Barcoded Bulk QTL mapping reveals highly polygenic and epistatic architecture of complex traits in yeast

  1. Alex N Nguyen Ba
  2. Katherine R Lawrence
  3. Artur Rego-Costa
  4. Shreyas Gopalakrishnan
  5. Daniel Temko
  6. Franziska Michor
  7. Michael M Desai  Is a corresponding author
  1. Harvard University, United States
  2. Massachusetts Institute of Technology, United States
  3. Dana-Farber Cancer Institute, United States

Abstract

Mapping the genetic basis of complex traits is critical to uncovering the biological mechanisms that underlie disease and other phenotypes. Genome-wide association studies (GWAS) in humans and quantitative trait locus (QTL) mapping in model organisms can now explain much of the observed heritability in many traits, allowing us to predict phenotype from genotype. However, constraints on power due to statistical confounders in large GWAS and smaller sample sizes in QTL studies still limit our ability to resolve numerous small-effect variants, map them to causal genes, identify pleiotropic effects across multiple traits, and infer non-additive interactions between loci (epistasis). Here, we introduce barcoded bulk quantitative trait locus (BB-QTL) mapping, which allows us to construct, genotype, and phenotype 100,000 offspring of a budding yeast cross, two orders of magnitude larger than the previous state of the art. We use this panel to map the genetic basis of eighteen complex traits, finding that the genetic architecture of these traits involves hundreds of small-effect loci densely spaced throughout the genome, many with widespread pleiotropic effects across multiple traits. Epistasis plays a central role, with thousands of interactions that provide insight into genetic networks. By dramatically increasing sample size, BB-QTL mapping demonstrates the potential of natural variants in high-powered QTL studies to reveal the highly polygenic, pleiotropic, and epistatic architecture of complex traits.

Data availability

Code used for this study is available at https://github.com/arturrc/bbqtl_inference. FASTQ files from high-throughput sequencing have been deposited in the NCBI BioProject database with accession number PRJNA767876. Inferred genotype and phenotype data is deposited in Dryad (doi: 10.5061/dryad.1rn8pk0vd).

The following data sets were generated

Article and author information

Author details

  1. Alex N Nguyen Ba

    Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, 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-1357-6386
  2. Katherine R Lawrence

    Department of Physics, Massachusetts Institute of Technology, Cambridge, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Artur Rego-Costa

    Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Shreyas Gopalakrishnan

    Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Daniel Temko

    Department of Biostatistics, Harvard University, Cambridge, United States
    Competing interests
    The authors declare that no competing interests exist.
  6. Franziska Michor

    Department of Data Science, Dana-Farber Cancer Institute, Boston, United States
    Competing interests
    The authors declare that no competing interests exist.
  7. Michael M Desai

    Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, United States
    For correspondence
    mdesai@oeb.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-9581-1150

Funding

Natural Sciences and Engineering Research Council of Canada (RGPIN-2021-02716)

  • Alex N Nguyen Ba

National Science Foundation (#1764269)

  • Katherine R Lawrence

National Institutes of Health (U54CA193461)

  • Michael M Desai

National Science Foundation (PHY-1914916)

  • Franziska Michor

National Institutes of Health (GM104239)

  • Michael M Desai

Natural Sciences and Engineering Research Council of Canada (DGECR-2021-00117)

  • Alex N Nguyen Ba

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

Copyright

© 2022, Nguyen Ba 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

  • 3,842
    views
  • 507
    downloads
  • 45
    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. Alex N Nguyen Ba
  2. Katherine R Lawrence
  3. Artur Rego-Costa
  4. Shreyas Gopalakrishnan
  5. Daniel Temko
  6. Franziska Michor
  7. Michael M Desai
(2022)
Barcoded Bulk QTL mapping reveals highly polygenic and epistatic architecture of complex traits in yeast
eLife 11:e73983.
https://doi.org/10.7554/eLife.73983

Share this article

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

Further reading

    1. Computational and Systems Biology
    2. Evolutionary Biology
    Pierre Barrat-Charlaix, Richard A Neher
    Research Article

    As pathogens spread in a population of hosts, immunity is built up, and the pool of susceptible individuals are depleted. This generates selective pressure, to which many human RNA viruses, such as influenza virus or SARS-CoV-2, respond with rapid antigenic evolution and frequent emergence of immune evasive variants. However, the host’s immune systems adapt, and older immune responses wane, such that escape variants only enjoy a growth advantage for a limited time. If variant growth dynamics and reshaping of host-immunity operate on comparable time scales, viral adaptation is determined by eco-evolutionary interactions that are not captured by models of rapid evolution in a fixed environment. Here, we use a Susceptible/Infected model to describe the interaction between an evolving viral population in a dynamic but immunologically diverse host population. We show that depending on strain cross-immunity, heterogeneity of the host population, and durability of immune responses, escape variants initially grow exponentially, but lose their growth advantage before reaching high frequencies. Their subsequent dynamics follows an anomalous random walk determined by future escape variants and results in variant trajectories that are unpredictable. This model can explain the apparent contradiction between the clearly adaptive nature of antigenic evolution and the quasi-neutral dynamics of high-frequency variants observed for influenza viruses.

    1. Ecology
    2. Evolutionary Biology
    Rebecca D Tarvin, Jeffrey L Coleman ... Richard W Fitch
    Research Article

    Understanding the origins of novel, complex phenotypes is a major goal in evolutionary biology. Poison frogs of the family Dendrobatidae have evolved the novel ability to acquire alkaloids from their diet for chemical defense at least three times. However, taxon sampling for alkaloids has been biased towards colorful species, without similar attention paid to inconspicuous ones that are often assumed to be undefended. As a result, our understanding of how chemical defense evolved in this group is incomplete. Here, we provide new data showing that, in contrast to previous studies, species from each undefended poison frog clade have measurable yet low amounts of alkaloids. We confirm that undefended dendrobatids regularly consume mites and ants, which are known sources of alkaloids. Thus, our data suggest that diet is insufficient to explain the defended phenotype. Our data support the existence of a phenotypic intermediate between toxin consumption and sequestration — passive accumulation — that differs from sequestration in that it involves no derived forms of transport and storage mechanisms yet results in low levels of toxin accumulation. We discuss the concept of passive accumulation and its potential role in the origin of chemical defenses in poison frogs and other toxin-sequestering organisms. In light of ideas from pharmacokinetics, we incorporate new and old data from poison frogs into an evolutionary model that could help explain the origins of acquired chemical defenses in animals and provide insight into the molecular processes that govern the fate of ingested toxins.