Principal components analysis of population admixture

With the availability of high-density genotype information, principal components analysis (PCA) is now routinely used to detect and quantify the genetic structure of populations in both population genetics and genetic epidemiology. An important issue is how to make appropriate and correct inferences about population relationships from the results of PCA, especially when admixed individuals are included in the analysis. We extend our recently developed theoretical formulation of PCA to allow for admixed populations. Because the sampled individuals are treated as features, our generalized formulation of PCA directly relates the pattern of the scatter plot of the top eigenvectors to the admixture proportions and parameters reflecting the population relationships, and thus can provide valuable guidance on how to properly interpret the results of PCA in practice. Using our formulation, we theoretically justify the diagnostic of two-way admixture. More importantly, our theoretical investigations based on the proposed formulation yield a diagnostic of multi-way admixture. For instance, we found that admixed individuals with three parental populations are distributed inside the triangle formed by their parental populations and divide the triangle into three smaller triangles whose areas have the same proportions in the big triangle as the corresponding admixture proportions. We tested and illustrated these findings using simulated data and data from HapMap III and the Human Genome Diversity Project.     

Ethnicity determination by names among the Aymara of Chile and Bolivia

The importance of surnames in genetic studies has been recognized for a century or so. While the ethnic affiliations of individuals are ordinarily established in genetic studies by admixture analysis based on gene frequencies, often there are implicit assumptions in these attempts that are difficult to validate in the absence of detailed ethnohistories. In northern Chile and western Bolivia, where genetic admixture has been known to occur among the Aymara Indians and Spanish Caucasoids, the naming pattern (parental patriand matrinyms) allowed us to classify individuals on the basis of the frequency of Aymara names into 9 'ethnic' groups. From a sample of 2525 individuals it is shown that admixture occurred in lineages nonrandomly, implying assortative mating of surnames. Admixture and genetic distance analysis on the basis of 31 genetic markers on approximately 1700 of these individuals reveals that there is a reasonable agreement of ethnic classification of individuals by name and phenotype data on genetic markers. The Aymara-named groups are shown to be predominantly Amerindian (89%) in their genetic profiles. Individuals whose current naming pattern is basically Spanish also exhibit a substantial fraction of genes of Amerindian origin (67%). Presence of some rare alleles not found in Amerindian or Spanish Caucasoids in the admixed groups suggest infiltration of Negroid genes in the past.     

The Developmental Basis of Quantitative Craniofacial Variation in Humans and Mice

The human skull is a complex and highly integrated structure that has long held the fascination of anthropologists and evolutionary biologists. Recent studies of the genetics of craniofacial variation reveal a very complex and multifactorial picture. These findings contrast with older ideas that posit much simpler developmental bases for variation in cranial morphology such as the growth of the brain or the growth of the chondrocranium relative to the dermatocranium. Such processes have been shown to have major effects on cranial morphology in mice. It is not known, however, whether they are relevant to explaining normal phenotypic variation in humans. To answer this question, we obtained vectors of shape change from mutant mouse models in which the developmental basis for the craniofacial phenotype is known to varying degrees, and compared these to a homologous dataset constructed from human crania obtained from a single population with a known genealogy. Our results show that the shape vectors associated with perturbations to chondrocranial growth, brain growth, and body size in mice do largely correspond to axes of covariation in humans. This finding supports the view that the developmental basis for craniofacial variation funnels down to a relatively small number of key developmental processes that are similar across mice and humans. Understanding these processes and how they influence craniofacial shape provides fundamental insights into the developmental basis for evolutionary change in the human skull as well as the developmental-genetic basis for normal phenotypic variation in craniofacial form.     

The Admixture Structure and Genetic Variation of the Archipelago of Cape Verde and Its Implications for Admixture Mapping Studies

Recently admixed populations offer unique opportunities for studying human history and for elucidating the genetic basis of complex traits that differ in prevalence between human populations. Historical records, classical protein markers, and preliminary genetic data indicate that the Cape Verde islands in West Africa are highly admixed and primarily descended from European males and African females. However, little is known about the variation in admixture levels, admixture dynamics and genetic diversity across the islands, or about the potential of Cape Verde for admixture mapping studies. We have performed a detailed analysis of phenotypic and genetic variation in Cape Verde based on objective skin color measurements, socio-economic status (SES) evaluations and data for 50 autosomal, 34 X-chromosome, and 21 non-recombinant Y-chromosome (NRY) markers in 845 individuals from six islands of the archipelago. We find extensive genetic admixture between European and African ancestral populations (mean West African ancestry = 0.57, sd = 0.08), with individual African ancestry proportions varying considerably among the islands. African ancestry proportions calculated with X and Y-chromosome markers confirm that the pattern of admixture has been sex-biased. The high-resolution NRY-STRs reveal additional patterns of variation among the islands that are most consistent with differentiation after admixture. The differences in the autosomal admixture proportions are clearly evident in the skin color distribution across the islands (Pearson r = 0.54, P-value<2e–16). Despite this strong correlation, there are significant interactions between SES and skin color that are independent of the relationship between skin color and genetic ancestry. The observed distributions of admixture, genetic variation and skin color and the relationship of skin color with SES relate to historical and social events taking place during the settlement history of Cape Verde, and have implications for the design of association studies using this population.     

Human genetic admixture

Throughout human history, large-scale migrations have facilitated the formation of populations with ancestry from multiple previously separated populations. This process leads to subsequent shuffling of genetic ancestry through recombination, producing variation in ancestry between populations, among individuals in a population, and along the genome within an individual. Recent methodological and empirical developments have elucidated the genomic signatures of this admixture process, bringing previously understudied admixed populations to the forefront of population and medical genetics. Under this theme, we present a collection of recent PLOS Genetics publications that exemplify recent progress in human genetic admixture studies, and we discuss potential areas for future work.     

A genomewide admixture mapping study for yield factors and morphological traits in a cultivated cocoa (Theobroma cacao L.) population

The selection of productive varieties of modern Criollo cocoa, showing fine aromatic qualities in their beans, is of major interest for some producing countries, such as Venezuela. Cultivated populations of Modern Criollo or Trinitario varieties may be suitable for admixture mapping analysis, as large blocks of alleles derived from two identified divergent ancestors, recently admixed, are still preserved, after a few generations of recombination, similar to experimental mapping progenies. Two hundred and fifty-seven individuals from a cultivated population of Modern Criollo were selected and analysed with 92 microsatellite markers distributed along the genome. This population exhibited a wide range of variability for yield factors and morphological features. Population structure analysis identified two main subgroups corresponding to the admixture from the two ancestors Criollo and Forastero. Several significant associations between markers and phenotypic data (yield factors and morphological traits) were identified by a least squares general linear model (GLM) taking into account the population structure and the percentage of admixture of each individual. Results were compared with classical QTL analyses previously reported for other cacao populations. Most markers associated to quantitative traits were very close to QTLs detected formerly for the same traits. Associations were also identified between markers and several qualitative traits including the red pigmentation observed in different organs, mainly associated to common markers in linkage group 4.     

Detecting directional selection in the presence of recent admixture in African-Americans

We investigate the performance of tests of neutrality in admixed populations using plausible demographic models for African-American history as well as resequencing data from African and African-American populations. The analysis of both simulated and human resequencing data suggests that recent admixture does not result in an excess of false-positive results for neutrality tests based on the frequency spectrum after accounting for the population growth in the parental African population. Furthermore, when simulating positive selection, Tajima's D, Fu and Li's D, and haplotype homozygosity have lower power to detect population-specific selection using individuals sampled from the admixed population than from the nonadmixed population. Fay and Wu's H test, however, has more power to detect selection using individuals from the admixed population than from the nonadmixed population, especially when the selective sweep ended long ago. Our results have implications for interpreting recent genome-wide scans for positive selection in human populations.     

Dermatoglyphic changes during the population admixture between Kam and Han Chinese

Genetic studies and gene localization for human dermatoglyphs are currently ongoing. However, the inheritance modes of various genetic traits are not well understood because of the complexity of dermatoglyph genetics. The study of admixed populations can contribute to the understanding of population genetic traits of dermatoglyphs. Here, we present the dermatoglyphic characteristics of Kam and Liujia Han, and the admixed population consisting of these two parent populations.The characteristics of the admixed population do not always fall in the same ranges as the parent population characters but do seem to be biased to Kam or Liujia parent populations, depending on sex and ethnicity of parents. The total frequencies of different fingerprint types do not differ among these populations, but several of the quantitative traits and the palm true pattern frequencies do significantly differ between admixed and parent populations. The simple arch fingerprint frequency decreases significantly in the admixed population in comparison with parent populations while both simple whorl fingerprint frequency and finger ridge counts increase significantly. True pattern frequency of the span area of interdigital III and IV areas on right hands and the radial-loop frequency of the right index fingers in the admixed populations are consistent with their matrilineal population. These dermatoglyphic changes may result from increased heterozygosity in the admixed population. The genetic modes of these changes may be relatively simple and will be useful for future dermatoglyph genetic studies.      

Mathematical properties of Fst between admixed populations and their parental source populations

We consider the properties of the F(st) measure of genetic divergence between an admixed population and its parental source populations. Among all possible populations admixed among an arbitrary set of parental populations, we show that the value of F(st) between an admixed population and a specific source population is maximized when the admixed population is simply the most distant of the other source populations. For the case with only two parental populations, as a function of the admixture fraction, we further demonstrate that this F(st) value is monotonic and convex, so that F(st) is informative about the admixture fraction. We illustrate our results using example human population-genetic data, showing how they provide a framework in which to interpret the features of F(st) in admixed populations.     

Autosomal Admixture Levels Are Informative About Sex Bias in Admixed Populations

Sex-biased admixture has been observed in a wide variety of admixed populations. Genetic variation in sex chromosomes and functions of quantities computed from sex chromosomes and autosomes have often been examined to infer patterns of sex-biased admixture, typically using statistical approaches that do not mechanistically model the complexity of a sex-specific history of admixture. Here, expanding on a model of Verdu and Rosenberg (2011) that did not include sex specificity, we develop a model that mechanistically examines sex-specific admixture histories. Under the model, multiple source populations contribute to an admixed population, potentially with their male and female contributions varying over time. In an admixed population descended from two source groups, we derive the moments of the distribution of the autosomal admixture fraction from a specific source population as a function of sex-specific introgression parameters and time. Considering admixture processes that are constant in time, we demonstrate that surprisingly, although the mean autosomal admixture fraction from a specific source population does not reveal a sex bias in the admixture history, the variance of autosomal admixture is informative about sex bias. Specifically, the long-term variance decreases as the sex bias from a contributing source population increases. This result can be viewed as analogous to the reduction in effective population size for populations with an unequal number of breeding males and females. Our approach suggests that it may be possible to use the effect of sex-biased admixture on autosomal DNA to assist with methods for inference of the history of complex sex-biased admixture processes.     

History shaped the geographic distribution of genomic admixture on the island of Puerto Rico

Contemporary genetic variation among Latin Americans human groups reflects population migrations shaped by complex historical, social and economic factors. Consequently, admixture patterns may vary by geographic regions ranging from countries to neighborhoods. We examined the geographic variation of admixture across the island of Puerto Rico and the degree to which it could be explained by historic and social events. We analyzed a census-based sample of 642 Puerto Rican individuals that were genotyped for 93 ancestry informative markers (AIMs) to estimate African, European and Native American ancestry. Socioeconomic status (SES) data and geographic location were obtained for each individual. There was significant geographic variation of ancestry across the island. In particular, African ancestry demonstrated a decreasing East to West gradient that was partially explained by historical factors linked to the colonial sugar plantation system. SES also demonstrated a parallel decreasing cline from East to West. However, at a local level, SES and African ancestry were negatively correlated. European ancestry was strongly negatively correlated with African ancestry and therefore showed patterns complementary to African ancestry. By contrast, Native American ancestry showed little variation across the island and across individuals and appears to have played little social role historically. The observed geographic distributions of SES and genetic variation relate to historical social events and mating patterns, and have substantial implications for the design of studies in the recently admixed Puerto Rican population. More generally, our results demonstrate the importance of incorporating social and geographic data with genetics when studying contemporary admixed populations.     

Genetic admixture, self-reported ethnicity, self-estimated admixture, and skin pigmentation among Hispanics and Native Americans

The relationship between ethnicity and biology is of interest to anthropologists, biomedical scientists, and historians in understanding how human groups are constructed. Ethnic self-identification in recently admixed groups such as Hispanics, African Americans, and Native Americans (NA) is likely to be complex due to the heterogeneity in individual admixture proportions and social environments within these groups. This study examines the relationships between self-identified ethnicity, self-estimated admixture proportions, skin pigmentation, and genetic marker estimated admixture proportions. These measures were assessed using questionnaires, skin color measurements, and genotyping of a panel of 76 ancestry informative markers, among 170 Hispanics and NAs from New Mexico, a state known for its complex history of interactions between people of NA and European (EU) ancestry. Results reveal that NAs underestimate their degree of EU admixture, and that Hispanics underestimate their degree of NA admixture. Within Hispanics, genetic-marker estimated admixture is better predicted by forehead skin pigmentation than by self-estimated admixture. We also find that Hispanic individuals self-identified as "half-White, half Hispanic" and "Spanish" have lower levels of NA admixture than those self-identified as "Mexican" and "Mexican American." Such results highlight the interplay between culture and biology in how individuals identify and view themselves, and have implications for how ethnicity and disease risk are assessed in a medical setting.     

A general population-genetic model for the production by population structure of spurious genotype-phenotype associations in discrete, admixed or spatially distributed populations

In linkage disequilibrium mapping of genetic variants causally associated with phenotypes, spurious associations can potentially be generated by any of a variety of types of population structure. However, mathematical theory of the production of spurious associations has largely been restricted to population structure models that involve the sampling of individuals from a collection of discrete subpopulations. Here, we introduce a general model of spurious association in structured populations, appropriate whether the population structure involves discrete groups, admixture among such groups, or continuous variation across space. Under the assumptions of the model, we find that a single common principle--applicable to both the discrete and admixed settings as well as to spatial populations--gives a necessary and sufficient condition for the occurrence of spurious associations. Using a mathematical connection between the discrete and admixed cases, we show that in admixed populations, spurious associations are less severe than in corresponding mixtures of discrete subpopulations, especially when the variance of admixture across individuals is small. This observation, together with the results of simulations that examine the relative influences of various model parameters, has important implications for the design and analysis of genetic association studies in structured populations.     

Disentangling Signatures of Selection Before and After European Colonization in Latin Americans

Throughout human evolutionary history, large-scale migrations have led to intermixing (i.e., admixture) between previously separated human groups. Although classical and recent work have shown that studying admixture can yield novel historical insights, the extent to which this process contributed to adaptation remains underexplored. Here, we introduce a novel statistical model, specific to admixed populations, that identifies loci under selection while determining whether the selection likely occurred post-admixture or prior to admixture in one of the ancestral source populations. Through extensive simulations, we show that this method is able to detect selection, even in recently formed admixed populations, and to accurately differentiate between selection occurring in the ancestral or admixed population. We apply this method to genome-wide SNP data of ∼4,000 individuals in five admixed Latin American cohorts from Brazil, Chile, Colombia, Mexico, and Peru. Our approach replicates previous reports of selection in the human leukocyte antigen region that are consistent with selection post-admixture. We also report novel signals of selection in genomic regions spanning 47 genes, reinforcing many of these signals with an alternative, commonly used local-ancestry-inference approach. These signals include several genes involved in immunity, which may reflect responses to endemic pathogens of the Americas and to the challenge of infectious disease brought by European contact. In addition, some of the strongest signals inferred to be under selection in the Native American ancestral groups of modern Latin Americans overlap with genes implicated in energy metabolism phenotypes, plausibly reflecting adaptations to novel dietary sources available in the Americas.     

Maximum-likelihood estimation of admixture proportions from genetic data

For an admixed population, an important question is how much genetic contribution comes from each parental population. Several methods have been developed to estimate such admixture proportions, using data on genetic markers sampled from parental and admixed populations. In this study, I propose a likelihood method to estimate jointly the admixture proportions, the genetic drift that occurred to the admixed population and each parental population during the period between the hybridization and sampling events, and the genetic drift in each ancestral population within the interval between their split and hybridization. The results from extensive simulations using various combinations of relevant parameter values show that in general much more accurate and precise estimates of admixture proportions are obtained from the likelihood method than from previous methods. The likelihood method also yields reasonable estimates of genetic drift that occurred to each population, which translate into relative effective sizes (N(e)) or absolute average N(e)'s if the times when the relevant events (such as population split, admixture, and sampling) occurred are known. The proposed likelihood method also has features such as relatively low computational requirement compared with previous ones, flexibility for admixture models, and marker types. In particular, it allows for missing data from a contributing parental population. The method is applied to a human data set and a wolflike canids data set, and the results obtained are discussed in comparison with those from other estimators and from previous studies.     

The influence of stochastic and selective forces in the population divergence of female colour polymorphism in damselflies of the genus Ischnura

Disentangling the relative importance and potential interactions of selection and genetic drift in driving phenotypic divergence of species is a classical research topic in population genetics and evolutionary biology. Here, we evaluate the role of stochastic and selective forces on population divergence of a colour polymorphism in seven damselfly species of the genus Ischnura, with a particular focus on I. elegans and I. graellsii. Colour-morph frequencies in Spanish I. elegans populations varied greatly, even at a local scale, whereas more similar frequencies were found among populations in eastern Europe. In contrast, I. graellsii and the other five Ischnura species showed little variation in colour-morph frequencies between populations. F(ST)-outlier analyses revealed that the colour locus deviated strongly from neutral expectations in Spanish populations of I. elegans, contrasting the pattern found in eastern European populations, and in I. graellsii, where no such discrepancy between morph divergence and neutral divergence could be detected. This suggests that divergent selection has been operating on the colour locus in Spanish populations of I. elegans, whereas processes such as genetic drift, possibly in combination with other forms of selection (such as negative frequency-dependent selection), appear to have been present in other regions, such as eastern Europe. Overall, the results indicate that both selective and stochastic processes operate on these colour polymorphisms, and suggest that the relative importance of factors varies between geographical regions.     

Characterization of mitochondrial haplogroups in a large population-based sample from the United States

Mitochondrial DNA (mtDNA) haplogroups are valuable for investigations in forensic science, molecular anthropology, and human genetics. In this study, we developed a custom panel of 61 mtDNA markers for high-throughput classification of European, African, and Native American/Asian mitochondrial haplogroup lineages. Using these mtDNA markers, we constructed a mitochondrial haplogroup classification tree and classified 18,832 participants from the National Health and Nutrition Examination Surveys (NHANES). To our knowledge, this is the largest study to date characterizing mitochondrial haplogroups in a population-based sample from the United States, and the first study characterizing mitochondrial haplogroup distributions in self-identified Mexican Americans separately from Hispanic Americans of other descent. We observed clear differences in the distribution of maternal genetic ancestry consistent with proposed admixture models for these subpopulations, underscoring the genetic heterogeneity of the United States Hispanic population. The mitochondrial haplogroup distributions in the other self-identified racial/ethnic groups within NHANES were largely comparable to previous studies. Mitochondrial haplogroup classification was highly concordant with self-identified race/ethnicity (SIRE) in non-Hispanic whites (94.8 %), but was considerably lower in admixed populations including non-Hispanic blacks (88.3 %), Mexican Americans (81.8 %), and other Hispanics (61.6 %), suggesting SIRE does not accurately reflect maternal genetic ancestry, particularly in populations with greater proportions of admixture. Thus, it is important to consider inconsistencies between SIRE and genetic ancestry when performing genetic association studies. The mitochondrial haplogroup data that we have generated, coupled with the epidemiologic variables in NHANES, is a valuable resource for future studies investigating the contribution of mtDNA variation to human health and disease.     

Identification of interspecific hybrids among domesticated apple and its wild relatives

Potential interspecific hybrids are usually identified in natural populations by their proximity to interbreeding species or their intermediate phenotypes; hybridization can then be confirmed by comparing the genetic makeup of putative hybrids to pure species. In contrast, detecting interspecific hybridization and misclassifications in ex situ collections can be difficult because fine-scale geographic locations and species-specific phenotypic data are generally unavailable. Thus, there is little a priori information available to suggest which individuals might be hybrids. Instead, hybrids or misclassified individuals must be identified based on molecular data via population assignment and admixture detection programs. We have applied a variety of population assignment and admixture detection programs to over 400 samples of four closely related Malus species held in the US Department of Agriculture?CAgricultural Research Service National Plant Germplasm System that were genotyped at 19 simple sequence repeat loci. Our findings indicate that over 10?% of the samples of the wild species Malus sieversii and Malus orientalis and nearly 20?% of the samples of the wild species Malus sylvestris may be admixed or misclassified. The percentage of admixed or misclassified samples of the domesticated species, Malus × domestica, was much lower, at <5?%. These findings provide an illustration of how to detect hybridization and misclassification in ex situ collections using molecular data and, ultimately, should help to maximize the utility of the collections. In addition, the presence of wild-collected samples that show admixture with domesticated apple suggests that gene flow may be occurring from the crop into natural populations of the wild species.     

Molecular control of facial morphology

We present a developmental perspective on the concept of phylotypic and phenotypic stages of craniofacial development. Within orders of avians and mammals, a phylotypic period exists when the morphology of the facial prominences is minimally divergent. We postulate that species-specific facial variations arise as a result of subtle shifts in the timing and the duration of molecular pathway activity (e.g., heterochrony), and present evidence demonstrating a critical role for Wnt and FGF signaling in this process. The same molecular pathways that shape the vertebrate face are also implicated in craniofacial deformities, indicating that comparisons between and among animal species may represent a novel method for the identification of human craniofacial disease genes.     

Unequal contributions of male and female gene pools from parental populations in the African descendants of the city of Melo, Uruguay

In admixed populations, genetic contributions from males and females of specific parental populations can be of different proportions due to past directional mating during the process of genetic admixture. In this research paper, we provide evidence of such male- and female-specific differential admixture components of African, European, and American Indian origin in an admixed population from the city of Melo, in the northeastern region of Uruguay. From data on 11 autosomal markers from a sample of 41 individuals of mixed African descent, we estimated 47% African, 38% European, and 15% Amerindian contributions. In contrast, 6 mtDNA site-specific polymorphic markers showed that the mtDNA genome of these individuals was 52% African, 19% European, and 29% Amerindian, while from 3 Y-specific polymorphic sites, we estimated 30% African, 64% European, and 6% Amerindian contributions. We argue that this heterogeneity of admixture estimates results from disproportionate unions of European males with African and American Indian females from which this mixed African population was formed. Also, we argue that the asymmetry of the admixture estimates from the three sets of markers (autosomal, mtDNA, and Y-linked) is a result of the changes in the direction of mating during the history of the population. Implications of such evidence of directional mating are discussed, indicating the need of further demographic data for a quantitative assessment of the impact of directional mating on genetic structure of admixed populations.