Archaic admixture in the human genome

Estimating bonobo (Pan paniscus) and chimpanzee (Pan troglodytes) evolutionary history from nucleotide site patterns

There is genomic evidence of widespread admixture in deep time between many closely related species, including humans. Our closest living relatives, bonobos and chimpanzees, may also exhibit such patterns. However, assessing the exact degree of interbreeding remains challenging because previous studies have resulted in multiple inconsistent demographic models. We use an approach that addresses these gaps by analyzing all lineages, simultaneously estimating parameters, and comparing previously models. We find evidence of considerable introgression from western into eastern chimpanzees. We also show more breeding females than males and evidence of male-biased dispersal in western chimpanzees. These findings highlight the extent of admixture in bonobo and chimpanzee evolutionary history and are consistent with substantial differences between past and present chimpanzee biogeography.

Admixture appears increasingly ubiquitous in the evolutionary history of various taxa, including humans. Such gene flow likely also occurred among our closest living relatives: bonobos (Pan paniscus) and chimpanzees (Pan troglodytes). However, our understanding of their evolutionary history has been limited by studies that do not consider all Pan lineages or do not analyze all lineages simultaneously, resulting in conflicting demographic models. Here, we investigate this gap in knowledge using nucleotide site patterns calculated from whole-genome sequences from the autosomes of 71 bonobos and chimpanzees, representing all five extant Pan lineages. We estimated demographic parameters and compared all previously proposed demographic models for this clade. We further considered sex bias in Pan evolutionary history by analyzing the site patterns from the X chromosome. We show that 1) 21% of autosomal DNA in eastern chimpanzees derives from western chimpanzee introgression and that 2) all four chimpanzee lineages share a common ancestor about 987,000 y ago, much earlier than previous estimates. In addition, we suggest that 3) there was male reproductive skew throughout Pan evolutionary history and find evidence of 4) male-biased dispersal from western to eastern chimpanzees. Collectively, these results offer insight into bonobo and chimpanzee evolutionary history and suggest considerable differences between current and historic chimpanzee biogeography.

Our Tangled Family Tree: New Genomic Methods Offer Insight into the Legacy of Archaic Admixture

The archaic ancestry present in the human genome has captured the imagination of both scientists and the wider public in recent years. This excitement is the result of new studies pushing the envelope of what we can learn from the archaic genetic information that has survived for over 50,000 years in the human genome. Here, we review the most recent ten years of literature on the topic of archaic introgression, including the current state of knowledge on Neanderthal and Denisovan introgression, as well as introgression from other as-yet unidentified archaic populations. We focus this review on four topics: 1) a reimagining of human demographic history, including evidence for multiple admixture events between modern humans, Neanderthals, Denisovans, and other archaic populations; 2) state-of-the-art methods for detecting archaic ancestry in population-level genomic data; 3) how these novel methods can detect archaic introgression in modern African populations; and 4) the functional consequences of archaic gene variants, including how those variants were co-opted into novel function in modern human populations. The goal of this review is to provide a simple-to-access reference for the relevant methods and novel data, which has changed our understanding of the relationship between our species and its siblings. This body of literature reveals the large degree to which the genetic legacy of these extinct hominins has been integrated into the human populations of today.

Recent effective population size in Eastern European plain Russians correlates with the key historical events

Effective population size reflects the history of population growth, contraction, and structuring. When the effect of structuring is negligible, the inferred trajectory of the effective population size can be informative about the key events in the history of a population. We used the IBDNe and DoRIS approaches, which exploit the data on IBD sharing between genomes, to reconstruct the recent effective population size in two population datasets of Russians from Eastern European plain: (1) ethnic Russians sampled from the westernmost part of Russia; (2) ethnic Russians, Bashkirs, and Tatars sampled from the Volga-Ural region. In this way, we examined changes in effective population size among ethnic Russians that reside in their historical area at the West of the plain, and that expanded eastward to come into contact with the indigenous peoples at the East of the plain. We compared the inferred demographic trajectories of each ethnic group to written historical data related to demographic events such as migration, war, colonization, famine, establishment, and collapse of empires. According to IBDNe estimations, 200 generations (~6000 years) ago, the effective size of the ancestral populations of Russians, Bashkirs, and Tatars hovered around 3,000, 30,000, and 8,000 respectively. Then, the ethnic Russians exponentially grew with increasing rates for the last 115 generations and become the largest ethnic group of the plain. Russians do not show any drop in effective population size after the key historical conflicts, including the Mongol invasion. The only exception is a moderate drop in the 17th century, which is well known in Russian history as The Smuta. Our analyses suggest a more eventful recent population history for the two small ethnic groups that came into contact with ethnic Russians in the Volga-Ural region. We found that the effective population size of Bashkirs and Tatars started to decrease during the time of the Mongol invasion. Interestingly, there is an even stronger drop in the effective population size that coincides with the expansion of Russians to the East. Thus, 15–20 generations ago, i.e. in the 16–18th centuries in the trajectories of Bashkirs and Tatars, we observe the bottlenecks of four and twenty thousand, respectively. Our results on the recent effective population size correlate with the key events in the history of populations of the Eastern European plain and have importance for designing biomedical studies in the region.

Fine scale mapping of genomic introgressions within the Drosophila yakuba clade

The process of speciation involves populations diverging over time until they are genetically and reproductively isolated. Hybridization between nascent species was long thought to directly oppose speciation. However, the amount of interspecific genetic exchange (introgression) mediated by hybridization remains largely unknown, although recent progress in genome sequencing has made measuring introgression more tractable. A natural place to look for individuals with admixed ancestry (indicative of introgression) is in regions where species co-occur. In west Africa, D. santomea and D. yakuba hybridize on the island of São Tomé, while D. yakuba and D. teissieri hybridize on the nearby island of Bioko. In this report, we quantify the genomic extent of introgression between the three species of the Drosophila yakuba clade (D. yakuba, D. santomea), D. teissieri). We sequenced the genomes of 86 individuals from all three species. We also developed and applied a new statistical framework, using a hidden Markov approach, to identify introgression. We found that introgression has occurred between both species pairs but most introgressed segments are small (on the order of a few kilobases). After ruling out the retention of ancestral polymorphism as an explanation for these similar regions, we find that the sizes of introgressed haplotypes indicate that genetic exchange is not recent (>1,000 generations ago). We additionally show that in both cases, introgression was rarer on X chromosomes than on autosomes which is consistent with sex chromosomes playing a large role in reproductive isolation. Even though the two species pairs have stable contemporary hybrid zones, providing the opportunity for ongoing gene flow, our results indicate that genetic exchange between these species is currently rare.

CoMuS: simulating coalescent histories and polymorphic data from multiple species

The simultaneous analysis of intra- and interspecies variation is challenging mainly because our knowledge about patterns of polymorphisms where both intra- and interspecies samples coexist is limited. In this study, we present CoMuS (Coalescent of Multiple Species), a multispecies coalescent software that can simulate intra- and interspecies polymorphisms. CoMuS supports a variety of speciation models and demographic scenarios related to the history of each species. In CoMuS, speciation can be accompanied by either instant or gradual isolation between sister species. Sampling may also occur in the past, and thus, we can study simultaneously extinct and extant species. Our software supports both the infinite- and the finite-site model, with substitution rate heterogeneity among sites and a user-defined proportion of invariable sites. We demonstrate the usage of CoMuS in various applications: species delimitation, software testing, model selection and parameter inference involving present-day and ancestral samples, comparison between gradual and instantaneous isolation models, estimation of speciation time between human and chimpanzee using both intra- and interspecies variation. We expect that CoMuS will be particularly useful for studies where species have been separated recently from their common ancestor and phenomena such as incomplete lineage sorting or introgression still occur.

Model-based analyses of whole-genome data reveal a complex evolutionary history involving archaic introgression in Central African Pygmies

Comparisons of whole-genome sequences from ancient and contemporary samples have pointed to several instances of archaic admixture through interbreeding between the ancestors of modern non-Africans and now extinct hominids such as Neanderthals and Denisovans. One implication of these findings is that some adaptive features in contemporary humans may have entered the population via gene flow with archaic forms in Eurasia. Within Africa, fossil evidence suggests that anatomically modern humans (AMH) and various archaic forms coexisted for much of the last 200,000 yr; however, the absence of ancient DNA in Africa has limited our ability to make a direct comparison between archaic and modern human genomes. Here, we use statistical inference based on high coverage whole-genome data (greater than 60×) from contemporary African Pygmy hunter-gatherers as an alternative means to study the evolutionary history of the genus Homo. Using whole-genome simulations that consider demographic histories that include both isolation and gene flow with neighboring farming populations, our inference method rejects the hypothesis that the ancestors of AMH were genetically isolated in Africa, thus providing the first whole genome-level evidence of African archaic admixture. Our inferences also suggest a complex human evolutionary history in Africa, which involves at least a single admixture event from an unknown archaic population into the ancestors of AMH, likely within the last 30,000 yr.

The peopling of the African continent and the diaspora into the new world

Africa is the birthplace of anatomically modern humans, and is the geographic origin of human migration across the globe within the last 100,000 years. The history of African populations has consisted of a number of demographic events that have influenced patterns of genetic and phenotypic variation across the continent. With the increasing amount of genomic data and corresponding developments in computational methods, researchers are able to explore long-standing evolutionary questions, expanding our understanding of human history within and outside of Africa. This review will summarize some of the recent findings regarding African demographic history, including the African Diaspora, and will briefly explore their implications for disease susceptibility in populations of African descent.

Neanderthal ancestry drives evolution of lipid catabolism in contemporary Europeans

Although Neanderthals are extinct, fragments of their genomes persist in contemporary humans. Here we show that while the genome-wide frequency of Neanderthal-like sites is approximately constant across all contemporary out-of-Africa populations, genes involved in lipid catabolism contain more than threefold excess of such sites in contemporary humans of European descent. Evolutionally, these genes show significant association with signatures of recent positive selection in the contemporary European, but not Asian or African populations. Functionally, the excess of Neanderthal-like sites in lipid catabolism genes can be linked with a greater divergence of lipid concentrations and enzyme expression levels within this pathway, seen in contemporary Europeans, but not in the other populations. We conclude that sequence variants that evolved in Neanderthals may have given a selective advantage to anatomically modern humans that settled in the same geographical areas.

Modern human genomes contain Neanderthal sequences, but it is unclear whether these were selected. Here, Khrameeva et al. show that Neanderthal sequences associated with lipid catabolism are three times more frequent in Europe, suggesting that these sequences might have been beneficial to Europeans.

Study of Modern Human Evolution via Comparative Analysis with the Neanderthal Genome

Many other human species appeared in evolution in the last 6 million years that have not been able to survive to modern times and are broadly known as archaic humans, as opposed to the extant modern humans. It has always been considered fascinating to compare the modern human genome with that of archaic humans to identify modern human-specific sequence variants and figure out those that made modern humans different from their predecessors or cousin species. Neanderthals are the latest humans to become extinct, and many factors made them the best representatives of archaic humans. Even though a number of comparisons have been made sporadically between Neanderthals and modern humans, mostly following a candidate gene approach, the major breakthrough took place with the sequencing of the Neanderthal genome. The initial genome-wide comparison, based on the first draft of the Neanderthal genome, has generated some interesting inferences regarding variations in functional elements that are not shared by the two species and the debated admixture question. However, there are certain other genetic elements that were not included or included at a smaller scale in those studies, and they should be compared comprehensively to better understand the molecular make-up of modern humans and their phenotypic characteristics. Besides briefly discussing the important outcomes of the comparative analyses made so far between modern humans and Neanderthals, we propose that future comparative studies may include retrotransposons, pseudogenes, and conserved non-coding regions, all of which might have played significant roles during the evolution of modern humans.

Length distributions of identity by descent reveal fine-scale demographic history

Data-driven studies of identity by descent (IBD) were recently enabled by high-resolution genomic data from large cohorts and scalable algorithms for IBD detection. Yet, haplotype sharing currently represents an underutilized source of information for population-genetics research. We present analytical results on the relationship between haplotype sharing across purportedly unrelated individuals and a population's demographic history. We express the distribution of IBD sharing across pairs of individuals for segments of arbitrary length as a function of the population's demography, and we derive an inference procedure to reconstruct such demographic history. The accuracy of the proposed reconstruction methodology was extensively tested on simulated data. We applied this methodology to two densely typed data sets: 500 Ashkenazi Jewish (AJ) individuals and 56 Kenyan Maasai (MKK) individuals (HapMap 3 data set). Reconstructing the demographic history of the AJ cohort, we recovered two subsequent population expansions, separated by a severe founder event, consistent with previous analysis of lower-throughput genetic data and historical accounts of AJ history. In the MKK cohort, high levels of cryptic relatedness were detected. The spectrum of IBD sharing is consistent with a demographic model in which several small-sized demes intermix through high migration rates and result in enrichment of shared long-range haplotypes. This scenario of historically structured demographies might explain the unexpected abundance of runs of homozygosity within several populations.

Inference of human population history from individual whole-genome sequences

The history of human population size is important for understanding human evolution. Various studies have found evidence for a founder event (bottleneck) in East Asian and European populations, associated with the human dispersal out-of-Africa event around 60 thousand years (kyr) ago. However, these studies have had to assume simplified demographic models with few parameters, and they do not provide a precise date for the start and stop times of the bottleneck. Here, with fewer assumptions on population size changes, we present a more detailed history of human population sizes between approximately ten thousand and a million years ago, using the pairwise sequentially Markovian coalescent model applied to the complete diploid genome sequences of a Chinese male (YH), a Korean male (SJK), three European individuals (J. C. Venter, NA12891 and NA12878 (ref. 9)) and two Yoruba males (NA18507 (ref. 10) and NA19239). We infer that European and Chinese populations had very similar population-size histories before 10-20 kyr ago. Both populations experienced a severe bottleneck 10-60 kyr ago, whereas African populations experienced a milder bottleneck from which they recovered earlier. All three populations have an elevated effective population size between 60 and 250 kyr ago, possibly due to population substructure. We also infer that the differentiation of genetically modern humans may have started as early as 100-120 kyr ago, but considerable genetic exchanges may still have occurred until 20-40 kyr ago.

Evolutionary genetics evidence of an essential, nonredundant role of the IFN-γ pathway in protective immunity

Identifying how natural selection has affected immunity-related genes can provide insights into the mechanisms that have been crucial for our survival against infection. Rare disorders of either chain of the IFN-γ receptor, but not of IFN-γ itself, have been shown to confer predisposition to mycobacterial disease in patients otherwise normally resistant to most viruses. Here, we defined the levels of naturally occurring variation in the three specific genes controlling the IFN-γ pathway (IFNG, IFNGR1, IFNGR2) and assessed whether and how natural selection has acted on them. To this end, we resequenced the three genes in 186 individuals from sub-Saharan Africa, Europe, and East-Asia. Our results show that IFNG is subject to strong purifying selection against nonsynonymous variants. Conversely, IFNGR1 and IFNGR2 evolve under more relaxed selective constraints, although they are not completely free to accumulate amino acid variation having a major impact on protein function. In addition, we have identified signatures of population-specific positive selection, including at one intronic variant known to be associated with higher production of IFN-γ. The integration of our population genetic data into a clinical framework demonstrates that the IFN-γ pathway is essential and nonredundant in host defense, probably because of its role in protective immunity against mycobacteria.

Testing for ancient admixture between closely related populations

One enduring question in evolutionary biology is the extent of archaic admixture in the genomes of present-day populations. In this paper, we present a test for ancient admixture that exploits the asymmetry in the frequencies of the two nonconcordant gene trees in a three-population tree. This test was first applied to detect interbreeding between Neandertals and modern humans. We derive the analytic expectation of a test statistic, called the D statistic, which is sensitive to asymmetry under alternative demographic scenarios. We show that the D statistic is insensitive to some demographic assumptions such as ancestral population sizes and requires only the assumption that the ancestral populations were randomly mating. An important aspect of D statistics is that they can be used to detect archaic admixture even when no archaic sample is available. We explore the effect of sequencing error on the false-positive rate of the test for admixture, and we show how to estimate the proportion of archaic ancestry in the genomes of present-day populations. We also investigate a model of subdivision in ancestral populations that can result in D statistics that indicate recent admixture.

Neandertal genome: the ins and outs of African genetic diversity

Analysis of the Neandertal genome indicates gene flow between Neandertals and modern humans of Eurasia but not Africa. This surprising result is difficult to reconcile with current models of human origins and might have to do with insufficient African sampling.

A draft sequence of the Neandertal genome

Neandertals, the closest evolutionary relatives of present-day humans, lived in large parts of Europe and western Asia before disappearing 30,000 years ago. We present a draft sequence of the Neandertal genome composed of more than 4 billion nucleotides from three individuals. Comparisons of the Neandertal genome to the genomes of five present-day humans from different parts of the world identify a number of genomic regions that may have been affected by positive selection in ancestral modern humans, including genes involved in metabolism and in cognitive and skeletal development. We show that Neandertals shared more genetic variants with present-day humans in Eurasia than with present-day humans in sub-Saharan Africa, suggesting that gene flow from Neandertals into the ancestors of non-Africans occurred before the divergence of Eurasian groups from each other.

Formulating a historical and demographic model of recent human evolution based on resequencing data from noncoding regions

BACKGROUND

Estimating the historical and demographic parameters that characterize modern human populations is a fundamental part of reconstructing the recent history of our species. In addition, the development of a model of human evolution that can best explain neutral genetic diversity is required to identify confidently regions of the human genome that have been targeted by natural selection.

METHODOLOGY/PRINCIPAL FINDINGS

We have resequenced 20 independent noncoding autosomal regions dispersed throughout the genome in 213 individuals from different continental populations, corresponding to a total of approximately 6 Mb of diploid resequencing data. We used these data to explore and co-estimate an extensive range of historical and demographic parameters with a statistical framework that combines the evaluation of multiple models of human evolution via a best-fit approach, followed by an Approximate Bayesian Computation (ABC) analysis. From a methodological standpoint, evaluating the accuracy of the parameter co-estimation allowed us to identify the most accurate set of statistics to be used for the estimation of each of the different historical and demographic parameters characterizing recent human evolution.

CONCLUSIONS/SIGNIFICANCE

Our results support a model in which modern humans left Africa through a single major dispersal event occurring approximately 60,000 years ago, corresponding to a drastic reduction of approximately 5 times the effective population size of the ancestral African population of approximately 13,800 individuals. Subsequently, the ancestors of modern Europeans and East Asians diverged much later, approximately 22,500 years ago, from the population of ancestral migrants. This late diversification of Eurasians after the African exodus points to the occurrence of a long maturation phase in which the ancestral Eurasian population was not yet diversified.

Detecting ancient admixture and estimating demographic parameters in multiple human populations

We analyze patterns of genetic variation in extant human polymorphism data from the National Institute of Environmental Health Sciences single nucleotide polymorphism project to estimate human demographic parameters. We update our previous work by considering a larger data set (more genes and more populations) and by explicitly estimating the amount of putative admixture between modern humans and archaic human groups (e.g., Neandertals, Homo erectus, and Homo floresiensis). We find evidence for this ancient admixture in European, East Asian, and West African samples, suggesting that admixture between diverged hominin groups may be a general feature of recent human evolution.

Positive selection in the human genome: from genome scans to biological significance

Here we review the evidence for positive selection in the human genome and its role in human evolution and population differentiation. In recent years, there has been a dramatic increase in the use of genome-wide scans to identify adaptively evolving loci in the human genome. Attention is now turning to understanding the biological relevance and adaptive significance of the regions identified as being subject to recent positive selection. Examples of adaptively evolving loci are discussed, specifically LCT and FOXP2. Comprehensive studies of these loci also provide information about the functional relevance of the selected alleles. We discuss current studies examining the role of positive selection in shaping copy number variation and noncoding genomic regions and highlight challenges presented by the study of positive selection in the human genome.

Subdivision in an ancestral species creates asymmetry in gene trees

We consider gene trees in three species for which the species tree is known. We show that population subdivision in ancestral species can lead to asymmetry in the frequencies of the two gene trees not concordant with the species tree and, if subdivision is extreme, cause the one of the nonconcordant gene trees to be more probable than the concordant gene tree. Although published data for the human-chimp-gorilla clade and for three species of Drosophila show asymmetry consistent with our model, sequencing error could also account for observed patterns. We show that substantial levels of persistent ancestral subdivision are needed to account for the observed levels of asymmetry found in these two studies.

Testing for archaic hominin admixture on the X chromosome: model likelihoods for the modern human RRM2P4 region from summaries of genealogical topology under the structured coalescent

A 2.4-kb stretch within the RRM2P4 region of the X chromosome, previously sequenced in a sample of 41 globally distributed humans, displayed both an ancient time to the most recent common ancestor (e.g., a TMRCA of approximately 2 million years) and a basal clade composed entirely of Asian sequences. This pattern was interpreted to reflect a history of introgressive hybridization from archaic hominins (most likely Asian Homo erectus) into the anatomically modern human genome. Here, we address this hypothesis by resequencing the 2.4-kb RRM2P4 region in 131 African and 122 non-African individuals and by extending the length of sequence in a window of 16.5 kb encompassing the RRM2P4 pseudogene in a subset of 90 individuals. We find that both the ancient TMRCA and the skew in non-African representation in one of the basal clades are essentially limited to the central 2.4-kb region. We define a new summary statistic called the minimum clade proportion (pmc), which quantifies the proportion of individuals from a specified geographic region in each of the two basal clades of a binary gene tree, and then employ coalescent simulations to assess the likelihood of the observed central RRM2P4 genealogy under two alternative views of human evolutionary history: recent African replacement (RAR) and archaic admixture (AA). A molecular-clock-based TMRCA estimate of 2.33 million years is a statistical outlier under the RAR model; however, the large variance associated with this estimate makes it difficult to distinguish the predictions of the human origins models tested here. The pmc summary statistic, which has improved power with larger samples of chromosomes, yields values that are significantly unlikely under the RAR model and fit expectations better under a range of archaic admixture scenarios.

No evidence of a Neanderthal contribution to modern human diversity

Recent studies of genomic DNA from Neanderthal fossils support earlier studies showing no contribution of Neanderthals to the gene pool of modern humans.

The relationship between Neanderthals and modern humans is contentious, but recent advances in Neanderthal genomics have shed new light on their evolutionary history. Here we review the available evidence and find no indication of any Neanderthal contribution to modern genetic diversity.

African genetic diversity: implications for human demographic history, modern human origins, and complex disease mapping

Comparative studies of ethnically diverse human populations, particularly in Africa, are important for reconstructing human evolutionary history and for understanding the genetic basis of phenotypic adaptation and complex disease. African populations are characterized by greater levels of genetic diversity, extensive population substructure, and less linkage disequilibrium (LD) among loci compared to non-African populations. Africans also possess a number of genetic adaptations that have evolved in response to diverse climates and diets, as well as exposure to infectious disease. This review summarizes patterns and the evolutionary origins of genetic diversity present in African populations, as well as their implications for the mapping of complex traits, including disease susceptibility.