On the limits of fitting complex models of population history to f-statistics

  1. Robert Maier  Is a corresponding author
  2. Pavel Flegontov  Is a corresponding author
  3. Olga Flegontova
  4. Ulas Isildak
  5. Piya Changmai
  6. David Reich  Is a corresponding author
  1. Harvard University, United States
  2. University of Ostrava, Czech Republic
  3. Broad Institute, United States

Abstract

Our understanding of population history in deep time has been assisted by fitting admixture graphs ('AGs') to data: models that specify the ordering of population splits and mixtures, which along with the amount of genetic drift on each lineage and the proportions of mixture, is the only information needed to predict the patterns of allele frequency correlation among populations. Not needing to specify population size changes, split times, or whether admixture events were sudden or drawn out simplifies the space of models that need to be searched. However, the space of possible AGs relating populations is vast and cannot be sampled fully, and thus most published studies have identified fitting AGs through a manual process driven by prior hypotheses, leaving the vast majority of alternative models unexplored. Here, we develop a method for systematically searching the space of all AGs that can incorporate non-genetic information in the form of topology constraints. We implement this findGraphs tool within a software package, ADMIXTOOLS 2, which is a reimplementation of the ADMIXTOOLS software with new features and large performance gains. We apply this methodology to identify alternative models to AGs that played key roles in eight published studies and find that graphs modeling more than six populations and two or three admixture events are often not unique, with many alternative models fitting nominally or significantly better than the published one. Our results suggest that strong claims about population history from AGs should only be made when all well-fitting and temporally plausible models share common topological features. Our re-evaluation of published data also provides insight into the population histories of humans, dogs, and horses, identifying features that are stable across the models we explored, as well as scenarios of populations relationships that differ in important ways from models that have been highlighted in the literature, that fit the allele frequency correlation data, and that are not obviously wrong.

Data availability

As indicated in the 'Materials availability statement', the ancient human genome newly reported in this manuscript (Table S2) is freely available at the European Nucleotide Archive in the form of an alignment of reads to the hg19 human reference genome (project accession number PRJEB58199. All the other data we analyze are previously reported. As we state in the 'Materials availability statement', the exact versions of the published archaeogenetic datasets re-analyzed in this manuscript were kindly shared by the corresponding authors of the following publications upon our requests:1.Bergström A, Frantz L, Schmidt R, et al. Initial Upper Palaeolithic humans in Europe had recent Neanderthal ancestry. Nature. 2021 Apr;592(7853):253-257. doi: 10.1038/s41586-021-03335-3.2.Lazaridis I, Patterson N, Mittnik A, et al. Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature. 2014 Sep 18;513(7518):409-13. doi: 10.1038/nature13673.3.Librado P, Khan N, Fages A, et al. The origins and spread of domestic horses from the Western Eurasian steppes. Nature. 2021 Oct;598(7882):634-640. doi: 10.1038/s41586-021-04018-9.4.Lipson M, Ribot I, Mallick S, et al. Ancient West African foragers in the context of African population history. Nature. 2020 Jan;577(7792):665-670. doi: 10.1038/s41586-020-1929-1.5.Shinde V, Narasimhan VM, Rohland N, et al. An Ancient Harappan Genome Lacks Ancestry from Steppe Pastoralists or Iranian Farmers. Cell. 2019 Oct 17;179(3):729-735.e10. doi: 10.1016/j.cell.2019.08.048.6.Sikora M, Pitulko VV, Sousa VC, et al. The population history of northeastern Siberia since the Pleistocene. Nature. 2019 Jun;570(7760):182-188. Doi: 10.1038/s41586-019-1279-z.7.Wang CC, Yeh HY, Popov AN, et al. Genomic insights into the formation of human populations in East Asia. Nature. 2021 Mar;591(7850):413-419. Doi: 10.1038/s41586-021-03336-2.Various statistics for these re-used datasets are summarized in Table S1.

Article and author information

Author details

  1. Robert Maier

    Harvard University, Cambridge, United States
    For correspondence
    rmaier@broadinstitute.org
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3044-090X
  2. Pavel Flegontov

    Department of Biology and Ecology, University of Ostrava, Ostrava, Czech Republic
    For correspondence
    Pavel_Flegontov@hms.harvard.edu
    Competing interests
    The authors declare that no competing interests exist.
  3. Olga Flegontova

    Department of Biology and Ecology, University of Ostrava, Ostrava, Czech Republic
    Competing interests
    The authors declare that no competing interests exist.
  4. Ulas Isildak

    Department of Biology and Ecology, University of Ostrava, Ostrava, Czech Republic
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6497-6254
  5. Piya Changmai

    Department of Biology and Ecology, University of Ostrava, University of Ostrava, Czech Republic
    Competing interests
    The authors declare that no competing interests exist.
  6. David Reich

    Program in Medical and Population Genetics, Broad Institute, Cambridge, United States
    For correspondence
    reich@genetics.med.harvard.edu
    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

Funding

Czech Ministry of Education, Youth and Sports (project no. LL2103)

  • Pavel Flegontov
  • Olga Flegontova
  • Piya Changmai

Czech Ministry of Education, Youth and Sports (LM2015070)

  • Pavel Flegontov
  • Piya Changmai

Czech Ministry of Education, Youth and Sports (project no. LTAUSA18153)

  • Pavel Flegontov
  • Piya Changmai

National Institutes of Health (GM100233)

  • Robert Maier
  • David Reich

National Institutes of Health (HG012287)

  • Robert Maier
  • David Reich

John Templeton Foundation (grant 61220)

  • Robert Maier
  • David Reich

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

Copyright

© 2023, Maier 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

  • 5,756
    views
  • 931
    downloads
  • 103
    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. Robert Maier
  2. Pavel Flegontov
  3. Olga Flegontova
  4. Ulas Isildak
  5. Piya Changmai
  6. David Reich
(2023)
On the limits of fitting complex models of population history to f-statistics
eLife 12:e85492.
https://doi.org/10.7554/eLife.85492

Share this article

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

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.