Argentine population genetic structure: large variance in Amerindian contribution

Argentine population genetic structure was examined using a set of 78 ancestry informative markers (AIMs) to assess the contributions of European, Amerindian, and African ancestry in 94 individuals members of this population. Using the Bayesian clustering algorithm STRUCTURE, the mean European contribution was 78%, the Amerindian contribution was 19.4%, and the African contribution was 2.5%. Similar results were found using weighted least mean square method: European, 80.2%; Amerindian, 18.1%; and African, 1.7%. Consistent with previous studies the current results showed very few individuals (four of 94) with greater than 10% African admixture. Notably, when individual admixture was examined, the Amerindian and European admixture showed a very large variance and individual Amerindian contribution ranged from 1.5 to 84.5% in the 94 individual Argentine subjects. These results indicate that admixture must be considered when clinical epidemiology or case control genetic analyses are studied in this population. Moreover, the current study provides a set of informative SNPs that can be used to ascertain or control for this potentially hidden stratification. In addition, the large variance in admixture proportions in individual Argentine subjects shown by this study suggests that this population is appropriate for future admixture mapping studies.     

Admixture analysis of a rural population of the state of Guerrero, Mexico

We studied 156 individuals of Native American descent from the city of Tlapa in the state of Guerrero in western Mexico. Most individuals' ethnicity was either Nahua, Mixtec, or Tlapanec, but self-identified Mestizos and individuals of mixed ethnicities were also included in the sample. We typed 24 autosomal, one Y-chromosome, and four mitochondrial ancestry-informative markers (AIMs) to estimate group and individual admixture proportions, and determine whether the admixture process involved directional gene flow between parental groups. When genetically defined (GD) Mestizos were excluded from the analysis, Native American ancestry represented approximately 98% of the population's gene pool, while European and West African ancestry represented approximately 1% each. Maternally inherited markers also showed an exceptionally high Native American contribution (98.5%), as did the paternally inherited marker, DYS199 (90.7%). We did not detect genetic structure in this population using these AIMs, which appears consistent with the homogeneity of the sample in terms of admixture proportions. The addition of GD Mestizos to the sample did not produce a considerable change in admixture estimates, but it had a major effect on population structure. These results show that the population of Tlapa in Guerrero, Mexico, has experienced little admixture with Europeans and/or West Africans. They also show that the impact of a small number of admixed individuals on an otherwise homogeneous population might have profound implications on subsequent ancestry/phenotype analysis and mapping strategies. We suggest that heterogeneity is a major characteristic of Mexican populations and, as a consequence, should not be disregarded when designing epidemiological studies of Mexican and Mexican American populations.     

Association of African Genetic Admixture with Resting Metabolic Rate and Obesity Among Women

Objective: To investigate the role of genetic admixture in explaining phenotypic variation in obesity‐related traits in a sample of African‐American women (n = 145) and to determine significant associations between obesity traits and admixture genetic markers. Research Methods and Procedures: Associations between genetic admixture and BMI, resting metabolic rate, fat mass, fat‐free mass, and bone mineral density were tested using linear regression considering the estimation of admixture by 1) a maximum‐likelihood approach (MLA) and 2) a Bayesian analysis. Results: Both the conservative MLA and the Bayesian approach support an association between African genetic admixture and BMI. Evidence for the associations of African genetic admixture with fat mass and fat‐free mass was supported by the Bayesian analysis; the MLA supported an association with bone mineral density. When the individual ancestry informative markers that were used to estimate admixture were tested for associations with BMI, significant associations were identified in chromosomes 1, 11, and 12. Discussion: These results provide evidence supporting the application of admixture mapping methods to the identification of genes that result in higher levels of obesity among African‐American women. Further research is needed to replicate and further explore these findings.     

Clumpak: a program for identifying clustering modes and packaging population structure inferences across K

The identification of the genetic structure of populations from multilocus genotype data has become a central component of modern population‐genetic data analysis. Application of model‐based clustering programs often entails a number of steps, in which the user considers different modelling assumptions, compares results across different predetermined values of the number of assumed clusters (a parameter typically denoted K), examines multiple independent runs for each fixed value of K, and distinguishes among runs belonging to substantially distinct clustering solutions. Here, we present Clumpak (Cluster Markov Packager Across K), a method that automates the postprocessing of results of model‐based population structure analyses. For analysing multiple independent runs at a single K value, Clumpak identifies sets of highly similar runs, separating distinct groups of runs that represent distinct modes in the space of possible solutions. This procedure, which generates a consensus solution for each distinct mode, is performed by the use of a Markov clustering algorithm that relies on a similarity matrix between replicate runs, as computed by the software Clumpp . Next, Clumpak identifies an optimal alignment of inferred clusters across different values of K, extending a similar approach implemented for a fixed K in Clumpp and simplifying the comparison of clustering results across different K values. Clumpak incorporates additional features, such as implementations of methods for choosing K and comparing solutions obtained by different programs, models, or data subsets. Clumpak , available at http://clumpak.tau.ac.il , simplifies the use of model‐based analyses of population structure in population genetics and molecular ecology.     

A PCA-based method for ancestral informative markers selection in structured populations

Identification of population structure can help trace population histories and identify disease genes. Structured association (SA) is a commonly used approach for population structure identification and association mapping. A major issue with SA is that its performance greatly depends on the informativeness and the numbers of ancestral informative markers (AIMs). Present major AIM selection methods mostly require prior individual ancestry information, which is usually not available or uncertain in practice. To address this potential weakness, we herein develop a novel approach for AIM selection based on principle component analysis (PCA), which does not require prior ancestry information of study subjects. Our simulation and real genetic data analysis results suggest that, with equivalent AIMs, PCA-based selected AIMs can significantly increase the accuracy of inferred individual ancestries compared with traditionally randomly selected AIMs. Our method can easily be applied to whole genome data to select a set of highly informative AIMs in population structure, which can then be used to identify potential population structure and correct possible statistical biases caused by population stratification.     

Fine-Scale Genetic Structure in the United Arab Emirates Reflects Endogamous and Consanguineous Culture,Population History,and Geography

The indigenous population of the United Arab Emirates (UAE) has a unique demographic and cultural history. Its tradition of endogamy and consanguinity is expected to produce genetic homogeneity and partitioning of gene pools while population movements and intercontinental trade are likely to have contributed to genetic diversity. Emiratis and neighboring populations of the Middle East have been underrepresented in the population genetics literature with few studies covering the broader genetic history of the Arabian Peninsula. Here, we genotyped 1,198 individuals from the seven Emirates using 1.7 million markers and by employing haplotype-based algorithms and admixture analyses, we reveal the fine-scale genetic structure of the Emirati population. Shared ancestry and gene flow with neighboring populations display their unique geographic position while increased intra- versus inter-Emirati kinship and sharing of uniparental haplogroups, reflect the endogamous and consanguineous cultural traditions of the Emirates and their tribes.     

随机SNP在全基因组关联研究人群分层分析中的应用

在复杂疾病的全基因组关联研究中,人群分层现象会增加结果的假阳性率,因此考虑人群遗传结构、控制人群分层是很有必要的。而在人群分层研究中,使用随机选择的SNP的效果还有待进一步探讨。文章利用HapMap Phase2人群中无关个体的Affymetrix SNP 6.0芯片分型数据,在全基因组上随机均匀选择不同数量的SNP,同时利用f值和Fisher精确检验方法筛选祖先信息标记（Ancestry Informative Markers,AIMs）。然后利用HapMap Phase3中的无关个体的数据,以F-statistics和STRUCTURE分析两种方法评估所选出的不同SNP组合对人群的区分效果。研究发现,随机均匀分布于全基因组的SNP可用于识别人群内部存在的遗传结构。文章进一步提示,在全基因组关联研究中,当没有针对特定人群的AIMs时,可在全基因组上随机选择3000以上均匀分布的SNP来控制人群分层。     

Population genetics models of local ancestry

Migrations have played an important role in shaping the genetic diversity of human populations. Understanding genomic data thus requires careful modeling of historical gene flow. Here we consider the effect of relatively recent population structure and gene flow and interpret genomes of individuals that have ancestry from multiple source populations as mosaics of segments originating from each population. This article describes general and tractable models for local ancestry patterns with a focus on the length distribution of continuous ancestry tracts and the variance in total ancestry proportions among individuals. The models offer improved agreement with Wright-Fisher simulation data when compared to the state-of-the art and can be used to infer time-dependent migration rates from multiple populations. Considering HapMap African-American (ASW) data, we find that a model with two distinct phases of "European" gene flow significantly improves the modeling of both tract lengths and ancestry variances.     

distruct: a program for the graphical display of population structure

In analysis of multilocus genotypes from structured populations, individual coefficients of membership in subpopulations are often estimated using programs such as structure . distruct provides a general method for visualizing these estimated membership coefficients. Subpopulations are represented as colours, and individuals are depicted as bars partitioned into coloured segments that correspond to membership coefficients in the subgroups. distruct , available at http://www.cmb.usc.edu/~noahr/distruct.html , can also be used to display subpopulation assignment probabilities when individuals are assumed to have ancestry in only one group.     

Evaluation of model fit of inferred admixture proportions

Model based methods for genetic clustering of individuals, such as those implemented in structure or ADMIXTURE, allow the user to infer individual ancestries and study population structure. The underlying model makes several assumptions about the demographic history that shaped the analysed genetic data. One assumption is that all individuals are a result of K homogeneous ancestral populations that are all well represented in the data, while another assumption is that no drift happened after the admixture event. The histories of many real world populations do not conform to that model, and in that case taking the inferred admixture proportions at face value might be misleading. We propose a method to evaluate the fit of admixture models based on estimating the correlation of the residual difference between the true genotypes and the genotypes predicted by the model. When the model assumptions are not violated, the residuals from a pair of individuals are not correlated. In the case of a bad fitting admixture model, individuals with similar demographic histories have a positive correlation of their residuals. Using simulated and real data, we show how the method is able to detect a bad fit of inferred admixture proportions due to using an insufficient number of clusters K or to demographic histories that deviate significantly from the admixture model assumptions, such as admixture from ghost populations, drift after admixture events and nondiscrete ancestral populations. We have implemented the method as an open source software that can be applied to both unphased genotypes and low depth sequencing data.     

Impact of cross-ancestry genetic architecture on GWASs in admixed populations

1. Download : Download high-res image (231KB)
2. Download : Download full-size image

     

The computer program structure for assigning individuals to populations: easy to use but easier to misuse

The computer program Structure implements a Bayesian method, based on a population genetics model, to assign individuals to their source populations using genetic marker data. It is widely applied in the fields of ecology, evolutionary biology, human genetics and conservation biology for detecting hidden genetic structures, inferring the most likely number of populations (K), assigning individuals to source populations and estimating admixture and migration rates. Recently, several simulation studies repeatedly concluded that the program yields erroneous inferences when samples from different populations are highly unbalanced in size. Analysing both simulated and empirical data sets, this study confirms that Structure indeed yields poor individual assignments to source populations and gives frequently incorrect estimates of K when sampling is unbalanced. However, this poor performance is mainly caused by the adoption of the default ancestry prior, which assumes all source populations contribute equally to the pooled sample of individuals. When the alternative ancestry prior, which allows for unequal representations of the source populations by the sample, is adopted, accurate individual assignments could be obtained even if sampling is highly unbalanced. The alternative prior also improves the inference of K by two estimators, albeit the improvement is not as much as that in individual assignments to populations. For the difficult case of many populations and unbalanced sampling, a rarely used parameter combination of the alternative ancestry prior, an initial ALPHA value much smaller than the default and the uncorrelated allele frequency model is required for Structure to yield accurate inferences. I conclude that Structure is easy to use but is easier to misuse because of its complicated genetic model and many parameter (prior) options which may not be obvious to choose, and suggest using multiple plausible models (parameters) and K estimators in conducting comparative and exploratory Structure analysis.     

Population structure and genetic diversity in tristylous Narcissus triandrus: insights from microsatellite and chloroplast DNA variation

We investigated cpDNA sequence and nuclear microsatellite variation among populations of the wild daffodil Narcissus triandrus to examine the role of historical vs. contemporary forces in shaping population structure, morphological differentiation and sexual-system evolution. This wide-ranging heterostylous species of the Iberian Peninsula is largely composed of two allopatric varieties (vars. cernuus and triandrus), and populations with either stylar trimorphism or dimorphism. Dimorphic populations only occur in var. triandrus, are mainly restricted to the northwestern portion of the species range, and uniformly lack the mid-styled morph (M-morph). Chloroplast DNA (cpDNA) sequence variation revealed strong geographical structuring and evidence for a fragmentation event associated with differentiation of the two varieties. In var. triandrus, population fragmentation, restricted gene flow and isolation-by-distance were also inferred. Significant differences in genetic diversity and population structure between the two varieties likely reflect historical and contemporary differences in demography and gene flow among populations. Discordance between cpDNA markers and both microsatellite and morphological variation indicate that hybridization has occurred between the two varieties at contact zones. There were no differences in genetic diversity or population structure between dimorphic and trimorphic populations, and chloroplast haplotypes were not associated with either sexual system, indicating transitions in morph structure within each maternal lineage. M-morph frequencies were positively correlated with differentiation at microsatellite loci, indicating that the evolutionary processes influencing these neutral markers also influence alleles controlling the style morphs.     

Indirect evidence for the genetic determination of short stature in African Pygmies

Central African Pygmy populations are known to be the shortest human populations worldwide. Many evolutionary hypotheses have been proposed to explain this short stature: adaptation to food limitations, climate, forest density, or high mortality rates. However, such hypotheses are difficult to test given the lack of long-term surveys and demographic data. Whether the short stature observed nowadays in African Pygmy populations as compared to their Non-Pygmy neighbors is determined by genetic factors remains widely unknown. Here, we study a uniquely large new anthropometrical dataset comprising more than 1,000 individuals from 10 Central African Pygmy and neighboring Non-Pygmy populations, categorized as such based on cultural criteria rather than height. We show that climate, or forest density may not play a major role in the difference in adult stature between existing Pygmies and Non-Pygmies, without ruling out the hypothesis that such factors played an important evolutionary role in the past. Furthermore, we analyzed the relationship between stature and neutral genetic variation in a subset of 213 individuals and found that the Pygmy individuals' stature was significantly positively correlated with levels of genetic similarity with the Non-Pygmy gene-pool for both men and women. Overall, we show that a Pygmy individual exhibiting a high level of genetic admixture with the neighboring Non-Pygmies is likely to be taller. These results show for the first time that the major morphological difference in stature found between Central African Pygmy and Non-Pygmy populations is likely determined by genetic factors.     

Microsatellite analysis of genetic diversity and population structure of Chinese mitten crab （Eriocheir sinensis）

Chinese mitten crab (Eriocheir sinensis) has higher commercial value as food source than any other species of Eriocheir in China.To evaluate the germplasm resources and characterize the genetic diversity and population structure of the crabs in different water systems,two stocks and two farming populations were assessed with 25 polymorphic microsallite loci available in public GenBank.Basic statistics showed that the average observed heterozygosity (Ho) amongst populations ranged from 0.5789 to 0.6824.However,a remarkable presence of inbreeding and heterozygote deficiencies were observed.To analyze population structure,pairwise FST coefficients explained only ～10.3% variability from the subdivision of mitten crab populations,the remaining variability stems from the subdivision within subpopulations.Although the four populations had slight differentiation,different allelic frequencies resulted in distinct population structures.Two stocks and one farming population were clustered together to the phylogenetic branch of Yangtze crab,with an approximate membership of 95%.Whereas,another fanning population was clustered singly to the phylogenetic branch of the Liaohe crab,with a membership of 97.1%.The tests for individual admixture showed that Yangtze crab had probably been contaminated with individuals from other water systems.Genetic relationships between populations also supported the conclusion that Yangtze crab and Liaohe crab had different gene pools in spite of the origins of the same species.