A genome scan to detect candidate regions influenced by local natural selection in human populations

As human populations dispersed throughout the world, they were subjected to new selective forces, which must have led to local adaptation via natural selection and hence altered patterns of genetic variation. Yet, there are very few examples known in which such local selection has clearly influenced human genetic variation. A potential approach for detecting local selection is to screen random loci across the genome; those loci that exhibit unusually large genetic distances between human populations are then potential markers of genomic regions under local selection. We investigated this approach by genotyping 332 short tandem repeat (STR) loci in Africans and Europeans and calculating the genetic differentiation for each locus. Patterns of genetic diversity at these loci were consistent with greater variation in Africa and with local selection operating on populations as they moved out of Africa. For 11 loci exhibiting the largest genetic differences, we genotyped an additional STR locus located nearby; the genetic distances for these nearby loci were significantly larger than average. These genomic regions therefore reproducibly exhibit larger genetic distances between populations than the "average" genomic region, consistent with local selection. Our results demonstrate that genome scans are a promising means of identifying candidate regions that have been subjected to local selection.     

Approximate Bayesian inference reveals evidence for a recent, severe bottleneck in a Netherlands population of Drosophila melanogaster

Genome-wide nucleotide variation in non-African populations of Drosophila melanogaster is a subset of variation found in East sub-Saharan African populations, suggesting a bottleneck in the history of the former. We implement an approximate Bayesian approach to infer the timing, duration, and severity of this putative bottleneck and ask whether this inferred model is sufficient to account for patterns of variability observed at 115 loci scattered across the X chromosome. We estimate a recent bottleneck 0.006N(e) generations ago, somewhat further in the past than suggested by biogeographical evidence. Using various proposed statistical tests, we find that this bottleneck model is able to predict the majority of observed features of diversity and linkage disequilibrium in the data. Thus, while precise estimates of bottleneck parameters (like inferences of selection) are sensitive to model assumptions, our results imply that it may be unnecessary to invoke frequent selective sweeps associated with the dispersal of D. melanogaster from Africa to explain patterns of variability in non-African populations.     

Natural selection at genomic regions associated with obesity and type-2 diabetes: East Asians and sub-Saharan Africans exhibit high levels of differentiation at type-2 diabetes regions

Different populations suffer from different rates of obesity and type-2 diabetes (T2D). Little is known about the genetic or adaptive component, if any, that underlies these differences. Given the cultural, geographic, and dietary variation that accumulated among humans over the last 60,000 years, we examined whether loci identified by genome-wide association studies for these traits have been subject to recent selection pressures. Using genome-wide SNP data on 938 individuals in 53 populations from the Human Genome Diversity Panel, we compare population differentiation and haplotype patterns at these loci to the rest of the genome. Using an “expanding window” approach (100–1,600 kb) for the individual loci as well as the loci as ensembles, we find a high degree of differentiation for the ensemble of T2D loci. This differentiation is most pronounced for East Asians and sub-Saharan Africans, suggesting that these groups experienced natural selection at loci associated with T2D. Haplotype analysis suggests an excess of obesity loci with evidence of recent positive selection among South Asians and Europeans, compared to sub-Saharan Africans and Native Americans. We also identify individual loci that may have been subjected to natural selection, such as the T2D locus, HHEX, which displays both elevated differentiation and extended haplotype homozygosity in comparisons of East Asians with other groups. Our findings suggest that there is an evolutionary genetic basis for population differences in these traits, and we have identified potential group-specific genetic risk factors.     

A microsatellite variability screen for positive selection associated with the "out of Africa" habitat expansion of Drosophila melanogaster

We report a "hitchhiking mapping" study in D. melanogaster, which searches for genomic regions with reduced variability. The study's aim was to identify selective sweeps associated with the "out of Africa" habitat expansion. We scanned 103 microsatellites on chromosome 3 and 102 microsatellites on the X chromosome for reduced variability in non-African populations. When the chromosomes were analyzed separately, the number of loci with a significant reduction in variability only slightly exceeded the expectation under neutrality--six loci on the third chromosome and four loci on the X chromosome. However, non-African populations also have a more pronounced average loss in variability on the X chromosomes as compared to the third chromosome, which suggests the action of selection. Therefore, comparing the X chromosome to the autosome yields a higher number of significantly reduced loci. However, a more pronounced loss of variability on the X chromosome may be caused by demographic events rather than by natural selection. We therefore explored a range of demographic scenarios and found that some of these captured most, but not all aspects of our data. More theoretical work is needed to evaluate how demographic events might differentially affect X chromosomes and autosomes and to estimate the most likely scenario associated with the out of Africa expansion of D. melanogaster.     

Limited Evidence for Classic Selective Sweeps in African Populations

While hundreds of loci have been identified as reflecting strong-positive selection in human populations, connections between candidate loci and specific selective pressures often remain obscure. This study investigates broader patterns of selection in African populations, which are underrepresented despite their potential to offer key insights into human adaptation. We scan for hard selective sweeps using several haplotype and allele-frequency statistics with a data set of nearly 500,000 genome-wide single-nucleotide polymorphisms in 12 highly diverged African populations that span a range of environments and subsistence strategies. We find that positive selection does not appear to be a strong determinant of allele-frequency differentiation among these African populations. Haplotype statistics do identify putatively selected regions that are shared across African populations. However, as assessed by extensive simulations, patterns of haplotype sharing between African populations follow neutral expectations and suggest that tails of the empirical distributions contain false-positive signals. After highlighting several genomic regions where positive selection can be inferred with higher confidence, we use a novel method to identify biological functions enriched among populations’ empirical tail genomic windows, such as immune response in agricultural groups. In general, however, it seems that current methods for selection scans are poorly suited to populations that, like the African populations in this study, are affected by ascertainment bias and have low levels of linkage disequilibrium, possibly old selective sweeps, and potentially reduced phasing accuracy. Additionally, population history can confound the interpretation of selection statistics, suggesting that greater care is needed in attributing broad genetic patterns to human adaptation.     

Signatures of Convergent Evolution and Natural Selection at the Alcohol Dehydrogenase Gene Region are Correlated with Agriculture in Ethnically Diverse Africans

The alcohol dehydrogenase (ADH) family of genes encodes enzymes that catalyze the metabolism of ethanol into acetaldehyde. Nucleotide variation in ADH genes can affect the catalytic properties of these enzymes and is associated with a variety of traits, including alcoholism and cancer. Some ADH variants, including the ADH1B*48His (rs1229984) mutation in the ADH1B gene, reduce the risk of alcoholism and are under positive selection in multiple human populations. The advent of Neolithic agriculture and associated increase in fermented foods and beverages is hypothesized to have been a selective force acting on such variants. However, this hypothesis has not been tested in populations outside of Asia. Here, we use genome-wide selection scans to show that the ADH gene region is enriched for variants showing strong signals of positive selection in multiple Afroasiatic-speaking, agriculturalist populations from Ethiopia, and that this signal is unique among sub-Saharan Africans. We also observe strong selection signals at putatively functional variants in nearby lipid metabolism genes, which may influence evolutionary dynamics at the ADH region. Finally, we show that haplotypes carrying these selected variants were introduced into Northeast Africa from a West-Eurasian source within the last ∼2,000 years and experienced positive selection following admixture. These selection signals are not evident in nearby, genetically similar populations that practice hunting/gathering or pastoralist subsistence lifestyles, supporting the hypothesis that the emergence of agriculture shapes patterns of selection at ADH genes. Together, these results enhance our understanding of how adaptations to diverse environments and diets have influenced the African genomic landscape.     

Worldwide patterns of genetic differentiation imply multiple 'domestications' of Aedes aegypti, a major vector of human diseases

Understanding the processes by which species colonize and adapt to human habitats is particularly important in the case of disease-vectoring arthropods. The mosquito species Aedes aegypti, a major vector of dengue and yellow fever viruses, probably originated as a wild, zoophilic species in sub-Saharan Africa, where some populations still breed in tree holes in forested habitats. Many populations of the species, however, have evolved to thrive in human habitats and to bite humans. This includes some populations within Africa as well as almost all those outside Africa. It is not clear whether all domestic populations are genetically related and represent a single 'domestication' event, or whether association with human habitats has developed multiple times independently within the species. To test the hypotheses above, we screened 24 worldwide population samples of Ae. aegypti at 12 polymorphic microsatellite loci. We identified two distinct genetic clusters: one included all domestic populations outside of Africa and the other included both domestic and forest populations within Africa. This suggests that human association in Africa occurred independently from that in domestic populations across the rest of the world. Additionally, measures of genetic diversity support Ae. aegypti in Africa as the ancestral form of the species. Individuals from domestic populations outside Africa can reliably be assigned back to their population of origin, which will help determine the origins of new introductions of Ae. aegypti.     

History and structure of sub-Saharan populations of Drosophila melanogaster

Drosophila melanogaster is an important model organism in evolutionary genetics, yet little is known about the population structure and the demographic history of this species within sub-Saharan Africa, which is thought to contain its ancestral range. We surveyed nucleotide variation at four 1-kb fragments in 240 individual lines representing 21 sub-Saharan and 4 Palearctic population samples of D. melanogaster. In agreement with recent studies, we find a small but significant level of genetic differentiation within sub-Saharan Africa. A clear geographic pattern is observed, with eastern and western African populations composing two genetically distinct groups. This pattern may have resulted from a relatively recent establishment of D. melanogaster in western Africa. Eastern populations show greater evidence for long-term stability, consistent with the hypothesis that eastern Africa contains the ancestral range of the species. Three sub-Saharan populations show evidence for cosmopolitan introgression. Apart from those cases, the closest relationships between Palearctic and sub-Saharan populations involve a sample from the rift zone (Uganda), suggesting that the progenitors of Palearctic D. melanogaster might have come from this region. Finally, we find a large excess of singleton polymorphisms in the full data set, which is best explained by a combination of population growth and purifying selection.     

Nucleotide variation at Msn and Alas2, two genes flanking the centromere of the X chromosome in humans

The centromeric region of the X chromosome in humans experiences low rates of recombination over a considerable physical distance. In such a region, the effects of selection may extend to linked sites that are far away. To investigate the effects of this recombinational environment on patterns of nucleotide variability, we sequenced 4581 bp at Msn and 4697 bp at Alas2, two genes situated on either side of the X chromosome centromere, in a worldwide sample of 41 men, as well as in one common chimpanzee and one orangutan. To investigate patterns of linkage disequilibrium (LD) across the centromere, we also genotyped several informative sites from each gene in 120 men from sub-Saharan Africa. By studying X-linked loci in males, we were able to recover haplotypes and study long-range patterns of LD directly. Overall patterns of variability were remarkably similar at these two loci. Both loci exhibited (i) very low levels of nucleotide diversity (among the lowest seen in the human genome); (ii) a strong skew in the distribution of allele frequencies, with an excess of both very-low and very-high-frequency derived alleles in non-African populations; (iii) much less variation in the non-African than in the African samples; (iv) very high levels of population differentiation; and (v) complete LD among all sites within loci. We also observed significant LD between Msn and Alas2 in Africa, despite the fact that they are separated by approximately 10 Mb. These observations are difficult to reconcile with a simple demographic model but may be consistent with positive and/or purifying selection acting on loci within this large region of low recombination.     

The extent of linkage disequilibrium caused by selection on G6PD in humans

The gene coding for glucose-6-phosphate dehydrogenase (G6PD) is subject to positive selection by malaria in some human populations. The G6PD A- allele, which is common in sub-Saharan Africa, is associated with deficient enzyme activity and protection from severe malaria. To delimit the impact of selection on patterns of linkage disequilibrium (LD) and nucleotide diversity, we resequenced 5.1 kb at G6PD and approximately 2-3 kb at each of eight loci in a 2.5-Mb region roughly centered on G6PD in a diverse sub-Saharan African panel of 51 unrelated men (including 20 G6PD A-, 11 G6PD A+, and 20 G6PD B chromosomes). The signature of selection is evident in the absence of genetic variation at G6PD and at three neighboring loci within 0.9 Mb from G6PD among all individuals bearing G6PD A- alleles. A genomic region of approximately 1.6 Mb around G6PD was characterized by long-range LD associated with the A- alleles. These patterns of nucleotide variability and LD suggest that G6PD A- is younger than previous age estimates and has increased in frequency in sub-Saharan Africa due to strong selection (0.1 < s < 0.2). These results also show that selection can lead to nonrandom associations among SNPs over great physical and genetic distances, even in African populations.     

Effects of natural selection and gene conversion on the evolution of human glycophorins coding for MNS blood polymorphisms in malaria-endemic African populations

Malaria has been a very strong selection pressure in recent human evolution, particularly in Africa. Of the one million deaths per year due to malaria, more than 90% are in sub-Saharan Africa, a region with high levels of genetic variation and population substructure. However, there have been few studies of nucleotide variation at genetic loci that are relevant to malaria susceptibility across geographically and genetically diverse ethnic groups in Africa. Invasion of erythrocytes by Plasmodium falciparum parasites is central to the pathology of malaria. Glycophorin A (GYPA) and B (GYPB), which determine MN and Ss blood types, are two major receptors that are expressed on erythrocyte surfaces and interact with parasite ligands. We analyzed nucleotide diversity of the glycophorin gene family in 15 African populations with different levels of malaria exposure. High levels of nucleotide diversity and gene conversion were found at these genes. We observed divergent patterns of genetic variation between these duplicated genes and between different extracellular domains of GYPA. Specifically, we identified fixed adaptive changes at exons 3-4 of GYPA. By contrast, we observed an allele frequency spectrum skewed toward a significant excess of intermediate-frequency alleles at GYPA exon 2 in many populations; the degree of spectrum distortion is correlated with malaria exposure, possibly because of the joint effects of gene conversion and balancing selection. We also identified a haplotype causing three amino acid changes in the extracellular domain of glycophorin B. This haplotype might have evolved adaptively in five populations with high exposure to malaria.     

Genetic differentiation of Anopheles gambiae s.s. populations in Mali,West Africa,using microsatellite loci

Anopheles gambiae sensu stricto is a principal vector of malaria through much of sub-Saharan Africa, where this disease is a major cause of morbidity and mortality in human populations. Accordingly, population sizes and gene flow in this species have received special attention, as these parameters are important in attempts to control malaria by impacting its mosquito vector. Past measures of genetic differentiation have sometimes yielded conflicting results, in some cases suggesting that gene flow is extensive over vast distances (6000 km) and is disrupted only by major geological disturbances and/or barriers. Using microsatellite DNA loci from populations in Mali, West Africa, we measured genetic differentiation over uniform habitats favorable to the species across distances ranging from 62 to 536 km. Gene flow was strongly correlated with distance (r(2) = 0.77), with no major differences among chromosomes. We conclude that in this part of Africa, at least, genetic differentiation for microsatellite DNA loci is consistent with traditional models of isolation by distance.     

Contrasting patterns of X-linked and autosomal nucleotide variation in Drosophila melanogaster and Drosophila simulans

Surveys of molecular variation in Drosophila melanogaster and Drosophila simulans have suggested that diversity outside of Africa is a subset of that within Africa. It has been argued that reduced levels of diversity in non-African populations reflect a population bottleneck, adaptation to temperate climates, or both. Here, I summarize the available single-nucleotide polymorphism data for both species. A simple "out of Africa" bottleneck scenario is consistent with geographic patterns for loci on the X chromosome but not with loci on the autosomes. Interestingly, there is a trend toward lower nucleotide diversity on the X chromosome relative to autosomes in non-African populations of D. melanogaster, but the opposite trend is seen in African populations. In African populations, autosomal inversion polymorphisms in D. melanogaster may contribute to reduced autosome diversity relative to the X chromosome. To elucidate the role that selection might play in shaping patterns of variability, I present a summary of within- and between-species patterns of synonymous and replacement variation in both species. Overall, D. melanogaster autosomes harbor an excess of amino acid replacement polymorphisms relative to D. simulans. Interestingly, range expansion from Africa appears to have had little effect on synonymous-to-replacement polymorphism ratios.     

Pleiotropic Effects of Immune Responses Explain Variation in the Prevalence of Fibroproliferative Diseases

Many diseases are differentially distributed among human populations. Differential selection on genetic variants in ancestral environments that coincidentally predispose to disease can be an underlying cause of these unequal prevalence patterns. Selected genes may be pleiotropic, affecting multiple phenotypes and resulting in more than one disease or trait. Patterns of pleiotropy may be helpful in understanding the underlying causes of an array of conditions in a population. For example, several fibroproliferative diseases are more prevalent and severe in populations of sub-Saharan ancestry. We propose that this disparity is due to selection for an enhanced Th2 response that confers resistance to helminthic infections, and concurrently increases susceptibility to fibrosis due to the profibrotic action of Th2 cytokines. Many studies on selection of Th2-related genes for host resistance to helminths have been reported, but the pleiotropic impact of this selection on the distribution of fibrotic disorders has not been explicitly investigated. We discuss the disproportionate occurrence of fibroproliferative diseases in individuals of African ancestry and provide evidence that adaptation of the immune system has shaped the genetic structure of these human populations in ways that alter the distribution of multiple fibroproliferative diseases.     

Population structure of African buffalo inferred from mtDNA sequences and microsatellite loci: high variation but low differentiation

The African buffalo ( Syncerus caffer ) is widespread throughout sub-Saharan Africa and is found in most major vegetation types, wherever permanent sources of water are available, making it physically able to disperse through a wide range of habitats. Despite this, the buffalo has been assumed to be strongly philopatric and to form large aggregations that remain within separate home ranges with little interchange between units, but the level of differentiation within the species is unknown. Genetic differences between populations were assessed using mitochondrial DNA (control region) sequence data and analysis of variation at six microsatellite loci among 11 localities in eastern and southern Africa. High levels of genetic variability were found, suggesting that reported severe population bottlenecks due to outbreak of rinderpest during the last century did not strongly reduce the genetic variability within the species. The high level of genetic variation within the species was found to be evenly distributed among populations and only at the continental level were we able to consistently detect significant differentiation, contrasting with the assumed philopatric behaviour of the buffalo. Results of mtDNA and microsatellite data were found to be congruent, disagreeing with the alleged male-biased dispersal. We propose that the observed pattern of the distribution of genetic variation between buffalo populations at the regional level can be caused by fragmentation of a previous panmictic population due to human activity, and at the continental level, reflects an effect of geographical distance between populations.     

High Genetic Variability of Schistosoma haematobium in Mali and Nigeria

Schistosoma haematobium is one of the most prevalent parasitic flatworms, infecting over 112 million people in Africa. However, little is known about the genetic diversity of natural S. haematobium populations from the human host because of the inaccessible location of adult worms in the host. We used 4 microsatellite loci to genotype individually pooled S. haematobium eggs directly from each patient sampled at 4 endemic locations in Africa. We found that the average allele number of individuals from Mali was significantly higher than that from Nigeria. In addition, no significant difference in allelic composition was detected among the populations within Nigeria; however, the allelic composition was significantly different between Mali and Nigeria populations. This study demonstrated a high level of genetic variability of S. haematobium in the populations from Mali and Nigeria, the 2 major African endemic countries, suggesting that geographical population differentiation may occur in the regions.     

Identifying genes underlying skin pigmentation differences among human populations

Skin pigmentation is a human phenotype that varies greatly among human populations and it has long been speculated that this variation is adaptive. We therefore expect the genes that contribute to these large differences in phenotype to show large allele frequency differences among populations and to possibly harbor signatures of positive selection. To identify the loci that likely contribute to among-population human skin pigmentation differences, we measured allele frequency differentiation among Europeans, Chinese and Africans for 24 human pigmentation genes from 2 publicly available, large scale SNP data sets. Several skin pigmentation genes show unusually large allele frequency differences among these populations. To determine whether these allele frequency differences might be due to selection, we employed a within-population test based on long-range haplotype structure and identified several outliers that have not been previously identified as putatively adaptive. Most notably, we identify the DCT gene as a candidate for recent positive selection in the Chinese. Moreover, our analyses suggest that it is likely that different genes are responsible for the lighter skin pigmentation found in different non-African populations. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Genomic ancestry of North Africans supports back-to-Africa migrations

North African populations are distinct from sub-Saharan Africans based on cultural, linguistic, and phenotypic attributes; however, the time and the extent of genetic divergence between populations north and south of the Sahara remain poorly understood. Here, we interrogate the multilayered history of North Africa by characterizing the effect of hypothesized migrations from the Near East, Europe, and sub-Saharan Africa on current genetic diversity. We present dense, genome-wide SNP genotyping array data (730,000 sites) from seven North African populations, spanning from Egypt to Morocco, and one Spanish population. We identify a gradient of likely autochthonous Maghrebi ancestry that increases from east to west across northern Africa; this ancestry is likely derived from "back-to-Africa" gene flow more than 12,000 years ago (ya), prior to the Holocene. The indigenous North African ancestry is more frequent in populations with historical Berber ethnicity. In most North African populations we also see substantial shared ancestry with the Near East, and to a lesser extent sub-Saharan Africa and Europe. To estimate the time of migration from sub-Saharan populations into North Africa, we implement a maximum likelihood dating method based on the distribution of migrant tracts. In order to first identify migrant tracts, we assign local ancestry to haplotypes using a novel, principal component-based analysis of three ancestral populations. We estimate that a migration of western African origin into Morocco began about 40 generations ago (approximately 1,200 ya); a migration of individuals with Nilotic ancestry into Egypt occurred about 25 generations ago (approximately 750 ya). Our genomic data reveal an extraordinarily complex history of migrations, involving at least five ancestral populations, into North Africa.