Nonneutral admixture of immigrant genotypes in African Drosophila melanogaster populations from Zimbabwe

Drosophila melanogaster originated in Africa and colonized the rest of the world only recently (approximately 10,000 to 15,000 years ago). Using 151 microsatellite loci, we investigated patterns of gene flow between African D. melanogaster populations representing presumptive ancestral variation and recently colonized European populations. Although we detected almost no evidence for alleles of non-African ancestry in a rural D. melanogaster population from Zimbabwe, an urban population from Zimbabwe showed evidence for admixture. Interestingly, the degree of admixture differed among chromosomes. X chromosomes of both rural and urban populations showed almost no non-African ancestry, but the third chromosome in the urban population showed up to 70% of non-African alleles. When chromosomes were broken into contingent microsatellite blocks, even higher estimates of admixture and significant heterogeneity in admixture was observed among these blocks. The discrepancy between the X chromosome and the third chromosome is not consistent with a neutral admixture hypothesis. The higher number of European alleles on the third chromosome could be due to stronger selection against foreign alleles on the X chromosome or to more introgression of (beneficial) alleles on the third chromosome.     

Population genomics as a new tool for wildlife management

Admixture and introgression have varied effects on population viability and fitness. Admixture might be an important source of new alleles, particularly for small, geographically isolated populations. However, admixture might also cause outbreeding depression if populations are adapted to different ecological or climatic conditions. Because of the emerging use of translocation and admixture as a conservation and wildlife management strategy to reduce genetic load (termed genetic rescue), the possible effects of admixture have practical consequences ( Bouzat et al. 2009 ; Hedrick & Fredrickson 2010 ). Importantly, genetic load and local adaptation are properties of individual loci and epistatic interactions among loci rather than properties of genomes. Likewise, the outcome and consequences of genetic rescue depend on the fitness effects of individual introduced alleles. In this issue of Molecular Ecology, Miller et al. (2012) use model‐based, population genomic analyses to document locus‐specific effects of a recent genetic rescue in the bighorn sheep population within the National Bison Range wildlife refuge (NBR; Montana, USA). They find a subset of introduced alleles associated with increased fitness in NBR bighorn sheep, some of which experienced accelerated introgression following their introduction. These loci mark regions of the genome that could constitute the genetic basis of the successful NBR bighorn sheep genetic rescue. Although population genomic analyses are frequently used to study local adaptation and selection (e.g. Hohenlohe et al. 2010 ; Lawniczak et al. 2010 ), this study constitutes a novel application of this analytical framework for wildlife management. Moreover, the detailed demographic data available for the NBR bighorn sheep population provide a rare and powerful source of information and allow more robust population genomic inference than is often possible.     

Surprising Differences in the Variability of Y Chromosomes in African and Cosmopolitan Populations of Drosophila melanogaster

The nonrecombining Drosophila melanogaster Y chromosome is heterochromatic and has few genes. Despite these limitations, there remains ample opportunity for natural selection to act on the genes that are vital for male fertility and on Y factors that modulate gene expression elsewhere in the genome. Y chromosomes of many organisms have low levels of nucleotide variability, but a formal survey of D. melanogaster Y chromosome variation had yet to be performed. Here we surveyed Y-linked variation in six populations of D. melanogaster spread across the globe. We find surprisingly low levels of variability in African relative to Cosmopolitan (i.e., non-African) populations. While the low levels of Cosmopolitan Y chromosome polymorphism can be explained by the demographic histories of these populations, the staggeringly low polymorphism of African Y chromosomes cannot be explained by demographic history. An explanation that is entirely consistent with the data is that the Y chromosomes of Zimbabwe and Uganda populations have experienced recent selective sweeps. Interestingly, the Zimbabwe and Uganda Y chromosomes differ: in Zimbabwe, a European Y chromosome appears to have swept through the population.     

Inferring the demographic history and rate of adaptive substitution in Drosophila

An important goal of population genetics is to determine the forces that have shaped the pattern of genetic variation in natural populations. We developed a maximum likelihood method that allows us to infer demographic changes and detect recent positive selection (selective sweeps) in populations of varying size from DNA polymorphism data. Applying this approach to single nucleotide polymorphism data at more than 250 noncoding loci on the X chromosome of Drosophila melanogaster from an (ancestral) African population and a (derived) European, we found that the African population expanded about 60,000 y ago and that the European population split off from the African lineage about 15,800 y ago, thereby suffering a severe population size bottleneck. We estimated that about 160 beneficial mutations (with selection coefficients s between 0.05% and 0.5%) were fixed in the euchromatic portion of the X in the African population since population size expansion, and about 60 mutations (with s around 0.5%) in the diverging European lineage.     

North African Populations Carry the Signature of Admixture with Neandertals

One of the main findings derived from the analysis of the Neandertal genome was the evidence for admixture between Neandertals and non-African modern humans. An alternative scenario is that the ancestral population of non-Africans was closer to Neandertals than to Africans because of ancient population substructure. Thus, the study of North African populations is crucial for testing both hypotheses. We analyzed a total of 780,000 SNPs in 125 individuals representing seven different North African locations and searched for their ancestral/derived state in comparison to different human populations and Neandertals. We found that North African populations have a significant excess of derived alleles shared with Neandertals, when compared to sub-Saharan Africans. This excess is similar to that found in non-African humans, a fact that can be interpreted as a sign of Neandertal admixture. Furthermore, the Neandertal''s genetic signal is higher in populations with a local, pre-Neolithic North African ancestry. Therefore, the detected ancient admixture is not due to recent Near Eastern or European migrations. Sub-Saharan populations are the only ones not affected by the admixture event with Neandertals.     

Extent of linkage disequilibrium between the androgen receptor gene CAG and GGC repeats in human populations: implications for prostate cancer risk

While studies have implicated alleles at the CAG and GGC trinucleotide repeats of the androgen receptor gene with high-grade, aggressive prostate cancer disease, little is known about the normal range of variation for these two loci, which are separated by about 1.1 kb. More importantly, few data exist on the extent of linkage disequilibrium (LD) between the two loci in different human populations. Here we present data on CAG and GGC allelic variation and LD in six diverse populations. Alleles at the CAG and GGC repeat loci of the androgen receptor were typed in over 1000 chromosomes from Africa, Asia, and North America. Levels of linkage disequilibrium between the two loci were compared between populations. Haplotype variation and diversity were estimated for each population. Our results reveal that populations of African descent possess significantly shorter alleles for the two loci than non-African populations (P<0.0001). Allelic diversity for both markers was higher among African Americans than any other population, including indigenous Africans from Sierra Leone and Nigeria. Analysis of molecular variance revealed that approx. 20% of CAG and GGC repeat variance could be attributed to differences between the populations. All non-African populations possessed the same common haplotype while the three populations of African descent possessed three divergent common haplotypes. Significant LD was observed in our sample of healthy African Americans. The LD observed in the African American population may be due to several reasons; recent migration of African Americans from diverse rural communities following urbanization, recurrent gene flow from diverse West African populations, and admixture with European Americans. This study represents the largest genotyping effort to be performed on the two androgen receptor trinucleotide repeat loci in diverse human populations.     

Early Back-to-Africa Migration into the Horn of Africa

Genetic studies have identified substantial non-African admixture in the Horn of Africa (HOA). In the most recent genomic studies, this non-African ancestry has been attributed to admixture with Middle Eastern populations during the last few thousand years. However, mitochondrial and Y chromosome data are suggestive of earlier episodes of admixture. To investigate this further, we generated new genome-wide SNP data for a Yemeni population sample and merged these new data with published genome-wide genetic data from the HOA and a broad selection of surrounding populations. We used multidimensional scaling and ADMIXTURE methods in an exploratory data analysis to develop hypotheses on admixture and population structure in HOA populations. These analyses suggested that there might be distinct, differentiated African and non-African ancestries in the HOA. After partitioning the SNP data into African and non-African origin chromosome segments, we found support for a distinct African (Ethiopic) ancestry and a distinct non-African (Ethio-Somali) ancestry in HOA populations. The African Ethiopic ancestry is tightly restricted to HOA populations and likely represents an autochthonous HOA population. The non-African ancestry in the HOA, which is primarily attributed to a novel Ethio-Somali inferred ancestry component, is significantly differentiated from all neighboring non-African ancestries in North Africa, the Levant, and Arabia. The Ethio-Somali ancestry is found in all admixed HOA ethnic groups, shows little inter-individual variance within these ethnic groups, is estimated to have diverged from all other non-African ancestries by at least 23 ka, and does not carry the unique Arabian lactase persistence allele that arose about 4 ka. Taking into account published mitochondrial, Y chromosome, paleoclimate, and archaeological data, we find that the time of the Ethio-Somali back-to-Africa migration is most likely pre-agricultural.     

Inferring human population sizes, divergence times and rates of gene flow from mitochondrial, X and Y chromosome resequencing data

We estimate parameters of a general isolation-with-migration model using resequence data from mitochondrial DNA (mtDNA), the Y chromosome, and two loci on the X chromosome in samples of 25-50 individuals from each of 10 human populations. Application of a coalescent-based Markov chain Monte Carlo technique allows simultaneous inference of divergence times, rates of gene flow, as well as changes in effective population size. Results from comparisons between sub-Saharan African and Eurasian populations estimate that 1500 individuals founded the ancestral Eurasian population approximately 40 thousand years ago (KYA). Furthermore, these small Eurasian founding populations appear to have grown much more dramatically than either African or Oceanian populations. Analyses of sub-Saharan African populations provide little evidence for a history of population bottlenecks and suggest that they began diverging from one another upward of 50 KYA. We surmise that ancestral African populations had already been geographically structured prior to the founding of ancestral Eurasian populations. African populations are shown to experience low levels of mitochondrial DNA gene flow, but high levels of Y chromosome gene flow. In particular, Y chromosome gene flow appears to be asymmetric, i.e., from the Bantu-speaking population into other African populations. Conversely, mitochondrial gene flow is more extensive between non-African populations, but appears to be absent between European and Asian populations.     

Demographic Inference Reveals African and European Admixture in the North American Drosophila melanogaster Population

Drosophila melanogaster spread from sub-Saharan Africa to the rest of the world colonizing new environments. Here, we modeled the joint demography of African (Zimbabwe), European (The Netherlands), and North American (North Carolina) populations using an approximate Bayesian computation (ABC) approach. By testing different models (including scenarios with continuous migration), we found that admixture between Africa and Europe most likely generated the North American population, with an estimated proportion of African ancestry of 15%. We also revisited the demography of the ancestral population (Africa) and found—in contrast to previous work—that a bottleneck fits the history of the population of Zimbabwe better than expansion. Finally, we compared the site-frequency spectrum of the ancestral population to analytical predictions under the estimated bottleneck model.     

Non-African populations of Drosophila melanogaster have a unique origin

Drosophila melanogaster is widely used as a model in DNA variation studies. Patterns of polymorphism have, however, been affected by the history of this species, which is thought to have recently spread out of Africa to the rest of the world. We analyzed DNA sequence variation in 11 populations, including four continental African and seven non-African samples (including Madagascar), at four independent X-linked loci. Variation patterns at all four loci followed neutral expectations in all African populations, but departed from it in all non-African ones due to a marked haplotype dimorphism at three out of four loci. We also found that all non-African populations show the same major haplotypes, though in various frequencies. A parsimonious explanation for these observations is that all non-African populations are derived from a single ancestral population having undergone a substantial reduction of polymorphism, probably through a bottleneck. Less likely alternatives involve either selection at all four loci simultaneously (including balancing selection at three of them), or admixture between two divergent populations. Small but significant structure was observed among African populations, and there were indications of differentiation across Eurasia for non-African ones. Since population history may result in non-equilibrium variation patterns, our study confirms that the search for footprints of selection in the D. melanogaster genome must include a sufficient understanding of its history.     

Inferring the population structure of the Maghreb sheep breeds using a medium‐density SNP chip

North Africa has a great diversity of indigenous sheep breeds whose origin is linked to its environmental characteristics and to certain historical events that took place in the region. To date, few genome‐wide studies have been conducted to investigate the population structure of North African indigenous sheep. The objective of the present study was to provide a detailed assessment of the genetic structure and admixture patterns of six Maghreb sheep populations using the Illumina 50K Ovine BeadChip and comparisons with 22 global populations of sheep and mouflon. Regardless of the method of analysis used, patterns of multiple hybridization events were observed within all North African populations, leading to a heterogeneous genetic architecture that varies according to the breed. The Barbarine population showed the lowest genetic heterogeneity and major southwest Asian ancestry, providing additional support to the Asian origin of the North African fat‐tailed sheep. All other breeds presented substantial Merino introgression ranging from 15% for D'man to 31% for Black Thibar. We highlighted several signals of ancestral introgression between North African and southern European sheep. In addition, we identified two opposite gradients of ancestry, southwest Asian and central European, occurring between North Africa and central Europe. Our results provide further evidence of the weak global population structure of sheep resulting from high levels of gene flow among breeds occurring worldwide. At the regional level, signs of recent admixture among North African populations, resulting in a change of the original genomic architecture of minority breeds, were also detected.     

Secondary contact and local adaptation contribute to genome‐wide patterns of clinal variation in Drosophila melanogaster

Populations arrayed along broad latitudinal gradients often show patterns of clinal variation in phenotype and genotype. Such population differentiation can be generated and maintained by both historical demographic events and local adaptation. These evolutionary forces are not mutually exclusive and can in some cases produce nearly identical patterns of genetic differentiation among populations. Here, we investigate the evolutionary forces that generated and maintain clinal variation genome‐wide among populations of Drosophila melanogaster sampled in North America and Australia. We contrast patterns of clinal variation in these continents with patterns of differentiation among ancestral European and African populations. Using established and novel methods we derive here, we show that recently derived North America and Australia populations were likely founded by both European and African lineages and that this hybridization event likely contributed to genome‐wide patterns of parallel clinal variation between continents. The pervasive effects of admixture mean that differentiation at only several hundred loci can be attributed to the operation of spatially varying selection using an F_ST outlier approach. Our results provide novel insight into the well‐studied system of clinal differentiation in D. melanogaster and provide a context for future studies seeking to identify loci contributing to local adaptation in a wide variety of organisms, including other invasive species as well as temperate endemics.     

Estimates of African, European and Native American ancestry in Afro-Caribbean men on the island of Tobago

BACKGROUND/AIMS: The Tobago Afro-Caribbean population is a valuable resource for studying the genetics of diseases that show significant differences in prevalence between populations of African descent and populations of other ancestries. Empirical confirmation of low European and Native American admixture may help in clarifying the ethnic variation in risk for such diseases. We hypothesize that the degree of European and Native American admixture in the Tobago population is low. METHODS: Admixture was estimated in a random sample of 220 men, from a population-based prostate cancer screening survey of 3,082 Tobago males, aged 40 to 79 years. We used a set of six autosomal markers with large allele frequency differences between the major ethnic populations involved in the admixture process, Europeans, Native Americans and West Africans. RESULTS: The ancestral proportions of Tobago population are estimated as 94.0+/-1.2% African, 4.6+/-3.4% European and 1.4+/-3.6% Native American. CONCLUSIONS: We conclude that Tobago Afro-Caribbean men are predominantly of West African ancestry, with minimal European and Native American admixture. The Tobago population, thus, may carry a higher burden of high-risk alleles of African origin for certain diseases than the more admixed African-American population. Conversely, this population may benefit from a higher prevalence of protective alleles of African origin.     

Human Genetic Ancestral Composition Correlates with the Origin of Mycobacterium leprae Strains in a Leprosy Endemic Population

Recent reports have suggested that leprosy originated in Africa, extended to Asia and Europe, and arrived in the Americas during European colonization and the African slave trade. Due to colonization, the contemporary Colombian population is an admixture of Native-American, European and African ancestries. Because microorganisms are known to accompany humans during migrations, patterns of human migration can be traced by examining genomic changes in associated microbes. The current study analyzed 118 leprosy cases and 116 unrelated controls from two Colombian regions endemic for leprosy (Atlantic and Andean) in order to determine possible associations of leprosy with patient ancestral background (determined using 36 ancestry informative markers), Mycobacterium leprae genotype and/or patient geographical origin. We found significant differences between ancestral genetic composition. European components were predominant in Andean populations. In contrast, African components were higher in the Atlantic region. M. leprae genotypes were then analyzed for cluster associations and compared with the ancestral composition of leprosy patients. Two M. leprae principal clusters were found: haplotypes C54 and T45. Haplotype C54 associated with African origin and was more frequent in patients from the Atlantic region with a high African component. In contrast, haplotype T45 associated with European origin and was more frequent in Andean patients with a higher European component. These results suggest that the human and M. leprae genomes have co-existed since the African and European origins of the disease, with leprosy ultimately arriving in Colombia during colonization. Distinct M. leprae strains followed European and African settlement in the country and can be detected in contemporary Colombian populations.     

Genomewide analysis of admixture and adaptation in the Africanized honeybee

Genetic exchange by hybridization or admixture can make an important contribution to evolution, and introgression of favourable alleles can facilitate adaptation to new environments. A small number of honeybees (Apis mellifera) with African ancestry were introduced to Brazil ~60 years ago, which dispersed and hybridized with existing managed populations of European origin, quickly spreading across much of the Americas in an example of a massive biological invasion. Here, we analyse whole‐genome sequences of 32 Africanized honeybees sampled from throughout Brazil to study the effect of this process on genome diversity. By comparison with ancestral populations from Europe and Africa, we infer that these samples have 84% African ancestry, with the remainder from western European populations. However, this proportion varies across the genome and we identify signals of positive selection in regions with high European ancestry proportions. These observations are largely driven by one large gene‐rich 1.4‐Mbp segment on chromosome 11 where European haplotypes are present at a significantly elevated frequency and likely confer an adaptive advantage in the Africanized honeybee population. This region has previously been implicated in reproductive traits and foraging behaviour in worker bees. Finally, by analysing the distribution of ancestry tract lengths in the context of the known time of the admixture event, we are able to infer an average generation time of 2.0 years. Our analysis highlights the processes by which populations of mixed genetic ancestry form and adapt to new environments.     

Genetic Variation on Chromosome 6 Influences F Cell Levels in Healthy Individuals of African Descent and HbF Levels in Sickle Cell Patients

Fetal haemoglobin (HbF) is a major ameliorating factor in sickle cell disease. We investigated if a quantitative trait locus on chromosome 6q23 was significantly associated with HbF and F cell levels in individuals of African descent. Single nucleotide polymorphisms (SNPs) in a 24-kb intergenic region, 33-kb upstream of the HBS1L gene and 80-kb upstream of the MYB gene, were typed in 177 healthy Afro-Caribbean subjects (AC) of approximately 7% European admixture, 631 healthy Afro-Germans (AG, a group of African and German descendents located in rural Jamaica with about 20% European admixture), 87 West African and Afro-Caribbean individuals with sickle cell anaemia (HbSS), as well as 75 Northern Europeans, which served as a contrasting population. Association with a tag SNP for the locus was detected in all four groups (AC, P = 0.005, AG, P = 0.002, HbSS patients, P = 0.019, Europeans, P = 1.5×10⁻⁷). The association signal varied across the interval in the African-descended groups, while it is more uniform in Europeans. The 6q QTL for HbF traits is present in populations of African origin and is also acting in sickle cell anaemia patients. We have started to distinguish effects originating from European and African ancestral populations in our admixed study populations.     

Selective sorting of ancestral introgression in maize and teosinte along an elevational cline

While often deleterious, hybridization can also be a key source of genetic variation and pre-adapted haplotypes, enabling rapid evolution and niche expansion. Here we evaluate these opposing selection forces on introgressed ancestry between maize (Zea mays ssp. mays) and its wild teosinte relative, mexicana (Zea mays ssp. mexicana). Introgression from ecologically diverse teosinte may have facilitated maize’s global range expansion, in particular to challenging high elevation regions (> 1500 m). We generated low-coverage genome sequencing data for 348 maize and mexicana individuals to evaluate patterns of introgression in 14 sympatric population pairs, spanning the elevational range of mexicana, a teosinte endemic to the mountains of Mexico. While recent hybrids are commonly observed in sympatric populations and mexicana demonstrates fine-scale local adaptation, we find that the majority of mexicana ancestry tracts introgressed into maize over 1000 generations ago. This mexicana ancestry seems to have maintained much of its diversity and likely came from a common ancestral source, rather than contemporary sympatric populations, resulting in relatively low F_ST between mexicana ancestry tracts sampled from geographically distant maize populations.Introgressed mexicana ancestry in maize is reduced in lower-recombination rate quintiles of the genome and around domestication genes, consistent with pervasive selection against introgression. However, we also find mexicana ancestry increases across the sampled elevational gradient and that high introgression peaks are most commonly shared among high-elevation maize populations, consistent with introgression from mexicana facilitating adaptation to the highland environment. In the other direction, we find patterns consistent with adaptive and clinal introgression of maize ancestry into sympatric mexicana at many loci across the genome, suggesting that maize also contributes to adaptation in mexicana, especially at the lower end of its elevational range. In sympatric maize, in addition to high introgression regions we find many genomic regions where selection for local adaptation maintains steep gradients in introgressed mexicana ancestry across elevation, including at least two inversions: the well-characterized 14 Mb Inv4m on chromosome 4 and a novel 3 Mb inversion Inv9f surrounding the macrohairless1 locus on chromosome 9. Most outlier loci with high mexicana introgression show no signals of sweeps or local sourcing from sympatric populations and so likely represent ancestral introgression sorted by selection, resulting in correlated but distinct outcomes of introgression in different contemporary maize populations.     

Population genomic analysis uncovers African and European admixture in Drosophila melanogaster populations from the south‐eastern United States and Caribbean Islands

Drosophila melanogaster is postulated to have colonized North America in the past several 100 years in two waves. Flies from Europe colonized the east coast United States while flies from Africa inhabited the Caribbean, which if true, make the south‐east US and Caribbean Islands a secondary contact zone for African and European D. melanogaster. This scenario has been proposed based on phenotypes and limited genetic data. In our study, we have sequenced individual whole genomes of flies from populations in the south‐east US and Caribbean Islands and examined these populations in conjunction with population sequences from the west coast US, Africa, and Europe. We find that west coast US populations are closely related to the European population, likely reflecting a rapid westward expansion upon first settlements into North America. We also find genomic evidence of African and European admixture in south‐east US and Caribbean populations, with a clinal pattern of decreasing proportions of African ancestry with higher latitude. Our genomic analysis of D. melanogaster populations from the south‐east US and Caribbean Islands provides more evidence for the Caribbean Islands as the source of previously reported novel African alleles found in other east coast US populations. We also find the border between the south‐east US and the Caribbean island to be the admixture hot zone where distinctly African‐like Caribbean flies become genomically more similar to European‐like south‐east US flies. Our findings have important implications for previous studies examining the generation of east coast US clines via selection.     

Latitudinal variability of Drosophila melanogaster: Allozyme frequencies divergence between European and Afrotropical populations

Allelic frequencies at five polymorphic loci were determined in seven European and six Afrotropical populations of Drosophila melanogaster. African populations, which may be considered as ancestral for the species, showed a greater genetic diversity as measured by the number of alleles found. Within each geographic group (Europe or tropical Africa) genetic distances between local populations were very small (D=0.027). By contrast, the average distance between European and African populations (D=0.389) was more than 12 times bigger. It was previously known that various morphological or physiological differences, which probably reflect genetic adaptations to different environments, exist between these temperate and tropical populations. Data presented here suggest that the divergence in allozyme frequencies also reflects some selective mechanisms.     

X-Linked MTMR8 Diversity and Evolutionary History of Sub-Saharan Populations

The genetic diversity within an 11 kb segment of the MTMR8 gene in a sample of 111 sub-Saharan and 49 non-African X chromosomes was investigated to assess the early evolutionary history of sub-Saharan Africans and the out-of-Africa expansion. The analyses revealed a complex genetic structure of the Africans that contributed to the emergence of modern humans. We observed partitioning of two thirds of old lineages among southern, west/central and east African populations indicating ancient population stratification predating the out of Africa migration. Age estimates of these lineages, older than coalescence times of uniparentally inherited markers, raise the question whether contemporary humans originated from a single population or as an amalgamation of different populations separated by years of independent evolution, thus suggesting a greater antiquity of our species than generally assumed. While the oldest sub-Saharan lineages, ∼500 thousand years, are found among Khoe-San from southern-Africa, a distinct haplotype found among Biaka is likely due to admixture from an even older population. An East African population that gave rise to non-Africans underwent a selective sweep affecting the subcentromeric region where MTMR8 is located. This and similar sweeps in four other regions of the X chromosome, documented in the literature, effectively reduced genetic diversity of non-African chromosomes and therefore may have exacerbated the effect of the demographic bottleneck usually ascribed to the out of Africa migration. Our data is suggestive, however, that a bottleneck, occurred in Africa before range expansion.