A community-maintained standard library of population genetic models

  1. Jeffrey R Adrion
  2. Christopher B Cole
  3. Noah Dukler
  4. Jared G Galloway
  5. Ariella L Gladstein
  6. Graham Gower
  7. Christopher C Kyriazis
  8. Aaron P Ragsdale
  9. Georgia Tsambos
  10. Franz Baumdicker
  11. Jedidiah Carlson
  12. Reed A Cartwright
  13. Arun Durvasula
  14. Ilan Gronau
  15. Bernard Y Kim
  16. Patrick McKenzie
  17. Philipp W Messer
  18. Ekaterina Noskova
  19. Diego Ortega Del Vecchyo
  20. Fernando Racimo
  21. Travis J Struck
  22. Simon Gravel
  23. Ryan N Gutenkunst
  24. Kirk E Lohmueller
  25. Peter L Ralph
  26. Daniel R Schrider
  27. Adam Siepel
  28. Jerome Kelleher  Is a corresponding author
  29. Andrew D Kern  Is a corresponding author
  1. University of Oregon, United States
  2. University of Oxford, United Kingdom
  3. Cold Spring Harbor Laboratory, United States
  4. University of North Carolina at Chapel Hill, United States
  5. The University of Adelaide, Australia
  6. University of California, Los Angeles, United States
  7. McGill University, Canada
  8. University of Melbourne, Australia
  9. University of Freiburg, Germany
  10. University of Washington, United States
  11. Arizona State University, United States
  12. IDC Herzliya, Israel
  13. Stanford University, United States
  14. Columbia University, United States
  15. Cornell University, United States
  16. ITMO University, Russian Federation
  17. National Autonomous University of Mexico, Mexico
  18. University of Copenhagen, Denmark
  19. University of Arizona, United States

Abstract

The explosion in population genomic data demands ever more complex modes of analysis, and increasingly these analyses depend on sophisticated simulations. Re-cent advances in population genetic simulation have made it possible to simulate large and complex models, but specifying such models for a particular simulation engine remains a difficult and error-prone task. Computational genetics researchers currently re-implement simulation models independently, leading to inconsistency and duplication of effort. This situation presents a major barrier to empirical researchers seeking to use simulations for power analyses of upcoming studies or sanity checks on existing genomic data. Population genetics, as a field, also lacks standard benchmarks by which new tools for inference might be measured. Here we describe a new resource, stdpopsim, that attempts to rectify this situation. Stdpopsim is a community-driven open source project, which provides easy access to a growing catalog of published simulation models from a range of organisms and supports multiple simulation engine backends. This resource is available as a well-documented python library with a simple command-line interface. We share some examples demonstrating how stdpopsim can be used to systematically compare demographic inference methods, and we encourage a broader community of developers to contribute to this growing resource.

Data availability

All resources are available from https://github.com/popsim-consortium/stdpopsim

Article and author information

Author details

  1. Jeffrey R Adrion

    Department of Biology, University of Oregon, Eugene, United States
    Competing interests
    No competing interests declared.
  2. Christopher B Cole

    Weatherall Institute of Molecular Medicine, University of Oxford, Oxford, United Kingdom
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6733-633X
  3. Noah Dukler

    Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8739-8052
  4. Jared G Galloway

    Department of Biology and Institute of Ecology and Evolution, University of Oregon, Eugene, United States
    Competing interests
    No competing interests declared.
  5. Ariella L Gladstein

    Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, United States
    Competing interests
    No competing interests declared.
  6. Graham Gower

    Australian Centre for Ancient DNA, School of Biological Sciences, The University of Adelaide, Adelaide, Australia
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6197-3872
  7. Christopher C Kyriazis

    Department of Ecology and Evolutionary Biology, University of California, Los Angeles, Los Angeles, United States
    Competing interests
    No competing interests declared.
  8. Aaron P Ragsdale

    Human Genetics, McGill University, Montreal, Canada
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0715-3432
  9. Georgia Tsambos

    Melbourne Integrative Genomics, School of Mathematics and Statistics, University of Melbourne, Melbourne, Australia
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7001-2275
  10. Franz Baumdicker

    Department of Mathematical Stochastics, University of Freiburg, Freiburg, Germany
    Competing interests
    No competing interests declared.
  11. Jedidiah Carlson

    Department of Genome Sciences, University of Washington, Seattle, United States
    Competing interests
    No competing interests declared.
  12. Reed A Cartwright

    The Biodesign Institute and The School of Life Sciences, Arizona State University, Tempe, United States
    Competing interests
    No competing interests declared.
  13. Arun Durvasula

    Department of Human Genetics, University of California, Los Angeles, Los Angeles, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0631-3238
  14. Ilan Gronau

    IDC Herzliya, Herzliya, Israel
    Competing interests
    No competing interests declared.
  15. Bernard Y Kim

    Department of Biology, Stanford University, Stanford, United States
    Competing interests
    No competing interests declared.
  16. Patrick McKenzie

    Department of Ecology, Evolution, and Environmental Biology, Columbia University, New York, United States
    Competing interests
    No competing interests declared.
  17. Philipp W Messer

    Department of Biological Statistics and Computational Biology, Cornell University, Ithaca, United States
    Competing interests
    Philipp W Messer, Reviewing editor, eLife.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8453-9377
  18. Ekaterina Noskova

    Computer Technologies Laboratory, ITMO University, Saint Petersburg, Russian Federation
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1168-0497
  19. Diego Ortega Del Vecchyo

    International Laboratory for Human Genome Research, National Autonomous University of Mexico, Juriquilla, Mexico
    Competing interests
    No competing interests declared.
  20. Fernando Racimo

    Lundbeck GeoGenetics Centre, Globe Institute, University of Copenhagen, Copenhagen, Denmark
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5025-2607
  21. Travis J Struck

    Molecular and Cellular Biology, University of Arizona, Tucson, United States
    Competing interests
    No competing interests declared.
  22. Simon Gravel

    Human Genetics, McGill University, Montreal, Canada
    Competing interests
    No competing interests declared.
  23. Ryan N Gutenkunst

    Molecular and Cellular Biology, University of Arizona, Tucson, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8659-0579
  24. Kirk E Lohmueller

    Department of Ecology and Evolutionary Biology, University of California, Los Angeles, Los Angeles, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3874-369X
  25. Peter L Ralph

    Institute of Ecology and Evolution, University of Oregon, Eugene, United States
    Competing interests
    No competing interests declared.
  26. Daniel R Schrider

    Department of Genetics, University of North Carolina at Chapel Hill, Chapel Hill, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5249-4151
  27. Adam Siepel

    Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, United States
    Competing interests
    No competing interests declared.
  28. Jerome Kelleher

    Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, University of Oxford, Oxford, United Kingdom
    For correspondence
    jerome.kelleher@bdi.ox.ac.uk
    Competing interests
    No competing interests declared.
  29. Andrew D Kern

    Institute of Ecology and Evolution, University of Oregon, Eugene, United States
    For correspondence
    adkern@uoregon.edu
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4381-4680

Funding

National Institute of General Medical Sciences (R35GM119856)

  • Christopher C Kyriazis
  • Kirk E Lohmueller

National Institute of General Medical Sciences (R01GM117241)

  • Jeffrey R Adrion
  • Andrew D Kern

National Institute of General Medical Sciences (R01GM127348)

  • Travis J Struck
  • Ryan N Gutenkunst

National Institute of General Medical Sciences (R00HG008696)

  • Ariella L Gladstein
  • Daniel R Schrider

National Institute of General Medical Sciences (R35GM127070)

  • Noah Dukler
  • Adam Siepel

National Human Genome Research Institute (R01HG010346)

  • Noah Dukler
  • Adam Siepel

Villum Fonden (00025300)

  • Graham Gower
  • Fernando Racimo

UC MEXUS-CONACYT

  • Diego Ortega Del Vecchyo

Robertson Foundation

  • Jerome Kelleher

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

Copyright

© 2020, Adrion 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

  • 6,544
    views
  • 663
    downloads
  • 149
    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. Jeffrey R Adrion
  2. Christopher B Cole
  3. Noah Dukler
  4. Jared G Galloway
  5. Ariella L Gladstein
  6. Graham Gower
  7. Christopher C Kyriazis
  8. Aaron P Ragsdale
  9. Georgia Tsambos
  10. Franz Baumdicker
  11. Jedidiah Carlson
  12. Reed A Cartwright
  13. Arun Durvasula
  14. Ilan Gronau
  15. Bernard Y Kim
  16. Patrick McKenzie
  17. Philipp W Messer
  18. Ekaterina Noskova
  19. Diego Ortega Del Vecchyo
  20. Fernando Racimo
  21. Travis J Struck
  22. Simon Gravel
  23. Ryan N Gutenkunst
  24. Kirk E Lohmueller
  25. Peter L Ralph
  26. Daniel R Schrider
  27. Adam Siepel
  28. Jerome Kelleher
  29. Andrew D Kern
(2020)
A community-maintained standard library of population genetic models
eLife 9:e54967.
https://doi.org/10.7554/eLife.54967

Share this article

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

Further reading

    1. Computational and Systems Biology
    Yang Tan, Bingxin Zhou ... Liang Hong
    Research Article

    Protein engineering is a pivotal aspect of synthetic biology, involving the modification of amino acids within existing protein sequences to achieve novel or enhanced functionalities and physical properties. Accurate prediction of protein variant effects requires a thorough understanding of protein sequence, structure, and function. Deep learning methods have demonstrated remarkable performance in guiding protein modification for improved functionality. However, existing approaches predominantly rely on protein sequences, which face challenges in efficiently encoding the geometric aspects of amino acids’ local environment and often fall short in capturing crucial details related to protein folding stability, internal molecular interactions, and bio-functions. Furthermore, there lacks a fundamental evaluation for developed methods in predicting protein thermostability, although it is a key physical property that is frequently investigated in practice. To address these challenges, this article introduces a novel pre-training framework that integrates sequential and geometric encoders for protein primary and tertiary structures. This framework guides mutation directions toward desired traits by simulating natural selection on wild-type proteins and evaluates variant effects based on their fitness to perform specific functions. We assess the proposed approach using three benchmarks comprising over 300 deep mutational scanning assays. The prediction results showcase exceptional performance across extensive experiments compared to other zero-shot learning methods, all while maintaining a minimal cost in terms of trainable parameters. This study not only proposes an effective framework for more accurate and comprehensive predictions to facilitate efficient protein engineering, but also enhances the in silico assessment system for future deep learning models to better align with empirical requirements. The PyTorch implementation is available at https://github.com/ai4protein/ProtSSN.

    1. Computational and Systems Biology
    2. Microbiology and Infectious Disease
    Saugat Poudel, Jason Hyun ... Bernhard O Palsson
    Research Article

    The Staphylococcus aureus clonal complex 8 (CC8) is made up of several subtypes with varying levels of clinical burden; from community-associated methicillin-resistant S. aureus USA300 strains to hospital-associated (HA-MRSA) USA500 strains and ancestral methicillin-susceptible (MSSA) strains. This phenotypic distribution within a single clonal complex makes CC8 an ideal clade to study the emergence of mutations important for antibiotic resistance and community spread. Gene-level analysis comparing USA300 against MSSA and HA-MRSA strains have revealed key horizontally acquired genes important for its rapid spread in the community. However, efforts to define the contributions of point mutations and indels have been confounded by strong linkage disequilibrium resulting from clonal propagation. To break down this confounding effect, we combined genetic association testing with a model of the transcriptional regulatory network (TRN) to find candidate mutations that may have led to changes in gene regulation. First, we used a De Bruijn graph genome-wide association study to enrich mutations unique to the USA300 lineages within CC8. Next, we reconstructed the TRN by using independent component analysis on 670 RNA-sequencing samples from USA300 and non-USA300 CC8 strains which predicted several genes with strain-specific altered expression patterns. Examination of the regulatory region of one of the genes enriched by both approaches, isdH, revealed a 38-bp deletion containing a Fur-binding site and a conserved single-nucleotide polymorphism which likely led to the altered expression levels in USA300 strains. Taken together, our results demonstrate the utility of reconstructed TRNs to address the limits of genetic approaches when studying emerging pathogenic strains.