Striking differentiation of sub-populations within a genetically homogeneous isolate (Ogliastra) in Sardinia as revealed by mtDNA analysis

Since the reduced genetic diversity found in isolates should simplify the study of complex traits, analyses of patterns of homogeneity within populations are of particular interest. We analysed the mtDNA haplogroups and hypervariable segment I (HVS-I) sequences of 475 individuals from a geographically restricted and isolated area (Ogliastra) within Sardinia, comprehending 175 random samples from 20 out of 23 villages. The remaining 300 subjects were chosen from the other three villages, Talana, Urzulei and Perdasdefogu, by sampling all maternal lineages. A comparison with other European populations reveals that Ogliastra ranks among the most genetically homogenous population and that it has been small and isolated throughout its history. The lack of variation and the high genetic homogeneity indicate that an important founder event and a demographic expansion took place during the Neolithic (~ 7,700 years before present) in Ogliastra's mtDNA gene pool. We present highly resolved phylogenetic networks for Ogliastra and for the three sub-isolates. MtDNA differentiation in the sub-populations versus Ogliastra is revealed by a strong demarcation in their genetic pools due to distinctive founder effects and genetic drift. We found that genetic homogeneity strictly depends on a scale factor in population size and on sampling methodology. The outstanding homogeneity and the reduced female gene pool observed in Ogliastra, in the European context, hide an extremely marked differentiation in sub-isolates originated from the same archaic population. Although Ogliastra can be considered a genetically homogeneous isolate, small villages' divergent genetic histories underline the importance of more systematic analysis of DNA variation between and within populations.C. Fraumene and E. Petretto contributed equally to this work     

High Differentiation among Eight Villages in a Secluded Area of Sardinia Revealed by Genome-Wide High Density SNPs Analysis

To better design association studies for complex traits in isolated populations it''s important to understand how history and isolation moulded the genetic features of different communities. Population isolates should not “a priori” be considered homogeneous, even if the communities are not distant and part of a small region. We studied a particular area of Sardinia called Ogliastra, characterized by the presence of several distinct villages that display different history, immigration events and population size. Cultural and geographic isolation characterized the history of these communities. We determined LD parameters in 8 villages and defined population structure through high density SNPs (about 360 K) on 360 unrelated people (45 selected samples from each village). These isolates showed differences in LD values and LD map length. Five of these villages show high LD values probably due to their reduced population size and extreme isolation. High genetic differentiation among villages was detected. Moreover population structure analysis revealed a high correlation between genetic and geographic distances. Our study indicates that history, geography and biodemography have influenced the genetic features of Ogliastra communities producing differences in LD and population structure. All these data demonstrate that we can consider each village an isolate with specific characteristics. We suggest that, in order to optimize the study design of complex traits, a thorough characterization of genetic features is useful to identify the presence of sub-populations and stratification within genetic isolates.     

High resolution analysis and phylogenetic network construction using complete mtDNA sequences in sardinian genetic isolates

For mitochondrial phylogenetic analysis, the best result comes from complete sequences. We therefore decided to sequence the entire mitochondrial DNA (mtDNA) (coding and D-loop regions) of 63 individuals selected in 3 small Ogliastra villages, an isolated area of eastern Sardinia: Talana, Urzulei, and Perdasdefogu. We studied at least one individual for each of the most frequent maternal genealogical lineages belonging to haplogroups H, V, J, K, T, U, and X. We found in our 63 samples, 172 and 69 sequence changes in the coding and in the D-loop region, respectively. Thirteen out of 172 sequence changes in the coding region are novel. It is our hypothesis that some of them are characteristic of the Ogliastra region and/or Sardinia. We reconstructed the phylogenetic network of the 63 complete mtDNA sequences for the 3 villages. We also drew a network including a large number of European sequences and calculated various indices of genetic diversity in Ogliastra. It appears that these small populations remained extremely isolated and genetically differentiated compared with other European populations. We also identified in our samples a never previously described subhaplogroup, U5b3, which seems peculiar to the Ogliastra region.     

Analysis of genome-wide structure, diversity and fine mapping of Mendelian traits in traditional and village chickens

Extensive phenotypic variation is a common feature among village chickens found throughout much of the developing world, and in traditional chicken breeds that have been artificially selected for traits such as plumage variety. We present here an assessment of traditional and village chicken populations, for fine mapping of Mendelian traits using genome-wide single-nucleotide polymorphism (SNP) genotyping while providing information on their genetic structure and diversity. Bayesian clustering analysis reveals two main genetic backgrounds in traditional breeds, Kenyan, Ethiopian and Chilean village chickens. Analysis of linkage disequilibrium (LD) reveals useful LD (r(2) ≥ 0.3) in both traditional and village chickens at pairwise marker distances of ~10 Kb; while haplotype block analysis indicates a median block size of 11-12 Kb. Association mapping yielded refined mapping intervals for duplex comb (Gga 2:38.55-38.89 Mb) and rose comb (Gga 7:18.41-22.09 Mb) phenotypes in traditional breeds. Combined mapping information from traditional breeds and Chilean village chicken allows the oocyan phenotype to be fine mapped to two small regions (Gga 1:67.25-67.28 Mb, Gga 1:67.28-67.32 Mb) totalling ~75 Kb. Mapping the unmapped earlobe pigmentation phenotype supports previous findings that the trait is sex-linked and polygenic. A critical assessment of the number of SNPs required to map simple traits indicate that between 90 and 110K SNPs are required for full genome-wide analysis of haplotype block structure/ancestry, and for association mapping in both traditional and village chickens. Our results demonstrate the importance and uniqueness of phenotypic diversity and genetic structure of traditional chicken breeds for fine-scale mapping of Mendelian traits in the species, with village chicken populations providing further opportunities to enhance mapping resolutions.     

Linkage disequilibrium between microsatellite markers in the Swedish Sami relative to a worldwide selection of populations

The pattern of linkage disequilibrium (LD) is affected by a number of factors, including population demography. High LD is seen in populations with a relatively limited and constant size, presumably because of genetic drift. We have examined the extent of LD among over 300 genome-wide pattern microsatellite loci in 29 populations from around the world. The pattern of LD varied between populations, with a larger extent of LD in populations with limited size relative to larger populations. In addition, the LD between 88 less well-spaced microsatellite markers from 10 different genomic regions was examined in the Sami compared with the general Swedish population. For these markers, increased LD extending up to 5 Mb was detected in the Sami. The amount of LD also differed between the chromosomal regions. The amount of LD in the Sami makes this population suitable for the mapping of complex genetic traits.Åsa Johansson, Veronika Vavruch-Nilsson contributed equally to the report     

Archival,demographic and genetic studies define a Sardinian sub-isolate as a suitable model for mapping complex traits

Genetic isolates represent exceptional resources for the mapping of complex traits but not all isolates are similar. We have selected a genetic and cultural isolate, the village of Talana from an isolated area of Sardinia, and propose that this population is suitable for the mapping of complex traits. A wealth of historical and archive data allowed the reconstruction of the demographic and genealogical history of the village. Key features of the population, which has grown slowly with no significant immigration, were defined by using a combination of historical, demographic and genetic studies. The genealogy of each Talana inhabitant was reconstructed and the main maternal and paternal lineages of the village were defined. Haplotype and phylogenetic analyses of the Y chromosome and characterisation of mitochondrial DNA haplogroups were used to determine the number of ancestral village founders. The extent of linkage disequilibrium (LD) was evaluated by the analysis of several microsatellites in chromosomal region Xq13.3, which was previously used to asses the extension of LD. Genealogical reconstructions were confirmed and reinforced by the genetic analyses, since some lineages were found to have merged prior to the beginning of the archival records, suggesting an even smaller number of founders than initially predicted. About 80% of the present-day population appears to derive from eight paternal and eleven maternal ancestral lineages. LD was found to span, on average, a 5-Mb region in Xq13.3. This suggests the possibility of identifying identical-by-descent regions associated with complex traits in a genome-wide search by using a low-density marker map. The present study emphasises the importance of combining genetic studies with genealogical and historical information.     

Gm and Km allotypes in four Sardinian population samples

Serum samples of 683 unrelated male and female individuals of four Sardinian population samples (Sassari, Nuoro, Oristano and Cagliari) were typed for G 1 m (1, 2, 3, 17), G 3 m (5, 6, 10, 11, 13, 14, 15, 16, 21, 26), and Km (1). Phenotype, haplotype (Gm), and allele frequencies (Km), respectively, show a remarkable variability between these four population samples. Comparisons with other Italian populations reveal the considerable genetic difference of the Sardinians, which is in particular caused by the presence of the haplotype Gm^{1, 3, 5, 10, 11, 13, 14, 26} in them. This haplotype is quite uncommon in Europeans and may reflect gene flow from Eastern populations (Phoenicians?) who came to this island in ancient-history times.     

Moroccan Leishmania infantum: Genetic Diversity and Population Structure as Revealed by Multi-Locus Microsatellite Typing

Leishmania infantum causes Visceral and cutaneous leishmaniasis in northern Morocco. It predominantly affects children under 5 years with incidence of 150 cases/year. Genetic variability and population structure have been investigated for 33 strains isolated from infected dogs and humans in Morocco. A multilocus microsatellite typing (MLMT) approach was used in which a MLMtype based on size variation in 14 independent microsatellite markers was compiled for each strain. MLMT profiles of 10 Tunisian, 10 Algerian and 21 European strains which belonged to zymodeme MON-1 and non-MON-1 according to multilocus enzyme electrophoresis (MLEE) were included for comparison. A Bayesian model-based approach and phylogenetic analysis inferred two L.infantum sub-populations; Sub-population A consists of 13 Moroccan strains grouped with all European strains of MON-1 type; and sub-population B consists of 15 Moroccan strains grouped with the Tunisian and Algerian MON-1 strains. Theses sub-populations were significantly different from each other and from the Tunisian, Algerian and European non MON-1 strains which constructed one separate population. The presence of these two sub-populations co-existing in Moroccan endemics suggests multiple introduction of L. infantum from/to Morocco; (1) Introduction from/to the neighboring North African countries, (2) Introduction from/to the Europe. These scenarios are supported by the presence of sub-population B and sub-population A respectively. Gene flow was noticed between sub-populations A and B. Five strains showed mixed A/B genotypes indicating possible recombination between the two populations. MLMT has proven to be a powerful tool for eco-epidemiological and population genetic investigations of Leishmania.     

Linkage disequilibrium analysis in the genetically isolated Norfolk Island population

Norfolk Island is a human genetic isolate, possessing unique population characteristics that could be utilized for complex disease gene localization. Our intention was to evaluate the extent and strength of linkage disequilibrium (LD) in the Norfolk isolate by investigating markers within Xq13.3 and the NOS2A gene encoding the inducible nitric oxide synthase. A total of six microsatellite markers spanning approximately 11 Mb were assessed on chromosome Xq13.3 in a group of 56 men from Norfolk Island. Additionally, three single nucleotide polymorphisms (SNPs) localizing to the NOS2A gene were analyzed in a subset of the complex Norfolk pedigree. With the exception of two of the marker pairs, one of which is the most distantly spaced marker, all the Xq13.3 marker pairs were found to be in significant LD indicating that LD extends up to 9.5-11.5 Mb in the Norfolk Island population. Also, all SNPs studied showed significant LD in both Norfolk Islanders and Australian Caucasians, with two of the marker pairs in complete LD in the Norfolk population only. The Norfolk Island study population possesses a unique set of characteristics including founder effect, geographical isolation, exhaustive genealogical information and phenotypic data of use to cardiovascular disease risk traits. With LD extending up to 9.5-11 Mb, the Norfolk isolate should be a powerful resource for the localization of complex disease genes.     

Linkage disequilibrium,SNP frequency change due to selection,and association mapping in popcorn chromosome regions containing QTLs for quality traits

The objectives of this study were to assess linkage disequilibrium (LD) and selection-induced changes in single nucleotide polymorphism (SNP) frequency, and to perform association mapping in popcorn chromosome regions containing quantitative trait loci (QTLs) for quality traits. Seven tropical and two temperate popcorn populations were genotyped for 96 SNPs chosen in chromosome regions containing QTLs for quality traits. The populations were phenotyped for expansion volume, 100-kernel weight, kernel sphericity, and kernel density. The LD statistics were the difference between the observed and expected haplotype frequencies (D), the proportion of D relative to the expected maximum value in the population, and the square of the correlation between the values of alleles at two loci. Association mapping was based on least squares and Bayesian approaches. In the tropical populations, D-values greater than 0.10 were observed for SNPs separated by 100-150 Mb, while most of the D-values in the temperate populations were less than 0.05. Selection for expansion volume indirectly led to increase in LD values, population differentiation, and significant changes in SNP frequency. Some associations were observed for expansion volume and the other quality traits. The candidate genes are involved with starch, storage protein, lipid, and cell wall polysaccharides synthesis.     

SNP identification, linkage disequilibrium, and haplotype analysis for a 200-kb genomic region in a Korean population

Understanding patterns of linkage disequilibrium (LD) across genomes may facilitate association mapping studies to localize genetic variants influencing complex diseases, a recognition that led to the International Haplotype Mapping Project (HapMap). Divergent patterns of haplotype frequency and LD across global populations require that the HapMap database be supplemented with haplotype and LD data from additional populations. We conducted a pilot study of the LD and haplotype structure of a genomic region in a Korean population. A total of 165 SNPs were identified in a 200-kb region of 22q13.2 by direct sequencing. Unphased genotype data were generated for 76 SNPs in 90 unrelated Korean individuals. LD, haplotype diversity, and recombination rates were assessed in this region and compared with the HapMap database. The pattern of LD and haplotype frequencies of Korean samples showed a high degree of similarity with Japanese data. There was a strong correlation between high LD and low recombination frequency in this region. We found considerable similarities in local LD patterns between three Asian populations (Han Chinese, Japanese, and Korean) and the CEPH population. Haplotype frequencies were, however, significantly different between them. Our results should further the understanding of distinctive Korean genomic features and assist in designing appropriate association studies.     

Genomewide linkage disequilibrium mapping of severe bipolar disorder in a population isolate

Genomewide association studies may offer the best promise for genetic mapping of complex traits. Such studies in outbred populations require very densely spaced single-nucleotide polymorphisms. In recently founded population isolates, however, extensive linkage disequilibrium (LD) may make these studies feasible with currently available sets of short tandem repeat markers, spaced at intervals as large as a few centimorgans. We report the results of a genomewide association study of severe bipolar disorder (BP-I), using patients from the isolated population of the central valley of Costa Rica. We observed LD with BP-I on several chromosomes; the most striking results were in proximal 8p, a region that has previously shown linkage to schizophrenia. This region could be important for severe psychiatric disorders, rather than for a specific phenotype.     

Genome association studies of complex diseases by case-control designs

One way to perform linkage-disequilibrium (LD) mapping of genetic traits is to use single markers. Since dense marker maps-such as single-nucleotide polymorphism and high-resolution microsatellite maps-are available, it is natural and practical to generalize single-marker LD mapping to high-resolution haplotype or multiple-marker LD mapping. This article investigates high-resolution LD-mapping methods, for complex diseases, based on haplotype maps or microsatellite marker maps. The objective is to explore test statistics that combine information from haplotype blocks or multiple markers. Based on two coding methods, genotype coding and haplotype coding, Hotelling's T2 statistics TG and TH are proposed to test the association between a disease locus and two haplotype blocks or two markers. The validity of the two T2 statistics is proved by theoretical calculations. A statistic TC, an extension of the traditional chi2 method of comparing haplotype frequencies, is introduced by simply adding the chi2 test statistics of the two haplotype blocks together. The merit of the three methods is explored by calculation and comparison of power and of type I errors. In the presence of LD between the two blocks, the type I error of TC is higher than that of TH and TG, since TC ignores the correlation between the two blocks. For each of the three statistics, the power of using two haplotype blocks is higher than that of using only one haplotype block. By power comparison, we notice that TC has higher power than that of TH, and TH has higher power than that of TG. In the absence of LD between the two blocks, the power of TC is similar to that of TH and higher than that of TG. Hence, we advocate use of TH in the data analysis. In the presence of LD between the two blocks, TH takes into account the correlation between the two haplotype blocks and has a lower type I error and higher power than TG. Besides, the feasibility of the methods is shown by sample-size calculation.     

Genetic diversity of three donkey populations in the Croatian coastal region

Genetic variation in three Croatian donkey populations, Istrian (IS), North Adriatic (NA) and Littoral-Dinaric (LD), was analysed using eight microsatellite loci and by sequence and SSCP analysis of the proximal portion of the mtDNA D-loop region. The analysis of microsatellite loci revealed observed heterozygosities in the range of 0.37 (MPZ002 in LD) to 0.85 (AHT21 in LD) and polymorphic information content values in the range of 0.36 (MPZ002 in NA) to 0.78 (AHT21 in LD). The overall probability of exclusion was 0.991. Two populations (IS and NA) were closely related (Fst=0.0034), whereas genetic distances between IS and LD (Fst=0.021) and NA and LD (Fst=0.027) were higher. Using AMOVA, 97.6% of the total genetic variance was portioned within populations, while 2.7% was portioned between the Littoral-Dinaric population and the Istrian/North Adriatic population group. Sequencing of the proximal part of the mtDNA D-loop region revealed 36 polymorphic sites representing 19 haplotypes which clustered into three haplotype groups (Y, W, Ws). Only the Y haplotype was found in the IS population which is characterized by a large body size. Haplotypes W and Ws were found in the NA and LD populations which include smaller animals. All three haplotypes were found in the LD population, indicating sporadic migration events from the IS into LD donkey population.     

On harvesting a structured ungulate population

Variation in demographic rates within a spatially structured population could have important consequences for management decisions, harvesting strategies and offtake rates. Although there is a growing body of evidence suggesting that demographic rates vary within populations over a range of spatial scales, there has been little research investigating the consequences of this variation for population management. In this paper, data on the dynamics of two female red deer sub-populations on Rum are analysed, and evidence is presented for differences between the fecundity and mortality rates of the two sub-populations. A simple harvesting model is developed to represent the dynamics of the two sub-populations, including density-independent migration between sub-populations and spatially correlated environmental variability. The highest monetary yield in the model is obtained by harvesting the more resilient sub-population at a higher rate. Surprisingly the losses involved in harvesting both sub-populations at the same rate are insignificant. However, if migration were density-dependent, the size of one sub-population would be more relevant to harvesting policy for the other sub-population. The results of this empirical study are compared to theoretical work on spatially structured populations; it is shown that when a species has complex age- and sex-structured population dynamics, previous theoretical results may not hold.     

Population differences in levels of linkage disequilibrium in the wild

Information about the levels of linkage disequilibrium (LD) in wild animal populations is still limited, and this is true particularly with respect to possible interpopulation variation in the levels of LD. We compared the levels and extent of LD at the genome‐wide scale in three Siberian jay (Perisoreus infaustus) populations, two of which (Kuusamo and Ylläs) represented outbred populations within the main distribution area of the species, whereas the third (Suupohja) was a semi‐isolated, partially inbred population at the margin of the species’ distribution area. Although extensive long‐range LD (>20 cM) was observed in all three populations, LD generally decayed to background levels at a distance of 1–5 cM or c. 200–600 kb. The degree and extent of LD differed markedly between populations but aligned closely with both observed levels of within‐population genetic variation and expectations based on population history. The levels of LD were highest in the most inbred population with strong population substructure (Suupohja), compared with the two outbred populations. Furthermore, the decay of LD with increasing distance was slower in Suupohja, compared with the other two populations. By demonstrating that levels of LD can vary greatly over relatively short geographical distances within a species, these results suggest that prospects for association mapping differ from population to population. In this example, the prospects are best in the Suupohja population, given that minimized marker genotyping and a minimum marker spacing of 1–5 cM (c. 200–600 kb) would be sufficient for a whole genome scan for detecting QTL.     

Red cell and serum protein polymorphisms in Sardinia

Four Sardinian population samples, from the provinces of Cagliari, Oristano, Nuoro and Sassari, were studied with regard to the erythrocyte enzyme systems ACP, ESD, PGM1, ADA, AK, 6PGD and Dia, and to the serum protein systems GC and C3. The findings showed a rather high degree of genetic heterogeneity of the Sardinians compared to the other populations from the Mediterranean area (Continental Italy, Sicily, Spain, North Africa).     

Environmental and Genetic Contribution to Hypertension Prevalence: Data from an Epidemiological Survey on Sardinian Genetic Isolates

Background and Objectives

Hypertension represents a major cause of cardiovascular morbidity and mortality worldwide but its prevalence has been shown to vary in different countries. The reasons for such differences are still matter of debate, the relative contributions given by environmental and genetic factors being still poorly defined. We estimated the current prevalence, distribution and determinants of hypertension in isolated Sardinian populations and also investigated the environmental and genetic contribution to hypertension prevalence taking advantage of the characteristics of such populations.

Methods and Results

An epidemiological survey with cross-sectional design was carried out measuring blood pressure in 9845 inhabitants of 10 villages of Ogliastra region between 2002 and 2008. Regression analysis for assessing blood pressure determinants and variance component models for estimating heritability were performed. Overall 38.8% of this population had hypertension, its prevalence varying significantly by age, sex and among villages taking into account age and sex structure of their population. About 50% of hypertensives had prior cardiovascular disease. High blood pressure was independently associated with age, obesity related factors, heart rate, total cholesterol, alcohol consumption, low education and smoking status, all these factors contributing more in women than in men. Heritability was 27% for diastolic and 36% for systolic blood pressure, its contribution being significantly higher in men (57%) than in women (46%). Finally, the genetic correlation between systolic and diastolic blood pressure was 0.74, indicating incomplete pleiotropy.

Conclusion

Genetic factors involved in the expression of blood pressure traits account for about 30% of the phenotypic variance, but seem to play a larger role in men; comorbidities and environmental factors remain of predominant importance, but seem to contribute much more in women.     

Identification of positive selection signatures in pigs by comparing linkage disequilibrium variances

Selection affects the patterns of linkage disequilibrium (LD) around the site of a beneficial allele with an increase in LD among the hitchhiking alleles. Comparing the differences in regional LD between pig populations could help to identify putative genomic regions with potential adaptations for economic traits. In this study, using Illumina Porcine SNP60K BeadChip genotyping data from 207 Chinese indigenous, 117 South American village and 408 Large White pigs, we estimated the variation of genome‐wide LD between populations using the varld program. The top 0.1% standardized VarLD scores were used as a criterion for all comparisons, and compared with LD blocks, a total of four selection signatures on Sus scrofa chromosome (SSC) 7, 9, 13 and 14 were identified in all populations. These signatures overlapped with quantitative trait loci for linoleic acid content, age at puberty, number of muscle fibers per unit area, hip structure and body weight traits in pigs. Among them, one of the signatures (56.5–56.6 Mb on SSC7) in Large White pigs harbored the ADAMTSL3 gene, which is known to affect body length. The findings of this study seem to point toward recent selection in different pig populations. Further investigations are encouraged to confirm the selection signatures detected by varld in the present study.     

Linkage Disequilibrium and Demographic History of Wild and Domestic Canids

Assessing the extent of linkage disequilibrium (LD) in natural populations of a nonmodel species has been difficult due to the lack of available genomic markers. However, with advances in genotyping and genome sequencing, genomic characterization of natural populations has become feasible. Using sequence data and SNP genotypes, we measured LD and modeled the demographic history of wild canid populations and domestic dog breeds. In 11 gray wolf populations and one coyote population, we find that the extent of LD as measured by the distance at which r² = 0.2 extends <10 kb in outbred populations to >1.7 Mb in populations that have experienced significant founder events and bottlenecks. This large range in the extent of LD parallels that observed in 18 dog breeds where the r² value varies from ~20 kb to >5 Mb. Furthermore, in modeling demographic history under a composite-likelihood framework, we find that two of five wild canid populations exhibit evidence of a historical population contraction. Five domestic dog breeds display evidence for a minor population contraction during domestication and a more severe contraction during breed formation. Only a 5% reduction in nucleotide diversity was observed as a result of domestication, whereas the loss of nucleotide diversity with breed formation averaged 35%.