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

  • 4,001
    views
  • 518
    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. Evolutionary Biology
    2. Microbiology and Infectious Disease
    Zach Hensel
    Short Report

    Accurate estimation of the effects of mutations on SARS-CoV-2 viral fitness can inform public-health responses such as vaccine development and predicting the impact of a new variant; it can also illuminate biological mechanisms including those underlying the emergence of variants of concern. Recently, Lan et al. reported a model of SARS-CoV-2 secondary structure and its underlying dimethyl sulfate reactivity data (Lan et al., 2022). I investigated whether base reactivities and secondary structure models derived from them can explain some variability in the frequency of observing different nucleotide substitutions across millions of patient sequences in the SARS-CoV-2 phylogenetic tree. Nucleotide basepairing was compared to the estimated ‘mutational fitness’ of substitutions, a measurement of the difference between a substitution’s observed and expected frequency that is correlated with other estimates of viral fitness (Bloom and Neher, 2023). This comparison revealed that secondary structure is often predictive of substitution frequency, with significant decreases in substitution frequencies at basepaired positions. Focusing on the mutational fitness of C→U, the most common type of substitution, I describe C→U substitutions at basepaired positions that characterize major SARS-CoV-2 variants; such mutations may have a greater impact on fitness than appreciated when considering substitution frequency alone.

    1. Evolutionary Biology
    Yiheng Zhang, Xing Wang ... Xiaoguang Yang
    Research Article

    Although fossil evidence suggests the existence of an early muscular system in the ancient cnidarian jellyfish from the early Cambrian Kuanchuanpu biota (ca. 535 Ma), south China, the mechanisms underlying the feeding and respiration of the early jellyfish are conjectural. Recently, the polyp inside the periderm of olivooids was demonstrated to be a calyx-like structure, most likely bearing short tentacles and bundles of coronal muscles at the edge of the calyx, thus presumably contributing to feeding and respiration. Here, we simulate the contraction and expansion of the microscopic periderm-bearing olivooid Quadrapyrgites via the fluid-structure interaction computational fluid dynamics (CFD) method to investigate their feeding and respiratory activities. The simulations show that the rate of water inhalation by the polyp subumbrella is positively correlated with the rate of contraction and expansion of the coronal muscles, consistent with the previous feeding and respiration hypothesis. The dynamic simulations also show that the frequent inhalation/exhalation of water through the periderm polyp expansion/contraction conducted by the muscular system of Quadrapyrgites most likely represents the ancestral feeding and respiration patterns of Cambrian sedentary medusozoans that predated the rhythmic jet-propelled swimming of the modern jellyfish. Most importantly for these Cambrian microscopic sedentary medusozoans, the increase of body size and stronger capacity of muscle contraction may have been indispensable in the stepwise evolution of active feeding and subsequent swimming in a higher flow (or higher Reynolds number) environment.