Divergent haplotypes and human history as revealed in a worldwide survey of X-linked DNA sequence variation

The population genetic history of a 10.1-kbp noncoding region of the human X chromosome was studied using the males of the HGDP-CEPH Human Genome Diversity Panel (672 individuals from 52 populations). The geographic distribution of patterns of variation was roughly consistent with previous studies, with the major exception that 1 highly divergent haplotype (haplotype X, hX) was observed at low frequency in widely scattered non-African populations and not at all observed in sub-Saharan African populations. Microsatellite (short tandem repeat) variation within the sequenced region was low among copies of hX, even though the estimated time of ancestry of hX and other sequences was 1.44 Myr. The estimated age of the common ancestor of all hX copies was 5,230 years (95% consistency index: 2,000-75,480 years). To further address the presence of hX in Africa, additional samples from Chad and Tanzania were screened. Five additional copies of hX were observed, consistent with a history in which hX was present in Africa prior to the migration of modern humans out of Africa and with eastern Africa being the source of non-African modern human populations. Taken together, these features of hX-that it is much older than other haplotypes and uncommon and patchily distributed throughout Africa, Europe, and Asia-present a cautionary tale for interpretations of human history.     

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.     

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.     

Spatial and temporal distribution of the neutral polymorphisms in the last ZFX intron: analysis of the haplotype structure and genealogy.

With 10 segregating sites (simple nucleotide polymorphisms) in the last intron (1089 bp) of the ZFX gene we have observed 11 haplotypes in 336 chromosomes representing a worldwide array of 15 human populations. Two haplotypes representing 77% of all chromosomes were distributed almost evenly among four continents. Five of the remaining haplotypes were detected in Africa and 4 others were restricted to Eurasia and the Americas. Using the information about the ancestral state of the segregating positions (inferred from human-great ape comparisons), we applied coalescent analysis to estimate the age of the polymorphisms and the resulting haplotypes. The oldest haplotype, with the ancestral alleles at all the sites, was observed at low frequency only in two groups of African origin. Its estimated age of 740 to 1100 kyr corresponded to the time to the most recent common ancestor. The two most frequent worldwide distributed haplotypes were estimated at 550 to 840 and 260 to 400 kyr, respectively, while the age of the continentally restricted polymorphisms was 120 to 180 kyr and smaller. Comparison of spatial and temporal distribution of the ZFX haplotypes suggests that modern humans diverged from the common ancestral stock in the Middle Paleolithic era. Subsequent range expansion prevented substantial gene flow among continents, separating African groups from populations that colonized Eurasia and the New World.     

High diversity of alpha-globin haplotypes in a Senegalese population, including many previously unreported variants.

RFLP haplotypes at the alpha-globin gene complex have been examined in 190 individuals from the Niokolo Mandenka population of Senegal: haplotypes were assigned unambiguously for 210 chromosomes. The Mandenka share with other African populations a sample size-independent haplotype diversity that is much greater than that in any non-African population: the number of haplotypes observed in the Mandenka is typically twice that seen in the non-African populations sampled to date. Of these haplotypes, 17.3% had not been observed in any previous surveys, and a further 19.1% have previously been reported only in African populations. The haplotype distribution shows clear differences between African and non-African peoples, but this is on the basis of population-specific haplotypes combined with haplotypes common to all. The relationship of the newly reported haplotypes to those previously recorded suggests that several mutation processes, particularly recombination as homologous exchange or gene conversion, have been involved in their production. A computer program based on the expectation-maximization (EM) algorithm was used to obtain maximum-likelihood estimates of haplotype frequencies for the entire data set: good concordance between the unambiguous and EM-derived sets was seen for the overall haplotype frequencies. Some of the low-frequency haplotypes reported by the estimation algorithm differ greatly, in structure, from those haplotypes known to be present in human populations, and they may not represent haplotypes actually present in the sample.     

Natural selection and population history in the human angiotensinogen gene (AGT): 736 complete AGT sequences in chromosomes from around the world

Several lines of evidence suggest that patterns of genetic variability in the human angiotensinogen gene (AGT) contribute to phenotypic variability in human hypertension. The A(-6) promoter variant of AGT is associated with higher plasma angiotensinogen levels and increased risk of essential hypertension. The geographic distribution of the A(-6) variant leads to the intriguing hypothesis that the G(-6) promoter variant has been selectively advantageous outside Africa. To test these hypotheses, we investigated the roles of population history and natural selection in shaping patterns of genetic diversity in AGT, by sequencing the entire AGT gene (14400 bp) in 736 chromosomes from Africa, Asia, and Europe. We found that the A(-6) variant is present at higher frequency in African populations than in non-African populations. Neutrality tests found no evidence of a departure from selective neutrality, when whole AGT sequences were compared. However, tests restricted to sites in the vicinity of the A(-6)G polymorphism found evidence of a selective sweep. Sliding-window analyses showed that evidence of the sweep is restricted to sites in tight linkage disequilibrium (LD) with the A(-6)G polymorphism. Further, haplotypes carrying the G(-6) variant showed elevated levels of LD, suggesting that they have risen recently to high frequency. Departures from neutral expectation in some but not all regions of AGT indicate that patterns of diversity in the gene cannot be accounted for solely by population history, which would affect all regions equally. Taken together, patterns of genetic diversity in AGT suggest that natural selection has generally favored the G(-6) variant over the A(-6) variant in non-African populations. However, important localized effects may also be present.     

Tracing the Route of Modern Humans out of Africa by Using 225 Human Genome Sequences from Ethiopians and Egyptians

The predominantly African origin of all modern human populations is well established, but the route taken out of Africa is still unclear. Two alternative routes, via Egypt and Sinai or across the Bab el Mandeb strait into Arabia, have traditionally been proposed as feasible gateways in light of geographic, paleoclimatic, archaeological, and genetic evidence. Distinguishing among these alternatives has been difficult. We generated 225 whole-genome sequences (225 at 8× depth, of which 8 were increased to 30×; Illumina HiSeq 2000) from six modern Northeast African populations (100 Egyptians and five Ethiopian populations each represented by 25 individuals). West Eurasian components were masked out, and the remaining African haplotypes were compared with a panel of sub-Saharan African and non-African genomes. We showed that masked Northeast African haplotypes overall were more similar to non-African haplotypes and more frequently present outside Africa than were any sets of haplotypes derived from a West African population. Furthermore, the masked Egyptian haplotypes showed these properties more markedly than the masked Ethiopian haplotypes, pointing to Egypt as the more likely gateway in the exodus to the rest of the world. Using five Ethiopian and three Egyptian high-coverage masked genomes and the multiple sequentially Markovian coalescent (MSMC) approach, we estimated the genetic split times of Egyptians and Ethiopians from non-African populations at 55,000 and 65,000 years ago, respectively, whereas that of West Africans was estimated to be 75,000 years ago. Both the haplotype and MSMC analyses thus suggest a predominant northern route out of Africa via Egypt.     

PKLR- GBA region shows almost complete linkage disequilibrium over 70 kb in a set of worldwide populations

Haplotype diversity in a genomic region of approximately 70 kb in 1q21 between genes PKLR and GBA was characterized by typing one single nucleotide polymorphism (SNP) in PKLR, two SNPs in GBA and one short tandem repeat polymorphism (STRP) in PKLR in 1792 chromosomes from 17 worldwide populations. Two other SNPs in GBA were typed in three African populations. Most chromosomes carried one of either two phylogenetically distinct haplotypes with different alleles at each site. Allele diversity at the STRP was tightly linked to haplotype background. Linkage disequilibrium (LD) was highly significant for all SNP pairs in all populations, although it was, on average, slightly higher in non-African populations than in sub-Saharan Africans. Variation at PKLR-GBA was also tightly linked to that at the GBA pseudogene, 16 kb downstream from GBA. Thus, a 90 kb-long LD block was observed, which points to a low recombination rate in this region. Detailed haplotype phylogeny suggests that the chimpanzee GBA haplotype is not one of the two most frequent haplotypes. Based on variability at the PKLR STRP and on the geographical distribution of LD, the expansion of the two main haplotypes may have predated the "Out of Africa" expansion of anatomically modern humans. LD and STRP variability in non-Africans are approximately 87% of those in Africans, in contrast with other loci; this implies that the "out of Africa" bottleneck may have had a broad distribution of effects across loci.     

Short tandem-repeat polymorphism/alu haplotype variation at the PLAT locus: implications for modern human origins

Two dinucleotide short tandem-repeat polymorphisms (STRPs) and a polymorphic Alu element spanning a 22-kb region of the PLAT locus on chromosome 8p12-q11.2 were typed in 1,287-1,420 individuals originating from 30 geographically diverse human populations, as well as in 29 great apes. These data were analyzed as haplotypes consisting of each of the dinucleotide repeats and the flanking Alu insertion/deletion polymorphism. The global pattern of STRP/Alu haplotype variation and linkage disequilibrium (LD) is informative for the reconstruction of human evolutionary history. Sub-Saharan African populations have high levels of haplotype diversity within and between populations, relative to non-Africans, and have highly divergent patterns of LD. Non-African populations have both a subset of the haplotype diversity present in Africa and a distinct pattern of LD. The pattern of haplotype variation and LD observed at the PLAT locus suggests a recent common ancestry of non-African populations, from a small population originating in eastern Africa. These data indicate that, throughout much of modern human history, sub-Saharan Africa has maintained both a large effective population size and a high level of population substructure. Additionally, Papua New Guinean and Micronesian populations have rare haplotypes observed otherwise only in African populations, suggesting ancient gene flow from Africa into Papua New Guinea, as well as gene flow between Melanesian and Micronesian populations.     

Fixation of the human-specific CMP-N-acetylneuraminic acid hydroxylase pseudogene and implications of haplotype diversity for human evolution

The human CMP-N-acetylneuraminic acid hydroxylase gene (CMAH) suffered deletion of an exon that encodes an active center for the enzyme approximately 3.2 million years ago (MYA). We analyzed a 7.3-kb intronic region of 132 CMAH genes to explore the fixation process of this pseudogene and the demographic implication of its haplotype diversity. Fifty-six variable sites were sorted into 18 different haplotypes with significant linkage disequilibrium. Despite the rather low nucleotide diversity, the most recent common ancestor at CMAH dates to 2.9 MYA. This deep genealogy follows shortly after the original exon deletion, indicating that the deletion has fixed in the population, although whether this fixation was facilitated by natural selection remains to be resolved. Remarkable features are exceptionally long persistence of two lineages and the confinement of one lineage in Africa, implying that some African local populations were in relative isolation while others were directly involved in multiple African exoduses of the genus Homo. Importantly, haplotypes found in Eurasia suggest interbreeding between then-contemporaneous human species. Although population structure within Africa complicates the interpretation of phylogeographic information of haplotypes, the data support a single origin of modern humans, but not with complete replacement of archaic inhabitants by modern humans.     

Features of evolution and expansion of modern humans,inferred from genomewide microsatellite markers

We study data on variation in 52 worldwide populations at 377 autosomal short tandem repeat loci, to infer a demographic history of human populations. Variation at di-, tri-, and tetranucleotide repeat loci is distributed differently, although each class of markers exhibits a decrease of within-population genetic variation in the following order: sub-Saharan Africa, Eurasia, East Asia, Oceania, and America. There is a similar decrease in the frequency of private alleles. With multidimensional scaling, populations belonging to the same major geographic region cluster together, and some regions permit a finer resolution of populations. When a stepwise mutation model is used, a population tree based on TD estimates of divergence time suggests that the branches leading to the present sub-Saharan African populations of hunter-gatherers were the first to diverge from a common ancestral population (approximately 71-142 thousand years ago). The branches corresponding to sub-Saharan African farming populations and those that left Africa diverge next, with subsequent splits of branches for Eurasia, Oceania, East Asia, and America. African hunter-gatherer populations and populations of Oceania and America exhibit no statistically significant signature of growth. The features of population subdivision and growth are discussed in the context of the ancient expansion of modern humans.     

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.     

Differential structuring of human populations for homologous X and Y microsatellite loci.

The global pattern of variation at the homologous microsatellite loci DYS413 (Yq11) and DXS8174 and DXS8175 (Xp22) was analyzed by examination of 30 world populations from four continents, accounting for more than 1,100 chromosomes per locus. The data showed discordant patterns of among- and within-population gene diversity for the Y-linked and the X-linked microsatellites. For the Y-linked polymorphism, all groups of populations displayed high FST values (the correlation between random haplotypes within subpopulations, relative to haplotypes of the total population) and showed a general trend for the haplotypes to cluster in a population-specific way. This was especially true for sub-Saharan African populations. The data also indicated that a large fraction of the variation among populations was due to the accumulation of new variants associated with the radiation process. Europeans exhibited the highest level of within-population haplotype diversity, whereas sub-Saharan Africans showed the lowest. In contrast, data for the two X-linked polymorphisms were concordant in showing lower FST values, as compared with those for DYS413, but higher within-population variances, for African versus non-African populations. Whereas the results for the X-linked loci agreed with a model of greater antiquity for the African populations, those for DYS413 showed a confounding pattern that is apparently at odds with such a model. Possible factors involved in this differential structuring for homologous X and Y microsatellite polymorphisms are discussed.     

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.     

Population Genomics of Sub-Saharan Drosophila melanogaster: African Diversity and Non-African Admixture

Drosophila melanogaster has played a pivotal role in the development of modern population genetics. However, many basic questions regarding the demographic and adaptive history of this species remain unresolved. We report the genome sequencing of 139 wild-derived strains of D. melanogaster, representing 22 population samples from the sub-Saharan ancestral range of this species, along with one European population. Most genomes were sequenced above 25X depth from haploid embryos. Results indicated a pervasive influence of non-African admixture in many African populations, motivating the development and application of a novel admixture detection method. Admixture proportions varied among populations, with greater admixture in urban locations. Admixture levels also varied across the genome, with localized peaks and valleys suggestive of a non-neutral introgression process. Genomes from the same location differed starkly in ancestry, suggesting that isolation mechanisms may exist within African populations. After removing putatively admixed genomic segments, the greatest genetic diversity was observed in southern Africa (e.g. Zambia), while diversity in other populations was largely consistent with a geographic expansion from this potentially ancestral region. The European population showed different levels of diversity reduction on each chromosome arm, and some African populations displayed chromosome arm-specific diversity reductions. Inversions in the European sample were associated with strong elevations in diversity across chromosome arms. Genomic scans were conducted to identify loci that may represent targets of positive selection within an African population, between African populations, and between European and African populations. A disproportionate number of candidate selective sweep regions were located near genes with varied roles in gene regulation. Outliers for Europe-Africa F_ST were found to be enriched in genomic regions of locally elevated cosmopolitan admixture, possibly reflecting a role for some of these loci in driving the introgression of non-African alleles into African populations.     

New founder haplotypes at the myotonic dystrophy locus in southern Africa.

The association between normal alleles at the CTG repeat and two nearby polymorphisms in the myotonin protein kinase gene, the Alu insertion/deletion polymorphism and the myotonic dystrophy kinase (DMK)(G/T) intron 9/HinfI polymorphism, has been analyzed in South African Negroids, a population in which myotonic dystrophy (DM) has not been described. South African Negroids have a CTG allelic distribution that is significantly different from that in Caucasoids and Japanese: the CTG repeat lengths of > or = 19 are very rare. The striking linkage disequilibrium between specific alleles at the Alu polymorphism (Alu(ins) and Alu(del)), the HinfI polymorphism (HinfI-1 and HinfI-2), and the CTG repeat polymorphism seen in Caucasoid (Europeans and Canadians) populations was also found in the South African Negroid population. Numerous haplotypes, not previously described in Europeans, were, however, found. It thus seems likely that only a small number of these "African" chromosomes were present in the progenitors of all non-African peoples. These data provide support for the "out of Africa" model for the origin of modern humans and suggest that the rare ancestral DM mutation event may have occurred after the migration from Africa, hence the absence of DM in sub-Saharan Negroid peoples.     

Geography predicts neutral genetic diversity of human populations

A leading theory for the origin of modern humans, the ‘recent African origin’ (RAO) model [1], postulates that the ancestors of all modern humans originated in East Africa and that, around 100,000 years ago, some modern humans left the African continent and subsequently colonised the entire world, displacing previously established human species such as Neanderthals in Europe 2., 3.. This scenario is supported by the observation that human populations from Africa are genetically the most diverse [2] and that the genetic diversity of non-African populations is negatively correlated with their genetic differentiation towards populations from Africa [3].     

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.     

Chromosomal patterns of microsatellite variability contrast sharply in African and non-African populations of Drosophila melanogaster.

Levels of neutral variation are influenced by background selection and hitchhiking. The relative contribution of these evolutionary forces to the distribution of neutral variation is still the subject of ongoing debates. Using 133 microsatellites, we determined levels of variability on X chromosomes and autosomes in African and non-African D. melanogaster populations. In the ancestral African populations microsatellite variability was higher on X chromosomes than on autosomes. In non-African populations X-linked polymorphism is significantly more reduced than autosomal variation. In non-African populations we observed a significant positive correlation between X chromosomal polymorphism and recombination rate. These results are consistent with the interpretation that background selection shapes levels of neutral variability in the ancestral populations, while the pattern in derived populations is determined by multiple selective sweeps during the colonization process. Further research, however, is required to investigate the influence of inversion polymorphisms and unequal sex ratios.     

Natural selection and molecular evolution in PTC, a bitter-taste receptor gene

The ability to taste phenylthiocarbamide (PTC) is a classic phenotype that has long been known to vary in human populations. This phenotype is of genetic, epidemiologic, and evolutionary interest because the ability to taste PTC is correlated with the ability to taste other bitter substances, many of which are toxic. Thus, variation in PTC perception may reflect variation in dietary preferences throughout human history and could correlate with susceptibility to diet-related diseases in modern populations. To test R. A. Fisher's long-standing hypothesis that variability in PTC perception has been maintained by balancing natural selection, we examined patterns of DNA sequence variation in the recently identified PTC gene, which accounts for up to 85% of phenotypic variance in the trait. We analyzed the entire coding region of PTC (1,002 bp) in a sample of 330 chromosomes collected from African (n=62), Asian (n=138), European (n=110), and North American (n=20) populations by use of new statistical tests for natural selection that take into account the potentially confounding effects of human population growth. Two intermediate-frequency haplotypes corresponding to "taster" and "nontaster" phenotypes were found. These haplotypes had similar frequencies across Africa, Asia, and Europe. Genetic differentiation between the continental population samples was low (FST=0.056) in comparison with estimates based on other genes. In addition, Tajima's D and Fu and Li's D and F statistics demonstrated a significant deviation from neutrality because of an excess of intermediate-frequency variants when human population growth was taken into account (P<.01). These results combine to suggest that balancing natural selection has acted to maintain "taster" and "nontaster" alleles at the PTC locus in humans.