Stress induces a rapid recruitment of RNA Pol II at stress-responsive loci
Analyses of gene expression have shown there is a rapid and strong induction of a set of stress-responsive genes in response to stress [2–9]. We quantified the increased fold induction of gene expression of 662 stress-responsive genes from microarray analysis (Materials and methods) and found an overall 6.4-fold increase of gene expression upon osmostress (0.4 M NaCl for 10 minutes). The induction pattern of these osmoresponsive genes in other stress conditions, such as heat shock (15 minutes at 37°C), oxidative stress (320 mM H2O2, 30 minutes), protein folding (250 mM dithiothreitol, 60 minutes) and amino acid starvation (30 minutes) [2] showed that osmoresponsive genes display a different expression pattern depending on each stress. In general, there is a poor overlap among the different stresses, with heat and osmostress overlapping the most (32%; Figure S1 in Additional file 1).
Whilst osmostress-induced genes showed a clear induction upon stress, the overall transcription of the whole genome, excluding the set of osmostress-induced genes, showed a 0.16-fold reduction in gene expression upon stress (Figure 1a). These data are consistent with previous reports [9] and indicate there must be a specialized mechanism that permits specific gene expression when global down-regulation of gene expression occurs.
To characterize how the changes of gene expression in response to osmostress are accomplished, we analyzed genome-wide binding of RNA Pol II in response to osmostress by ChIP-Seq in wild-type and hog1 cells. Association of RNA Pol II with ORFs is reduced when the overall genome is considered, whereas it clearly increases for stress-responsive genes (Figure 1b). Earlier studies showed that some housekeeping genes suffered a strong reduction of RNA Pol II occupancy at early time points in response to stress [22]. This is exemplified by the increase of RNA Pol II at the STL1 osmoresponsive gene in contrast to the reduction of overall RNA Pol II observed at the housekeeping gene PMA1, which encodes an essential H-ATPase (Figure 1c). Thus, STL1 and PMA1 genes are clear examples that represent the trend in osmoresponsive versus housekeeping genes. It also shows that while RNA Pol II is lost at the housekeeping genes in both wild-type and hog1 strains, the wild-type strain shows a faster recovery of RNA Pol II. Of note, the down-regulation of RNA Pol II in house-keeping genes precedes the recruitment of RNA Pol II at stress-responsive genes, indicating that the overall reduction of RNA Pol II occupancy cannot be due to a decrease in its availability. Taken together, genome-wide RNA Pol II localization suggests a strong bias for its localization towards stress-responsive genes.
Stress-responsive genes can be classified into two groups, Hog1-dependent and Hog1-independent, on the basis of gene expression data (Materials and methods; Additional file 2). When the 100 most responsive genes of each group were analyzed we found that indeed there was a clear difference in the degree of induction; the Hog1-dependent genes displayed a fold change that was almost four times higher compared to Hog1-independent genes (Figure S2 in Additional file 1). We then analyzed each group with regard to RNA Pol II association and found that RNA Pol II was recruited significantly to Hog1-dependent (green dots) and Hog1-independent (red dots) responsive genes in response to stress (Figure 1d, left-hand panels). By contrast, recruitment of RNA Pol II to Hog1-dependent genes was significantly different between wild-type and hog1 strains (green dots), while no differences were observed between a wild-type and a hog1 strain with regard to Hog1-independent genes (red dots) (Figure 1d, right-hand panels). The association of RNA Pol II with down-regulated genes or housekeeping genes in a wild-type strain was similar to that in a hog1 strain, indicating that Hog1 does not play a role in the initial changes observed in non-stress-dependent genes (Figure S3 in Additional file 1). Thus, Hog1 has a crucial role in the genome-wide redistribution of RNA Pol II to stress-responsive genes upon stress.
Hog1 associates with the chromatin of RNA Pol II and Pol III genes
Genome-wide studies using ChIP and microarray analysis have been instrumental in uncovering the presence of Hog1 associated with a number of stress-responsive genes as well as its localization at both promoter and coding regions of stress-responsive genes. However, the number of genes uncovered by these approaches has been rather limited and never totaled more than 70 genes [15–17]. The relevance of Hog1 in gene expression and RNA Pol II recruitment suggested that the number of genes with Hog1 association could have been underestimated. We undertook ChIP-Seq analysis to improve the sensitivity of detection and found that Hog1 is present in at least 340 genome loci upon osmostress (0.4 M NaCl for 5 minutes; Figure S4a in Additional file 1). We analyzed binding at 5 minutes because this was the peak of Hog1 association with STL1 and CTT1 [19] (Figure S5 in Additional file 1). Albeit ChIP experiments generate data that are population averages, previous single cell analyses showed that Hog1 is activated in all cells similarly upon osmostress and that transcriptional induction correlates very well with the localization of Hog1 to stress-responsive genes [23]. Recruitment of Hog1 was not restricted to RNA Pol II transcribed genes but was present, albeit to a lesser extent, at RNA Pol III transcribed genes as well as long terminal repeat (LTR) DNA regions.
When the presence of Hog1 was analyzed on RNA Pol II transcribed genes, we found that Hog1 was associated with approximately 80% of genes, with expression described to be highly dependent on the SAPK, confirming that Hog1 is widely associated with Hog1-regulated genes. By contrast, only 30% of the genes induced upon osmostress are Hog1-independent and showed Hog1 associated with their loci (Figure 2a; Figure S4b in Additional file 1, green and red dots). Hog1 was not present at down-regulated or house-keeping genes (Figure S4b in Additional file 1, blue and yellow dots); therefore, Hog1 is associated with stress-responsive RNA Pol II genes.
Hog1 has been shown to associate with promoters and the ORFs of stress-responsive genes. We asked whether there was a biased interaction towards promoters or ORFs and found that Hog1 was associated with both promoters and ORFs in a significant number of genes (more than 41%). Also, for those genes where only one region was above the threshold, association of Hog1 with ORFs was more prominent than with promoters (Figure S6 in Additional file 1). When we analyzed association of Hog1 with Hog1-dependent or Hog1-independent genes, we found that Hog1 binding is biased slightly towards the ORFs in Hog1-dependent genes, whereas for genes for which Hog1 was less relevant, Hog1 localization was associated mainly with promoters (Figure 2b). Several scenarios could explain the presence of the SAPK at Hog1-independent loci, such as the use of a too stringent threshold for Hog1 dependency, their presence close to a Hog1-dependent gene, or the induction of the gene is mediated by redundant pathways, including Hog1. Thus, stress-induced RNA Pol II transcribed genes appear to have strong enrichment of Hog1 at their promoters and ORFs.
Remarkably, Hog1 was also associated with RNA Pol III transcribed genes, including at least 16 tRNA genes as well as the two reference genes SCR1 and RPR1 (Figure 2c). ChIP experiments showed similar kinetics of association of Hog1 with tF(GAA)D, RPR1 and SCR1 as with RNA Pol II transcribed genes. We then investigated the association of RNA Pol III (Rpc82 subunit) with two tRNA loci (tF(GAA)D and tP(UGG)O3) and found that, albeit RNA Pol III disassociated rapidly from chromatin in a stress-dependent manner, a rapid recovery of RNA Pol III levels occurred in wild-type that was not observed in hog1 cells (Figure S7 in Additional file 1).
It has been reported that Hog1 interacts with RNA Pol II (most likely through Rpb1), which facilitates gene expression in RNA Pol II transcribed genes [13, 17]. Therefore, we used co-precipitation experiments in extracts from cells expressing endogenously tagged HA-Hog1 and Myc-Rpc82 (a subunit of the RNA Pol III complex not shared with RNA Pol II) to assess whether Hog1 is able to interact with RNA Pol III. We found Hog1 was able to interact with endogenous tagged-Rpc82 and vice versa (Figure 2d). It is worth noting that this interaction was observed only when cells were subjected to osmostress. Thus, Hog1 is targeted to RNA Pol III loci in response to stress and associates physically with RNA Pol III, as it does with RNA Pol II transcribed genes.
Efficient recruitment of RNA Pol II and maximal gene expression requires Hog1
To assess the relevance of the association of Hog1 with RNA Pol II, we compared the degree of gene expression of stress-responsive genes with the presence of RNA Pol II and Hog1 (Figure 3a). Several groups of osmoresponsive genes can be identified, depending on the presence of Hog1 and/or RNA Pol II upon stress. A group of genes showed no significant association of Hog1 with RNA Pol II but were induced upon osmostress. This group of genes correlated quite well with genes that showed stabilization of mRNAs upon stress (40 out of the 43 analyzed were found to be stabilized) [8, 9]. We found a significant number of genes with increased RNA Pol II association that did not have Hog1 present on them (91 out of 391); these genes correspond to Hog1-independent genes. There was, however, a prominent overlapping group of genes that showed increased expression, increased recruitment of RNA Pol II and association with Hog1 (a total of 144 genes; Figure 3a). We then compared the degree of gene induction in those groups of genes and found that there was a strong correlation between the presence of both Hog1 and RNA Pol II with high expression rates when compared to genes that were enriched only with RNA Pol II (Figure 3b). If the presence of Hog1 improved the recruitment of RNA Pol II and transcription, it should be possible to establish a quantitative relationship between weak, or strong, Hog1 binding with RNA Pol II and transcription. Co-localization studies of Hog1 with RNA Pol II showed that RNA Pol II association with stress-responsive genes was more efficient for genes with higher Hog1 association (Figure 3c). Therefore, a high level of induction in stress-responsive genes is accomplished by strong association with Hog1 and increased RNA Pol II recruitment.
It has been reported that very low salt (0.1 M NaCl) stress results in maximal Hog1 activation. When cells are exposed to higher concentrations of NaCl, however, activation of Hog1 remains associated with stress-responsive loci for an extended period of time [19, 23]. If the presence of Hog1 improves the recruitment of RNA Pol II and transcription efficiency, it would be expected that an increase of Hog1 at specific promoters will result in enhanced recruitment of RNA Pol II and transcription. We followed STL1 expression in response to 0.1, 0.2 and 0.4 M NaCl together with the association of Hog1 and RNA Pol II. Weak expression of STL1 was observed at 0.1 and 0.2 M NaCl, in clear contrast to the induction observed with 0.4 M NaCl (similar results were observed for CTT1 and ALD3). Remarkably, the initial recruitment of RNA Pol II at STL1 was similar at 0.2 and 0.4 M NaCl; however, the residence time of Hog1 at the loci was clearly shorter (Figure 3d; Figure S8 in Additional file 1). Thus, the association time of Hog1 with stress loci appears to be crucial for determining the degree of gene induction upon osmostress.
Hog1 mediates chromatin changes at stress-responsive loci
Hog1 stimulates chromatin remodeling at specific stress-responsive loci by recruiting the RSC chromatin remodeler [21]. We investigated whether all stress-responsive genes were subjected to changes on chromatin organization and the relevance of the SAPK to those changes. We used genome-wide MNase digestion of chromatin and deep sequencing (MNase-Seq) before and after stress. Wild-type and hog1 strains were subjected (or not) to osmostress and cells were fixed before digestion of chromatin by MNase to prevent Hog1 activation during the preparation of spheroplasts (see Materials and methods). The nucleosomal profile around the transcription start site (TSS) in genes that were not regulated upon stress did not change when cells were subjected to osmostress (Figure 4a, upper panel). We then analyzed nucleosome positioning in stress-induced genes with expression that does not depend on Hog1. Upon osmostress there were slight changes of nucleosome occupancy, especially around the TSS and those changes were similar in hog1-deficient cells (Figure 4a, middle panel). In clear contrast, when Hog1-dependent genes were analyzed, a dramatic change of nucleosome occupancy occurred upon stress at both the promoter and ORF regions. These changes on nucleosomes were completely abolished in hog1 cells. It is noteworthy that the +1 nucleosome in stress-responsive genes appears to be shifted slightly when compared to localization of the genome-wide +1 nucleosome, suggesting a particular chromatin structure for stress-responsive genes. Color maps of -1,000 to +1,000 of each ORF aligned by TSS have also been included to illustrate nucleosome organization genome-wide and in stress-responsive genes (Figure S9 in Additional file 1). Taken together, efficient nucleosome re-organization at stress-responsive genes is completely dependent on the SAPK (Figure 4a, lower panel). When changes on the chromatin structure were quantified (percentage of nucleosome occupancy), we found that osmostress induces a 25% decrease in nucleosome occu-pancy in Hog1-independent genes (similar to that found in wild-type and hog1 cells), whereas it was decreased 51% in Hog1-dependent genes (Figure 4b). Thus, Hog1 is crucial to inducing major changes in chromatin structure.
In summary, genome-wide binding studies in combination with analysis of chromatin structure have shown that Hog1 serves to bypass the general down-regulation of gene expression that occurs in response to stress. Hog1 permits efficient targeting of the RNA Pol II machinery and, furthermore, induces major changes of chromatin structure at stress-responsive loci. The combination of targeted recruitment of RNA Pol II with chromatin remodeling is essential to maximize gene expression in response to external stimuli.