Universal temporal profile of replication origin activation in eukaryotes

Genome replication in asynchronously growing microbial populations

Biological cells replicate their genomes in a well-planned manner. The DNA replication program of an organism determines the timing at which different genomic regions are replicated, with fundamental consequences for cell homeostasis and genome stability. In a growing cell culture, genomic regions that are replicated early should be more abundant than regions that are replicated late. This abundance pattern can be experimentally measured using deep sequencing. However, a general quantitative theory linking this pattern to the replication program is still lacking. In this paper, we predict the abundance of DNA fragments in asynchronously growing cultures from any given stochastic model of the DNA replication program. As key examples, we present stochastic models of the DNA replication programs in budding yeast and Escherichia coli. In both cases, our model results are in excellent agreement with experimental data and permit to infer key information about the replication program. In particular, our method is able to infer the locations of known replication origins in budding yeast with high accuracy. These examples demonstrate that our method can provide insight into a broad range of organisms, from bacteria to eukaryotes.

Rif1 restrains the rate of replication origin firing in Xenopus laevis

Metazoan genomes are duplicated by the coordinated activation of clusters of replication origins at different times during S phase, but the underlying mechanisms of this temporal program remain unclear during early development. Rif1, a key replication timing factor, inhibits origin firing by recruiting protein phosphatase 1 (PP1) to chromatin counteracting S phase kinases. We have previously described that Rif1 depletion accelerates early Xenopus laevis embryonic cell cycles. Here, we find that in the absence of Rif1, patterns of replication foci change along with the acceleration of replication cluster activation. However, initiations increase only moderately inside active clusters. Our numerical simulations suggest that the absence of Rif1 compresses the temporal program towards more homogeneity and increases the availability of limiting initiation factors. We experimentally demonstrate that Rif1 depletion increases the chromatin-binding of the S phase kinase Cdc7/Drf1, the firing factors Treslin, MTBP, Cdc45, RecQL4, and the phosphorylation of both Treslin and MTBP. We show that Rif1 globally, but not locally, restrains the replication program in early embryos, possibly by inhibiting or excluding replication factors from chromatin.

DNA replication timing: Biochemical mechanisms and biological significance

The regulation of DNA replication is a fascinating biological problem both from a mechanistic angle-How is replication timing regulated?-and from an evolutionary one-Why is replication timing regulated? Recent work has provided significant insight into the first question. Detailed biochemical understanding of the mechanism and regulation of replication initiation has made possible robust hypotheses for how replication timing is regulated. Moreover, technical progress, including high-throughput, single-molecule mapping of replication initiation and single-cell assays of replication timing, has allowed for direct testing of these hypotheses in mammalian cells. This work has consolidated the conclusion that differential replication timing is a consequence of the varying probability of replication origin initiation. The second question is more difficult to directly address experimentally. Nonetheless, plausible hypotheses can be made and one-that replication timing contributes to the regulation of chromatin structure-has received new experimental support.

Organization of DNA Replication Origin Firing in Xenopus Egg Extracts: The Role of Intra-S Checkpoint

During cell division, the duplication of the genome starts at multiple positions called replication origins. Origin firing requires the interaction of rate-limiting factors with potential origins during the S(ynthesis)-phase of the cell cycle. Origins fire as synchronous clusters which is proposed to be regulated by the intra-S checkpoint. By modelling the unchallenged, the checkpoint-inhibited and the checkpoint protein Chk1 over-expressed replication pattern of single DNA molecules from Xenopus sperm chromatin replicated in egg extracts, we demonstrate that the quantitative modelling of data requires: (1) a segmentation of the genome into regions of low and high probability of origin firing; (2) that regions with high probability of origin firing escape intra-S checkpoint regulation and (3) the variability of the rate of DNA synthesis close to replication forks is a necessary ingredient that should be taken in to account in order to describe the dynamic of replication origin firing. This model implies that the observed origin clustering emerges from the apparent synchrony of origin firing in regions with high probability of origin firing and challenge the assumption that the intra-S checkpoint is the main regulator of origin clustering.

Genome-wide mapping of human DNA replication by optical replication mapping supports a stochastic model of eukaryotic replication

The heterogeneous nature of eukaryotic replication kinetics and the low efficiency of individual initiation sites make mapping the location and timing of replication initiation in human cells difficult. To address this challenge, we have developed optical replication mapping (ORM), a high-throughput single-molecule approach, and used it to map early-initiation events in human cells. The single-molecule nature of our data and a total of >2,500-fold coverage of the human genome on 27 million fibers averaging ∼300 kb in length allow us to identify initiation sites and their firing probability with high confidence. We find that the distribution of human replication initiation is consistent with inefficient, stochastic activation of heterogeneously distributed potential initiation complexes enriched in accessible chromatin. These observations are consistent with stochastic models of initiation-timing regulation and suggest that stochastic regulation of replication kinetics is a fundamental feature of eukaryotic replication, conserved from yeast to humans.

Chromatin structure restricts origin utilization when quiescent cells re-enter the cell cycle

Quiescent cells reside in G0 phase, which is characterized by the absence of cell growth and proliferation. These cells remain viable and re-enter the cell cycle when prompted by appropriate signals. Using a budding yeast model of cellular quiescence, we investigated the program that initiated DNA replication when these G0 cells resumed growth. Quiescent cells contained very low levels of replication initiation factors, and their entry into S phase was delayed until these factors were re-synthesized. A longer S phase in these cells correlated with the activation of fewer origins of replication compared to G1 cells. The chromatin structure around inactive origins in G0 cells showed increased H3 occupancy and decreased nucleosome positioning compared to the same origins in G1 cells, inhibiting the origin binding of the Mcm4 subunit of the MCM licensing factor. Thus, quiescent yeast cells are under-licensed during their re-entry into S phase.

High-resolution Repli-Seq defines the temporal choreography of initiation, elongation and termination of replication in mammalian cells

BACKGROUND

DNA replication in mammalian cells occurs in a defined temporal order during S phase, known as the replication timing (RT) programme. Replication timing is developmentally regulated and correlated with chromatin conformation and local transcriptional potential. Here, we present RT profiles of unprecedented temporal resolution in two human embryonic stem cell lines, human colon carcinoma line HCT116, and mouse embryonic stem cells and their neural progenitor derivatives.

RESULTS

Fine temporal windows revealed a remarkable degree of cell-to-cell conservation in RT, particularly at the very beginning and ends of S phase, and identified 5 temporal patterns of replication in all cell types, consistent with varying degrees of initiation efficiency. Zones of replication initiation (IZs) were detected throughout S phase and interacted in 3D space preferentially with other IZs of similar firing time. Temporal transition regions were resolved into segments of uni-directional replication punctuated at specific sites by small, inefficient IZs. Sites of convergent replication were divided into sites of termination or large constant timing regions consisting of many synchronous IZs in tandem. Developmental transitions in RT occured mainly by activating or inactivating individual IZs or occasionally by altering IZ firing time, demonstrating that IZs, rather than individual origins, are the units of developmental regulation. Finally, haplotype phasing revealed numerous regions of allele-specific and allele-independent asynchronous replication. Allele-independent asynchronous replication was correlated with the presence of previously mapped common fragile sites.

CONCLUSIONS

Altogether, these data provide a detailed temporal choreography of DNA replication in mammalian cells.

On the Interplay of the DNA Replication Program and the Intra-S Phase Checkpoint Pathway

DNA replication in eukaryotes is achieved by the activation of multiple replication origins which needs to be precisely coordinated in space and time. This spatio-temporal replication program is regulated by many factors to maintain genome stability, which is frequently threatened through stresses of exogenous or endogenous origin. Intra-S phase checkpoints monitor the integrity of DNA synthesis and are activated when replication forks are stalled. Their activation leads to the stabilization of forks, to the delay of the replication program by the inhibition of late firing origins, and the delay of G2/M phase entry. In some cell cycles during early development these mechanisms are less efficient in order to allow rapid cell divisions. In this article, we will review our current knowledge of how the intra-S phase checkpoint regulates the replication program in budding yeast and metazoan models, including early embryos with rapid S phases. We sum up current models on how the checkpoint can inhibit origin firing in some genomic regions, but allow dormant origin activation in other regions. Finally, we discuss how numerical and theoretical models can be used to connect the multiple different actors into a global process and to extract general rules.

The eukaryotic bell-shaped temporal rate of DNA replication origin firing emanates from a balance between origin activation and passivation

The time-dependent rate I(t) of origin firing per length of unreplicated DNA presents a universal bell shape in eukaryotes that has been interpreted as the result of a complex time-evolving interaction between origins and limiting firing factors. Here, we show that a normal diffusion of replication fork components towards localized potential replication origins (p-oris) can more simply account for the I(t) universal bell shape, as a consequence of a competition between the origin firing time and the time needed to replicate DNA separating two neighboring p-oris. We predict the I(t) maximal value to be the product of the replication fork speed with the squared p-ori density. We show that this relation is robustly observed in simulations and in experimental data for several eukaryotes. Our work underlines that fork-component recycling and potential origins localization are sufficient spatial ingredients to explain the universality of DNA replication kinetics.

Before a cell can divide, it must duplicate its DNA. In eukaryotes – organisms such as animals and fungi, which store their DNA in the cell’s nucleus – DNA replication starts at specific sites in the genome called replication origins. At each origin sits a protein complex that will activate when it randomly captures an activating protein that diffuses within the nucleus. Once a replication origin activates or “fires”, the complex then splits into two new complexes that move away from each other as they duplicate the DNA. If an active complex collides with an inactive one at another origin, the latter is inactivated – a phenomenon known as origin passivation. When two active complexes meet, they release the activating proteins, which diffuse away and eventually activate other origins in unreplicated DNA.

The number of origins that activate each minute divided by the length of unreplicated DNA is referred to as the “rate of origin firing”. In all eukaryotes, this rate – also known as I(t) – follows the same pattern. First, it increases until more than half of the DNA is duplicated. Then it decreases until everything is duplicated. This means that, if plotted out, the graph of origin firing rate would always be a bell-shaped curve, even for organisms with genomes of different sizes that have different numbers of origins. The reason for this universal shape remained unclear.

Scientists had tried to create numerical simulations that model the rate of origin firing. However, for these simulations to reproduce the bell-shape curve, a number of untested assumptions had to be made about how DNA replication takes place. In addition, these models ignored the fact that it takes time to replicate the DNA between origins.

To take this time into account, Arbona et al. instead decided to model the replication origins as discrete and distinct entities. This way of building the mathematical model succeeded in reproducing the universal bell curve shape without additional assumptions. With this simulation, the balance between origin activation and passivation is enough to achieve the observed pattern.

The new model also predicts that the maximum rate of origin firing is determined by the speed of DNA replication and the density of origins in the genome. Arbona et al. verified this prediction in yeast, fly, frog and human cells – organisms with different sized genomes that take between 20 minutes and 8 hours to replicate their DNA. Lastly, the prediction also held true in yeast treated with hydroxyurea, an anticancer drug that slows DNA replication.

A better understanding of DNA replication can help scientists to understand how this process is perturbed in cancers and how drugs that target DNA replication can treat these diseases. Future work will explore how the 3D organization of the genome affects the diffusion of activating proteins within the cell nucleus.

Replication fork slowing and stalling are distinct, checkpoint-independent consequences of replicating damaged DNA

In response to DNA damage during S phase, cells slow DNA replication. This slowing is orchestrated by the intra-S checkpoint and involves inhibition of origin firing and reduction of replication fork speed. Slowing of replication allows for tolerance of DNA damage and suppresses genomic instability. Although the mechanisms of origin inhibition by the intra-S checkpoint are understood, major questions remain about how the checkpoint regulates replication forks: Does the checkpoint regulate the rate of fork progression? Does the checkpoint affect all forks, or only those encountering damage? Does the checkpoint facilitate the replication of polymerase-blocking lesions? To address these questions, we have analyzed the checkpoint in the fission yeast Schizosaccharomyces pombe using a single-molecule DNA combing assay, which allows us to unambiguously separate the contribution of origin and fork regulation towards replication slowing, and allows us to investigate the behavior of individual forks. Moreover, we have interrogated the role of forks interacting with individual sites of damage by using three damaging agents-MMS, 4NQO and bleomycin-that cause similar levels of replication slowing with very different frequency of DNA lesions. We find that the checkpoint slows replication by inhibiting origin firing, but not by decreasing fork rates. However, the checkpoint appears to facilitate replication of damaged templates, allowing forks to more quickly pass lesions. Finally, using a novel analytic approach, we rigorously identify fork stalling events in our combing data and show that they play a previously unappreciated role in shaping replication kinetics in response to DNA damage.

Model-based analysis of DNA replication profiles: predicting replication fork velocity and initiation rate by profiling free-cycling cells

Eukaryotic cells initiate DNA synthesis by sequential firing of hundreds of origins. This ordered replication is described by replication profiles, which measure the DNA content within a cell population. Here, we show that replication dynamics can be deduced from replication profiles of free-cycling cells. While such profiles lack explicit temporal information, they are sensitive to fork velocity and initiation capacity through the passive replication pattern, namely the replication of origins by forks emanating elsewhere. We apply our model-based approach to a compendium of profiles that include most viable budding yeast mutants implicated in replication. Predicted changes in fork velocity or initiation capacity are verified by profiling synchronously replicating cells. Notably, most mutants implicated in late (or early) origin effects are explained by global modulation of fork velocity or initiation capacity. Our approach provides a rigorous framework for analyzing DNA replication profiles of free-cycling cells.

How MCM loading and spreading specify eukaryotic DNA replication initiation sites

DNA replication origins strikingly differ between eukaryotic species and cell types. Origins are localized and can be highly efficient in budding yeast, are randomly located in early fly and frog embryos, which do not transcribe their genomes, and are clustered in broad (10-100 kb) non-transcribed zones, frequently abutting transcribed genes, in mammalian cells. Nonetheless, in all cases, origins are established during the G1-phase of the cell cycle by the loading of double hexamers of the Mcm 2-7 proteins (MCM DHs), the core of the replicative helicase. MCM DH activation in S-phase leads to origin unwinding, polymerase recruitment, and initiation of bidirectional DNA synthesis. Although MCM DHs are initially loaded at sites defined by the binding of the origin recognition complex (ORC), they ultimately bind chromatin in much greater numbers than ORC and only a fraction are activated in any one S-phase. Data suggest that the multiplicity and functional redundancy of MCM DHs provide robustness to the replication process and affect replication time and that MCM DHs can slide along the DNA and spread over large distances around the ORC. Recent studies further show that MCM DHs are displaced along the DNA by collision with transcription complexes but remain functional for initiation after displacement. Therefore, eukaryotic DNA replication relies on intrinsically mobile and flexible origins, a strategy fundamentally different from bacteria but conserved from yeast to human. These properties of MCM DHs likely contribute to the establishment of broad, intergenic replication initiation zones in higher eukaryotes.

3D replicon distributions arise from stochastic initiation and domino-like DNA replication progression

DNA replication dynamics in cells from higher eukaryotes follows very complex but highly efficient mechanisms. However, the principles behind initiation of potential replication origins and emergence of typical patterns of nuclear replication sites remain unclear. Here, we propose a comprehensive model of DNA replication in human cells that is based on stochastic, proximity-induced replication initiation. Critical model features are: spontaneous stochastic firing of individual origins in euchromatin and facultative heterochromatin, inhibition of firing at distances below the size of chromatin loops and a domino-like effect by which replication forks induce firing of nearby origins. The model reproduces the empirical temporal and chromatin-related properties of DNA replication in human cells. We advance the one-dimensional DNA replication model to a spatial model by taking into account chromatin folding in the nucleus, and we are able to reproduce the spatial and temporal characteristics of the replication foci distribution throughout S-phase.

Deciphering DNA replication dynamics in eukaryotic cell populations in relation with their averaged chromatin conformations

We propose a non-local model of DNA replication that takes into account the observed uncertainty on the position and time of replication initiation in eukaryote cell populations. By picturing replication initiation as a two-state system and considering all possible transition configurations, and by taking into account the chromatin’s fractal dimension, we derive an analytical expression for the rate of replication initiation. This model predicts with no free parameter the temporal profiles of initiation rate, replication fork density and fraction of replicated DNA, in quantitative agreement with corresponding experimental data from both S. cerevisiae and human cells and provides a quantitative estimate of initiation site redundancy. This study shows that, to a large extent, the program that regulates the dynamics of eukaryotic DNA replication is a collective phenomenon that emerges from the stochastic nature of replication origins initiation.

Replication landscape of the human genome

Despite intense investigation, human replication origins and termini remain elusive. Existing data have shown strong discrepancies. Here we sequenced highly purified Okazaki fragments from two cell types and, for the first time, quantitated replication fork directionality and delineated initiation and termination zones genome-wide. Replication initiates stochastically, primarily within non-transcribed, broad (up to 150 kb) zones that often abut transcribed genes, and terminates dispersively between them. Replication fork progression is significantly co-oriented with the transcription. Initiation and termination zones are frequently contiguous, sometimes separated by regions of unidirectional replication. Initiation zones are enriched in open chromatin and enhancer marks, even when not flanked by genes, and often border ‘topologically associating domains' (TADs). Initiation zones are enriched in origin recognition complex (ORC)-binding sites and better align to origins previously mapped using bubble-trap than λ-exonuclease. This novel panorama of replication reveals how chromatin and transcription modulate the initiation process to create cell-type-specific replication programs.

The physical origin and termination sites of DNA replication in human cells have remained elusive. Here the authors use Okazaki fragment sequencing to reveal global replication patterns and show how chromatin and transcription modulate the process.

Tight Chk1 Levels Control Replication Cluster Activation in Xenopus

DNA replication in higher eukaryotes initiates at thousands of origins according to a spatio-temporal program. The ATR/Chk1 dependent replication checkpoint inhibits the activation of later firing origins. In the Xenopus in vitro system initiations are not sequence dependent and 2-5 origins are grouped in clusters that fire at different times despite a very short S phase. We have shown that the temporal program is stochastic at the level of single origins and replication clusters. It is unclear how the replication checkpoint inhibits late origins but permits origin activation in early clusters. Here, we analyze the role of Chk1 in the replication program in sperm nuclei replicating in Xenopus egg extracts by a combination of experimental and modelling approaches. After Chk1 inhibition or immunodepletion, we observed an increase of the replication extent and fork density in the presence or absence of external stress. However, overexpression of Chk1 in the absence of external replication stress inhibited DNA replication by decreasing fork densities due to lower Cdk2 kinase activity. Thus, Chk1 levels need to be tightly controlled in order to properly regulate the replication program even during normal S phase. DNA combing experiments showed that Chk1 inhibits origins outside, but not inside, already active clusters. Numerical simulations of initiation frequencies in the absence and presence of Chk1 activity are consistent with a global inhibition of origins by Chk1 at the level of clusters but need to be combined with a local repression of Chk1 action close to activated origins to fit our data.

The spatial and temporal organization of origin firing during the S-phase of fission yeast

Eukaryotes duplicate their genomes using multiple replication origins, but the organization of origin firing along chromosomes and during S-phase is not well understood. Using fission yeast, we report the first genome-wide analysis of the spatial and temporal organization of replication origin firing, analyzed using single DNA molecules that can approach the full length of chromosomes. At S-phase onset, origins fire randomly and sparsely throughout the chromosomes. Later in S-phase, clusters of fired origins appear embedded in the sparser regions, which form the basis of nuclear replication foci. The formation of clusters requires proper histone methylation and acetylation, and their locations are not inherited between cell cycles. The rate of origin firing increases gradually, peaking just before mid S-phase. Toward the end of S-phase, nearly all the available origins within the unreplicated regions are fired, contributing to the timely completion of genome replication. We propose that the majority of origins do not fire as a part of a deterministic program. Instead, origin firing, both individually and as clusters, should be viewed as being mostly stochastic.

Embryonic stem cell specific "master" replication origins at the heart of the loss of pluripotency

Epigenetic regulation of the replication program during mammalian cell differentiation remains poorly understood. We performed an integrative analysis of eleven genome-wide epigenetic profiles at 100 kb resolution of Mean Replication Timing (MRT) data in six human cell lines. Compared to the organization in four chromatin states shared by the five somatic cell lines, embryonic stem cell (ESC) line H1 displays (i) a gene-poor but highly dynamic chromatin state (EC4) associated to histone variant H2AZ rather than a HP1-associated heterochromatin state (C4) and (ii) a mid-S accessible chromatin state with bivalent gene marks instead of a polycomb-repressed heterochromatin state. Plastic MRT regions (≲ 20% of the genome) are predominantly localized at the borders of U-shaped timing domains. Whereas somatic-specific U-domain borders are gene-dense GC-rich regions, 31.6% of H1-specific U-domain borders are early EC4 regions enriched in pluripotency transcription factors NANOG and OCT4 despite being GC poor and gene deserts. Silencing of these ESC-specific “master” replication initiation zones during differentiation corresponds to a loss of H2AZ and an enrichment in H3K9me3 mark characteristic of late replicating C4 heterochromatin. These results shed a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and lineage commitment.

During development, embryonic stem cell (ESC) enter a program of cell differentiation eventually leading to all the necessary differentiated cell types. Understanding the mechanisms responsible for the underlying modifications of the gene expression program is of fundamental importance, as it will likely have strong impact on the development of regenerative medicine. We show that besides some epigenetic regulation, ubiquitous master replication origins at replication timing U-domain borders shared by 6 human cell types are transcriptionally active open chromatin regions specified by a local enrichment in nucleosome free regions encoded in the DNA sequence suggesting that they have been selected during evolution. In contrast, ESC specific master replication origins bear a unique epigenetic signature (enrichment in CTCF, H2AZ, NANOG, OCT4, …) likely contributing to maintain ESC chromatin in a highly dynamic and accessible state that is refractory to polycomb and HP1 heterochromatin spreading. These ESC specific master origins thus appear as key genomic regions where epigenetic control of chromatin organization is at play to maintain pluripotency of stem cell lineages and to guide lineage commitment to somatic cell types.

Possible role of H1 histone in replication timing

AT-rich repetitive DNA sequences become late replicating during cell differentiation. Replication timing is not correlated with LINE density in human cells (Ryba et al. 2010). However, short and properly spaced runs of oligo dA or dT present in nuclear matrix attachment regions (MARs) of the genome are good candidates for elements of AT-rich repetitive late replicating DNA. MAR attachment to the nuclear matrix is negatively regulated by chromatin binding of H1 histone, but this is counteracted by H1 phosphorylation, high mobility group proteins or, indirectly, core histone acetylation. Fewer MAR attachments correlates positively with longer average DNA loop size, longer replicons and an increase of late replicating DNA.

Large replication skew domains delimit GC-poor gene deserts in human

Besides their large-scale organization in isochores, mammalian genomes display megabase-sized regions, spanning both genes and intergenes, where the strand nucleotide composition asymmetry decreases linearly, possibly due to replication activity. These so-called skew-N domains cover about a third of the human genome and are bordered by two skew upward jumps that were hypothesized to compose a subset of "master" replication origins active in the germline. Skew-N domains were shown to exhibit a particular gene organization. Genes with CpG-rich promoters likely expressed in the germline are over represented near the master replication origins, with large genes being co-oriented with replication fork progression, which suggests some coordination of replication and transcription. In this study, we describe another skew structure that covers ∼13% of the human genome and that is bordered by putative master replication origins similar to the ones flanking skew-N domains. These skew-split-N domains have a shape reminiscent of a N, but split in half, leaving in the center a region of null skew whose length increases with domain size. These central regions (median size ∼860 kb) have a homogeneous composition, i.e. both a null and constant skew and a constant and low GC content. They correspond to heterochromatin gene deserts found in low-GC isochores with an average gene density of 0.81 promoters/Mb as compared to 7.73 promoters/Mb genome wide. The analysis of epigenetic marks and replication timing data confirms that, in these late replicating heterochomatic regions, the initiation of replication is likely to be random. This contrasts with the transcriptionally active euchromatin state found around the bordering well positioned master replication origins. Altogether skew-N domains and skew-split-N domains cover about 50% of the human genome.

DNA replication timing: Coordinating genome stability with genome regulation on the X chromosome and beyond

Recent studies based on next-generation DNA sequencing have revealed that the female inactive X chromosome is replicated in a rapid, unorganized manner, and undergoes increased rates of mutation. These observations link the organization of DNA replication timing to gene regulation on one hand, and to the generation of mutations on the other hand. More generally, the exceptional biology of the inactive X chromosome highlights general principles of genome replication. Cells may control replication timing by a combination of intrinsic replication origin properties, local chromatin states and global levels of replication factors, leading to a functional separation between the activity of genes and their mutation.

A chromatin structure-based model accurately predicts DNA replication timing in human cells

The metazoan genome is replicated in precise cell lineage‐specific temporal order. However, the mechanism controlling this orchestrated process is poorly understood as no molecular mechanisms have been identified that actively regulate the firing sequence of genome replication. Here, we develop a mechanistic model of genome replication capable of predicting, with accuracy rivaling experimental repeats, observed empirical replication timing program in humans. In our model, replication is initiated in an uncoordinated (time‐stochastic) manner at well‐defined sites. The model contains, in addition to the choice of the genomic landmark that localizes initiation, only a single adjustable parameter of direct biological relevance: the number of replication forks. We find that DNase‐hypersensitive sites are optimal and independent determinants of DNA replication initiation. We demonstrate that the DNA replication timing program in human cells is a robust emergent phenomenon that, by its very nature, does not require a regulatory mechanism determining a proper replication initiation firing sequence.

Human genome replication proceeds through four chromatin states

Advances in genomic studies have led to significant progress in understanding the epigenetically controlled interplay between chromatin structure and nuclear functions. Epigenetic modifications were shown to play a key role in transcription regulation and genome activity during development and differentiation or in response to the environment. Paradoxically, the molecular mechanisms that regulate the initiation and the maintenance of the spatio-temporal replication program in higher eukaryotes, and in particular their links to epigenetic modifications, still remain elusive. By integrative analysis of the genome-wide distributions of thirteen epigenetic marks in the human cell line K562, at the 100 kb resolution of corresponding mean replication timing (MRT) data, we identify four major groups of chromatin marks with shared features. These states have different MRT, namely from early to late replicating, replication proceeds though a transcriptionally active euchromatin state (C1), a repressive type of chromatin (C2) associated with polycomb complexes, a silent state (C3) not enriched in any available marks, and a gene poor HP1-associated heterochromatin state (C4). When mapping these chromatin states inside the megabase-sized U-domains (U-shaped MRT profile) covering about 50% of the human genome, we reveal that the associated replication fork polarity gradient corresponds to a directional path across the four chromatin states, from C1 at U-domains borders followed by C2, C3 and C4 at centers. Analysis of the other genome half is consistent with early and late replication loci occurring in separate compartments, the former correspond to gene-rich, high-GC domains of intermingled chromatin states C1 and C2, whereas the latter correspond to gene-poor, low-GC domains of alternating chromatin states C3 and C4 or long C4 domains. This new segmentation sheds a new light on the epigenetic regulation of the spatio-temporal replication program in human and provides a framework for further studies in different cell types, in both health and disease.

Previous studies revealed spatially coherent and biological-meaningful chromatin mark combinations in human cells. Here, we analyze thirteen epigenetic mark maps in the human cell line K562 at 100 kb resolution of MRT data. The complexity of epigenetic data is reduced to four chromatin states that display remarkable similarities with those reported in fly, worm and plants. These states have different MRT: (C1) is transcriptionally active, early replicating, enriched in CTCF; (C2) is Polycomb repressed, mid-S replicating; (C3) lacks of marks and replicates late and (C4) is a late-replicating gene-poor HP1 repressed heterochromatin state. When mapping these states inside the 876 replication U-domains of K562, the replication fork polarity gradient observed in these U-domains comes along with a remarkable epigenetic organization from C1 at U-domain borders to C2, C3 and ultimately C4 at centers. The remaining genome half displays early replicating, gene rich and high GC domains of intermingled C1 and C2 states segregating from late replicating, gene poor and low GC domains of concatenated C3 and/or C4 states. This constitutes the first evidence of epigenetic compartmentalization of the human genome into replication domains likely corresponding to autonomous units in the 3D chromatin architecture.

From simple bacterial and archaeal replicons to replication N/U-domains

The Replicon Theory proposed 50 years ago has proven to apply for replicons of the three domains of life. Here, we review our knowledge of genome organization into single and multiple replicons in bacteria, archaea and eukarya. Bacterial and archaeal replicator/initiator systems are quite specific and efficient, whereas eukaryotic replicons show degenerate specificity and efficiency, allowing for complex regulation of origin firing time. We expand on recent evidence that ~50% of the human genome is organized as ~1,500 megabase-sized replication domains with a characteristic parabolic (U-shaped) replication timing profile and linear (N-shaped) gradient of replication fork polarity. These N/U-domains correspond to self-interacting segments of the chromatin fiber bordered by open chromatin zones and replicate by cascades of origin firing initiating at their borders and propagating to their center, possibly by fork-stimulated initiation. The conserved occurrence of this replication pattern in the germline of mammals has resulted over evolutionary times in the formation of megabase-sized domains with an N-shaped nucleotide compositional skew profile due to replication-associated mutational asymmetries. Overall, these results reveal an evolutionarily conserved but developmentally plastic organization of replication that is driving mammalian genome evolution.

Random replication of the inactive X chromosome

In eukaryotic cells, genomic DNA replicates in a defined temporal order. The inactive X chromosome (Xi), the most extensive instance of facultative heterochromatin in mammals, replicates later than the active X chromosome (Xa), but the replication dynamics of inactive chromatin are not known. By profiling human DNA replication in an allele-specific, chromosomally phased manner, we determined for the first time the replication timing along the active and inactive chromosomes (Xa and Xi) separately. Replication of the Xi was different from that of the Xa, varied among individuals, and resembled a random, unstructured process. The Xi replicated rapidly and at a time largely separable from that of the euchromatic genome. Late-replicating, transcriptionally inactive regions on the autosomes also replicated in an unstructured manner, similar to the Xi. We conclude that DNA replication follows two strategies: slow, ordered replication associated with transcriptional activity, and rapid, random replication of silent chromatin. The two strategies coexist in the same cell, yet are segregated in space and time.

A pool of peptides extracted from wheat bud chromatin inhibits tumor cell growth by causing defective DNA synthesis

Background

We previously reported that a pool of low molecular weight peptides can be extracted by alkali treatment of DNA preparations obtained from prokaryotic and eukaryotic cells after intensive deproteinization. This class of peptides, isolated from wheat bud chromatin, induces growth inhibition, DNA damage, G2 checkpoint activation and apoptosis in HeLa cells. In this work we studied their mechanism of action by investigating their ability to interfere with DNA synthesis.

Methods

BrdUrd comet assays were used to detect DNA replication defects during S phase. DNA synthesis, cell proliferation, cell cycle progression and DNA damage response pathway activation were assessed using 3H-thymidine incorporation, DNA flow cytometry and Western blotting, respectively.

Results

BrdUrd labelling close to DNA strand discontinuities (comet tails) detects the number of active replicons. This number was significantly higher in treated cells (compared to controls) from entry until mid S phase, but markedly lower in late S phase, indicating the occurrence of defective DNA synthesis. In mid S phase the treated cells showed less 3H-thymidine incorporation with respect to the controls, which supports an early arrest of DNA synthesis. DNA damage response activation was also shown in both p53-defective HeLa cells and p53-proficient U2OS cells by the detection of the phosphorylated form of H2AX after peptide treatment. These events were accompanied in both cell lines by an increase in p21 levels and, in U2OS cells, of phospho-p53 (Ser15) levels. At 24 h of recovery after peptide treatment the cell cycle phase distribution was similar to that seen in controls and CDK1 kinase accumulation was not detected.

Conclusion

The data reported here show that the antiproliferative effect exhibited by these chromatin peptides results from their ability to induce genomic stress during DNA synthesis. This effect seems to be S-phase specific since surviving cells are able to progress through their normal cell cycle when the peptide fraction is removed from the culture medium. It is likely that the subsequent apoptosis is a consequence of the failed attempt of the tumour cells to repair the DNA damage induced by the peptides.

DNA replication timing

Patterns of replication within eukaryotic genomes correlate with gene expression, chromatin structure, and genome evolution. Recent advances in genome-scale mapping of replication kinetics have allowed these correlations to be explored in many species, cell types, and growth conditions, and these large data sets have allowed quantitative and computational analyses. One striking new correlation to emerge from these analyses is between replication timing and the three-dimensional structure of chromosomes. This correlation, which is significantly stronger than with any single histone modification or chromosome-binding protein, suggests that replication timing is controlled at the level of chromosomal domains. This conclusion dovetails with parallel work on the heterogeneity of origin firing and the competition between origins for limiting activators to suggest a model in which the stochastic probability of individual origin firing is modulated by chromosomal domain structure to produce patterns of replication. Whether these patterns have inherent biological functions or simply reflect higher-order genome structure is an open question.

A variable fork rate affects timing of origin firing and S phase dynamics in Saccharomyces cerevisiae

Activation (in the following referred to as firing) of replication origins is a continuous and irreversible process regulated by availability of DNA replication molecules and cyclin-dependent kinase activities, which are often altered in human cancers. The temporal, progressive origin firing throughout S phase appears as a characteristic replication profile, and computational models have been developed to describe this process. Although evidence from yeast to human indicates that a range of replication fork rates is observed experimentally in order to complete a timely S phase, those models incorporate velocities that are uniform across the genome. Taking advantage of the availability of replication profiles, chromosomal position and replication timing, here we investigated how fork rate may affect origin firing in budding yeast. Our analysis suggested that patterns of origin firing can be observed from a modulation of the fork rate that strongly correlates with origin density. Replication profiles of chromosomes with a low origin density were fitted with a variable fork rate, whereas for the ones with a high origin density a constant fork rate was appropriate. This indeed supports the previously reported correlation between inter-origin distance and fork rate changes. Intriguingly, the calculated correlation between fork rate and timing of origin firing allowed the estimation of firing efficiencies for the replication origins. This approach correctly retrieved origin efficiencies previously determined for chromosome VI and provided testable prediction for other chromosomal origins. Our results gain deeper insights into the temporal coordination of genome duplication, indicating that control of the replication fork rate is required for the timely origin firing during S phase.

Multiscale analysis of genome-wide replication timing profiles using a wavelet-based signal-processing algorithm

In this protocol, we describe the use of the LastWave open-source signal-processing command language (http://perso.ens-lyon.fr/benjamin.audit/LastWave/) for analyzing cellular DNA replication timing profiles. LastWave makes use of a multiscale, wavelet-based signal-processing algorithm that is based on a rigorous theoretical analysis linking timing profiles to fundamental features of the cell's DNA replication program, such as the average replication fork polarity and the difference between replication origin density and termination site density. We describe the flow of signal-processing operations to obtain interactive visual analyses of DNA replication timing profiles. We focus on procedures for exploring the space-scale map of apparent replication speeds to detect peaks in the replication timing profiles that represent preferential replication initiation zones, and for delimiting U-shaped domains in the replication timing profile. In comparison with the generally adopted approach that involves genome segmentation into regions of constant timing separated by timing transition regions, the present protocol enables the recognition of more complex patterns of the spatio-temporal replication program and has a broader range of applications. Completing the full procedure should not take more than 1 h, although learning the basics of the program can take a few hours and achieving full proficiency in the use of the software may take days.

Mathematical modeling of genome replication

Eukaryotic DNA replication is initiated from multiple sites on the chromosome, but little is known about the global and local regulation of replication. We present a mathematical model for the spatial dynamics of DNA replication, which offers insight into the kinetics of replication in different types of organisms. Most biological experiments involve average quantities over large cell populations (typically >10(7) cells) and therefore can mask the cell-to-cell variability present in the system. Although the model is formulated in terms of a population of cells, using mathematical analysis we show that one can obtain signatures of stochasticity in individual cells from averaged quantities. This work generalizes the result by Retkute et al. [Phys. Rev. Lett. 107, 068103 (2011)] to a broader set of parameter regimes.

Regulation of DNA replication within the immunoglobulin heavy-chain locus during B cell commitment

The temporal order of replication of mammalian chromosomes appears to be linked to their functional organization, but the process that establishes and modifies this order during cell differentiation remains largely unknown. Here, we studied how the replication of the Igh locus initiates, progresses, and terminates in bone marrow pro-B cells undergoing B cell commitment. We show that many aspects of DNA replication can be quantitatively explained by a mechanism involving the stochastic firing of origins (across the S phase and the Igh locus) and extensive variations in their firing rate (along the locus). The firing rate of origins shows a high degree of coordination across Igh domains that span tens to hundreds of kilobases, a phenomenon not observed in simple eukaryotes. Differences in domain sizes and firing rates determine the temporal order of replication. During B cell commitment, the expression of the B-cell-specific factor Pax5 sharply alters the temporal order of replication by modifying the rate of origin firing within various Igh domains (particularly those containing Pax5 binding sites). We propose that, within the Igh C_H-3′RR domain, Pax5 is responsible for both establishing and maintaining high rates of origin firing, mostly by controlling events downstream of the assembly of pre-replication complexes.

Each time a mammalian cell duplicates its genome in preparation for cell division it activates thousands of so called “DNA origins of replication.” The timely and complete duplication of the genome depends on careful orchestration of origin activation, which is modified when cells differentiate to perform a specific function. We currently lack a universally accepted model of origin regulation that can explain the replication dynamics in complex eukaryotes. Here, we studied the mouse immunoglobulin heavy-chain locus, one of the antibody-encoding portions of the genome, where origins change activity when antibody-producing B cells differentiate in the bone marrow. We show that multiple aspects of DNA replication initiation, progression, and termination can be explained mathematically by the interplay between randomly firing origins and two independent variables: the speed of progression of replication forks and the firing rate of origins along the locus. The rate of origin firing varies extensively along the locus during B cell differentiation and, thus, is a dominant factor in establishing the temporal order of replication. A differentiation factor called Pax5 can alter the temporal order of replication by modifying the rate of origin firing across various parts of the locus.

A quantitative model of the initiation of DNA replication in Saccharomyces cerevisiae predicts the effects of system perturbations

Background

Eukaryotic cell proliferation involves DNA replication, a tightly regulated process mediated by a multitude of protein factors. In budding yeast, the initiation of replication is facilitated by the heterohexameric origin recognition complex (ORC). ORC binds to specific origins of replication and then serves as a scaffold for the recruitment of other factors such as Cdt1, Cdc6, the Mcm2-7 complex, Cdc45 and the Dbf4-Cdc7 kinase complex. While many of the mechanisms controlling these associations are well documented, mathematical models are needed to explore the network’s dynamic behaviour. We have developed an ordinary differential equation-based model of the protein-protein interaction network describing replication initiation.

Results

The model was validated against quantified levels of protein factors over a range of cell cycle timepoints. Using chromatin extracts from synchronized Saccharomyces cerevisiae cell cultures, we were able to monitor the in vivo fluctuations of several of the aforementioned proteins, with additional data obtained from the literature. The model behaviour conforms to perturbation trials previously reported in the literature, and accurately predicts the results of our own knockdown experiments. Furthermore, we successfully incorporated our replication initiation model into an established model of the entire yeast cell cycle, thus providing a comprehensive description of these processes.

Conclusions

This study establishes a robust model of the processes driving DNA replication initiation. The model was validated against observed cell concentrations of the driving factors, and characterizes the interactions between factors implicated in eukaryotic DNA replication. Finally, this model can serve as a guide in efforts to generate a comprehensive model of the mammalian cell cycle in order to explore cancer-related phenotypes.

Replication timing and its emergence from stochastic processes

The temporal organization of DNA replication has puzzled cell biologists since before the mechanism of replication was understood. The realization that replication timing correlates with important features, such as transcription, chromatin structure and genome evolution, and is misregulated in cancer and aging has only deepened the fascination. Many ideas about replication timing have been proposed, but most have been short on mechanistic detail. However, recent work has begun to elucidate basic principles of replication timing. In particular, mathematical modeling of replication kinetics in several systems has shown that the reproducible replication timing patterns seen in population studies can be explained by stochastic origin firing at the single-cell level. This work suggests that replication timing need not be controlled by a hierarchical mechanism that imposes replication timing from a central regulator, but instead results from simple rules that affect individual origins.

Replication fork polarity gradients revealed by megabase-sized U-shaped replication timing domains in human cell lines

In higher eukaryotes, replication program specification in different cell types remains to be fully understood. We show for seven human cell lines that about half of the genome is divided in domains that display a characteristic U-shaped replication timing profile with early initiation zones at borders and late replication at centers. Significant overlap is observed between U-domains of different cell lines and also with germline replication domains exhibiting a N-shaped nucleotide compositional skew. From the demonstration that the average fork polarity is directly reflected by both the compositional skew and the derivative of the replication timing profile, we argue that the fact that this derivative displays a N-shape in U-domains sustains the existence of large-scale gradients of replication fork polarity in somatic and germline cells. Analysis of chromatin interaction (Hi-C) and chromatin marker data reveals that U-domains correspond to high-order chromatin structural units. We discuss possible models for replication origin activation within U/N-domains. The compartmentalization of the genome into replication U/N-domains provides new insights on the organization of the replication program in the human genome.

Evidence for sequential and increasing activation of replication origins along replication timing gradients in the human genome

Genome-wide replication timing studies have suggested that mammalian chromosomes consist of megabase-scale domains of coordinated origin firing separated by large originless transition regions. Here, we report a quantitative genome-wide analysis of DNA replication kinetics in several human cell types that contradicts this view. DNA combing in HeLa cells sorted into four temporal compartments of S phase shows that replication origins are spaced at 40 kb intervals and fire as small clusters whose synchrony increases during S phase and that replication fork velocity (mean 0.7 kb/min, maximum 2.0 kb/min) remains constant and narrowly distributed through S phase. However, multi-scale analysis of a genome-wide replication timing profile shows a broad distribution of replication timing gradients with practically no regions larger than 100 kb replicating at less than 2 kb/min. Therefore, HeLa cells lack large regions of unidirectional fork progression. Temporal transition regions are replicated by sequential activation of origins at a rate that increases during S phase and replication timing gradients are set by the delay and the spacing between successive origin firings rather than by the velocity of single forks. Activation of internal origins in a specific temporal transition region is directly demonstrated by DNA combing of the IGH locus in HeLa cells. Analysis of published origin maps in HeLa cells and published replication timing and DNA combing data in several other cell types corroborate these findings, with the interesting exception of embryonic stem cells where regions of unidirectional fork progression seem more abundant. These results can be explained if origins fire independently of each other but under the control of long-range chromatin structure, or if replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA. These findings shed a new light on the replication timing program of mammalian genomes and provide a general model for their replication kinetics.

Eukaryotic chromosomes replicate from multiple replication origins that fire at different times in S phase. The mechanisms that specify origin position and firing time and coordinate origins to ensure complete genome duplication are unclear. Previous studies proposed either that origins are arranged in temporally coordinated groups or fire independently of each other in a stochastic manner. Here, we have performed a quantitative analysis of human genome replication kinetics using a combination of DNA combing, which reveals local patterns of origin firing and replication fork progression on single DNA molecules, and massive sequencing of newly replicated DNA, which reveals the population-averaged replication timing profile of the entire genome. We show that origins are activated synchronously in large regions of uniform replication timing but more gradually in temporal transition regions and that the rate of origin firing increases as replication progresses. Large regions of unidirectional fork progression are abundant in embryonic stem cells but rare in differentiated cells. We propose a model in which replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA in a manner that explains the shape of the replication timing profile. These results provide a fundamental insight into the temporal regulation of mammalian genome replication.

Do replication forks control late origin firing in Saccharomyces cerevisiae?

Recent studies of eukaryotic DNA replication timing profiles suggest that the time-dependent rate of origin firing, I(t), has a universal shape, which ensures a reproducible replication completion time. However, measurements of I(t) are based on population averages, which may bias the shape of the I(t) because of imperfect cell synchrony and cell-to-cell variability. Here, we measure the population-averaged I(t) profile from synchronized Saccharomyces cerevisiae cells using DNA combing and we extract the single-cell I(t) profile using numerical deconvolution. The single cell I(t) and the population-averaged I(t) extracted from DNA combing and replication timing profiles are similar, indicating a genome scale invariance of the replication process, and excluding cell-to-cell variability in replication time as an explanation for the shape of I(t). The single cell I(t) correlates with fork density in wild-type cells, which is specifically loosened in late S phase in the clb5Δ mutant. A previously proposed numerical model that reproduces the wild-type I(t) profile, could also describe the clb5Δ mutant I(t) once modified to incorporate the decline in CDK activity and the looser dependency of initiation on fork density in the absence of Clb5p. Overall, these results suggest that the replication forks emanating from early fired origins facilitate origin firing in later-replicating regions.

Genetic variation and DNA replication timing, or why is there late replicating DNA?

Mutation rates vary significantly within the genome and across species. Recent studies revealed a long suspected replication-timing effect on mutation rate, but the mechanisms that regulate the increase in mutation rate as the genome is replicated remain unclear. Evidence is emerging, however, that DNA repair systems, in general, are less efficient in late replicating heterochromatic regions compared to early replicating euchromatic regions of the genome. At the same time, mutation rates in both vertebrates and invertebrates have been shown to vary with generation time (GT). GT is correlated with genome size, which suggests a possible nucleotypic effect on species-specific mutation rates. These and other observations all converge on a role for DNA replication checkpoints in modulating generation times and mutation rates during the DNA synthetic phase (S phase) of the cell cycle. The following will examine the potential role of the intra-S checkpoint in regulating cell cycle times (GT) and mutation rates in eukaryotes. This article was published online on August 5, 2011. An error was subsequently identified. This notice is included in the online and print versions to indicate that both have been corrected October 4, 2011.

Modeling genome-wide replication kinetics reveals a mechanism for regulation of replication timing

We developed analytical models of DNA replication that include probabilistic initiation of origins, fork progression, passive replication, and asynchrony.

We fit the model to budding yeast genome-wide microarray data probing the replication fraction and found that initiation times correlate with the precision of timing.

We extracted intrinsic origin properties, such as potential origin efficiency and firing-time distribution, which cannot be done using phenomenological approaches.

We propose that origin timing is controlled by stochastically activated initiators bound to origin sites rather than explicit time-measuring mechanisms.

The kinetics of DNA replication must be controlled for cells to develop properly. Although the biochemical mechanisms of origin initiations are increasingly well understood, the organization of initiation timing as a genome-wide program is still a mystery. With the advance of technology, researchers have been able to generate large amounts of data revealing aspects of replication kinetics. In particular, the use of microarrays to probe the replication fraction of budding yeast genome wide has been a successful first step towards unraveling the details of the replication program (Raghuraman et al, 2001; Alvino et al, 2007; McCune et al, 2008). On the surface, the microarray data shows apparent patterns of early and late replicating regions and seems to support the prevailing picture of eukaryotic replication—origins are positioned at defined sites and initiated at defined, preprogrammed times (Donaldson, 2005). Molecular combing, a single-molecule technique, however, showed that the initiation of origins is stochastic (Czajkowsky et al, 2008). Motivated by these conflicting viewpoints, we developed a model that is flexible enough to describe both deterministic and stochastic initiation.

We modeled origin initiation as probabilistic events. We first propose a model where each origin is allowed to have its distinct ‘firing-time distribution.' Origins that have well-determined initiation times have narrow distributions, whereas more stochastic origins have wider distributions. Similar models based on simulations have previously been proposed (Lygeros et al, 2008; Blow and Ge, 2009; de Moura et al, 2010); however, our model is novel in that it is analytic. It is much faster than simulations and allowed us, for the first time, to fit genome-wide microarray data and extract parameters that describe the replication program in unprecedented detail (Figure 2).

Our main result is this: origins that fire early, on average, have precisely defined initiation times, whereas origins that fire late, on average, do not have a well-defined initiation time and initiate throughout S phase. What kind of global controlling mechanism can account for this trend? We propose a second model where an origin is composed of multiple initiators, each of which fires independently and identically. A good candidate for the initiator is the minichromosome maintenance (MCM) complex, as it is found to be associated with origin firing and loaded in abundance (Hyrien et al, 2003). We show that the aforementioned relationship can be explained quantitatively if the earlier-firing origins have more MCM complexes. This model offers a new view of replication: controlled origin timing can emerge from stochastic firing and does not need an explicit time-measuring mechanism, a ‘clock.' This model provides a new, detailed, plausible, and testable mechanism for replication timing control.

Our models also capture the effects of passive replication, which is often neglected in phenomenological approaches (Eshaghi et al, 2007). There are two ways an origin site can be replicated. The site can be replicated by the origin binding to it but can also be passively replicated by neighboring origins. This complication makes it difficult to extract the intrinsic properties of origins. By modeling passive replication, we can separate the contribution from each origin and extract the potential efficiency of origins, i.e., the efficiency of the origin given that there is no passive replication. We found that while most origins are potentially highly efficient, their observed efficiency varies greatly. This implies that many origins, though capable of initiating, are often passively replicated and appear dormant. Such a design makes the replication process robust against replication stress such as fork stalling (Blow and Ge, 2009). If two approaching forks stall, normally dormant origins in the region, not being passively replicated, will initiate to replicate the gap.

With the advance of the microarray and molecular-combing technology, experiments have been done to probe many different types of cells, and large amounts of replication fraction data have been generated. Our model can be applied to spatiotemporally resolved replication fraction data for any organism, as the model is flexible enough to capture a wide range of replication kinetics. The analytical model is also much faster than simulation-based models. For these reasons, we believe that the model is a powerful tool for analyzing these large datasets. This work opens the possibility for understanding the replication program across species in more rigor and detail (Goldar et al, 2009).

Microarrays are powerful tools to probe genome-wide replication kinetics. The rich data sets that result contain more information than has been extracted by current methods of analysis. In this paper, we present an analytical model that incorporates probabilistic initiation of origins and passive replication. Using the model, we performed least-squares fits to a set of recently published time course microarray data on Saccharomyces cerevisiae. We extracted the distribution of firing times for each origin and found that the later an origin fires on average, the greater the variation in firing times. To explain this trend, we propose a model where earlier-firing origins have more initiator complexes loaded and a more accessible chromatin environment. The model demonstrates how initiation can be stochastic and yet occur at defined times during S phase, without an explicit timing program. Furthermore, we hypothesize that the initiators in this model correspond to loaded minichromosome maintenance complexes. This model is the first to suggest a detailed, testable, biochemically plausible mechanism for the regulation of replication timing in eukaryotes.

Arabidopsis thaliana chromosome 4 replicates in two phases that correlate with chromatin state

DNA replication programs have been studied extensively in yeast and animal systems, where they have been shown to correlate with gene expression and certain epigenetic modifications. Despite the conservation of core DNA replication proteins, little is known about replication programs in plants. We used flow cytometry and tiling microarrays to profile DNA replication of Arabidopsis thaliana chromosome 4 (chr4) during early, mid, and late S phase. Replication profiles for early and mid S phase were similar and encompassed the majority of the euchromatin. Late S phase exhibited a distinctly different profile that includes the remaining euchromatin and essentially all of the heterochromatin. Termination zones were consistent between experiments, allowing us to define 163 putative replicons on chr4 that clustered into larger domains of predominately early or late replication. Early-replicating sequences, especially the initiation zones of early replicons, displayed a pattern of epigenetic modifications specifying an open chromatin conformation. Late replicons, and the termination zones of early replicons, showed an opposite pattern. Histone H3 acetylated on lysine 56 (H3K56ac) was enriched in early replicons, as well as the initiation zones of both early and late replicons. H3K56ac was also associated with expressed genes, but this effect was local whereas replication time correlated with H3K56ac over broad regions. The similarity of the replication profiles for early and mid S phase cells indicates that replication origin activation in euchromatin is stochastic. Replicon organization in Arabidopsis is strongly influenced by epigenetic modifications to histones and DNA. The domain organization of Arabidopsis is more similar to that in Drosophila than that in mammals, which may reflect genome size and complexity. The distinct patterns of association of H3K56ac with gene expression and early replication provide evidence that H3K56ac may be associated with initiation zones and replication origins.

During growth and development, all plants and animals must replicate their DNA. This process is regulated to ensure that all sequences are completely and accurately replicated and is limited to S phase of the cell cycle. In the cell, DNA is packaged with histone proteins into chromatin, and both DNA and histones are subject to epigenetic modifications that affect chromatin state. Euchromatin and heterochromatin are chromatin states marked by epigenetic modifications specifying open and closed conformations, respectively. Using the model plant Arabidopsis thaliana, we show that the time at which a DNA sequence replicates is influenced by the epigenetic modifications to the surrounding chromatin. DNA replication occurs in two phases, with euchromatin replicating in early and mid S phase and heterochromatin replicating late. DNA replication time has been linked to gene expression in other organisms, and this is also true in Arabidopsis because more genes are active in euchromatin when compared to heterochromatin. The earliest replicating DNA sequences are associated with acetylation of histone H3 on lysine 56 (H3K56ac). H3K56ac is also abundant in active genes, but the patterns of association of H3K56ac with gene expression and DNA replication are distinct, suggesting that H3K56ac is independently linked to both processes.

Molecular analysis of the replication program in unicellular model organisms

Eukaryotes have long been reported to show temporal programs of replication, different portions of the genome being replicated at different times in S phase, with the added possibility of developmentally regulated changes in this pattern depending on species and cell type. Unicellular model organisms, primarily the budding yeast Saccharomyces cerevisiae, have been central to our current understanding of the mechanisms underlying the regulation of replication origins and the temporal program of replication in particular. But what exactly is a temporal program of replication, and how might it arise? In this article, we explore this question, drawing again on the wealth of experimental information in unicellular model organisms.

S-phase progression in mammalian cells: modelling the influence of nuclear organization

The control of DNA replication is of fundamental importance as cell proliferation demands that identical copies of the genetic material are passed to the two daughter cells that form during mitosis. These genetic copies are generated in the preceding S phase, where the entire DNA complement of the mother cell must be copied exactly once. As part of this process, it is known that different regions of mammalian genomes are replicated at specific times of a temporally defined replication programme. The key feature of this programme is that active genes in euchromatin are replicated before inactive ones in heterochromatin. This separation of S phase into periods where different classes of chromatin are duplicated is important in maintaining changes in gene expression that define individual cell types. Recent attempts to understand the structure of the S-phase timing programme have focused on the use of genome-wide strategies that inevitably use DNA isolated from large cell populations for analysis. However, this approach provides a composite view of events that occur within a population without knowledge of the cell-to-cell variability across the population. In this review, we attempt to combine information generated using genome-wide and single cell strategies in order to develop a coherent molecular understanding of S-phase progression. During this integration, we have explored how available information can be introduced into a modelling environment that best describes S-phase progression in mammalian cells.

Mathematical modelling of eukaryotic DNA replication

Eukaryotic DNA replication is a complex process. Replication starts at thousand origins that are activated at different times in S phase and terminates when converging replication forks meet. Potential origins are much more abundant than actually fire within a given S phase. The choice of replication origins and their time of activation is never exactly the same in any two cells. Individual origins show different efficiencies and different firing time probability distributions, conferring stochasticity to the DNA replication process. High-throughput microarray and sequencing techniques are providing increasingly huge datasets on the population-averaged spatiotemporal patterns of DNA replication in several organisms. On the other hand, single-molecule replication mapping techniques such as DNA combing provide unique information about cell-to-cell variability in DNA replication patterns. Mathematical modelling is required to fully comprehend the complexity of the chromosome replication process and to correctly interpret these data. Mathematical analysis and computer simulations have been recently used to model and interpret genome-wide replication data in the yeast Saccharomyces cerevisiae and Schizosaccharomyces pombe, in Xenopus egg extracts and in mammalian cells. These works reveal how stochasticity in origin usage confers robustness and reliability to the DNA replication process.

Reconciling stochastic origin firing with defined replication timing

Eukaryotic chromosomes replicate with defined timing patterns. However, the mechanism that regulates the timing of replication is unknown. In particular, there is an apparent conflict between population experiments, which show defined average replication times, and single-molecule experiments, which show that origins fire stochastically. Here, we provide a simple simulation that demonstrates that stochastic origin firing can produce defined average patterns of replication firing if two criteria are met. The first is that origins must have different relative firing probabilities, with origins that have relatively high firing probability being likely to fire in early S phase and origins with relatively low firing probability being unlikely to fire in early S phase. The second is that the firing probability of all origins must increase during S phase to ensure that origins with relatively low firing probability, which are unlikely to fire in early S phase, become likely to fire in late S phase. In addition, we propose biochemically plausible mechanisms for these criteria and point out how stochastic and defined origin firing can be experimentally distinguished in population experiments.

Mathematical modelling of whole chromosome replication

All chromosomes must be completely replicated prior to cell division, a requirement that demands the activation of a sufficient number of appropriately distributed DNA replication origins. Here we investigate how the activity of multiple origins on each chromosome is coordinated to ensure successful replication. We present a stochastic model for whole chromosome replication where the dynamics are based upon the parameters of individual origins. Using this model we demonstrate that mean replication time at any given chromosome position is determined collectively by the parameters of all origins. Combining parameter estimation with extensive simulations we show that there is a range of model parameters consistent with mean replication data, emphasising the need for caution in interpreting such data. In contrast, the replicated-fraction at time points through S phase contains more information than mean replication time data and allowed us to use our model to uniquely estimate many origin parameters. These estimated parameters enable us to make a number of predictions that showed agreement with independent experimental data, confirming that our model has predictive power. In summary, we demonstrate that a stochastic model can recapitulate experimental observations, including those that might be interpreted as deterministic such as ordered origin activation times.