Skip to main content

Global nucleosome occupancy in yeast

Abstract

Background

Although eukaryotic genomes are generally thought to be entirely chromatin-associated, the activated PHO5 promoter in yeast is largely devoid of nucleosomes. We systematically evaluated nucleosome occupancy in yeast promoters by immunoprecipitating nucleosomal DNA and quantifying enrichment by microarrays.

Results

Nucleosome depletion is observed in promoters that regulate active genes and/or contain multiple evolutionarily conserved motifs that recruit transcription factors. The Rap1 consensus was the only binding motif identified in a completely unbiased search of nucleosome-depleted promoters. Nucleosome depletion in the vicinity of Rap1 consensus sites in ribosomal protein gene promoters was also observed by real-time PCR and micrococcal nuclease digestion. Nucleosome occupancy in these regions was increased by the small molecule rapamycin or, in the case of the RPS11B promoter, by removing the Rap1 consensus sites.

Conclusions

The presence of transcription factor-binding motifs is an important determinant of nucleosome depletion. Most motifs are associated with marked depletion only when they appear in combination, consistent with a model in which transcription factors act collaboratively to exclude nucleosomes and gain access to target sites in the DNA. In contrast, Rap1-binding sites cause marked depletion under steady-state conditions. We speculate that nucleosome depletion enables Rap1 to define chromatin domains and alter them in response to environmental cues.

Background

Global gene-expression patterns are established and maintained by the concerted actions of transcription factors and the proteins that constitute chromatin. The global network of interactions between transcription factors and promoters in yeast is increasingly being characterized [1]. The role of chromatin in gene regulation is less clear, however. For example, the distribution of nucleosomes, the fundamental units of chromatin, is poorly understood on a gene-specific basis, much less a global basis [2].

The nucleosome consists of approximately 146 base-pairs (bp) of DNA wrapped around an octamer of histone proteins - two each of histones H2A, H2B, H3 and H4. Eukaryotic genomes are packaged into repeating units of nucleosomes separated by around 10-80 bp of linker DNA. High occupancy by nucleosomes is thought to be generally repressive [3], and extensive remodeling (and loss) of nucleosomes occurs in the promoters of genes undergoing activation [4]. In the case of the PHO5 promoter in yeast, this remodeling proceeds until essentially no nucleosomes are detected across a region of several hundred base-pairs [5, 6].

Transcription factors and chromatin proteins each form complex regulatory networks that interact in a variety of ways [1, 7]. Transcription factors modify chromatin structure by recruiting enzymes that remodel nucleosomes or posttranslationally modify histones (by acetylation or methylation, for example) [8–10]. The modifications can be maintained through cell division and propagated to proximal nucleosomes by positive-feedback mechanisms [7, 11, 12]. Hence, a signal such as the activation of a transcription factor can be temporally and spatially transmitted through chromatin. Conversely, chromatin can influence transcription factor function by modulating the accessibility of target binding sites in the DNA [13, 14].

We used chromatin immunoprecipitation (ChIP) and DNA microarrays to evaluate nucleosome occupancy levels for essentially all promoters in yeast. Promoters that regulate active genes, contain multiple conserved motifs or recruit Rap1 tend to be relatively nucleosome-depleted. We also used real-time PCR and micrococcal nuclease digestion to show that nucleosomes are depleted in the vicinity of Rap1 consensus sites. This depletion can be partially reversed by the actions of the small molecule rapamycin or by removing Rap1-binding sites. We suggest that other transcription factors have less robust nucleosome-depleting activities than Rap1 and must therefore act collaboratively to gain access to their cognate sites in the DNA.

Results

ChIP-based assay for nucleosome occupancy

Histones are essential components of the nucleosome and efficiently cross-link to nucleosomal DNA. Antibodies against invariant portions of histones have been used previously in ChIP assays to follow nucleosome loss at the yeast PHO5 promoter [5, 6]. We extended this approach to evaluate relative nucleosome occupancy at essentially all promoters and other intergenic regions in yeast. DNA associated in vivo with histone H3 was isolated by ChIP using antibody against the carboxy terminus of histone H3 (no posttranslational modifications are thought to occur in this region). ChIP DNA and unenriched control DNA were amplified by in vitro transcription and evaluated using microarrays. DNA associated with histone H2B was evaluated in a similar fashion using anti-FLAG antibody and a FLAG-H2B strain. H3 and H2B datasets were compiled by averaging four and three independent biological experiments, respectively. These datasets are remarkably similar as shown by a genome-wide correlation of 0.83 (Figure 1a-c). This correlation is comparable to that observed when comparing replicate H3 datasets (or H2B datasets), and suggests that both assays measure similar phenomena. In the H3 and H2B datasets, respectively, there are 347 and 214 regions depleted at least 1.5-fold relative to the average over all intergenics. In contrast, there are just 84 and 6 regions in the respective datasets enriched at least 1.5-fold relative to this average. The relatively narrow range of ChIP enrichment and the negative skew of the data (Figure 1b) are consistent with the conventional view that the majority of the genome is packaged into nucleosomes with intervening stretches of free DNA such as the activated PHO5 promoter [5, 6].

Figure 1
figure 1

Correlation between H3 and FLAG-H2B ChIP datasets. DNA associated with histones in vivo was enriched in ChIP assays using antibodies against histone H3 or FLAG-H2B, and quantified by microarrays. (a) Relative enrichment of promoters and other non-coding regions in the H3 and H2B ChIP assays is shown. (b) Histogram showing distributions of enrichment for promoter regions in the H3 and H2B ChIP assays. (c) Overlap between regions depleted in the H3 and FLAG-H2B assays is shown. Overall, there is an 0.83 correlation between these ChIP datasets. (d) Overlap between regions depleted in the H3 ChIP assay and regions enriched by aqueous extraction is shown [62].

Despite these consistencies, a possible caveat to using ChIP to evaluate nucleosome occupancy is that immunoprecipitation efficiency can depend on epitope accessibility. Rather than having low occupancies, genomic regions depleted in the H3 ChIP might be inaccessible as a result of association with large protein complexes in chromatin. To investigate this possibility, we examined a published chromatin fractionation dataset in which cross-linked chromatin fragments were subjected to phenol-chloroform extraction and DNA that partitioned into the aqueous phase was quantified by microarrays [15]. Given the polar nature of DNA and the hydrophobic nature of denatured protein, aqueous extraction should generally enrich for free DNA. We found that regions depleted in the H3 ChIP assay overlap extensively with regions enriched by aqueous extraction, but not with regions depleted by aqueous extraction (Figure 1d). Overall, there is a negative correlation of -0.54 between the H3 ChIP and aqueous-extraction datasets. Although the fractionation data may partially reflect differential cross-linking of lysines in the histone tails [15], this analysis suggests that regions depleted in the H3 ChIP experiment are relatively protein-free, as would be expected of non-nucleosomal DNA.

Nucleosome occupancy correlates inversely with promoter strength

As previous studies show that PHO5 activation is accompanied by marked nucleosome loss in the promoter region [5, 6], we sought to determine whether nucleosome depletion is a general attribute of active promoters. A total of 4,365 intergenic regions that reside immediately upstream of one or more validated yeast genes were assigned as promoters. Relative transcription rates were determined for each yeast gene from transcript levels measured by array and previously collected mRNA half-life data [16]. We found an inverse correlation of -0.39 between the enrichment of promoters in the H3 and H2B ChIP assays and the transcription rates of downstream genes (Figure 2a). Under the conditions examined, PHO5 is not induced and its promoter has an average nucleosome occupancy according to these datasets. To evaluate further the relationship between nucleosome depletion and transcription, we collated a set of 308 nucleosome-depleted promoters on the basis of their relative depletion across the replicate H3 and H2B experiments. Of these nucleosome-depleted promoters, 42% regulate highly active genes (Figure 2b). These data suggest that there is a systematic relationship between promoter strength and nucleosome depletion. However, as this correspondence is not complete there are likely to be other determinants of nucleosome occupancy.

Figure 2
figure 2

Inverse association between nucleosome occupancy and promoter strength. (a) Relative enrichment of promoter regions in the H3 and FLAG-H2B ChIP assays plotted against transcription rate of downstream genes (moving average, window 50). (b) Overlap between promoters upstream of active genes and the set of nucleosome-depleted promoters defined on the basis of depletion across the replicate H3 and FLAG-H2B experiments.

Transcription factor binding motifs are over-represented in nucleosome-depleted promoters

To identify additional determinants of occupancy, we sought sequence elements associated with nucleosome depletion. Specifically, we carried out an unbiased search for elements up to 10 bp in length that occur with higher frequency in nucleosome-depleted promoters. Two distinct categories of sequences emerged (Figure 3a). The first includes poly(dA.dT) elements. Stretches of 10 or more dA.dT nucleotides appear in 38% of depleted promoters, compared with 26% of promoters overall (hypergeometric p < 10-5). dA.dT stretches destabilize nucleosome formation in vitro and in vivo [17, 18]. The enrichment of poly(dA.dT) elements in nucleosome-depleted promoters probably reflects, at least in part, this destabilizing influence. As a high proportion of the poly(dA.dT) elements identified in nucleosome-depleted promoters are more than 10 bp long (30% are at least 14 bp), these data do not address the minimum length required for destabilization. However, in vitro studies show that a 16-bp insertion leads to a 1.7-fold increase in accessibility of nucleosomal target sites [18].

Figure 3
figure 3

Sequence motifs over-represented in nucleosome-depleted promoters. (a) An unbiased search for sequences up to 10 bp in length over-represented in nucleosome-depleted promoters (relative to promoters overall) identified the poly(dA.dT) sequence element and variants of the Rap1 consensus motif ACACCCATACAT [21]. (b) Overlap between nucleosome-depleted promoters and promoters that contain multiple conserved motifs is shown [21]. (c) Histogram showing average numbers of motifs in 1,000 randomly generated promoter sets. Nucleosome-depleted promoters contain an average of 6.1 conserved motifs, significantly higher than in these randomly generated sets.

The second sequence element enriched in nucleosome-depleted promoters corresponds to the consensus motif for the Rap1 transcription factor. This motif commonly occurs in the promoters of ribosomal proteins genes and is required for Rap1 binding in vitro and in vivo [19, 20]. Some variant of this motif appears in 22% of nucleosome-depleted promoters, compared with just 8% of promoters overall (hypergeometric p < 10-5). Furthermore, multiple Rap1 sites are found in 19% of nucleosome-depleted promoters with Rap1 sites, compared to 8% of promoters with Rap1 sites overall (hypergeometric p < 10-3). These data suggest that Rap1 recruitment may lead to nucleosome loss.

Because only the Rap1 consensus site was identified in an unbiased search, we sought to identify additional sequence motifs by incorporating species conservation data. Specifically, we evaluated a set of 71 conserved motifs identified by Kellis and colleagues, a majority of which function in transcription factor recruitment [21]. Nearly half of these 71 motifs are over-represented in nucleosome-depleted promoters relative to promoters overall, as defined by a hypergeometric p < 0.001. However, many of the implicated motifs appear in the same promoters. For example, nine of the over-represented motifs are associated with filamentation gene promoters [21]. We therefore considered the possibility that the total number of conserved motifs might be a more relevant predictor of nucleosome depletion. Indeed, we found that 31% of nucleosome-depleted promoters contain at least eight motifs, compared with 11% of promoters overall (hypergeometric p <10-5; Figure 3b). Furthermore, nucleosome-depleted promoters contain an average of 6.1 motifs, whereas the average promoter contains 3.1 (permutation p < 0.001; Figure 3c). Next, we sought motifs associated with nucleosome depletion in the absence of multiple motifs, by confining our analysis to promoters containing a maximum of four motifs. This analysis identified just two over-represented motifs, which correspond to the Rap1 and Swi4 binding sites. Hence, although a large number of conserved motifs are enriched in nucleosome-depleted promoters, most appear to be relevant mainly when occurring in combination.

Functionally cooperative transcription factors associate with nucleosome-depleted promoters

As a majority of the conserved motifs recruit transcription factors [21], we examined the relationship between transcription factor binding and nucleosome occupancy more directly. Lee and colleagues combined ChIP and microarrays to identify target promoters for essentially all yeast transcription factors under the same conditions used here to evaluate nucleosome occupancy [1]. For each factor, we determined the significance of overlap between its target promoters and the set of nucleosome-depleted promoters. Of the 113 transcription factors in their database, 31 tend to associate with nucleosome-depleted promoters as defined by a hypergeometric p < 0.001. Rap1 has the most significant association (Figure 4a), consistent with the enrichment of its binding motif (see above). Other top-ranked factors include Fhl1, which associates with many Rap1-bound promoters, and Swi4, whose binding motif is also enriched (Table 1).

Figure 4
figure 4

Nucleosome depletion in the vicinity of Rap1-binding sites. (a) Overlap between the 308 most nucleosome-depleted promoters and promoters found to recruit Rap1 in a global ChIP study [1]. (b) Nucleosome depletion in the vicinity of Rap1-binding sites in ribosomal gene promoters evaluated by ChIP. Fold-enrichment was determined by real-time PCR using primers that span Rap1-binding motifs in the RPS22A, RPS15, RPS11B and RPL23A promoters. (c) Southern blots showing DNA from yeast spheroplasts digested with increasing concentrations of micrococcal nuclease probed with labeled PCR products spanning the TUB2 promoter and the Rap1 sites in the RPS11B and RPS15 promoters. (d) Nucleosome occupancy for a mutant RPS11B promoter lacking Rap1 consensus sites was determined by H3 ChIP and real-time PCR. The mutant promoter is enriched 2.1-fold relative to wild type. (e) Nucleosome occupancy at Rap1-binding sites in ribosomal protein gene promoters after treatment with rapamycin evaluated by H3 ChIP and real-time PCR.

Table 1 Transcription factors that tend to associate with nucleosome-depleted promoters

We sought an underlying binding mechanism or function common to the transcription factors we had identified. However, these factors utilize a variety of binding domains, regulate different pathways, and only a minority have significant associations with promoters of highly active genes. Nonetheless, a commonality does emerge when transcription factor cooperativity is considered. A recent informatics study by Banerjee and Zhang identified 31 functionally cooperative transcription factor pairs (representing a total of 33 factors) on the basis of comprehensive binding and expression data [22]. Only a fraction of these are known to interact physically, suggesting that other mechanisms also confer cooperative function. There is a remarkable correspondence between these functionally cooperative factors and those that preferentially associate with nucleosome-depleted promoters (see Table 1). Of the 31 factors we found to associate with nucleosome-depleted promoters, 17 were found to be functionally cooperative by Banerjee and Zhang (p < 10-5). Furthermore, an evaluation of nucleosome occupancy at promoters bound by both members of a cooperative pair revealed a significant association with nucleosome-depletion for 18 of the 31 pairs (hypergeometric p < 0.01). Together, these findings suggest that binding motifs and transcription factors act in combination to deplete nucleosomes and suggest a role for nucleosomes in transcription factor cooperativity [23–25].

Conditional nucleosome depletion at Rap1 consensus motifs

Although a number of transcription factors appear to act in defining promoter nucleosome occupancy, only the Rap1 consensus motif was identified in an unbiased search of nucleosome-depleted promoters. Furthermore, there is a highly significant association between nucleosome-depleted promoters and promoters bound by this factor in vivo [1] (Figure 4a). To investigate the relationship between Rap1 recruitment and nucleosome depletion further, we used ChIP and real-time PCR to evaluate nucleosome occupancy at several Rap1 binding sites in ribosomal protein promoters. We found that these regions are depleted 3- to 10-fold in H3 and FLAG-H2B ChIP assays, relative to a control promoter (TUB2) with average occupancy by global analysis (Figure 4b). We also used an orthogonal approach in which micrococcal nuclease digestion [26] was used to probe for nucleosomes at the TUB2, RPS11B and RPS15 promoters (Figure 4c). A pattern of nuclease protection indicative of a regular nucleosome array is evident at the TUB2 promoter, consistent with the average nucleosome occupancy attributed to this promoter by global ChIP analysis. In contrast, nuclease protection is not evident at the RAP1 sites in the RPS15 promoter, consistent with the marked nucleosome-depletion attributed to this region by global ChIP and real-time PCR analysis. The region surrounding the RAP1 sites in RPS11B exhibits weak nuclease protection, consistent with the modest nucleosome-depletion attributed to this region by global ChIP and real-time PCR. Although these focused analyses specifically addressed Rap1 sites in ribosomal protein genes, our global analyses indicate that approximately 30% of nucleosome-depleted promoters containing Rap1 motifs do not regulate ribosomal protein genes. Together these data confirm that nucleosomes are markedly depleted in the vicinity of Rap1 consensus sites in vivo, and thus extend previous studies showing that Rap1 induces local alterations in chromatin structure that, for example, result in increased nuclease sensitivity [27–29].

To gain further insight into the relationship between Rap1 and nucleosome depletion, we examined a mutant RPS11B promoter lacking its Rap1 consensus sites. We found that removal of these sites, which completely abrogates Rap1 binding [30], causes nucleosomes to return to the region, as reflected by a greater than twofold change in H3 ChIP enrichment (Figure 4d). We also examined the effect of rapamycin treatment on nucleosome occupancy in the vicinity of these consensus sites. Although ribosomal protein gene expression is dramatically reduced by rapamycin [31, 32], Rap1 remains bound to its target promoters ([30, 33], and B.B., E.P. and S.S., unpublished results). We found that rapamycin treatment causes nucleosomes to return to the vicinity of Rap1 sites, as reflected by twofold and greater increases in H3 ChIP enrichment (Figure 4e). Together these data show that Rap1 consensus sites are required for conditional nucleosome depletion at ribosomal protein gene promoters.

Discussion

To gain further insight into the role of nucleosomes in gene regulation, we systematically evaluated promoter nucleosome occupancy in yeast by immunoprecipitating nucleosomal DNA and quantifying enrichment with microarrays. Promoters that are inefficiently immunoprecipitated by general anti-histone antibodies, and are therefore presumed to be relatively nucleosome-depleted, tend to regulate active genes (Figure 2). This is consistent with the previous observation that the activated PHO5 promoter is largely devoid of nucleosomes [5, 6]. However, as not all nucleosome-depleted promoters regulate active genes, there are most likely to be additional determinants of depletion. An unbiased search for sequence elements enriched in nucleosome-depleted promoters revealed poly(dA.dT) elements, previously shown to destabilize nucleosome formation [17, 18], and the Rap1 consensus motif. By incorporating sequence conservation data [21], more than 30 other enriched motifs could be identified. However, most of these appear to be relevant mainly when occurring in combination. When we limited this analysis to promoters containing four or fewer motifs, all but two of these additional motifs drop out (only the Rap1 and Swi4 consensus sites remain). As the majority of conserved motifs incorporated in this analysis recruit transcription factors [21], these data suggest that multiple transcription factors act in combination to deplete nucleosomes. This possibility is further supported by our finding that functionally cooperative transcription factors tend to bind nucleosome-depleted promoters. These associations may reflect a mechanistic model in which transcription factors compete collaboratively to displace nucleosomes in order to gain access to target sites in the DNA [23]. This model was formulated to explain why certain pairs of transcription factors bind cooperatively to proximal target sites in vivo and on a chromatin template, but not to naked DNA [23–25]. This view invokes a broad role for nucleosomes as ubiquitous negative regulators of transcription factor binding and function. We speculate that by promoting synergy among multiple transcription factors and impeding the activities of individual ones, nucleosomes facilitate threshold behavior and filter noise (for example, genetic variation in motif sequence) in the transcriptional regulatory network.

Although many factors appear to act in defining promoter nucleosome occupancy, our data indicate that Rap1 has a uniquely important role. Rap1 and its consensus motif are both markedly enriched in nucleosome-depleted promoters. Follow-up studies using real-time PCR and micrococcal nuclease digestion also demonstrate marked nucleosome depletion in the vicinity of Rap1 sites in the promoters of ribosomal protein genes. Moreover, nucleosomes appeared to return when the Rap1 consensus sites in one of these promoters were removed. These findings are consistent with previously described roles for Rap1 in opening chromatin and altering nucleosome positioning [27, 28]. However, Rap1 recruitment is not equally associated with nucleosome depletion under all conditions. We find that nucleosomes partially return to the vicinity of Rap1 sites during a rapamycin-induced starvation response [34], even though Rap1 remains bound ([30, 33], and B.B, E.P. and S.S., unpublished results). Hence, the nucleosome loss associated with Rap1 recruitment is most likely to require additional proteins, such as Esa1, a histone acetyltransferase recruited by Rap1 under exponential growth conditions but released in stress [30].

These findings may also offer insight into the barrier activity previously documented for Rap1 [35]. Heterochromatin propagation involves the sequential modification of histones in adjacent nucleosomes through positive-feedback mechanisms [7, 11]. Certain factors such as Rap1 are able to block this propagation by largely unknown mechanisms [36]. One model speculates that these barriers create nucleosome-free 'holes' lacking the histone substrate required for heterochromatin propagation [29, 35]. By identifying such a 'hole' in the vicinity of Rap1-binding sites in vivo our data support this model. Remarkably, the nucleosomal hole and the barrier function ascribed to Rap1 may be conditional, as nucleosomes return following treatment with the small molecule rapamycin, which activates a starvation response. Heterochromatic silencing has been shown previously to moderate under these conditions [37]. Hence, we speculate that dynamic influences on nucleosome occupancy may enable Rap1 to define chromatin domains and vary them in response to environmental cues.

More broadly, the widespread nucleosome loss observed in the promoters of active genes provides a general caveat for ChIP studies examining posttranslational histone modifications, as a decrease in signal for a histone modification at a promoter undergoing activation may actually reflect nucleosome loss. Similarly, regions that appear relatively hypo-modified by ChIP may actually be nucleosome-depleted. However, this is not the case for low levels of acetylation [38] and H3 lysine 4 methylation [39] observed at yeast telomeres, as these regions have high occupancy. The data also provide insight into the maintenance of epigenetic information by histone modifications. Whereas epigenetic memory of a repressed state can be maintained on histones in promoters, memory of an activated state must be maintained on histones outside the promoters, for example in transcribed regions, which may not undergo significant nucleosome loss during activation [5, 6]. Methylation of histone H3 at lysines 4 and 36, targeted to transcribed regions in yeast via interactions between RNA polymerase and the methylases [39–47], may represent such 'activating' marks.

Materials and methods

Chromatin immunoprecipitation (ChIP)

DNA associated with histone H3 in vivo was immunoprecipitated with antibodies against the invariant H3 carboxy terminus using a ChIP protocol described previously [39, 48, 49]. Briefly, 45 ml log-phase w303a yeast (OD600 ~ 1.0) growing in yeast extract/peptone/dextrose (YPD) were cross-linked in 1% formaldehyde for 15 min, washed twice in PBS, resuspended in 400 μl lysis buffer (50 mM Hepes-KOH pH 7.5, 140 mM NaCl, 1 mM EDTA, 1% Triton X-100, 0.1% sodium deoxycholate) and lysed with glass beads. The resulting extract was sonicated to fragment chromatin (4 × 20 sec burst/30 sec rest with a Branson Sonifier 250 at 70% duty, power 3) and centrifuged for 15 min. Solubilized chromatin was then immunoprecipitated with polyclonal antibodies against the carboxy terminus of histone H3 (Abcam or Cell Signaling). A unenriched whole-cell extract sample (WCE) was also retained as a control. After enrichment, cross-links were reversed by incubating samples in 10 mM Tris-HCl pH 8.0, 1 mM EDTA, 1.0% SDS, 150 mM NaCl at 65°C overnight. DNA was purified from ChIP and WCE samples by proteinase K treatment, phenol/chloroform extraction, ethanol precipitation, and incubation with RNAse. DNA associated with histone H2B in vivo was isolated in a similar manner from yeast containing epitope-tagged H2B [50] using anti-FLAG M2 monoclonal antibodies (Sigma).

DNA amplification and hybridization

To obtain sufficient quantities for hybridization, immunoprecipitated DNA (from approximately 108 cells) and whole-cell extract DNA (unenriched control) were amplified in a linear fashion as described [51]. Briefly, terminal transferase was used to add poly(T) tails to DNA fragment and a T7-poly(A) adaptor primer was used to incorporate T7 promoters. The reaction products were used as template for an in vitro transcription reaction carried out with the T7 Megascript Kit (Ambion) and RNA samples were purified using an RNeasy Mini Kit (Qiagen). Amplified RNA was reverse-transcribed, incorporating amino-allyl dUTP, and the resulting DNA was fluorescently labeled by incubation with monofunctional reactive Cy5 (enriched sample) or Cy3 (unenriched control) dye as described [52]. Microarrays containing 6,438 PCR-amplified intergenic regions were prepared as described previously [39, 53, 54]. Mixed Cy5-/Cy3-labeled probe was hybridized to intergenic microarrays for 12-14 h at 60°C, washed and then scanned using a GenePix 4000A scanner with GenePix Pro software (Axon Instruments) as described [55]. In addition, transcript levels were determined by hybridizing Cy5-labeled mRNA extracted from log phase w303a yeast against Cy3-labeled genomic DNA on microarrays containing 6,218 open reading frames (ORFs), as described previously [16].

Microarray data processing

Cy5 and Cy3 fluorescence were integrated for each feature using GenePix Pro Software (Axon). Data were processed and composite Cy5:Cy3 ratios determined according to protocols at the Stanford Microarray Database [56]. Correlations between replicate datasets were ~0.8 for all experiments. Composite datasets were log2 transformed and zero centered before further analysis. The histone H3 ChIP dataset was determined from four independent immunoprecipitations and hybridizations (two each using antibodies from Cell Signaling or Abcam). The FLAG-H2B ChIP dataset was determined from three independent immunoprecipitations and hybridizations. The mRNA dataset was determined from three independent extractions and hybridizations of mRNA against genomic DNA. Relative transcription rates were determined by dividing transcript levels by half-life data collected by Wang and colleagues [16]. A set of activated promoters was defined as those in the top 10% by mRNA expression level of associated gene, with divergent promoters assigned to the more highly expressed gene. Complete datasets are available online [57].

Analysis of nucleosome-depleted promoters

Z-scores were assigned to each intergenic that reflect depletion across the four H3 and three H2B ChIP experiments, using the formula Z = (x - μ)/σ where x is the average of the replicate measurements, μ is the average of all intergenics and σ is the standard error of the replicate measurements. We defined as nucleosome-depleted the 410 features with the highest Z-scores. This set, which includes 308 promoters, contains nucleosome-depleted outliers and is not inclusive of all promoters that immunoprecipitate with average or lower efficiency. The average aqueous enrichment ratio [15] for these 308 depleted promoters is 1.7-fold, significantly higher than expected by chance (permutation p < 0.001), consistent with the premise that these promoters are relatively free of nucleosomes.

Sequence elements common to nucleosome-depleted promoters were identified by searching between 10 and 500 bp upstream of gene start sites for over-represented sequences up to 10 bp in length using the GeneSpring program suite (Silicon Genetics). Enrichment was confirmed by evaluating the significance of overlap between the set of nucleosome-depleted promoters and the set of promoters containing Rap1 consensus motifs (ACACCCATACAT with up to two mismatches) or poly dA.dT stretches at least 10 bp in length (identified using PatMatch, Saccharomyces Genome Database [58]). Statistical significances of overlaps between sets are expressed as P-values calculated by a hypergeometric probability model. The P-values reflect the extent to which observed overlaps exceed that expected under the null hypothesis that there is no relationship between the sets [59]. Where specified, permutation analyses were carried out by generating 1,000 random but representative promoter sets with an Excel macro and used to confirm statistical significance. Lists of promoters containing the 71 conserved motifs [21] were collated from gene sets available online [60]. Lists of promoters bound by transcription factors at a significance of p < 0.001 [1] were collated from data available at [61].

Real-time PCR

Regions approximately 200 bp in size that span one or more Rap1 consensus sites in ribosomal protein gene promoters were amplified from ChIP and unenriched control samples using SYBR green PCR mix (Qiagen) in an MJ Research real-time PCR machine according to the manufacturers' instructions. Fold-ratios that reflect relative enrichment or depletion of a given region in the H3 or FLAG-H2B ChIP assays were determined using the 2-ΔΔCT method described in the Applied Biosystems User Bulletin. For each region examined, the TUB2 promoter was used as the normalizer (this promoter is used as a control because its occupancy approximates that of the average promoter by global analysis), and the unenriched control sample was used as the calibrator. Each reported ratio represents the average of three independent ChIP experiments analyzed in duplicate by real-time PCR. The following primer pairs were used:

RPS22A promoter: 5'-GCCTAAAACGCCCATAAGTT-3' and 5'-ACTGCAAACCCATATTCAAGA-3'

RPS15 promoter: 5'-TACACCGCGCGTATAAATCA-3' and 5'-CCCAGCAAGGAGTTTCTCAG-3'

RPS11B promoter: 5'-GAAGAAATATTTCCTTGCTGCACC-3' and 5'-AAGGGAAACGTAAAGCTATTGGAC-3'

RPL23A promoter: 5'-ATTAACATCTGTACACCCCCAACT-3' and 5'-TACAGTTCGTTTCCTGCC ATATTA-3'

TUB2 promoter: 5'-GGCCTAACAGTAAAGATATCCTCC-3' and 5'-GTTGTAGTAGCTGCTATGT CACTC-3'

Centromeric vectors containing either a mutant RPS11B promoter lacking the two Rap1 consensus motifs [30] or an essentially wild-type allele were transformed into wild-type yeast and used in an H3 ChIP assay to evaluate the consequence of removing Rap1 binding sites on nucleosome occupancy. Enrichment was evaluated by real-time PCR using the following primer pair that selectively amplifies the plasmid alleles but not the endogenous RPS11B promoter: 5'-CTGGAAGAAATATTTCCTT GCTCTAG-3' and 5'-AAGGGAAACGTAAAGCTATTGGAC-3'.

Micrococcal nuclease assay

Log-phase cultures of W303a yeast grown in 450 ml YPD to OD600 of 1.0 were spheroplasted with zymolase (10 mg in 40 ml volume of 1 M sorbitol, 50 mM Tris pH 7.4, 10 mM β-mercaptoethanol (β-ME), at 30°C for 38 min shaking at 300 rpm), divided into five aliquots, and digested with increasing concentrations (20 U to 320 U) of micrococcal nuclease (Worthington Biochem) in 600 μl 0.5 mM spermidine, 1 mM β-ME, 0.075% NP-40. DNA from digested samples was extracted with phenol twice and chloroform once and precipitated in ethanol. Samples were washed, resuspended in 10 mM Tris pH 7.5, subjected to RNAse treatment, cleaned up with the MinElute kit (Qiagen) and run out in a 1% agarose gel. Following depurination, denaturation and neutralization of the gel, DNA was transferred onto nylon membranes by capillary action and covalently linked to the membranes by UV irradiation. Southern blotting was carried out using a DIG Luminescent Detection Kit (Roche) and DIG-labeled probe generated by PCR using the TUB2, RPS11B and RPS15 primers described above.

References

  1. Lee TI, Rinaldi NJ, Robert F, Odom DT, Bar-Joseph Z, Gerber GK, Hannett NM, Harbison CT, Thompson CM, Simon I, et al: Transcriptional regulatory networks in Saccharomyces cerevisiae. Science. 2002, 298: 799-804. 10.1126/science.1075090.

    Article  PubMed  CAS  Google Scholar 

  2. Kornberg RD, Lorch Y: Twenty-five years of the nucleosome, fundamental particle of the eukaryote chromosome. Cell. 1999, 98: 285-294. 10.1016/S0092-8674(00)81958-3.

    Article  PubMed  CAS  Google Scholar 

  3. Han M, Grunstein M: Nucleosome loss activates yeast downstream promoters in vivo. Cell. 1988, 55: 1137-1145. 10.1016/0092-8674(88)90258-9.

    Article  PubMed  CAS  Google Scholar 

  4. Narlikar GJ, Fan HY, Kingston RE: Cooperation between complexes that regulate chromatin structure and transcription. Cell. 2002, 108: 475-487. 10.1016/S0092-8674(02)00654-2.

    Article  PubMed  CAS  Google Scholar 

  5. Boeger H, Griesenbeck J, Strattan JS, Kornberg RD: Nucleosomes unfold completely at a transcriptionally active promoter. Mol Cell. 2003, 11: 1587-1598. 10.1016/S1097-2765(03)00231-4.

    Article  PubMed  CAS  Google Scholar 

  6. Reinke H, Horz W: Histones are first hyperacetylated and then lose contact with the activated PHO5 promoter. Mol Cell. 2003, 11: 1599-1607. 10.1016/S1097-2765(03)00186-2.

    Article  PubMed  CAS  Google Scholar 

  7. Schreiber SL, Bernstein BE: Signaling network model of chromatin. Cell. 2002, 111: 771-778. 10.1016/S0092-8674(02)01196-0.

    Article  PubMed  CAS  Google Scholar 

  8. Hassan AH, Neely KE, Vignali M, Reese JC, Workman JL: Promoter targeting of chromatin-modifying complexes. Front Biosci. 2001, 6: D1054-D1064.

    Article  PubMed  CAS  Google Scholar 

  9. Kouzarides T: Histone methylation in transcriptional control. Curr Opin Genet Dev. 2002, 12: 198-209. 10.1016/S0959-437X(02)00287-3.

    Article  PubMed  CAS  Google Scholar 

  10. Kurdistani SK, Grunstein M: Histone acetylation and deacetylation in yeast. Nat Rev Mol Cell Biol. 2003, 4: 276-284. 10.1038/nrm1075.

    Article  PubMed  CAS  Google Scholar 

  11. Grewal SI, Elgin SC: Heterochromatin: new possibilities for the inheritance of structure. Curr Opin Genet Dev. 2002, 12: 178-187. 10.1016/S0959-437X(02)00284-8.

    Article  PubMed  CAS  Google Scholar 

  12. Turner BM: Cellular memory and the histone code. Cell. 2002, 111: 285-291. 10.1016/S0092-8674(02)01080-2.

    Article  PubMed  CAS  Google Scholar 

  13. Gregory PD, Horz W: Chromatin and transcription - how transcription factors battle with a repressive chromatin environment. Eur J Biochem. 1998, 251: 9-18. 10.1046/j.1432-1327.1998.2510009.x.

    Article  PubMed  CAS  Google Scholar 

  14. Morse RH: Getting into chromatin: how do transcription factors get past the histones?. Biochem Cell Biol. 2003, 81: 101-112. 10.1139/o03-039.

    Article  PubMed  CAS  Google Scholar 

  15. Nagy PL, Cleary ML, Brown PO, Lieb JD: Genomewide demarcation of RNA polymerase II transcription units revealed by physical fractionation of chromatin. Proc Natl Acad Sci USA. 2003, 100: 6364-6369. 10.1073/pnas.1131966100.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  16. Wang Y, Liu CL, Storey JD, Tibshirani RJ, Herschlag D, Brown PO: Precision and functional specificity in mRNA decay. Proc Natl Acad Sci USA. 2002, 99: 5860-5865. 10.1073/pnas.092538799.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Iyer V, Struhl K: Poly(dA:dT), a ubiquitous promoter element that stimulates transcription via its intrinsic DNA structure. EMBO J. 1995, 14: 2570-2579.

    PubMed  CAS  PubMed Central  Google Scholar 

  18. Anderson JD, Widom J: Poly(dA-dT) promoter elements increase the equilibrium accessibility of nucleosomal DNA target sites. Mol Cell Biol. 2001, 21: 3830-3839. 10.1128/MCB.21.11.3830-3839.2001.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Graham IR, Chambers A: Use of a selection technique to identify the diversity of binding sites for the yeast RAP1 transcription factor. Nucleic Acids Res. 1994, 22: 124-130.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  20. Buchman AR, Kimmerly WJ, Rine J, Kornberg RD: Two DNA-binding factors recognize specific sequences at silencers, upstream activating sequences, autonomously replicating sequences, and telomeres in Saccharomyces cerevisiae. Mol Cell Biol. 1988, 8: 210-225.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  21. Kellis M, Patterson N, Endrizzi M, Birren B, Lander ES: Sequencing and comparison of yeast species to identify genes and regulatory elements. Nature. 2003, 423: 241-254. 10.1038/nature01644.

    Article  PubMed  CAS  Google Scholar 

  22. Banerjee N, Zhang MQ: Identifying cooperativity among transcription factors controlling the cell cycle in yeast. Nucleic Acids Res. 2003, 31: 7024-7031. 10.1093/nar/gkg894.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  23. Miller JA, Widom J: Collaborative competition mechanism for gene activation in vivo. Mol Cell Biol. 2003, 23: 1623-1632. 10.1128/MCB.23.5.1623-1632.2003.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  24. Vashee S, Melcher K, Ding WV, Johnston SA, Kodadek T: Evidence for two modes of cooperative DNA binding in vivo that do not involve direct protein-protein interactions. Curr Biol. 1998, 8: 452-458. 10.1016/S0960-9822(98)70179-4.

    Article  PubMed  CAS  Google Scholar 

  25. Adams CC, Workman JL: Binding of disparate transcriptional activators to nucleosomal DNA is inherently cooperative. Mol Cell Biol. 1995, 15: 1405-1421.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. Gregory PD, Horz W: Mapping chromatin structure in yeast. Methods Enzymol. 1999, 304: 365-76. 10.1016/S0076-6879(99)04022-7.

    Article  PubMed  CAS  Google Scholar 

  27. Yu L, Morse RH: Chromatin opening and transactivator potentiation by RAP1 in Saccharomyces cerevisiae. Mol Cell Biol. 1999, 19: 5279-88.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. De Sanctis V, La Terra S, Bianchi A, Shore D, Burderi L, Di Mauro E, Negri R: In vivo topography of Rap1p-DNA complex at Saccharomyces cerevisiae TEF2 UAS(RPG) during transcriptional regulation. J Mol Biol. 2002, 318: 333-349. 10.1016/S0022-2836(02)00110-9.

    Article  PubMed  CAS  Google Scholar 

  29. Yu Q, Qiu R, Foland TB, Griesen D, Galloway CS, Chiu YH, Sandmeier J, Broach JR, Bi X: Rap1p and other transcriptional regulators can function in defining distinct domains of gene expression. Nucleic Acids Res. 2003, 31: 1224-1233. 10.1093/nar/gkg200.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  30. Reid JL, Iyer VR, Brown PO, Struhl K: Coordinate regulation of yeast ribosomal protein genes is associated with targeted recruitment of Esa1 histone acetylase. Mol Cell. 2000, 6: 1297-1307. 10.1016/S1097-2765(00)00128-3.

    Article  PubMed  CAS  Google Scholar 

  31. Powers T, Walter P: Regulation of ribosome biogenesis by the rapamycin-sensitive TOR-signaling pathway in Saccharomyces cerevisiae. Mol Biol Cell. 1999, 10: 987-1000.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  32. Hardwick JS, Kuruvilla FG, Tong JK, Shamji AF, Schreiber SL: Rapamycin-modulated transcription defines the subset of nutrient-sensitive signaling pathways directly controlled by the Tor proteins. Proc Natl Acad Sci USA. 1999, 96: 14866-14870. 10.1073/pnas.96.26.14866.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  33. Rohde JR, Cardenas ME: The tor pathway regulates gene expression by linking nutrient sensing to histone acetylation. Mol Cell Biol. 2003, 23: 629-635. 10.1128/MCB.23.2.629-635.2003.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Shamji AF, Kuruvilla FG, Schreiber SL: Partitioning the transcriptional program induced by rapamycin among the effectors of the Tor proteins. Curr Biol. 2000, 10: 1574-1581. 10.1016/S0960-9822(00)00866-6.

    Article  PubMed  CAS  Google Scholar 

  35. Bi X, Broach JR: UASrpg can function as a heterochromatin boundary element in yeast. Genes Dev. 1999, 13: 1089-1101.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  36. Donze D, Kamakaka RT: Braking the silence: how heterochromatic gene repression is stopped in its tracks. BioEssays. 2002, 24: 344-349. 10.1002/bies.10072.

    Article  PubMed  CAS  Google Scholar 

  37. Ai W, Bertram PG, Tsang CK, Chan TF, Zheng XF: Regulation of subtelomeric silencing during stress response. Mol Cell. 2002, 10: 1295-1305. 10.1016/S1097-2765(02)00695-0.

    Article  PubMed  CAS  Google Scholar 

  38. Braunstein M, Rose AB, Holmes SG, Allis CD, Broach JR: Transcriptional silencing in yeast is associated with reduced nucleosome acetylation. Genes Dev. 1993, 7: 592-604.

    Article  PubMed  CAS  Google Scholar 

  39. Bernstein BE, Humphrey EL, Erlich RL, Schneider R, Bouman P, Liu JS, Kouzarides T, Schreiber SL: Methylation of histone H3 Lys 4 in coding regions of active genes. Proc Natl Acad Sci USA. 2002, 99: 8695-8700. 10.1073/pnas.112318199.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  40. Ng HH, Robert F, Young RA, Struhl K: Targeted recruitment of Set1 histone methylase by elongating Pol II provides a localized mark and memory of recent transcriptional activity. Mol Cell. 2003, 11: 709-719. 10.1016/S1097-2765(03)00092-3.

    Article  PubMed  CAS  Google Scholar 

  41. Li J, Moazed D, Gygi SP: Association of the histone methyltransferase Set2 with RNA polymerase II plays a role in transcription elongation. J Biol Chem. 2002, 277: 49383-49388. 10.1074/jbc.M209294200.

    Article  PubMed  CAS  Google Scholar 

  42. Krogan NJ, Kim M, Tong A, Golshani A, Cagney G, Canadien V, Richards DP, Beattie BK, Emili A, Boone C, et al: Methylation of histone H3 by Set2 in Saccharomyces cerevisiae is linked to transcriptional elongation by RNA polymerase II. Mol Cell Biol. 2003, 23: 4207-4218. 10.1128/MCB.23.12.4207-4218.2003.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  43. Schaft D, Roguev A, Kotovic KM, Shevchenko A, Sarov M, Neugebauer KM, Stewart AF: The histone 3 lysine 36 methyltransferase, SET2, is involved in transcriptional elongation. Nucleic Acids Res. 2003, 31: 2475-2482. 10.1093/nar/gkg372.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  44. Xiao T, Hall H, Kizer KO, Shibata Y, Hall MC, Borchers CH, Strahl BD: Phosphorylation of RNA polymerase II CTD regulates H3 methylation in yeast. Genes Dev. 2003, 17: 654-663. 10.1101/gad.1055503.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  45. Li B, Howe L, Anderson S, Yates JR, Workman JL: The Set2 histone methyltransferase functions through the phosphorylated carboxyl-terminal domain of RNA polymerase II. J Biol Chem. 2003, 278: 8897-8903. 10.1074/jbc.M212134200.

    Article  PubMed  CAS  Google Scholar 

  46. Gerber M, Shilatifard A: Transcriptional elongation by RNA polymerase II and histone methylation. J Biol Chem. 2003, 278: 26303-26306. 10.1074/jbc.R300014200.

    Article  PubMed  CAS  Google Scholar 

  47. Hampsey M, Reinberg D: Tails of intrigue: phosphorylation of RNA polymerase II mediates histone methylation. Cell. 2003, 113: 429-432. 10.1016/S0092-8674(03)00360-X.

    Article  PubMed  CAS  Google Scholar 

  48. Kuo MH, Allis CD: In vivo cross-linking and immunoprecipitation for studying dynamic Protein:DNA associations in a chromatin environment. Methods. 1999, 19: 425-433. 10.1006/meth.1999.0879.

    Article  PubMed  CAS  Google Scholar 

  49. Bernstein BE, Humphrey EL, Liu CL, Schreiber SL: The use of chromatin immunoprecipitation assays in genome-wide analyses of histone modifications. Methods Enzymol. 2004, 376: 349-360.

    Article  PubMed  CAS  Google Scholar 

  50. Robzyk K, Recht J, Osley MA: Rad6-dependent ubiquitination of histone H2B in yeast. Science. 2000, 287: 501-504. 10.1126/science.287.5452.501.

    Article  PubMed  CAS  Google Scholar 

  51. Liu CL, Schreiber SL, Bernstein BE: Development and validation of a T7 based linear amplification for genomic DNA. BMC Genomics. 2003, 4: 19-10.1186/1471-2164-4-19.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  52. Carroll AS, Bishop AC, DeRisi JL, Shokat KM, O'Shea EK: Chemical inhibition of the Pho85 cyclin-dependent kinase reveals a role in the environmental stress response. Proc Natl Acad Sci USA. 2001, 98: 12578-12583. 10.1073/pnas.211195798.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  53. Iyer VR, Horak CE, Scafe CS, Botstein D, Snyder M, Brown PO: Genomic binding sites of the yeast cell-cycle transcription factors SBF and MBF. Nature. 2001, 409: 533-538. 10.1038/35054095.

    Article  PubMed  CAS  Google Scholar 

  54. Ren B, Robert F, Wyrick JJ, Aparicio O, Jennings EG, Simon I, Zeitlinger J, Schreiber J, Hannett N, Kanin E, et al: Genome-wide location and function of DNA binding proteins. Science. 2000, 290: 2306-2309. 10.1126/science.290.5500.2306.

    Article  PubMed  CAS  Google Scholar 

  55. Microarrays.org. [http://www.microarrays.org]

  56. Gollub J, Ball CA, Binkley G, Demeter J, Finkelstein DB, Hebert JM, Hernandez-Boussard T, Jin H, Kaloper M, Matese JC, et al: The Stanford Microarray Database: data access and quality assessment tools. Nucleic Acids Res. 2003, 31: 94-96. 10.1093/nar/gkg078.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  57. Stuart Schreiber group website. [http://www.broad.mit.edu/chembio/lab_schreiber/pubs/protocols/global%20nucleosome/globalnuc.html]

  58. Saccharomyces Genome Database. [http://www.yeastgenome.org]

  59. Tavazoie S, Hughes JD, Campbell MJ, Cho RJ, Church GM: Systematic determination of genetic network architecture. Nat Genet. 1999, 22: 281-285. 10.1038/10343.

    Article  PubMed  CAS  Google Scholar 

  60. Yeast comparative genomics. [http://www.broad.mit.edu/annotation/fungi/comp_yeasts]

  61. Young lab. [http://web.wi.mit.edu/young]

  62. Lieb JD, Liu X, Botstein D, Brown PO: Promoter-specific binding of Rap1 revealed by genome-wide maps of protein-DNA association. Nat Genet. 2001, 28: 327-334. 10.1038/ng569.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We thank Mary Ann Osley and Kevin Struhl for generously providing yeast strains, and Jay Bradner, Jeff McMahon, Aly Shamji, Jianping Cui, Manolis Kellis, Vamsi Mootha, Mike Kamal and Oliver Rando for insightful discussions. This study was supported by a grant from NIGMS (GM38627, awarded to S.L.S.). B.E.B. is supported by a K08 Development Award from the National Cancer Institute. C.L.L. is supported by a Graduate Research Fellowship from the National Science Foundation. S.L.S. is an Investigator at the Howard Hughes Medical Institute.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Bradley E Bernstein or Stuart L Schreiber.

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Bernstein, B.E., Liu, C.L., Humphrey, E.L. et al. Global nucleosome occupancy in yeast. Genome Biol 5, R62 (2004). https://doi.org/10.1186/gb-2004-5-9-r62

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/gb-2004-5-9-r62

Keywords