A general mechanistic model for admixture histories of hybrid populations

Admixed populations have been used for inferring migrations, detecting natural selection, and finding disease genes. These applications often use a simple statistical model of admixture rather than a modeling perspective that incorporates a more realistic history of the admixture process. Here, we develop a general model of admixture that mechanistically accounts for complex historical admixture processes. We consider two source populations contributing to the ancestry of a hybrid population, potentially with variable contributions across generations. For a random individual in the hybrid population at a given point in time, we study the fraction of genetic admixture originating from a specific one of the source populations by computing its moments as functions of time and of introgression parameters. We show that very different admixture processes can produce identical mean admixture proportions, but that such processes produce different values for the variance of the admixture proportion. When introgression parameters from each source population are constant over time, the long-term limit of the expectation of the admixture proportion depends only on the ratio of the introgression parameters. The variance of admixture decreases quickly over time after the source populations stop contributing to the hybrid population, but remains substantial when the contributions are ongoing. Our approach will facilitate the understanding of admixture mechanisms, illustrating how the moments of the distribution of admixture proportions can be informative about the historical admixture processes contributing to the genetic diversity of hybrid populations.     

The Tamang: a case of hybridization.

The Tamang, a Mongoloid and hybrid population of the eastern part of the cis-Himalayan region of West Bengal, have been studied for an estimation of the Tibetan and Nepalese components. The degree of admixture of the two parental stocks is found to be about 59 and 41%, respectively. The deviation between the observed and expected gene frequencies of the hybrid is insignificant.     

Cherokee admixture and its estimation by the gene identity method: A critique

A previous application to a Cherokee hybrid population of two methods for measuring admixture proportions by use of gene identity probabilities was examined. The results suggest that neither procedure has been correctly applied, and that underlying statistical difficulties limit their reliability.     

The effects of admixture and population subdivision on cytonuclear disequilibria

We examine the generation of cytonuclear disequilibria by admixture and continued gene flow. General formulas analogous to the nuclear case are first derived showing that the allelic and genotypic disequilibria from admixture or population subdivision equal their expected value across the contributing (sub) populations plus the covariance across these sources between the cytoplasmic gene frequency and the relevant nuclear frequency. A detailed study is then presented of the cytonuclear dynamics, in a random-mating population under two different migration scenarios. In both cases closed-form solutions are given for all variables as a function of the initial conditions and relevant migration parameters. The dynamics of the gene frequencies and allelic disequilibria, which dominate each system, are the same as those involving two unlinked nuclear loci, while the dynamics of the genotypic disequilibria and cytonuclear frequencies have no nuclear counterpart. The continent-island formulation focuses on a population receiving continued immigration from a large source of constant composition. A major discovery is that cytonuclear disequilibria can transiently build up on the "island" to levels far exceeding those found at equilibrium. In contrast, the admixture formulation focuses on the dynamics within two populations undergoing continued intermigration. Although in this case all cytonuclear associations must ultimately decay to zero, long-term transient disequilibria can develop which are many times their initial admixture values. For both migration scenarios it is shown that the time of population censusing relative to migration and reproduction dramatically affects both the amount and pattern of the nonrandom associations produced. The empirical relevance of these models is discussed in light of nuclear-mitochondrial data from a hybrid zone between European and North American eels and from a zone of racial admixture in humans.     

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.     

Detecting signatures of inter‐regional and inter‐specific hybridization among the Chinese rhesus macaque specific pathogen‐free (SPF) population using single nucleotide polymorphic (SNP) markers

Background While rates of gene flow between rhesus and longtail macaque populations near their hybrid zone in Indochina have been quantified elsewhere, this study demonstrates that the inter‐specific introgression is not limited to the Indochinese hybrid zone but is more geographically widespread. Methods Twelve rhesus and longtail macaque populations were analyzed using single nucleotide polymorphic (SNP) loci. Results There is evidence for inter‐specific admixture between Chinese rhesus and mainland longtails, with implications for genetic diversity both in the Chinese super‐SPF population at the California National Primate Research Center and in other primate facilities. Eastern Chinese rhesus appeared more highly derived than western Chinese rhesus, and allele sharing between longtails and Chinese rhesus was not random with regard to geographic distance, but no significant nuclear genetic differences between eastern and western Chinese rhesus were detected among the 245 genic SNPs assayed. Conclusion The implications of this inter‐specific admixture for the use of Chinese rhesus and mainland longtail in biomedical research should be considered.     

Overcoming barriers to active interventions for genetic diversity

As a result of processes such as habitat loss and overharvest, many species persist in small, isolated populations that experience reduced fitness, decreased evolutionary potential, and increased extinction risk. The goal of species conservation is to restore genetic diversity and adaptive potential caused by isolation and small population size. For populations trapped in an extinction vortex, habitat protection may be inadequate for successful conservation. Alternative actions such as deliberate admixture by introducing individuals from related subspecies may be necessary to recover population fitness. While there is precedent for such actions, admixture temporarily disrupts the taxonomic integrity of a species. Concerns for the taxonomic integrity or “naturalness” of a species may prevent the use of active interventions that involve admixture and transient hybrid gene pools even though extinction may be imminent. We explore the cultural barriers to using tools such as genetic rescue and make suggestions for overcoming those barriers. We focus mainly on examples from animals, but the same evolutionary processes are ongoing in other life forms and are subject to the same cultural barriers.     

Exploring Population Admixture Dynamics via Empirical and Simulated Genome-wide Distribution of Ancestral Chromosomal Segments

The processes of genetic admixture determine the haplotype structure and linkage disequilibrium patterns of the admixed population, which is important for medical and evolutionary studies. However, most previous studies do not consider the inherent complexity of admixture processes. Here we proposed two approaches to explore population admixture dynamics, and we demonstrated, by analyzing genome-wide empirical and simulated data, that the approach based on the distribution of chromosomal segments of distinct ancestry (CSDAs) was more powerful than that based on the distribution of individual ancestry proportions. Analysis of 1,890 African Americans showed that a continuous gene flow model, in which the African American population continuously received gene flow from European populations over about 14 generations, best explained the admixture dynamics of African Americans among several putative models. Interestingly, we observed that some African Americans had much more European ancestry than the simulated samples, indicating substructures of local ancestries in African Americans that could have been caused by individuals from some particular lineages having repeatedly admixed with people of European ancestry. In contrast, the admixture dynamics of Mexicans could be explained by a gradual admixture model in which the Mexican population continuously received gene flow from both European and Amerindian populations over about 24 generations. Our results also indicated that recent gene flows from Sub-Saharan Africans have contributed to the gene pool of Middle Eastern populations such as Mozabite, Bedouin, and Palestinian. In summary, this study not only provides approaches to explore population admixture dynamics, but also advances our understanding on population history of African Americans, Mexicans, and Middle Eastern populations.     

The population genomics of multiple tsetse fly (Glossina fuscipes fuscipes) admixture zones in Uganda

Understanding the mechanisms that enforce, maintain or reverse the process of speciation is an important challenge in evolutionary biology. This study investigates the patterns of divergence and discusses the processes that form and maintain divergent lineages of the tsetse fly Glossina fuscipes fuscipes in Uganda. We sampled 251 flies from 18 sites spanning known genetic lineages and the four admixture zones between them. We apply population genomics, hybrid zone and approximate Bayesian computation to the analysis of three types of genetic markers: 55,267 double‐digest restriction site‐associated DNA (ddRAD) SNPs to assess genome‐wide admixture, 16 microsatellites to provide continuity with published data and accurate biogeographic modelling, and a 491‐bp fragment of mitochondrial cytochrome oxidase I and II to infer maternal inheritance patterns. Admixture zones correspond with regions impacted by the reorganization of Uganda's river networks that occurred during the formation of the West African Rift system over the last several hundred thousand years. Because tsetse fly population distributions are defined by rivers, admixture zones likely represent both old and new regions of secondary contact. Our results indicate that older hybrid zones contain mostly parental types, while younger zones contain variable hybrid types resulting from multiple generations of interbreeding. These findings suggest that reproductive barriers are nearly complete in the older admixture zones, while nearly absent in the younger admixture zones. Findings are consistent with predictions of hybrid zone theory: Populations in zones of secondary contact transition rapidly from early to late stages of speciation or collapse all together.     

On the formula for admixture linkage disequilibrium

Admixture linkage disequilibrium (ALD), a phenomenon created by gene flow between genetically distinct populations, has for some time been used as a tool in gene mapping. It is therefore important to analyze the pattern of ALD over generations. In this study we explore two models of admixture: the gradual admixture (GA) model, in which admixture occurs at a variable rate in every generation; and the immediate admixture (IA) model, a special case of the GA model, in which admixture occurs in a single generation. In the case of ALD, the well-known formula of linkage disequilibrium (Delta(t)=(1-r)t Delta(0)) is not applicable under these two models. We note the effect of a random mating population (RMP) to the gametic frequencies from the parental population to the offspring population, and provide the correct formula for ALD.     

Unexpected ancestry of Populus seedlings from a hybrid zone implies a large role for postzygotic selection in the maintenance of species

In the context of potential interspecific gene flow, the integrity of species will be maintained by reproductive barriers that reduce genetic exchange, including traits associated with prezygotic isolation or poor performance of hybrids. Hybrid zones can be used to study the importance of different reproductive barriers, particularly when both parental species and hybrids occur in close spatial proximity. We investigated the importance of barriers to gene flow that act early vs. late in the life cycle of European Populus by quantifying the prevalence of homospecific and hybrid matings within a mosaic hybrid zone. We obtained genotypic data for 11 976 loci from progeny and their maternal parents and constructed a Bayesian model to estimate individual admixture proportions and hybrid classes for sampled trees and for the unsampled pollen parent. Matings that included one or two hybrid parents were common, resulting in admixture proportions of progeny that spanned the whole range of potential ancestries between the two parental species. This result contrasts strongly with the distribution of admixture proportions in adult trees, where intermediate hybrids and each of the parental species are separated into three discrete ancestry clusters. The existence of the full range of hybrids in seedlings is consistent with weak reproductive isolation early in the life cycle of Populus. Instead, a considerable amount of selection must take place between the seedling stage and maturity to remove many hybrid seedlings. Our results highlight that high hybridization rates and appreciable hybrid fitness do not necessarily conflict with the maintenance of species integrity.     

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 相似文献   

Skin colorimetry and admixture measurement: Some further considerations

Several techniques for using skin colorimetric data to measure admixture contributions to a hybrid population have been reviewed by Lees and Relethford ('78), who find this general approach to be useful. This paper extends their evaluation, providing a method for obtaining an error term and considering at greater length certain potential sources of estimate bias. Questions are raised about their discussion of the appropriate employment of distance statistics in this context, and an alternative procedure for obtaining a single estimate of admixture proportions jointly over all reflectance measurements is offered. While it is shown that the hybrid reflectance curves do not support entirely an assumption of unlinked loci contributing equally and additively to skin color, preliminary findings indicate that admixture estimates made from these data are apparently unbiased and reasonably efficient.     

Frequencies of complex diseases in hybrid populations

Diseases of complex etiology demonstrate considerable variation in their frequencies in different ethnic populations. Noninsulin-dependent diabetes mellitus (NIDDM), rheumatoid arthritis, and several cardiovascular diseases constitute examples of such disorders. In genetic studies involving hybrid populations of known ancestry, it is of interest to compare and correlate disease prevalence with the admixture proportion, the latter estimated from a number of polymorphic genetic markers. Theoretical formulations are provided relating disease prevalence in a hybrid population to the admixture proportion under different models of disease transmission. It is shown that the relationship between admixture proportion and disease frequency provides discriminatory power regarding the mode of inheritance. This method is illustrated with an example comparing the proportion of Amerindian ancestry in Mexican-Americans and the prevalence of NIDDM. It is found that genetic factors are involved in susceptibility to NIDDM, but the mode of inheritance cannot be explained by any simple genetic model, and the role of sporadic events cannot be totally ruled out.     

Uncontrolled admixture and loss of genetic diversity in a local Vietnamese pig breed

The expansion of intensive livestock production systems in developing countries has increased the introduction of highly productive exotic breeds facilitating indiscriminate crossbreeding with local breeds. In this study, we set out to investigate the genetic status of the Vietnamese Black H'mong pig breed by evaluating (1) genetic diversity and (2) introgression from exotic breeds. Two exotic breeds, namely Landrace and Yorkshire used for crossbreeding, and the H'mong pig population from Ha Giang (HG) province were investigated using microsatellite markers. Within the province, three phenotypes were observed: a White, a Spotted and a Black phenotype. Genetic differentiation between phenotypes was low (0.5-6.1%). The White phenotypes showed intermediate admixture values between exotic breeds and the Black HG population (0.53), indicating a crossbreed status. Management practices were used to predict the rate of private diversity loss due to exotic gene introgressions. After 60 generations, 100% of Black private alleles will be lost. This loss is accelerated if the admixture rate is increased but can be slowed down if the mortality rate (e.g., recruitment rate) is decreased. Our study showed that a large number of markers are needed for accurately identifying hybrid classes for closely related populations. While our estimate of admixture still seems underestimated, genetic erosion can occur very fast even through indiscriminate crossbreeding.     

Development of two multiplex mini-sequencing panels of ancestry informative SNPs for studies in Latin Americans: an application to populations of the State of Minas Gerais (Brazil)

Admixture occurs when individuals from parental populations that have been isolated for hundreds of generations form a new hybrid population. Currently, interest in measuring biogeographic ancestry has spread from anthropology to forensic sciences, direct-to-consumers personal genomics, and civil rights issues of minorities, and it is critical for genetic epidemiology studies of admixed populations. Markers with highly differentiated frequencies among human populations are informative of ancestry and are called ancestry informative markers (AIMs). For tri-hybrid Latin American populations, ancestry information is required for Africans, Europeans and Native Americans. We developed two multiplex panels of AIMs (for 14 SNPs) to be genotyped by two mini-sequencing reactions, suitable for investigators of medium-small laboratories to estimate admixture of Latin American populations. We tested the performance of these AIMs by comparing results obtained with our 14 AIMs with those obtained using 108 AIMs genotyped in the same individuals, for which DNA samples is available for other investigators. We emphasize that this type of comparison should be made when new admixture/population structure panels are developed. At the population level, our 14 AIMs were useful to estimate European admixture, though they overestimated African admixture and underestimated Native American admixture. Combined with more AIMs, our panel could be used to infer individual admixture. We used our panel to infer the pattern of admixture in two urban populations (Montes Claros and Manhua?u) of the State of Minas Gerais (southeastern Brazil), obtaining a snapshot of their genetic structure in the context of their demographic history.     

Hybridization in contact zone between temperate European pine species

Hybridization studies are important to advance our understanding of the interspecific gene flow and its evolutionary consequences in closely related species. Hybridization and admixture patterns were assessed in a contact zone and reference populations of European pine species using sequence data from 26 nuclear genes and a species-diagnostic cpDNA marker. Reference populations formed three distinct genetic clusters comprising Pinus sylvestris, Pinus mugo/Pinus uliginosa, and Pinus uncinata. Evidence of population structure was found only in P. uliginosa. Based on phenotypic characteristics and molecular data, we identified five groups of individuals in the contact zone in Poland, comprising forms of the parental species and intermediates that were most probably the result of interspecific crosses. A combination of nuclear gene sequence data and a diagnostic organelle marker were used to show that hybridization is frequent in the contact zone and results in hybrid trees with distinct phenotypic identity. The influence of selection in maintaining hybrid phenotypes in environments unsuited to parental species was inferred from nucleotide polymorphism data. A lack of admixture in reference populations suggests that hybridization has not occurred during post-glacial migration and so the contact zone represents a distinct, active example of ongoing evolution. Pine populations in this zone will be a valuable system for studying the genetic basis of hybrid advantage in environmental conditions untypical of pure parental species.     

Inferring Demographic History from a Spectrum of Shared Haplotype Lengths

There has been much recent excitement about the use of genetics to elucidate ancestral history and demography. Whole genome data from humans and other species are revealing complex stories of divergence and admixture that were left undiscovered by previous smaller data sets. A central challenge is to estimate the timing of past admixture and divergence events, for example the time at which Neanderthals exchanged genetic material with humans and the time at which modern humans left Africa. Here, we present a method for using sequence data to jointly estimate the timing and magnitude of past admixture events, along with population divergence times and changes in effective population size. We infer demography from a collection of pairwise sequence alignments by summarizing their length distribution of tracts of identity by state (IBS) and maximizing an analytic composite likelihood derived from a Markovian coalescent approximation. Recent gene flow between populations leaves behind long tracts of identity by descent (IBD), and these tracts give our method power by influencing the distribution of shared IBS tracts. In simulated data, we accurately infer the timing and strength of admixture events, population size changes, and divergence times over a variety of ancient and recent time scales. Using the same technique, we analyze deeply sequenced trio parents from the 1000 Genomes project. The data show evidence of extensive gene flow between Africa and Europe after the time of divergence as well as substructure and gene flow among ancestral hominids. In particular, we infer that recent African-European gene flow and ancient ghost admixture into Europe are both necessary to explain the spectrum of IBS sharing in the trios, rejecting simpler models that contain less population structure.     

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.