Clustered regulatory elements at nucleosome-depleted regions punctuate a constant nucleosomal landscape in Schizosaccharomyces pombe
© Soriano et al.; licensee BioMed Central Ltd. 2013
Received: 22 July 2013
Accepted: 14 November 2013
Published: 21 November 2013
Nucleosomes facilitate the packaging of the eukaryotic genome and modulate the access of regulators to DNA. A detailed description of the nucleosomal organization under different transcriptional programmes is essential to understand their contribution to genomic regulation.
To visualize the dynamics of individual nucleosomes under different transcriptional programmes we have generated high-resolution nucleosomal maps in Schizosaccharomyces pombe. We show that 98.5% of the genome remains almost invariable during mitosis and meiosis while remodelling is limited to approximately 1100 nucleosomes in the promoters of a subset of meiotic genes. These inducible nucleosome-depleted regions (NDR) and also those constitutively present in the genome overlap precisely with clusters of binding sites for transcription factors (TF) specific for meiosis and for different functional classes of genes, respectively. Deletion of two TFs affects only a small fraction of all the NDRs to which they bind in vivo, indicating that TFs collectively contribute to NDR maintenance.
Our results show that the nucleosomal profile in S. pombe is largely maintained under different physiological conditions and patterns of gene expression. This relatively constant landscape favours the concentration of regulators in constitutive and inducible NDRs. The combinatorial analysis of binding motifs in this discrete fraction of the genome will facilitate the definition of the transcriptional regulatory networks.
Nucleosomes play an essential role in the management of the eukaryotic genome by facilitating its packaging inside the nucleus. They also regulate basic genomic processes such as transcription, replication and recombination, either directly by controlling the physical access of regulators to DNA or indirectly by modulating their binding through a complex repertoire of histone modifications . In recent years, DNA microarray and high-throughput sequencing technologies have enabled the mapping of nucleosomes at genome-wide scale in many organisms, including yeasts [2–6], Drosophila, C. elegans and mammalian cells [9, 10]. In general, nucleosome profiles conform to a pattern where nucleosome-depleted regions (NDR) are present immediately upstream from the transcription start site (TSS) of most genes in the genome . From these sites, regular nucleosome arrays are generated that extend into the transcribed regions. Since most of these studies have focused on the aggregated nucleosome profiles at the 5′ end of genes and transcribed units, the distribution of nucleosomes along intergenic regions has been comparatively less studied.
Recent analyses have revealed the extremely precise organization of nucleosomes at promoters to regulate the interaction between transcription factors and DNA. For example, periodic expression of the CLN2 and HO genes during the S. cerevisiae cell cycle depends on the binding of their regulators to NDR in their promoters. Experimental manipulation of these promoters showed that when the binding sites were embedded in nucleosomes, transcription was still active but the cell cycle periodic expression was lost . By contrast, other studies have shown that transcription factors such as Rap1 and Reb1 preferentially bind to their cognate sites on the surface of the nucleosome immediately upstream from the promoter NDR. The fact that the binding sites for Rap1 face outwards on the surface of the −1 nucleosome further highlights the relevance of their positioning and of their rotational symmetry in the regulation of gene expression .
How is this precise nucleosome organization generated and maintained? Poly (dA:dT) tracts can exclude nucleosomes [14, 15] and NDRs in S. cerevisiae are enriched in these elements . Other sequences have also been selected to favour low nucleosome occupancy, such as those at the promoters of respiration genes, which are usually expressed in aerobic yeasts . However, although some sequences favor nucleosome exclusion and others can promote nucleosome occupancy , the DNA sequence alone cannot specify the nucleosome profile observed in vivo at genome-wide scale [18, 19].
The statistical positioning model  proposed that nucleosome stacking against a physical barrier would passively generate regular arrays with an internucleosomal spacing inversely proportional to nucleosome density. Recent work has shown that nucleosome positioning in S. cerevisiae is an energy-dependent process, which in the presence of ATP and a cell extract can recapitulate the in vivo profile emanating from the 5′ NDRs, even at low nucleosome density and in the absence of transcription . The model still requires the existence of physical barriers to act as organizing centers or focal points to set the beginning and confer directionality to the array [21, 22].
Nucleosome distribution in S. pombe has been comparatively much less studied than in S. cerevisiae. Several studies have generated genome-wide nucleosome maps using tiling microarrays or next generation sequencing. Some of these analyses have focused primarily on promoters and transcribed regions and have derived their conclusions from the aggregated nucleosome profiles of hundreds or thousands of genes [4, 5, 23, 24], while others have focused on the NDR profiles of replication origins and recombination hotspots [6, 25].
We have sequenced mononucleosomal DNA to a depth that allows a precise description of nucleosome dynamics at the level of individual nucleosomes across the S. pombe genome during mitosis and meiosis. We show that the great majority of the genome is organized in a very stable pattern of positioned nucleosomes and that NDRs overlap precisely with clusters of binding sites for transcription factors, which, in turn, could contribute to maintaining a regular nucleosome pattern across the genome.
Widespread nucleosome positioning in the S. pombe genome
We have generated nucleosome maps in S. pombe by sequencing mononucleosomal DNA at a genome coverage ranging from 46- to 177-fold. This sequencing depth allows the high resolution mapping of individual nucleosomes by aligning sequence reads directly onto the reference genome followed by signal smoothing, with minimal mathematical modification of the raw data (see Methods and Additional file 1: Figure S1 and Additional file 2: Figure S2 for a comparison between raw and processed data and between the sequencing and hybridization analyses).
Alignment to the midpoint position of +1 N generates sharper profiles than alignment to the TSS since +1 nucleosomes from different genes are perfectly aligned. Alignment to the TSS is expected to generate comparable nucleosomal profiles but of lower resolution due to variations in the distance between the TSS and the midpoint position of +1 N in different genes. This is reflected in the lower amplitude of the nucleosome peaks in the aggregated profiles. Alignment to TSS, however, is useful for establishing where transcription initiates relative to +1 N. Comparative alignments to +1 N and to TSS of S. pombe are shown in Additional file 3: Figure S3A.
It has been reported that in S. pombe, unlike in S. cerevisiae, nucleosome positioning only occurs downstream from the 5′ NDRs [4, 24]. However, Figure 1A shows that alignment to nucleosome +1 generated a regular nucleosomal profile upstream from the NDR. The low amplitude of the peaks (indicative of higher variation in the internucleosomal distances) could be due to the different size of individual NDRs, which would cause the −1 nucleosome (−1 N, the closest upstream from the TSS) and those upstream from it to be out of phase between different genes. This possibility was confirmed by the finding that the alignment relative to -1 N generated a clear regular profile upstream from the NDR (Figure 1B). Alignment relative to -1 N was expected to blur the periodicity of nucleosomes downstream from the NDR generated when +1 N was used as a reference, as was indeed the case (compare Figure 1A and B). To test whether positioning was maintained downstream from transcribed regions, we used the terminal nucleosome (TN, the closest upstream from the TTS) as a reference for the alignment. Figure 1D shows that the resulting profile was virtually symmetric on both sides of TN indicating that nucleosome positioning extended beyond the TTS into the adjacent intergenic region. This symmetry, however, is not present in the nucleosomal profile of S. cerevisiae, due to the presence of a NDR at the 3′ end of many genes . The comparative profiles generated by alignment to the TN or to the TTS in both yeasts is shown in Additional file 3: Figure S3C. Taken together, these results revealed a very high degree of genome-wide nucleosome positioning in S. pombe and indicated that aggregated nucleosome profiles vary depending on the reference used to align them and do not result from different properties of specific nucleosomes or gene regions.
Nucleosome remodelling at meiosis-specific promoters
To address the point of how stable this pattern was under different transcriptional regimes, we searched for differences in the genomic distribution of nucleosomes in exponential vegetative cells and at different stages through meiosis. Specifically, we compared the nucleosome profiles of asynchronous haploid and diploid cells during mitosis and of diploid cells arrested in G0 (0 h) and at 3 and 5 hours after synchronous entry into meiosis, when chromosomes are recombining (3 h), and when the four haploid nuclei have separated and start maturing into spores (5 h) (see Methods). Previous studies have shown that several hundred genes are differentially transcribed during these early, medium and late stages [27–29] (see also Additional file 5: Table S1). Despite these differences, we found that the aggregated nucleosomal profile at these stages was indistinguishable from those seen during mitosis (See Additional file 6: Figure S5 for comparison of the profiles of mitosis and meiosis at 3 hours). In order to detect possible differences at the level of specific genes, we compared the nucleosomal profiles along the entire genome and found that the major differences between meiotic and vegetative cells were limited to 782 NDRs, which were present at 0, 3 or 5 h of meiosis but were absent in mitotic cells. We defined NDRs as regions spanning at least 150 nucleotides (corresponding to the eviction of at least one nucleosome. See Methods). As an example, Figure 2A shows a 15-kb region where a nucleosome at the 5′ end of the dni2 gene in meiosis at 0 hours was absent while another occupied the NDR at the 5′ end of the divergently transcribed mcm2 and mlo2 genes in the 5-hour sample (green arrowheads). By contrast, all the other nucleosomes displayed a pattern of sharp, regularly spaced peaks that remained virtually invariable despite the different physiological stages and the chromosomal processes undergone by the cells. Since the average length of the 782 meiosis-specific NDRs was 214.68 +/− 71.06 nucleotides and the average distance between the mid-point of adjacent nucleosomes is 152 bp (Figure 1) [4, 25], we estimated that these NDRs were generated by the eviction of around 1100 nucleosomes. We have detected 78188 nucleosome peaks of nucleosome occupancy above the average genome-wide number of reads. Of these, we have considered 62824 nucleosomes (80.3%) as well positioned since their central coordinate was more than 100 nucleotides away from the nearest flanking peaks. These figures are consistent with the 77796 estimated number of nucleosomes obtained by dividing the 11825 kb of the S. pombe genome in http://www.pombase.org that remain after excluding rDNA, centromeric regions and constitutive NDRs by 152 bp. This means that nucleosome remodelling in the three meiotic samples that we analyzed relative to mitosis was limited to only 1.5% of the approximately 78000 nucleosomes in the S. pombe genome. The genomic position of the 2046 constitutive, 782 meiosis-specific and 26 mitosis-specific NDRs is indicated in Additional file 7: Table S2.
To monitor whether the meiosis-specific induced NDRs were associated with the differential transcription of the adjacent genes, we generated strand-specific transcription maps of exponential mitotic cells and of cells at 0, 3 and 5 hours during meiosis using tiling microarrays. The RNA level at these stages for all the annotated S. pombe genes (http://www.pombase.org) is indicated in Additional file 5: Table S1. Figure 2B illustrates that the eviction or repositioning of the two nucleosomes in the region shown coincided with the up- or down-regulation of the promoters of the dni2 and mcm2/mlo2 genes at 0 and 5 hours during meiosis, respectively. The nucleosomal and transcriptional profiles for the entire genome can be visualized in a searchable genome browser at http://genomics.usal.es/cgi-bin/gb2/gbrowse/Sp_nucdyn.
How widespread is the association between nucleosome remodelling and transcriptional regulation? Altogether, of the 782 meiosis-specific NDRs, 607 (77.6%) mapped to sites of transcription initiation. No transcripts were detected associated with the remaining 175 NDRs (32.4%) possibly due to a very low level of transcription, as has recently been shown for many genes in S. pombe. Alternatively, these NDRs could be generated independently of transcription. They do not seem to be dependent on replication initiation either since only 20% of them colocalize with sites of binding for the Origin Recognition Complex (ORC) . Out of the 607 NDRs that colocalized with TSS, 287 were adjacent to genes overexpressed at least 1.5-fold during meiosis relative to mitosis [27–29] (Additional file 5: Table S1).
To analyze the extent to which transcriptional activation was associated with promoter remodelling, we selected a total of 352 genes showing at least 4-fold differential expression during the following stages: 102 genes overexpressed in meiosis at 0 h relative to mitosis; 178 genes in meiosis at 3 h relative to 0 h; and 72 genes in meiosis at 5 h relative to 3 h [27–29] (Additional file 8: Table S3). Comparison of the nucleosomal profiles at these stages revealed that in 58/102 (0 h), 22/178 (3 h) and 23/72 (5 h) of the genes, differential expression was associated with the eviction of 1–2 nucleosomes from each promoter. This means that, altogether, a meiosis-specific NDR was generated in 103 out of the 352 promoters analyzed (30%) while the other 248 remained virtually identical despite the meiosis-specific overexpression. Of these 248, 118 did and 130 did not harbour constitutive NDRs. This indicates that the constant nucleosome profile present in the great majority of the genome is also maintained in approximately 70% of the promoters of meiosis-specific genes regardless of their differential activity.
Nucleosome positioning along differentially expressed genes
Clusters of binding sites for transcription factors at nucleosome-depleted regions of meiosis-specific genes
As indicated above, 103 out of the 352 genes overexpressed over 4-fold during meiosis were associated with meiosis-specific NDRs while 118 genes harboured constitutive NDRs. Would binding sites for Ste11 and Mei4 also be overrepresented in such NDRs? MEME sequence analysis did not reveal enrichment in sites for Ste11 in any of them in the 0 h, 3 h and 5 h samples. The same analysis also showed that none of the genes overexpressed at 0 h relative to mitosis had any binding motifs for Mei4. However, 50% of the genes overexpressed at 3 h relative to 0 h had 1–2 motifs and 62.5% of the genes overexpressed at 5 h relative to 3 h had 1–3 motifs for Mei4. These results are consistent with those presented in Additional file 10: Figure S7 and imply that this factor is likely to be also involved in the specific expression of the middle and late meiotic genes through the binding to those NDRs that are constitutively present at the 5′ end of some genes specifically expressed during meiosis.
Clusters of binding sites for transcription factors at nucleosome-depleted regions of cell cycle and stress response genes
Transcription factors contribute redundantly to maintaining NDRs
The other 94 (81%) NDRs that remained unaltered in the atf1 Δ and pcr1 Δ mutants were probably maintained by other transcription factors that bind to the same NDRs, consistent with the presence of multiple binding sites for TFs in NDRs of meiosis-specific and cell cycle-regulated genes (Figures 4 and 5). A possible candidate to contribute to the maintenance of these NDRs could be the Php4 CCAAT-binding factor  since analysis of the 116 NDRs bound by Atf1/Pcr1 showed that the CCAATCA sequence was present in half of the 94 NDRs that remained unaffected but was absent in the 22 that dissappeared in atf1 Δ and pcr1 Δ cells. It is also possible that Atf21 and Atf31, another two members of the ATF CREB family of transcription factors, could functionally replace Atf1 and Pcr1 in some promoters . This is supported by the fact that of the 148 genes immediately downstream from the NDRs bound by Atf1 and Pcr1 in the genome, only 18 and 14 genes were overexpressed and 11 and 16 genes were underexpressed more than 2-fold in the atf1 Δ and pcr1 Δ mutants, respectively (Additional file 12: Table S5). The transcriptional profile for the entire S. pombe genome in wild-type, atf1 Δ and pcr1 Δ mutant cells is shown in the genome browser linked to this article.
Nucleosome dynamics of the S. pombe genome
As a consequence of the high resolution of the maps generated by us here, our results regarding the nucleosomal organization in S. pombe differ in several aspects from previous studies. For example, it has been proposed that nucleosomal positioning would coincide with the length of the transcribed units and that it is absent in inactive genes, suggesting an active role for transcription in positioning . However, we observed that nucleosomal positioning extended beyond the transcription termination site (Figure 1D) and that it was also present at genes showing low or undetectable transcription levels (Figure 3A, B and Additional file 9: Figure S6). Also, although previous analyses have failed to detect positioning upstream from the NDR at the 5′ position of the genes [4, 24, 25], we found that nucleosomal arrays emanate bidirectionally from this NDR (Figure 1B) in a fashion comparable to that in S. cerevisiae[2, 3, 26]. A comparison between the nucleosomal profiles of S. pombe and S. cerevisiae using the same criteria shows that profiles downstream from +1 N were similar (Additional file 3: Figure S3A). However, nucleosome positioning upstream from the NDR was slightly higher in S. cerevisiae (as indicated by the height of the peaks). To test whether this effect might be partially due to the fact that the size of NDRs is more homogeneous in S. cerevisiae[16, 26] than in S. pombe, we also aligned the nucleosomal profiles to -1 N and obtained comparable results in both yeasts (Additional file 3: Figure S3B). In the two cases, positioning was stricter downstream (+1 N alignment) than upstream (−1 N alignment) from the NDR, probably due to the variable distance of other upstream NDRs and genes that set the nucleosomal arrays out of phase in the aggregated profiles.
It is possible that the regular nucleosome pattern upstream from the NDR in S. pombe was not detected earlier owing to the limited resolution of the microarray analyses and the lower sequencing depth of previous studies, which made it difficult to define the position of the −1 nucleosome [4, 5, 24, 25].
Another contribution of our work is the possibility of a genome-wide quantitative assessment of nucleosome dynamics. We estimate that only 1.5% of the approximately 78000 nucleosomes in the genome are remodelled during the expression of the meiotic transcriptional programme. Remodelling, however, is not always associated with transcriptional activation, as illustrated by the invariable nucleosome profile of 70% of the promoters of the 352 genes overexpressed more than 4-fold at specific stages of meiosis relative to mitosis (Additional file 8: Table S3). A lack of correlation between transcriptional activity and the presence of NDRs at promoters has also been observed in S. cerevisiae.
Specification of genomic nucleosomal patterns
In S. cerevisiae, poly (dA:dT) elements are overrepresented in the NDRs associated with promoters [14, 15]. By contrast, these elements are not particularly enriched in S. pombe and, in agreement with previous data , we found that only 8.3% and 3.4% of the poly (dA:dT) elements 15 and 7 nucleotides long, respectively, colocalized with constitutive NDRs. Conversely, only 13% of these included poly (dA:dT) tracts 7 or 15 nucleotides long (data not shown). These results suggest that the contribution of these elements to the generation of NDRs is much smaller than in S. cerevisiae[14–16].
The statistical positioning model proposes the existence of physical barriers from which regular nucleosomal arrays are generated [20, 21] and transcriptional regulatory complexes are obvious candidates for this role. These large complexes include transcription factors, coactivators, histone modifiers and chromatin remodellers  and are targeted to promoters by transcription factors, which are their only components capable of recognizing specific sequences on DNA. Binding motifs for transcription factors are made up of degenerated sequences 6–10 bp long scattered across the genome, although only a small fraction of them are actually bound in vivo. For example, in S. cerevisiae the Leu3 transcription factor binds only a subset of all potential binding sites in the chromosomes even though the protein binds to all of them with comparable affinity on naked DNA. Sites bound in vivo were strongly correlated with a low nucleosome occupancy, which was not dependent on the presence of the Leu3 protein, suggesting an opportunistic use of accessible chromatin sites already available in the genome . In S. pombe, 148 (6.9%) out of the 2141 TGACGT hexamers representing potential binding sites for Atf1/Pcr1 are detectably bound in the genome . Of these 148 sites, 116 (78.4%) colocalized with NDRs, of which only 22 disappeared in the absence of the factors (Figure 6), pointing to the presence of redundant elements in their maintenance, as in the case of Leu3. This possibility is strongly supported by the recent finding that up to 8 proteins are involved in the maintenance of the NDR at the CLN2 promoter in S. cerevisiae. Along the same lines, NDRs at heterochromatic regions in S. pombe depend on the combinatorial contribution of at least 5 different proteins . However, NDRs can also be generated as a consequence of the binding of a single transcription factor, as in the case of the M26 mutation in S. pombe, which generates a binding site for Atf1 . Also, DNA sequences refractory to bending such as poly (dA:dT) [14, 15], poly G tracts  and sequences at some promoters  favour the exclusion of nucleosomes.
Previous studies have found that the binding sites for transcription factors tend to cluster at promoters [3, 36, 56, 57] and our results show that these clusters overlap precisely with NDRs at the promoters of meiosis-specific (Figure 4), cell cycle-regulated genes, and stress-response genes (Figure 5). The complexity of the clusters of binding sites for TFs suggests that each NDR and each promoter is probably unique, making it very difficult to predict the effect that the removal of specific elements will have in its maintenance or on the expression of the adjacent genes. For example, while the lack of Leu3 in S. cerevisiae and of Atf1/Pcr1 in S. pombe (Figure 6) has a limited impact on NDRs genome-wide, deletion of the Abf1 and Reb1 transcription factors in S. cerevisiae negatively affects a larger number of NDRs [33, 57]. Similarly, the Sap1 protein, which is involved in DNA replication and mating-type switching and binds preferentially to NDRs, is required for the maintenance of a large number of them in S. pombe.
This collaborative strategy of TFs to assemble transcription complexes at promoters has several immediate advantages for genome regulation. The first is that it favours transcription initiation from the 5′ end of genes and reduces spurious initiation at single binding sites dispersed along the genome that are usually occluded by nucleosomes. Second, the redundant contribution of several factors guarantees the maintenance of NDRs, regardless of fluctuations in their concentrations. Altogether, we have identified 2046 constitutive NDRs in mitotic and meiotic cells and 782 meiosis-specific NDRs. The constitutive presence of a large fraction of NDRs could facilitate the rapid transcriptional induction of, for example, stress response genes , and could also contribute to other processes such as the specification of meiotic recombination hotspots [6, 58]. A third function of NDRs and the complexes bound to them would be to act as barriers  or organizing centers from where nucleosome arrays are generated [16, 21, 22]. The targeting of these complexes by transcription factors to the same specific sites in every genome and the relatively close proximity between promoters could explain why a virtually identical nucleosomal arrangement is maintained in all the cells in the population under very different physiological conditions. A stable nucleosomal pattern during mitosis and meiosis has also been observed in S. cerevisiae. This is compatible with an active turnover of nucleosomes as has been described in S. cerevisiae and with variations in the epigenetic modifications of histones  that could provide flexible regulatory signalling while maintaining a constant nucleosomal framework.
Nucleosome dynamics in S. pombe genome is limited to a very small fraction of the genome that overlaps with regulatory regions. We have shown that NDRs encompass clusters of TFBS specific for different gene functions such as meiosis, cell-cycle regulation and stress response. Sequence analysis of these narrow and well-defined regions should help to define genomic regulatory networks based on the combinatorial and collective contribution of regulatory elements shared by different promoters. On the other hand, the high degree of order at nucleosome level probably underlies the ordered structure of chromatin organization at higher levels, evidenced by the specific pattern of interactions between different chromosomes  and by the global three-dimensional architecture of the nucleus, where each chromosome occupies a specific territory [63, 64]. Despite the many instances in which nucleosome remodelling at promoters is not associated with changes in transcriptional activity , the maintenance of a highly organized nucleosomal pattern is likely to be important for genomic stability, as illustrated by the gross alterations in recombination , and for preventing cryptic and unscheduled antisense transcription that results when the level of histones or the positioning of nucleosomes is altered [66–68].
S. pombe strains, growth conditions and meiosis synchronization
Wild-type (972 h-), atf1Δ (h + ura4.d18 leu1.32 ade6M210 atf1Δ::kanMX4) and pcr1Δ (h + ura4.d18 leu1.32 ade6M210 pcr1Δ::kanMX4) cells were grown in rich medium (YES) at 32°C up to a A595 = 0.8. Diploid pat1.114 (h-/h- pat1.114/pat1.114 leu1.32/leu1.32 ade6M210/ade6M216) asynchronous cells were grown in Minimal Medium (MM) supplemented with 0,1 g/l leucine at 25°C up to a A595 = 0.8. Synchronous meiosis of diploid pat1.114 cells was induced as described .
Preparation of mononucleosomal DNA
Mononucleosomal DNA for microarray hybridization and sequencing was isolated as described . The amount of Zymolyase 20 T used to prepare spheroplasts was optimized experimentally for each S. pombe strain and for the different physiological conditions to generate a 80:20 ratio of mononucleosomes to dinucleosomes, as described . The following amounts of Zymolyase were added to the cell suspension in 10 ml of sorbitol–Tris buffer: Meiosis 0 h (10 mg), meiosis 3 h (40 mg), meiosis 5 h (65 mg), mitosis (40 mg), atf1 Δ (8 mg), pcr1 Δ (15 mg).
Microarray analysis of transcription
For quantitative transcriptional analyses, Affymetrix GeneChip S. pombe 1.0FR tiling microarrays were used. Target labelling preserving the original polarity of RNAs was performed following the instructions of the GeneChip whole-transcript sense target-labelling assay manual from Affymetrix. Quantitative measurement of strand-specific differential transcription in Additional file 5: Table S1, Additional file 7: Table S2, Additional file 8: Table S3, Additional file 11: Table S4 and Additional file 12: Table S5 are indicated as the log2 value of the averaged hybridization signal from the probes spanning every ORF in the microarray. A detailed description of the method has been previously described . Normalized raw microarray signals without smoothing or denoising are shown in Figures 2, 3, 4, 5, 6 and Additional file 9: Figure S6 (green vertical lines). Hybridization signals from probes mapping to more than one position in the genome were normalized relative to the number of repeats. For all the genomic analyses, we used the S. pombe genome version of 23/08/07 and the annotation of 24/02/11 in (http://www.pombase.org) as a reference.
Affymetrix GeneChip S. pombe 1.0FR tiling microarrays were hybridized with mononucleosomal DNA from wild type, atf1 Δ and pcr1 Δ S. pombe strains (Figure 6) following the instructions of the Affymetrix GeneChip whole-transcript double-stranded target-labelling assay manual.
Sequencing analysis of mononucleosomal DNA and nucleosome depleted regions
Mononucleosomal DNA isolated as described above from haploid and diploid S. pombe strains and from cells at 0, 3 and 5 h into meiosis was sequenced in an Illumina Genome Analyzer IIx. 16588557 to 35703552 single reads 36 or 40 nucleotides long, depending on the sample, representing an average genome coverage ranging from 46- to 177-fold were aligned to the S. pombe reference genome described above. The alignment of reads generated two peaks (one on each strand) corresponding to the boundaries of each nucleosome. We used the smoothed signal generated by using the multilevel 1-D biorthogonal wavelet decomposition/reconstruction tool implemented in the Matlab “Wavelet Toolbox” to calculate the average spacing between boundary peaks for individual nucleosomes.
This parameter (which was estimated for every independent experiment) defined the distance that the individual profiles from each DNA strand had to be shifted to converge and define the midpoint position of each nucleosome. The resulting combined profile was wavelet-smoothed to generate the final nucleosome positioning profile. Comparison of the raw and wavelet-processed data in Additional file 1: Figure S1 shows that the mathematical modification of the raw sequence data was minimal. The wavelet-smoothed signal facilitated the straightforward detection of specific nucleosomes (−1, +1, central and terminal) relative to the transcription initiation or termination sites. Alignment to the +1 nucleosome (Figure 1) has been used in previous studies [68, 69] and generates sharper profiles than alignment to the transcription start site (TSS) (Additional file 3: Figure S3A).
The sequencing coverage for every nucleotide was divided by the average genomic coverage to normalize the different experiments. The Pearson Correlation Coefficient between all possible combinations of the complete raw sequencing datasets of mitosis and meiosis at 0, 3 and 5 hours, ranged between 0.74 and 0.88. We defined nucleosome-depleted regions (NDRs) as regions spanning at least 150 nucleotides (corresponding to the eviction of at least one nucleosome) with a normalized sequence coverage lower than 0.4. The resulting population of NDRs was largely coincident with those detected previously using tiling microarrays .
Sequence motif analysis
To identify DNA motifs, sequences from specific groups of NDRs were extracted and used as input for the MEME (Multiple EM for Motif Elicitation) algorithm  using the same settings as described , except that we used as a background set a fifth-order Markov model representing possible nucleotide pentuplets in all S. pombe 5′ NDRs.
The genomic data described in this work can be accessed from a searchable genome browser at http://genomics.usal.es/cgi-bin/gb2/gbrowse/Sp_nucdyn. All sequencing and microarray data are deposited in the Gene Expression Omnibus (GEO) database under the accession number GSE41773.
We are grateful to Rodrigo Bermejo for insightful comments on the manuscript. We also thank Encarnación Fermiñán and Orlando Domínguez for excellent technical support with microarray hybridization and high throughput sequencing, respectively. IS was supported by a postgraduate fellowship from the Spanish Ministerio de Economía y Competitividad (MINECO). This work was funded by grants BFU2011-28804 and Consolider-Ingenio CSD2007-00015 (from MINECO).
- Luger K, Dechassa ML, Tremethick DJ: New insights into nucleosome and chromatin structure: an ordered state or a disordered affair?. Nat Rev Mol Cell Biol. 2012, 13: 436-447.PubMed CentralView ArticlePubMedGoogle Scholar
- Yuan GC, Liu YJ, Dion MF, Slack MD, Wu LF, Altschuler SJ, Rando OJ: Genome-scale identification of nucleosome positions in S. cerevisiae. Science. 2005, 309: 626-630.View ArticlePubMedGoogle Scholar
- Lee W, Tillo D, Bray N, Morse RH, Davis RW, Hughes TR, Nislow C: A high-resolution atlas of nucleosome occupancy in yeast. Nat Genet. 2007, 39: 1235-1244.View ArticlePubMedGoogle Scholar
- Lantermann AB, Straub T, Stralfors A, Yuan GC, Ekwall K, Korber P: Schizosaccharomyces pombe genome-wide nucleosome mapping reveals positioning mechanisms distinct from those of Saccharomyces cerevisiae. Nat Struct Mol Biol. 2010, 17: 251-257.View ArticlePubMedGoogle Scholar
- Tsankov A, Yanagisawa Y, Rhind N, Regev A, Rando OJ: Evolutionary divergence of intrinsic and trans-regulated nucleosome positioning sequences reveals plastic rules for chromatin organization. Genome Res. 2011, 21: 1851-1862.PubMed CentralView ArticlePubMedGoogle Scholar
- de Castro E, Soriano I, Marin L, Serrano R, Quintales L, Antequera F: Nucleosomal organization of replication origins and meiotic recombination hotspots in fission yeast. EMBO J. 2012, 31: 124-137.PubMed CentralView ArticlePubMedGoogle Scholar
- Mavrich TN, Jiang C, Ioshikhes IP, Li X, Venters BJ, Zanton SJ, Tomsho LP, Qi J, Glaser RL, Schuster SC, Gilmour DS, Albert I, Pugh BF: Nucleosome organization in the Drosophila genome. Nature. 2008, 453: 358-362.PubMed CentralView ArticlePubMedGoogle Scholar
- Valouev A, Ichikawa J, Tonthat T, Stuart J, Ranade S, Peckham H, Zeng K, Malek JA, Costa G, McKernan K, Sidow A, Fire A, Johnson SM: A high-resolution, nucleosome position map of C. elegans reveals a lack of universal sequence-dictated positioning. Genome Res. 2008, 18: 1051-1063.PubMed CentralView ArticlePubMedGoogle Scholar
- Schones DE, Cui K, Cuddapah S, Roh TY, Barski A, Wang Z, Wei G, Zhao K: Dynamic regulation of nucleosome positioning in the human genome. Cell. 2008, 132: 887-898.View ArticlePubMedGoogle Scholar
- Valouev A, Johnson SM, Boyd SD, Smith CL, Fire AZ, Sidow A: Determinants of nucleosome organization in primary human cells. Nature. 2011, 474: 516-520.PubMed CentralView ArticlePubMedGoogle Scholar
- Whitehouse I, Tsukiyama T: Opening windows to the genome. Cell. 2009, 137: 400-402.View ArticlePubMedGoogle Scholar
- Bai L, Charvin G, Siggia ED, Cross FR: Nucleosome-depleted regions in cell-cycle-regulated promoters ensure reliable gene expression in every cell cycle. Dev Cell. 2010, 18: 544-555.PubMed CentralView ArticlePubMedGoogle Scholar
- Koerber RT, Rhee HS, Jiang C, Pugh BF: Interaction of transcriptional regulators with specific nucleosomes across the Saccharomyces genome. Mol Cell. 2009, 35: 889-902.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 CentralPubMedGoogle Scholar
- Segal E, Widom J: Poly(dA:dT) tracts: major determinants of nucleosome organization. Curr Opin Struct Biol. 2009, 19: 65-71.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Z, Wippo CJ, Wal M, Ward E, Korber P, Pugh BF: A packing mechanism for nucleosome organization reconstituted across a eukaryotic genome. Science. 2011, 332: 977-980.View ArticlePubMedGoogle Scholar
- Field Y, Fondufe-Mittendorf Y, Moore IK, Mieczkowski P, Kaplan N, Lubling Y, Lieb JD, Widom J, Segal E: Gene expression divergence in yeast is coupled to evolution of DNA-encoded nucleosome organization. Nat Genet. 2009, 41: 438-445.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Y, Moqtaderi Z, Rattner BP, Euskirchen G, Snyder M, Kadonaga JT, Liu XS, Struhl K: Intrinsic histone-DNA interactions are not the major determinant of nucleosome positions in vivo. Nat Struct Mol Biol. 2009, 16: 847-852.PubMed CentralView ArticlePubMedGoogle Scholar
- Struhl K, Segal E: Determinants of nucleosome positioning. Nat Struct Mol Biol. 2013, 20: 267-273.PubMed CentralView ArticlePubMedGoogle Scholar
- Kornberg RD, Stryer L: Statistical distributions of nucleosomes: nonrandom locations by a stochastic mechanism. Nucleic Acids Res. 1988, 16: 6677-6690.PubMed CentralView ArticlePubMedGoogle Scholar
- Mavrich TN, Ioshikhes IP, Venters BJ, Jiang C, Tomsho LP, Qi J, Schuster SC, Albert I, Pugh BF: A barrier nucleosome model for statistical positioning of nucleosomes throughout the yeast genome. Genome Res. 2008, 18: 1073-1083.PubMed CentralView ArticlePubMedGoogle Scholar
- Yen K, Vinayachandran V, Batta K, Koerber RT, Pugh BF: Genome-wide nucleosome specificity and directionality of chromatin remodelers. Cell. 2012, 149: 1461-1473.PubMed CentralView ArticlePubMedGoogle Scholar
- Kristell C, Orzechowski Westholm J, Olsson I, Ronne H, Komorowski J, Bjerling P: Nitrogen depletion in the fission yeast Schizosaccharomyces pombe causes nucleosome loss in both promoters and coding regions of activated genes. Genome Res. 2010, 20: 361-371.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu J, Yanagisawa Y, Tsankov AM, Hart C, Aoki K, Kommajosyula N, Steinmann KE, Bochicchio J, Russ C, Regev A, Rando OJ, Nusbaum C, Niki H, Milos P, Weng Z, Rhind N: Genome-wide identification and characterization of replication origins by deep sequencing. Genome Biol. 2012, 13: R27-PubMed CentralView ArticlePubMedGoogle Scholar
- Givens RM, Lai WK, Rizzo JM, Bard JE, Mieczkowski PA, Leatherwood J, Huberman JA, Buck MJ: Chromatin architectures at fission yeast transcriptional promoters and replication origins. Nucleic Acids Res. 2012, 40: 7176-7189.PubMed CentralView ArticlePubMedGoogle Scholar
- Shivaswamy S, Bhinge A, Zhao Y, Jones S, Hirst M, Iyer VR: Dynamic remodeling of individual nucleosomes across a eukaryotic genome in response to transcriptional perturbation. PLoS Biol. 2008, 6: e65-PubMed CentralView ArticlePubMedGoogle Scholar
- Mata J, Lyne R, Burns G, Bähler J: The transcriptional program of meiosis and sporulation in fission yeast. Nat Genet. 2002, 32: 143-147.View ArticlePubMedGoogle Scholar
- Mata J, Bähler J: Global roles of Ste11p, cell type, and pheromone in the control of gene expression during early sexual differentiation in fission yeast. Proc Natl Acad Sci USA. 2006, 103: 15517-15522.PubMed CentralView ArticlePubMedGoogle Scholar
- Bitton DA, Grallert A, Scutt PJ, Yates T, Li Y, Bradford JR, Hey Y, Pepper SD, Hagan IM, Miller CJ: Programmed fluctuations in sense/antisense transcript ratios drive sexual differentiation in S. pombe. Mol Syst Biol. 2011, 7: 559-PubMed CentralView ArticlePubMedGoogle Scholar
- Marguerat S, Schmidt A, Codlin S, Chen W, Aebersold R, Bähler J: Quantitative analysis of fission yeast transcriptomes and proteomes in proliferating and quiescent cells. Cell. 2012, 151: 671-683.PubMed CentralView ArticlePubMedGoogle Scholar
- Belotserkovskaya R, Oh S, Bondarenko VA, Orphanides G, Studitsky VM, Reinberg D: FACT facilitates transcription-dependent nucleosome alteration. Science. 2003, 301: 1090-1093.View ArticlePubMedGoogle Scholar
- Wilhelm BT, Marguerat S, Aligianni S, Codlin S, Watt S, Bahler J: Differential patterns of intronic and exonic DNA regions with respect to RNA polymerase II occupancy, nucleosome density and H3K36me3 marking in fission yeast. Genome Biol. 2011, 12: R82-PubMed CentralView ArticlePubMedGoogle Scholar
- Hartley PD, Madhani HD: Mechanisms that specify promoter nucleosome location and identity. Cell. 2009, 137: 445-458.PubMed CentralView ArticlePubMedGoogle Scholar
- Charoensawan V, Janga SC, Bulyk ML, Babu MM, Teichmann SA: DNA sequence preferences of transcriptional activators correlate more strongly than repressors with nucleosomes. Mol Cell. 2012, 47: 183-192.PubMed CentralView ArticlePubMedGoogle Scholar
- Bailey TL, Elkan C: Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol. 1994, 2: 28-36.PubMedGoogle Scholar
- Oliva A, Rosebrock A, Ferrezuelo F, Pyne S, Chen H, Skiena S, Futcher B, Leatherwood J: The cell cycle-regulated genes of Schizosaccharomyces pombe. PLoS Biol. 2005, 3: e225-PubMed CentralView ArticlePubMedGoogle Scholar
- van Beest M, Dooijes D, van De Wetering M, Kjaerulff S, Bonvin A, Nielsen O, Clevers H: Sequence-specific high mobility group box factors recognize 10-12-base pair minor groove motifs. J Biol Chem. 2000, 275: 27266-27273.PubMedGoogle Scholar
- Abe H, Shimoda C: Autoregulated expression of Schizosaccharomyces pombe meiosis-specific transcription factor Mei4 and a genome-wide search for its target genes. Genetics. 2000, 154: 1497-1508.PubMed CentralPubMedGoogle Scholar
- Mata J, Wilbrey A, Bähler J: Transcriptional regulatory network for sexual differentiation in fission yeast. Genome Biol. 2007, 8: R217-PubMed CentralView ArticlePubMedGoogle Scholar
- Murakami-Tonami Y, Yamada-Namikawa C, Tochigi A, Hasegawa N, Kojima H, Kunimatsu M, Nakanishi M, Murakami H: Mei4p coordinates the onset of meiosis I by regulating cdc25+ in fission yeast. Proc Natl Acad Sci USA. 2007, 104: 14688-14693.PubMed CentralView ArticlePubMedGoogle Scholar
- Hao Z, Furunobu A, Nagata A, Okayama H: A zinc finger protein required for stationary phase viability in fission yeast. J Cell Sci. 1997, 110: 2557-2566.PubMedGoogle Scholar
- Horie S, Watanabe Y, Tanaka K, Nishiwaki S, Fujioka H, Abe H, Yamamoto M, Shimoda C: The Schizosaccharomyces pombe mei4+ gene encodes a meiosis-specific transcription factor containing a forkhead DNA-binding domain. Mol Cell Biol. 1998, 18: 2118-2129.PubMed CentralView ArticlePubMedGoogle Scholar
- Quintales L, Sanchez M, Antequera F: Analysis of DNA strand-specific differential expression with high density tiling microarrays. BMC Bioinformatics. 2010, 11: 136-PubMed CentralView ArticlePubMedGoogle Scholar
- Sansó M, Vargas-Perez I, Quintales L, Antequera F, Ayte J, Hidalgo E: Gcn5 facilitates Pol II progression, rather than recruitment to nucleosome-depleted stress promoters, in Schizosaccharomyces pombe. Nucleic Acids Res. 2011, 39: 6369-6379.PubMed CentralView ArticlePubMedGoogle Scholar
- Benbrook DM, Jones NC: Different binding specificities and transactivation of variant CRE’s by CREB complexes. Nucleic Acids Res. 1994, 22: 1463-1469.PubMed CentralView ArticlePubMedGoogle Scholar
- Kanoh J, Watanabe Y, Ohsugi M, Iino Y, Yamamoto M: Schizosaccharomyces pombe gad7+ encodes a phosphoprotein with a bZIP domain, which is required for proper G1 arrest and gene expression under nitrogen starvation. Genes Cells. 1996, 1: 391-408.View ArticlePubMedGoogle Scholar
- Watanabe Y, Yamamoto M: Schizosaccharomyces pombe pcr1+ encodes a CREB/ATF protein involved in regulation of gene expression for sexual development. Mol Cell Biol. 1996, 16: 704-711.PubMed CentralView ArticlePubMedGoogle Scholar
- Morita T, Yamada T, Yamada S, Matsumoto K, Ohta K: Fission yeast ATF/CREB family protein Atf21 plays important roles in production of normal spores. Genes Cells. 2011, 16: 217-230.View ArticlePubMedGoogle Scholar
- Eshaghi M, Lee JH, Zhu L, Poon SY, Li J, Cho KH, Chu Z, Karuturi RK, Liu J: Genomic binding profiling of the fission yeast stress-activated MAPK Sty1 and the bZIP transcriptional activator Atf1 in response to H2O2. PLoS One. 2010, 5: e11620-PubMed CentralView ArticlePubMedGoogle Scholar
- Mercier A, Watt S, Bahler J, Labbe S: Key function for the CCAAT-binding factor Php4 to regulate gene expression in response to iron deficiency in fission yeast. Eukaryot Cell. 2008, 7: 493-508.PubMed CentralView ArticlePubMedGoogle Scholar
- Weake VM, Workman JL: Inducible gene expression: diverse regulatory mechanisms. Nat Rev Genet. 2010, 11: 426-437.View ArticlePubMedGoogle Scholar
- Liu X, Lee CK, Granek JA, Clarke ND, Lieb JD: Whole-genome comparison of Leu3 binding in vitro and in vivo reveals the importance of nucleosome occupancy in target site selection. Genome Res. 2006, 16: 1517-1528.PubMed CentralView ArticlePubMedGoogle Scholar
- Bai L, Ondracka A, Cross FR: Multiple sequence-specific factors generate the nucleosome-depleted region on CLN2 promoter. Mol Cell. 2011, 42: 465-476.PubMed CentralView ArticlePubMedGoogle Scholar
- Garcia JF, Dumesic PA, Hartley PD, El-Samad H, Madhani HD: Combinatorial, site-specific requirement for heterochromatic silencing factors in the elimination of nucleosome-free regions. Genes Dev. 2010, 24: 1758-1771.PubMed CentralView ArticlePubMedGoogle Scholar
- Mizuno K, Emura Y, Baur M, Kohli J, Ohta K, Shibata T: The meiotic recombination hot spot created by the single-base substitution ade6-M26 results in remodeling of chromatin structure in fission yeast. Genes Dev. 1997, 11: 876-886.View ArticlePubMedGoogle Scholar
- Harbison CT, Gordon DB, Lee TI, Rinaldi NJ, Macisaac KD, Danford TW, Hannett NM, Tagne JB, Reynolds DB, Yoo J, Jennings EG, Zeitlinger J, Pokholok DK, Kellis M, Rolfe PA, Takusagawa KT, Lander ES, Gifford DK, Fraenkel E, Young RA: Transcriptional regulatory code of a eukaryotic genome. Nature. 2004, 431: 99-104.PubMed CentralView ArticlePubMedGoogle Scholar
- Badis G, Chan ET, van Bakel H, Pena-Castillo L, Tillo D, Tsui K, Carlson CD, Gossett AJ, Hasinoff MJ, Warren CL, Gebbia M, Talukder S, Yang A, Mnaimneh S, Terterov D, Coburn D, Li Yeo A, Yeo ZX, Clarke ND, Lieb JD, Ansari AZ, Nislow C, Hughes TR: A library of yeast transcription factor motifs reveals a widespread function for Rsc3 in targeting nucleosome exclusion at promoters. Mol Cell. 2008, 32: 878-887.PubMed CentralView ArticlePubMedGoogle Scholar
- Pan J, Sasaki M, Kniewel R, Murakami H, Blitzblau HG, Tischfield SE, Zhu X, Neale MJ, Jasin M, Socci ND, Hochwagen A, Keeney S: A hierarchical combination of factors shapes the genome-wide topography of yeast meiotic recombination initiation. Cell. 2011, 144: 719-731.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang L, Ma H, Pugh BF: Stable and dynamic nucleosome states during a meiotic developmental process. Genome Res. 2011, 122: 875-884.View ArticleGoogle Scholar
- Dion MF, Kaplan T, Kim M, Buratowski S, Friedman N, Rando OJ: Dynamics of replication-independent histone turnover in budding yeast. Science. 2007, 315: 1405-1408.View ArticlePubMedGoogle Scholar
- Cam HP, Sugiyama T, Chen ES, Chen X, FitzGerald PC, Grewal SI: Comprehensive analysis of heterochromatin- and RNAi-mediated epigenetic control of the fission yeast genome. Nat Genet. 2005, 37: 809-819.View ArticlePubMedGoogle Scholar
- Kalhor R, Tjong H, Jayathilaka N, Alber F, Chen L: Genome architectures revealed by tethered chromosome conformation capture and population-based modeling. Nat Biotechnol. 2012, 30: 90-98.View ArticleGoogle Scholar
- Misteli T: Beyond the sequence: cellular organization of genome function. Cell. 2007, 128: 787-800.View ArticlePubMedGoogle Scholar
- Rajapakse I, Groudine M: On emerging nuclear order. J Cell Biol. 2011, 192: 711-721.PubMed CentralView ArticlePubMedGoogle Scholar
- Prado F, Aguilera A: Partial depletion of histone H4 increases homologous recombination-mediated genetic instability. Mol Cell Biol. 2005, 25: 1526-1536.PubMed CentralView ArticlePubMedGoogle Scholar
- Gossett AJ, Lieb JD: In vivo effects of histone H3 depletion on nucleosome occupancy and position in Saccharomyces cerevisiae. PLoS Genet. 2012, 8: e1002771-PubMed CentralView ArticlePubMedGoogle Scholar
- Shim YS, Choi Y, Kang K, Cho K, Oh S, Lee J, Grewal SI, Lee D: Hrp3 controls nucleosome positioning to suppress non-coding transcription in eu- and heterochromatin. EMBO J. 2012, 31: 4375-4387.PubMed CentralView ArticlePubMedGoogle Scholar
- Hennig BP, Bendrin K, Zhou Y, Fischer T: Chd1 chromatin remodelers maintain nucleosome organization and repress cryptic transcription. EMBO Rep. 2012, 13: 997-1003.PubMed CentralView ArticlePubMedGoogle Scholar
- Vaillant C, Palmeira L, Chevereau G, Audit B, d’Aubenton-Carafa Y, Thermes C, Arneodo A: A novel strategy of transcription regulation by intragenic nucleosome ordering. Genome Res. 2010, 20: 59-67.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsui K, Dubuis S, Gebbia M, Morse RH, Barkai N, Tirosh I, Nislow C: Evolution of nucleosome occupancy: conservation of global properties and divergence of gene-specific patterns. Mol Cell Biol. 2011, 31: 4348-4355.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.