Transcription-factor-dependent enhancer transcription defines a gene regulatory network for cardiac rhythm

  1. Xinan H Yang
  2. Rangarajan D Nadadur  Is a corresponding author
  3. Catharina RE Hilvering  Is a corresponding author
  4. Valerio Bianchi  Is a corresponding author
  5. Michael Werner
  6. Stefan R Mazurek
  7. Margaret Gadek  Is a corresponding author
  8. Kaitlyn M Shen  Is a corresponding author
  9. Joseph Aaron Goldman  Is a corresponding author
  10. Leonid Tyan  Is a corresponding author
  11. Jenna Bekeny  Is a corresponding author
  12. Johnathan M Hall
  13. Nutishia Lee  Is a corresponding author
  14. Carlos Perez-Cervantes  Is a corresponding author
  15. Ozanna Burnicka-Turek  Is a corresponding author
  16. Kenneth D Poss
  17. Christopher R Weber  Is a corresponding author
  18. Wouter de Laat  Is a corresponding author
  19. Alexander J Ruthenburg  Is a corresponding author
  20. Ivan P Moskowitz  Is a corresponding author
  1. The University of Chicago, United States
  2. University Medical Center Utrecht, Netherlands
  3. Duke University School of Medicine, United States
  4. Duke University, United States
5 figures and 1 additional file

Figures

Figure 1 with 4 supplements
TF-dependent noncoding transcription defines regulatory elements.

(A) Venn-diagram of peak call overlaps for HL-1 ATAC-sequencing (ATAC.HL1), histone 3 lysine 27 acetylation in mouse heart (H3K27ac.Heart, GSE52123), DNase hypersensitivity in mouse heart …

https://doi.org/10.7554/eLife.31683.002
Figure 1—figure supplement 1
Genomic features of the identified TF-dependent ncRNAs recapitulate known features of annotated lincRNAs.

(A) Pie chart showing the composition for the identified TBX5-dependent ncRNAs using de novo assembly (Cufflinks). The 3067 ncRNAs were categorized into six groups. GENCODE (mm10)-annotated ncRNAs …

https://doi.org/10.7554/eLife.31683.003
Figure 1—figure supplement 2
Genomic view of nine TF-dependent ncRNAs (mm9).

Long noncoding RNAs (lncRNAs) are labeled by nearest TBX5-dependent coding gene. Top track is chromosomal location, followed by the ncRNA read density from R26CreErt2 control and Tbx5fl/fl;R26CreERt2

https://doi.org/10.7554/eLife.31683.004
Figure 1—figure supplement 3
Identifying TF-dependent ncRNA targets from TF-dependent expressed genes.

(A) TF-dependent genes are coherent with TF-dependent ncRNAs. Volcano-plot of TF-dependent genes (FDR < 0.05, |FC| > 1.5). Genes residing within 2 Mb of coherently TBX5-dependent ncRNAs are colored …

https://doi.org/10.7554/eLife.31683.005
Figure 1—figure supplement 4
Identifying TF-dependent ncRNA targets from open chromatin regions.

161 Tbx5-dependent ncRNAs at open chromatin regions were determined by the overlap between ncRNAs and ATAC peaks (allowing up to 500-bp edge-to-edge distance). The empirical significance of …

https://doi.org/10.7554/eLife.31683.006
TF-dependent ncRNAs predict functional enhancers.

(A) Meta-gene plot showing reads per million (RPM) of Tbx5-dependent ncRNAs, being divided by enhancer makers DNase Hypersensitivity (DHS), ATAC Sequencing in HL-1 cells, TBX5 ChIP-Seq, and H3K27Ac, …

https://doi.org/10.7554/eLife.31683.007
Figure 3 with 1 supplement
Enhancer transcription mediates target gene expression.

(A) Genomic loci showing circularized chromatin conformation capture (4C) from the viewpoint of four identified regulatory elements (near Atp2a2, Wrn, Ryr2, and Ubc). Total reads (top) and …

https://doi.org/10.7554/eLife.31683.008
Figure 3—figure supplement 1
A TBX5-dependent enhancer at Sln that was not identified by TF ChIP-seq.

(A) Genomic view of a Tbx5-activated ncRNA region adjacent to gene Sln. This locus shows significant Tbx5-dependence as determined by transcription expression in mouse atrium, open chromatin …

https://doi.org/10.7554/eLife.31683.009
Figure 4 with 2 supplements
Enhancer-associated ncRNAs at Ryr2 and Atp2a2 are necessary for calcium-handling gene expression and cellular phenotype.

(A) Gene Ontology functional enrichment for gene targets within a 2-Mb window of the TBX5-activated enhancers. Odds Ratio (OR) and P-value of overlap with GO terms ‘Calcium Ion Transport Genes’ and …

https://doi.org/10.7554/eLife.31683.010
Figure 4—figure supplement 1
Knockdown of control Hprt locus does not change Ryr2 or RACER expression.

Gene expression of Ryr2 mRNA and RACER transcript after knockdown of Hprt mRNA. Relative transcript expression (RTE) after mRNA knockdown for RACER 1.23 ± 0.25, p = 0.77; for Ryr2 RTE 0.84 ± 0.09, p …

https://doi.org/10.7554/eLife.31683.011
Figure 4—figure supplement 2
Ryr2 enhancer shows cardiac expression in 3 days post fertilization (dpf) zebrafish.

(A) Representative F0 (left) and F1 (right) 3 dpf zebrafish embryos surviving injection with a reporter construct containing the wildtype Ryr2 enhancer reporter, showing cardiac EGFP expression. …

https://doi.org/10.7554/eLife.31683.012
TF-dependent transcription of regulatory elements.

mRNA and enhancer transcription occur in the presence of transcription factor (TBX5) through the Transcription Preinitiation Complex (PIC) (left). Loss of the TBX5 results in decreased PIC complex, …

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

Additional files

Download links