DNA combing reveals intrinsic temporal disorder in the replication of yeast chromosome VI

Cdc45 is a critical effector of myc-dependent DNA replication stress

c-Myc oncogenic activity is thought to be mediated in part by its ability to generate DNA replication stress and subsequent genomic instability when deregulated. Previous studies have demonstrated a nontranscriptional role for c-Myc in regulating DNA replication. Here, we analyze the mechanisms by which c-Myc deregulation generates DNA replication stress. We find that overexpression of c-Myc alters the spatiotemporal program of replication initiation by increasing the density of early-replicating origins. We further show that c-Myc deregulation results in elevated replication-fork stalling or collapse and subsequent DNA damage. Notably, these phenotypes are independent of RNA transcription. Finally, we demonstrate that overexpression of Cdc45 recapitulates all c-Myc-induced replication and damage phenotypes and that Cdc45 and GINS function downstream of Myc.

Back to the origin: reconsidering replication, transcription, epigenetics, and cell cycle control

In bacteria, replication is a carefully orchestrated event that unfolds the same way for each bacterium and each cell division. The process of DNA replication in bacteria optimizes cell growth and coordinates high levels of simultaneous replication and transcription. In metazoans, the organization of replication is more enigmatic. The lack of a specific sequence that defines origins of replication has, until recently, severely limited our ability to define the organizing principles of DNA replication. This question is of particular importance as emerging data suggest that replication stress is an important contributor to inherited genetic damage and the genomic instability in tumors. We consider here the replication program in several different organisms including recent genome-wide analyses of replication origins in humans. We review recent studies on the role of cytosine methylation in replication origins, the role of transcriptional looping and gene gating in DNA replication, and the role of chromatin's 3-dimensional structure in DNA replication. We use these new findings to consider several questions surrounding DNA replication in metazoans: How are origins selected? What is the relationship between replication and transcription? How do checkpoints inhibit origin firing? Why are there early and late firing origins? We then discuss whether oncogenes promote cancer through a role in DNA replication and whether errors in DNA replication are important contributors to the genomic alterations and gene fusion events observed in cancer. We conclude with some important areas for future experimentation.

Quantitative, genome-wide analysis of eukaryotic replication initiation and termination

Many fundamental aspects of DNA replication, such as the exact locations where DNA synthesis is initiated and terminated, how frequently origins are used, and how fork progression is influenced by transcription, are poorly understood. Via the deep sequencing of Okazaki fragments, we comprehensively document replication fork directionality throughout the S. cerevisiae genome, which permits the systematic analysis of initiation, origin efficiency, fork progression, and termination. We show that leading-strand initiation preferentially occurs within a nucleosome-free region at replication origins. Using a strain in which late origins can be induced to fire early, we show that replication termination is a largely passive phenomenon that does not rely on cis-acting sequences or replication fork pausing. The replication profile is predominantly determined by the kinetics of origin firing, allowing us to reconstruct chromosome-wide timing profiles from an asynchronous culture.

Location, location, location: it's all in the timing for replication origins

The differential replication timing of eukaryotic replication origins has long been linked with epigenetic regulation of gene expression and more recently with genome stability and mutation rates; however, the mechanism has remained obscure. Recent studies have shed new light by identifying novel factors that determine origin timing in yeasts and mammalian cells and implicate the spatial organization of origins within nuclear territories in the mechanism. These new insights, along with recent findings that several initiation factors are limiting relative to licensed origins, support and shape an emerging model for replication timing control. The mechanisms that control the spatial organization of replication origins have potential impacts for genome regulation beyond replication.

The spatiotemporal program of replication in the genome of Lachancea kluyveri

We generated a genome-wide replication profile in the genome of Lachancea kluyveri and assessed the relationship between replication and base composition. This species diverged from Saccharomyces cerevisiae before the ancestral whole genome duplication. The genome comprises eight chromosomes among which a chromosomal arm of 1 Mb has a G + C-content much higher than the rest of the genome. We identified 252 active replication origins in L. kluyveri and found considerable divergence in origin location with S. cerevisiae and with Lachancea waltii. Although some global features of S. cerevisiae replication are conserved: Centromeres replicate early, whereas telomeres replicate late, we found that replication origins both in L. kluyveri and L. waltii do not behave as evolutionary fragile sites. In L. kluyveri, replication timing along chromosomes alternates between regions of early and late activating origins, except for the 1 Mb GC-rich chromosomal arm. This chromosomal arm contains an origin consensus motif different from other chromosomes and is replicated early during S-phase. We showed that precocious replication results from the specific absence of late firing origins in this chromosomal arm. In addition, we found a correlation between GC-content and distance from replication origins as well as a lack of replication-associated compositional skew between leading and lagging strands specifically in this GC-rich chromosomal arm. These findings suggest that the unusual base composition in the genome of L. kluyveri could be linked to replication.

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.

Spatio-temporal organization of replication in bacteria and eukaryotes (nucleoids and nuclei)

Here we discuss the spatio-temporal organization of replication in eubacteria and eukaryotes. Although there are significant differences in how replication is organized in cells that contain nuclei from those that do not, you will see that organization of replication in all organisms is principally dictated by the structured arrangement of the chromosome. We will begin with how replication is organized in eubacteria with particular emphasis on three well studied model organisms. We will then discuss spatial and temporal organization of replication in eukaryotes highlighting the similarities and differences between these two domains of life.

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.

Separation of DNA replication from the assembly of break-competent meiotic chromosomes

The meiotic cell division reduces the chromosome number from diploid to haploid to form gametes for sexual reproduction. Although much progress has been made in understanding meiotic recombination and the two meiotic divisions, the processes leading up to recombination, including the prolonged pre-meiotic S phase (meiS) and the assembly of meiotic chromosome axes, remain poorly defined. We have used genome-wide approaches in Saccharomyces cerevisiae to measure the kinetics of pre-meiotic DNA replication and to investigate the interdependencies between replication and axis formation. We found that replication initiation was delayed for a large number of origins in meiS compared to mitosis and that meiotic cells were far more sensitive to replication inhibition, most likely due to the starvation conditions required for meiotic induction. Moreover, replication initiation was delayed even in the absence of chromosome axes, indicating replication timing is independent of the process of axis assembly. Finally, we found that cells were able to install axis components and initiate recombination on unreplicated DNA. Thus, although pre-meiotic DNA replication and meiotic chromosome axis formation occur concurrently, they are not strictly coupled. The functional separation of these processes reveals a modular method of building meiotic chromosomes and predicts that any crosstalk between these modules must occur through superimposed regulatory mechanisms.

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.

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.

Replication origins and timing of temporal replication in budding yeast: how to solve the conundrum?

Similarly to metazoans, the budding yeast Saccharomyces cereviasiae replicates its genome with a defined timing. In this organism, well-defined, site-specific origins, are efficient and fire in almost every round of DNA replication. However, this strategy is neither conserved in the fission yeast Saccharomyces pombe, nor in Xenopus or Drosophila embryos, nor in higher eukaryotes, in which DNA replication initiates asynchronously throughout S phase at random sites. Temporal and spatial controls can contribute to the timing of replication such as Cdk activity, origin localization, epigenetic status or gene expression. However, a debate is going on to answer the question how individual origins are selected to fire in budding yeast. Two opposing theories were proposed: the "replicon paradigm" or "temporal program" vs. the "stochastic firing". Recent data support the temporal regulation of origin activation, clustering origins into temporal blocks of early and late replication. Contrarily, strong evidences suggest that stochastic processes acting on origins can generate the observed kinetics of replication without requiring a temporal order. In mammalian cells, a spatiotemporal model that accounts for a partially deterministic and partially stochastic order of DNA replication has been proposed. Is this strategy the solution to reconcile the conundrum of having both organized replication timing and stochastic origin firing also for budding yeast? In this review we discuss this possibility in the light of our recent study on the origin activation, suggesting that there might be a stochastic component in the temporal activation of the replication origins, especially under perturbed conditions.

Preferential localization of human origins of DNA replication at the 5'-ends of expressed genes and at evolutionarily conserved DNA sequences

BACKGROUND

Replication of mammalian genomes requires the activation of thousands of origins which are both spatially and temporally regulated by as yet unknown mechanisms. At the most fundamental level, our knowledge about the distribution pattern of origins in each of the chromosomes, among different cell types, and whether the physiological state of the cells alters this distribution is at present very limited.

METHODOLOGY/PRINCIPAL FINDINGS

We have used standard λ-exonuclease resistant nascent DNA preparations in the size range of 0.7-1.5 kb obtained from the breast cancer cell line MCF-7 hybridized to a custom tiling array containing 50-60 nt probes evenly distributed among genic and non-genic regions covering about 1% of the human genome. A similar DNA preparation was used for high-throughput DNA sequencing. Array experiments were also performed with DNA obtained from BT-474 and H520 cell lines. By determining the sites showing nascent DNA enrichment, we have localized several thousand origins of DNA replication. Our major findings are: (a) both array and DNA sequencing assay methods produced essentially the same origin distribution profile; (b) origin distribution is largely conserved (>70%) in all cell lines tested; (c) origins are enriched at the 5'ends of expressed genes and at evolutionarily conserved intergenic sequences; and (d) ChIP on chip experiments in MCF-7 showed an enrichment of H3K4Me3 and RNA Polymerase II chromatin binding sites at origins of DNA replication.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that the program for origin activation is largely conserved among different cell types. Also, our work supports recent studies connecting transcription initiation with replication, and in addition suggests that evolutionarily conserved intergenic sequences have the potential to participate in origin selection. Overall, our observations suggest that replication origin selection is a stochastic process significantly dependent upon local accessibility to replication factors.

Mrc1 marks early-firing origins and coordinates timing and efficiency of initiation in fission yeast

How early- and late-firing origins are selected on eukaryotic chromosomes is largely unknown. Here, we show that Mrc1, a conserved factor required for stabilization of stalled replication forks, selectively binds to the early-firing origins in a manner independent of Cdc45 and Hsk1 kinase in the fission yeast Schizosaccharomyces pombe. In mrc1Δ cells (and in swi1Δ cells to some extent), efficiency of firing is stimulated, and its timing is advanced selectively at those origins that are normally bound by Mrc1. In contrast, the late or inefficient origins which are not bound by Mrc1 are not activated in mrc1Δ cells. The enhanced firing and precocious Cdc45 loading at Mrc1-bound early-firing origins are not observed in a checkpoint mutant of mrc1, suggesting that non-checkpoint function is involved in maintaining the normal program of early-firing origins. We propose that prefiring binding of Mrc1 is an important marker of early-firing origins which are precociously activated by the absence of this protein.

pH- and salt-dependent molecular combing of DNA: experiments and phenomenological model

λ-DNA as well as plasmids can be successfully deposited by molecular combing on hydrophobic surfaces, for pH values ranging from 4 to 10. On polydimethylsiloxane (PDMS) substrates, the deposited DNA molecules are overstretched by about 60-100%. There is a significant influence of sodium ions (NaCl) on the surface density of the deposited DNA, with a maximum near to 100 mM NaCl for a DNA solution (28 ng µl(-1)) at pH 8. The combing process can be described by a micromechanical model including: (i) the adsorption of free moving coiled DNA at the substrate; (ii) the stretching of the coiled DNA by the preceding meniscus; (iii) the relaxation of the deposited DNA to the final length. The sticky ends of λ-DNA cause an adhesion force in the range of about 400 pN which allows a stable overstretching of the DNA by the preceding meniscus. The exposing of hidden hydrophobic bonds of the overstretched DNA leads to a stable deposition on the hydrophobic substrate. The pH-dependent density of deposited DNA as well as the observed influence of sodium ions can be explained by their screening of the negatively charged DNA backbone and sticky ends, respectively. The final DNA length can be derived from a balance of the stored elastic energy of the overstretched molecules and the energy of adhesion.

The DNA damage response pathway contributes to the stability of chromosome III derivatives lacking efficient replicators

In eukaryotic chromosomes, DNA replication initiates at multiple origins. Large inter-origin gaps arise when several adjacent origins fail to fire. Little is known about how cells cope with this situation. We created a derivative of Saccharomyces cerevisiae chromosome III lacking all efficient origins, the 5ORIΔ-ΔR fragment, as a model for chromosomes with large inter-origin gaps. We used this construct in a modified synthetic genetic array screen to identify genes whose products facilitate replication of long inter-origin gaps. Genes identified are enriched in components of the DNA damage and replication stress signaling pathways. Mrc1p is activated by replication stress and mediates transduction of the replication stress signal to downstream proteins; however, the response-defective mrc1^AQ allele did not affect 5ORIΔ-ΔR fragment maintenance, indicating that this pathway does not contribute to its stability. Deletions of genes encoding the DNA-damage-specific mediator, Rad9p, and several components shared between the two signaling pathways preferentially destabilized the 5ORIΔ-ΔR fragment, implicating the DNA damage response pathway in its maintenance. We found unexpected differences between contributions of components of the DNA damage response pathway to maintenance of ORIΔ chromosome derivatives and their contributions to DNA repair. Of the effector kinases encoded by RAD53 and CHK1, Chk1p appears to be more important in wild-type cells for reducing chromosomal instability caused by origin depletion, while Rad53p becomes important in the absence of Chk1p. In contrast, RAD53 plays a more important role than CHK1 in cell survival and replication fork stability following treatment with DNA damaging agents and hydroxyurea. Maintenance of ORIΔ chromosomes does not depend on homologous recombination. These observations suggest that a DNA-damage-independent mechanism enhances ORIΔ chromosome stability. Thus, components of the DNA damage response pathway contribute to genome stability, not simply by detecting and responding to DNA template damage, but also by facilitating replication of large inter-origin gaps.

Loss of genome integrity underlies aspects of aging and human disease. During DNA replication, two parallel signaling pathways play important roles in the maintenance of genome integrity. One pathway detects DNA damage, while the other senses replication stress. Both pathways activate responses that include arrest of cell cycle progression, giving cells time to cope with the problem. These pathways have been defined by treating cells with compounds that induce either replication stress or DNA damage, but little is known about their roles during unperturbed DNA replication. They may be important when several adjacent replication origins fail to initiate and forks from flanking origins must replicate longer regions of DNA than normal to complete replication. We have used a derivative of budding yeast chromosome III lacking all efficient replication origins to identify mutants that preferentially destabilize this chromosome fragment, which mimics a chromosome with a large inter-origin gap. We found that the DNA damage response pathway, but not the replication stress response pathway, plays an important role in maintaining this fragment. The signal recognized in this case may be replisome failure rather than forks stalled at endogenous DNA damage.

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

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.

Eukaryotic chromosome DNA replication: where, when, and how?

DNA replication is central to cell proliferation. Studies in the past six decades since the proposal of a semiconservative mode of DNA replication have confirmed the high degree of conservation of the basic machinery of DNA replication from prokaryotes to eukaryotes. However, the need for replication of a substantially longer segment of DNA in coordination with various internal and external signals in eukaryotic cells has led to more complex and versatile regulatory strategies. The replication program in higher eukaryotes is under a dynamic and plastic regulation within a single cell, or within the cell population, or during development. We review here various regulatory mechanisms that control the replication program in eukaryotes and discuss future directions in this dynamic field.

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.

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.

Defining replication origin efficiency using DNA fiber assays

The timely duplication of eukaryotic genomes depends on the coordinated activation of thousands of replication origins distributed along the chromosomes. Origin activation follows a temporal program that is imposed by the chromosomal context and is under the control of S-phase checkpoints. Although the general mechanisms regulating DNA replication are now well-understood at the level of individual origins, little is known on the coordination of thousands of initiation events at a genome-wide level. Recent studies using DNA combing and other single-molecule assays have shown that eukaryotic genomes contain a large excess of replication origins. Most of these origins remain "dormant" in normal growth conditions but are activated when fork progression is impeded. In this review, we discuss how DNA fiber technologies have changed our view of eukaryotic replication programs and how origin redundancy contributes to the maintenance of genome integrity in eukaryotic cells.

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.

Universal temporal profile of replication origin activation in eukaryotes

Although replication proteins are conserved among eukaryotes, the sequence requirements for replication initiation differ between species. In all species, however, replication origins fire asynchronously throughout S phase. The temporal program of origin firing is reproducible in cell populations but largely probabilistic at the single-cell level. The mechanisms and the significance of this program are unclear. Replication timing has been correlated with gene activity in metazoans but not in yeast. One potential role for a temporal regulation of origin firing is to minimize fluctuations in replication end time and avoid persistence of unreplicated DNA in mitosis. Here, we have extracted the population-averaged temporal profiles of replication initiation rates for S. cerevisiae, S. pombe, D. melanogaster, X. laevis and H. sapiens from genome-wide replication timing and DNA combing data. All the profiles have a strikingly similar shape, increasing during the first half of S phase then decreasing before its end. A previously proposed minimal model of stochastic initiation modulated by accumulation of a recyclable, limiting replication-fork factor and fork-promoted initiation of new origins, quantitatively described the observed profiles without requiring new implementations.

The selective pressure for timely completion of genome replication and optimal usage of replication proteins that must be imported into the cell nucleus can explain the generic shape of the profiles. We have identified a universal behavior of eukaryotic replication initiation that transcends the mechanisms of origin specification. The population-averaged efficiency of replication origin usage changes during S phase in a strikingly similar manner in a highly diverse set of eukaryotes. The quantitative model previously proposed for origin activation in X. laevis can be generalized to explain this evolutionary conservation.

A model for the spatiotemporal organization of DNA replication in Saccharomyces cerevisiae

DNA replication in eukaryotes is considered to proceed according to a precise program in which each chromosomal region is duplicated in a defined temporal order. However, recent studies reveal an intrinsic temporal disorder in the replication of yeast chromosome VI. Here we provide a model of the chromosomal duplication to study the temporal sequence of origin activation in budding yeast. The model comprises four parameters that influence the DNA replication system: the lengths of the chromosomes, the explicit chromosomal positions for all replication origins as well as their distinct initiation times and the replication fork migration rate. The designed model is able to reproduce the available experimental data in form of replication profiles. The dynamics of DNA replication was monitored during simulations of wild type and randomly perturbed replication conditions. Severe loss of origin function showed only little influence on the replication dynamics, so systematic deletions of origins (or loss of efficiency) were simulated to provide predictions to be tested experimentally. The simulations provide new insights into the complex system of DNA replication, showing that the system is robust to perturbation, and giving hints about the influence of a possible disordered firing.

Introduction to molecular combing: genomics, DNA replication, and cancer

The sequencing of the human genome inaugurated a new era in both fundamental and applied genetics. At the same time, the emergence of new technologies for probing the genome has transformed the field of pharmaco-genetics and made personalized genomic profiling and high-throughput screening of new therapeutic agents all but a matter of routine. One of these technologies, molecular combing, has served to bridge the technical gap between the examination of gross chromosomal abnormalities and sequence-specific alterations. Molecular combing provides a new perspective on the structure and dynamics of the human genome at the whole genome and sub-chromosomal levels with a resolution ranging from a few kilobases up to a megabase and more. Originally developed to study genetic rearrangements and to map genes for positional cloning, recent advances have extended the spectrum of its applications to studying the real-time dynamics of the replication of the genome. Understanding how the genome is replicated is essential for elucidating the mechanisms that both maintain genome integrity and result in the instabilities leading to human genetic disease and cancer. In the following, we will examine recent discoveries and advances due to the application of molecular combing to new areas of research in the fields of molecular cytogenetics and cancer genomics.

Determining the replication dynamics of specific gene loci by single-molecule analysis of replicated DNA

In metazoans, development and cell differentiation are known to affect various aspects of chromosomal organization at developmentally regulated gene loci (e.g., nuclear localization, locus accessibility, chromatin modifications, etc.). Recent evidence also indicates that cell differentiation can have dramatic effects on DNA replication initiation and replication fork progression, at several of these loci. Hence, origin selection and activation are no longer viewed as constitutive events, but as actively regulated processes, the deregulation of which could potentially affect gene expression and genomic stability at particular chromosomal locations. For these reasons, the study of DNA replication across large genomic regions has recently received increasing attention. Here we describe an assay based on the analysis of single DNA molecules that can be used to determine the steady-state distribution of replication forks across large genomic regions. This assay can be used to study how DNA replication initiates, progresses, pauses, and terminates at single-copy loci in mammalian cells.

The temporal program of chromosome replication: genomewide replication in clb5{Delta} Saccharomyces cerevisiae

Temporal regulation of origin activation is widely thought to explain the pattern of early- and late-replicating domains in the Saccharomyces cerevisiae genome. Recently, single-molecule analysis of replication suggested that stochastic processes acting on origins with different probabilities of activation could generate the observed kinetics of replication without requiring an underlying temporal order. To distinguish between these possibilities, we examined a clb5Delta strain, where origin firing is largely limited to the first half of S phase, to ask whether all origins nonspecifically show decreased firing (as expected for disordered firing) or if only some origins ("late" origins) are affected. Approximately half the origins in the mutant genome show delayed replication while the remainder replicate largely on time. The delayed regions can encompass hundreds of kilobases and generally correspond to regions that replicate late in wild-type cells. Kinetic analysis of replication in wild-type cells reveals broad windows of origin firing for both early and late origins. Our results are consistent with a temporal model in which origins can show some heterogeneity in both time and probability of origin firing, but clustering of temporally like origins nevertheless yields a genome that is organized into blocks showing different replication times.

A revisionist replicon model for higher eukaryotic genomes

The replicon model devised to explain replication control in bacteria has served as the guiding paradigm in the search for origins of replication in the more complex genomes of eukaryotes. In Saccharomyces cerevisiae, this model has proved to be extremely useful, leading to the identification of specific genetic elements (replicators) and the interacting initiator proteins that activate them. However, replication control in organisms ranging from Schizosaccharomyces pombe to mammals is far more fluid: only a small number of origins seem to represent classic replicators, while the majority correspond to zones of inefficient, closely spaced start sites none of which are indispensable for origin activity. In addition, it is apparent that the epigenetic state of a given sequence largely determines its ability to be used as a replication initiation site. These conclusions were arrived at over a period of three decades, and required the development of several novel replicon mapping techniques, as well as new ways of examining the chromatin architecture of any sequence of interest. Recently, methods have been elaborated for isolating all of the active origins in the genomes of higher eukaryotes en masse. Microarray analyses and more recent high-throughput sequencing technology will allow all the origins to be mapped onto the chromosomes of any organism whose genome has been sequenced. With the advent of whole-genome studies on gene expression and chromatin composition, the field is now positioned to define both the genetic and epigenetic rules that govern origin activity.

DNA replication timing is deterministic at the level of chromosomal domains but stochastic at the level of replicons in Xenopus egg extracts

Replication origins in Xenopus egg extracts are located at apparently random sequences but are activated in clusters that fire at different times during S phase under the control of ATR/ATM kinases. We investigated whether chromosomal domains and single sequences replicate at distinct times during S phase in egg extracts. Replication foci were found to progressively appear during early S phase and foci labelled early in one S phase colocalized with those labelled early in the next S phase. However, the distribution of these two early labels did not coincide between single origins or origin clusters on single DNA fibres. The 4 Mb Xenopus rDNA repeat domain was found to replicate later than the rest of the genome and to have a more nuclease-resistant chromatin structure. Replication initiated more frequently in the transcription unit than in the intergenic spacer. These results suggest for the first time that in this embryonic system, where transcription does not occur, replication timing is deterministic at the scale of large chromatin domains (1–5 Mb) but stochastic at the scale of replicons (10 kb) and replicon clusters (50–100 kb).

A dynamic stochastic model for DNA replication initiation in early embryos

Background

Eukaryotic cells seem unable to monitor replication completion during normal S phase, yet must ensure a reliable replication completion time. This is an acute problem in early Xenopus embryos since DNA replication origins are located and activated stochastically, leading to the random completion problem. DNA combing, kinetic modelling and other studies using Xenopus egg extracts have suggested that potential origins are much more abundant than actual initiation events and that the time-dependent rate of initiation, I(t), markedly increases through S phase to ensure the rapid completion of unreplicated gaps and a narrow distribution of completion times. However, the molecular mechanism that underlies this increase has remained obscure.

Methodology/Principal Findings

Using both previous and novel DNA combing data we have confirmed that I(t) increases through S phase but have also established that it progressively decreases before the end of S phase. To explore plausible biochemical scenarios that might explain these features, we have performed comparisons between numerical simulations and DNA combing data. Several simple models were tested: i) recycling of a limiting replication fork component from completed replicons; ii) time-dependent increase in origin efficiency; iii) time-dependent increase in availability of an initially limiting factor, e.g. by nuclear import. None of these potential mechanisms could on its own account for the data. We propose a model that combines time-dependent changes in availability of a replication factor and a fork-density dependent affinity of this factor for potential origins. This novel model quantitatively and robustly accounted for the observed changes in initiation rate and fork density.

Conclusions/Significance

This work provides a refined temporal profile of replication initiation rates and a robust, dynamic model that quantitatively explains replication origin usage during early embryonic S phase. These results have significant implications for the organisation of replication origins in higher eukaryotes.

Global regulation of genome duplication in eukaryotes: an overview from the epifluorescence microscope

In eukaryotes, DNA replication is initiated along each chromosome at multiple sites called replication origins. Locally, each replication origin is "licensed" or specified at the end of the M and the beginning of the G1 phases of the cell cycle. During the S phase when DNA synthesis takes place, origins are activated in stages corresponding to early and late-replicating domains. The staged and progressive activation of replication origins reflects the need to maintain a strict balance between the number of active replication forks and the rate at which DNA synthesis proceeds. This suggests that origin densities (frequency of initiation) and replication fork movement (rates of elongation) must be coregulated to guarantee the efficient and complete duplication of each subchromosomal domain. Emerging evidence supports this proposal and suggests that the ATM/ATR intra-S phase checkpoint plays an important role in the coregulation of initiation frequencies and rates of elongation. In this paper, we review recent results concerning the mechanisms governing the global regulation of DNA replication and discuss the roles these mechanisms play in maintaining genome stability during both a normal and perturbed S phase.