Email updates

Keep up to date with the latest news and content from Genome Biology and BioMed Central.

Open Access Highly Accessed Research

Cell-cycle regulated transcription associates with DNA replication timing in yeast and human

Hunter B Fraser

Author Affiliations

Department of Biology, Stanford University, Stanford, CA 94305, USA

Genome Biology 2013, 14:R111  doi:10.1186/gb-2013-14-10-r111

The electronic version of this article is the complete one and can be found online at: http://genomebiology.com/2013/14/10/R111


Received:7 June 2013
Accepted:7 October 2013
Published:7 October 2013

© 2013 Fraser; licensee 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.

Abstract

Background

Eukaryotic DNA replication follows a specific temporal program, with some genomic regions consistently replicating earlier than others, yet what determines this program is largely unknown. Highly transcribed regions have been observed to replicate in early S-phase in all plant and animal species studied to date, but this relationship is thought to be absent from both budding yeast and fission yeast. No association between cell-cycle regulated transcription and replication timing has been reported for any species.

Results

Here I show that in budding yeast, fission yeast, and human, the genes most highly transcribed during S-phase replicate early, whereas those repressed in S-phase replicate late. Transcription during other cell-cycle phases shows either the opposite correlation with replication timing, or no relation. The relationship is strongest near late-firing origins of replication, which is not consistent with a previously proposed model—that replication timing may affect transcription—and instead suggests a potential mechanism involving the recruitment of limiting replication initiation factors during S-phase.

Conclusions

These results suggest that S-phase transcription may be an important determinant of DNA replication timing across eukaryotes, which may explain the well-established association between transcription and replication timing.

Background

The timing of DNA replication during S-phase of the cell cycle plays an important role in genome integrity, the mutational spectrum, and a wide range of human diseases [1]. Despite many recent advances in our ability to measure the time of replication (Trep) across entire genomes [2-7], our understanding of what regulates this timing remains far from complete [1,8-11]. The time at which origins of replication (ORIs) fire is thought to be determined in M-phase [12] or G1 [13,14], at which point factors such as Cdc45 and Sld3 bind to ORIs that will fire early in the following S-phase [15,16]. These and several other proteins critical for replication initiation are present at copy-numbers lower than the number of ORIs [17-19], and their over-expression advances Trep for many late-firing ORIs in both budding and fission yeast [12,17-20], suggesting that their re-use may be a key step in regulating ORI firing time. However, what determines the relative affinities of different ORIs for these limiting factors - and hence their temporal order of initiation - is largely unknown [19].

Among the strongest correlates (and potential determinants) of Trep in metazoans are transcriptional activity and chromatin state. Although transcriptionally active euchromatin has been known to replicate earlier than repressive heterochromatin for over 50 years [11,21], the reason - and even the direction of causation - has remained elusive. The two major models [8,11], not mutually exclusive, are that 1) the euchromatic chromatin structure is more permissive both to transcription and to DNA replication initiation, or 2) Trep itself affects chromatin structure and transcription as a result of changes in the nuclear milieu during S-phase. The former is most directly supported by experiments altering ORI firing time via manipulation of histone modifications [8-10,18,22-24], whereas the latter is supported by differences in chromatin and transcription of DNA templates injected into cells during either early or late S-phase [8,9,25,26].

Measuring Trep genome-wide in the budding yeast Saccharomyces cerevisiae (Sc), Raghuraman et al. [2] reported a surprising lack of association between transcription and Trep (with the exception of the eight histone genes, which are highly transcribed in S phase and are replicated early). However, this analysis only involved clusters of co-expressed genes, and did not actually compare the highest- versus lowest-expressed genes. Nevertheless, it has been widely interpreted in the literature as indicating the absence of any association, and many authors have speculated as to why budding yeast lacks this relationship [5,8-11]. Similarly, the fission yeast Schizosaccharomyces pombe (Sp) is thought to lack any association between transcription and replication timing [8], though again no systematic comparison has been reported.

Results

Because DNA replication is confined to a specific period during the cell cycle, I reasoned that the relationship between Trep and transcription may depend on when in the cell cycle transcription is occurring. The transcription of most genes does not vary greatly throughout the cell cycle, so cannot be used to determine phase-dependent effects. However, several hundred genes have been identified in both Sc and Sp for which transcription does vary consistently during the cell cycle [27,28]. I compared the expression levels of these cell cycle-regulated genes measured in synchronized cells [27,28] with the Trep for each gene to determine if any relationship exists. For both Sc and Sp expression levels measured in G2 phase, higher expression associated with earlier Trep (Figure  1A). However, at other points in the cell cycle the relationship was quite different; mostly notably in M/G1 (Sc) or G1 (Sp), the relationship reversed, such that highly expressed genes were replicated late (Figure  1A).

thumbnailFigure 1. The transcription/Trep association varies by cell-cycle stage. (A) Comparing mean Trep of the top decile (10%) of most-induced versus most-repressed cell cycle-regulated genes reveals that genes highly expressed in G2 replicate early in both Sc and Sp, whereas those highly expressed in M/G1 (Sc) or G1 (Sp) replicate late (error bars represent the standard error). (B) The correlation between Trep and expression levels of known cell cycle-regulated genes was calculated separately for expression levels from each time point of cell cycle-synchronized time courses [27,28]. An oscillation of the correlation coefficient (Pearson’s r) was observed for both budding yeast (all |r| > 0.107 are significant at P < 0.0025) and fission yeast (all |r| > 0.177 are at P < 0.0025). The approximate cell-cycle phase of each time point is shown [27,28]. Similar oscillations are observed for other methods of synchronization as well (Additional file 1: Figures S1 and S2). (C) A moving average of Trep is shown for all cell cycle-regulated genes, arranged in order of their time of maximal expression, beginning immediately following mitosis. A similar pattern is observed for both yeast species, with the latest Trep for genes with maximal transcript levels in G1, and the earliest Trep for genes with maximal transcript levels in G2.

Additional file 1: Figure S1. Correlation analysis as in Figure  1B (left), but using a different method for synchronization of Sc (a temperature-sensitive cdc28 mutant). Figure S2. correlation analysis as in Figure  1B (right), but using a different method for synchronization of Sp (a temperature-sensitive cdc25 mutant). Figure S3. repeating the moving average analysis from Figure  1C, with standard errors shown for each point (grey lines) in (A)Sc and (B)Sp. Results suggest the differences between high and low windows are unlikely to be due to random fluctuations. Figure S4.(A) The mean expression level of 100 genes comprising the window with earliest Trep in Figure  1C (left) is plotted as a function of time in the cell cycle. The genes that reach a maximum mRNA level in G2 have their maximum rate of increase (and likely maximum rate of transcription) in S phase. (B) As for part (A), but showing the mean expression for the 100 genes in the window with the latest Trep in Figure  1C (left). The genes that reach a minimum mRNA level in G2 have their maximum rate of decrease (and likely minimum rate of transcription) in S phase. Figure S5. repeating the ORI proximal/distal analysis from Figure  3C using a cutoff of 7.5 kb to define ORI proximal in (A)Sc and (B)Sp. Results are qualitatively identical to Figure  3C.

Format: PDF Size: 207KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

To more systematically visualize these patterns, I calculated the correlation between the expression levels of all cell cycle-regulated genes measured in synchronized cultures [27,28] with their Trep, separately for each expression data time-point (see Materials and methods). Plotting these correlation coefficients as a function of the time at which the expression data were sampled, I found a striking relationship: both the strength and direction of the correlation oscillate as a function of cell-cycle stage (Figure  1B). In these plots, positive r values represent time-points at which up-regulated genes tend to be replicated late in S phase; negative r values indicate times when up-regulated genes are replicated early. Consistent with the results in Figure  1A, in both species of yeast, genes highly expressed in G2 phase are replicated early, while those expressed in late M/G1 are replicated late. The oscillation is observed regardless of the method used to achieve cell-cycle synchronization (Additional file 1: Figures S1 and S2).

To further characterize this relationship, I plotted a moving average of Trep for the cell cycle-regulated genes in each species, ordered by their time of maximal expression. If expression in certain cell-cycle phases correlates with early or late replication, this should be reflected by troughs or peaks in such a plot. Again in both species a similar trend emerged: Trep reaches a maximum for genes expressed in G1, and a minimum for those expressed in G2 (Figure  1C; Additional file 1: Figure S3), consistent with the correlation analysis (Figure  1B). The strong conservation of this pattern was surprising, considering how much the regulation of DNA replication has diverged in the hundreds of millions of years separating these two yeast lineages [29].

Although the strongest association between high mRNA levels and early replication was observed for G2-phase expression levels, it is important to note that this does not imply these genes are maximally transcribed in G2. Rather, one would expect maximal transcription to occur in the time leading up to the maximal transcript level, that is, in S phase. Indeed, plotting mRNA levels for G2-upregulated genes (those with early Trep in Figure  1C), it is clear that their transcript levels show the greatest increase - likely reflecting active transcription - in S phase (Additional file 1: Figure S4A). Likewise, genes with late Trep show the opposite pattern: maximal decrease in mRNA levels during S phase (Additional file 1: Figure S4B).

The oscillating relationships shown in Figure  1 do not establish whether Trep is more directly associated with transcription in S phase or in M phase. For example, if M-phase repression led to early Trep, S-phase induction could be associated with early Trep simply as an indirect consequence, because genes repressed in M phase are typically induced in S phase (Additional file 1: Figure S4A). To disentangle the effects of S and M phases, I examined the Trep of genes that are expressed at similar levels throughout the cell cycle. If M-phase repression leads to early Trep, then genes repressed throughout the cell cycle would be expected to have early Trep, as a result of their repression in M phase (in this scenario, their S-phase expression levels are not relevant). However, if the association is instead due to S-phase induction, genes with constitutive high expression would have earlier Trep because of their active transcription in S phase (in which case M-phase expression levels would be irrelevant). This analysis showed a clear trend: highly expressed genes replicate 5.9 minutes earlier in Sc and 3.0 minutes earlier in Sp (Figure  2). Therefore, the results shown in Figure  1 can be entirely, and most parsimoniously, explained by the association of Trep with S-phase transcription; the M-phase relationship is likely to be an indirect side effect of this. This result also suggests a more general association between transcription and Trep in yeast that extends beyond cell cycle-regulated genes.

thumbnailFigure 2. Asynchronous gene expression associates with Trep in budding and fission yeast. Comparison of the 100 highest-expressed genes with the 100 lowest-expressed shows that highly expressed genes are replicated earlier in both budding yeast and fission yeast. Error bars represent the standard error.

To further investigate the connection between S-phase transcription and Trep, I tested whether the relationship differed for genes replicated during early versus late S phase. In this analysis I separated all cell cycle-regulated genes into 10 bins (that is, deciles) by their Trep, and plotted the median G2-phase transcript level (the time-point most closely reflecting S-phase transcription; Additional file 1: Figure S4A) for each. Across all 10 Trep bins, I observed the expected relationship: decreasing expression of bins with increasing Trep (Figure  3A). However, closer examination revealed that, for both yeast species, this pattern was almost entirely driven by late-replicating genes. In other words, there was no correlation between expression levels and Trep for genes in the first five bins (replicated in early S phase), while in late S phase the relationship was quite strong (Figure  3A). Consistent with this, applying the correlation analysis from Figure  1A to just early- or late-replicating genes revealed that the oscillation is entirely driven by replication in late S phase; genes with early Trep showed no oscillation, and only a weak correlation at nearly all time-points (Figure  3B). These results parallel the finding in mouse that genes replicated in the second half of S phase show the strongest association between transcript levels and Trep[7].

thumbnailFigure 3. Factors affecting the strength of the transcription/Trep association. (A) Median G2-phase transcript levels (representing S-phase transcription; Additional file 1: Figure S4A) are shown for all cell cycle-regulated genes separated into 10 equally sized bins (deciles) by their Trep. For both yeast species, no correlation is observed for the first five bins, whereas a strong relationship is present for later Trep. (B) Consistent with the decile analysis, no oscillation is observed in the correlation between expression level and Trep for early Trep genes, while a strong oscillation is observed for late Trep genes. (C) Only weak oscillation is observed in the correlation between expression level and Trep for ORI-distal genes (>5 kb from the nearest ORI in budding yeast, or 10 kb in fission yeast), while a strong oscillation is observed for ORI-proximal genes.

Another factor that may influence the relationship between S-phase transcription and replication timing is a gene’s distance from the nearest ORI. Under the model where chromatin affects both transcription and Trep, the strongest association would be expected for genes near ORIs, whereas if instead Trep affects a gene’s level of S-phase transcription, the relationship should be independent of distance to the nearest ORI [3]. Separating genes into two classes, ORI-proximal or ORI-distal, the ORI-proximal class showed far stronger oscillations (Figure  3C; ORI distance cutoffs, chosen to result in approximately equal-sized lists, were 5 kb from the nearest ORI in Sc and 10 kb in Sp, due to the higher density of known ORIs in Sc; results from equal distance cutoffs are shown in Additional file 1: Figure S5). Because ORI-proximal genes tend to be replicated earlier than ORI-distal genes, this result could not be an indirect effect of the stronger association for late Trep genes, as it acts in the opposite direction. This result suggests that the relationship is unlikely to be caused by an effect of Trep on S-phase transcription, which is one of the two major classes of models that have been proposed to explain the transcription/Trep association [8,25,26].

To test whether the relationship between S-phase transcription and replication timing is conserved outside of fungi, I applied the same correlation analysis to cell-cycle gene expression and Trep data from human HeLa cells [6,30]. Analyzing all known HeLa cell-cycle regulated genes [30], I found no significant relationship of any kind (Figure  4A). However, applying each of the two filters identified from yeast - late Trep and ORI proximity (within 10 kb) - resulted in clear and significant oscillations, of a magnitude similar to that observed for both yeast species (Figure  4B,C). As observed for yeast, the minimum correlation (indicating early Trep of up-regulated genes) occurred in G2, and the maximum in late M/G1. The fact that the same oscillating relationship exists in human, and that its strength is influenced by the same two factors, suggests that it is likely to be caused by a mechanism conserved between fungi and metazoans.

thumbnailFigure 4. Transcription and Trep in human. (A) No oscillation is observed when comparing the Trep versus expression levels of all cell cycle-regulated genes in HeLa cells (all |r| > 0.063 are significant at P < 0.05; the four time-points that exceed this are within the range expected by chance, given that 47 time-points were analyzed). (B) Significant oscillation is observed when comparing Trep versus expression levels of cell cycle-regulated genes with late Trep (red line; the final 50% of S phase; all |r| > 0.195 are significant), but not early Trep (blue line). (C) Significant oscillation is observed when comparing Trep versus expression levels of cell cycle-regulated genes within 10 kb of an ORI (blue line; all |r| > 0.197 are significant), but not further than 10 kb from an ORI (red line).

To put into perspective the strength of the relationship between Trep and cell cycle-regulated gene expression in human, I compared it to the well-established association between Trep and average (asynchronous) gene expression. The latter provides a useful benchmark because it is regarded as a strong relationship that has been observed in numerous studies across diverse metazoans [5,6,8,9]. To facilitate a direct comparison with the results in Figure  4, I used the same Trep data [6] for the same genes, but replaced the cell cycle-synchronized gene expression data [30] with high-coverage RNA-seq data from asynchronous HeLa cells [31]. The correlation between asynchronous expression and Trep was r = -0.16 for late Trep genes (the genes represented by the red line in Figure  4B) and r = -0.15 for ORI-proximal genes (represented by the blue line in Figure  4C). In both cases, the asynchronous data explained less than a third of the variance in Trep that is explained by S-phase transcription (see Materials and methods). Differing quality of the two gene expression data sets [30,31] could contribute to this difference; however, because RNA-seq is of far higher precision than spotted cDNA microarrays [32], any difference would likely underestimate the strength of the cell-cycle oscillations (Figure  4). These results suggest that the relationship between Trep and S-phase transcription in human is substantially stronger than the well-established association with asynchronous expression.

Discussion

These results suggest that 1) S-phase transcription is associated with DNA replication timing in budding yeast, fission yeast, and human; 2) the association is strongest for genomic regions near ORIs, excluding the causal model in which Trep affects transcription [8,9,25,26]; 3) it is also strongest for regions replicated in late S phase, implying that early-firing ORIs are not affected by this relationship; and 4) this association explains at least three times more of the variability in Trep than the well-known association with (asynchronous) gene expression in human.

Although the replication of these patterns across three species (and across multiple data sets within species; Additional file 1: Figures S1 and S2) lends confidence to their robustness, several caveats should be considered. First, gene expression was represented by transcript abundances, which is a function of both transcription and mRNA decay; therefore, the correlations reported here may underestimate the relationship between transcription and Trep. This prediction can be tested once rates of transcription have been measured throughout the cell cycle. Second, data quality is critical in any analysis; poor-quality data can reduce, or entirely mask, a real relationship. However, in most analyses reported here this is not a major concern, because it could only make the current results conservative (one exception to this is the ORI-proximal versus distal analyses (Figures  3C and 4C): if Trep was measured more accurately near ORIs, this would lead to stronger ORI-proximal correlations, although additional analysis suggests this is not the case (see Materials and methods)). Third, correlation does not imply causation. Although the evidence does not support a model where Trep affects transcription (Figures  3C and 4C), I cannot determine whether transcription itself is affecting Trep, or whether unobserved (latent) factors may be involved. With this caveat in mind, I believe there is still sufficient evidence to propose a testable model to account for these data.

A plausible mechanism explaining these observations draws from the finding that the firing of ORIs in late S phase is governed by recruitment of limiting replication initiation factors [12,17-20]. These factors are sequestered by early-firing ORIs from G1 until early S phase, and are reused at late-firing ORIs after their release from early-firing ORIs. I propose that the level of S-phase transcription near a late-firing ORI reflects local chromatin accessibility and/or subnuclear positioning, and in turn the ability of ORIs to recruit these limiting factors during S phase (Figure  5). This model accounts for the relationship of Trep with S-phase transcription (and the differing relationships in other phases); for the relationship being strongest near late-firing ORIs; and for the inferred direction of causality (that is, Trep not being causal).

thumbnailFigure 5. A model to explain these observations. Components: ORC and MCM2-7 are protein complexes comprising the pre-replicative complex. Blue cylinders represent nucleosomes, with dark blue indicating closed/repressive chromatin and light blue indicating open/accessible chromatin. Red proteins are limiting replication initiation factors (such as Cdc45 and Sld3). Txn = transcription. Sequence of events: in G1 (not depicted), the limiting replication initiation factors (red circles) associate with the earliest-firing ORIs (top row). When S phase begins, these early ORIs fire and release the factors, which are then free to associate with other ORIs (though note that Cdc45 is a component of the replication fork, so can only be recycled after fork termination). The relative affinities of the remaining ORIs for these factors - and thus their relative firing times - are determined by the chromatin state near the ORI during S-phase. ORIs near genes highly transcribed in S phase (middle row) have an accessible chromatin structure and thus high affinity, so will tend to fire earlier than those with little nearby S-phase transcription and thus less accessible chromatin (bottom row). Although not shown here, subnuclear positioning could help determine ORI accessibility, either by influencing chromatin structure or through other mechanisms. Figure adapted from [19].

The proposed mechanism likely acts in concert with other factors determining Trep, and thus is not inconsistent with evidence supporting these other factors. For example, although the determination of early versus late-firing ORIs is completed during M/G1 [12-14], S-phase transcription may still influence firing time specifically at late-firing ORIs (Figure  5).

Conclusions

Future work integrating these results with other (non-mutually-exclusive) mechanisms affecting Trep - for example, Forkhead transcription factors [33] and subnuclear positioning [8,29,34,35] - may lead to a unified framework for understanding the causes, and consequences, of the temporal program of DNA replication across eukaryotes.

Materials and methods

Data sources

Genome-wide Trep values were downloaded for all three species [2,4,6], and mapped onto genes by linear interpolation to the gene’s midpoint. Asynchronous yeast expression levels (used in Figure  2) were taken from [36,37], using the poly-A data for Sc and the median of wild-type replicates for Sp. Asynchronous HeLa RNA-seq data were from the ENCODE project [31]. Identities of cell cycle-regulated genes, their expression levels, and the cell-cycle phase of each expression time-point were acquired from [27,28,30]. All cell-cycle expression data were measured as mRNA levels relative to asynchronous levels of each gene, as opposed to absolute mRNA abundances that can be measured by RNA-seq; therefore, these expression levels represent the relative induction or repression of each gene throughout the cell cycle. The order of maximum expression levels was obtained from [38] for Sp and [27] for Sc. ORI locations were downloaded from ORIdb [39] for both yeasts (using only 'confirmed' or 'likely' ORIs), and from [6] for human (see below).

Yeast data analysis

All correlations were Pearson’s (significance cutoffs given in each figure legend). Trep moving averages (Figure  1C) were calculated for windows of 100 genes for Sc and 60 genes for Sp (due to the smaller number of cell cycle-regulated genes in Sp). For Figures  1A and 3A, the G2 expression data were represented by the 42 minute time-point for Sc and 135 minute time-point for Sp; for Figure  1A, Sc M/G1 was represented by the 70 minute time-point, and Sp G1 was represented by the 225 minute time-point. For Figure  3B, the early/late S-phase cutoff was chosen at halfway through S phase of each Trep data set (39.6 minutes after release from hydroxyurea arrest in Sp, and 26.8 minutes after release from cdc7 arrest in Sc). The cutoff for ORI-proximal versus ORI-distal (5 kb from each gene’s 5' end in Sc and 10 kb in Sp) was chosen in each yeast to result in gene lists of approximately equal size.

P-values in Figure  2 were calculated with a two-tailed Student’s t-test. Because the Sc expression levels were calculated as a ratio of mRNA/genomic DNA from asynchronous cells [37], they represent the number of mRNAs per DNA copy, and thus account for the fact that genes with early Trep spend a greater portion of the cell cycle with two copies. Although the Sp expression data [36] do not account for this, correcting for the effect by subtracting a fraction of each expression level proportional to the time each gene spends with two copies had only a minimal effect.

All code and data are available at [40].

Human data analysis

Human ORIs were defined as Orc1 binding sites [6] located within 1 mb of early-replicating peaks in the HeLa Trep profile, which indicate active ORIs (this window size was necessitated by the low resolution of the Trep profile) [6]. The early/late Trep cutoff was the first 50% of S phase and the ORI-proximal/distal cutoff was 10 kb from each gene’s 5' end. Due to the higher number of expression data points per cell cycle in human (approximately 15 in human versus approximately 9 for both yeasts), a two-point moving average was used for plotting human correlation coefficients.

To compare asynchronous expression versus S-phase transcription in HeLa cells, I compared high-coverage RNA-seq data from HeLa cells [31] with Trep[6] for the same genes analyzed in Figure  4B,C. The fraction of variance in Trep explained by the expression data is simply the r2 value from the Pearson’s correlation. Comparing these values for the asynchronous data with the strongest G2-phase (used to represent S-phase transcription, as described above) correlations, among the late-replicating genes (represented by the red line in Figure  4B) 2.7% of the variance in Trep was explained by the asynchronous data, versus 8.1% for S-phase transcription. Likewise for ORI-proximal genes (represented by the blue line in Figure  4C), the asynchronous data explained 2.3% of the variance in Trep, versus 7.6% for S-phase transcription.

To determine whether Trep is measured with greater accuracy near ORIs, I compared the Trep data used in Figure  4[6] with an independent Trep data set from HeLa cells [41]. Restricting the analysis to the cell cycle-regulated genes analyzed in Figure  4C, I found that ORI-distal genes actually showed better agreement between Trep data sets than did ORI-proximal genes (r = 0.59 and 0.46, respectively). This implies that, if anything, Trep is measured less accurately in ORI-proximal regions, which would lead to an underestimate of the strength of the oscillating correlation (blue line in Figure  4C).

Abbreviations

ORI: Origin of replication; Sc: Saccharomyces cerevisiae; Sp: Schizosaccharomyces pombe; Trep: Time of replication.

Competing interests

The author declares no competing interest.

Acknowledgements

I would like to thank M Botchan, N Cozzarelli, A Donaldson, D Gilbert, M Kobor, and J Rine for helpful advice. This work was supported by NIH grant 1R01GM097171-01A1. HBF is an Alfred P Sloan Fellow and a Pew Scholar in the Biomedical Sciences.

References

  1. Watanabe Y, Maekawa M: Spatiotemporal regulation of DNA replication in the human genome and its association with genomic instability and disease.

    Curr Med Chem 2010, 17:222. OpenURL

  2. 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. OpenURL

  3. Farkash-Amar S, Lipson D, Polten A, Goren A, Helmstetter C, Yakhini Z, Simon I: Global organization of replication time zones of the mouse genome.

    Genome Res 2008, 18:1562. OpenURL

  4. Eshaghi M, Karuturi RK, Li J, Chu Z, Liu ET, Liu J: Global profiling of DNA replication timing and efficiency reveals that efficient replication/firing occurs late during S-Phase in S. pombe.

    PLoS ONE 2007, 2:e722. OpenURL

  5. Schübeler D, Scalzo D, Kooperberg C, van Steensel B, Delrow J, Groudine M: Genome-wide DNA replication profile for Drosophila melanogaster: a link between transcription and replication timing.

    Nat Genet 2002, 32:438. OpenURL

  6. Dellino GI, Cittaro D, Piccioni R, Luzi L, Banfi S, Segalla S, Cesaroni M, Mendoza-Maldonado R, Giacca M, Pelicci PG: Genome-wide mapping of human DNA-replication origins: Levels of transcription at ORC1 sites regulate origin selection and replication timing.

    Genome Res 2013, 23:1. OpenURL

  7. Hiratani I, Ryba T, Itoh M, Yokochi T, Schwaiger M, Chang CW, Lyou Y, Townes TM, Schübeler D, Gilbert DM: Global reorganization of replication domains during embryonic stem cell differentiation.

    PLoS Biol 2008, 6:e245. OpenURL

  8. Hiratani I, Takebayashi S, Lu J, Gilbert DM: Replication timing and transcriptional control: beyond cause and effect - part II.

    Curr Opin Genet Dev 2009, 19:142. OpenURL

  9. Donaldson AD: Shaping time: chromatin structure and the DNA replication programme.

    Trends Genet 2005, 21:444. OpenURL

  10. McCune HJ, Donaldson AD: DNA replication: telling time with microarrays.

    Genome Biol 2003, 4:204. OpenURL

  11. Gilbert DM: Replication timing and transcriptional control: beyond cause and effect.

    Curr Opin Cell Biol 2002, 14:377. OpenURL

  12. Wu PYJ, Nurse P: Establishing the program of origin firing during S phase in fission yeast.

    Cell 2009, 136:852. OpenURL

  13. Raghuraman MK, Brewer BJ, Fangman WL: Cell cycle–dependent establishment of a late replication program.

    Science 1997, 276:806. OpenURL

  14. Dimitrova DS, Gilbert DM: The spatial position and replication timing of chromosomal domains are both established in early G1 phase.

    Mol Cell 1999, 4:983. OpenURL

  15. 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. OpenURL

  16. 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. OpenURL

  17. 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. OpenURL

  18. 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. OpenURL

  19. Douglas ME, Diffley JFX: Replication timing: the early bird catches the worm.

    Curr Biol 2012, 22:R81. OpenURL

  20. Patel PK, Kommajosyula N, Rosebrock A, Bensimon A, Leatherwood J, Bechhoefer J, Rhind N: The Hsk1(Cdc7) replication kinase regulates origin efficiency.

    Mol Biol Cell 2008, 19:5550. OpenURL

  21. Lima de Faria A: Incorporation of tritiated thymidine into meiotic chromosomes.

    Science 1959, 130:503. OpenURL

  22. Vogelauer M, Rubbi L, Lucas I, Brewer BJ, Grunstein M: Histone acetylation regulates the time of replication origin firing.

    Mol Cell 2002, 10:1223. OpenURL

  23. Goren A, Tabib A, Hecht M, Cedar H: DNA replication timing of the human beta-globin domain is controlled by histone modification at the origin.

    Genes Dev 2008, 22:1319. OpenURL

  24. Knott SR, Viggiani CJ, Tavaré 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. OpenURL

  25. Lande-Diner L, Zhang J, Cedar H: Shifts in replication timing actively affect histone acetylation during nucleosome reassembly.

    Mol Cell 2009, 34:767. OpenURL

  26. Zhang J, Xu F, Hashimshony T, Keshet I, Cedar H: Establishment of transcriptional competence in early and late S phase.

    Nature 2002, 420:198. OpenURL

  27. Spellman PT, Sherlock G, Zhang MQ, Iyer VR, Anders K, Eisen MB, Brown PO, Botstein D, Futcher B: Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization.

    Mol Biol Cell 1998, 9:3273. OpenURL

  28. Rustici G, Mata J, Kivinen K, Lió P, Penkett CJ, Burns G, Hayles J, Brazma A, Nurse P, Bähler J: Periodic gene expression program of the fission yeast cell cycle.

    Nat Genet 2004, 36:809. OpenURL

  29. Natsume T, Tanaka TU: Spatial regulation and organization of DNA replication within the nucleus.

    Chromosome Res 2010, 18:7. OpenURL

  30. Whitfield ML, Sherlock G, Saldanha AJ, Murray JI, Ball CA, Alexander KE, Matese JC, Perou CM, Hurt MM, Brown PO, Botstein D: Identification of genes periodically expressed in the human cell cycle and their expression in tumors.

    Mol Biol Cell 2002, 13:1977. OpenURL

  31. Djebali S, Davis CA, Merkel A, Dobin A, Lassmann T, Mortazavi A, Tanzer A, Lagarde J, Lin W, Schlesinger F, Xue C, Marinov GK, Khatun J, Williams BA, Zaleski C, Rozowsky J, Röder M, Kokocinski F, Abdelhamid RF, Alioto T, Antoshechkin I, Baer MT, Bar NS, Batut P, Bell K, Bell I, Chakrabortty S, Chen X, Chrast J, Curado J, et al.: Landscape of transcription in human cells.

    Nature 2012, 489:101-108. OpenURL

  32. Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics.

    Nat Rev Genet 2009, 10:57-63. OpenURL

  33. Knott SR, Peace JM, Ostrow AZ, Gan Y, Rex AE, Viggiani CJ, Tavaré S, Aparicio OM: Forkhead transcription factors establish origin timing and long-range clustering in S. cerevisiae.

    Cell 2012, 148:99-111. OpenURL

  34. Yaffe E, Farkash-Amar S, Polten A, Yakhini Z, Tanay A, Simon I: Comparative analysis of DNA replication timing reveals conserved large-scale chromosomal architecture.

    PLoS Genet 2010, 6:e1001011. OpenURL

  35. Ryba T, Hiratani I, Lu J, Itoh M, Kulik M, Zhang J, Schulz TC, Robins AJ, Dalton S, Gilbert DM: Evolutionarily conserved replication timing profiles predict long-range chromatin interactions and distinguish closely related cell types.

    Genome Res 2010, 20:761-770. OpenURL

  36. Wells ML, Huang W, Li L, Gerrish KE, Fargo DC, Ozsolak F, Blackshear PJ: Posttranscriptional regulation of cell-cell interaction protein-encoding transcripts by Zfs1p in Schizosaccharomyces pombe.

    Mol Cell Biol 2012, 32:4206. OpenURL

  37. Wang Y, Liu CL, Storey JD, Tibshirani RJ, Herschlag D, Brown PO: Precision and functional specificity in mRNA decay.

    Proc Natl Acad Sci U S A 2002, 99:5860. OpenURL

  38. Gauthier NP, Jensen LJ, Wernersson R, Brunak S, Jensen TS: Cyclebase.org: version 2.0, an updated comprehensive, multi-species repository of cell cycle experiments and derived analysis results.

    Nucleic Acids Res 2010, 38:D699. OpenURL

  39. Siow CC, Nieduszynska SR, Müller CA, Nieduszynski CA: OriDB, the DNA replication origin database updated and extended.

    Nucleic Acids Res 2012, 40:D682. OpenURL

  40. txn-rep.zip [http://www.stanford.edu/group/fraserlab/txn-rep webcite]

  41. Chen CL, Rappailles A, Duquenne L, Huvet M, Guilbaud G, Farinelli L, Audit B, D'Aubenton-Carafa Y, Arneodo A, Hyrien O, Thermes C: Impact of replication timing on non-CpG and CpG substitution rates in mammalian genomes.

    Genome Res 2010, 20:447-457. OpenURL