Mammalian Heat Shock Response and Mechanisms Underlying Its Genome-wide Transcriptional Regulation

A bi-stable feedback loop between GDNF, EGR1, and ERα contribute to endocrine resistant breast cancer

Discovering regulatory interactions between genes that specify the behavioral properties of cells remains an important challenge. We used the dynamics of transcriptional changes resolved by PRO-seq to identify a regulatory network responsible for endocrine resistance in breast cancer. We show that GDNF leads to endocrine resistance by switching the active state in a bi-stable feedback loop between GDNF, EGR1, and the master transcription factor ERα. GDNF stimulates MAP kinase, activating the transcription factors SRF and AP-1. SRF initiates an immediate transcriptional response, activating EGR1 and suppressing ERα. Newly translated EGR1 protein activates endogenous GDNF, leading to constitutive GDNF and EGR1 up-regulation, and the sustained down-regulation of ERα. Endocrine resistant MCF-7 cells are constitutively in the GDNF-high/ ERα-low state, suggesting that the state in the bi-stable feedback loop may provide a 'memory' of endocrine resistance. Thus, we identified a regulatory network switch that contributes to drug resistance in breast cancer.

ER-positive breast cancer cells are poised for RET-mediated endocrine resistance

The RET tyrosine kinase signaling pathway is involved in the development of endocrine resistant ER+ breast cancer. However, we know little about how ER+ cells activate RET signaling and initiate an endocrine resistant phenotype. Here we show that both ER+ endocrine resistant and sensitive breast cancers have a functional RET tyrosine kinase signaling pathway, but that endocrine sensitive breast cancer cells lack RET ligands that are necessary to drive endocrine resistance. Transcription of one RET ligand, GDNF, is necessary and sufficient to confer resistance in the ER+ MCF-7 cell line. Endogenous GDNF produced by endocrine resistant cells is translated, secreted into the media, and activates RET signaling in nearby cells. In patients, RET ligand expression predicts responsiveness to endocrine therapies and correlates with survival. Collectively, our findings show that ER+ tumor cells are "poised" for RET mediated endocrine resistance, expressing all components of the RET signaling pathway, but endocrine sensitive cells lack high expression of RET ligands that are necessary to initiate the resistance phenotype.

Deregulated Mucosal Immune Surveillance through Gut-Associated Regulatory T Cells and PD-1+ T Cells in Human Colorectal Cancer

Disturbed balance between immune surveillance and tolerance may lead to poor clinical outcomes in some malignancies. In paired analyses of adenocarcinoma and normal mucosa from 142 patients, we found a significant increase of the CD4/CD8 ratio and accumulation of regulatory T cells (Tregs) within the adenocarcinoma. The increased frequency of Tregs correlated with the local infiltration and extension of the tumor. There was concurrent maturation arrest, upregulation of programmed death-1 expression, and functional impairment in CD8⁺ T cells (CTLs) isolated from the adenocarcinoma. Adenocarcinoma-associated Tregs directly inhibit the function of normal human CTLs in vitro. With histopathological analysis, Foxp3⁺ Tregs were preferentially located in stroma. Concurrent transcriptome analysis of epithelial cells, stromal cells, and T cell subsets obtained from carcinomatous and normal intestinal samples from patients revealed a distinct gene expression signature in colorectal adenocarcinoma-associated Tregs, with overexpression of CCR1, CCR8, and TNFRSF9, whereas their ligands CCL4 and TNFSF9 were found upregulated in cancerous epithelium. Overexpression of WNT2 and CADM1, associated with carcinogenesis and metastasis, in cancer-associated stromal cells suggests that both cancer cells and stromal cells play important roles in the development and progression of colorectal cancer through the formation of a tumor microenvironment. The identification of CTL anergy by Tregs and the unique gene expression signature of human Tregs and stromal cells in colorectal cancer patients may facilitate the development of new therapeutics against malignancies.

TimeLapse-seq: adding a temporal dimension to RNA sequencing through nucleoside recoding

RNA sequencing (RNA-seq) offers a snapshot of cellular RNA populations, but not temporal information about the sequenced RNA. Here we report TimeLapse-seq, which uses oxidative-nucleophilic-aromatic substitution to convert 4-thiouridine into cytidine analogs, yielding apparent U-to-C mutations that mark new transcripts upon sequencing. TimeLapse-seq is a single-molecule approach that is adaptable to many applications and reveals RNA dynamics and induced differential expression concealed in traditional RNA-seq.

Cdk9 regulates a promoter-proximal checkpoint to modulate RNA polymerase II elongation rate in fission yeast

Post-translational modifications of the transcription elongation complex provide mechanisms to fine-tune gene expression, yet their specific impacts on RNA polymerase II regulation remain difficult to ascertain. Here, in Schizosaccharomyces pombe, we examine the role of Cdk9, and related Mcs6/Cdk7 and Lsk1/Cdk12 kinases, on transcription at base-pair resolution with Precision Run-On sequencing (PRO-seq). Within a minute of Cdk9 inhibition, phosphorylation of Pol II-associated factor, Spt5 is undetectable. The effects of Cdk9 inhibition are more severe than inhibition of Cdk7 and Cdk12, resulting in a shift of Pol II toward the transcription start site (TSS). A time course of Cdk9 inhibition reveals that early transcribing Pol II can escape promoter-proximal regions, but with a severely reduced elongation rate of only ~400 bp/min. Our results in fission yeast suggest the existence of a conserved global regulatory checkpoint that requires Cdk9 kinase activity.

Hsf1 and Hsp70 constitute a two-component feedback loop that regulates the yeast heat shock response

Models for regulation of the eukaryotic heat shock response typically invoke a negative feedback loop consisting of the transcriptional activator Hsf1 and a molecular chaperone. Previously we identified Hsp70 as the chaperone responsible for Hsf1 repression and constructed a mathematical model that recapitulated the yeast heat shock response (Zheng et al., 2016). The model was based on two assumptions: dissociation of Hsp70 activates Hsf1, and transcriptional induction of Hsp70 deactivates Hsf1. Here we validate these assumptions. First, we severed the feedback loop by uncoupling Hsp70 expression from Hsf1 regulation. As predicted by the model, Hsf1 was unable to efficiently deactivate in the absence of Hsp70 transcriptional induction. Next, we mapped a discrete Hsp70 binding site on Hsf1 to a C-terminal segment known as conserved element 2 (CE2). In vitro, CE2 binds to Hsp70 with low affinity (9 µM), in agreement with model requirements. In cells, removal of CE2 resulted in increased basal Hsf1 activity and delayed deactivation during heat shock, while tandem repeats of CE2 sped up Hsf1 deactivation. Finally, we uncovered a role for the N-terminal domain of Hsf1 in negatively regulating DNA binding. These results reveal the quantitative control mechanisms underlying the heat shock response.

Acute Heat Stress Alters the Expression of Orexin System in Quail Muscle

Accumulating evidences indicate that the hypothalamic neuropeptide orexins or hypocretins are involved in stress-induced responses in mammals. Recently, we found that orexin is expressed and secreted from avian muscle cells, however its regulation is still unknown. In this study, we investigated the effect of heat and oxidative stress, the most challenging stressors in poultry production, on the expression of orexin system in quail muscle tissues and myoblast cell lines. Four week-old genetically selected susceptible and resistant Japanese quail (Coturnix coturnix Japonica) lines were exposed to acute heat stress (HS, 37°C for 1.5 h) or maintained at thermoneutral conditions (24°C). Quail myoblast (QM7) cell line was exposed to heat stress (45°C) for 0.5, 1, 2, or 4 h. The control cells were maintained at 37°C. The cells were also treated with several doses of hydrogen peroxide (H₂O₂, 10-200 μM) or 4-Hydroxynonenal (4-HNE, 10-30 μM) as oxidative stress. Untreated cells were used as controls. Acute HS significantly induced the expression of HSP70 and down-regulated orexin system in both quail muscle tissue and QM7 cells. Similarly, H₂O₂ but not 4-HNE treatment significantly increased HSP70 protein levels and dysregulated the expression of orexin and its related receptors in a dose-dependent manner in QM7 cells. Transient overexpression of HSP70 down-regulated the expression of orexin system in QM7 cells. Taken together, these data indicate that orexin may be a key player in stress response in avian muscle by demonstrating that heat and oxidative stress alter the expression of orexin system in quail muscle. This effect might be mediated through HSP70. Unraveling the upstream regulators and downstream effectors of orexin in avian muscle merits further in depth investigations.

The role of declining adaptive homeostasis in ageing

Adaptive homeostasis is "the transient expansion or contraction of the homeostatic range for any given physiological parameter in response to exposure to sub-toxic, non-damaging, signalling molecules or events, or the removal or cessation of such molecules or events" (Davies, 2016). Adaptive homeostasis enables biological systems to make continuous short-term adjustments for optimal functioning despite ever-changing internal and external environments. Initiation of adaptation in response to an appropriate signal allows organisms to successfully cope with much greater, normally toxic, stresses. These short-term responses are initiated following effective signals, including hypoxia, cold shock, heat shock, oxidative stress, exercise-induced adaptation, caloric restriction, osmotic stress, mechanical stress, immune response, and even emotional stress. There is now substantial literature detailing a decline in adaptive homeostasis that, unfortunately, appears to manifest with ageing, especially in the last third of the lifespan. In this review, we present the hypothesis that one hallmark of the ageing process is a significant decline in adaptive homeostasis capacity. We discuss the mechanistic importance of diminished capacity for short-term (reversible) adaptive responses (both biochemical and signal transduction/gene expression-based) to changing internal and external conditions, for short-term survival and for lifespan and healthspan. Studies of cultured mammalian cells, worms, flies, rodents, simians, apes, and even humans, all indicate declining adaptive homeostasis as a potential contributor to age-dependent senescence, increased risk of disease, and even mortality. Emerging work points to Nrf2-Keap1 signal transduction pathway inhibitors, including Bach1 and c-Myc, both of whose tissue concentrations increase with age, as possible major causes for age-dependent loss of adaptive homeostasis.

Fibronectin is a stress responsive gene regulated by HSF1 in response to geldanamycin

Fibronectin is an extracellular matrix glycoprotein with key roles in cell adhesion and migration. Hsp90 binds directly to fibronectin and Hsp90 depletion regulates fibronectin matrix stability. Where inhibition of Hsp90 with a C-terminal inhibitor, novobiocin, reduced the fibronectin matrix, treatment with an N-terminal inhibitor, geldanamycin, increased fibronectin levels. Geldanamycin treatment induced a stress response and a strong dose and time dependent increase in fibronectin mRNA via activation of the fibronectin promoter. Three putative heat shock elements (HSEs) were identified in the fibronectin promoter. Loss of two of these HSEs reduced both basal and geldanamycin-induced promoter activity, as did inhibition of the stress-responsive transcription factor HSF1. Binding of HSF1 to one of the putative HSE was confirmed by ChIP under basal conditions, and occupancy shown to increase with geldanamycin treatment. These data support the hypothesis that fibronectin is stress-responsive and a functional HSF1 target gene. COLA42 and LAMB3 mRNA levels were also increased with geldanamycin indicating that regulation of extracellular matrix (ECM) genes by HSF1 may be a wider phenomenon. Taken together, these data have implications for our understanding of ECM dynamics in stress-related diseases in which HSF1 is activated, and where the clinical application of N-terminal Hsp90 inhibitors is intended.

Global unleashing of transcription elongation waves in response to genotoxic stress restricts somatic mutation rate

Complex molecular responses preserve gene expression accuracy and genome integrity in the face of environmental perturbations. Here we report that, in response to UV irradiation, RNA polymerase II (RNAPII) molecules are dynamically and synchronously released from promoter-proximal regions into elongation to promote uniform and accelerated surveillance of the whole transcribed genome. The maximised influx of de novo released RNAPII correlates with increased damage-sensing, as confirmed by RNAPII progressive accumulation at dipyrimidine sites and by the average slow-down of elongation rates in gene bodies. In turn, this transcription elongation ‘safe’ mode guarantees efficient DNA repair regardless of damage location, gene size and transcription level. Accordingly, we detect low and homogenous rates of mutational signatures associated with UV exposure or cigarette smoke across all active genes. Our study reveals a novel advantage for transcription regulation at the promoter-proximal level and provides unanticipated insights into how active transcription shapes the mutagenic landscape of cancer genomes.

Precise orchestration of gene expression regulation upon DNA damage is essential for genome integrity. Here the authors identify a novel widespread stress-triggered defence mechanism that promotes rapid transcription-driven genomic surveillance thus limiting mutagenesis and shaping cancer genomes.

Molecular Genetic Analysis of Human Endometrial Mesenchymal Stem Cells That Survived Sublethal Heat Shock

High temperature is a critical environmental and personal factor. Although heat shock is a well-studied biological phenomenon, hyperthermia response of stem cells is poorly understood. Previously, we demonstrated that sublethal heat shock induced premature senescence in human endometrial mesenchymal stem cells (eMSC). This study aimed to investigate the fate of eMSC-survived sublethal heat shock (SHS) with special emphasis on their genetic stability and possible malignant transformation using methods of classic and molecular karyotyping, next-generation sequencing, and transcriptome functional analysis. G-banding revealed random chromosome breakages and aneuploidy in the SHS-treated eMSC. Molecular karyotyping found no genomic imbalance in these cells. Gene module and protein interaction network analysis of mRNA sequencing data showed that compared to untreated cells, SHS-survived progeny revealed some difference in gene expression. However, no hallmarks of cancer were found. Our data identified downregulation of oncogenic signaling, upregulation of tumor-suppressing and prosenescence signaling, induction of mismatch, and excision DNA repair. The common feature of heated eMSC is the silence of MYC, AKT1/PKB oncogenes, and hTERT telomerase. Overall, our data indicate that despite genetic instability, SHS-survived eMSC do not undergo transformation. After long-term cultivation, these cells like their unheated counterparts enter replicative senescence and die.

The HSF1-PARP13-PARP1 complex facilitates DNA repair and promotes mammary tumorigenesis

Poly(ADP-ribose) polymerase 1 (PARP1) is involved in DNA repair, chromatin structure, and transcription. However, the mechanisms that regulate PARP1 distribution on DNA are poorly understood. Here, we show that heat shock transcription factor 1 (HSF1) recruits PARP1 through the scaffold protein PARP13. In response to DNA damage, activated and auto-poly-ADP-ribosylated PARP1 dissociates from HSF1–PARP13, and redistributes to DNA lesions and DNA damage-inducible gene loci. Histone deacetylase 1 maintains PARP1 in the ternary complex by inactivating PARP1 through deacetylation. Blocking ternary complex formation impairs redistribution of PARP1 during DNA damage, which reduces gene expression and DNA repair. Furthermore, ternary complex formation and PARP1 redistribution protect cells from DNA damage by promoting DNA repair, and support growth of BRCA1-null mammary tumors, which are sensitive to PARP inhibitors. Our findings identify HSF1 as a regulator of genome integrity and define this function as a guarding mechanism for a specific type of mammary tumorigenesis.

PARP1 recruitment to DNA lesions is critical for DNA damage repair but how DNA damage induces PARP1 redistribution is largely unknown. Here, the authors provide evidence that PARP1 redistribution and DNA repair in tumor cells depend on the formation of a HSF1–PARP13–PARP1 complex.

Nascent RNA sequencing reveals a dynamic global transcriptional response at genes and enhancers to the natural medicinal compound celastrol

Most studies of responses to transcriptional stimuli measure changes in cellular mRNA concentrations. By sequencing nascent RNA instead, it is possible to detect changes in transcription in minutes rather than hours and thereby distinguish primary from secondary responses to regulatory signals. Here, we describe the use of PRO-seq to characterize the immediate transcriptional response in human cells to celastrol, a compound derived from traditional Chinese medicine that has potent anti-inflammatory, tumor-inhibitory, and obesity-controlling effects. Celastrol is known to elicit a cellular stress response resembling the response to heat shock, but the transcriptional basis of this response remains unclear. Our analysis of PRO-seq data for K562 cells reveals dramatic transcriptional effects soon after celastrol treatment at a broad collection of both coding and noncoding transcription units. This transcriptional response occurred in two major waves, one within 10 min, and a second 40-60 min after treatment. Transcriptional activity was generally repressed by celastrol, but one distinct group of genes, enriched for roles in the heat shock response, displayed strong activation. Using a regression approach, we identified key transcription factors that appear to drive these transcriptional responses, including members of the E2F and RFX families. We also found sequence-based evidence that particular transcription factors drive the activation of enhancers. We observed increased polymerase pausing at both genes and enhancers, suggesting that pause release may be widely inhibited during the celastrol response. Our study demonstrates that a careful analysis of PRO-seq time-course data can disentangle key aspects of a complex transcriptional response, and it provides new insights into the activity of a powerful pharmacological agent.

Heat Shock Factor 1: From Fire Chief to Crowd-Control Specialist

HSF1 is the supposed master regulator of the heat shock response. In this issue of Molecular Cell, Solís et al. reveal that it has a much narrower job description: organizing a small team of molecular chaperones that keep the proteome moving.

Transcriptional Regulation of the Ambient Temperature Response by H2A.Z Nucleosomes and HSF1 Transcription Factors in Arabidopsis

Temperature influences the distribution, range, and phenology of plants. The key transcriptional activators of heat shock response in eukaryotes, the heat shock factors (HSFs), have undergone large-scale gene amplification in plants. While HSFs are central in heat stress responses, their role in the response to ambient temperature changes is less well understood. We show here that the warm ambient temperature transcriptome is dependent upon the HSFA1 clade of Arabidopsis HSFs, which cause a rapid and dynamic eviction of H2A.Z nucleosomes at target genes. A transcriptional cascade results in the activation of multiple downstream stress-responsive transcription factors, triggering large-scale changes to the transcriptome in response to elevated temperature. H2A.Z nucleosomes are enriched at temperature-responsive genes at non-inducible temperature, and thus likely confer inducibility of gene expression and higher responsive dynamics. We propose that the antagonistic effects of H2A.Z and HSF1 provide a mechanism to activate gene expression rapidly and precisely in response to temperature, while preventing leaky transcription in the absence of an activation signal.

HSF1-dependent and -independent regulation of the mammalian in vivo heat shock response and its impairment in Huntington's disease mouse models

The heat shock response (HSR) is a mechanism to cope with proteotoxic stress by inducing the expression of molecular chaperones and other heat shock response genes. The HSR is evolutionarily well conserved and has been widely studied in bacteria, cell lines and lower eukaryotic model organisms. However, mechanistic insights into the HSR in higher eukaryotes, in particular in mammals, are limited. We have developed an in vivo heat shock protocol to analyze the HSR in mice and dissected heat shock factor 1 (HSF1)-dependent and -independent pathways. Whilst the induction of proteostasis-related genes was dependent on HSF1, the regulation of circadian function related genes, indicating that the circadian clock oscillators have been reset, was independent of its presence. Furthermore, we demonstrate that the in vivo HSR is impaired in mouse models of Huntington's disease but we were unable to corroborate the general repression of transcription that follows a heat shock in lower eukaryotes.

New inhibitor targeting human transcription factor HSF1: effects on the heat shock response and tumor cell survival

Comparative modeling of the DNA-binding domain of human HSF1 facilitated the prediction of possible binding pockets for small molecules and definition of corresponding pharmacophores. In silico screening of a large library of lead-like compounds identified a set of compounds that satisfied the pharmacophoric criteria, a selection of which compounds was purchased to populate a biased sublibrary. A discriminating cell-based screening assay identified compound 001, which was subjected to systematic analysis of structure–activity relationships, resulting in the development of compound 115 (I_HSF115). I_HSF115 bound to an isolated HSF1 DNA-binding domain fragment. The compound did not affect heat-induced oligomerization, nuclear localization and specific DNA binding but inhibited the transcriptional activity of human HSF1, interfering with the assembly of ATF1-containing transcription complexes. I_HSF115 was employed to probe the human heat shock response at the transcriptome level. In contrast to earlier studies of differential regulation in HSF1-naïve and -depleted cells, our results suggest that a large majority of heat-induced genes is positively regulated by HSF1. That I_HSF115 effectively countermanded repression in a significant fraction of heat-repressed genes suggests that repression of these genes is mediated by transcriptionally active HSF1. I_HSF115 is cytotoxic for a variety of human cancer cell lines, multiple myeloma lines consistently exhibiting high sensitivity.

Comparative analysis reveals genomic features of stress-induced transcriptional readthrough

Transcription is a highly regulated process, and stress-induced changes in gene transcription have been shown to play a major role in stress responses and adaptation. Genome-wide studies reveal prevalent transcription beyond known protein-coding gene loci, generating a variety of RNA classes, most of unknown function. One such class, termed downstream of gene-containing transcripts (DoGs), was reported to result from transcriptional readthrough upon osmotic stress in human cells. However, how widespread the readthrough phenomenon is, and what its causes and consequences are, remain elusive. Here we present a genome-wide mapping of transcriptional readthrough, using nuclear RNA-Seq, comparing heat shock, osmotic stress, and oxidative stress in NIH 3T3 mouse fibroblast cells. We observe massive induction of transcriptional readthrough, both in levels and length, under all stress conditions, with significant, yet not complete, overlap of readthrough-induced loci between different conditions. Importantly, our analyses suggest that stress-induced transcriptional readthrough is not a random failure process, but is rather differentially induced across different conditions. We explore potential regulators and find a role for HSF1 in the induction of a subset of heat shock-induced readthrough transcripts. Analysis of public datasets detected increases in polymerase II occupancy in DoG regions after heat shock, supporting our findings. Interestingly, DoGs tend to be produced in the vicinity of neighboring genes, leading to a marked increase in their antisense-generating potential. Finally, we examine genomic features of readthrough transcription and observe a unique chromatin signature typical of DoG-producing regions, suggesting that readthrough transcription is associated with the maintenance of an open chromatin state.

Rethinking HSF1 in Stress, Development, and Organismal Health

The heat shock response (HSR) was originally discovered as a transcriptional response to elevated temperature shock and led to the identification of heat shock proteins and heat shock factor 1 (HSF1). Since then HSF1 has been shown to be important for combating other forms of environmental perturbations as well as genetic variations that cause proteotoxic stress. The HSR has long been thought to be an absolute response to conditions of cell stress and the primary mechanism by which HSF1 promotes organismal health by preventing protein aggregation and subsequent proteome imbalance. Accumulating evidence now shows that HSF1, the central player in the HSR, is regulated according to specific cellular requirements through cell-autonomous and non-autonomous signals, and directs transcriptional programs distinct from the HSR during development and in carcinogenesis. We discuss here these 'non-canonical' roles of HSF1, its regulation in diverse conditions of development, reproduction, metabolism, and aging, and posit that HSF1 serves to integrate diverse biological and pathological responses.

Chaperone co-inducer BGP-15 inhibits histone deacetylases and enhances the heat shock response through increased chromatin accessibility

Defects in cellular protein homeostasis are associated with many severe and prevalent pathological conditions such as neurodegenerative diseases, muscle dystrophies, and metabolic disorders. One way to counteract these defects is to improve the protein homeostasis capacity through induction of the heat shock response. Despite numerous attempts to develop strategies for chemical activation of the heat shock response by heat shock transcription factor 1 (HSF1), the underlying mechanisms of drug candidates' mode of action are poorly understood. To lower the threshold for the heat shock response activation, we used the chaperone co-inducer BGP-15 that was previously shown to have beneficial effects on several proteinopathic disease models. We found that BGP-15 treatment combined with heat stress caused a substantial increase in HSF1-dependent heat shock protein 70 (HSPA1A/B) expression already at a febrile range of temperatures. Moreover, BGP-15 alone inhibited the activity of histone deacetylases (HDACs), thereby increasing chromatin accessibility at multiple genomic loci including the stress-inducible HSPA1A. Intriguingly, treatment with well-known potent HDAC inhibitors trichostatin A and valproic acid enhanced the heat shock response and improved cytoprotection. These results present a new pharmacological strategy for restoring protein homeostasis by inhibiting HDACs, increasing chromatin accessibility, and lowering the threshold for heat shock response activation.

PAF1 regulation of promoter-proximal pause release via enhancer activation

Gene expression in metazoans is regulated by RNA polymerase II (Pol II) promoter-proximal pausing and its release. Previously, we showed that Pol II-associated factor 1 (PAF1) modulates the release of paused Pol II into productive elongation. Here, we found that PAF1 occupies transcriptional enhancers and restrains hyperactivation of a subset of these enhancers. Enhancer activation as the result of PAF1 loss releases Pol II from paused promoters of nearby PAF1 target genes. Knockout of PAF1-regulated enhancers attenuates the release of paused Pol II on PAF1 target genes without major interference in the establishment of pausing at their cognate promoters. Thus, a subset of enhancers can primarily modulate gene expression by controlling the release of paused Pol II in a PAF1-dependent manner.

A Model of Exposure to Extreme Environmental Heat Uncovers the Human Transcriptome to Heat Stress

The molecular mechanisms by which individuals subjected to environmental heat stress either recover or develop heat-related complications are not well understood. We analysed the changes in blood mononuclear gene expression patterns in human volunteers exposed to extreme heat in a sauna (temperature of 75.7 ± 0.86 °C). Our analysis reveals that expression changes occur rapidly with no significant increase in core temperature and continue to amplify one hour after the end of heat stress. The reprogramed transcriptome was predominantly inhibitory, as more than two-thirds of the expressed genes were suppressed. The differentially expressed genes encoded proteins that function in stress-associated pathways; including proteostasis, energy metabolism, cell growth and proliferation, and cell death, and survival. The transcriptome also included mitochondrial dysfunction, altered protein synthesis, and reduced expression of genes -related to immune function. The findings reveal the human transcriptomic response to heat and highlight changes that might underlie the health outcomes observed during heat waves.

The heat shock response and humoral immune response are mutually antagonistic in honey bees

The honey bee is of paramount importance to humans in both agricultural and ecological settings. Honey bee colonies have suffered from increased attrition in recent years, stemming from complex interacting stresses. Defining common cellular stress responses elicited by these stressors represents a key step in understanding potential synergies. The proteostasis network is a highly conserved network of cellular stress responses involved in maintaining the homeostasis of protein production and function. Here, we have characterized the Heat Shock Response (HSR), one branch of this network, and found that its core components are conserved. In addition, exposing bees to elevated temperatures normally encountered by honey bees during typical activities results in robust HSR induction with increased expression of specific heat shock proteins that was variable across tissues. Surprisingly, we found that heat shock represses multiple immune genes in the abdomen and additionally showed that wounding the cuticle of the abdomen results in decreased expression of multiple HSR genes in proximal and distal tissues. This mutually antagonistic relationship between the HSR and immune activation is unique among invertebrates studied to date and may promote understanding of potential synergistic effects of disparate stresses in this critical pollinator and social insects more broadly.

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

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.Heat Shock Factor 1 (HSF1) is a regulator of stress-induced transcription. Here, the authors investigate changes to transcription and chromatin organization upon stress and find that activated HSF1 binds to transcription-primed promoters and enhancers, and to CTCF occupied, untranscribed chromatin.

The Aging Epigenome

During aging, the mechanisms that normally maintain health and stress resistance strikingly decline, resulting in decrepitude, frailty, and ultimately death. Exactly when and how this decline occurs is unknown. Changes in transcriptional networks and chromatin state lie at the heart of age-dependent decline. These epigenomic changes are not only observed during aging but also profoundly affect cellular function and stress resistance, thereby contributing to the progression of aging. We propose that the dysregulation of transcriptional and chromatin networks is a crucial component of aging. Understanding age-dependent epigenomic changes will yield key insights into how aging begins and progresses and should lead to the development of new therapeutics that delay or even reverse aging and age-related diseases.

Multisubunit DNA-Dependent RNA Polymerases from Vaccinia Virus and Other Nucleocytoplasmic Large-DNA Viruses: Impressions from the Age of Structure

The past 17 years have been marked by a revolution in our understanding of cellular multisubunit DNA-dependent RNA polymerases (MSDDRPs) at the structural level. A parallel development over the past 15 years has been the emerging story of the giant viruses, which encode MSDDRPs. Here we link the two in an attempt to understand the specialization of multisubunit RNA polymerases in the domain of life encompassing the large nucleocytoplasmic DNA viruses (NCLDV), a superclade that includes the giant viruses and the biochemically well-characterized poxvirus vaccinia virus. The first half of this review surveys the recently determined structural biology of cellular RNA polymerases for a microbiology readership. The second half discusses a reannotation of MSDDRP subunits from NCLDV families and the apparent specialization of these enzymes by virus family and by subunit with regard to subunit or domain loss, subunit dissociability, endogenous control of polymerase arrest, and the elimination/customization of regulatory interactions that would confer higher-order cellular control. Some themes are apparent in linking subunit function to structure in the viral world: as with cellular RNA polymerases I and III and unlike cellular RNA polymerase II, the viral enzymes seem to opt for speed and processivity and seem to have eliminated domains associated with higher-order regulation. The adoption/loss of viral RNA polymerase proofreading functions may have played a part in matching intrinsic mutability to genome size.

Genome-wide characterization of mammalian promoters with distal enhancer functions

Gene expression in mammals is precisely regulated by the combination of promoters and gene-distal regulatory regions, known as enhancers. Several studies have suggested that some promoters might have enhancer functions. However, the extent of this type of promoters and whether they actually function to regulate the expression of distal genes have remained elusive. Here, by exploiting a high-throughput enhancer reporter assay, we unravel a set of mammalian promoters displaying enhancer activity. These promoters have distinct genomic and epigenomic features and frequently interact with other gene promoters. Extensive CRISPR-Cas9 genomic manipulation demonstrated the involvement of these promoters in the cis regulation of expression of distal genes in their natural loci. Our results have important implications for the understanding of complex gene regulation in normal development and disease.

Conserved Transcription Factors Steer Growth-Related Genomic Programs in Daphnia

Ecological genomics aims to understand the functional association between environmental gradients and the genes underlying adaptive traits. Many genes that are identified by genome-wide screening in ecologically relevant species lack functional annotations. Although gene functions can be inferred from sequence homology, such approaches have limited power. Here, we introduce ecological regulatory genomics by presenting an ontology-free gene prioritization method. Specifically, our method combines transcriptome profiling with high-throughput cis-regulatory sequence analysis in the water fleas Daphnia pulex and Daphnia magna. It screens coexpressed genes for overrepresented DNA motifs that serve as transcription factor binding sites, thereby providing insight into conserved transcription factors and gene regulatory networks shaping the expression profile. We first validated our method, called Daphnia-cisTarget, on a D. pulex heat shock data set, which revealed a network driven by the heat shock factor. Next, we performed RNA-Seq in D. magna exposed to the cyanobacterium Microcystis aeruginosa. Daphnia-cisTarget identified coregulated gene networks that associate with the moulting cycle and potentially regulate life history changes in growth rate and age at maturity. These networks are predicted to be regulated by evolutionary conserved transcription factors such as the homologues of Drosophila Shavenbaby and Grainyhead, nuclear receptors, and a GATA family member. In conclusion, our approach allows prioritising candidate genes in Daphnia without bias towards prior knowledge about functional gene annotation and represents an important step towards exploring the molecular mechanisms of ecological responses in organisms with poorly annotated genomes.

Systematic comparison of the response properties of protein and RNA mediated gene regulatory motifs

We present a framework enabling the dissection of the effects of motif structure (feedback or feedforward), the nature of the controller (RNA or protein), and the regulation mode (transcriptional, post-transcriptional or translational) on the response to a step change in the input. We have used a common model framework for gene expression where both motif structures have an activating input and repressing regulator, with the same set of parameters, to enable a comparison of the responses. We studied the global sensitivity of the system properties, such as steady-state gain, overshoot, peak time, and peak duration, to parameters. We find that, in all motifs, overshoot correlated negatively whereas peak duration varied concavely with peak time. Differences in the other system properties were found to be mainly dependent on the nature of the controller rather than the motif structure. Protein mediated motifs showed a higher degree of adaptation i.e. a tendency to return to baseline levels; in particular, feedforward motifs exhibited perfect adaptation. RNA mediated motifs had a mild regulatory effect; they also exhibited a lower peaking tendency and mean overshoot. Protein mediated feedforward motifs showed higher overshoot and lower peak time compared to the corresponding feedback motifs.

Destabilization of B2 RNA by EZH2 Activates the Stress Response

More than 98% of the mammalian genome is noncoding, and interspersed transposable elements account for ∼50% of noncoding space. Here, we demonstrate that a specific interaction between the polycomb protein EZH2 and RNA made from B2 SINE retrotransposons controls stress-responsive genes in mouse cells. In the heat-shock model, B2 RNA binds stress genes and suppresses their transcription. Upon stress, EZH2 is recruited and triggers cleavage of B2 RNA. B2 degradation in turn upregulates stress genes. Evidence indicates that B2 RNA operates as a "speed bump" against advancement of RNA polymerase II, and temperature stress releases the brakes on transcriptional elongation. These data attribute a new function to EZH2 that is independent of its histone methyltransferase activity and reconcile how EZH2 can be associated with both gene repression and activation. Our study reveals that EZH2 and B2 together control activation of a large network of genes involved in thermal stress.

Chromatin SUMOylation in heat stress: To protect, pause and organise?: SUMO stress response on chromatin

Post-translational modifications, e.g. SUMO modifications (SUMOylation), provide a mechanism for swiftly changing a protein's activity. Various stress conditions trigger a SUMO stress response (SSR) - a stress-induced rapid change in the conjugation of SUMO to multiple proteins, which predominantly targets nuclear proteins. The SSR has been postulated to protect stressed cells by preserving the functionality of crucial proteins. However, it is unclear how it exerts its protective functions. Interestingly, heat stress (HS) increases SUMOylation of proteins at active promoters and enhancers. In promoters, HS-induced SUMOylation correlates with gene transcription and stress-induced RNA polymerase II (Pol2) pausing. Conversely, a disappearance of SUMOylation in HS occurs at chromatin anchor points that maintain chromatin-looping structures and the spatial organisation of chromatin. In reviewing the literature, we hypothesise that the SSR regulates Pol2 pausing by modulating the interactions of pausing-regulating proteins, whereas deSUMOylation alters the function of chromatin anchors.

The origins of developmental gene regulation

The leap from simple unicellularity to complex multicellularity remains one of life's major enigmas. The origins of metazoan developmental gene regulatory mechanisms are sought by analyzing gene regulation in extant eumetazoans, sponges, and unicellular organisms. The main hypothesis of this manuscript is that, developmental enhancers evolved from unicellular inducible promoters that diversified the expression of regulatory genes during metazoan evolution. Promoters and enhancers are functionally similar; both can regulate the transcription of distal promoters and both direct local transcription. Additionally, enhancers have experimentally characterized structural features that reveal their origin from inducible promoters. The distal co-operative regulation among promoters identified in unicellular opisthokonts possibly represents the precursor of distal regulation of promoters by enhancers. During metazoan evolution, constitutive-type promoters of regulatory genes would have acquired novel receptivity to distal regulatory inputs from promoters of inducible genes that eventually specialized as enhancers. The novel regulatory interactions would have caused constitutively expressed genes controlling differential gene expression in unicellular organisms to become themselves differentially expressed. The consequence of the novel regulatory interactions was that regulatory pathways of unicellular organisms became interlaced and ultimately evolved into the intricate developmental gene regulatory networks (GRNs) of extant metazoans.

Use of conditioned media is critical for studies of regulation in response to rapid heat shock

Heat shock response (HSR) maintains and restores protein homeostasis when cells are exposed to proteotoxic heat stress. Heat shock (HS) triggers a rapid and robust change in genome-wide transcription, protein synthesis, and chaperone activity; and therefore, the HSR has been widely used as a model system in these studies. The conventional method of performing instantaneous HS in the laboratory uses heated fresh media to induce HSR when added to cells. However, addition of fresh media to cells may evoke additional cellular responses and signaling pathways. Here, we compared the change in global transcription profile when HS is performed with either heated fresh media or heated conditioned media. We found that the use of heated fresh media induces transcription of hundreds of genes that HS alone does not induce, and masks or partially masks HS-mediated downregulation of thousands of genes. The fresh-media-dependent upregulated genes encode ribosomal subunit proteins involved in translation and RNA processing factors. More importantly, fresh media also induce transcription of several heat shock protein genes (Hsps) in a heat shock factor 1 (HSF1)-independent manner. Thus, we conclude that a conventional method of HS with heated fresh media causes changes in transcription regulation that confound the actual change caused solely by elevated temperature of cells.

Dynamic control of Hsf1 during heat shock by a chaperone switch and phosphorylation

Heat shock factor (Hsf1) regulates the expression of molecular chaperones to maintain protein homeostasis. Despite its central role in stress resistance, disease and aging, the mechanisms that control Hsf1 activity remain unresolved. Here we show that in budding yeast, Hsf1 basally associates with the chaperone Hsp70 and this association is transiently disrupted by heat shock, providing the first evidence that a chaperone repressor directly regulates Hsf1 activity. We develop and experimentally validate a mathematical model of Hsf1 activation by heat shock in which unfolded proteins compete with Hsf1 for binding to Hsp70. Surprisingly, we find that Hsf1 phosphorylation, previously thought to be required for activation, in fact only positively tunes Hsf1 and does so without affecting Hsp70 binding. Our work reveals two uncoupled forms of regulation - an ON/OFF chaperone switch and a tunable phosphorylation gain - that allow Hsf1 to flexibly integrate signals from the proteostasis network and cell signaling pathways.

DOI:http://dx.doi.org/10.7554/eLife.18638.001

Proteins are strings of amino acids that carry out crucial activities inside cells, such as harvesting energy and generating the building blocks that cells need to grow. In order to carry out their specific roles inside the cell, the proteins need to “fold” into precise three-dimensional shapes.

Protein folding is critical for life, and cells don’t leave it up to chance. Cells employ “molecular chaperones” to help proteins to fold properly. However, under some conditions – such as high temperature – proteins are more difficult to fold and the chaperones can become overwhelmed. In these cases, unfolded proteins can pile up in the cell. This leads not only to the cell being unable to work properly, but also to the formation of toxic “aggregates”. These aggregates are tangles of unfolded proteins that are hallmarks of many neurodegenerative diseases such as Alzheimer’s, Parkinson’s and amyotrophic lateral sclerosis (ALS).

Protein aggregates can be triggered by high temperature in a condition termed “heat shock”. A sensor named heat shock factor 1 (Hsf1 for short) increases the amount of chaperones following heat shock. But what controls the activity of Hsf1?

To answer this question, Zheng, Krakowiak et al. combined mathematical modelling and experiments in yeast cells. The most important finding is that the ‘on/off switch’ that controls Hsf1 is based on whether Hsf1 is itself bound to a chaperone. When bound to the chaperone, Hsf1 is turned ‘off’; when the chaperone falls off, Hsf1 turns ‘on’ and makes more chaperones; when there are enough chaperones, they once again bind to Hsf1 and turn it back ‘off’. In this way, Hsf1 and the chaperones form a feedback loop that ensures that there are always enough chaperones to keep the cell’s proteins folded.

Now that we know how Hsf1 is controlled, can we harness this understanding to tune the activity of Hsf1 without disrupting how the chaperones work? If we can activate Hsf1, we can provide cells with more chaperones. This could be a therapeutic strategy to combat neurodegenerative diseases.

DOI:http://dx.doi.org/10.7554/eLife.18638.002

High-Resolution Mapping of RNA Polymerases Identifies Mechanisms of Sensitivity and Resistance to BET Inhibitors in t(8;21) AML

Bromodomain and extra-terminal domain (BET) family inhibitors offer an approach to treating hematological malignancies. We used precision nuclear run-on transcription sequencing (PRO-seq) to create high-resolution maps of active RNA polymerases across the genome in t(8;21) acute myeloid leukemia (AML), as these polymerases are exceptionally sensitive to BET inhibitors. PRO-seq identified over 1,400 genes showing impaired release of promoter-proximal paused RNA polymerases, including the stem cell factor receptor tyrosine kinase KIT that is mutated in t(8;21) AML. PRO-seq also identified an enhancer 3' to KIT. Chromosome conformation capture confirmed contacts between this enhancer and the KIT promoter, while CRISPRi-mediated repression of this enhancer impaired cell growth. PRO-seq also identified microRNAs, including MIR29C and MIR29B2, that target the anti-apoptotic factor MCL1 and were repressed by BET inhibitors. MCL1 protein was upregulated, and inhibition of BET proteins sensitized t(8:21)-containing cells to MCL1 inhibition, suggesting a potential mechanism of resistance to BET-inhibitor-induced cell death.

Size doesn't matter in the heat shock response

Heat shock factor 1 (Hsf1) is a transcription factor that is often described as the master regulator of the heat shock response in all eukaryotes. However, due to its essentiality in yeast, Hsf1's contribution to the transcriptome under basal and heat shock conditions has never been directly determined. Using a chemical genetics approach that allowed rapid Hsf1 inactivation, my colleagues and I have recently shown that the bulk of the heat shock response is Hsf1 independent. Rather than inducing genes responsible for carrying out the various cellular processes required for adaptation to thermal stress, Hsf1 controls a dedicated set of chaperone protein genes devoted to restoring protein-folding homeostasis. The limited scope of the Hsf1 regulon belies its outsize importance in cellular fitness.

Base-pair-resolution genome-wide mapping of active RNA polymerases using precision nuclear run-on (PRO-seq)

We provide a protocol for precision nuclear run-on sequencing (PRO-seq) and its variant, PRO-cap, which map the location of active RNA polymerases (PRO-seq) or transcription start sites (TSSs) (PRO-cap) genome-wide at high resolution. The density of RNA polymerases at a particular genomic locus directly reflects the level of nascent transcription at that region. Nuclei are isolated from cells and, under nuclear run-on conditions, transcriptionally engaged RNA polymerases incorporate one or, at most, a few biotin-labeled nucleotide triphosphates (biotin-NTPs) into the 3' end of nascent RNA. The biotin-labeled nascent RNA is used to prepare sequencing libraries, which are sequenced from the 3' end to provide high-resolution positional information for the RNA polymerases. PRO-seq provides much higher sensitivity than ChIP-seq, and it generates a much larger fraction of usable sequence reads than ChIP-seq or NET-seq (native elongating transcript sequencing). Similarly to NET-seq, PRO-seq maps the RNA polymerase at up to base-pair resolution with strand specificity, but unlike NET-seq it does not require immunoprecipitation. With the protocol provided here, PRO-seq (or PRO-cap) libraries for high-throughput sequencing can be generated in 4-5 working days. The method has been applied to human, mouse, Drosophila melanogaster and Caenorhabditis elegans cells and, with slight modifications, to yeast.

Transcription factors GAF and HSF act at distinct regulatory steps to modulate stress-induced gene activation

The coordinated regulation of gene expression at the transcriptional level is fundamental to development and homeostasis. Inducible systems are invaluable when studying transcription because the regulatory process can be triggered instantaneously, allowing the tracking of ordered mechanistic events. Here, we use precision run-on sequencing (PRO-seq) to examine the genome-wide heat shock (HS) response in Drosophila and the function of two key transcription factors on the immediate transcription activation or repression of all genes regulated by HS. We identify the primary HS response genes and the rate-limiting steps in the transcription cycle that GAGA-associated factor (GAF) and HS factor (HSF) regulate. We demonstrate that GAF acts upstream of promoter-proximally paused RNA polymerase II (Pol II) formation (likely at the step of chromatin opening) and that GAF-facilitated Pol II pausing is critical for HS activation. In contrast, HSF is dispensable for establishing or maintaining Pol II pausing but is critical for the release of paused Pol II into the gene body at a subset of highly activated genes. Additionally, HSF has no detectable role in the rapid HS repression of thousands of genes.

Defining the Essential Function of Yeast Hsf1 Reveals a Compact Transcriptional Program for Maintaining Eukaryotic Proteostasis

Despite its eponymous association with the heat shock response, yeast heat shock factor 1 (Hsf1) is essential even at low temperatures. Here we show that engineered nuclear export of Hsf1 results in cytotoxicity associated with massive protein aggregation. Genome-wide analysis revealed that Hsf1 nuclear export immediately decreased basal transcription and mRNA expression of 18 genes, which predominately encode chaperones. Strikingly, rescuing basal expression of Hsp70 and Hsp90 chaperones enabled robust cell growth in the complete absence of Hsf1. With the exception of chaperone gene induction, the vast majority of the heat shock response was Hsf1 independent. By comparative analysis of mammalian cell lines, we found that only heat shock-induced but not basal expression of chaperones is dependent on the mammalian Hsf1 homolog (HSF1). Our work reveals that yeast chaperone gene expression is an essential housekeeping mechanism and provides a roadmap for defining the function of HSF1 as a driver of oncogenesis.