Systematic identification of cis-regulatory variants that cause gene expression differences in a yeast cross

  1. Kaushik Renganaath
  2. Rockie Chong
  3. Laura Day
  4. Sriram Kosuri
  5. Leonid Kruglyak  Is a corresponding author
  6. Frank Wolfgang Albert  Is a corresponding author
  1. University of Minnesota, United States
  2. University of California, Los Angeles, United States

Abstract

Sequence variation in regulatory DNA alters gene expression and shapes genetically complex traits. However, the identification of individual, causal regulatory variants is challenging. Here, we used a massively parallel reporter assay to measure the cis-regulatory consequences of 5,832 natural DNA variants in the promoters of 2,503 genes in the yeast Saccharomyces cerevisiae. We identified 451 causal variants, which underlie genetic loci known to affect gene expression. Several promoters harbored multiple causal variants. In five promoters, pairs of variants showed non-additive, epistatic interactions. Causal variants were enriched at conserved nucleotides, tended to have low derived allele frequency, and were depleted from promoters of essential genes, which is consistent with the action of negative selection. Causal variants were also enriched for alterations in transcription factor binding sites. Models integrating these features provided modest, but statistically significant, ability to predict causal variants. This work revealed a complex molecular basis for cis-acting regulatory variation.

Data availability

Raw data and barcode assignments to oligos are available under GEO accession GSE155944. Source Data is provided for Figures 2, 3, 4, 5, and 6. Additional processed data and the MPRA design are available as Supplementary Files.

The following data sets were generated
The following previously published data sets were used

Article and author information

Author details

  1. Kaushik Renganaath

    Department of Genetics, Cell Biology, & Development, University of Minnesota, Minneapolis, 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-1010-3604
  2. Rockie Chong

    Department of Chemistry & Biochemistry, University of California, Los Angeles, Los Angeles, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Laura Day

    Department of Human Genetics, University of California, Los Angeles, Los Angeles, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Sriram Kosuri

    Chemistry and Biochemistry, 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-4661-0600
  5. Leonid Kruglyak

    Department of Human Genetics, University of California, Los Angeles, Los Angeles, United States
    For correspondence
    LKruglyak@mednet.ucla.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8065-3057
  6. Frank Wolfgang Albert

    Department of Genetics, Cell Biology, and Development, University of Minnesota, Minneapolis, United States
    For correspondence
    falbert@umn.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1380-8063

Funding

National Institutes of Health (R35GM124676)

  • Frank Wolfgang Albert

Howard Hughes Medical Institute

  • Leonid Kruglyak

Pew Charitable Trusts

  • Frank Wolfgang Albert

Alfred P. Sloan Foundation

  • Frank Wolfgang Albert

Kinship Foundation

  • Sriram Kosuri

Department of Energy, Labor and Economic Growth (DE-FC02-02ER63421)

  • Sriram Kosuri

National Institutes of Health (R01GM102308)

  • Leonid Kruglyak

National Institutes of Health (DP2GM114829)

  • Sriram Kosuri

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

Copyright

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

  • 2,990
    views

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

Download links

Share this article

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

Further reading

    1. Genetics and Genomics
    Yi Li, Long Gong ... Shangbang Gao
    Research Article

    Resistance to anthelmintics, particularly the macrocyclic lactone ivermectin (IVM), presents a substantial global challenge for parasite control. We found that the functional loss of an evolutionarily conserved E3 ubiquitin ligase, UBR-1, leads to IVM resistance in Caenorhabditis elegans. Multiple IVM-inhibiting activities, including viability, body size, pharyngeal pumping, and locomotion, were significantly ameliorated in various ubr-1 mutants. Interestingly, exogenous application of glutamate induces IVM resistance in wild-type animals. The sensitivity of all IVM-affected phenotypes of ubr-1 is restored by eliminating proteins associated with glutamate metabolism or signaling: GOT-1, a transaminase that converts aspartate to glutamate, and EAT-4, a vesicular glutamate transporter. We demonstrated that IVM-targeted GluCls (glutamate-gated chloride channels) are downregulated and that the IVM-mediated inhibition of serotonin-activated pharynx Ca2+ activity is diminished in ubr-1. Additionally, enhancing glutamate uptake in ubr-1 mutants through ceftriaxone completely restored their IVM sensitivity. Therefore, UBR-1 deficiency-mediated aberrant glutamate signaling leads to ivermectin resistance in C. elegans.

    1. Genetics and Genomics
    Minsoo Noh, Xiangguo Che ... Sihoon Lee
    Research Article

    Osteoporosis, characterized by reduced bone density and strength, increases fracture risk, pain, and limits mobility. Established therapies of parathyroid hormone (PTH) analogs effectively promote bone formation and reduce fractures in severe osteoporosis, but their use is limited by potential adverse effects. In the pursuit of safer osteoporosis treatments, we investigated R25CPTH, a PTH variant wherein the native arginine at position 25 is substituted by cysteine. These studies were prompted by our finding of high bone mineral density in a hypoparathyroidism patient with the R25C homozygous mutation, and we explored its effects on PTH type-1 receptor (PTH1R) signaling in cells and bone metabolism in mice. Our findings indicate that R25CPTH(1–84) forms dimers both intracellularly and extracellularly, and the synthetic dimeric peptide, R25CPTH(1–34), exhibits altered activity in PTH1R-mediated cyclic AMP (cAMP) response. Upon a single injection in mice, dimeric R25CPTH(1–34) induced acute calcemic and phosphaturic responses comparable to PTH(1–34). Furthermore, repeated daily injections increased calvarial bone thickness in intact mice and improved trabecular and cortical bone parameters in ovariectomized (OVX) mice, akin to PTH(1–34). The overall results reveal a capacity of a dimeric PTH peptide ligand to activate the PTH1R in vitro and in vivo as PTH, suggesting a potential path of therapeutic PTH analog development.