CoChat
Sign Up
Cited by 165Genomics and Chromatin DynamicsChromatin Remodeling and CancerHeat shock proteins research

The DOI matches its Crossref record; 81 of the 81 references it cites were found in a registry or an index.

Transcriptional response to stress is pre-wired by promoter and enhancer architecture

Open access
Vihervaara, Anniina, Mahat, Dig Bijay, Guertin, Michael J., Chu, Tinyi, Danko, Charles G., Lis, John T., Sistonen, Lea·2017·Nature Communications
DOI 10.1038/s41467-017-00151-0PMC PMC5557961
Open this paper in CoChat — free accountOpen PDF

Abstract

Programs of gene expression are executed by a battery of transcription factors that coordinate divergent transcription from a pair of tightly linked core initiation regions of promoters and enhancers. Here, to investigate how divergent transcription is reprogrammed upon stress, we measured nascent RNA synthesis at nucleotide-resolution, and profiled histone H4 acetylation in human cells. Our results globally show that the release of promoter-proximal paused RNA polymerase into elongation functions as a critical switch at which a gene’s response to stress is determined. Highly transcribed and highly inducible genes display strong transcriptional directionality and selective assembly of general transcription factors on the core sense promoter. Heat-induced transcription at enhancers, instead, correlates with prior binding of cell-type, sequence-specific transcription factors. Activated Heat Shock Factor 1 (HSF1) binds to transcription-primed promoters and enhancers, and CTCF-occupied, non-transcribed chromatin. These results reveal chromatin architectural features that orient transcription at divergent regulatory elements and prime transcriptional responses genome-wide.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41467-017-00151-0

Cites

81 of 81 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Rapid transcriptional reprogramming of genes is defined at the step of promoter-proximal pause-release. a Schematic representation (upper panel) of promoter-proximal region (−100 to +400 nt from TSS) from which the Pol II pause site was scored as a 50-nt window, and gene body (+500 nt from TSS to −500 nt from polyA site) from which the average Pol II density across the coding region was measured. Divergent transcription is depicted upstream of the gene, and the transcription initiation sites towards the sense and anti-sense are indicated with arrows. The lower panel illustrates the strand-specific scanning of Pol II density at promoter–proximal (dashed line) and gene body (grey box) regions at BAG3 gene prior to (NHS) and upon (HS) heat shock. b MA-plot (upper panel) showing the heat-induced transcriptional change at the coding regions of individual genes. Genes with significantly upregulated (Up) or downregulated (Down) transcription upon HS are colored red and blue, respectively. The lower panel indicates the number and transcriptional change of genes that were significantly upregulated, downregulated or remained unchanged (UnCh) upon acute heat stress, or that were not transcribed (UnExp) prior to or upon HS in human K562 cells. c Strand-specific average intensity of transcriptionally engaged Pol II at the TSS of upregulated (Up), downregulated (Down), unchanged (UnCh) and unexpressed (UnExp) genes. Coding strand is indicated with solid, divergent strand with dashed line. d Heatmap depicting the change in the Pol II density at the coding strand of significantly changed genes upon acute stress. e The change in the pausing index at individual upregulated, downregulated and unchanged genes. f Comparison of PRO-seq reads prior to (NHS) and upon heat shock (HS) at promoter-proximal and gene body regions of each gene. The density of genes in the scatter plot is indicated with the color scale
Fig. 1 Rapid transcriptional reprogramming of genes is defined at the step of promoter-proximal pause-release. a Schematic representation (upper panel) of promoter-proximal region (−100 to +400 nt from TSS) from which the Pol II pause site was scored as a 50-nt window, and gene body (+500 nt from TSS to −500 nt from polyA site) from which the average Pol II density across the coding region was measured. Divergent transcription is depicted upstream of the gene, and the transcription initiation sites towards the sense and anti-sense are indicated with arrows. The lower panel illustrates the strand-specific scanning of Pol II density at promoter–proximal (dashed line) and gene body (grey box) regions at BAG3 gene prior to (NHS) and upon (HS) heat shock. b MA-plot (upper panel) showing the heat-induced transcriptional change at the coding regions of individual genes. Genes with significantly upregulated (Up) or downregulated (Down) transcription upon HS are colored red and blue, respectively. The lower panel indicates the number and transcriptional change of genes that were significantly upregulated, downregulated or remained unchanged (UnCh) upon acute heat stress, or that were not transcribed (UnExp) prior to or upon HS in human K562 cells. c Strand-specific average intensity of transcriptionally engaged Pol II at the TSS of upregulated (Up), downregulated (Down), unchanged (UnCh) and unexpressed (UnExp) genes. Coding strand is indicated with solid, divergent strand with dashed line. d Heatmap depicting the change in the Pol II density at the coding strand of significantly changed genes upon acute stress. e The change in the pausing index at individual upregulated, downregulated and unchanged genes. f Comparison of PRO-seq reads prior to (NHS) and upon heat shock (HS) at promoter-proximal and gene body regions of each gene. The density of genes in the scatter plot is indicated with the color scaleCC BY · PubMed Central
Fig. 2 Immediate emergence of a stress-specific repertoire of active regulatory elements. a Schematic representation (upper panel) and a browser shot example (lower panel) of a transcribed distal regulatory element (dTRE) showing the characteristic pattern of short divergent transcription. The strand-specific measurement of Pol II density along the length of the dTRE is indicated with the red (plus strand) and blue (minus strand) boxes. The numbers of dTREs in K562 cells that occur prior to (NHS), upon acute heat shock (HS), or in both conditions are shown. b MA-plot showing dTREs with significantly up- or downregulated Pol II density across the length of the regulatory element. c Density of transcriptionally engaged Pol II at individual up- or downregulated dTREs. The dTREs are sorted by the increasing distance between the Pol II pause sites at the sense and the anti-sense strands, and the signal is centered to the middle coordinate between the pause sites
Fig. 2 Immediate emergence of a stress-specific repertoire of active regulatory elements. a Schematic representation (upper panel) and a browser shot example (lower panel) of a transcribed distal regulatory element (dTRE) showing the characteristic pattern of short divergent transcription. The strand-specific measurement of Pol II density along the length of the dTRE is indicated with the red (plus strand) and blue (minus strand) boxes. The numbers of dTREs in K562 cells that occur prior to (NHS), upon acute heat shock (HS), or in both conditions are shown. b MA-plot showing dTREs with significantly up- or downregulated Pol II density across the length of the regulatory element. c Density of transcriptionally engaged Pol II at individual up- or downregulated dTREs. The dTREs are sorted by the increasing distance between the Pol II pause sites at the sense and the anti-sense strands, and the signal is centered to the middle coordinate between the pause sitesCC BY · PubMed Central
Fig. 3 The state of histone H4 acetylation changes in cells exposed to acute stress. a Genome browser images of heat-induced (DNAJB1) and heat-repressed (SNHG8) genes, showing the histone H4 acetylation (black) and transcription from plus (red) and minus (blue) strands prior to (NHS) and upon (HS) heat stress. b Average ChIP-seq intensity of histone H4 acetylation at promoters (upper panels) and gene bodies (lower panels) of genes grouped by their transcriptional response. c The change in histone H4 acetylation upon stress at individual promoters of up- and downregulated genes. The grey dashed line in (b, c) marks the +500 nt position from TSS. d Genome browser images of histone H4 acetylation at dTREs with upregulated (left) and downregulated (right) Pol II density upon stress. e The change in histone H4 acetylation upon stress at individual up- or downregulated dTREs
Fig. 3 The state of histone H4 acetylation changes in cells exposed to acute stress. a Genome browser images of heat-induced (DNAJB1) and heat-repressed (SNHG8) genes, showing the histone H4 acetylation (black) and transcription from plus (red) and minus (blue) strands prior to (NHS) and upon (HS) heat stress. b Average ChIP-seq intensity of histone H4 acetylation at promoters (upper panels) and gene bodies (lower panels) of genes grouped by their transcriptional response. c The change in histone H4 acetylation upon stress at individual promoters of up- and downregulated genes. The grey dashed line in (b, c) marks the +500 nt position from TSS. d Genome browser images of histone H4 acetylation at dTREs with upregulated (left) and downregulated (right) Pol II density upon stress. e The change in histone H4 acetylation upon stress at individual up- or downregulated dTREsCC BY · PubMed Central
Fig. 4 Directionality and rapid transcriptional induction is pre-wired in the promoter architecture. a PRO-seq profile of highly upregulated, highly transcribed, moderately upregulated and all transcribed genes, centered on the Pol II pause site at the coding strand. The intensity of transcriptionally engaged Pol II at the coding strand is indicated above the value zero, respective intensity at the non-coding strand is depicted with negative values. The right panel schematically depicts the Pol II profile at divergent regulatory elements, indicating TSSs to sense and anti-sense directions with arrows, Pol II pause sites at sense and anti-sense strands with dotted lines, and the mid coordinate between the Pol II pause sites with a bracket. b, c Average intensity of (b) TBP ChIP-nexus and (c) NELF-E ChIP-seq at highly upregulated, highly transcribed and all transcribed genes. The corresponding heatmaps are shown in Supplementary Fig. 4c. d Browser shot of HSPA8 gene indicating the positioning of NELF-E and TBP with respect to the GRO-cap mapping of transcription start sites at the sense and the anti-sense strands, and the PRO-seq profile of transcriptionally engaged Pol II prior to (NHS) and upon (HS) heat stress. ChIP-nexus data set for TBP was obtained from He et al.54, ChIP-seq for NELF-E is from the ENCODE55, and GRO-cap from Core et al.19
Fig. 4 Directionality and rapid transcriptional induction is pre-wired in the promoter architecture. a PRO-seq profile of highly upregulated, highly transcribed, moderately upregulated and all transcribed genes, centered on the Pol II pause site at the coding strand. The intensity of transcriptionally engaged Pol II at the coding strand is indicated above the value zero, respective intensity at the non-coding strand is depicted with negative values. The right panel schematically depicts the Pol II profile at divergent regulatory elements, indicating TSSs to sense and anti-sense directions with arrows, Pol II pause sites at sense and anti-sense strands with dotted lines, and the mid coordinate between the Pol II pause sites with a bracket. b, c Average intensity of (b) TBP ChIP-nexus and (c) NELF-E ChIP-seq at highly upregulated, highly transcribed and all transcribed genes. The corresponding heatmaps are shown in Supplementary Fig. 4c. d Browser shot of HSPA8 gene indicating the positioning of NELF-E and TBP with respect to the GRO-cap mapping of transcription start sites at the sense and the anti-sense strands, and the PRO-seq profile of transcriptionally engaged Pol II prior to (NHS) and upon (HS) heat stress. ChIP-nexus data set for TBP was obtained from He et al.54, ChIP-seq for NELF-E is from the ENCODE55, and GRO-cap from Core et al.19CC BY · PubMed Central
Fig. 5 Local chromatin architecture is permissive or restrictive for HSF1-mediated trans-activation. a Transcriptional change of up and downregulated genes, plotted against the transcriptional level upon stress. The HSF1-bound upregulated genes are indicated with red, the HSF1-bound downregulated genes with blue closed circles. b Localization of HSF1 summit point from the Pol II pause site at the coding strand, plotted against the gene’s transcriptional change upon heat shock. c The average occupancy of HSF2 and SP2 at HSF1 target promoters. d Left panel: ENCODE binding score (proportional number between 0–1000; https://genome.ucsc.edu/FAQ/FAQformat.html#format12) of SP2 at the HSF1-bound and HSF1-unbound upregulated (Up) and downregulated (Down) genes. Significant P-values (Mann–Whitney U-test) are shown; the three asterisks indicating values lower than 0.0005. The position of HSF1 peak summit from the peak summit of SP2 is depicted with respect to the directionality of the divergent promoter (right panel). e Scaled ChIP-seq, PRO-seq and GRO-cap intensities at indicated gene groups. The highest average signal intensity for each factor in any bin across the gene groups is used as normalization constant, bringing the maximum signal to value 1. ChIP-seq data sets for NELF-E and SP2 were obtained from the ENCODE55, GRO-cap is from Core et al.19
Fig. 5 Local chromatin architecture is permissive or restrictive for HSF1-mediated trans-activation. a Transcriptional change of up and downregulated genes, plotted against the transcriptional level upon stress. The HSF1-bound upregulated genes are indicated with red, the HSF1-bound downregulated genes with blue closed circles. b Localization of HSF1 summit point from the Pol II pause site at the coding strand, plotted against the gene’s transcriptional change upon heat shock. c The average occupancy of HSF2 and SP2 at HSF1 target promoters. d Left panel: ENCODE binding score (proportional number between 0–1000; https://genome.ucsc.edu/FAQ/FAQformat.html#format12) of SP2 at the HSF1-bound and HSF1-unbound upregulated (Up) and downregulated (Down) genes. Significant P-values (Mann–Whitney U-test) are shown; the three asterisks indicating values lower than 0.0005. The position of HSF1 peak summit from the peak summit of SP2 is depicted with respect to the directionality of the divergent promoter (right panel). e Scaled ChIP-seq, PRO-seq and GRO-cap intensities at indicated gene groups. The highest average signal intensity for each factor in any bin across the gene groups is used as normalization constant, bringing the maximum signal to value 1. ChIP-seq data sets for NELF-E and SP2 were obtained from the ENCODE55, GRO-cap is from Core et al.19CC BY · PubMed Central
Fig. 6 Lineage-specific transcription factors prime enhancers for transcriptional induction. a Transcriptional change of significantly changed dTREs plotted as function of transcription level upon stress. The HSF1-bound upregulated dTREs are indicated with red, the HSF1-bound downregulated dTREs with blue closed circles. The numbers of HSF1-targeted or un-targeted dTREs are indicated. b PRO-seq profile of indicated dTREs in NHS and HS conditions. c HSF1 binding intensity as the function of heat-induced change in the histone H4 acetylation at individual dTREs. The coefficients (rho) and P-values are according to Spearman’s rank correlation. The correlation lines are fitted for linear regression. d Average histone H4 acetylation at indicated dTREs upon HS. e Left panels: Average ChIP-seq intensity of GTF2B, TAF7, TAL1 and GATA2 at up- and downregulated dTREs. Right panels show heatmaps of TAL1 and GATA2 at up- and downregulated dTREs. ChIP-seq data sets for GTF2B, TAF7, GATA2 and TAL1 are from the ENCODE55
Fig. 6 Lineage-specific transcription factors prime enhancers for transcriptional induction. a Transcriptional change of significantly changed dTREs plotted as function of transcription level upon stress. The HSF1-bound upregulated dTREs are indicated with red, the HSF1-bound downregulated dTREs with blue closed circles. The numbers of HSF1-targeted or un-targeted dTREs are indicated. b PRO-seq profile of indicated dTREs in NHS and HS conditions. c HSF1 binding intensity as the function of heat-induced change in the histone H4 acetylation at individual dTREs. The coefficients (rho) and P-values are according to Spearman’s rank correlation. The correlation lines are fitted for linear regression. d Average histone H4 acetylation at indicated dTREs upon HS. e Left panels: Average ChIP-seq intensity of GTF2B, TAF7, TAL1 and GATA2 at up- and downregulated dTREs. Right panels show heatmaps of TAL1 and GATA2 at up- and downregulated dTREs. ChIP-seq data sets for GTF2B, TAF7, GATA2 and TAL1 are from the ENCODE55CC BY · PubMed Central
Fig. 7 HSF1 binds to regions that are devoid of initiating Pol II but gain hyper-acetylation at histone H4 upon stress. a Left panels: Average Pol II density at the HSF1-bound promoters, dTREs and untranscribed genomic regions. N indicates the number, and percentage the share, of HSF1 target sites in respective category. The right panel shows the number and percentage of HSF1 target sites that contain a MEME-ChIP identified HSE (Supplementary Methods). b–e Average intensities of (b) DNaseI sensitivity, (c) CTCF ChIP-seq, (d) HSF1 ChIP-seq and (e) histone H4 acetylation ChIP-seq at HSF1-targeted promoters (purple), dTREs (green) and untranscribed regions (black). f Genome browser examples of a promoter and a dTRE (upper panel), and an untranscribed genomic region (lower panel) that are bound by HSF1 upon heat shock. HSF1 binding (gray), histone H4 acetylation (black) and transcriptionally engaged Pol II (red and blue) are shown in NHS and HS conditions. Data sets of DNaseI sensitivity and CTCF ChIP-seq (Broad institute) were obtained from the ENCODE55
Fig. 7 HSF1 binds to regions that are devoid of initiating Pol II but gain hyper-acetylation at histone H4 upon stress. a Left panels: Average Pol II density at the HSF1-bound promoters, dTREs and untranscribed genomic regions. N indicates the number, and percentage the share, of HSF1 target sites in respective category. The right panel shows the number and percentage of HSF1 target sites that contain a MEME-ChIP identified HSE (Supplementary Methods). b–e Average intensities of (b) DNaseI sensitivity, (c) CTCF ChIP-seq, (d) HSF1 ChIP-seq and (e) histone H4 acetylation ChIP-seq at HSF1-targeted promoters (purple), dTREs (green) and untranscribed regions (black). f Genome browser examples of a promoter and a dTRE (upper panel), and an untranscribed genomic region (lower panel) that are bound by HSF1 upon heat shock. HSF1 binding (gray), histone H4 acetylation (black) and transcriptionally engaged Pol II (red and blue) are shown in NHS and HS conditions. Data sets of DNaseI sensitivity and CTCF ChIP-seq (Broad institute) were obtained from the ENCODE55CC BY · PubMed Central
Fig. 8 Rapid and coordinated reprogramming of genes and distal regulatory elements in stressed human cells. Model showing the rapid heat-induced response of the human genome, including reprogramming of gene transcription and establishment of a stress-specific repertoire of distal regulatory elements. Inhibition of the pause-release of promoter–proximal Pol II clears transcription complexes from the downregulated genes, elevating the concentration of free Pol II in the cell. The heat-induced genes are primed for directionality and rapid activation by the pre-assembled PIC on the core promoter of the coding strand, and by the paused Pol II at the 5′ region of the gene. Upon stress, trans-activators, such as HSF1, launch Pol II from the primed genes into productive elongation. The elevated levels of free Pol II and the highly positioned PIC on the coding strand enable instant loading of Pol II to the freed pause sites, efficiently launching rounds of transcript synthesis from the activated genes. The free Pol II allows also tuning up of the enhancer repertoire, increasing eRNA-production at dTREs that are marked by lineage-specific transcription factors, for example GATA1, GATA2 and TAL1 as shown in this study. In addition to reprogramming TREs, heat stress causes emergence of putative untranscribed regulatory elements, as demonstrated by the acetylation of histone H4 and recruitment of HSF1 to CTCF-rich loci that do not contain components of Pol II complex. Only the key regulatory factors discussed in the text are shown. ac, acetylation; CTCF, CCCTC-binding factor; Down, downregulated genes upon acute heat stress, dTRE, distal transcription regulatory element; dURE, distal untranscribed regulatory element; eRNA, enhancer RNA; GATA, GATA-binding protein; HSF1, heat shock factor 1; PIC, pre-initiation complex; P-TEFb, positive transcription elongation factor b; TAL1, T-cell acute lymphoid leukemia 1; TF, transcription factor, Up, upregulated genes or dTREs upon acute heat stress
Fig. 8 Rapid and coordinated reprogramming of genes and distal regulatory elements in stressed human cells. Model showing the rapid heat-induced response of the human genome, including reprogramming of gene transcription and establishment of a stress-specific repertoire of distal regulatory elements. Inhibition of the pause-release of promoter–proximal Pol II clears transcription complexes from the downregulated genes, elevating the concentration of free Pol II in the cell. The heat-induced genes are primed for directionality and rapid activation by the pre-assembled PIC on the core promoter of the coding strand, and by the paused Pol II at the 5′ region of the gene. Upon stress, trans-activators, such as HSF1, launch Pol II from the primed genes into productive elongation. The elevated levels of free Pol II and the highly positioned PIC on the coding strand enable instant loading of Pol II to the freed pause sites, efficiently launching rounds of transcript synthesis from the activated genes. The free Pol II allows also tuning up of the enhancer repertoire, increasing eRNA-production at dTREs that are marked by lineage-specific transcription factors, for example GATA1, GATA2 and TAL1 as shown in this study. In addition to reprogramming TREs, heat stress causes emergence of putative untranscribed regulatory elements, as demonstrated by the acetylation of histone H4 and recruitment of HSF1 to CTCF-rich loci that do not contain components of Pol II complex. Only the key regulatory factors discussed in the text are shown. ac, acetylation; CTCF, CCCTC-binding factor; Down, downregulated genes upon acute heat stress, dTRE, distal transcription regulatory element; dURE, distal untranscribed regulatory element; eRNA, enhancer RNA; GATA, GATA-binding protein; HSF1, heat shock factor 1; PIC, pre-initiation complex; P-TEFb, positive transcription elongation factor b; TAL1, T-cell acute lymphoid leukemia 1; TF, transcription factor, Up, upregulated genes or dTREs upon acute heat stressCC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

1 retracted

Show all 81 references

Cited by

Sorted by citations

See all 165

Related

Metadata

Authors
Vihervaara, A., Mahat, D. B., Guertin, M. J., Chu, T., Danko, C. G., Lis, J. T., & Sistonen, L.
Affiliations
  • Åbo Akademi University (FI)
  • Cornell University (US)
  • University of Virginia (US)
Corresponding author
  • Lis, John T.
  • Sistonen, Lea
ORCID
Venue
Nature Communications
Publisher
Springer Science and Business Media LLC
Published
2017-08-15
Volume
8
Issue
1
Pages
255-255
ISSN
2041-1723 (ISSN-L)
DOI
10.1038/s41467-017-00151-0
PMC
PMC5557961
OpenAlex
W2743018412
Cited by
164
References
81
Citations by year
2026: 13 · 2025: 19 · 2024: 9 · 2023: 14 · 2022: 24 · 2021: 29 · 2020: 25 · 2019: 15 · 2018: 14 · 2017: 1
FWCI
6.15
Type
journal article
Language
en
Open access
Licence
Locations
  • Nature Communications (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • RePEc: Research Papers in Economics (submittedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • Åbo Akademi University Research Portal (submittedVersion)
Funders
  • Magnus Ehrnroothin Säätiö
  • Sigrid Juséliuksen Säätiö
  • Åbo Akademi
  • National Institutes of Health (GM25232)
  • National Institute of General Medical Sciences (R01GM025232)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Genomics and Chromatin Dynamics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Chromatin Remodeling and Cancer
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Heat shock proteins research
Keywords
divergent transcription, promoter architecture, enhancer architecture, histone H4 acetylation, RNA polymerase pausing, heat shock factor 1, CTCF, chromatin organization, transcriptional reprogramming, heat shock response
MeSH terms
Acetylation, Cell Line, Enhancer Elements, Genetic, Gene Expression Regulation, Heat-Shock Proteins, Humans, Promoter Regions, Genetic, Stress, Physiological, Transcription Factors, Transcription, Genetic, Heat-Shock Response

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review