- Open Access
The histone deacetylase Rpd3p is required for transient changes in genomic expression in response to stress
Genome Biology volume 10, Article number: R57 (2009)
Yeast responding to stress activate a large gene expression program called the Environmental Stress Response that consists of approximately 600 repressed genes and approximately 300 induced genes. Numerous factors are implicated in regulating subsets of Environmental Stress Response genes; however, a complete picture of Environmental Stress Response regulation remains unclear. We investigated the role of the histone deacetylase Rpd3p, previously linked to the upstream regions of many Environmental Stress Response genes, in producing Environmental Stress Response gene expression changes in response to stress.
We found that the Rpd3-Large complex is required for proper expression of both induced and repressed Environmental Stress Response genes under multiple stress conditions. Cells lacking RPD3 or the Rpd3-Large subunit PHO23 had a major defect in Environmental Stress Response initiation, particularly during the transient phase of expression immediately after stress exposure. Chromatin-immunoprecipitation showed a direct role for Rpd3-Large at representative genes; however, there were different effects on nucleosome occupancy and histone deacetylation at different promoters. Computational analysis implicated regulators that may act with Rpd3p at Environmental Stress Response genes. We provide genetic and biochemical evidence that Rpd3p is required for binding and action of the stress-activated transcription factor Msn2p, although the contribution of these factors differs for different genes.
Our results implicate Rpd3p as an important co-factor in the Environmental Stress Response regulatory network, and suggest the importance of histone modification in producing transient changes in gene expression triggered by stress.
Sudden environmental changes can trigger rapid and dramatic changes in genomic expression. This involves coordinated expression of hundreds to thousands of genes, whose expression is precisely modulated in timing and magnitude. Many different transcription factors function in the cell at any given time and respond to distinct upstream signals. Therefore, cells must integrate the action of numerous signals and regulatory factors to produce a coherent genomic expression program customized for each new environment.
Yeast respond to stress in part by initiating the Environmental Stress Response (ESR), which consists of approximately 600 genes whose expression is repressed and approximately 300 genes whose expression is induced by diverse stresses [1, 2]. The repressed genes include approximately 130 ribosomal protein ('RP') genes and a distinct group of approximately 450 genes more broadly related to protein synthesis ('PS genes'). Both groups are highly expressed in actively growing cells but sharply repressed, with slightly different expression profiles, in response to stress. Genes induced in the ESR ('iESR genes') are involved in varied aspects of stress defense, including redox regulation, protein folding, osmo-tolerance, cell signaling, and other functions. Initiation of the ESR is not required to survive the offending stress but rather helps to protect cells against subsequent severe doses of the same or different stress (although it cannot fully explain acquired stress resistance in all cases) .
Although activated by many different stresses, the ESR is regulated differently depending on the environment. Numerous upstream signaling pathways have been implicated in condition-specific ESR regulation, including the high osmolarity glycerol (HOG)  (Jessica Clarke and APG, unpublished data), MEC , and protein kinase C (Scott Topper and APG, unpublished) pathways in response to osmotic shock, DNA damage, or reducing agents, respectively, and the protein kinase A and target of rapamycin (TOR) pathways upon stress relief [6–10] (reviewed in ). Furthermore, different subsets of iESR genes can be induced by stress-specific transcription factors, such as the oxidative-stress factor Yap1p , the heat shock factor Hsf1p [12–14], Sko1p and Hot1p upon osmotic stress [15–18], and the 'general-stress' transcription factors Msn2p and Msn4p in response to diverse stresses (reviewed in ). However, little is known about how these signals are integrated to mediate ESR initiation, or how genes repressed in the ESR are coordinated with genes induced in the program.
One mechanism of altering gene expression is through changes in chromatin state. The histone deacetylase Rpd3p deacetylates histones in both coding and noncoding regions, where it is thought to function in at least two distinct complexes (reviewed in [19, 20]). A small complex (Rpd3S) suppresses cryptic transcription initiation by deacetylating histones after elongating polymerase [21–23]. Rpd3S is recruited via the combined action of the Eaf3p and Rco1p subunits to histone H3 methylated by Set2p during transcription of the open reading frame [21–23]. In contrast, a large complex (Rpd3L) is recruited to promoters by site-specific DNA binding proteins, including the Ume6p subunit of Rpd3L, where it is thought to function in transcription initiation [23–27]. Rpd3p is known to bind different promoters under different conditions, such as cold shock and rapamycin treatment [28–30]. In fact, many promoters to which Rpd3p relocalizes are of genes repressed in the ESR. The effects of Rpd3p at these promoters have not been shown on a global scale, but the result suggests Rpd3p is required for stress-dependent repression of ESR genes [11, 30].
Although traditionally linked to repression, histone deacetylases can also function during gene activation [31–36]. Induction of several different yeast genes requires Rpd3p following salt treatment, hypoxia, or DNA damage [32–34]. The precise mechanism is not clear but requires Rpd3p for recruitment of RNA polymerase to promoters of genes (including iESR genes) induced by osmotic shock and DNA damage [32, 34]. Furthermore, induction of hypoxic genes requires Rpd3p-dependent histone deacetylation for nucleosome displacement and stable binding of the Upc2p transcription factor within the genes' regulatory regions . That Rpd3p has been linked to stress-dependent gene induction and repression raised the possibility that Rpd3p participates in regulating both induced and repressed genes within the ESR.
Indeed, here we show that Rpd3p is required for proper initiation of the ESR, including normal regulation of both induced and repressed genes, in yeast responding to multiple stresses. Cells lacking RPD3 or the Rpd3L subunit PHO23 had a major defect, specifically during the transient phase immediately after H2O2 treatment, while cells lacking the Rpd3S subunit RCO1 did not. Chromatin-immunoprecipitation (ChIP) at candidate ESR genes revealed that Rpd3p moves to numerous promoters upon stress to mediate histone deacetylation; however, the precise pattern of chromatin change was different for different nucleosomes and genes investigated. We show that Rpd3p binds directly to genes induced by stress and is required for normal binding of Msn2p to numerous promoters. Together, this work implicates Rpd3L as an important co-factor in the ESR regulatory network.
Rpd3p is required for the full dynamic range of stress-activated gene expression changes
We followed genomic expression in wild-type and rpd3Δ cells responding over time to a 25°C to 37°C heat shock, 0.4 mM H2O2, and 0.75 M NaCl. A large fraction (56 to 80%) of the gene expression changes seen in wild-type cells was affected by RPD3 deletion, and this included both repressed and induced genes (Table 1). In particular, Rpd3p was required for normal expression of the vast majority of ESR genes (Figure 1). Repression of PS genes was heavily dependent on Rpd3p in response to all stresses, whereas repression of RP genes required Rpd3p for full repression in response to heat and H2O2 stress but not salt treatment. Normal induction of iESR genes also required Rpd3p, since the rpd3Δ strain displayed more than twofold decreased induction levels at the peak of the response. Interestingly, a subset of iESR genes (approximately 50% at a false discovery rate of 0.05) showed slight derepression (approximately 1.5-fold) in the rpd3Δ mutant in the absence of stress (Figure 1; Figure S1 in Additional data file 1). The defect in stress-dependent induction was not due to an already activated stress response in mutant cells, indicated by normal cytosolic localization of Msn2p before stress but substantial Msn2p nuclear accumulation after stress, similar to wild-type (Figure S2 in Additional data file 1). Furthermore, these iESR genes (as well as those with no significant difference in basal expression) still had a defect in induction beyond what could be accounted for by basal expression differences (Figure S1 in Additional data file 1). Thus, Rpd3p is required for the induction and repression of ESR genes during stress, although each ESR subgroup shows a qualitatively different dependence on the protein.
Stress-dependent gene expression changes are often transient, in that large changes immediately after stress subsequently relax to new 'steady-state' levels as cells acclimate (reviewed in ). We found that Rpd3p is particularly important for this transient phase of expression (Figure 1b). PS genes showed almost no transient expression changes, while iESR genes showed reduced expression levels specifically at the peak of the transient phase. RP genes also showed diminished expression differences at the peak of the response to heat shock and H2O2 treatment. Despite the defect in transient ESR expression, the rpd3Δ mutant eventually reached near-wild-type expression changes by the end of these time courses. This indicates that Rpd3p is not necessarily required to maintain new steady-state levels of expression in cells acclimated to high temperature or H2O2, but is critical in producing a large, rapid response to stress.
ESR regulation requires histone deacetylase activity through the Rpd3L complex
We found that the catalytic activity of Rpd3p, as well as modifiable histones and subunits of the Rpd3L complex, were required for proper ESR regulation. Cells harboring the catalytically inactive rpd3-H150:151A protein  or treated with the Rpd3p inhibitor trichostatin A displayed the same widespread defect as the rpd3Δ strain (Figure S3 in Additional data file 1). A similar defect was observed in cells harboring a mutant histone H4 (H4KQ), in which amino-terminal lysines were changed to glutamine to mimic the acetylated histone state  (Figure S3 in Additional data file 1). This effect was particularly clear for PS and iESR genes, although there was only a subtle defect in repression of the RP genes in the H4 mutant strain.
To distinguish between the effects of the different Rpd3p complexes, we characterized the H2O2 response in cells lacking Pho23p or Rco1p, exclusive members of the Rpd3L and Rpd3S complexes, respectively [21, 23, 39]. The expression defect seen in the pho23Δ mutant, but not the rco1Δ cells, was highly similar to that in the rpd3Δ mutant. Over 80% of Rpd3p-affected genes were equally dependent on Pho23p (R = 0.94, m = 0.98), whereas less than 12% of Rpd3-affected genes showed a partial expression defect in cells lacking RCO1. Furthermore, the pho23Δ strain showed the same defect in transient expression as the rpd3Δ cells (Figure S4 in Additional data file 1). In contrast, the rco1Δ cells showed large changes in expression similar to wild type, albeit with a slightly delayed response that is difficult to interpret due to spurious internal transcripts in this mutant [21, 23]. Nonetheless, these data show that defects in the magnitude and transience of gene expression can be accounted for by the Rpd3L complex. Consistent with previous studies [28, 40, 41], we found few of the Rpd3L-dependent expression changes were dependent on the Ume6p subunit (data not shown), which is thought to recruit the complex to specific loci [24, 25, 27]. This suggests that other DNA binding proteins may be required for Rpd3L-dependent gene expression changes (see below).
Representative ESR genes show Rpd3p-dependent changes in chromatin following stress
Previous studies showed Rpd3p physically bound to many ESR-gene promoters during times of stress [28–30]. Global studies probing Rpd3p binding after cold shock (inadvertently inflicted by ) and rapamycin treatment  showed that promoters of 60% of PS genes (P < 10-32) and 90% of RP genes (P < 10-20) were bound by Rpd3p. Few of these regions are bound under standard conditions [29, 30]. Roughly 20% of iESR-gene upstream regions were bound by Rpd3p under stress conditions, though this may be an underestimate, since chromatin-remodeling enzymes are difficult to ChIP, particularly during dynamic responses . Consistent with these studies, we found Rpd3p bound upstream of four representative ESR genes (including one PS, one RP, and two iESR genes) after H2O2 treatment (Figure 2). Three of the targets also showed some Rpd3p binding before stress, and all but the UBC5 promoter showed increased Rpd3p binding after H2O2 treatment. These results were similar to those seen in cold-shock (Figure 2), suggesting that many of the previously observed binding events from  also occur during H2O2 stress.
We therefore characterized changes in nucleosome occupancy and H4 acetylation at nucleosomes spanning the same four ESR genes in wild-type, rpd3Δ or pho23Δ strains using mononucleosome digestion and ChIP of acetylated H4 before and after H2O2 exposure. The results showed different trends at different genes. Nucleosomes at repressed ESR genes GAR1 and RPL16A showed Rpd3L-dependent changes in histone deacetylation following H2O2 treatment. Though wild-type cells showed an approximately three- to eightfold decrease (depending on the gene and nucleosome) in the fraction of acetylated nucleosomes (Figure 3b), both the rpd3Δ and pho23Δ mutants had a major defect in histone deacetylation across both repressed ESR genes. This defect correlated with the defect in their H2O2-dependent repression (Figure 3c). Interestingly, the rpd3Δ mutant, and to some extent the pho23Δ strain, also had a defect in nucleosome repositioning at these repressed genes: whereas wild-type cells responding to H2O2 showed a dramatic increase in nucleosome occupancy upstream of RPL16A, the rpd3Δ mutant showed a major defect in this response (Figure 3a). The pho23Δ mutant displayed a weaker defect than the rpd3Δ strain, indicating that Pho23p is only partially required for the stress-dependent increase in nucleosome occupancy at this locus. Together with results in Figure 2, this indicates that Rpd3L-dependent histone deacetylation is required for repression of these PS and RP genes.
The two representative iESR genes each displayed a unique profile in chromatin change. Nucleosomes surrounding the transcription start site of the induced gene UBC5 displayed decreased histone acetylation in wild-type cells but not the rpd3Δ or pho23Δ mutants responding to H2O2. In addition, nucleosome occupancy at these loci increased in wild-type cells, but not the mutants. In contrast, both the promoter and open reading frame of iESR gene XKS1 showed increased histone acetylation and nucleosome loss in wild-type cells, with no significant defect in either mutant. Nonetheless, this gene showed approximately threefold weaker induction in the rpd3Δ and pho23Δ mutants, specifically during the transient phase of expression. This reveals a decoupling of chromatin changes upstream of XKS1 and XKS1 gene induction in the mutant strains responding to stress, in a manner dependent on direct Rpd3p binding to the region (see Discussion).
Implication of Rpd3p-dependent and -independent transcriptional regulators
The above results indicate that Rpd3p has different effects at different ESR genes, perhaps due to different regulators functioning at those genes. To identify additional stress-dependent regulators, we systematically analyzed clustered expression data for enrichment of known transcription factor targets or functional gene groups. We manually identified gene clusters in the hierarchically clustered dataset and scored enrichment of Gene Ontology annotations , targets of known transcription factors , and genes with different upstream cis-regulatory elements . This analysis pointed to transcription factors involved in the Rpd3p-dependent and Rpd3p-independent regulation of gene expression (Table S1 in Additional data file 2).
Multiple clusters of Rpd3p-dependent induced genes were enriched for genes with upstream Msn2p and Msn4p binding sites (CCCCT [45, 46]), consistent with the known role of Msn2/4p in regulating iESR genes [1, 46, 47]. Another cluster of Rpd3-dependent repressed genes was heavily enriched for genes with upstream Polymerase A and C (PAC; GCGATGAG) elements and Ribosomal RNA Processing Elements (RRPEs; AAAAWTTTT), known to be enriched in PS genes and previously linked to promoters bound by Rpd3p [1, 28, 41, 48]. Another cluster was enriched for proteasome genes and genes containing binding sites of the proteasome regulator Rpn4p. These associations raise the possibility that Rpd3p may work with these factors to mediate the observed gene expression changes (see more below).
Interestingly, we identified some genes whose expression was conditionally dependent on Rpd3p. Targets of the heat shock transcription factor Hsf1p or the oxidative stress transcription factor Yap1p were only dependent on Rpd3p in response to specific conditions (Figure 4). The majority of Hsf1p targets did not require Rpd3p for induction following heat shock but showed Rpd3-dependent induction in response to H2O2 and NaCl treatment (Figure 4a). Similarly, induction of Yap1p targets (Figure 4b) was independent of Rpd3p in response to H2O2, while a subset induced with the ESR required Rpd3p for full induction following heat shock and salt stress only. Hsf1p and Yap1p are known to be condition-specific regulators of subsets of iESR genes, functioning during heat shock and oxidative stress, respectively (reviewed in ). Under other conditions, many of these genes are regulated by Msn2/4p. Our observations are consistent with the model that Hsf1p and Yap1p function independently of Rpd3p to regulate gene induction, whereas Msn2/4p act in an Rpd3p-dependent manner.
Rpd3p is required for normal Msn2p binding and transcription initiation
To investigate the link between Msn2/4p and Rpd3p function, we measured genomic expression in strains lacking RPD3, MSN2/MSN4, or MSN2/MSN4/RPD3 as cells responded to H2O2. Interestingly, genes fell into different categories depending on their expression defect (Figure 5). One class of genes was equally dependent on Rpd3p and Msn2/4p for induction, with no additional defect in the triple mutant (Figure 5a). A second class required both sets of factors but was more dependent on Msn2/4p (Figure 5b), while a third class suggests redundant function of Rpd3p and Msn2/4p at these genes (Figure 5c). The latter group was enriched for genes involved in carbohydrate metabolism (P < 10-8) and trehalose synthesis (P < 10-5), suggesting functional relevance of the categorization. A fourth class of genes was dependent only on Rpd3p (data not shown), indicating that additional Rpd3p-dependent transcription factors are required for proper initiation of the ESR (including Rpn4p and others). Importantly, a fifth group of genes was dependent only on Msn2/Msn4p (data not shown), which underscores that the Rpd3p-dependent defect in iESR-gene induction is not simply caused by failure to activate Msn2/4p, consistent with microscopy data (Figure S2 in Additional data file 1). Thus, most but not all of Msn2/4p-dependent genes require Rpd3p for full induction, and these targets show qualitative differences in their dependence.
These results suggest Rpd3p may be required for Msn2/4p action during gene induction. We therefore measured Msn2p binding upstream of genes representing each category above in wild-type and rpd3Δ cells responding to H2O2 (Figure 5d). While none of the promoters tested showed Msn2p bound before stress, as expected, wild-type cells showed an increase in Msn2p promoter binding that was defective in the rpd3Δ strain at most targets, regardless of class. The exception was YPR127W, an Msn2p-dependent but Rpd3p-independent target, which showed no significant defect in Msn2p binding in the rpd3Δ strain. Thus, Rpd3p was required for Msn2p binding upstream of targets that showed dependence on Rpd3p for induction.
It is important to note that over half the H2O2-induced gene expression changes were not affected by RPD3 deletion or MSN2/4 deletion. This underscores that Rpd3p is not universally required for all gene expression changes in response to stress, and shows that the defect in expression is not due to a gross alteration in the rpd3Δ mutant's response.
Rpd3p is required for ESR suppression following stress relief
That Rpd3p is implicated in both gene induction and repression following stress treatment raised the possibility that Rpd3p participates in the reciprocal regulation of the same genes during stress relief, when the ESR is suppressed. To test the role of Rpd3p in ESR suppression, we measured gene expression in wild-type and rpd3Δ cells acclimated to 37°C as cells were returned to 25°C. Strikingly, the rpd3Δ strain had a significant defect in ESR suppression during stress relief (Figure S5 in Additional data file 1): whereas wild-type cells rapidly repressed expression of iESR genes in response to stress relief, rpd3Δ cells displayed a significantly weaker response. Similarly, induction levels of PS genes were significantly smaller in the rpd3Δ strain compared to the wild-type cells recovering from stress. Consistent with results presented above, the RP genes were distinct in that induction upon stress relief was only mildly affected by RPD3 deletion. These results suggest that Rpd3p is not exclusively required for the repression or for the induction of the ESR genes but instead is required for proper changes in the genes' expression regardless of the directionality of the change.
Our results reveal that Rpd3p is required for many stress-dependent gene expression changes, particularly genes in the yeast ESR. We show that Rpd3p and the Rpd3L subunit Pho23p (but not the Rpd3S component Rco1p), as well as Rpd3p catalytic activity and modifiable histones, are required to produce these effects. Rpd3p binds directly to promoters of representative ESR genes, indicating that the Rpd3-dependent changes in chromatin structure that we see are direct at these promoters. Furthermore, the observed defects in iESR induction correlate with decreased Msn2p binding at candidate promoters in the rpd3Δ strain. Together with previous global studies of Rpd3p localization [28–30], these results indicate that Rpd3L acts directly at many ESR genes to mediate transient changes in gene expression. The defect in stress-activated expression leads to a corresponding defect in acquired stress resistance (Figure S6 in Additional data file 1), similar to that we have previously shown in cells lacking Msn2p and/or Msn4p . Thus, Rpd3p is an important cofactor in initiating the ESR. Models for how Rpd3p fits into the ESR regulatory network are discussed below.
Role of Rpd3p in ESR initiation under diverse stress conditions
Rpd3p likely acts with distinct transcription factors at different classes of ESR promoters. PS genes are heavily enriched for upstream PAC elements (GCGATGAG) and RRPEs (AAAAWTTTT) [1, 48], which have also been linked to Rpd3p binding . Recently, the binding proteins of both elements have been identified and linked to PS expression. PAC is bound by Dot6p and Pbs1p [49, 50], and deletion of the two genes leads to defective PS gene repression in response to heat shock . The RRPE binding factor was recently identified as Stb3p, which interacts with the Sin3p subunit of Rpd3p complexes [51, 52] and is required for PS gene induction upon starvation relief but represses PS gene transcription when overexpressed (D Liko and W Heideman, personal communication). Although we found no expression defect in an stb3Δ mutant responding to stress (data not shown), the link between Stb3p, Sin3p/Rpd3p, and RRPEs suggests that the proteins function together at this regulatory motif to affect PS gene expression.
Rpd3p has a distinct role in repressing RP genes, since their expression was mildly Rpd3L-dependent under certain conditions only. Nonetheless, we found that Rpd3p moves to the promoter of RPL16A upon H2O2 treatment (Figure 2), as previously found in response to cold shock [28, 30], and is required for normal histone deacetylation and nucleosome deposition/repositioning (Figure 3). Rpd3p has previously been linked to RP gene repression after rapamycin treatment [29, 53, 54], although we found no requirement for the proposed repressor Crf1p (data not shown) . We have, however, found a requirement for the ATP-dependent nucleosome-remodeling complex, RSC, which is important for proper nucleosome organization upstream of many genes . RSC mutants have increased RP expression in the absence of stress , while cells lacking Rsc1p fail to fully repress RP expression and, to some extent, PS gene expression upon H2O2 treatment (our unpublished data). Like Rpd3p, RSC binds RP promoters in a condition-specific manner . Thus, Rpd3p and RSC may function in parallel pathways at these genes. Interestingly, stress-dependent changes in nucleosome occupancy at RPL16A were only partially dependent on Pho23p, raising the possibility that Rpd3L functions partially independently of Pho23p or that Rpd3p is acting through multiple complexes, at least one of which does not require Pho23p [20–23].
The role of Rpd3L at iESR genes is less clear; however, our ChIP experiments suggest four general models for how Rpd3p may affect gene induction. The first is that some iESR genes may be indirectly affected by Rpd3L activity, particularly those for which there is no evidence of Rpd3p binding in response to stress. The second model is that Rpd3p plays an important and direct role in repressing iESR expression in the absence of stress, since Rpd3p binds directly to the promoters of UBC5 and XKS1 before stress (Figure 2) and these genes (plus nearly half of iESR genes) show slight derepression under normal conditions (Figure 1). This model is not incompatible with separate roles for Rpd3p in regulating stress-dependent expression changes, demonstrated by UBC5 and XKS1. At the UBC5 promoter, Rpd3p directly deacetylates promoter-based histones to mediate gene induction. This is consistent with results of De Nadal et al. , who showed Rpd3-dependent histone deacetylation is required for polymerase recruitment. In contrast, H2O2-dependent chromatin changes at XKS1 were not detectibly dependent on Rpd3L, despite increased Rpd3p binding upon treatment. The rpd3Δ mutant ultimately induced XKS1 to levels higher than wild type, but with a major defect in the normal transient burst of expression. Thus, the changes in histone acetylation did not lead to normal gene induction. One possibility is that gene induction triggered by H2O2 requires proper Rpd3-dependent promoter architecture before stress; alternatively, Rpd3p may play a role late in gene induction, after active nucleosome acetylation, as previously proposed for DNA damage-responsive genes .
We also show that Rpd3p activity is required for normal Msn2p binding to representative promoters. This is reminiscent of the requirement of Rpd3p for nucleosome displacement and Upc2p binding at the promoters of hypoxia-regulated genes . The exact mechanism of Rpd3p involvement at Msn2/4p targets is unclear; however, Lindstrom et al.  recently showed that Msn2/4p activity is inhibited by NuA4-dependent histone acetylation. This raises the possibility that histone deacetylation by Rpd3p counteracts the inhibitory effects of NuA4-dependent acetylation to allow Msn2p binding and gene induction. That different targets of Msn2/4p and Rpd3p show distinct sensitivities to the factors' deletion again implies distinct regulatory mechanisms for the different subclasses of targets. Understanding the differences in regulation will be an interesting area of future investigation.
Rpd3p functions as a 'general-stress' co-factor in the ESR regulatory network
The ESR regulatory network consists of condition-specific regulators - those that only regulate ESR expression under specific circumstances - as well as 'general-stress' factors (such as Msn2/4p) that function under a wide variety of conditions. Our results suggest Rpd3p acts with the 'general-stress' set of ESR regulators at iESR and PS genes. Rpd3L is required for proper expression of these genes in response to numerous stresses (Figure 1). Furthermore, Msn2/4-dependent induction, but not condition-specific regulation by Hsf1p and Yap1p, requires Rpd3p (Figures 4 and 5). Like Msn2/4p, the 'general stress' role of Rpd3p persists despite the involvement of different upstream regulators under different conditions. For example, De Nadal et al.  showed that Rpd3p is recruited to numerous iESR promoters in a manner dependent on the Hog1p kinase following salt stress but independent of Hog1p after heat shock. Thus, the involvement of Rpd3p, and the transcription factors it interacts with at these promoters, is controlled by different upstream signaling pathways under different environments. It will be interesting to decipher the mechanisms by which Rpd3p associates with stress-activated transcription factors despite distinct, condition-specific upstream pathways.
Rpd3p is required for the transient phase of stress-activated gene expression changes
This study also demonstrates the importance of histone modification in mediating rapid and transient responses to environmental changes. The Rpd3L complex is particularly important in producing the large, rapid expression changes during the period of stress acclimation. The transient expression changes produced by acute stress treatment are qualitatively distinct from continuous expression changes seen under different nutrients. However, Rpd3p can affect the rapid kinetics of both types of expression responses. Upon phosphate limitation, cells lacking RPD3 showed delayed induction of PHO5 but eventually altered expression similar to wild-type cells . Interestingly, a similar effect was reported in cells lacking the histone acetyltransferase Gcn5p, which also showed delayed induction of metabolic genes . These results reflect that changes in chromatin states, mediated by both deacetylases and acetyltransferases, are particularly important for rapid kinetics of gene-expression changes in response to variable environments. Consistently, we found that rpd3Δ cells display defects in reciprocal expression changes of the same genes upon stress exposure as well as stress relief. Dynamic and successive alterations in histone modification are crucial in producing proper transcriptional changes (for example, [61–67]). Elucidating the dynamics of chromatin changes upon stress treatment will continue to shed light on the dynamics of stress-dependent gene expression changes.
Rpd3p is an important co-factor in the regulatory network that controls ESR gene expression in response to stress, working with different factors at different subsets of ESR genes. Many questions remain about the mechanistic details of Rpd3p action at these promoters. While future studies will be required to dissect the precise mechanism of Rpd3p in regulating these genes, this work contributes to our understanding of the ESR regulatory network and provides an avenue for identifying additional factors that work with Rpd3p in regulating the ESR.
Materials and methods
Strains and growth conditions
Strains used in this study are listed in Table S2 in Additional data file 3. PHO23 and RCO1 deletion strains were purchased from Open Biosystems (Huntsville, AL, USA), and each deletion was verified by PCR. The rpd3Δ and msn2Δ msn4Δ rpd3Δ strains were constructed by homologous recombination to replace RPD3 with KANMX or LEU2 in BY4741 or AGY0249, respectively. Unless otherwise noted, cells were grown at 30°C in YPD medium. Although the growth rate of the rpd3Δ strain is approximately 1.5-fold slower than wild type, this cannot explain the observed expression defects, since the mutant phenotypes are recapitulated by the pho23Δ mutant, whose doubling rate is indistinguishable from wild type.
Cell collection for microarray analysis
Cells were grown approximately three doublings to an optical density (OD600) of approximately 0.6 to 0.8 and a sample was collected for the unstressed control, as previously described . Basal expression in rpd3Δ versus wild type was measured in triplicate. For heat shock time courses, cells were grown at 25°C, filtered and resuspended in 37°C YPD. Aliquots were collected at 5, 15, 30, 45, and 60 minutes (time course HS_1) or at 5, 10, 20, 30, and 60 minutes (time course HS_2) as previously described . For the H2O2 experiments, peroxide was added to 0.4 mM and cells were collected at 10, 20, 30, 40, and 60 minutes (time course H2O2_1) or at 30 minutes for single-time point experiments, done in triplicate. For sodium chloride (NaCl) time courses, NaCl was added to 0.75 M and cells were collected at 15, 30, 60, and 90 minutes (time course NaCl_1) or at 30, 45, and 60 minutes (time course NaCl_2). Experiments probing the catalytically inactive rpd3  were done in SC-leucine. The catalytically inactive rpd3 plasmid and the histone H4KQ mutant  strain were generously provided by F Posas and R Morse, respectively.
Wild-type cells were also exposed to heat shock with and without exposure to 10 μM trichostatinA (Sigma-Aldrich, St Louis, MO, USA), added 15 minutes before and throughout shock. For stress relief, cells grown at 37°C were collected by centrifugation, resuspended in 25°C YPD, and collected at 5, 10, 20, and 40 minutes (time course RH_1) or 10, 40, and 60 minutes (time course RH_2).
Microarrays and genomic analysis
Total RNA extraction, cDNA synthesis and labeling were performed as previously described [3, 68], using Superscript RT III (Invitrogen, Carlsbad, CA, USA), amino-allyl dUTP (Ambion, Austin, TX, USA) and NHS-ester cyanine dyes (Flownamics, Madison, WI, USA). Microarray data are available in the NIH Gene Expression Omnibus database with the access number [GEO:GSE9108].
Microarray data were analyzed by average-linkage hierarchical clustering, using the programs Cluster and Java-Treeview  as previously described . Genes affected in wild-type cells were defined based on triplicate single-time-point measurements [70, 71] or based on time courses  if q < 0.01 or if expression was altered more than 1.5-fold in at least two time points from replicate experiments. Genes affected in deletion strains were identified similarly, except the q-value cutoff was relaxed to 0.05.
Chromatin immunopreciptation and quantitative PCR
Rpd3-myc and Msn2p ChIP experiments were done as previously described . Briefly, cells were grown as described above and were either untreated or exposed to 0.4 mM H2O2 for 10 minutes, or washed twice with cold phosphate-buffered saline for the cold-shock control then exposed to 1% formaldehyde for 30 minutes (Rpd3-myc) or 45 minutes (Msn2 ChIPs) at 25°C. Cells were flash frozen, resuspended, and lysed; isolated chromatin was sonicated to an average size of approximately 400 bp. Protein (2.0 mg) was incubated with 5 μl anti-c-myc (9E11, Abcam (Cambridge, MA, USA) ab-56) or 15 μl anti-Msn2 (y-300, Santa Cruz (Santa Cruz, CA, USA) sc-33631) antibody overnight at 4°C. For chromatin ChIP, cells were exposed to 0.4 mM H2O2 for 20 minutes, cross-linked as above, then digested to spheroplasts with zymolyase (Seikagaku Biosystems, Tokyo, Japan) for 60 minutes at 30°C and treated with micrococcal nuclease (Worthington Biochemical, Lakewood, NJ, USA) for 20 minutes at 37°C to isolate mononucleosomes. This sample measured total nucleosome occupancy; in addition, 1.5 mg protein was mixed with 3 μl anti-acetylated H4 (Upstate 06-866 (Millipore, Billerica, MA, USA)) to immunoprecipitate acetylated histone H4. DNA purified from each sample was amplified  and converted to cDNA using SuperScript III (Invitrogen). All ChIPs were done in triplicate and quantified by real-time quantitative PCR reactions, using Sybrgreen Jumpstart Taq (Sigma-Aldrich, St Louis, MO, USA) and an Applied Biosystems 7500 detector (Foster City, CA, USA). Each ChIP PCR was normalized to a control fragment between YEL073C and YEL072W on chromosome V as previously described . Apparent histone acetylation levels were normalized to nucleosome occupancy at each locus to report the fraction of acetylated nucleosomes. Primers were designed to span approximately 75 bp regions within positioned nucleosomes  and data not shown) and were validated by amplifying genomic DNA; primer sequences are available upon request.
Additional data files
The following additional data are available with the online version of this paper: six supplemental figures (Figures S1 to S6; Additional data file 1); Table S1, showing enrichment of functional categories and transcription factor targets in gene groups taken from clustered expression data (Additional data file 2); Table S2, listing strains used in this study (Additional data file 3).
environmental stress response
Polymerase A and C
Ribosomal RNA Processing Element.
Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, Botstein D, Brown PO: Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000, 11: 4241-4257.
Causton HC, Ren B, Koh SS, Harbison CT, Kanin E, Jennings EG, Lee TI, True HL, Lander ES, Young RA: Remodeling of yeast genome expression in response to environmental changes. Mol Biol Cell. 2001, 12: 323-337.
Berry DB, Gasch AP: Stress-activated genomic expression changes serve a preparative role for impending stress in yeast. Mol Biol Cell. 2008, 19: 4580-4587.
O'Rourke SM, Herskowitz I: Unique and redundant roles for HOG MAPK pathway components as revealed by whole-genome expression analysis. Mol Biol Cell. 2004, 15: 532-542.
Gasch AP, Huang M, Metzner S, Botstein D, Elledge SJ, Brown PO: Genomic expression responses to DNA-damaging agents and the regulatory role of the yeast ATR homolog Mec1p. Mol Biol Cell. 2001, 12: 2987-3003.
Gorner W, Durchschlag E, Martinez-Pastor MT, Estruch F, Ammerer G, Hamilton B, Ruis H, Schuller C: Nuclear localization of the C2H2 zinc finger protein Msn2p is regulated by stress and protein kinase A activity. Genes Dev. 1998, 12: 586-597.
Smith A, Ward MP, Garrett S: Yeast PKA represses Msn2p/Msn4p-dependent gene expression to regulate growth, stress response and glycogen accumulation. EMBO J. 1998, 17: 3556-3564.
Beck T, Hall MN: The TOR signaling pathway controls nuclear localization of nutrient-regulated transcription factors. Nature. 1999, 402: 689-692.
Schmelzle T, Beck T, Martin DE, Hall MN: Activation of the RAS/cyclic AMP pathway suppresses a TOR deficiency in yeast. Mol Cell Biol. 2004, 24: 338-351.
Zurita-Martinez SA, Cardenas ME: Tor and cyclic AMP-protein kinase A: two parallel pathways regulating expression of genes required for cell growth. Eukaryotic Cell. 2005, 4: 63-71.
Gasch AP: The Environmental Stress Response: a Common Yeast Response to Environmental Stresses. 2002, Heidelberg: Springer-Verlag
Amoros M, Estruch F: Hsf1p and Msn2/4p cooperate in the expression of Saccharomyces cerevisiae genes HSP26 and HSP104 in a gene- and stress type-dependent manner. Mol Microbiol. 2001, 39: 1523-1532.
Boy-Marcotte E, Lagniel G, Perrot M, Bussereau F, Boudsocq A, Jacquet M, Labarre J: The heat shock response in yeast: differential regulations and contributions of the Msn2p/Msn4p and Hsf1p regulons. Mol Microbiol. 1999, 33: 274-283.
Hahn JS, Hu Z, Thiele DJ, Iyer VR: Genome-wide analysis of the biology of stress responses through heat shock transcription factor. Mol Cell Biol. 2004, 24: 5249-5256.
Rep M, Krantz M, Thevelein JM, Hohmann S: The transcriptional response of Saccharomyces cerevisiae to osmotic shock. Hot1p and Msn2p/Msn4p are required for the induction of subsets of high osmolarity glycerol pathway-dependent genes. J Biol Chem. 2000, 275: 8290-8300.
Rep M, Proft M, Remize F, Tamas M, Serrano R, Thevelein JM, Hohmann S: The Saccharomyces cerevisiae Sko1p transcription factor mediates HOG pathway-dependent osmotic regulation of a set of genes encoding enzymes implicated in protection from oxidative damage. Mol Microbiol. 2001, 40: 1067-1083.
Rep M, Reiser V, Gartner U, Thevelein JM, Hohmann S, Ammerer G, Ruis H: Osmotic stress-induced gene expression in Saccharomyces cerevisiae requires Msn1p and the novel nuclear factor Hot1p. Mol Cell Biol. 1999, 19: 5474-5485.
Proft M, Serrano R: Repressors and upstream repressing sequences of the stress-regulated ENA1 gene in Saccharomyces cerevisiae: bZIP protein Sko1p confers HOG-dependent osmotic regulation. Mol Cell Biol. 1999, 19: 537-546.
Yang XJ, Seto E: Lysine acetylation: codified crosstalk with other posttranslational modifications. Mol Cell. 2008, 31: 449-461.
Shevchenko A, Roguev A, Schaft D, Buchanan L, Habermann B, Sakalar C, Thomas H, Krogan NJ, Shevchenko A, Stewart AF: Chromatin Central: towards the comparative proteome by accurate mapping of the yeast proteomic environment. Genome Biol. 2008, 9: R167-
Carrozza MJ, Li B, Florens L, Suganuma T, Swanson SK, Lee KK, Shia WJ, Anderson S, Yates J, Washburn MP, Workman JL: Histone H3 methylation by Set2 directs deacetylation of coding regions by Rpd3S to suppress spurious intragenic transcription. Cell. 2005, 123: 581-592.
Joshi AA, Struhl K: Eaf3 chromodomain interaction with methylated H3-K36 links histone deacetylation to Pol II elongation. Mol Cell. 2005, 20: 971-978.
Keogh MC, Kurdistani SK, Morris SA, Ahn SH, Podolny V, Collins SR, Schuldiner M, Chin K, Punna T, Thompson NJ, Boone C, Emili A, Weissman JS, Hughes TR, Strahl BD, Grunstein M, Greenblatt JF, Buratowski S, Krogan NJ: Cotranscriptional set2 methylation of histone H3 lysine 36 recruits a repressive Rpd3 complex. Cell. 2005, 123: 593-605.
Rundlett SE, Carmen AA, Suka N, Turner BM, Grunstein M: Transcriptional repression by UME6 involves deacetylation of lysine 5 of histone H4 by RPD3. Nature. 1998, 392: 831-835.
Kadosh D, Struhl K: Targeted recruitment of the Sin3-Rpd3 histone deacetylase complex generates a highly localized domain of repressed chromatin in vivo. Mol Cell Biol. 1998, 18: 5121-5127.
Carrozza MJ, Florens L, Swanson SK, Shia WJ, Anderson S, Yates J, Washburn MP, Workman JL: Stable incorporation of sequence specific repressors Ash1 and Ume6 into the Rpd3L complex. Biochim Biophys Acta. 2005, 1731: 77-87. discussion 75-76.
Kadosh D, Struhl K: Repression by Ume6 involves recruitment of a complex containing Sin3 corepressor and Rpd3 histone deacetylase to target promoters. Cell. 1997, 89: 365-371.
Kurdistani SK, Robyr D, Tavazoie S, Grunstein M: Genome-wide binding map of the histone deacetylase Rpd3 in yeast. Nat Genet. 2002, 31: 248-254.
Humphrey EL, Shamji AF, Bernstein BE, Schreiber SL: Rpd3p relocation mediates a transcriptional response to rapamycin in yeast. Chem Biol. 2004, 11: 295-299.
Robert F, Pokholok DK, Hannett NM, Rinaldi NJ, Chandy M, Rolfe A, Workman JL, Gifford DK, Young RA: Global position and recruitment of HATs and HDACs in the yeast genome. Mol Cell. 2004, 16: 199-209.
Borecka-Melkusova S, Kozovska Z, Hikkel I, Dzugasova V, Subik J: RPD3 and ROM2 are required for multidrug resistance in Saccharomyces cerevisiae. FEMS Yeast Res. 2008, 8: 414-424.
De Nadal E, Zapater M, Alepuz PM, Sumoy L, Mas G, Posas F: The MAPK Hog1 recruits Rpd3 histone deacetylase to activate osmoresponsive genes. Nature. 2004, 427: 370-374.
Sertil O, Vemula A, Salmon SL, Morse RH, Lowry CV: Direct role for the Rpd3 complex in transcriptional induction of the anaerobic DAN/TIR genes in yeast. Mol Cell Biol. 2007, 27: 2037-2047.
Sharma VM, Tomar RS, Dempsey AE, Reese JC: Histone deacetylases RPD3 and HOS2 regulate the transcriptional activation of DNA damage-inducible genes. Mol Cell Biol. 2007, 27: 3199-3210.
Wang A, Kurdistani SK, Grunstein M: Requirement of Hos2 histone deacetylase for gene activity in yeast. Science. 2002, 298: 1412-1414.
Xin X, Lan C, Lee HC, Zhang L: Regulation of the HAP1 gene involves positive actions of histone deacetylases. Biochem Biophys Res Comm. 2007, 362: 120-125.
Gasch AP, Werner-Washburne M: The genomics of yeast responses to environmental stress and starvation. Funct Integr Genomics. 2002, 2: 181-192.
Sabet N, Volo S, Yu C, Madigan JP, Morse RH: Genome-wide analysis of the relationship between transcriptional regulation by Rpd3p and the histone H3 and H4 amino termini in budding yeast. Mol Cell Biol. 2004, 24: 8823-8833.
Loewith R, Smith JS, Meijer M, Williams TJ, Bachman N, Boeke JD, Young D: Pho23 is associated with the Rpd3 histone deacetylase and is required for its normal function in regulation of gene expression and silencing in Saccharomyces cerevisiae. J Biol Chem. 2001, 276: 24068-24074.
Bernstein BE, Tong JK, Schreiber SL: Genomewide studies of histone deacetylase function in yeast. Proc Natl Acad Sci USA. 2000, 97: 13708-13713.
Fazzio TG, Kooperberg C, Goldmark JP, Neal C, Basom R, Delrow J, Tsukiyama T: Widespread collaboration of Isw2 and Sin3-Rpd3 chromatin remodeling complexes in transcriptional repression. Mol Cell Biol. 2001, 21: 6450-6460.
Cherry JM, Adler C, Ball C, Chervitz SA, Dwight SS, Hester ET, Jia Y, Juvik G, Roe T, Schroeder M, Weng S, Botstein D: SGD: Saccharomyces Genome Database. Nucleic Acids Res. 1998, 26: 73-79.
Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, Jennings EG, Zeitlinger J, Pokholok DK, Kellis M, Rolfe PA, Takusagawa KT, Lander ES, Gifford DK, Fraenkel E, Young RA: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431: 99-104.
Gasch AP, Moses AM, Chiang DY, Fraser HB, Berardini M, Eisen MB: Conservation and evolution of cis-regulatory systems in ascomycete fungi. PLoS Biol. 2004, 2: e398-
Martinez-Pastor MT, Marchler G, Schuller C, Marchler-Bauer A, Ruis H, Estruch F: The Saccharomyces cerevisiae zinc finger proteins Msn2p and Msn4p are required for transcriptional induction through the stress response element (STRE). EMBO J. 1996, 15: 2227-2235.
Schmitt AP, McEntee K: Msn2p, a zinc finger DNA-binding protein, is the transcriptional activator of the multistress response in Saccharomyces cerevisiae. Proc Natl Acad Sci USA. 1996, 93: 5777-5782.
Estruch F, Carlson M: Two homologous zinc finger genes identified by multicopy suppression in a SNF1 protein kinase mutant of Saccharomyces cerevisiae. Mol Cell Biol. 1993, 13: 3872-3881.
Tavazoie S, Hughes JD, Campbell MJ, Cho RJ, Church GM: Systematic determination of genetic network architecture. Nat Genet. 1999, 22: 281-285.
Badis G, Chan ET, van Bakel H, Pena-Castillo L, Tillo D, Tsui K, Carlson CD, Gossett AJ, Hasinoff MJ, Warren CL, Gebbia M, Talukder S, Yang A, Mnaimneh S, Terterov D, Coburn D, Li Yeo A, Yeo ZX, Clarke ND, Lieb JD, Ansari AZ, Nislow C, Hughes TR: A library of yeast transcription factor motifs reveals a widespread function for Rsc3 in targeting nucleosome exclusion at promoters. Mol Cell. 2008, 32: 878-887.
Zhu C, Byers K, McCord R, Shi Z, Berger M, Newburger D, Saulrieta K, Smith Z, Shah M, Radhakrishnan M, Philippakis AA, Hu Y, De Masi F, Pacek M, Rolfs A, Murthy T, Labaer J, Bulyk ML: High-resolution DNA binding specificity analysis of yeast transcription factors. Genome Res. 2009, 19: 556-566.
Liko D, Slattery MG, Heideman W: Stb3 binds to ribosomal RNA processing element motifs that control transcriptional responses to growth in Saccharomyces cerevisiae. J Biol Chem. 2007, 282: 26623-26628.
Kasten MM, Dorland S, Stillman DJ: A large protein complex containing the yeast Sin3p and Rpd3p transcriptional regulators. Mol Cell Biol. 1997, 17: 4852-4858.
Mager WH, Planta RJ: Multifunctional DNA-binding proteins mediate concerted transcription activation of yeast ribosomal protein genes. Biochim Biophys Acta. 1990, 1050: 351-355.
Rohde JR, Cardenas ME: The tor pathway regulates gene expression by linking nutrient sensing to histone acetylation. Mol Cell Biol. 2003, 23: 629-635.
Zhao Y, McIntosh KB, Rudra D, Schawalder S, Shore D, Warner JR: Fine-structure analysis of ribosomal protein gene transcription. Mol Cell Biol. 2006, 26: 4853-4862.
Angus-Hill ML, Schlichter A, Roberts D, Erdjument-Bromage H, Tempst P, Cairns BR: A Rsc3/Rsc30 zinc cluster dimer reveals novel roles for the chromatin remodeler RSC in gene expression and cell cycle control. Mol Cell. 2001, 7: 741-751.
Damelin M, Simon I, Moy TI, Wilson B, Komili S, Tempst P, Roth FP, Young RA, Cairns BR, Silver PA: The genome-wide localization of Rsc9, a component of the RSC chromatin-remodeling complex, changes in response to stress. Mol Cell. 2002, 9: 563-573.
Lindstrom KC, Vary JC, Parthun MR, Delrow J, Tsukiyama T: Isw1 functions in parallel with the NuA4 and Swr1 complexes in stress-induced gene repression. Mol Cell Biol. 2006, 26: 6117-6129.
Rundlett SE, Carmen AA, Kobayashi R, Bavykin S, Turner BM, Grunstein M: HDA1 and RPD3 are members of distinct yeast histone deacetylase complexes that regulate silencing and transcription. Proc Natl Acad Sci USA. 1996, 93: 14503-14508.
Barbaric S, Walker J, Schmid A, Svejstrup JQ, Horz W: Increasing the rate of chromatin remodeling and gene activation--a novel role for the histone acetyltransferase Gcn5. EMBO J. 2001, 20: 4944-4951.
Clayton AL, Hazzalin CA, Mahadevan LC: Enhanced histone acetylation and transcription: a dynamic perspective. Mol Cell. 2006, 23: 289-296.
Geng F, Laurent BC: Roles of SWI/SNF and HATs throughout the dynamic transcription of a yeast glucose-repressible gene. EMBO J. 2004, 23: 127-137.
Hazzalin CA, Mahadevan LC: Dynamic acetylation of all lysine 4-methylated histone H3 in the mouse nucleus: analysis at c-fos and c-jun. PLoS Biol. 2005, 3: e393-
Metivier R, Penot G, Hubner MR, Reid G, Brand H, Kos M, Gannon F: Estrogen receptor-alpha directs ordered, cyclical, and combinatorial recruitment of cofactors on a natural target promoter. Cell. 2003, 115: 751-763.
Morillon A, Karabetsou N, Nair A, Mellor J: Dynamic lysine methylation on histone H3 defines the regulatory phase of gene transcription. Mol Cell. 2005, 18: 723-734.
Reinke H, Gregory PD, Horz W: A transient histone hyperacetylation signal marks nucleosomes for remodeling at the PHO8 promoter in vivo. Mol Cell. 2001, 7: 529-538.
Yamagoe S, Kanno T, Kanno Y, Sasaki S, Siegel RM, Lenardo MJ, Humphrey G, Wang Y, Nakatani Y, Howard BH, Ozato K: Interaction of histone acetylases and deacetylases in vivo. Mol Cell Biol. 2003, 23: 1025-1033.
Gasch AP: Yeast genomic expression studies using DNA microarrays. Methods Enzymol. 2002, 350: 393-414.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868.
Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-
Storey JD, Tibshirani R: Statistical significance for genomewide studies. Proc Natl Acad Sci USA. 2003, 100: 9440-9445.
Storey JD, Xiao W, Leek JT, Tompkins RG, Davis RW: Significance analysis of time course microarray experiments. Proc Natl Acad Sci USA. 2005, 102: 12837-12842.
Drouin S, Robert F: Genome-wide location analysis of chromatin-associated proteins by ChIP on CHIP: Controls matter. Methods.
Liu CL, Schreiber SL, Bernstein BE: Development and validation of a T7 based linear amplification for genomic DNA. BMC Genomics. 2003, 4: 19-
Lee CK, Shibata Y, Rao B, Strahl BD, Lieb JD: Evidence for nucleosome depletion at active regulatory regions genome-wide. Nat Genet. 2004, 36: 900-905.
Lee W, Tillo D, Bray N, Morse RH, Davis RW, Hughes TR, Nislow C: A high-resolution atlas of nucleosome occupancy in yeast. Nat Genet. 2007, 39: 1235-1244.
We thank A Ansari, J Lieb, R Morse, F Posas, F Robert, and T Tsukiyama and the UW-Madison Gene Expression Center for providing reagents and equipment, D Berry and E Kelley for experimental support, and W Heideman, A Ansari, J Lieb, O Rando, P Provenzano, and members of the Gasch lab for helpful discussions and critical reading of the manuscript. This work was supported by an NSF CAREER Award (#0447887) to APG and NIH NIGMS R01GM083989-01; support was provided to AAO through an NHGRI grant to the Genomic Sciences Training Program (5T32HG002760) and to DJH through an NIGMS grant to the Molecular Biosciences Training Program (T32GM007215).
AAO conducted microarray analysis, data analysis, and wrote the manuscript. DJH conducted microarray analysis, ChIP studies, microscopy, data analysis, and wrote the manuscript. MS conducted reciprocal heat shift time-courses (Figure S5 in Additional data file 1), and SN carried out acquired stress experiments (Figure S6 in Additional data file 1). DP and JLW assisted in experimental setup, RNA preparation, and strain construction. APG carried out data analysis and wrote the manuscript.
Adriana L Alejandro-Osorio, Dana J Huebert contributed equally to this work.
Electronic supplementary material
Additional data file 1: Figure S1: Basal expression in rpd3Δ cells does not account for stress-dependent expression defects. Gene expression in unstressed rpd3Δ versus wild type (average of triplicate experiments) and for wild-type and rpd3Δ cells responding to stress is shown as in Figure 1. The middle panel shows the differences between wild-type and rpd3Δ expression where yellow represents higher transcript abundance (that is, weaker repression) and blue indicates lower transcript abundance (that is, weaker induction) in the rpd3Δ mutant. The right panel shows the difference between transcript abundance in wild-type and rpd3Δ cells after adjusting for the basal expression differences in unstressed cells. These data show that the observed defect in ESR initiation in the rpd3Δ strain is not due to the possibility that basal expression of the ESR genes already reflects the 'ON' state of the program. First, genes repressed in the ESR show little discernable difference in basal expression in rpd3Δ versus wild-type cells. Second, although a subset of the iESR genes are subtly derepressed (approximately 1.5-fold) in untreated rpd3Δ cells, the mutant still displays lower absolute transcript levels relative to wild-type at the peak of the expression response (right panel). Third, the rpd3Δ mutant ultimately alters expression comparable to the expression changes in wild-type cells, resulting in higher absolute transcript levels for these iESR genes in acclimated rpd3Δ cells compared to acclimated wild-type cells, particularly in response to heat shock and H2O2 treatment (right panel and data not shown). Finally, almost half the iESR genes show no significant difference in basal expression (within 1.3-fold of wild-type and P > 0.01) but are still induced to lower peak levels than in wild-type cells. Thus, the defect in stress-dependent expression changes seen in the rpd3Δ strain is not accounted for by basal expression differences across all ESR genes, although these results suggest that Rpd3p is required to suppress some iESR genes in the absence of stress. Figure S2: Msn2 localizes to the nucleus upon stress treatment in cells lacking Rpd3p. Msn2p localization in wild-type and rpd3Δ cells responding to H2O2 treatment was scored, in cells transformed with a plasmid constitutively expressing Msn2-green fluorescent (GFP) protein obtained from T Tsukiyama [6, 58]. (a) Percent of cells with nuclear Msn2p localization at several time points after H2O2 treatment in wild-type and rpd3Δ cells according to the key. (b) Examples of cytoplasmic and nuclear Msn2-GFP before and after stress. Nuclear Msn2p is indicated with an arrow. Figure S3: Rpd3p catalytic activity and modifiable histone H4 are required for ESR expression in response to stress. The average expression and standard deviation of genes in the PS, RP, and iESR gene groups is shown as cells responded to (a) 25°C to 37°C heat shock or (b) 0.4 mM H2O2 treatment. Plots represent the response of rpd3Δ cells harboring plasmid-borne RPD3 (left), the blank vector (middle), or the catalytically inactive allele rpd3-H150:151A (right) according to the key. (c) Gene expression was also measured in wild-type cells responding to 25°C to 37°C heat shock, with and without pretreatment with 10 μM trichostatin A. The average log2 expression change of genes in each group is plotted for treated and untreated cells. Time points with smaller expression changes in the trichostatin A-treated cells (P < 0.01, paired t-test) are indicated with an asterisk. The effect of trichostatin A is less severe than RPD3 deletion or mutation (a, b), although still statistically significant, likely due to incomplete inhibition of Rpd3p by the drug. Each plot represents the standard deviation of expression of genes in the group (not error). (d) Expression of the PS, RP, and iESR genes is shown for wild-type, rpd3Δ, pho23Δ, and H4KQ cells responding to 0.4 mM H2O2 treatment, as described in Figure 1. Each column shows the average of at least biological triplicates. The difference in expression between the isogenic wild-type and each mutant strain is shown in the right panel as described in Figure 1. Figure S4: The Rpd3L subunit PHO23 is required for transient initiation of the ESR. Average expression of iESR, PS, and RP genes in wild-type, rpd3Δ, pho23Δ, and rco1Δ cells responding to 0.4 mM H2O2 treatment at 10, 20, 30, 40, and 60 minutes. Figure S5: Rpd3p is required for suppressed ESR expression during stress relief. (a) Gene expression diagrams and (b) average expression plots are shown as described in Figure 1 for wild-type and rpd3Δ cells responding to 25°C to 37°C heat shock (left panels) or the reciprocal shift from 37°C to 25°C (right panels). The starting point of the recovery experiment is indicated by a dashed line and was adjusted to the acclimated expression levels seen in the heat shock experiment for clarity. Figure S6: PHO23 mutant cells display a defect in acquired stress resistance. Previous work from our lab showed that the gene expression response to a single dose of stress is not required to survive that condition, but rather protects cells against future stress . Here, H2O2 tolerance in wild-type and pho23Δ cells was measured after a 60-minute pretreatment with 0.7 M NaCl, similarly to that previously described . Briefly, log-phase cells were treated with 0.7 M NaCl (or YPD for mock-treated cells) for 60 minutes, then washed once with YPD to remove salt and exposed to 23 doses of H2O2 (ranging from 0 to 8 mM) for 2 h. Cell viability at each dose was measured using live-dead staining (Molecular Probes (Invitrogen, Carlsbad, CA, USA) and a Guava flow cytometer. The plot shows the growth score (the sum of all viability scores normalized to the starting viability score in untreated cells) for wild-type and the pho23Δ mutant. Error bars represent the standard deviation from triplicate experiments. The pho23Δ strain was not sensitive to the mild dose of NaCl in this assay (P = 0.4) but had a significant defect in acquired tolerance to H2O2 after mild NaCl treatment (P = 0.029, paired t-test). The defect is similar to that seen for cells lacking MSN2 and/or MSN4 . (PDF 531 KB)
Additional data file 2: Enrichment of functional categories and transcription factor targets in gene groups taken from clustered expression data. (XLS 124 KB)
About this article
Cite this article
Alejandro-Osorio, A.L., Huebert, D.J., Porcaro, D.T. et al. The histone deacetylase Rpd3p is required for transient changes in genomic expression in response to stress. Genome Biol 10, R57 (2009). https://doi.org/10.1186/gb-2009-10-5-r57
- Additional Data File
- Gene Expression Change
- Stress Relief
- H2O2 Treatment
- Transient Phase