- Research article
- Open Access
Different nucleosomal architectures at early and late replicating origins in Saccharomyces cerevisiae
- Ignacio Soriano†1,
- Esther C Morafraile†1,
- Enrique Vázquez1,
- Francisco Antequera1 and
- Mónica Segurado1, 2Email author
© Soriano et al.; licensee BioMed Central Ltd. 2014
- Received: 14 April 2014
- Accepted: 10 September 2014
- Published: 13 September 2014
Eukaryotic genomes are replicated during S phase according to a temporal program. Several determinants control the timing of origin firing, including the chromatin environment and epigenetic modifications. However, how chromatin structure influences the timing of the activation of specific origins is still poorly understood.
By performing high-resolution analysis of genome-wide nucleosome positioning we have identified different chromatin architectures at early and late replication origins. These different patterns are already established in G1 and are tightly correlated with the organization of adjacent transcription units. Moreover, specific early and late nucleosomal patterns are fixed robustly, even in rpd3 mutants in which histone acetylation and origin timing have been significantly altered. Nevertheless, higher histone acetylation levels correlate with the local modulation of chromatin structure, leading to increased origin accessibility. In addition, we conducted parallel analyses of replication and nucleosome dynamics that revealed that chromatin structure at origins is modulated during origin activation.
Our results show that early and late replication origins present distinctive nucleosomal configurations, which are preferentially associated to different genomic regions. Our data also reveal that origin structure is dynamic and can be locally modulated by histone deacetylation, as well as by origin activation. These data offer novel insight into the contribution of chromatin structure to origin selection and firing in budding yeast.
- Histone Acetylation
- Replication Origin
- Replication Timing
- Early Origin
- Late Origin
Eukaryotic chromosomes are replicated from multiple replication origins. In budding yeast, replication origins are termed ARS (autonomously replicating sequences) and contain the ARS consensus sequence (ACS), which is essential for origin activity . During the G1 phase, pre-replicative complexes (pre-RCs), composed of the origin recognition complex (ORC), Cdc6, Cdt1, and an inactive form of the replicative helicase Mcm2-7 complex are assembled into ACS-containing origin regions, a process called ‘licensing’ . During S phase, origin activation is induced by the activation of the MCM2-7 helicase and by the recruitment of further initiation factors, including Mcm10, Cdc45, GINS, Sld2, Sld3 and Dpb11, to form pre-initiation complexes (pre-IC) . Although all replication origins are ‘licensed’ during G1 phase [4, 5], not all of them are activated during the subsequent replication cycle and, moreover, origins are activated at different times along S phase. The temporal order of origin activation is called the replication timing program, and it is established in G1 phase .
The mechanisms orchestrating the timing program are not fully understood, but many studies have indicated that the firing of origins can be influenced by different factors involving replication initiation mechanisms, epigenetic modifications, chromosomal position and the chromatin context [7–17].
Limiting levels of initiation factors have been reported [15, 16], and constitute a point of control of replication initiation. The amount of Sld3, Sld2, Sld7, Dpb11, Cdc45 and Dbf4 factors is below the number of “licensed” origins, and hence activation is necessarily restricted to a subset of origins. Some of these factors, i.e. Sld3 and Cdc45, can be detected specifically at some early firing origins during G1 [11, 18, 19], which is critical for promoting their activation at the beginning of S phase. However, even though the sequential activation of origins can be explained by limiting levels of essential initiation proteins, an open question is why limiting factors are preferentially recruited by early origins.
Recent reports have indicated that the Forkhead transcription factors Fkh1 and Fkh2 are involved in early origin firing [11, 17]. Although FKH1 and FKH2 are not required for origin licensing, the binding of Cdc45 to early origins in G1 seems to be Fkh1/2-dependent, and it has been suggested that the role of Fkh1/2 in origin firing would involve the establishment of replication timing domains and the spatial organization of origins . Fkh1/2 binding sites are enriched in a subset of origins, but ChIP–chip analyses indicate that there is no strict correlation between the presence of consensus motifs and Fkh1/2 binding [20, 21]. In addition, the presence of Fkh1/2 binding sites is not sufficient to confer early firing, and the number and distribution of Fkh motifs within the origin region seem to be important for conferring early origin activation [11, 17], arguing for an additional role of the chromatin environment in the regulation of origin firing.
Chromatin modifications also play a role in the regulation of origin timing and several studies have demonstrated that increases in histone acetylation, either by deletion of the deacetylase RPD3 or by artificial recruitment of the histone acetylase Gcn5, correlate with advanced origin firing [9, 10, 12, 22]. However, the explanation as to how histone acetylation modifies replication timing has remained obscure.
Previous studies have shown that replication origins in budding yeast exhibit a specific nucleosomal configuration, consisting of a nucleosome-depleted region (NDR) next to the ACS element, flanked by well-positioned nucleosomes [23, 24]. The correct organization of these origin-associated chromatin features is important for efficient replication initiation  and it would also be important for the timing of origin activation.
In light of the above, here we examined the contribution of the chromatin structure to the timing of origin firing by genome-wide analysis of the nucleosome organization at early and late origins and its relationship with epigenetic modifications, in particular with histone acetylation. Our findings show that early and late origins exhibit different nucleosomal patterns, which are influenced by the genomic context. We found that higher histone acetylation in rpd3 mutants was linked to local changes in chromatin structure, leading to a more open conformation at the affected replication origins. Finally, our analysis indicates that origin structure is dynamic during the cell cycle.
Early and late replication origins exhibit different nucleosome profiles
Since replication origins displayed different types of nucleosomal configuration as well as a different timing of activation, we wondered whether early and late origins might be preferentially linked to some of the nucleosome profiles identified above. First, we independently aligned nucleosome maps corresponding to early and late origins groups relative to the ACS. The assignment of replication origins within early or late categories for our analysis was made according to timing data available from at least three independent studies (OriDB; [12, 14, 27–29] (Additional file 1: Table S2). Comparison of their average nucleosome profiles showed that early origins displayed a higher occupancy of nucleosomes immediately upstream (-1) and downstream (+1) from a broader NDR feature, along with a better positioning of the adjacent nucleosomes (Figure 1C).
Second, we examined the nucleosomal arrangement of early and late replication origins within the six types of nucleosome patterns described above. Although the 6 different classes of origin organization were found within early and late origin groups (Additional file 2: Figure S1A), quantification of the distribution of early and late origins between the different classes revealed a marked preferential distribution of early origins for classes 1 and 2 and of late origins for classes 5 and 6 (Figure 1D). Indeed, 60% of early origins belonged to classes 1 and 2, while 65% of late origins belonged to classes 5 and 6. Moreover, nearly 40% of early or late origins were clustered in class 1 or 6 respectively (asterisks in Figure 1D). The overlap between the profiles of classes 1 and 6, which include the highest proportion of early and late origins, shows a striking difference in chromatin architecture between them (Figure 1E).
Next, we wondered whether early and late origins, or their different chromatin structure, could be influenced by the genomic context in which origins are placed. In S. cerevisiae, replication origins are mainly located at intergenic regions (IGR), which may be flanked by genes transcribed co-directionally, convergently or divergently. Replication origins in budding yeast appeared in these three types of IGR, but when we quantified the number of early and late origins located in each type of IGR, we found that they were not randomly distributed between them. In fact, while early origins can be found in all types of IGR with nearly the expected probability, late origins showed a strong bias for convergent IGR (Additional file 2: Figure S1B).
We asked whether this preferential localization would be related to their different chromatin architecture. To address this question we quantified the number of C1&2-early and C5&6-late origins located in each type of IGR, and we found that 60% of origins in each group mapped to IGR between tandem and convergent transcripts, respectively (Figure 1F). Moreover, class 6 was particularly enriched in convergent IGR, and the percentage of C6-late origins at this type of IGRs reached 71%. Our results indicated that class 6 was linked to this type of gene organization (convergent transcription), and class 6 is highly enriched with late origins, therefore late origins are necessarily over-represented at convergent IGR. Origin-free NDR in convergent IGR were on average shorter and presented lower nucleosome occupancy than NDR located in IGR with gene units organized in tandem or divergently (Additional file 2: Figure S1C), resembling class 6-origin patterns. This might suggest that IGR with convergent transcription are prone to have lower nucleosome occupancy and shorter NDR.
In sum, although a nucleosome-depleted region next to the ACS element flanked by positioned nucleosomes is a common feature shared by all origins, the above analysis identified previously unknown separate modes of nucleosome organization in early and late replication origins. In addition, these early and late-nucleosomal patterns were preferentially located in different kinds of intergenic regions, which might condition their different chromatin structure and timing.
Fkh-binding sites co-localize with the NDR in early origins
These results advance our previous knowledge by showing that Fhk1/2 motifs are often positioned at the boundaries of the NDR, which could explain why critical Fkh1/2 motifs are restricted to the close vicinity of ACS and why they are not further apart than 120 bp. Therefore, the localization and distribution of Fkh1/2 motifs required for early origin function might be imposed or determined by the NDR element.
These data also suggest that the location of Fkh motifs within accessible chromatin is important for early origin firing. It would be possible that the higher accessibility of Fkh1/2 motifs located within the NDR may favour the binding of Fkh1/2 factors, or/and the interaction between Fkh and replication factors such as Cdc45 in early origins [11, 17]. However, not all early origins contain Fkh-binding sites, and only 25.6% of the early origins considered in our analysis had this specific configuration of 2 Fkh sites close to the ACS, indicating that other factors must influence early firing too.
An increase in histone acetylation correlates with local nucleosomal remodelling in rpd3mutants
Histone modifications also play a role in determining origin timing, and an increase in histone acetylation around origins induces earlier firing . Our data revealed a correlation between origin structure and replication timing and we therefore wondered whether histone acetylation might modify replication timing by altering the chromatin architecture of replication origins. To examine this possibility, we generated genome-wide nucleosome maps in an rpd3Δ mutant, where the RPD3 histone deacetylase gen had been deleted.
Deletion of the deacetylase RPD3increases origin accessibility
To explore the above possibility, we wondered whether chromatin structure in origin regions would also be modified locally in an Rpd3-dependent manner. Analysis of the average nucleosomal organization of replication origins from G1-synchronized cells showed that origins presented a NDR flanked by well-positioned nucleosomes in rpd3Δ cells too and, also, that early and late origins displayed different nucleosome profiles when analyzed separately (Additional file 4: Figure S3A). This result indicated that the different profiles observed at early and late origins were Rpd3-independent.
To discriminate whether the G1-specific nucleosomal remodelling observed was affecting early or late replication origins or both, we analyzed the average nucleosome profiles of early and late origin groups independently. Figure 4D-E shows that chromatin structure was remodelled in both types of origins, although changes in the NDR width and occupancy in late origins were slightly more pronounced. According to this, the NDR size increased in rpd3Δ mutants in both early and late origins (Figure 4F). Local modifications in the NDR region also included the aperture of an NDR that was occupied by a nucleosome in a wild-type strain. This NDR “opening” was taken place in the absence of Rpd3 in 21% of early and 39.2% of late origins with these characteristics, and an example is shown in Additional file 4: Figure S3D for ARS1410.
Chromatin remodelling in the de-acetylation mutant also affected the occupancy of nucleosomes flanking the NDR. These local modifications were examined in individual origin regions in the absence of Rpd3, and some origins such as ARS603 and ARS423 showed a noticeable decrease in occupancy in the +1 nucleosome (arrows in Figure 4G) while in others like ARS1413 the nucleosomes -1, -2 and -3 were also affected (Additional file 4: Figure S3E). When acetylation data were available [9, 31] for the origin analyzed, we found that the chromatin-remodelled region was hyper-acetylated (Figure 4G).
Our data indicate that Rpd3 affects the nucleosomal configuration of many early and late origins. Similarly, Rpd3 affects the replication timing of numerous origins, including early and late ones [9, 12]. We confirmed that origin replication timing was altered in our rpd3Δ sample by examining origin activation of the early origin ARS305 and the late origin ARS603 during S phase as described in Additional file 5: Figure S4. Thus, 2D gel analysis indicated that when RPD3 was deleted the bubble arc for ARS305 was most prominent at 30 min, whereas in the wild-type strain the maximum was observed at 45 min (Additional file 5: Figure S4A). In addition, we examined the position of the replication forks by monitoring the association of the Replication Protein A (RPA), a single-strand DNA-binding protein that has been shown to match the location of the replisome . Chromatin immunoprecipitation (ChIP) experiments showed that the association of Rpa1 peaked at 5 kb away from ARS305 in rpd3Δ cells while in wild-type cells Rpa1-association was restricted to the origin position (Additional file 5: Figure S4B). Since origin efficiency seems comparable in both strains, we concluded that ARS305 fires earlier in the mutant strain and therefore that more DNA will be synthesized before dNTP levels are depleted by HU. These data are in agreement with previous observations showing that the firing of ARS305 is advanced in rpd3 mutants . When we examined the activation of ARS603 we detected initiation structures (Additional file 5: Figure S4A), indicating that in rpd3Δ mutants late origins can escape checkpoint inhibition to a certain extent, as described . However, the lower bubble arc signal in ARS603, suggested that late origin firing before critical dNTP depletion was taking place in a small proportion of cells. Therefore, replication timing was altered in rpd3 mutants, and so was the nucleosomal organization around origin regions.
Then we took the opposite approach, and we analyzed the nucleosome profile of a described set of origins (Additional file 1: Table S2) with advanced replication in the absence of Rpd3 . Additional file 4: Figure S3C showed similar variations in the NDR-width and occupancy of flanking nucleosomes as described in Figures 4A-D-E. These results strongly suggest that nucleosome remodelling in rpd3 mutants is related with changes in replication timing.
Taken together, our results indicate that there is a correlation between changes in histone acetylation and replication timing with local modifications in origin architecture. Despite the fact that early and late specific chromatin patterns are conserved in rpd3Δ mutants, the remodelling of key origin features such as the NDR and -1 /+1 nucleosomes leads to a more open conformation. Although it has been speculated that changes in histone acetylation can influence replication timing by modifying chromatin structure [9, 10], our results demonstrate for the first time that higher levels of histone acetylation caused by the deletion of the deacetylase RPD3 increase NDR opening.
Origin architecture is dynamic during activation
As mentioned in the Introduction, origin activation involves the modification of chromatin-bound replicative complexes and we wondered whether this process might entail modifications in origin structure during S phase. To explore this possibility, we performed a genome-wide analysis of the dynamics of chromatin organization at replication origins before and after origin activation by comparing the nucleosomal pattern of synchronized populations in G1 and cells proceeding into S-phase in the presence of hydroxyurea (HU). The slow movement of replication forks in the presence of HU facilitates the analysis of replication dynamics and nucleosome organization during S-phase. Additionally, late origin firing is inhibited in the presence of HU by the replication checkpoint [33, 34], providing a system to compare chromatin organization in active or inactive origins within the same replication cycle.
The above observations revealed that the structure of the NDR is dynamic during replication but we then wondered why NDR-modulation would be different at early and late origins during S phase. We reasoned that the observed changes might be linked to origin activation, which in HU should take place exclusively at early origins. To confirm the unique activation of early origins in our experimental conditions, we performed 2D gel analysis of DNA replication intermediates on ARS305 and ARS603 as previously described . In ARS305, the replication bubble arc indicated origin activation (Figure 5E). By contrast, when the membranes were re-probed to analyze ARS603 activation, no initiation structures at this late origin were observed. Therefore, one possibility was that the NDR-conformation would be different in origins before and after origin firing (early origin = post-firing state or late origin = pre-firing state in S phase in the presence of HU).
We conclude that origin nucleosomal configuration is dynamic during the cell cycle and propose that origin activation is associated with local changes in chromatin organization. Changes in chromatin architecture at origins seem to be intimately linked to their timing of activation and thus chromatin modulation in early and late origins is separated temporally. It is possible that changes in chromatin organization in origins might be required to allow the correct assembly/disassembly of replication factors before and after origin activation.
To elucidate the role of chromatin structure in the timing of origin firing in budding yeast, we conducted a high-resolution analysis of genome-wide nucleosomal organization at origins. Our results revealed distinctive nucleosomal configurations in early and late replication origins. The high resolution of our nucleosome maps, together with the independent analysis of early and late origins clusters were critical to disclose these different nucleosomal architectures, emphasizing the notion that the average nucleosome profile of all replication origins hides important variations in origin structure.
Previous studies have suggested that genomic regions upstream from genes tend to replicate early, while regions containing a gene end tend to replicate late . In agreement with this, our analysis indicates that there is a strong coincidence of the major late patterns (5&6) at convergent IGR that contain two gene ends, and moreover, they seem to be almost completely excluded from regions that do not contain any gene end (divergent IGR). The situation seems to vary further for early origin patterns, because these can be found in regions that coincide with transcription initiation or termination, or both. However, it is unclear why late origin patterns are highly represented at convergent IGR. It is possible that a higher rate of transcription termination at convergent IGR may interfere with nucleosome stability, leading to the lower occupancy observed in late origin patterns. This possibility is supported by the fact that other convergent IGR, non-associated with replication origins, also presented low nucleosome occupancy.
Alternatively, the co-localization of origins and transcription initiation sites (TSS) may promote early firing. In this sense, replication origins from several organisms frequently contain transcription factor-binding sites , arguing for a role of TF in origin regulation. In S. cerevisiae, the transcription factor Abf1 seems to help origin activation , and importantly, the transcription factor Fkh is involved in early origin firing [11, 17]. Therefore, it is possible that origins located in the vicinity of a promoter would be more likely to bind Fkh1/2 transcription factors, which can strongly influence the selection of these origins for early firing. Supporting this possibility, 81% of early origins enriched with Fkh1/2-binding sites were in IGRs containing at least one promoter, while only 19% were in promoter-free regions. However, only 25.6% of the early origins analyzed contain Fkh-binding sites, suggesting that the role of TF in origin activation might not be exclusive to Fkh and that other TFs could also contribute to early firing in other origins. In fact, binding sites for other transcription factors have been found close to replication origins and have been implicated in origin activation in both yeast and mammals [39–42]. Future examination of the interactions between transcription and replication factors should provide important information about this issue.
Regarding the distribution of Fkh-binding motifs, the overlapping of predominant sites with the NDR origin feature suggests that the location of Forkhead sites within an open/accessible chromatin region is important for origin regulation. In addition, the particular localization and spacing of Fkh sites within NDR raise the possibility that the binding of TFs such as Fkh1/2 might also contribute to the maintenance of a nucleosome-free region and to the positioning of proximal nucleosomes, similar to the previously described role of ORC-binding .
The fact that both ORC and Fkh1/2-binding sites are restricted to the NDR suggests that the NDR is critical for accommodating functional complexes in the origins. The average NDR size from origin patterns 1–6 varies from 193 to 129 bp, indicating that different NDR widths are valid for origin activity. However, classes 1 and 2, with the broadest NDRs, and classes 5 and 6, with the narrowest NDRs, generally correlate with early and late replication respectively, indicating that NDR width might influence the timing of origin replication.
We propose that the different affinities of initiation factors for early and late origins would be determined by their different origin structure. Thus, initiation factors might have a higher affinity for early origins that generally adopt class 1 and 2 conformations and show broader and better defined NDRs. The importance of an appropriate NDR width in origin timing is highlighted by the observation that previously described dormant origins such as ARS301, ARS302, ARS303 and ARS320 show either a closed NDR or an extremely small NDR (97 bp) (Additional file 8: Figure S7). Dormant origins bind a pre-RC and are licensed in G1 , and if given enough time they are competent to initiate replication [43, 44]. This small NDR size probably imposes a critically low affinity for regulatory factors in these origins, being excluded for firing in a normal replication cycle. The higher occupancy by -1/+1 nucleosomes in early origins suggests that the interaction of the initiation factors with the origin region or with some components of the pre-RC complex might be facilitated or mediated by interaction with these nucleosomes. In fact, previous studies have reported an interaction between transcriptional regulators and +1 and +2 nucleosomes [45, 46].
Our results also revealed that origin architecture was actively modulated during origin firing. To visualize replication-dependent origin remodelling it is important to compare high-resolution maps from synchronized G1 and S samples, perhaps explaining why it passed unnoticed in previous studies performed on asynchronous cultures and/or with tiling microarrays [23, 24]. An interesting possibility is that origin structure might be influenced by the different replication complexes or factors bound to the origin before and after origin firing. This is in agreement with the fact that some pre-RC components and/or initiation factors are bound and released from chromatin before and after origin activation [2, 3]. Interestingly, some initiation factors such as Sld3 and Cdc45 can be detected specifically at early firing origins during G1 [11, 18, 19]. These essential initiation factors are present in the cell in limiting amounts [15, 16], and it has been suggested that they are recycled from early to late origins during S phase. Therefore, the dynamic of nucleosomes and replication proteins during the process of origin activation might be related, and we are currently examining this possibility.
In addition, we also examined the contribution of histone acetylation to early origin firing at the level of chromatin organization. Our data indicate that increased levels of histone acetylation around origins induce local modifications in chromatin organization, creating a more accessible structure. One possibility is that histone acetylation might weaken histone/DNA interactions, diminishing the chromatin association of ACS-flanking nucleosomes and altering NDR formation. The resulting open conformation in origins might facilitate, either directly or indirectly (through interaction with other chromatin-bound factors), the recruitment of critical replication factors to origins.
Taken together, the results reported argue for an important contribution of the chromatin structure to origin selection and firing, and help to understand the interplay between replication regulators such as Fkh1/2, epigenetic modifications, and chromatin structure.
Coordinated origin firing is important to ensure complete replication of the entire genome and recent evidences suggest that the replication temporal program can influence genome stability [47, 48]. In this work, we have analyzed the role of chromatin structure in replication timing at the level of nucleosome organization. Our data highlight the importance of chromatin features, such as the NDR and nucleosomes +1 and -1, in the regulation of origin activation and timing. Thus, early origins present broader NDR features and a high occupancy of flanking nucleosomes while tighter NDR and a lower occupancy appear in late origins. Besides, factors involved in the establishment of origin timing such as histone deacetylation induce local modifications in these key origin features, leading to a more open origin conformation. Similarly, Fkh transcription factors seem to be in close relationship with origin structure, as suggested by the overlap of predominant Fkh-binding sites with the NDR origin feature.
Origin activation is also associated with local changes in chromatin organization. Although the specific chromatin conformations of origins are already established in G1, origin structure is dynamic during origin activation. It has been described that chromatin structure affects origin function , and it would be possible that origin architecture can also affect the affinity and efficiency of the formation of replication complexes and, therefore, influence the timing of origin firing.
Yeast strains and growth conditions
The strains used in this study are described in Additional file 9: Table S1 and derived from W303-1a (MAT a ade2-1 ura3-1 his3-11 15 trp1-1 leu2-3 112 can1-100). Yeast cultures were grown in YP medium (1% yeast extract, 2% bacto peptone) supplemented with 2% glucose (YPD) at 23°C to a final concentration of 0.7 × 107 cells/ml. Cells were synchronized in G1 phase by adding 5 μg/ml α factor (αF) for one generation time and released into fresh medium in the presence/absence of 0.2 M hydroxyurea (HU) depending on the experiment.
Preparation of mononucleosomal DNA
Previously published nucleosome maps were mostly obtained from asynchronous cultures, which would make difficult to detect nucleosome dynamics during different stages of the cell cycle. To overcome this limitation, we generated nucleosome maps from G1 and S-phase synchronized samples. Thus, 2 × 109 cells were collected at G1 and at 60 minutes in S-phase in the presence of 0.2 M HU for preparation of mononucleosomal DNA and for ChIP experiments (see below). Preparation of mononucleosomal DNA was adapted from . Spheroplasts were induced by treatment with 10 mg of zymolyase 20 T for 10 min at 30°C, and mononucleosomal fragments were generated by digesting DNA with 600 or 450 units/ml (WT and rpd3 strains respectively) of micrococcal nuclease (MNAse) at 37°C during 10 minutes. The amount of MNase was optimized experimentally for each strain to generate an 80:20 ratio of mononucleosomes to dinucleosomes, as described .
Sequencing and nucleosome maps representation
Mononucleosomal DNA was sequenced using the single-end sequencing protocol in an Illumina Genome Analyzer IIx. All sequencing data are deposited in the Gene Expression Omnibus (GEO) database under the accession number [GSE54377]. 23117530 to 28511688 single reads 40 nucleotides long were aligned with the reference genome of S. cerevisiae S288C (R64-1-1 assembly 3/2/2011). For the alignment we used Bowtie with no mismatches allowed. and the average genome coverage obtained was from 76 to 95-fold. This high level of resolution is a critical factor to detect subtle local modifications at individual nucleosome positions. The sequencing coverage for every nucleotide was normalized by the average genomic coverage. The nucleosome positioning profile was generated as described . Briefly, the average spacing between boundary peaks for individual nucleosomes was calculated from the smoothed signal generated with Python “Pywavelets” (with the multilevel 1-D biorthogonal wavelet decomposition/reconstruction tool). Then, the midpoint position of each nucleosome was defined and the resulting combined profile was wavelet-smoothed. For the copy number correction shown in Additional file 6: Figure S5, the ratio of mapped reads from the G1 and S-phase nucleosomal samples to the naked-DNA samples from the same time points (SRR398641 and SRR398652 files deposited in SRA database) was calculated. The number of reads from the G1 and S naked-DNA samples were corrected to obtain an average number value across the genome of one, as previously described in replication timing profiles analysis . The resulting signal was then smoothed with a 1000 bp window, step 1 bp.
Clustering, analysis of NDRs, TSS and search for Fkh motif
A list of replication origins with annotated ACS in the +/-20 nt limits of each origin was obtained from OriDB and those origins that mapped on a repetitive region were excluded. To generate the average nucleosome profiles, the nucleosome maps were aligned relative to the ACS and were oriented with respect to the T-rich strand. The nucleosome profiles of the [-600, 600] nt, centered on each ACS sequence, from the WT G1 sample were clustered using a k-mean clustering function with 10.000 iterations. We used 4 to 8 groups in our analysis, and six groups were finally chosen because they accomplished the best consensus between cluster size and information. Nucleosome-depleted regions (NDRs) were defined as regions spanning at least 80 nucleotides as described , with a normalized sequence coverage lower than 0.4. The annotations of  were used as a TSS reference. The whole genome was searched for the FKH motif (RTAAAYA) previously described  by using the tool FIMO with a maximum threshold of 0.001. We cured it by deleting all non-strict coincidences with the motif.
Two-dimensional gel analysis
DNA was isolated from 2 × 109 cells using a G-20 column (QIAGEN) according to the manufacturer’s instructions. DNA fragments digested using XbaI and EcoRI (Fermentas) or BamHI and NcoI (Fermentas) were separated by neutral/neutral two-dimensional agarose gel electrophoresis as described .
Chromatin Immunoprecipitation (ChIP) and quantitative real-time PCR (qPCR)
ChIP extracts were prepared from 7.5 × 108 cells as described . Immunoprecipitation was carried out overnight with 5 mg of anti-PK antibody (SV5-Pk1 Fisher Scientific Pierce) coupled to magnetic beads (Dynabeads, Invitrogen). Uncoupled Dynabeads were used for background controls. The DNA fragments recovered from the whole-cell extract and ChIP samples were quantified by real-time PCR using the BioRad-CFX96 system. The q-PCR signals obtained from the ChIP samples were normalized by the signals obtained from an input sample (percent IP). Primer sequences and amplification reactions are indicated in Additional file 9: Table S1. Primers were tested to achieve amplification curves with a slope of 3.3 ± 10% and R2 value >0.99. All real-time PCRs were performed in duplicate with SYBR Premix Ex Taq (TaKaRa) in two independent experiments.
Digestion with MNase and indirect end-labelling analyses
2×109 cells were permeabilized as described above and spheroplasts were split into six fractions in NP-buffer and treated with increasing amounts of MNase (Fermentas) (0, 0.5, 1, 1.5, 2 and 3 units/ml) for 10 minutes at 37°C. DNA was purified and quantified. For each sample, 2.5 μg of DNA were digested with HindIII and PvuII enzymes (Fermentas), separated by electrophoresis in a 1.5% agarose gel, blotted and hybridized to an specific end-terminal probe.
Availability of supporting data
The data sets supporting the results of this article are available in the [GSE54377] repository, [http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE54377].
We thank Rodrigo Bermejo for critical reading of the manuscript and for providing strains. We also thank Orlando Domínguez for technical support with high throughput sequencing. This work was supported by the Spanish Ministry of Economy and Competitiveness [RYC-2009-04407 and BFU2010-20034 to M.S., BFU2011-28804 and Consolider-Ingenio CSD2007-00015 to F.A.]; Spanish Ministry of Economy and Competitiveness (Ramón y Cajal contract to M.S.); Spanish Ministry of Economy and Competitiveness (predoctoral fellowships to I.S. and E.C.); Consejo Superior de Investigaciones Cientificas (contract to E.V.). Funding for open access charge: Spanish Ministry of Economy and Competitiveness.
- Bell SP, Stillman B: ATP-dependent recognition of eukaryotic origins of DNA replication by a multiprotein complex. Nature. 1992, 357: 128-134. 10.1038/357128a0.PubMedView ArticleGoogle Scholar
- Remus D, Diffley JF: Eukaryotic DNA replication control: lock and load, then fire. Curr Opin Cell Biol. 2009, 21: 771-777. 10.1016/j.ceb.2009.08.002.PubMedView ArticleGoogle Scholar
- Labib K: How do Cdc7 and cyclin-dependent kinases trigger the initiation of chromosome replication in eukaryotic cells?. Genes Dev. 2010, 24: 1208-1219. 10.1101/gad.1933010.PubMed CentralPubMedView ArticleGoogle Scholar
- Santocanale C, Sharma K, Diffley JF: Activation of dormant origins of DNA replication in budding yeast. Genes Dev. 1999, 13: 2360-2364. 10.1101/gad.13.18.2360.PubMed CentralPubMedView ArticleGoogle Scholar
- Wyrick JJ, Aparicio JG, Chen T, Barnett JD, Jennings EG, Young RA, Bell SP, Aparicio OM: Genome-wide distribution of ORC and MCM proteins in S. cerevisiae: high-resolution mapping of replication origins. Science. 2001, 294: 2357-2360. 10.1126/science.1066101.PubMedView ArticleGoogle Scholar
- Raghuraman MK, Brewer BJ, Fangman WL: Cell cycle-dependent establishment of a late replication program. Science. 1997, 276: 806-809. 10.1126/science.276.5313.806.PubMedView ArticleGoogle Scholar
- Friedman KL, Brewer BJ, Fangman WL: Replication profile of Saccharomyces cerevisiae chromosome VI. Genes Cells. 1997, 2: 667-678. 10.1046/j.1365-2443.1997.1520350.x.PubMedView ArticleGoogle Scholar
- Ferguson BM, Fangman WL: A position effect on the time of replication origin activation in yeast. Cell. 1992, 68: 333-339. 10.1016/0092-8674(92)90474-Q.PubMedView ArticleGoogle Scholar
- Vogelauer M, Rubbi L, Lucas I, Brewer BJ, Grunstein M: Histone acetylation regulates the time of replication origin firing. Mol Cell. 2002, 10: 1223-1233. 10.1016/S1097-2765(02)00702-5.PubMedView ArticleGoogle Scholar
- Aparicio JG, Viggiani CJ, Gibson DG, Aparicio OM: The Rpd3-Sin3 histone deacetylase regulates replication timing and enables intra-S origin control in Saccharomyces cerevisiae. Mol Cell Biol. 2004, 24: 4769-4780. 10.1128/MCB.24.11.4769-4780.2004.PubMed CentralPubMedView ArticleGoogle Scholar
- Knott SR, Peace JM, Ostrow AZ, Gan Y, Rex AE, Viggiani CJ, Tavare S, Aparicio OM: Forkhead transcription factors establish origin timing and long-range clustering in S. cerevisiae. Cell. 2012, 148: 99-111. 10.1016/j.cell.2011.12.012.PubMed CentralPubMedView ArticleGoogle Scholar
- Knott SR, Viggiani CJ, Tavare S, Aparicio OM: Genome-wide replication profiles indicate an expansive role for Rpd3L in regulating replication initiation timing or efficiency, and reveal genomic loci of Rpd3 function in Saccharomyces cerevisiae. Genes Dev. 2009, 23: 1077-1090. 10.1101/gad.1784309.PubMed CentralPubMedView ArticleGoogle Scholar
- Pohl TJ, Brewer BJ, Raghuraman MK: Functional centromeres determine the activation time of pericentric origins of DNA replication in Saccharomyces cerevisiae. PLoS Genet. 2012, 8: e1002677-10.1371/journal.pgen.1002677.PubMed CentralPubMedView ArticleGoogle Scholar
- Raghuraman MK, Winzeler EA, Collingwood D, Hunt S, Wodicka L, Conway A, Lockhart DJ, Davis RW, Brewer BJ, Fangman WL: Replication dynamics of the yeast genome. Science. 2001, 294: 115-121. 10.1126/science.294.5540.115.PubMedView ArticleGoogle Scholar
- Mantiero D, Mackenzie A, Donaldson A, Zegerman P: Limiting replication initiation factors execute the temporal programme of origin firing in budding yeast. EMBO J. 2011, 30: 4805-4814. 10.1038/emboj.2011.404.PubMed CentralPubMedView ArticleGoogle Scholar
- Tanaka S, Nakato R, Katou Y, Shirahige K, Araki H: Origin association of Sld3, Sld7, and Cdc45 proteins is a key step for determination of origin-firing timing. Curr Biol. 2011, 21: 2055-2063. 10.1016/j.cub.2011.11.038.PubMedView ArticleGoogle Scholar
- Looke M, Kristjuhan K, Varv S, Kristjuhan A: Chromatin-dependent and -independent regulation of DNA replication origin activation in budding yeast. EMBO Rep. 2013, 14: 191-198. 10.1038/embor.2012.196.PubMed CentralPubMedView ArticleGoogle Scholar
- Kamimura Y, Tak YS, Sugino A, Araki H: Sld3, which interacts with Cdc45 (Sld4), functions for chromosomal DNA replication in Saccharomyces cerevisiae. EMBO J. 2001, 20: 2097-2107. 10.1093/emboj/20.8.2097.PubMed CentralPubMedView ArticleGoogle Scholar
- Aparicio OM, Stout AM, Bell SP: Differential assembly of Cdc45p and DNA polymerases at early and late origins of DNA replication. Proc Natl Acad Sci U S A. 1999, 96: 9130-9135. 10.1073/pnas.96.16.9130.PubMed CentralPubMedView ArticleGoogle 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. 10.1038/nature02800.PubMed CentralPubMedView ArticleGoogle Scholar
- MacIsaac KD, Wang T, Gordon DB, Gifford DK, Stormo GD, Fraenkel E: An improved map of conserved regulatory sites for Saccharomyces cerevisiae. BMC Bioinformatics. 2006, 7: 113-10.1186/1471-2105-7-113.PubMed CentralPubMedView ArticleGoogle Scholar
- Stevenson JB, Gottschling DE: Telomeric chromatin modulates replication timing near chromosome ends. Genes Dev. 1999, 13: 146-151. 10.1101/gad.13.2.146.PubMed CentralPubMedView ArticleGoogle Scholar
- Berbenetz NM, Nislow C, Brown GW: Diversity of eukaryotic DNA replication origins revealed by genome-wide analysis of chromatin structure. PLoS Genet. 2010, 6: e1001092-10.1371/journal.pgen.1001092.PubMed CentralPubMedView ArticleGoogle Scholar
- Eaton ML, Galani K, Kang S, Bell SP, MacAlpine DM: Conserved nucleosome positioning defines replication origins. Genes Dev. 2010, 24: 748-753. 10.1101/gad.1913210.PubMed CentralPubMedView ArticleGoogle Scholar
- Lipford JR, Bell SP: Nucleosomes positioned by ORC facilitate the initiation of DNA replication. Mol Cell. 2001, 7: 21-30. 10.1016/S1097-2765(01)00151-4.PubMedView ArticleGoogle Scholar
- Nieduszynski CA, Knox Y, Donaldson AD: Genome-wide identification of replication origins in yeast by comparative genomics. Genes Dev. 2006, 20: 1874-1879. 10.1101/gad.385306.PubMed CentralPubMedView ArticleGoogle Scholar
- Feng W, Collingwood D, Boeck ME, Fox LA, Alvino GM, Fangman WL, Raghuraman MK, Brewer BJ: Genomic mapping of single-stranded DNA in hydroxyurea-challenged yeasts identifies origins of replication. Nat Cell Biol. 2006, 8: 148-155. 10.1038/ncb1358.PubMed CentralPubMedView ArticleGoogle Scholar
- Yabuki N, Terashima H, Kitada K: Mapping of early firing origins on a replication profile of budding yeast. Genes Cells. 2002, 7: 781-789. 10.1046/j.1365-2443.2002.00559.x.PubMedView ArticleGoogle Scholar
- Alvino GM, Collingwood D, Murphy JM, Delrow J, Brewer BJ, Raghuraman MK: Replication in hydroxyurea: it’s a matter of time. Mol Cell Biol. 2007, 27: 6396-6406. 10.1128/MCB.00719-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Vogelauer M, Wu J, Suka N, Grunstein M: Global histone acetylation and deacetylation in yeast. Nature. 2000, 408: 495-498. 10.1038/35044127.PubMedView ArticleGoogle Scholar
- Robyr D, Suka Y, Xenarios I, Kurdistani SK, Wang A, Suka N, Grunstein M: Microarray deacetylation maps determine genome-wide functions for yeast histone deacetylases. Cell. 2002, 109: 437-446. 10.1016/S0092-8674(02)00746-8.PubMedView ArticleGoogle Scholar
- De Piccoli G, Katou Y, Itoh T, Nakato R, Shirahige K, Labib K: Replisome stability at defective DNA replication forks is independent of S phase checkpoint kinases. Mol Cell. 2012, 45: 696-704. 10.1016/j.molcel.2012.01.007.PubMedView ArticleGoogle Scholar
- Santocanale C, Diffley JF: A Mec1- and Rad53-dependent checkpoint controls late-firing origins of DNA replication. Nature. 1998, 395: 615-618. 10.1038/27001.PubMedView ArticleGoogle Scholar
- Shirahige K, Hori Y, Shiraishi K, Yamashita M, Takahashi K, Obuse C, Tsurimoto T, Yoshikawa H: Regulation of DNA-replication origins during cell-cycle progression. Nature. 1998, 395: 618-621. 10.1038/27007.PubMedView ArticleGoogle Scholar
- Bell SP, Dutta A: DNA replication in eukaryotic cells. Annu Rev Biochem. 2002, 71: 333-374. 10.1146/annurev.biochem.71.110601.135425.PubMedView ArticleGoogle Scholar
- Friedman KL, Brewer BJ: Analysis of replication intermediates by two-dimensional agarose gel electrophoresis. Methods Enzymol. 1995, 262: 613-627.PubMedView ArticleGoogle Scholar
- Kohzaki H, Murakami Y: Transcription factors and DNA replication origin selection. Bioessays. 2005, 27: 1107-1116. 10.1002/bies.20316.PubMedView ArticleGoogle Scholar
- Diffley JF, Stillman B: Similarity between the transcriptional silencer binding proteins ABF1 and RAP1. Science. 1989, 246: 1034-1038. 10.1126/science.2511628.PubMedView ArticleGoogle Scholar
- Dominguez-Sola D, Ying CY, Grandori C, Ruggiero L, Chen B, Li M, Galloway DA, Gu W, Gautier J, Dalla-Favera R: Non-transcriptional control of DNA replication by c-Myc. Nature. 2007, 448: 445-451. 10.1038/nature05953.PubMedView ArticleGoogle Scholar
- Minami H, Takahashi J, Suto A, Saitoh Y, Tsutsumi K: Binding of AlF-C, an Orc1-binding transcriptional regulator, enhances replicator activity of the rat aldolase B origin. Mol Cell Biol. 2006, 26: 8770-8780. 10.1128/MCB.00949-06.PubMed CentralPubMedView ArticleGoogle Scholar
- Irlbacher H, Franke J, Manke T, Vingron M, Ehrenhofer-Murray AE: Control of replication initiation and heterochromatin formation in Saccharomyces cerevisiae by a regulator of meiotic gene expression. Genes Dev. 2005, 19: 1811-1822. 10.1101/gad.334805.PubMed CentralPubMedView ArticleGoogle Scholar
- Morse RH: RAP, RAP, open up! New wrinkles for RAP1 in yeast. Trends Genet. 2000, 16: 51-53. 10.1016/S0168-9525(99)01936-8.PubMedView ArticleGoogle Scholar
- Vujcic M, Miller CA, Kowalski D: Activation of silent replication origins at autonomously replicating sequence elements near the HML locus in budding yeast. Mol Cell Biol. 1999, 19: 6098-6109.PubMed CentralPubMedView ArticleGoogle Scholar
- Doksani Y, Bermejo R, Fiorani S, Haber JE, Foiani M: Replicon dynamics, dormant origin firing, and terminal fork integrity after double-strand break formation. Cell. 2009, 137: 247-258. 10.1016/j.cell.2009.02.016.PubMed CentralPubMedView ArticleGoogle 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. 10.1016/j.molcel.2009.09.011.PubMed CentralPubMedView ArticleGoogle Scholar
- Son J, Shen SS, Margueron R, Reinberg D: Nucleosome-binding activities within JARID2 and EZH1 regulate the function of PRC2 on chromatin. Genes Dev. 2013, 27: 2663-2677. 10.1101/gad.225888.113.PubMed CentralPubMedView ArticleGoogle Scholar
- Ibarra A, Schwob E, Mendez J: Excess MCM proteins protect human cells from replicative stress by licensing backup origins of replication. Proc Natl Acad Sci U S A. 2008, 105: 8956-8961. 10.1073/pnas.0803978105.PubMed CentralPubMedView ArticleGoogle Scholar
- Shima N, Alcaraz A, Liachko I, Buske TR, Andrews CA, Munroe RJ, Hartford SA, Tye BK, Schimenti JC: A viable allele of Mcm4 causes chromosome instability and mammary adenocarcinomas in mice. Nat Genet. 2007, 39: 93-98. 10.1038/ng1936.PubMedView ArticleGoogle 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-10.1371/journal.pbio.0060065.PubMed CentralPubMedView ArticleGoogle Scholar
- Lantermann A, Stralfors A, Fagerstrom-Billai F, Korber P, Ekwall K: Genome-wide mapping of nucleosome positions in Schizosaccharomyces pombe. Methods. 2009, 48: 218-225. 10.1016/j.ymeth.2009.02.004.PubMedView ArticleGoogle Scholar
- Soriano I, Quintales L, Antequera F: Clustered regulatory elements at nucleosome-depleted regions punctuate a constant nucleosomal landscape in Schizosaccharomyces pombe. BMC Genomics. 2013, 14: 813-10.1186/1471-2164-14-813.PubMed CentralPubMedView ArticleGoogle Scholar
- Muller CA, Hawkins M, Retkute R, Malla S, Wilson R, Blythe MJ, Nakato R, Komata M, Shirahige K, de Moura AP, Nieduszynski CA: The dynamics of genome replication using deep sequencing. Nucleic Acids Res. 2014, 42: e3-10.1093/nar/gkt878.PubMed CentralPubMedView ArticleGoogle Scholar
- Jiang C, Pugh BF: A compiled and systematic reference map of nucleosome positions across the Saccharomyces cerevisiae genome. Genome Biol. 2009, 10: R109-10.1186/gb-2009-10-10-r109.PubMed CentralPubMedView ArticleGoogle 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. 10.1038/ng2117.PubMedView ArticleGoogle Scholar
- Brewer BJ, Lockshon D, Fangman WL: The arrest of replication forks in the rDNA of yeast occurs independently of transcription. Cell. 1992, 71: 267-276. 10.1016/0092-8674(92)90355-G.PubMedView ArticleGoogle Scholar
- Bermejo R, Katou YM, Shirahige K, Foiani M: ChIP-on-chip analysis of DNA topoisomerases. Methods Mol Biol. 2009, 582: 103-118. 10.1007/978-1-60761-340-4_9.PubMedView ArticleGoogle 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.