Population structure in admixed populations: effect of admixture dynamics on the pattern of linkage disequilibrium

Gene flow between genetically distinct populations creates linkage disequilibrium (admixture linkage disequilibrium [ALD]) among all loci (linked and unlinked) that have different allele frequencies in the founding populations. We have explored the distribution of ALD by using computer simulation of two extreme models of admixture: the hybrid-isolation (HI) model, in which admixture occurs in a single generation, and the continuous-gene-flow (CGF) model, in which admixture occurs at a steady rate in every generation. Linkage disequilibrium patterns in African American population samples from Jackson, MS, and from coastal South Carolina resemble patterns observed in the simulated CGF populations, in two respects. First, significant association between two loci (FY and AT3) separated by 22 cM was detected in both samples. The retention of ALD over relatively large (>10 cM) chromosomal segments is characteristic of a CGF pattern of admixture but not of an HI pattern. Second, significant associations were also detected between many pairs of unlinked loci, as observed in the CGF simulation results but not in the simulated HI populations. Such a high rate of association between unlinked markers in these populations could result in false-positive linkage signals in an admixture-mapping study. However, we demonstrate that by conditioning on parental admixture, we can distinguish between true linkage and association resulting from shared ancestry. Therefore, populations with a CGF history of admixture not only are appropriate for admixture mapping but also have greater power for detection of linkage disequilibrium over large chromosomal regions than do populations that have experienced a pattern of admixture more similar to the HI model, if methods are employed that detect and adjust for disequilibrium caused by continuous admixture.     

Analysis of genomic admixture in Uyghur and its implication in mapping strategy

The Uyghur (UIG) population, settled in Xinjiang, China, is a population presenting a typical admixture of Eastern and Western anthropometric traits. We dissected its genomic structure at population level, individual level, and chromosome level by using 20,177 SNPs spanning nearly the entire chromosome 21. Our results showed that UIG was formed by two-way admixture, with 60% European ancestry and 40% East Asian ancestry. Overall linkage disequilibrium (LD) in UIG was similar to that in its parental populations represented in East Asia and Europe with regard to common alleles, and UIG manifested elevation of LD only within 500 kb and at a level of 0.1 相似文献   

Bayesian analyses of admixture in wild and domestic cats (Felis silvestris) using linked microsatellite loci

Methods recently developed to infer population structure and admixture mostly use individual genotypes described by unlinked neutral markers. However, Hardy-Weinberg and linkage disequilibria among independent markers decline rapidly with admixture time, and the admixture signals could be lost in a few generations. In this study, we aimed to describe genetic admixture in 182 European wild and domestic cats (Felis silvestris), which hybridize sporadically in Italy and extensively in Hungary. Cats were genotyped at 27 microsatellites, including 21 linked loci mapping on five distinct feline linkage groups. Genotypes were analysed with structure 2.1, a Bayesian procedure designed to model admixture linkage disequilibrium, which promises to assess efficiently older admixture events using tightly linked markers. Results showed that domestic and wild cats sampled in Italy were split into two distinct clusters with average proportions of membership Q > 0.90, congruent with prior morphological identifications. In contrast, free-living cats sampled in Hungary were assigned partly to the domestic and the wild cat clusters, with Q < 0.50. Admixture analyses of individual genotypes identified, respectively, 5/61 (8%), and 16-20/65 (25-31%) hybrids among the Italian wildcats and Hungarian free-living cats. Similar results were obtained in the past using unlinked loci, although the new linked markers identified additional admixed wildcats in Italy. Linkage analyses confirm that hybridization is limited in Italian, but widespread in Hungarian wildcats, a population that is threatened by cross-breeding with free-ranging domestic cats. The total panel of 27 loci performed better than the linked loci alone in the identification of domestic and known hybrid cats, suggesting that a large number of linked plus unlinked markers can improve the results of admixture analyses. Inferred recombination events led to identify the population of origin of chromosomal segments, suggesting that admixture mapping experiments can be designed also in wild populations.     

Sensitive Detection of Chromosomal Segments of Distinct Ancestry in Admixed Populations

Identifying the ancestry of chromosomal segments of distinct ancestry has a wide range of applications from disease mapping to learning about history. Most methods require the use of unlinked markers; but, using all markers from genome-wide scanning arrays, it should in principle be possible to infer the ancestry of even very small segments with exquisite accuracy. We describe a method, HAPMIX, which employs an explicit population genetic model to perform such local ancestry inference based on fine-scale variation data. We show that HAPMIX outperforms other methods, and we explore its utility for inferring ancestry, learning about ancestral populations, and inferring dates of admixture. We validate the method empirically by applying it to populations that have experienced recent and ancient admixture: 935 African Americans from the United States and 29 Mozabites from North Africa. HAPMIX will be of particular utility for mapping disease genes in recently admixed populations, as its accurate estimates of local ancestry permit admixture and case-control association signals to be combined, enabling more powerful tests of association than with either signal alone.     

Reconstructing genetic ancestry blocks in admixed individuals

A chromosome in an individual of recently admixed ancestry resembles a mosaic of chromosomal segments, or ancestry blocks, each derived from a particular ancestral population. We consider the problem of inferring ancestry along the chromosomes in an admixed individual and thereby delineating the ancestry blocks. Using a simple population model, we infer gene-flow history in each individual. Compared with existing methods, which are based on a hidden Markov model, the Markov-hidden Markov model (MHMM) we propose has the advantage of accounting for the background linkage disequilibrium (LD) that exists in ancestral populations. When there are more than two ancestral groups, we allow each ancestral population to admix at a different time in history. We use simulations to illustrate the accuracy of the inferred ancestry as well as the importance of modeling the background LD; not accounting for background LD between markers may mislead us to false inferences about mixed ancestry in an indigenous population. The MHMM makes it possible to identify genomic blocks of a particular ancestry by use of any high-density single-nucleotide-polymorphism panel. One application of our method is to perform admixture mapping without genotyping special ancestry-informative-marker panels.     

Consistent long-range linkage disequilibrium generated by admixture in a Bantu-Semitic hybrid population

Both the optimal marker density for genome scans in case-control association studies and the appropriate study design for the testing of candidate genes depend on the genomic pattern of linkage disequilibrium (LD). In this study, we provide the first conclusive demonstration that the diverse demographic histories of human populations have produced dramatic differences in genomewide patterns of LD. Using a panel of 66 markers spanning the X chromosome, we show that, in the Lemba, a Bantu-Semitic hybrid population, markers 2 cM. Moreover, analysis of Bantu and Ashkenazi populations as putative parental populations of the Lemba shows a significant relationship between allele-frequency differentials and the LD observed in the Lemba, which demonstrates that much of the excess LD is due to admixture. Our results suggest that demographic history has such a profound effect on LD that it will not be possible to predict patterns a priori but that it will be necessary to empirically evaluate the patterns in all populations of interest.     

Comparison of Linkage Disequilibrium in Elite European Maize Inbred Lines using AFLP and SSR Markers

Application of association mapping to plant breeding populations has the potential to revolutionize plant genetics. The main objectives of this study were to (i) investigate the extent and genomic distribution of linkage disequilibrium (LD) between pairs of amplified fragment length polymorphism (AFLP) markers, (ii) compare these results with those obtained with simple sequence repeat (SSR) markers, and (iii) compare the usefulness of AFLP and SSR markers for genomewide association mapping in plant breeding populations. We examined LD in a cross-section of 72 European elite inbred lines genotyped with 452 AFLP and 93 SSR markers. LD was significant (p < 0.05) for about 15% of the AFLP marker pairs and for about 49% of the SSR marker pairs in each of the two germplasm groups, flint and dent. In both germplasm groups the ratio of linked to unlinked loci pairs in LD was higher for AFLPs than for SSRs. The observation of LD due to linkage for both marker types suggested that genome-wide association mapping should be possible using either AFLPs or SSRs. The results of our study indicated that SSRs should be favored over AFLPs but the opposite applies to populations with a long history of recombination.     

Ethnic-difference markers for use in mapping by admixture linkage disequilibrium

Mapping by admixture linkage disequilibrium (MALD) is a potentially powerful technique for the mapping of complex genetic diseases. The practical requirements of this method include (a) a set of markers spanning the genome that have large allele-frequency differences between the parental ethnicities contributing to the admixed population and (b) an understanding of the extent of admixture in the study population. To this end, a DNA-pooling technique was used to screen microsatellite and diallelic insertion/deletion markers for allele-frequency differences between putative representatives of the parental populations of the admixed Mexican American (MA) and African American (AA) populations. Markers with promising pooled differences were then confirmed by individual genotyping in both the parental and admixed populations. For the MA population, screening of >600 markers identified 151 ethnic-difference markers (EDMs) with delta>0.30 (where delta is the absolute value of each allele-frequency difference between two populations, summed over all marker alleles and divided by two) that are likely to be useful for MALD analysis. For the AA population, analysis of >400 markers identified 97 EDMs. In addition, individual genotyping of these markers in Pima Amerindians, Yavapai Amerindians, European American (EA) individuals, Africans from Zimbabwe, MA individuals, and AA individuals, as well as comparison to the CEPH genotyping set, suggests that the differences between subpopulations of an ethnicity are small for many markers with large interethnic differences. Estimates of admixture that are based on individual genotyping of these markers are consistent with a 60% EA:40% Amerindian contribution to MA populations and with a 20% EA:80% African contribution to AA populations. Taken together, these data suggest that EDMs with large interpopulation and small intrapopulation differences can be readily identified for MALD studies in both AA and MA populations.     

Linkage disequilibrium in European elite maize germplasm investigated with SSRs

Information about the extent and genomic distribution of linkage disequilibrium (LD) is of fundamental importance for association mapping. The main objectives of this study were to (1) investigate genetic diversity within germplasm groups of elite European maize (Zea mays L.) inbred lines, (2) examine the population structure of elite European maize germplasm, and (3) determine the extent and genomic distribution of LD between pairs of simple sequence repeat (SSR) markers. We examined genetic diversity and LD in a cross section of European and US elite breeding material comprising 147 inbred lines genotyped with 100 SSR markers. For gene diversity within each group, significant (P<0.05) differences existed among the groups. The LD was significant (P<0.05) for 49% of the SSR marker pairs in the 80 flint lines and for 56% of the SSR marker pairs in the 57 dent lines. The ratio of linked to unlinked loci in LD was 1.1 for both germplasm groups. The high incidence of LD suggests that the extent of LD between SSR markers should allow the detection of marker-phenotype associations in a genome scan. However, our results also indicate that a high proportion of the observed LD is generated by forces, such as relatedness, population stratification, and genetic drift, which cause a high risk of detecting false positives in association mapping.     

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     

A genome-wide analysis of admixture in Uyghurs and a high-density admixture map for disease-gene discovery

Following up on our previous study, we conducted a genome-wide analysis of admixture for two Uyghur population samples (HGDP-UG and PanAsia-UG), collected from the northern and southern regions of Xinjiang in China, respectively. Both HGDP-UG and PanAsia-UG showed a substantial admixture of East-Asian (EAS) and European (EUR) ancestries, with an empirical estimation of ancestry contribution of 53:47 (EAS:EUR) and 48:52 for HGDP-UG and PanAsia-UG, respectively. The effective admixture time under a model with a single pulse of admixture was estimated as 110 generations and 129 generations, or admixture events occurred about 2200 and 2580 years ago for HGDP-UG and PanAsia-UG, respectively, assuming an average of 20 yr per generation. Despite Uyghurs' earlier history compared to other admixture populations, admixture mapping, holds promise for this population, because of its large size and its mixture of ancestry from different continents. We screened multiple databases and identified a genome-wide single-nucleotide polymorphism panel that can distinguish EAS and EUR ancestry of chromosomal segments in Uyghurs. The panel contains 8150 ancestry-informative markers (AIMs) showing large frequency differences between EAS and EUR populations (F_ST > 0.25, mean F_ST = 0.43) but small frequency differences (7999 AIMs validated) within both populations (F_ST < 0.05, mean F_ST < 0.01). We evaluated the effectiveness of this admixture map for localizing disease genes in two Uyghur populations. To our knowledge, our map constitutes the first practical resource for admixture mapping in Uyghurs, and it will enable studies of diseases showing differences in genetic risk between EUR and EAS populations.     

Evidence for a large-scale population structure of Arabidopsis thaliana from genome-wide single nucleotide polymorphism markers

Population-based methods for the genetic mapping of adaptive traits and the analysis of natural selection require that the population structure and demographic history of a species are taken into account. We characterized geographic patterns of genetic variation in the model plant Arabidopsis thaliana by genotyping 115 genome-wide single nucleotide polymorphism (SNP) markers in 351 accessions from the whole species range using a matrix-assisted laser desorption/ionization time-of-flight assay, and by sequencing of nine unlinked short genomic regions in a subset of 64 accessions. The observed frequency distribution of SNPs is not consistent with a constant-size neutral model of sequence polymorphism due to an excess of rare polymorphisms. There is evidence for a significant population structure as indicated by differences in genetic diversity between geographic regions. Accessions from Central Asia have a low level of polymorphism and an increased level of genome-wide linkage disequilibrium (LD) relative to accessions from the Iberian Peninsula and Central Europe. Cluster analysis with the structure program grouped Eurasian accessions into K=6 clusters. Accessions from the Iberian Peninsula and from Central Asia constitute distinct populations, whereas Central and Eastern European accessions represent admixed populations in which genomes were reshuffled by historical recombination events. These patterns likely result from a rapid postglacial recolonization of Eurasia from glacial refugial populations. Our analyses suggest that mapping populations for association or LD mapping should be chosen from regional rather than a species-wide sample or identified genetically as sets of individuals with similar average genetic distances. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Association mapping of morphological traits in wild and captive zebra finches: reliable within,but not between populations

Identifying causal genetic variants underlying heritable phenotypic variation is a long‐standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole‐genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild‐caught zebra finches. Then, we validated the most promising SNP–phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype–phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.     

Admixture and patterns of linkage disequilibrium in a free-living vertebrate population

Linkage disequilibrium (LD), a measure of nonrandom association of alleles at different loci, is of great interest to evolutionary geneticists as it can be used to help identify loci that explain phenotypic variation. Surveys of the extent of LD across genomes have been carried out in a number of systems, most notably humans and model organisms. However, studies of natural populations of vertebrates have rarely been performed. Here, we describe an investigation of LD in a free-living island population of red deer Cervus elaphus. Relatively high levels of LD extended several tens of centimorgans, and significant LD was frequently detected between unlinked markers. The magnitude of LD varied depending on how the population was sampled. It also varied across different chromosomes, and was shown to be a function of sample size, intermarker distance and marker heterozygosity. A recent admixture event in the population led to an ephemeral increase in LD. Association mapping may be possible in this population, although a high 'baseline' level of LD could lead to false positive associations between marker loci and a trait of interest.     

A covariance structure model for the admixture of binary genetic variation

I derive a covariance structure model for pairwise linkage disequilibrium (LD) between binary markers in a recently admixed population and use a generalized least-squares method to fit the model to two different data sets. Both linked and unlinked marker pairs are incorporated in the model. Under the model, a pairwise LD matrix is decomposed into two component matrices, one containing LD attributable to admixture, and another containing, in an aggregate form, LD specific to the populations forming the mixture. I use population genetics theory to show that the latter matrix has block-diagonal structure. For the data sets considered here, I show that the number of source populations can be determined by statistical inference on the canonical correlations of the sample LD matrix.     

A first insight into population structure and linkage disequilibrium in the US peanut minicore collection

Knowledge of genetic diversity, population structure, and degree of linkage disequilibrium (LD) in target association mapping populations is of great importance and is a prerequisite for LD-based mapping. In the present study, 96 genotypes comprising 92 accessions of the US peanut minicore collection, a component line of the tetraploid variety Florunner, diploid progenitors A. duranensis (AA) and A. ipaënsis (BB), and synthetic amphidiploid accession TxAG-6 were investigated with 392 simple sequence repeat (SSR) marker bands amplified using 32 highly-polymorphic SSR primer pairs. Both distance- and model-based (Bayesian) cluster analysis revealed the presence of structured diversity. In general, the wild-species accessions and the synthetic amphidiploid grouped separately from most minicore accessions except for COC155, and were eliminated from most subsequent analyses. UPGMA analysis divided the population into four subgroups, two major subgroups representing subspecies fastigiata and hypogaea, a third group containing individuals from each subspecies or possibly of mixed ancestry, and a fourth group, either consisting of COC155 alone if wild species were excluded, or of COC155, the diploid species, and the synthetic amphidiploid. Model-based clustering identified four subgroups- one each for fastigiata and hypogaea subspecies, a third consisting of individuals of both subspecies or of mixed ancestry predominantly from Africa or Asia, and a fourth group, consisting of individuals predominantly of var fastigiata, peruviana, and aequatoriana accessions from South America, including COC155. Analysis of molecular variance (AMOVA) revealed statistically-significant (P < 0.0001) genetic variance of 16.87% among subgroups. A total of 4.85% of SSR marker pairs revealed significant LD (at r² ≥ 0.1). Of the syntenic marker pairs separated by distances < 10 cM, 11–20 cM, 21–50 cM, and > 50 cM, 19.33, 5.19, 6.25 and 5.29% of marker pairs were found in strong LD (P ≤ 0.01), in accord with LD extending to great distances in self pollinated crops. A threshold value of r² > 0.035 was found to distinguish mean r² values of linkage distance groups statistically from the mean r² values of unlinked markers; LD was found to extend to 10 cM over the entire minicore collection by this criterion. However, there were large differences in r² values among marker pairs even among tightly-linked markers. The implications of these findings with regard to the possibility of using association mapping for detection of genome-wide SSR marker-phenotype association are discussed.     

Comparative AFLP mapping in two hexaploid oat populations

Amplified fragment length polymorphisms (AFLPs) can be used to quickly develop linkage maps in plant species and are especially useful for crops with large genomes like oat (Avena sativa L., 2n=6x=42). High reproducibility and consistency are crucial if AFLP linkage maps are employed for comparative mapping. We mapped AFLP markers in combination with restriction fragment length polymorphism (RFLP) markers in two recombinant inbred populations of hexaploid oat in two laboratories to test the consistency of AFLP markers in a polyploid crop. Eight primer combinations produced 102 and 121 scoreable AFLP markers in the respective populations. In a population from the cross Kanota×Ogle, AFLP markers were placed onto a RFLP reference map consisting of 32 linkage groups. Nineteen linkage groups from another population from the cross Kanota×Marion were assigned to the reference map using AFLP and RFLP markers homologous to those used in the Kanota× Ogle cross. Reproducibility of AFLP assays was high in both laboratories and between laboratories. The AFLP markers were well-distributed across the genome in both populations. Many AFLP markers tended to extend the distance between adjacent RFLP markers in linkage analysis. Of the 27 polymorphic AFLPs common in both populations, 20 mapped to homologous linkage groups, 4 were unlinked in at least one population, and 3 mapped to different linkage groups in the two crosses. We believe that 1 of the 3 markers that mapped to a different linkage group in the two populations mapped to homoeologous linkage groups. The linkage map of hexaploid oat is not yet complete, and genomic rearrangements such as translocations exist among cultivars and are likely to account for the remaining two non-syntenous mapping results. AFLPs provide not only a fast and powerful tool for mapping but could be useful in characterizing genomic structural variations among germplasms in hexaploid oat. Received: 17 December 1999 / Accepted: 28 July 2000     

Markers that discriminate between European and African ancestry show limited variation within Africa

Markers informative for ancestry are necessary for admixture mapping and improving case-control association analyses. In particular, African Americans are an admixed population for which genetic studies require accurately evaluating admixture. This will require markers that can be used in African Americans to determine if a given genomic region is of European or African ancestry. This report shows that, despite studies indicating high intra-African sequence variation, markers with large inter-ethnic differences have only small variations in allele distribution among divergent African populations and should be valuable for evaluating admixture in complex disease genetic studies.     

The extent of linkage disequilibrium in four populations with distinct demographic histories

The design and feasibility of whole-genome-association studies are critically dependent on the extent of linkage disequilibrium (LD) between markers. Although there has been extensive theoretical discussion of this, few empirical data exist. The authors have determined the extent of LD among 38 biallelic markers with minor allele frequencies >.1, since these are most comparable to the common disease-susceptibility polymorphisms that association studies aim to detect. The markers come from three chromosomal regions-1,335 kb on chromosome 13q12-13, 380 kb on chromosome 19q13.2, and 120 kb on chromosome 22q13.3-which have been extensively mapped. These markers were examined in approximately 1,600 individuals from four populations, all of European origin but with different demographic histories; Afrikaners, Ashkenazim, Finns, and East Anglian British. There are few differences, either in allele frequencies or in LD, among the populations studied. A similar inverse relationship was found between LD and distance in each genomic region and in each population. Mean D' is.68 for marker pairs <5 kb apart and is.24 for pairs separated by 10-20 kb, and the level of LD is not different from that seen in unlinked marker pairs separated by >500 kb. However, only 50% of marker pairs at distances <5 kb display sufficient LD (delta>.3) to be useful in association studies. Results of the present study, if representative of the whole genome, suggest that a whole-genome scan searching for common disease-susceptibility alleles would require markers spaced < or = 5 kb apart.