The lingering effects of Neanderthal introgression on human complex traits

  1. Xinzhu Wei
  2. Christopher R Robles
  3. Ali Pazokitoroudi
  4. Andrea Ganna
  5. Alexander Gusev
  6. Arun Durvasula
  7. Steven Gazal
  8. Po-Ru Loh
  9. David Reich
  10. Sriram Sankararaman  Is a corresponding author
  1. Cornell University, United States
  2. University of California, Los Angeles, United States
  3. Broad Institute, United States
  4. Harvard University, United States
  5. University of Southern California, United States

Abstract

The genetic variants introduced into the ancestors of modern humans from interbreeding with Neanderthals have been suggested to contribute an unexpected extent to complex human traits. However, testing this hypothesis has been challenging due to the idiosyncratic population genetic properties of introgressed variants. We developed rigorous methods to assess the contribution of introgressed Neanderthal variants to heritable trait variation relative to that of modern human variants. We applied these methods to analyze 235,592 introgressed Neanderthal variants and 96 distinct phenotypes measured in about 300,000 unrelated white British individuals in the UK Biobank. Introgressed Neanderthal variants have a significant contribution to trait variation consistent with the polygenic architecture of complex phenotypes (contributing 0.12% of heritable variation averaged across phenotypes). However, the contribution of introgressed variants tends to be significantly depleted relative to modern human variants matched for allele frequency and linkage disequilibrium (about 59% depletion on average), consistent with purifying selection on introgressed variants. Different from previous studies (McArthur 2021), we find no evidence for elevated heritability across the phenotypes examined. We identified 348 independent significant associations of introgressed Neanderthal variants with 64 phenotypes . Previous work (Skov 2020) has suggested that a majority of such associations are likely driven by statistical association with nearby modern human variants that are the true causal variants. We therefore developed a customized statistical fine-mapping methodology for introgressed variants that led us to identify 112 regions (at a false discovery proportion of 16%) across 47 phenotypes containing 4,303 unique genetic variants where introgressed variants are highly likely to have a phenotypic effect. Examination of these variants reveal their substantial impact on genes that are important for the immune system, development, and metabolism. Our results provide the first rigorous basis for understanding how Neanderthal introgression modulates complex trait variation in present-day humans.

Data availability

New software is deposited on GitHub: https://github.com/alipazokit/RHEmc-coeffData for figures (and supplement data) are deposited on GitHub: https://github.com/AprilWei001/NIM

Article and author information

Author details

  1. Xinzhu Wei

    Department of Computational Biology, Cornell University, New York, United States
    Competing interests
    The authors declare that no competing interests exist.
  2. Christopher R Robles

    Department of Human Genetics, University of California, Los Angeles, Los Angeles, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5667-7625
  3. Ali Pazokitoroudi

    Department of Computer Science, University of California, Los Angeles, Los Angeles, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2839-2291
  4. Andrea Ganna

    Program in Medical and Population Genetics, Broad Institute, Cambridge, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Alexander Gusev

    Harvard University, Boston, United States
    Competing interests
    The authors declare that no competing interests exist.
  6. Arun Durvasula

    Department of Genetics, University of California, Los Angeles, Los Angeles, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0631-3238
  7. Steven Gazal

    Department of Public and Population Health Sciences, University of Southern California, Los Angeles, United States
    Competing interests
    The authors declare that no competing interests exist.
  8. Po-Ru Loh

    Program in Medical and Population Genetics, Broad Institute, Boston, United States
    Competing interests
    The authors declare that no competing interests exist.
  9. David Reich

    Program in Medical and Population Genetics, Broad Institute, Boston, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7037-5292
  10. Sriram Sankararaman

    Department of Human Genetics, University of California, Los Angeles, Los Angeles, United States
    For correspondence
    sriram@cs.ucla.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1586-9641

Funding

National Institutes of Health (GM100233)

  • David Reich

Alfred P. Sloan Foundation

  • Sriram Sankararaman

Okawa Foundation

  • Sriram Sankararaman

Burroughs Wellcome Fund

  • Po-Ru Loh

Next Generation Fund at the Broad Institute of MIT and Harvard

  • Po-Ru Loh

National Science Foundation (BCS1032255)

  • David Reich

National Institutes of Health (HG006399)

  • David Reich

Paul G. Allen Frontiers Group (Allen Discovery Center Grant on Brain Evolution)

  • David Reich

John Templeton Foundation (61220)

  • David Reich

Howard Hughes Medical Institute

  • David Reich

National Institutes of Health (R35GM125055)

  • Sriram Sankararaman

National Science Foundation (III1705121)

  • Sriram Sankararaman

National Science Foundation (CAREER1943497)

  • Sriram Sankararaman

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

© 2023, Wei et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 12,288
    views
  • 920
    downloads
  • 24
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Xinzhu Wei
  2. Christopher R Robles
  3. Ali Pazokitoroudi
  4. Andrea Ganna
  5. Alexander Gusev
  6. Arun Durvasula
  7. Steven Gazal
  8. Po-Ru Loh
  9. David Reich
  10. Sriram Sankararaman
(2023)
The lingering effects of Neanderthal introgression on human complex traits
eLife 12:e80757.
https://doi.org/10.7554/eLife.80757

Share this article

https://doi.org/10.7554/eLife.80757

Further reading

    1. Evolutionary Biology
    Ljiljana Mihajlovic, Bharat Ravi Iyengar ... Yolanda Schaerli
    Research Article

    Gene duplication drives evolution by providing raw material for proteins with novel functions. An influential hypothesis by Ohno (1970) posits that gene duplication helps genes tolerate new mutations and thus facilitates the evolution of new phenotypes. Competing hypotheses argue that deleterious mutations will usually inactivate gene duplicates too rapidly for Ohno’s hypothesis to work. We experimentally tested Ohno’s hypothesis by evolving one or exactly two copies of a gene encoding a fluorescent protein in Escherichia coli through several rounds of mutation and selection. We analyzed the genotypic and phenotypic evolutionary dynamics of the evolving populations through high-throughput DNA sequencing, biochemical assays, and engineering of selected variants. In support of Ohno’s hypothesis, populations carrying two gene copies displayed higher mutational robustness than those carrying a single gene copy. Consequently, the double-copy populations experienced relaxed purifying selection, evolved higher phenotypic and genetic diversity, carried more mutations and accumulated combinations of key beneficial mutations earlier. However, their phenotypic evolution was not accelerated, possibly because one gene copy rapidly became inactivated by deleterious mutations. Our work provides an experimental platform to test models of evolution by gene duplication, and it supports alternatives to Ohno’s hypothesis that point to the importance of gene dosage.

    1. Cell Biology
    2. Evolutionary Biology
    Paul Richard J Yulo, Nicolas Desprat ... Heather L Hendrickson
    Research Article

    Maintenance of rod-shape in bacterial cells depends on the actin-like protein MreB. Deletion of mreB from Pseudomonas fluorescens SBW25 results in viable spherical cells of variable volume and reduced fitness. Using a combination of time-resolved microscopy and biochemical assay of peptidoglycan synthesis, we show that reduced fitness is a consequence of perturbed cell size homeostasis that arises primarily from differential growth of daughter cells. A 1000-generation selection experiment resulted in rapid restoration of fitness with derived cells retaining spherical shape. Mutations in the peptidoglycan synthesis protein Pbp1A were identified as the main route for evolutionary rescue with genetic reconstructions demonstrating causality. Compensatory pbp1A mutations that targeted transpeptidase activity enhanced homogeneity of cell wall synthesis on lateral surfaces and restored cell size homeostasis. Mechanistic explanations require enhanced understanding of why deletion of mreB causes heterogeneity in cell wall synthesis. We conclude by presenting two testable hypotheses, one of which posits that heterogeneity stems from non-functional cell wall synthesis machinery, while the second posits that the machinery is functional, albeit stalled. Overall, our data provide support for the second hypothesis and draw attention to the importance of balance between transpeptidase and glycosyltransferase functions of peptidoglycan building enzymes for cell shape determination.