Dynamic effects of genetic variation on gene expression revealed following hypoxic stress in cardiomyocytes

Abstract

One life-threatening outcome of cardiovascular disease is myocardial infarction, where cardiomyocytes are deprived of oxygen. To study inter-individual differences in response to hypoxia, we established an in vitro model of induced pluripotent stem cell-derived cardiomyocytes from 15 individuals. We measured gene expression levels, chromatin accessibility, and methylation levels in four culturing conditions that correspond to normoxia, hypoxia and short or long-term re-oxygenation. We characterized thousands of gene regulatory changes as the cells transition between conditions. Using available genotypes, we identified 1,573 genes with a cis expression quantitative locus (eQTL) in at least one condition, as well as 367 dynamic eQTLs, which are classified as eQTLs in at least one, but not in all conditions. A subset of genes with dynamic eQTLs is associated with complex traits and disease. Our data demonstrate how dynamic genetic effects on gene expression, which are likely relevant for disease, can be uncovered under stress.

Data availability

Sequencing data have been deposited in GEO under accession codesGSE144426

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

Article and author information

Author details

  1. Michelle C Ward

    Medicine, University of Chicago, Chicago, United States
    For correspondence
    miward@utmb.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1485-320X
  2. Nicholas E Banovich

    Human Genetics, University of Chicago, Chicago, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Abhishek Sarkar

    Human Genetics, University of Chicago, Chicago, United States
    Competing interests
    The authors declare that no competing interests exist.
  4. Matthew Stephens

    Department of Human Genetics, University of Chicago, Chicago, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Yoav Gilad

    Department of Medicine, University of Chicago, Chicago, United States
    For correspondence
    gilad@uchicago.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8284-8926

Funding

National Heart, Lung, and Blood Institute (HL092206)

  • Yoav Gilad

EMBO Long-Term Fellowship (ALTF 751-2014)

  • Michelle C Ward

National Institute on Aging (F31 AG044948)

  • Nicholas E Banovich

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

Copyright

© 2021, Ward 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,470
    views
  • 337
    downloads
  • 43
    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. Michelle C Ward
  2. Nicholas E Banovich
  3. Abhishek Sarkar
  4. Matthew Stephens
  5. Yoav Gilad
(2021)
Dynamic effects of genetic variation on gene expression revealed following hypoxic stress in cardiomyocytes
eLife 10:e57345.
https://doi.org/10.7554/eLife.57345

Share this article

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

Further reading

    1. Chromosomes and Gene Expression
    Carmina Lichauco, Eric J Foss ... Antonio Bedalov
    Research Article

    The association between late replication timing and low transcription rates in eukaryotic heterochromatin is well known, yet the specific mechanisms underlying this link remain uncertain. In Saccharomyces cerevisiae, the histone deacetylase Sir2 is required for both transcriptional silencing and late replication at the repetitive ribosomal DNA (rDNA) arrays. We have previously reported that in the absence of SIR2, a de-repressed RNA PolII repositions MCM replicative helicases from their loading site at the ribosomal origin, where they abut well-positioned, high-occupancy nucleosomes, to an adjacent region with lower nucleosome occupancy. By developing a method that can distinguish activation of closely spaced MCM complexes, here we show that the displaced MCMs at rDNA origins have increased firing propensity compared to the nondisplaced MCMs. Furthermore, we found that both activation of the repositioned MCMs and low occupancy of the adjacent nucleosomes critically depend on the chromatin remodeling activity of FUN30. Our study elucidates the mechanism by which Sir2 delays replication timing, and it demonstrates, for the first time, that activation of a specific replication origin in vivo relies on the nucleosome context shaped by a single chromatin remodeler.

    1. Chromosomes and Gene Expression
    2. Structural Biology and Molecular Biophysics
    Liza Dahal, Thomas GW Graham ... Xavier Darzacq
    Research Article

    Type II nuclear receptors (T2NRs) require heterodimerization with a common partner, the retinoid X receptor (RXR), to bind cognate DNA recognition sites in chromatin. Based on previous biochemical and overexpression studies, binding of T2NRs to chromatin is proposed to be regulated by competition for a limiting pool of the core RXR subunit. However, this mechanism has not yet been tested for endogenous proteins in live cells. Using single-molecule tracking (SMT) and proximity-assisted photoactivation (PAPA), we monitored interactions between endogenously tagged RXR and retinoic acid receptor (RAR) in live cells. Unexpectedly, we find that higher expression of RAR, but not RXR, increases heterodimerization and chromatin binding in U2OS cells. This surprising finding indicates the limiting factor is not RXR but likely its cadre of obligate dimer binding partners. SMT and PAPA thus provide a direct way to probe which components are functionally limiting within a complex TF interaction network providing new insights into mechanisms of gene regulation in vivo with implications for drug development targeting nuclear receptors.