![]() |
|
|
Vol. 20, Issue 15, 3503-3513, August 1, 2009
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

*Department of Molecular Biology, Princeton University, Princeton, NJ 08544; and
Department of Physics and Astronomy and BioMaPS Institute for Quantitative Biology, Rutgers University, Piscataway, NJ 08854
Submitted February 6, 2009;
Revised April 28, 2009;
Accepted May 26, 2009
Monitoring Editor: William P. Tansey
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Chromatin can also play a more passive role in gene regulation by simply controlling which transcription factor binding sites are available for participation in gene regulation. Recent genome-wide studies of nucleosome positioning found that most yeast promoters contain an extended nucleosome-depleted region (NDR) that is enriched for transcription factor binding sites (Lee et al., 2004
; Yuan et al., 2005
; Lee et al., 2007
; Mavrich et al., 2008
; Shivaswamy et al., 2008
). Moreover, sequences to which a transcription factor can bind are often more prevalent in the genome than the sites to which that factor is actually bound in vivo. For example, the yeast genome contains many more high-affinity binding motifs for Rap1 and Leu3 transcription factors than are bound by the cognate factor in vivo (Lieb et al., 2001
; Liu et al., 2006
). In addition, not all genes with a given transcription factor binding site show correlated expression patterns. These discrepancies may be explained by selective occlusion of transcription factor binding sites by positioned nucleosomes. For example, Leu3 sites to which Leu3 is not bound in vivo tend to lie under well positioned nucleosomes (Liu et al., 2006
). Thus, as Morse has proposed, chromatin may be "instructive" for transcription factor binding (Morse, 2007
).
To distinguish active versus passive roles of chromatin in regulating gene expression, Elgin and coworkers postulated two general categories of promoter architecture (Lu et al., 1994
; Wallrath et al., 1994
): "pre-set" promoters whose nucleosome structure does not significantly change with changes in transcription, such as the Drosophila hsp26 and hsp70 promoters; and "remodeling" promoters whose nucleosome structure changes significantly with gene expression changes, such as the mammalian mouse mammary tumor virus long terminal repeat promoter and the yeast PHO5 promoter. To assess whether this dichotomy applies to most yeast promoters and, if so, which yeast promoters fall into each category, we have explored the positions of nucleosomes across the entire yeast genome under conditions in which a majority of yeast genes undergo significant change in gene expression. Several recent studies have reported the precise nucleosome positions across the yeast genome (Yuan et al., 2005
; Lee et al., 2007
; Mavrich et al., 2008
; Shivaswamy et al., 2008
; Kaplan et al., 2009
). Most studies of genome-wide nucleosome positioning have been conducted at steady states by using cells grown in rich media. However, the Segal laboratory observed that most nucleosome positions remain unchanged in several alternative carbon sources (Kaplan et al., 2009
). Shivaswamy et al. (2008)
examined nucleosome positioning before and after heat shock and observed discrete local nucleosome remodeling events at gene promoters but found no general relationship between transcriptional changes and nucleosome remodeling events. In contrast, Lee et al. (2004)
also examined nucleosome occupancy changes during heat shock, by using lower resolution techniques, and concluded that promoter nucleosome occupancy changed inversely with transcriptional change.
We determined nucleosome occupancy before and after addition of glucose to cells grown in a poor carbon source, a nutrient upshift in which more than half of all yeast genes significantly change expression (Wang et al., 2004
; Zaman et al., 2009
). Accordingly, we have been able to correlate changes in nucleosome occupancy to changes in transcriptional activity in a large number of genes and relate those changes to the presence of previously predicted transcription factor binding motifs. We find that most genes that undergo significant expression changes do not exhibit any observable change in nucleosome positioning within their promoters, although those that do are significantly enriched for TATA boxes. Moreover, we find quite often that the accessibility of specific transcription factor binding motifs within a promoter before the nutrient upshift determines whether the corresponding gene will exhibit a change in expression. Thus, our results ratify Morse's postulate of a largely instructive role for chromatin in regulating gene expression.
| MATERIALS AND METHODS |
|---|
|
|
|---|
30 min). We washed spheroplasts twice with 40 ml of 0°C spheroplasting buffer, gently resuspended them in 5 ml of NP-buffer (50 mM NaCl, 10 mM Tris, pH 7.4, 1 M sorbitol, 5 mM MgCl2, 1 mM CaCl2, with 1 mM β-mercaptoethanol and 500 µM spermidine added fresh), and then added Nonidet-P40 at a final concentration of 0.025%. We added micrococcal nuclease (10 U; Sigma-Aldrich, St. Louis, MO) and incubated the lysate for 45 min at 37°C, digestion conditions in which chromatin was converted predominantly to mononucleosomes (Supplemental Figure S1). We stopped digestions by addition of 1 ml of 5% SDS and 50 mM EDTA, followed by addition of 10 µl of 10 mg/ml RNaseA and incubation at 37°C for 30 min. We then added 5 µl of 10 mg/ml proteinase K and reversed the cross-links by incubation at 65°C for at least 4 h. DNA was purified using a polymerase chain reaction (PCR) clean-up kit (QIAGEN, Valencia, CA) and then analyzed by gel electrophoresis to assess the extent of digestion. We performed this experiment for four biological replicates, all of which yielded essentially identical results. Data for one of the replicates are presented in the Results.
Microarray Hybridization and Data Analysis
DNA (8 µg) was labeled and hybridized to 1.0R yeast tiling arrays (Affymetrix, Santa Clara, CA) as described previously (Gresham et al., 2006
), except that fragmentation was accomplished by digestion with 0.05 U of DNaseI (Invitrogen, Carlsbad, CA) for 2 min at 37°C. Raw microarray data for individual time points was analyzed with Affymetrix Tiling Analysis Software, version 1.1 with settings as described previously (Lee et al., 2007
). The microarray data are available for download from http://puma.princeton.edu/publications.html. Yeast mRNA expression data have been published previously (Zaman et al., 2009
).
A Probabilistic Algorithm for Predicting Nucleosome Positions from Tiled Microarray Data
We developed a hidden Markov model (HMM) (Rabiner, 1989
) to infer genome-wide nucleosome positions from hybridization data (Supplemental Figure S2). The model specifies nucleosome positions as peaks of probability that DNA will enter a nucleosome at a given base pair in the genome. The model also returns the overall nucleosome occupancy of any base pair in the genome as the probability that it is covered by any nucleosome. To compute the nucleosome occupancy for a base pair, we sum over posterior probabilities for all nucleosomes that initiate in the appropriate direction within 148 base pairs (length of 37 probes separated by 4 base pairs) of that site (Morozov et al., 2008). Hybridization values from the tiled Affymatrix arrays were used as input to the HMM, and the parameters of the model were fit by maximum likelihood. Because consecutive probes are not always equidistant on the chip, all log intensities were first placed on a regular four-base pair grid by linear interpolation. Gaps of 50 base pairs or longer were assigned zero hybridization values for the HMM but omitted from all subsequent bioinformatic analysis. The parameters of the HMM are fit on a set of 19 genes that show no appreciable changes in expression between 0 and 20 min (POL4, CDC20, RNR3, RAM2, RSM22, CDC73, ASC1, SCP160, ILV3, NNF1, AKR1, SKI2, MDL1, SAS4, SNF5, URE2, YFB3, SRL2, and MEC3), and their associated 5' intergenic regions. The model uses a mixture of two Gaussians to represent both nucleosome state and linker state. The final set of fitting parameters includes means and widths of four one-dimensional Gaussians, two independent mixture coefficients, and 38 initial probabilities for the nucleosome and linker states. The topology of the Bayesian network (Supplemental Figure S2) reflects the fact that all 37 nucleosome states are forced to share the same mixture of two Gaussians. The HMM parameters inferred independently for all 19 loci were averaged and the resulting model was run genome-wide on each data set. We implemented the HMM in Matlab (The MathWorks, Natick, MA) as a dynamic Bayesian network using standard libraries from the Bayes Net Toolbox (BNT; http://www.cs.ubc.ca/
murphyk/Software/BNT/bnt.html). Our Matlab code is available upon request.
| RESULTS |
|---|
|
|
|---|
To facilitate interpretation of the nucleosome positioning data, we developed a novel HMM to predict nucleosome occupancy from hybridization intensity data. Although HMMs have been previously used to predict nucleosome positions (Yuan et al., 2005
; Lee et al., 2007
), our approach has several important distinctions. First, all nucleosomes have a canonical length of 148 base pairs (37 probes separated by 4 base pairs). Second, the probability of starting a new nucleosome (PLN in Supplemental Figure S2) is not fit by maximum likelihood but rather is adjusted manually to reproduce physiologically relevant nucleosome occupancies in bulk chromatin (Wolffe, 1998
). The manual adjustment reflects our lack of control over the zero intensity baseline, which can shift depending on the relative amounts of DNA used in the nucleosomal and control samples. Third, we model the probability of observing a nucleosome or linker state by a mixture of two Gaussians. We sought to account for the fact that a single Gaussian will cause the assignment of small nucleosome probabilities to high log intensity values. This problem is resolved by using a mixture model which can in principle fit functions of arbitrary shape, although in practice the number of Gaussian components has to be limited to two or three to reduce the total number of fitting parameters. In sum, our model avoids the need for post hoc assumptions of dubious biological accuracy, such as nucleosomes of noncanonical length, whereas, as described below, maintaining accurate predictions of experimentally determined in vivo nucleosome positions.
To test our HMM, we applied it to our hybridization data for a set of seven reference loci, independently of those used for training, with previously determined nucleosome positions. This stringent analysis indicated that interpretation of our data by this model places nucleosomes with statistically significant accuracy (Supplemental Figure S3) similar to that of a previously published study (Lee et al., 2007
). We also applied the HMM to our hybridization data at the CHA1 gene, whose nucleosome organization has been previously determined both by indirect end labeling and by global nucleosome positioning studies (Moreira and Holmberg, 1998
; Yuan et al., 2005
; Lee et al., 2007
). Our HMM places nucleosomes at locations determined by previous studies (Figure 1). CHA1 expression is glucose independent, and we find that our calculated nucleosome occupancy at this locus remains unchanged before glucose addition and 20 and 60 min post-glucose addition (Supplemental Figure S4). From these analyses, we conclude that our data sets, in conjunction with the HMM we developed to interpret those data, yield results highly consistent with prior nucleosome position determinations, and we conclude that both the model and the data are reliable.
|
To further validate our data and methodology, we examined nucleosome occupancy at SUC2 and ADH2, two glucose-regulated genes that undergo well-characterized glucose-induced changes in nucleosome positioning. Glucose addition represses SUC2 expression by fivefold. Twenty minutes after glucose addition, we detect a shift in a promoter nucleosome to a site further upstream from the ORF. Sixty minutes after glucose addition, we detect addition of a new nucleosome adjacent to the shifted nucleosome (Figure 2), as has been reported previously (Perez-Ortin et al., 1987
; Hirschhorn et al., 1992
; Gavin and Simpson, 1997
). Glucose addition represses ADH2 gene expression >50-fold. Our analysis indicates the addition of a well-positioned nucleosome to the ADH2 promoter 20 min after glucose addition (Supplemental Figure S6), as has been observed in previous studies (Verdone et al., 1996
). In addition, the array of nucleosomes extending into the ADH2 ORF shifts toward the promoter upon repression. This repressed promoter structure is retained 60 min after glucose addition. For both ADH2 and SUC2, we compared our nucleosome positions with previously published positions determined for cells grown in glucose (Verdone et al., 1996
; Gavin and Simpson, 1997
). The correlation between the published positions and our post-glucose nucleosome occupancy is quite high, and clearly higher than for our pre-glucose positions (Supplemental Figure S7), indicating that our nucleosome position determinations are consistent with those reported previously.
|
|
|
Of the genes that exhibit changes in promoter nucleosome occupancy, repressed genes tend to gain nucleosomes and induced genes tend to lose promoter nucleosomes. However, across all promoters the relationship between changes in nucleosome occupancy and transcription is surprisingly weak, with an overall correlation of 0.34. Similar results are obtained 60 min after glucose addition. Moreover, the same result is obtained by plotting changes in expression versus changes in nucleosome occupancy as predicted by our HMM rather than changes in hybridization intensity (Supplemental Figure S10). To visualize changes in nucleosome occupancy across all promoters, we calculated the difference in nucleosome occupancy between 0 and 20 min, relative to the transcription start site (TSS), at each promoter and then grouped these by K-means clustering into three groups. This clustering separated promoters into those that gained nucleosomes, those that lost nucleosomes, and those that remained unchanged (Figure 4B). Consistent with the results in Figure 4A, only 5% of all promoters gained nucleosomes and only 5% lost nucleosomes, despite that 50% of all genes changed expression by more than twofold. Moreover, for most of those promoters that gained or lost nucleosomes, only a single nucleosome was affected and most changes occurred between –550 and –100 base pairs from the transcriptional start site. In sum, 77 of 578 genes repressed by more than fourfold gained promoter nucleosomes, whereas 39 of 398 genes induced by more than fourfold lost promoter nucleosomes. However, the majority of promoters underwent no significant change in nucleosome occupancy profiles.
Our analysis of the relationship between changes in promoter nucleosome occupancy and gene expression measures the latter by mRNA microarrays. However, mRNA microarrays measure steady-state mRNA levels, which can be altered by changes in mRNA turnover as well as changes in transcription rate. To assess whether the dearth of nucleosome remodeling we observe on glucose upshift derives from a limited change in transcription rates of most genes, we performed chromatin immunoprecipitation (ChIP) against RNA polymerase II in ORFs (Supplemental Figure S11). At the SUC2 and HXT1 genes, changes in promoter nucleosome occupancy are associated with changes in mRNA levels (Figure 4A); after glucose upshift, we observe an increase in RNA polymerase II detected by ChIP at the HXT1 gene and a decrease at the SUC2 gene, which in both cases correlates with the changes in their mRNA levels. Examining several genes with glucose-regulated changes in mRNA levels but no changes in promoter nucleosome structure (HSP104, ALD3, REX4, and DHR2), we observed in every case a change in the occupancy of RNA polymerase II that correlated strongly with the change in mRNA level. These data indicate that genes with glucose-regulated changes in mRNA levels do undergo associated changes in transcription initiation rates, even in the absence of changes in promoter nucleosome occupancy.
Because we observed relatively few changes in nucleosome occupancy with transcriptional change, we considered the possibility that glucose-regulated promoters might alter their nucleosome positioning without altering their net nucleosome occupancy. To address this possibility, we determined the average nucleosome occupancy profiles across the promoters of genes whose transcription is either induced or repressed by glucose, before and after glucose addition (Supplemental Figure S12). The promoter nucleosome profiles of both glucose-induced and repressed genes showed little change, demonstrating that glucose-regulated genes do not undergo coherent nucleosome movements. Examination of nucleosome changes at individual promoters similarly failed to reveal examples of genes regulated by large-scale nucleosome sliding (Figure 4B). Because this level of analysis would not capture small shifts in nucleosome positioning that might alter the ability of transcription factors to bind, we specifically examined the relation between nucleosome position and putative transcription factor binding motifs. These results, described below, are consistent with the above observations that nucleosome repositioning plays a limited role in altering gene expression in response to glucose upshift.
We also examined changes in nucleosome occupancy at the 3' end of transcripts. Difference maps of nucleosome occupancy at the 3' end of all genes revealed a small number of genes that undergo nucleosome loss or addition at the 3' end (Supplemental Figure S13A). Genes that exhibit nucleosome occupancy changes at their 3' end are equally divided between those convergent with their downstream neighboring gene and those collinear with the downstream gene, indicating that the nucleosome remodeling at the 3' end cannot be explained solely as remodeling of a downstream promoter. However, these changes in nucleosome occupancy at the 3' end are not correlated with changes in gene expression (r = –0.01) (Supplemental Figure S13B). Thus, the function of 3' NDR and changes to that structure remain unresolved.
Promoters Containing TATA Boxes Are More Likely to Undergo Nucleosome Remodeling
We observed that genes containing TATA boxes, which constitute
20% of all genes in yeast (Basehoar et al., 2004
), were unusual in several respects. First, as noted previously (Ioshikhes et al., 2006
; Albert et al., 2007
; Mavrich et al., 2008
), we found that, although the 5' NDRs in non-TATA containing promoters were all 200 base pairs wide and stereotypically centered 100 base pairs upstream of the TSS, the NDRs in promoters with TATA boxes were of the same size as those at non-TATA promoters but were dispersed over a much wider region, with each gene containing an individually positioned NDR. This pattern is evident in the heat map of nucleosome positions in individual promoters, shown in Figure 3, with the TATA-containing genes significantly overrepresented in cluster a, and in plot of the average nucleosome density as a function of distance from the TSS (Supplemental Figure S14). Second, we observed that glucose-repressed genes were significantly enriched for TATA box-containing promoters (71% contain TATA boxes as opposed to a genome-wide average of 22%; binomial p = 1.3 x 10–36), whereas glucose-induced genes were significantly depleted for TATA box-containing promoters (15% contain TATA boxes; binomial p = 3.0 x 10–3). In sum, TATA box-containing genes behave differently from other genes both in nucleosome positioning and in expression changes after a glucose upshift.
Because the asymmetry in the distribution of TATA box containing genes between glucose-induced and repressed genes mirrored the asymmetry in nucleosome remodeling in induced versus repressed genes, we directly examined if promoters containing a TATA box were more likely to undergo nucleosome remodeling. The promoter nucleosome profile for glucose repressed genes was subdivided between genes with or without a TATA box (Figure 5A). Strikingly, TATA box-containing promoters gained considerable nucleosome occupancy during gene repression, whereas promoters lacking TATA boxes gained little nucleosome occupancy during repression. We also subdivided the promoter nucleosome profile of glucose-induced genes between those with or without a TATA box (Figure 5B). Most glucose-induced genes lack TATA boxes and showed little difference in their promoter nucleosome profile before or after induction, whereas the minority of glucose-induced genes with TATA boxes lost nucleosome occupancy after induction. Thus, promoters of genes that contain TATA boxes tend to undergo changes in nucleosome occupancy upon changes in transcriptional activity, whereas those genes lacking TATA boxes tend to retain the same nucleosome profile upon changes in transcriptional activity.
|
Functional Transcription Factor Binding Motifs Reside in Regions of Low Nucleosome Occupancy
Our data allowed us to examine the relationship among transcription factor binding motifs, nucleosome occupancy, and transcriptional changes. In particular, we explored the degree to which nucleosome occupancy over transcription factor binding motifs changed as a function of changes in expression or influenced their use in promoting transcriptional change. Genes regulated by the glucose-repressive transcription factors Adr1 and Cat8 (Tachibana et al., 2005
) gain promoter nucleosomes during glucose repression (Verdone et al., 1996
; Agricola et al., 2004
). Consistent with these prior observations, we observed that many Adr1/Cat8-regulated genes gained promoter nucleosomes on transcriptional repression (Figure 4A). However, the overall correlation between changes in promoter nucleosomes and transcription for these genes is weak (r = 0.2), although this correlation is higher if alternative measures of nucleosome occupancy are used (Supplemental Figure S10).
We extended our analysis to determine whether shifts in nucleosome positions expose or occlude transcription factor binding motifs in conjunction with changes in gene expression. We previously applied a mutual information based algorithm, FIRE (Elemento et al., 2007
), to identify promoter-localized sequence motifs that mediate glucose-induced changes in gene expression, many of which corresponded to binding sites for transcription factors known to be involved in the glucose response, such as Msn2/4 and Hap4 (Zaman et al., 2009
). Limiting our analysis to promoters, we measured differences in nucleosome occupancy in response to glucose addition directly over all instances of these individual motifs and observed little change in their nucleosome occupancy upon glucose addition (data not shown). To reduce the effect of spurious, nonfunctional representatives of these motifs, we calculated the difference in nucleosome occupancy over occurrences of each motif at 20 and 60 min after glucose addition only for promoters of those genes that exhibited transcriptional induction or repression (Figure 6, A and B). Most motifs undergo no significant change in average nucleosome occupancy. However, a small number of motifs, such as the Hap1-related CCG(N)5CC motif, become exposed on average during induction or, such as the Gln3 binding motif, become occluded during repression. This bias remains even after excluding those promoters that undergo overall changes in nucleosome occupancy during the upshift (data not shown). Examining the actual distribution of individual nucleosome occupancy changes over motifs, rather than the average changes, reveals that most motif occurrences undergo no significant changes in nucleosome occupancy, but for some motifs a subpopulation undergoes substantial changes in nucleosome occupancy (Supplemental Figure S16). Thus, at a limited number of promoters subtle shifts in nucleosome positioning can occasionally lead to the exposure or occlusion of transcription factor binding motifs.
|
Promoters of the glucose-induced ribosomal biogenesis genes are highly enriched for RRPE and PAC motifs (Beer and Tavazoie, 2004
; Jorgensen et al., 2004
; Wade et al., 2006
). Quite strikingly, in the promoters of those genes containing a PAC or RRPE motif (labeled Motif B in Figure 6) that were subsequently induced by glucose addition, the motifs were significantly more exposed than they were in those genes that were subsequently repressed or remained unchanged upon glucose addition, suggesting that likely physiologically relevant occurrences of RRPE/PAC motifs may be free of nucleosomes and accessible, whereas likely spurious instances of the RRPE/PAC motifs in promoters of repressed or unchanged genes are more likely to be contained within nucleosomes and perhaps inaccessible. The Msn2/4 transcription factors mediate portions of the general stress response which is alleviated by glucose addition (Gasch et al., 2000
; Zaman et al., 2009
), suggesting that physiologically relevant occurrences of the Msn2/4 motifs would be expected to occur in the promoters of glucose-repressed genes. We find that occurrences of the Msn2/4 motifs in promoters of glucose-repressed genes are significantly more exposed than in genes that were subsequently induced or transcriptionally unchanged. Several other motifs exhibited this dichotomy in occupancy, although many others did not. Nonetheless, these results suggest that accessibility dictated by nucleosome occupancy likely presages the functional activity of several transcription factor binding motifs.
| DISCUSSION |
|---|
|
|
|---|
Remodeling Promoters
The promoters of only 10% of all genes gain or lose nucleosomes after glucose addition despite that >50% of all genes exhibit a change in mRNA levels of twofold or more. Those genes that gain promoter nucleosomes are predominantly those that undergo repression, whereas those that lose nucleosomes are predominantly unchanged in expression or, to a lesser extent, undergo induction. Those promoters that undergo nucleosome remodeling—both those that lose as well as those that gain nucleosomes—are highly enriched for predicted TATA boxes, with approximately half of all remodeled promoters containing a TATA box. As a consequence, TATA box-containing promoters exhibit a much higher correlation between transcriptional change and nucleosome remodeling than do promoters lacking TATA boxes. Pugh and coworkers noted previously that TATA boxes are enriched in the promoters of genes induced by starvation and other cell stresses and that these genes are particularly dependent on chromatin modification activities (Basehoar et al., 2004
). We can now appreciate this dependency as a consequence of the significant proportion of TATA-containing promoters, relative to TATA-less promoters, that undergo nucleosome remodeling during transcriptional reprogramming. However, we still do not understand why TATA-containing boxes are associated with promoter remodeling. The TATA box itself does not seem to be required for nucleosome remodeling, as mutant PHO5, SUC2, and CUP1 promoters that lack TATA boxes still undergo wild type nucleosome remodeling (Hirschhorn et al., 1992
; Fascher et al., 1993
; Shen et al., 2001
). Moreover, artificial recruitment of TBP to several remodeling promoters is not sufficient to initiate nucleosome remodeling (Ryan et al., 2000
). Thus, although the correlation between TATA boxes and chromatin remodeling is clear, the mechanistic basis for that correlation is not.
In addition to remodeling at TATA-containing promoters, half the promoters that undergo remodeling under our conditions do not contain TATA boxes. What determines whether the promoter of a gene will undergo nucleosome remodeling during transcriptional reprogramming? We were unable to find any significant genome-wide correlation between binding motifs for particular transcription factors and nucleosome remodeling. Similarly, we were unable to detect specific motifs enriched at sites of nucleosome removal or addition (data not shown). One possibility is that remodeling is limited to those promoters in which the intrinsic energetics of nucleosome position, dictated by the ease at which segments of DNA can conform to the structure dictated by a nucleosome, does not strongly favor a single chromatin conformation. For many genes at which remodeling does not take place, the two nucleosomes bracketing the NDR (the +1 and –1 nucleosomes) assume their positions largely as a result of the intrinsic DNA sequence within the promoter (Ioshikhes et al., 2006
; Segal et al., 2006
). Such promoters may not be readily reconfigured, whereas other promoters in which multiple distinct sites are equally energetically favorable may be more easily remodeled. The availability from this study of an extensive list of remodeling versus nonremodeling promoters, in conjunction with improved computational methods to determine the free energy associated with nucleosome binding to any specific sequence of DNA, should allow us to directly test this hypothesis.
Preset Promoters
The majority of genes transcriptionally induced or repressed by glucose do not alter their promoter nucleosome structure. These preset promoters tend to be underrepresented for TATA boxes and tend to have a strong nucleosome-depleted region at a stereotypical position relative to the transcriptional start site, a promoter architecture that has been termed "type II" (Field et al., 2008
). Previous work has characterized TATA-less promoters as enriched in constitutively expressed housekeeping genes (Basehoar et al., 2004
), but the work reported here clearly demonstrates that individual TATA-less preset promoters are capable of dramatic gene expression changes.
The unaltered promoter nucleosome structure for most glucose-regulated genes implies the existence of constitutively accessible binding sites for the factors that control expression of these genes. In agreement with this inference, previous studies mapping genome-wide nucleosome positions at steady state have discovered a strong enrichment of predicted transcription factor binding sites in NDRs (Yuan et al., 2005
; Lee et al., 2007
; Morse, 2007
; Field et al., 2008
). Extending these studies, we found that gene promoters with an accessible transcription factor binding motif are significantly more likely to be regulated like known targets of that transcription factor than promoters where the motif is occluded by a nucleosome. This correlation was particularly noteworthy for induction through the RRPE/PAC elements and repression through Msn2/4 binding motifs. Because we did not directly measure transcription factor binding to these motifs, we cannot formally conclude that nucleosomes occlude transcription factor binding at any given motif; however, these results are consistent with previous work that found significant differences between the predicted binding sites and actual in vivo binding profile of the Leu3 and Rap1 transcription factors (Lieb et al., 2001
; Liu et al., 2006
). Given these observations, we would anticipate that coupling nucleosome position data with clustered gene expression data should improve computational approaches to predicting regulatory motifs and transcription factor binding sites associated with specific transcriptional responses (Beer and Tavazoie, 2004
; Elemento et al., 2007
). In sum, our results suggest that while nucleosome repositioning does not seem to be a major driving force in yeast for altering transcriptional activity, nucleosome positioning seems to play a significant role in determining which transcription factor binding sites are available for directing transcriptional changes.
Unexpected Promoter Nucleosome Stability in Yeast
Despite the massive transcriptional reprogramming induced by glucose addition we observed surprisingly little change in nucleosome positioning, particularly across coding regions. Moreover, a significant majority of yeast promoters show little or no change in nucleosome occupancy in association with significant change in expression. This stability was unanticipated given previous work on promoters for genes such as PHO5, ADH2, SUC2, and CHA1, where gene induction leads to substantial promoter nucleosome remodeling. We also observe the previously reported nucleosome remodeling in the promoters of the ADH2 and SUC2 genes, but our genome-wide analysis of promoter nucleosomes indicates that these remodeling promoters cannot be considered typical. Thus, the majority of yeast promoters fall into the category of promoters Elgin and colleagues described as preset.
Our analysis is limited to nucleosome positioning and does not examine many other changes that may occur at promoter nucleosomes, such as altered histone modifications, changing incorporation of histone variants into positioned nucleosomes, or even histone turnover. In addition, we have only examined nucleosome positioning changes in response to glucose addition, leaving the formal possibility that other perturbations might induce nucleosome remodeling to a greater degree. However, both our work and recent work of the Iyer laboratory, which examined global nucleosome positioning changes in response to heat shock (Shivaswamy et al., 2008
), demonstrate remarkably stable nucleosome structure for most genes in yeast.
| ACKNOWLEDGMENTS |
|---|
| Footnotes |
|---|
Address correspondence to: James R. Broach (jbroach{at}princeton.edu) or Alexandre V. Morozov (morozov{at}physics.rutgers.edu)
Abbreviations used: HMM, hidden Markov model; NDR, nucleosome-depleted region; ORF, open reading frame.
| REFERENCES |
|---|
|
|
|---|
Albert, I., Mavrich, T. N., Tomsho, L. P., Qi, J., Zanton, S. J., Schuster, S. C., and Pugh, B. F. (2007). Translational and rotational settings of H2A.Z nucleosomes across the Saccharomyces cerevisiae genome. Nature 446, 572–576.[CrossRef][Medline]
Almer, A., Rudolph, H., Hinnen, A., and Horz, W. (1986). Removal of positioned nucleosomes from the yeast PHO5 promoter upon PHO5 induction releases additional upstream activating DNA elements. EMBO J 5, 2689–2696.[Medline]
Basehoar, A. D., Zanton, S. J., and Pugh, B. F. (2004). Identification and distinct regulation of yeast TATA box-containing genes. Cell 116, 699–709.[CrossRef][Medline]
Beer, M. A., and Tavazoie, S. (2004). Predicting gene expression from sequence. Cell 117, 185–198.[CrossRef][Medline]
Durrin, L. K., Mann, R. K., and Grunstein, M. (1992). Nucleosome loss activates CUP1 and HIS3 promoters to fully induced levels in the yeast Saccharomyces cerevisiae. Mol. Cell. Biol 12, 1621–1629.
Elemento, O., Slonim, N., and Tavazoie, S. (2007). A universal framework for regulatory element discovery across all genomes and data types. Mol Cell 28, 337–350.[CrossRef][Medline]
Fascher, K. D., Schmitz, J., and Horz, W. (1993). Structural and functional requirements for the chromatin transition at the PHO5 promoter in Saccharomyces cerevisiae upon PHO5 activation. J. Mol. Biol 231, 658–667.[CrossRef][Medline]
Fedor, M. J., and Kornberg, R. D. (1989). Upstream activation sequence-dependent alteration of chromatin structure and transcription activation of the yeast GAL1-GAL10 genes. Mol. Cell. Biol 9, 1721–1732.
Field, Y., Kaplan, N., Fondufe-Mittendorf, Y., Moore, I. K., Sharon, E., Lubling, Y., Widom, J., and Segal, E. (2008). Distinct modes of regulation by chromatin encoded through nucleosome positioning signals. PLoS Comput. Biol 4, e1000216.[CrossRef][Medline]
Gasch, A. P., Spellman, P. T., Kao, C. M., Carmel-Harel, O., Eisen, M. B., Storz, G., Botstein, D., and Brown, P. O. (2000). Genomic expression programs in the response of yeast cells to environmental changes. Mol. Biol. Cell 11, 4241–4257.
Gavin, I. M., and Simpson, R. T. (1997). Interplay of yeast global transcriptional regulators Ssn6p-Tup1p and Swi-Snf and their effect on chromatin structure. EMBO J 16, 6263–6271.[CrossRef][Medline]
Gresham, D., Ruderfer, D. M., Pratt, S. C., Schacherer, J., Dunham, M. J., Botstein, D., and Kruglyak, L. (2006). Genome-wide detection of polymorphisms at nucleotide resolution with a single DNA microarray. Science 311, 1932–1936.
Han, M., and Grunstein, M. (1988). Nucleosome loss activates yeast downstream promoters in vivo. Cell 55, 1137–1145.[CrossRef][Medline]
Han, M., Kim, U. J., Kayne, P., and Grunstein, M. (1988). Depletion of histone H4 and nucleosomes activates the PHO5 gene in Saccharomyces cerevisiae. EMBO J 7, 2221–2228.[Medline]
Hirschhorn, J. N., Brown, S. A., Clark, C. D., and Winston, F. (1992). Evidence that SNF2/SWI2 and SNF5 activate transcription in yeast by altering chromatin structure. Genes Dev 6, 2288–2298.
Ioshikhes, I. P., Albert, I., Zanton, S. J., and Pugh, B. F. (2006). Nucleosome positions predicted through comparative genomics. Nat. Genet 38, 1210–1215.[CrossRef][Medline]
Jorgensen, P., Rupes, I., Sharom, J. R., Schneper, L., Broach, J. R., and Tyers, M. (2004). A dynamic transcriptional network communicates growth potential to ribosome synthesis and critical cell size. Genes Dev 18, 2491–2505.
Kaplan, N. et al. (2009). The DNA-encoded nucleosome organization of a eukaryotic genome. Nature 458, 362–366.[CrossRef][Medline]
Kent, N. A., Bird, L. E., and Mellor, J. (1993). Chromatin analysis in yeast using NP-40 permeabilised sphaeroplasts. Nucleic Acids Res 21, 4653–4654.
Kim, Y., McLaughlin, N., Lindstrom, K., Tsukiyama, T., and Clark, D. J. (2006). Activation of Saccharomyces cerevisiae HIS3 results in Gcn4p-dependent, SWI/SNF-dependent mobilization of nucleosomes over the entire gene. Mol. Cell. Biol 26, 8607–8622.
Lee, C. K., Shibata, Y., Rao, B., Strahl, B. D., and Lieb, J. D. (2004). Evidence for nucleosome depletion at active regulatory regions genome-wide. Nat. Genet 36, 900–905.[CrossRef][Medline]
Lee, W., Tillo, D., Bray, N., Morse, R. H., Davis, R. W., Hughes, T. R., and Nislow, C. (2007). A high-resolution atlas of nucleosome occupancy in yeast. Nat. Genet 39, 1235–1244.[CrossRef][Medline]
Li, B., Carey, M., and Workman, J. L. (2007). The role of chromatin during transcription. Cell 128, 707–719.[CrossRef][Medline]
Lieb, J. D., Liu, X., Botstein, D., and Brown, P. O. (2001). Promoter-specific binding of Rap1 revealed by genome-wide maps of protein-DNA association. Nat. Genet 28, 327–334.[CrossRef][Medline]
Liu, X., Lee, C. K., Granek, J. A., Clarke, N. D., and Lieb, J. D. (2006). Whole-genome comparison of Leu3 binding in vitro and in vivo reveals the importance of nucleosome occupancy in target site selection. Genome Res 16, 1517–1528.
Lu, Q., Wallrath, L. L., and Elgin, S. C. (1994). Nucleosome positioning and gene regulation. J. Cell. Biochem 55, 83–92.[CrossRef][Medline]
Mavrich, T. N., Ioshikhes, I. P., Venters, B. J., Jiang, C., Tomsho, L. P., Qi, J., Schuster, S. C., Albert, I., and Pugh, B. F. (2008). A barrier nucleosome model for statistical positioning of nucleosomes throughout the yeast genome. Genome Res 18, 1073–1083.
Moreira, J. M., and Holmberg, S. (1998). Nucleosome structure of the yeast CHA1 promoter: analysis of activation-dependent chromatin remodeling of an RNA-polymerase-II-transcribed gene in TBP and RNA pol II mutants defective in vivo in response to acidic activators. EMBO J 17, 6028–6038.[CrossRef][Medline]
Morozov, A. V., Fortney, K., Gaykalova, D. A., Studitsky, V. M., Widom, J., and Siggia, E. D. (2009). Using DNA mechanics to predict in vitro nucleosome positions and formation energies. Nucleic Acids Res in press.
Morse, R. H. (2003). Getting into chromatin: how do transcription factors get past the histones? Biochem. Cell Biol 81, 101–112.[CrossRef][Medline]
Morse, R. H. (2007). Transcription factor access to promoter elements. J. Cell. Biochem 102, 560–570.[CrossRef][Medline]
Nagalakshmi, U., Wang, Z., Waern, K., Shou, C., Raha, D., Gerstein, M., and Snyder, M. (2008). The transcriptional landscape of the yeast genome defined by RNA sequencing. Science 320, 1344–1349.
Perez-Ortin, J. E., Estruch, F., Matallana, E., and Franco, L. (1987). Fine analysis of the chromatin structure of the yeast SUC2 gene and of its changes upon derepression. Comparison between the chromosomal and plasmid-inserted genes. Nucleic Acids Res 15, 6937–6956.
Rabiner, L. (1989). A tutorial on Hidden Markov Models and selected applications in speech recognition. Proc. IEEE 77, 257–286.[CrossRef]
Ryan, M. P., Stafford, G. A., Yu, L., and Morse, R. H. (2000). Artificially recruited TATA-binding protein fails to remodel chromatin and does not activate three promoters that require chromatin remodeling. Mol. Cell. Biol 20, 5847–5857.
Segal, E., Fondufe-Mittendorf, Y., Chen, L., Thastrom, A., Field, Y., Moore, I. K., Wang, J. P., and Widom, J. (2006). A genomic code for nucleosome positioning. Nature 442, 772–778.[CrossRef][Medline]
Shen, C. H., Leblanc, B. P., Alfieri, J. A., and Clark, D. J. (2001). Remodeling of yeast CUP1 chromatin involves activator-dependent repositioning of nucleosomes over the entire gene and flanking sequences. Mol. Cell. Biol 21, 534–547.
Shivaswamy, S., Bhinge, A., Zhao, Y., Jones, S., Hirst, M., and Iyer, V. R. (2008). Dynamic remodeling of individual nucleosomes across a eukaryotic genome in response to transcriptional perturbation. PLoS Biol 6, e65.[CrossRef][Medline]
Tachibana, C., Yoo, J. Y., Tagne, J. B., Kacherovsky, N., Lee, T. I., and Young, E. T. (2005). Combined global localization analysis and transcriptome data identify genes that are directly coregulated by Adr1 and Cat8. Mol. Cell. Biol 25, 2138–2146.
Verdone, L., Camilloni, G., Di Mauro, E., and Caserta, M. (1996). Chromatin remodeling during Saccharomyces cerevisiae ADH2 gene activation. Mol. Cell. Biol 16, 1978–1988.
Wade, C. H., Umbarger, M. A., and McAlear, M. A. (2006). The budding yeast rRNA and ribosome biosynthesis (RRB) regulon contains over 200 genes. Yeast 23, 293–306.[CrossRef][Medline]
Wallrath, L. L., Lu, Q., Granok, H., and Elgin, S. C. (1994). Architectural variations of inducible eukaryotic promoters: preset and remodeling chromatin structures. Bioessays 16, 165–170.[CrossRef][Medline]
Wang, Y., Pierce, M., Schneper, L., Guldal, C. G., Zhang, X., Tavazoie, S., and Broach, J. R. (2004). Ras and Gpa2 mediate one branch of a redundant glucose signaling pathway in yeast. PLoS Biol 2, E128.[CrossRef][Medline]
Williams, S. K., and Tyler, J. K. (2007). Transcriptional regulation by chromatin disassembly and reassembly. Curr. Opin. Genet. Dev 17, 88–93.[CrossRef][Medline]
Wolffe, A. (1998). Chromatin: Structure and Function, San Diego, CA: Academic Press.
Yuan, G. C., Liu, Y. J., Dion, M. F., Slack, M. D., Wu, L. F., Altschuler, S. J., and Rando, O. J. (2005). Genome-scale identification of nucleosome positions in S. cerevisiae. Science 309, 626–630.
Zaman, S., Lippman, S. I., Schneper, L., Slonim, N., and Broach, J. R. (2009). Glucose regulates transcription in yeast through a network of signaling pathways. Mol. Syst. Biol 5, 245.[Medline]
This article has been cited by other articles:
![]() |
C. Kristell, J. Orzechowski Westholm, I. Olsson, H. Ronne, J. Komorowski, and P. Bjerling Nitrogen depletion in the fission yeast Schizosaccharomyces pombe causes nucleosome loss in both promoters and coding regions of activated genes Genome Res., March 1, 2010; 20(3): 361 - 371. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. I. Lippman and J. R. Broach Protein kinase A and TORC1 activate genes for ribosomal biogenesis by inactivating repressors encoded by Dot6 and its homolog Tod6 PNAS, November 24, 2009; 106(47): 19928 - 19933. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||