Admixture in Southern Africa (Petersen et al. 2013)

Related:

• The genetic prehistory of southern Africa
• Genomic Variation in Seven Khoe-San Groups Reveals Adaptation and Complex African History

PLoS Genet 9(3): e1003309. doi:10.1371/journal.pgen.1003309
Complex Patterns of Genomic Admixture within Southern Africa
Desiree C. Petersen et al.
Within-population genetic diversity is greatest within Africa, while between-population genetic diversity is directly proportional to geographic distance. The most divergent contemporary human populations include the click-speaking forager peoples of southern Africa, broadly defined as Khoesan. Both intra- (Bantu expansion) and inter-continental migration (European-driven colonization) have resulted in complex patterns of admixture between ancient geographically isolated Khoesan and more recently diverged populations. Using gender-specific analysis and almost 1 million autosomal markers, we determine the significance of estimated ancestral contributions that have shaped five contemporary southern African populations in a cohort of 103 individuals. Limited by lack of available data for homogenous Khoesan representation, we identify the Ju/'hoan (n = 19) as a distinct early diverging human lineage with little to no significant non-Khoesan contribution. In contrast to the Ju/'hoan, we identify ancient signatures of Khoesan and Bantu unions resulting in significant Khoesan- and Bantu-derived contributions to the Southern Bantu amaXhosa (n = 15) and Khoesan !Xun (n = 14), respectively. Our data further suggests that contemporary !Xun represent distinct Khoesan prehistories. Khoesan assimilation with European settlement at the most southern tip of Africa resulted in significant ancestral Khoesan contributions to the Coloured (n = 25) and Baster (n = 30) populations. The latter populations were further impacted by 170 years of East Indian slave trade and intra-continental migrations resulting in a complex pattern of genetic variation (admixture). The populations of southern Africa provide a unique opportunity to investigate the genomic variability from some of the oldest human lineages to the implications of complex admixture patterns including ancient and recently diverged human lineages.
Link

Extremely old (237–581 kya) root of human Y-chromosome phylogeny

I had mentioned this research before, and now it has officially been published. There are three things to be excited about this new paper:
First, it forces us to consider the possibility of either (i) archaic admixture in Africa, or (ii) a much more ancient time depth of modern humans than the first fossils from Ethiopia dated to about ~200 thousand years ago.
Second, it underscores the importance of collaboration between academia and regular folk, since it was the combined contributions of academics, genetic genealogists, and the owners of the new A00 basal Y-chromosomes that made this new discovery possible.
And, third, it shows that the extraordinary can be discovered without contacting isolated tribes or seeking human bones in remote regions, but rather through careful scrutiny of large volumes of data for the proverbial needle in the haystack.
The paper developed a model of Y-chromosome mutation based on the estimate of Kong et al. Significantly, though:

If we were to use the higher mutation rate (1.0x10^-9 per base per year6) rather than a realistic range derived from whole-genome sequencing (4.39x10^-9  -   7.07x10^- 9), the estimated TMRCA for the tree incorporating A00 as the basal lineage would be 209 kya, which is only slightly older than current estimates of the TMRCA of mtDNA and the age of the oldest AMH fossil remains. We note, however, that the higher mutation rate produces an estimate for the common ancestor of all non-African Y chromosome haplogroups (C through T) of ~39 kya6 (i.e., versus ~63 kya for the mutation rate used here).

A 39kya common ancestor for Eurasian Y-chromosomes makes no sense, since we now know for sure that by that time, the differentiation of Eurasians was already well on its way and modern humans in remote parts of the Old World have been documented much earlier than that time.
A ~63kya common ancestor, on the other hand, fits nicely with my "two deserts" theory of modern human origins, according to which the ancestors of Eurasians faced an ecological crisis in Arabia when it became much drier post-70kya; that seems like a most opportune time for the major Eurasian bottleneck and the corresponding coalescence of Eurasian Y-chromosomes to a single man. And, while there is no a priori reason for Y chromosomes and mtDNA to behave similarly, the age of the "older" Eurasian ur-mother, haplogroup N at 59 thousand years, with presumably an older ancestor within mtDNA haplogroup L3 founding the Eurasian population.
Also, if modern human-Neandertal admixture had occurred  "most likely 47,000–65,000 years ago", then the expansion of modern Proto-Eurasians within a 70-60kya timeframe north Out-of-Arabia would have brought them in contact with their northern Neandertal neighbors. On the other hand, it would be incredible if modern humans experienced admixture with Neandertals but were still much later a very small population (to allow for the coalescence of their Y-chromosomes to one man ~39kya).
So, in summary, the mutation rate used by the authors seems consistent with what we know about an important calibration point of the human story.
But, who were the people in Africa responsible for the introgression of A00 chromosomes? Mendez et al. used the haplotype of the African American A00 individual and discovered his patrilineal kin among theMbo of Cameroon, who are Bantu farmers.
I have observed before that Pygmies and Bushmen represent only a tiny fraction of pre-existing African genetic diversity, the part that had not yet been absorbed into the farmers' expanding population by the time that Africa came to the attention of of modern science. We see traces of Pygmy and Bushman ancestry in some African farmers, and there were probably other groups, no longer extant as distinct ethno-cultural entities, but, nonetheless, surviving as genetic fragments in the genomes of the farmers.
Thus, while it still makes sense to study the surviving hunter-gatherers of Africa who make up perhaps a percent or less of the population of Africa, it may be equally important to study different groups of African farmers who may possess a much richer treasure trove (albeit diluted) of such "Palaeoafrican" ancestry.
Finally:

Although the stochastic nature of the evolutionary process can explain the aforementioned incongruences, the extreme age and rarity of the A00 lineage point to the possibility of a highly structured ancestral population, consistent with recent work on the autosomes.40,41,43,44 This could take the form of long-standing population structure among AMH populations45 or archaic introgression from an archaic form into the ancestors of AMHs.46 Interestingly, the Mbo live less than 800 km away from a Nigerian site known as Iwo Eleru, where human skeletal remains with both archaic and modern features were found and dated to ~13 kya.47 Further surveys in sub- Saharan Africa and in the African Diaspora might uncover more diverged basal lineages, which will help to disentangle some of the complex evolutionary processes that shape patterns of Y chromosome diversity.

AJHG 10.1016/j.ajhg.2013.02.002
An African American Paternal Lineage Adds an Extremely Ancient Root to the Human Y Chromosome Phylogenetic Tree
Fernando L. Mendez et al.
We report the discovery of an African American Y chromosome that carries the ancestral state of all SNPs that defined the basal portion of the Y chromosome phylogenetic tree. We sequenced ∼240 kb of this chromosome to identify private, derived mutations on this lineage, which we named A00. We then estimated the time to the most recent common ancestor (TMRCA) for the Y tree as 338 thousand years ago (kya) (95% confidence interval = 237–581 kya). Remarkably, this exceeds current estimates of the mtDNA TMRCA, as well as those of the age of the oldest anatomically modern human fossils. The extremely ancient age combined with the rarity of the A00 lineage, which we also find at very low frequency in central Africa, point to the importance of considering more complex models for the origin of Y chromosome diversity. These models include ancient population structure and the possibility of archaic introgression of Y chromosomes into anatomically modern humans. The A00 lineage was discovered in a large database of consumer samples of African Americans and has not been identified in traditional hunter-gatherer populations from sub-Saharan Africa. This underscores how the stochastic nature of the genealogical process can affect inference from a single locus and warrants caution during the interpretation of the geographic location of divergent branches of the Y chromosome phylogenetic tree for the elucidation of human origins.
Link

Deep mtDNA substructure in southern Africa (Barbieri et al. 2013)

Deep mtDNA substructure in southern Africa (Barbieri et al. 2013)

The Khoisan have been used in many different ways in reconstructions of human history.

Being probably the most genetically diverse modern human population, they are occasionally viewed as akin to the ur-humans, with everyone else shedding diversity via founder effects as they moved away from a south African modern human urheimat.

They are also sometimes viewed as a basal branch of the human family tree, and they probably are -if modern humans are made to fit a tree model. But, modern humans didn't really evolve tree-like (some African farmers have Khoisan-like admixture, and the Khoisan themselves have relatively "shallow" common ancestry with other Africans and many Eurasians on account of their possession of a respectable frequency of Y-haplogroup E).

I have sometimes noted that in the case of South African groups were are lucky that the Khoisan exist as a discrete set of populations, making it easier to discern the legacy of South African hunter-gatherers in the genomes of immigrant farmers and pastoralists who converged southwards over the last few thousand years. This can be contrasted with the presumable situation in places like West Africa (the cradle of Sub-Saharan African farming), in which any indigenous hunter-gatherer groups have ceased to exist as distinct entities a long time ago.

A new AJHG paper sample south African genomes extensively and arrives at a startling conclusion. In the words of the authors:

Overall, the results of this analysis indicate that it is very unlikely that the highly divergent L0k1b/L0k2 lineages were incorporated into the Bantu-speaking populations via gene flow from a population that was ancestral to a Khoisan population in our sample but subsequently lost from the Khoisan population via drift. Instead, these results support the hypothesis that the ancestors of the Bantu-speaking populations carrying the divergent L0k lineages (who now live mainly in Zambia) experienced gene flow from a pre-Bantu population that is nowadays extinct. Alternatively, it is possible that descendants from this pre-Bantu population do exist but have not yet been included in population genetic studies; however, our extensive sampling of populations from Botswana, Namibia, andWest Zambia (which includes representatives of nearly all known Khoisan groups) makes it highly unlikely that this pre-Bantu Khoisan population has not yet been sampled.

In other words, we must resist the tendency to think of the Khoisan as representatives of all pre-Bantu south Africans. The Khoisan are certainly descendants of old south Africans, and represent a part of the pre-Bantu genetic landscape that retained its cultural distinctiveness (and hence can be nowadays sampled as a distinct population). But, there were other, now submerged, peaks in that landscape that are no longer extant in distinct form, but only in absorbed form in the gene pool of south African farmers.

This is fairly interesting in itself, and certainly ought to change our belief about what Africa looked like pre-Bantu expansion. We ought to think of, perhaps, a cornucopia of groups: many of them may have gone extinct; some may have been completely absorbed into more successful ones, and perhaps only a handful survive as distinct entities. Such a view would agree with the conclusions of physical anthropology about the persistence of archaic-leaning groups in parts of Africa down to the Holocene boundary.

The American Journal of Human Genetics, 17 January 2013 doi:10.1016/j.ajhg.2012.12.010

Ancient Substructure in Early mtDNA Lineages of Southern Africa

Chiara Barbieri et al.


Among the deepest-rooting clades in the human mitochondrial DNA (mtDNA) phylogeny are the haplogroups defined as L0d and L0k, which are found primarily in southern Africa. These lineages are typically present at high frequency in the so-called Khoisan populations of hunter-gatherers and herders who speak non-Bantu languages, and the early divergence of these lineages led to the hypothesis of ancient genetic substructure in Africa. Here we update the phylogeny of the basal haplogroups L0d and L0k with 500 full mtDNA genome sequences from 45 southern African Khoisan and Bantu-speaking populations. We find previously unreported subhaplogroups and greatly extend the amount of variation and time-depth of most of the known subhaplogroups. Our major finding is the definition of two ancient sublineages of L0k (L0k1b and L0k2) that are present almost exclusively in Bantu-speaking populations from Zambia; the presence of such relic haplogroups in Bantu speakers is most probably due to contact with ancestral pre-Bantu populations that harbored different lineages than those found in extant Khoisan. We suggest that although these populations went extinct after the immigration of the Bantu-speaking populations, some traces of their haplogroup composition survived through incorporation into the gene pool of the immigrants. Our findings thus provide evidence for deep genetic substructure in southern Africa prior to the Bantu expansion that is not represented in extant Khoisan populations.


Link

Ethiopian origins (Pagani et al. 2012)

Ethiopian origins (Pagani et al. 2012)

The study attempts to answer four questions:

Our current study is motivated by four questions. First, where do the Ethiopians stand in the African genetic landscape? Second, what is the extent of recent gene flow from outside Africa into Ethiopia, when did it occur, and is there evidence of selection effects? Third, do genomic data support a route for out-of-Africa migration of modern humans across the mouth of the Red Sea? Fourth, assuming temporal stability of current populations, what are the estimated ages of Ethiopian populations relative to other African groups?

Link to press release. Link the supplemental data.

The authors reiterate that modern humans left Africa 50-70kya, a hypothesis that seems to me pretty much dead in the light of recent archaeological evidence.

The lack of antiquity in the Ethiopian population, even in only the African component thereof argues against that population being ancestral to modern humans. Note that if the Out-of-East Africa hypothesis is correct, then skulls like Omo I represent ancestral modern humans and they are followed much later by modern humans anywhere else. So, while anatomical modernity may have emerged in East Africa --or maybe not; let's not forget that we have early modern skulls from the region in part because of the excellent preservation conditions and excess of scholarly interest-- there is no evidence that they spread from there.

I have little doubt that my own theory about substantial back-migration of Eurasians into Africa will eventually win the day. Of course, I am not referring to the recent (in the last 3,000 years) admixture with West Eurasians that the Ethiopian population has undergone, but rather to the more ancient migration that was probably associated with Y-haplogroup DE-YAP.

The fact that the African component of diverse African populations is more closely related to West than to East Eurasians is one piece of evidence among many for that scenario. Hopefully, it can be tested soon using whole genome data which may have enough density to detect much older admixture events.

UPDATE I: Since the dates in the paper are based on ROLLOFF, a piece of software that is not publicly available more than a year after its announcement, and which contradicts other software released by the same authors, I will take the Queen of Sheba stories circulated in the media with a huge grain of salt.

The American Journal of Human Genetics, 21 June 2012 doi:10.1016/j.ajhg.2012.05.015

Ethiopian Genetic Diversity Reveals Linguistic Stratification and Complex Influences on the Ethiopian Gene Pool

Luca Pagani et al.

Humans and their ancestors have traversed the Ethiopian landscape for millions of years, and present-day Ethiopians show great cultural, linguistic, and historical diversity, which makes them essential for understanding African variability and human origins. We genotyped 235 individuals from ten Ethiopian and two neighboring (South Sudanese and Somali) populations on an Illumina Omni 1M chip. Genotypes were compared with published data from several African and non-African populations. Principal-component and STRUCTURE-like analyses confirmed substantial genetic diversity both within and between populations, and revealed a match between genetic data and linguistic affiliation. Using comparisons with African and non-African reference samples in 40-SNP genomic windows, we identified “African” and “non-African” haplotypic components for each Ethiopian individual. The non-African component, which includes the SLC24A5 allele associated with light skin pigmentation in Europeans, may represent gene flow into Africa, which we estimate to have occurred ∼3 thousand years ago (kya). The African component was found to be more similar to populations inhabiting the Levant rather than the Arabian Peninsula, but the principal route for the expansion out of Africa ∼60 kya remains unresolved. Linkage-disequilibrium decay with genomic distance was less rapid in both the whole genome and the African component than in southern African samples, suggesting a less ancient history for Ethiopian populations.

Link

Introgression of archaic haplotype at OAS1 in Melanesians (Mendez et al. 2012)

Introgression of archaic haplotype at OAS1 in Melanesians (Mendez et al. 2012)

It seems that Michael Hammer was good on his promise that in 2012 "This year, we should be able to confirm what we found and go way beyond that."  In a new paper, conclusive evidence is presented about introgression of an archaic sequence into Melanesian populations. The argument is as follows:

• Melanesians are more diverse in that region than Africans.
• The common ancestor of the "archaic" and "African" haplotypes lived >3 million years ago.
• The "archaic" haplotype matches the ancient DNA from the Denisova hominin.
• Balancing selection (which can sometimes maintain extremely old polymorphism) is not reasonable in this case, because it would need to maintain both "archaic" and "African" haplotypes for a long time, but then (inexplicably) would continue to operate in Melanesia and cease to operate everywhere else.

Notice that once again, this is based on resequencing a small region of the genome. This is why I am all the more confident in my prediction that the advent of full genome sequencing will uncover more archaic admixture in humans. It may not always be able to use all the above listed criteria to confirm this admixture (since we do not and cannot have ancient DNA from all the archaic hominins that once roamed the planet), but all the remaining ones will suffice to make a very good case for introgression.

What I find particularly interesting, is that Mendez et al. re-iterate a few times that genomewide averages admit to different explanations:

Full genome comparisons of the Neandertal and Denisova draft genomes with modern human sequences have revealed different amounts of shared ancestry between each of these archaic forms and anatomically modern human (AMH) populations from different geographic regions. For example, a higher proportion of SNPs was shared between non-African and Neandertal, and between Melanesian and the Denisova genomes, than between either Neandertal or Denisova and extant African genomes (Green et al. 2010; Reich et al. 2010). An intriguing possibility is that these patterns result from introgression of archaic genes into AMH populations in Eurasia. However, this SNP sharing pattern could also be explained by ancestral population structure in Africa (i.e., without the need to posit introgression). For example, if non-Africans and the ancestors of Neandertals descend from the same deme in a subdivided African population, and this structure persisted with low levels of gene flow among African residents until the ancestors of non-Africans migrated into Eurasia, then we would expect more SNP sharing between non-Africans and Neandertals (Durand et al. 2011). 

... 

While genome-wide comparisons detect more sequence agreement between non-African and Neandertal genomes, and between Melanesian and Denisova genomes, the specific loci exhibiting these signals have not yet been identified. Furthermore, current analyses do not elucidate the relative roles of recent introgression versus long-term population structure in Africa in explaining these patterns.

The current paper does a good job at showing how in one particular region archaic introgression into Melanesians is indeed the best explanation for the evidence. But, the fact that the authors seem to re-iterate the possibility of African population structure and repeatedly caution against using patterns of genomewide sharing between modern and archaic humans is a strong hint that there are more things to come on the topic.

We should remember that the widely-circulated estimates of Neandertal->Eurasian introgression are based on genomewide averages. It is true that Reich et al. (2010) identified 13 regions of potential Neandertal introgression, which together make up a very small portion of the human genome. So, the jury is out on whether African population structure or Neandertal introgression is responsible for most of the genomewide pattern.

What you can be sure of is that many scientists are busy lining up full genomes from different human populations as we speak, and finding plenty of regions where haplotypes of extremely old divergence times co-exist in our species. We will probably learn more about such efforts during 2012.



Mol Biol Evol (2012)doi: 10.1093/molbev/msr301

Global genetic variation at OAS1 provides evidence of archaic admixture in Melanesian populations

Fernando L. Mendez, Joseph C. Watkins and Michael F. Hammer

Recent analysis of DNA extracted from two Eurasian forms of archaic human show that more genetic variants are shared with humans currently living in Eurasia than with anatomically modern humans in sub-Saharan Africa. While these genome-wide average measures of genetic similarity are consistent with the hypothesis of archaic admixture in Eurasia, analyses of individual loci exhibiting the signal of archaic introgression are needed to test alternative hypotheses and investigate the admixture process. Here, we provide a detailed sequence analysis of the innate immune gene, OAS1, a locus with a divergent Melanesian haplotype that is very similar to the Denisova sequence from the Altai region of Siberia. We re-sequenced a 7 kb region encompassing the OAS1 gene in 88 individuals from 6 Old World populations (San, Biaka, Mandenka, French Basque, Han Chinese, and Papua New Guineans) and discovered previously unknown and ancient genetic variation. The 5' region of this gene has unusual patterns of diversity, including 1) higher levels of nucleotide diversity in Papuans than in sub-Saharan Africans, 2) very deep ancestry with an estimated time to the most recent common ancestor of >3 million years, and 3) a basal branching pattern with Papuan individuals on either side of the rooted network. A global geographic survey of >1500 individuals showed that the divergent Papuan haplotype is nearly restricted to populations from eastern Indonesia and Melanesia. Polymorphic sites within this haplotype are shared with the draft Denisova genome over a span of ∼90 kb and are associated with an extended block of linkage disequilibrium, supporting the hypothesis that this haplotype introgressed from an archaic source that likely lived in Eurasia.

Link