Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Neanderthal introgression reintroduced functional ancestral alleles lost in Eurasian populations

Abstract

Neanderthal ancestry remains across modern Eurasian genomes and introgressed sequences influence diverse phenotypes. Here, we demonstrate that introgressed sequences reintroduced thousands of ancestral alleles that were lost in Eurasian populations before introgression. Our simulations and variant effect predictions argue that these reintroduced alleles (RAs) are more likely to be tolerated by modern humans than are introgressed Neanderthal-derived alleles (NDAs) due to their distinct evolutionary histories. Consistent with this, we show enrichment for RAs and depletion for NDAs on introgressed haplotypes with expression quantitative trait loci (eQTL) and phenotype associations. Analysis of available cross-population eQTLs and massively parallel reporter assay data show that RAs commonly influence gene expression independent of linked NDAs. We further validate these independent effects for one RA in vitro. Finally, we demonstrate that NDAs are depleted for regulatory activity compared to RAs, while RAs have activity levels similar to non-introgressed variants. In summary, our study reveals that Neanderthal introgression reintroduced thousands of lost ancestral variants with gene regulatory activity and that these RAs were more tolerated than NDAs. Thus, RAs and their distinct evolutionary histories must be considered when evaluating the effects of introgression.

This is a preview of subscription content, access via your institution

Access options

Rent or buy this article

Prices vary by article type

from$1.95

to$39.95

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Schematic of the reintroduction of lost ancestral alleles by Neanderthal introgression.
Fig. 2: Neanderthal introgression reintroduced thousands of lost ancestral alleles to Eurasian populations.
Fig. 3: RAs are more prevalent than NDAs among introgressed haplotypes with GWAS hits and eQTL.
Fig. 4: RAs restore regulatory functions lost in Eurasians.
Fig. 5: Hundreds of RAs have regulatory activity independent of NDAs.

Similar content being viewed by others

Data availability

All results reported in this paper are available in supplementary material and/or on the project’s github repository (https://github.com/DaRinker/Neanderthal.reintroduction).

Code availability

Full code used in this analysis are available on the project’s github (https://github.com/DaRinker/Neanderthal.reintroduction).

References

  1. The 1000 Genomes Project Consortium. A global reference for human genetic variation. Nature 526, 68–74 (2015).

    PubMed Central  Google Scholar 

  2. Mallick, S. et al. The Simons genome diversity project: 300 genomes from 142 diverse populations. Nature 538, 201–206 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  3. Pagani, L. et al. Genomic analyses inform on migration events during the peopling of Eurasia. Nature 538, 238–242 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Henn, B. M., Botigué, L. R., Bustamante, C. D., Clark, A. G. & Gravel, S. Estimating mutation load in human genomes. Nat. Rev. Genet. 16, 333–343 (2015).

    CAS  PubMed  PubMed Central  Google Scholar 

  5. Prüfer, K. et al. The complete genome sequence of a Neanderthal from the Altai mountains. Nature 505, 43–49 (2014).

    PubMed  Google Scholar 

  6. Prüfer, K. et al. A high-coverage Neandertal genome from Vindija cave in Croatia. Science 358, 655–658 (2017).

    PubMed  PubMed Central  Google Scholar 

  7. Meyer, M. et al. A high-coverage genome sequence from an archaic Denisovan individual. Science 338, 222–226 (2012).

    CAS  PubMed  PubMed Central  Google Scholar 

  8. Green, R. E. et al. A draft sequence of the Neandertal genome. Science 328, 710–722 (2010).

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Sankararaman, S., Mallick, S., Patterson, N. & Reich, D. The combined landscape of Denisovan and Neanderthal ancestry in present-day humans. Curr. Biol. 26, 1241–1247 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Sankararaman, S. et al. The genomic landscape of Neanderthal ancestry in present-day humans. Nature 507, 354–357 (2014).

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Vernot, B. & Akey, J. M. Resurrecting surviving Neandertal lineages from modern human genomes. Science 343, 1017–1021 (2014).

    CAS  PubMed  Google Scholar 

  12. Abi-Rached, L. et al. The shaping of modern human immune systems by multiregional admixture with archaic humans. Science 334, 89–94 (2011).

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Mendez, F. L., Watkins, J. C. & Hammer, M. F. A Haplotype at STAT2 introgressed from Neanderthals and serves as a candidate of positive selection in Papua New Guinea. Am. J. Hum. Genet. 91, 265–274 (2012).

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Dannemann, M., Prüfer, K. & Kelso, J. Functional implications of Neandertal introgression in modern humans. Genome Biol. 18, 61 (2017).

    PubMed  PubMed Central  Google Scholar 

  15. Racimo, F., Marnetto, D. & Huerta-Sánchez, E. Signatures of archaic adaptive introgression in present-day human populations. Mol. Biol. Evol. 34, 296–317 (2017).

  16. Racimo, F., Sankararaman, S., Nielsen, R. & Huerta-Sánchez, E. Evidence for archaic adaptive introgression in humans. Nat. Rev. Genet. 16, 359–371 (2015).

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Juric, I., Aeschbacher, S. & Coop, G. The strength of selection against Neanderthal introgression. PLoS Genet. 12, e1006340 (2016).

    PubMed  PubMed Central  Google Scholar 

  18. Harris, K. & Nielsen, R. The genetic cost of Neanderthal introgression. Genetics 3, 881–891 (2016).

    Google Scholar 

  19. Gittelman, R. M. et al. Archaic hominin admixture facilitated adaptation to out-of-Africa environments. Curr. Biol. 26, 3375–3382 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Petr, M., Pääbo, S., Kelso, J. & Vernot, B. Limits of long-term selection against Neandertal introgression. Proc. Natl Acad. Sci. USA 116, 1639–1644 (2019).

  21. Dannemann, M. & Kelso, J. The contribution of Neanderthals to phenotypic variation in modern humans. Am. J. Hum. Genet. 101, 578–589 (2017).

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Simonti, C. N. et al. The phenotypic legacy of admixture between modern humans and Neandertals. Science 351, 737–741 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Nédélec, Y. et al. Genetic ancestry and natural selection drive population differences in immune responses to pathogens. Cell 167, 657–669 (2016).

  24. Quach, H. et al. Genetic adaptation and Neandertal admixture shaped the immune system of human populations. Cell 167, 643–656 (2016).

  25. Sams, A. J. et al. Adaptively introgressed Neandertal haplotype at the OAS locus functionally impacts innate immune responses in humans. Genome Biol. 17, 246 (2016).

    PubMed  PubMed Central  Google Scholar 

  26. Hu, Y., Ding, Q., He, Y., Xu, S. & Jin, L. Reintroduction of a homocysteine level-associated allele into East Asians by Neanderthal introgression. Mol. Biol. Evol. 32, msv176 (2015).

    Google Scholar 

  27. Vernot, B. et al. Excavating Neandertal and Denisovan DNA from the genomes of melanesian individuals. Science 352, 235–239 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Kim, B. Y. & Lohmueller, K. E. Selection and reduced population size cannot explain higher amounts of Neandertal ancestry in East Asian than in European human populations. Am. J. Hum. Genet. 96, 454–461 (2015).

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Villanea, F. A. & Schraiber, J. G. Multiple episodes of interbreeding between Neanderthal and modern humans. Nat. Ecol. Evol. 3, 39–44 (2019).

    PubMed  Google Scholar 

  30. Wall, J. D. et al. Higher levels of Neanderthal ancestry in East Asians than in Europeans. Genetics 194, 199–209 (2013).

    PubMed  PubMed Central  Google Scholar 

  31. MacArthur, J. et al. The new NHGRI-EBI Catalog of published genome-wide association studies (GWAS Catalog). Nucleic Acids Res. 45, D896–D901 (2017).

    CAS  Google Scholar 

  32. Franke, A. et al. Genome-wide meta-analysis increases to 71 the number of confirmed Crohn’s disease susceptibility loci. Nat. Genet. 42, 1118–1125 (2010).

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Jostins, L. et al. Host–microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature 491, 119–124 (2012).

    CAS  PubMed  PubMed Central  Google Scholar 

  34. Lee, M. K. et al. Genome-wide association study of facial morphology reveals novel associations with FREM1 and PARK2. PLoS ONE 12, e0176566 (2017).

    PubMed  PubMed Central  Google Scholar 

  35. Park, S. L. et al. Mercapturic acids derived from the toxicants acrolein and crotonaldehyde in the urine of cigarette smokers from five ethnic groups with differing risks for lung cancer. PLoS ONE 10, e0124841 (2015).

    PubMed  PubMed Central  Google Scholar 

  36. Spada, J. et al. Genome-wide association analysis of actigraphic sleep phenotypes in the LIFE Adult Study. J. Sleep Res. 25, 690–701 (2016).

    PubMed  Google Scholar 

  37. Kulminski, A. M. et al. Strong impact of natural-selection-free heterogeneity in genetics of age-related phenotypes. Aging 10, 492–514 (2018).

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Lutz, S. M. et al. A genome-wide association study identifies risk loci for spirometric measures among smokers of European and African ancestry. BMC Genet. 16, 138 (2015).

    PubMed  PubMed Central  Google Scholar 

  39. McCoy, R. C., Wakefield, J. & Akey, J. M. Impacts of Neanderthal-introgressed sequences on the landscape of human gene expression. Cell 168, 916–927 (2017).

    CAS  PubMed  PubMed Central  Google Scholar 

  40. GTEx Consortium. Genetic effects on gene expression across human tissues. Nature 550, 204–213 (2017).

    PubMed Central  Google Scholar 

  41. Lappalainen, T. et al. Transcriptome and genome sequencing uncovers functional variation in humans. Nature 501, 506–511 (2013).

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Cat Eye Syndrome (CES) (OMIM, accessed 1 November 2018); https://www.omim.org/entry/115470

  43. Tewhey, R. et al. Direct identification of hundreds of expression-modulating variants using a multiplexed reporter assay. Cell 165, 1519–1529 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  44. van Arensbergen, J. et al. High-throughput identification of human SNPs affecting regulatory element activity. Nat. Genet. 51, 1160–1169 (2019).

    PubMed  PubMed Central  Google Scholar 

  45. Brawand, D. et al. The evolution of gene expression levels in mammalian organs. Nature 478, 343–348 (2011).

    CAS  PubMed  Google Scholar 

  46. Colbran, L. L. et al. Inferred divergent gene regulation in archaic hominins reveals potential phenotypic differences. Nat. Ecol. Evol. 3, 1598–1606 (2019).

    PubMed  PubMed Central  Google Scholar 

  47. Telis, N., Aguilar, R. & Harris, K. Selection against archaic DNA in human regulatory regions. Preprint at bioRxiv https://doi.org/10.1101/708230 (2019).

  48. Peyregne, S., Boyle, M. J., Dannemann, M. & Prufer, K. Detecting ancient positive selection in humans using extended lineage sorting. Genome Res. 27, 1563–1572 (2017).

    CAS  PubMed  PubMed Central  Google Scholar 

  49. Castellano, S. et al. Patterns of coding variation in the complete exomes of three Neandertals. Proc. Natl Acad. Sci. USA 111, 6666–6671 (2014).

    CAS  PubMed  Google Scholar 

  50. Henn, B. M. et al. Distance from sub-Saharan Africa predicts mutational load in diverse human genomes. Proc. Natl Acad. Sci. USA 113, E440–E449 (2016).

    CAS  PubMed  Google Scholar 

  51. Lohmueller, K. E. The distribution of deleterious genetic variation in human populations. Curr. Opin. Genet. Dev. 29, 139–146 (2014).

    CAS  PubMed  Google Scholar 

  52. Lohmueller, K. E. et al. Proportionally more deleterious genetic variation in European than in African populations. Nature 451, 994–997 (2008).

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Simons, Y. B., Turchin, M. C., Pritchard, J. K. & Sella, G. The deleterious mutation load is insensitive to recent population history. Nat. Genet. 46, 220–224 (2014).

    CAS  PubMed  PubMed Central  Google Scholar 

  54. Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156–2158 (2011).

    CAS  PubMed  PubMed Central  Google Scholar 

  55. Chen, L., Wolf, A. B., Fu, W., Li, L. & Akey, J. M. Identifying and interpreting apparent Neanderthal ancestry in African individuals. Cell 180, 677–687 (2020).

    CAS  PubMed  Google Scholar 

  56. Bergström, A. et al. Insights into human genetic variation and population history from 929 diverse genomes. Science 367, eaay5012 (2020).

    PubMed  PubMed Central  Google Scholar 

  57. Hodgkinson, A. & Eyre-Walker, A. Variation in the mutation rate across mammalian genomes. Nat. Rev. Genet. 12, 756–766 (2011).

    CAS  PubMed  Google Scholar 

  58. Kircher, M. et al. A general framework for estimating the relative pathogenicity of human genetic variants. Nat. Genet. 46, 310–315 (2014).

    CAS  PubMed  PubMed Central  Google Scholar 

  59. Pers, T. H., Timshel, P. & Hirschhorn, J. N. SNPsnap: a web-based tool for identification and annotation of matched SNPs. Bioinformatics 31, 418–420 (2015).

    CAS  PubMed  Google Scholar 

  60. Haller, B. C. & Messer, P. W. SLiM 2: Flexible, interactive forward genetic simulations. Mol. Biol. Evol. 34, 230–240 (2017).

  61. Eyre-Walker, A., Woolfit, M. & Phelps, T. The distribution of fitness effects of new deleterious amino acid mutations in humans. Genetics 173, 891–900 (2006).

    CAS  PubMed  PubMed Central  Google Scholar 

  62. Gravel, S. et al. Demographic history and rare allele sharing among human populations. Proc. Natl Acad. Sci. USA 108, 11983–11988 (2011).

    CAS  Google Scholar 

Download references

Acknowledgements

We thank B. Haller, P. Messer and K. Harris for advice on evolutionary simulations. We thank R. Tewhey for discussions of MPRA results. We thank L. Colbran and other members of the Capra Lab for helpful comments on the figures and manuscript. This work was supported by the National Institutes of Health grant nos. T32EY021453 to C.N.S., T32GM080178 to D.S., K22CA184308 to E.H. and R01GM115836 and R35GM127087 to J.A.C. This work was conducted in part using the resources of the Advanced Computing Center for Research and Education at Vanderbilt University, Nashville, TN, USA.

Author information

Authors and Affiliations

Authors

Contributions

D.C.R., C.N.S., E.M. and J.A.C. conceived and conducted the computational analyses. D.S. and E.H. performed the luciferase assays. D.C.R. and J.A.C. wrote the manuscript with input from all authors.

Corresponding author

Correspondence to John A. Capra.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 Introgressed allele class assignment decision tree and allele count summary.

Decision tree by which 1000 Genomes variants in perfect LD with Neanderthal tag SNPs were classified as RAs and NDAs. The counts of variants making it to each of the numbered steps (1–5) is summarized in the lower table.

Extended Data Fig. 2 Evolutionary simulations suggest that RAs are common and more tolerated than NDAs.

a, The demographic model used to simulate human–Neanderthal admixture and quantify the reintroduction of lost alleles. The model and effective population sizes (Ne) were based on previous simulations of Neanderthal admixture. We considered models in which mutations incurred a fitness cost (mildly purifying selection) or no fitness cost (strict neutrality). Two different admixture fractions (f = 0.02 and f = 0.04) and three mutation rates were used in the simulations (Methods). b, The ratios of RAs to NDAs over 100 simulated Eurasian populations. The simulations predict approximately one RA for every two NDAs, and these estimates are robust to changes in the simulated Neanderthal admixture fraction. Misclassification of non-RAs as RAs due to independent, convergent mutations is extremely rare and the overall false discovery rate for LD-based RA identification is below ~1% (Supplementary Tables 3). While these forward-time simulations only approximate the demographic histories of these populations, the observed RA/NDAs ratio are qualitatively consistent with the simulations (Fig. 2). c, Boxplots summarize the frequencies of these potentially confounding NDAs among all sites that would be called as RAs at the time of admixture (c.f. Figure 1). The incidence of these confounding mutations is slightly higher under a purely neutral model (left) than under a model in which new mutations can be deleterious (right). d, Comparison of the effect of elevated mutation rates on the incidence of potentially confounding variants. Under a neutral model, the false positive rate scales with the mutation rate. The highest rate (μ= 7e-7) provides an estimate for CpG sites and results in a 3% false positive rate. Each boxplot represents 100 simulated populations. e, Selection coefficients in Eurasians from SLiM simulations with high (0.04) and low (0.02) admixture fractions. Each boxplot summarizes the average selection coefficient of all alleles in each introgressed class in each of 100 simulated modern Eurasian populations. The differences between the selection coefficients between RAs and NDAs is large and not dependent upon admixture fraction (~2.5×, P ≈ 0, Mann–Whitney U-test test).

Extended Data Fig. 3 Introgressed allele sharing across three Eurasian populations.

Venn diagram showing the fractions of each introgressed variant class that are shared between populations.

Extended Data Fig. 4 Reintroduced alleles cluster within introgressed Neanderthal haplotypes.

a, Scatter plot of the numbers of RAs and NDAs contained on all introgressed haplotypes in EUR. The correlation between the NDA and RA content is moderate (Pearson’s r2 = 0.46), with 18% of the haplotypes containing no RAs and 10% having more RAs than NDAs. b, Scatter plot of the number of introgressed variants on each haplotype versus haplotype length. The NDA content of a haplotype is proportional to its length (r2 = 0.85), but the number of RAs in each haplotype is less strongly correlated with length (r2 = 0.56). c, Heatmap of the fraction of NDAs and RAs in density percentiles (high to low, left to right) averaged over all introgressed Eurasian haplotypes. This information is summarized in a cumulative density function (CDF) above the heatmaps. A higher fraction of all RAs are found in the most dense percentiles; this reflects the fact that RAs are often present in more dense clusters than are NDAs.

Extended Data Fig. 5 Reintroduced alleles have different predicted fitness effects than Neanderthal-derived alleles.

a, In modern European (EUR) populations, RAs are predicted to be significantly less deleterious than NDAs by CADD (median scaled CADD: NDA=2.67; RA=2.1; P≈0). The upper tail of highly deleterious mutations (CADD score >10) is highlighted in the inset. Results are similar for unscaled scores. b, At the haplotype level, the maximum RA CADD score per introgressed haplotype is significantly lower than for NDAs (median scaled max CADD: NDA=13.3; RA=5.8; P≈0). This is in part due to the overall difference demonstrated in (a) and to the greater number of NDAs per haplotype. RAs are rarely the most deleterious variant per haplotype. Results in East Asian and South Asian populations are similar (data not shown).

Extended Data Fig. 6 PolyPhen2 predicts RAs to be less damaging than NDAs.

PolyPhen2 is more likely to classify RAs as ‘benign’ in all three Eurasian populations. Conversely, NDAs are more likely to be classified as ‘damaging’ in both EUR and SAS populations. The y-axis reports the difference in PolyPhen category membership for each population (that is, the fraction all RAs in the population in the PolyPhen category minus the fraction of all NDAs in population in the PolyPhen category). Per population hypergeometric test is calculated on the enrichment (positive delta) or depletion (negative delta) for RA content within each PolyPhen category.

Extended Data Fig. 7 Distribution of RegulomeDB scores among all introgressed Eurasian SNPs.

Comparison of the composition of NDAs and RAs within in each functional category of RegulomeDB. Boxplots refer to 1000 independently samples sets of background variiants matched on the allele frequency and local LD structure of the highest frequency Neanderthal tag SNP per introgressed EUR haplotype (Vernot 2016).

Extended Data Fig. 8 Comparison of allele frequencies across Eurasian populations.

Stratified by presence/absence of allele in modern sub-Saharan African populations. Reintroduced Ancestral Alleles (RAAs) that are also present in modern African (AFR) populations segregate at higher allele frequencies (AF) in all Eurasian populations than RAAs for which the allele is absent in AFR. Intra-population median differences in AF are displayed along with P values (Mann–Whitney U-test). Outliers are not shown. Circles indicate mean AF.

Extended Data Fig. 9 RA sublass composition of introgressed haplotypes containing GWAS hits.

The fraction of RAs among introgressed alleles on introgressed haplotypes in EUR that contain GWAS Catalog associations in Europeans versus introgressed haplotypes having no reported GWAS associations.

Extended Data Fig. 10 RA fraction in introgressed haplotypes containing eQTL in GTEx tissues.

Summary of the RA fraction among introgressed variants in Neanderthal haplotypes in Europeans (EUR). Boxplots show the distributions RA fractions of all haplotypes containing at least one introgressed eQTL (RA or NDA) in the given GTEx tissue (grey box plots). These distributions are then compared pairwise with distribution for introgressed haplotypes that contain no introgressed GTEx eQTL (top, blue boxplot; n = 4237). Haplotypes containing GTEx eQTL have RA contents higher than non-eQTL containing haplotypes in 46 tissues, with 34 of the tissues (*) having a significantly higher the RA fraction (P < 0.05, Mann–Whitney U-Test).

Supplementary information

Supplementary Information

Supplementary Tables 1–11.

Reporting Summary

Supplementary Data 1

List of NDAs and RAs for each Eurasian population.

Supplementary Data 2

Table of GWAS hits for NDAs and RAs.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rinker, D.C., Simonti, C.N., McArthur, E. et al. Neanderthal introgression reintroduced functional ancestral alleles lost in Eurasian populations. Nat Ecol Evol 4, 1332–1341 (2020). https://doi.org/10.1038/s41559-020-1261-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41559-020-1261-z

This article is cited by

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing