Molecular analysis of the replication program in unicellular model organisms

The yeast Dbf4 Zn2+ finger domain suppresses single-stranded DNA at replication forks initiated from a subset of origins

Dbf4 is the cyclin-like subunit for the Dbf4-dependent protein kinase (DDK), required for activating the replicative helicase at DNA replication origin that fire during S phase. Dbf4 also functions as an adaptor, targeting the DDK to different groups of origins and substrates. Here we report a genome-wide analysis of origin firing in a budding yeast mutant, dbf4-zn, lacking the Zn²⁺ finger domain within the C-terminus of Dbf4. At one group of origins, which we call dromedaries, we observe an unanticipated DNA replication phenotype: accumulation of single-stranded DNA spanning ± 5kbp from the center of the origins. A similar accumulation of single-stranded DNA at origins occurs more globally in pri1-m4 mutants defective for the catalytic subunit of DNA primase and rad53 mutants defective for the S phase checkpoint following DNA replication stress. We propose the Dbf4 Zn²⁺ finger suppresses single-stranded gaps at replication forks emanating from dromedary origins. Certain origins may impose an elevated requirement for the DDK to fully initiate DNA synthesis following origin activation. Alternatively, dbf4-zn may be defective for stabilizing/restarting replication forks emanating from dromedary origins during replication stress.

Genomic methods for measuring DNA replication dynamics

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density, and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analog incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

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.

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.

The Eukaryotic Mismatch Recognition Complexes Track with the Replisome during DNA Synthesis

During replication, mismatch repair proteins recognize and repair mispaired bases that escape the proofreading activity of DNA polymerase. In this work, we tested the model that the eukaryotic mismatch recognition complex tracks with the advancing replisome. Using yeast, we examined the dynamics during replication of the leading strand polymerase Polε using Pol2 and the eukaryotic mismatch recognition complex using Msh2, the invariant protein involved in mismatch recognition. Specifically, we synchronized cells and processed samples using chromatin immunoprecipitation combined with custom DNA tiling arrays (ChIP-chip). The Polε signal was not detectable in G1, but was observed at active origins and replicating DNA throughout S-phase. The Polε signal provided the resolution to track origin firing timing and efficiencies as well as replisome progression rates. By detecting Polε and Msh2 dynamics within the same strain, we established that the mismatch recognition complex binds origins and spreads to adjacent regions with the replisome. In mismatch repair defective PCNA mutants, we observed that Msh2 binds to regions of replicating DNA, but the distribution and dynamics are altered, suggesting that PCNA is not the sole determinant for the mismatch recognition complex association with replicating regions, but may influence the dynamics of movement. Using biochemical and genomic methods, we provide evidence that both MutS complexes are in the vicinity of the replisome to efficiently repair the entire spectrum of mutations during replication. Our data supports the model that the proximity of MutSα/β to the replisome for the efficient repair of the newly synthesized strand before chromatin reassembles.

During replication, errors that escape the replication machinery are identified and repaired by DNA mismatch repair proteins. A mismatch in the helix is recognized by MutS homologs and subsequent events include excision of the error-containing strand followed by re-synthesis. A critical step in this process is directing repair to the newly synthesized strand. Current data suggest that transient discontinuities in the DNA backbone, known as nicks, generated during replication serve as the strand discrimination signals. Additionally, proteins that package DNA have the capacity to block mismatch recognition and are known to rapidly assemble behind the replication fork. Thus, there must be a short window of opportunity for the mismatch recognition complexes to scan for mismatches and access the strand discrimination signals. To address these issues, we tested the model that the mismatch recognition complexes track with the replisome. We employed high resolution genomic methods to determine that during replication, the mismatch recognition complexes bind origins of replication and advances with the replisome. The findings support the hypothesis that the mismatch recognition proteins track with the DNA replication machinery to accurately survey and repair the newly synthesized strands while the DNA is unpackaged and strand specificity signals are accessible.

Chromosomal context and replication properties of ARS plasmids in Schizosaccharomyces pombe

Short, specific DNA sequences called as Autonomously Replicating Sequence (ARS) elements function as plasmid as well as chromosomal replication origins in yeasts. As compared to ARSs, different chromosomal origins vary greatly in their efficiency and timing of replication probably due to their wider chromosomal context. The two Schizosaccharomyces pombe ARS elements, ars727 and ars2004, represent two extremities in their chromosomal origin activity - ars727 is inactive and late replicating, while ars2004 is a highly active, early-firing origin. To determine the effect of chromosomal context on the activity of these ARS elements, we have cloned them with their extended chromosomal context as well as in the context of each other in both orientations and analysed their replication efficiency by ARS and plasmid stability assays. We found that these ARS elements retain their origin activity in their extended/altered context. However, deletion of a 133-bp region of the previously reported ars727- associated late replication enforcing element (LRE) caused advancement in replication timing of the resulting plasmid. These results confirm the role of LRE in directing plasmid replication timing and suggest that the plasmid origin efficiency of ars2004 or ars727 remains unaltered by the extended chromosomal context.

Peaks cloaked in the mist: the landscape of mammalian replication origins

Replication of mammalian genomes starts at sites termed replication origins, which historically have been difficult to locate as a result of large genome sizes, limited power of genetic identification schemes, and rareness and fragility of initiation intermediates. However, origins are now mapped by the thousands using microarrays and sequencing techniques. Independent studies show modest concordance, suggesting that mammalian origins can form at any DNA sequence but are suppressed by read-through transcription or that they can overlap the 5' end or even the entire gene. These results require a critical reevaluation of whether origins form at specific DNA elements and/or epigenetic signals or require no such determinants.

Genome architecture and its roles in human copy number variation

Besides single-nucleotide variants in the human genome, large-scale genomic variants, such as copy number variations (CNVs), are being increasingly discovered as a genetic source of human diversity and the pathogenic factors of diseases. Recent experimental findings have shed light on the links between different genome architectures and CNV mutagenesis. In this review, we summarize various genomic features and discuss their contributions to CNV formation. Genomic repeats, including both low-copy and high-copy repeats, play important roles in CNV instability, which was initially known as DNA recombination events. Furthermore, it has been found that human genomic repeats can also induce DNA replication errors and consequently result in CNV mutations. Some recent studies showed that DNA replication timing, which reflects the high-order information of genomic organization, is involved in human CNV mutations. Our review highlights that genome architecture, from DNA sequence to high-order genomic organization, is an important molecular factor in CNV mutagenesis and human genomic instability.

Checkpoint-independent scaling of the Saccharomyces cerevisiae DNA replication program

Background

In budding yeast, perturbations that prolong S phase lead to a proportionate delay in the activation times of most origins. The DNA replication checkpoint was implicated in this scaling phenotype, as an intact checkpoint was shown to be required for the delayed activation of late origins in response to hydroxyurea treatment. In support of that, scaling is lost in cells deleted of mrc1, a mediator of the replication checkpoint signal. Mrc1p, however, also plays a role in normal replication.

Results

To examine whether the replication checkpoint is required for scaling the replication profile with S phase duration we measured the genome-wide replication profile of different MRC1 alleles that separate its checkpoint function from its role in normal replication, and further analyzed the replication profiles of S phase mutants that are checkpoint deficient. We found that the checkpoint is not required for scaling; rather the unique replication phenotype of mrc1 deleted cells is attributed to the role of Mrc1 in normal replication. This is further supported by the replication profiles of tof1Δ which functions together with Mrc1p in normal replication, and by the distinct replication profiles of specific POL2 alleles which differ in their interaction with Mrc1p.

Conclusions

We suggest that the slow fork progression in mrc1 deleted cells reduces the likelihood of passive replication leading to the activation of origins that remain mostly dormant in wild-type cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-014-0079-z) contains supplementary material, which is available to authorized users.

High-resolution replication profiles define the stochastic nature of genome replication initiation and termination

Eukaryotic genome replication is stochastic, and each cell uses a different cohort of replication origins. We demonstrate that interpreting high-resolution Saccharomyces cerevisiae genome replication data with a mathematical model allows quantification of the stochastic nature of genome replication, including the efficiency of each origin and the distribution of termination events. Single-cell measurements support the inferred values for stochastic origin activation time. A strain, in which three origins were inactivated, confirmed that the distribution of termination events is primarily dictated by the stochastic activation time of origins. Cell-to-cell variability in origin activity ensures that termination events are widely distributed across virtually the whole genome. We propose that the heterogeneity in origin usage contributes to genome stability by limiting potentially deleterious events from accumulating at particular loci.

•

Deep sequencing reveals a high-resolution view of genome replication dynamics

•

Genome-wide modeling and single-cell imaging reveal stochastic origin activity

•

Origin activity determines the location of replication termination events

•

Termination events are widely distributed across the whole genome

Time to be versatile: regulation of the replication timing program in budding yeast

Eukaryotic replication origins are activated at different times during the S phase of the cell cycle, following a temporal program that is stably transmitted to daughter cells. Although the mechanisms that control initiation at the level of individual origins are now well understood, much less is known on how cells coordinate replication at hundreds of origins distributed on the chromosomes. In this review, we discuss recent advances shedding new light on how this complex process is regulated in the budding yeast Saccharomyces cerevisiae. The picture that emerges from these studies is that replication timing is regulated in cis by mechanisms modulating the chromatin structure and the subnuclear organization of origins. These mechanisms do not affect the licensing of replication origins but determine their ability to compete for limiting initiation factors, which are recycled from early to late origins throughout the length of the S phase.

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.

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.

Analysis of DNA replication profiles in budding yeast and mammalian cells using DNA combing

DNA combing is a powerful method developed by Bensimon and colleagues to stretch DNA molecules on silanized glass coverslips. This technique provides a unique way to monitor the activation of replication origins and the progression of replication forks at the level of single DNA molecules, after incorporation of thymidine analogs, such as 5-bromo-2'-deoxyuridine (BrdU), 5-iodo-2'-deoxyuridine (IdU) and 5-chloro-2'-deoxyuridine (CldU) in newly-synthesized DNA. Unlike microarray-based approaches, this assay gives access to the variability of replication profiles in individual cells. It can also be used to monitor the effect of DNA lesions on fork progression, arrest and restart. In this review, we propose standard DNA combing methods to analyze DNA replication in budding yeast and in human cells. We also show that 5-ethynyl-2'-deoxyuridine (EdU) can be used as a good alternative to BrdU for DNA combing analysis, as unlike halogenated nucleotides, it can be detected without prior denaturation of DNA.

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

Principles of chromosomal organization: lessons from yeast

The spatial organization of genes and chromosomes plays an important role in the regulation of several DNA processes. However, the principles and forces underlying this nonrandom organization are mostly unknown. Despite its small dimension, and thanks to new imaging and biochemical techniques, studies of the budding yeast nucleus have led to significant insights into chromosome arrangement and dynamics. The dynamic organization of the yeast genome during interphase argues for both the physical properties of the chromatin fiber and specific molecular interactions as drivers of nuclear order.

Genome-scale analysis of replication timing: from bench to bioinformatics

Replication timing profiles are cell type-specific and reflect genome organization changes during differentiation. In this protocol, we describe how to analyze genome-wide replication timing (RT) in mammalian cells. Asynchronously cycling cells are pulse labeled with the nucleotide analog 5-bromo-2-deoxyuridine (BrdU) and sorted into S-phase fractions on the basis of DNA content using flow cytometry. BrdU-labeled DNA from each fraction is immunoprecipitated, amplified, differentially labeled and co-hybridized to a whole-genome comparative genomic hybridization microarray, which is currently more cost effective than high-throughput sequencing and equally capable of resolving features at the biologically relevant level of tens to hundreds of kilobases. We also present a guide to analyzing the resulting data sets based on methods we use routinely. Subjects include normalization, scaling and data quality measures, LOESS (local polynomial) smoothing of RT values, segmentation of data into domains and assignment of timing values to gene promoters. Finally, we cover clustering methods and means to relate changes in the replication program to gene expression and other genetic and epigenetic data sets. Some experience with R or similar programming languages is assumed. All together, the protocol takes ∼3 weeks per batch of samples.

Initiation of DNA replication: functional and evolutionary aspects

BACKGROUND

The initiation of DNA replication is a very important and highly regulated step in the cell division cycle. It is of interest to compare different groups of eukaryotic organisms (a) to identify the essential molecular events that occur in all eukaryotes, (b) to start to identify higher-level regulatory mechanisms that are specific to particular groups and (c) to gain insights into the evolution of initiation mechanisms.

SCOPE

This review features a wide-ranging literature survey covering replication origins, origin recognition and usage, modification of origin usage (especially in response to plant hormones), assembly of the pre-replication complex, loading of the replisome, genomics, and the likely origin of these mechanisms and proteins in Archaea.

CONCLUSIONS

In all eukaryotes, chromatin is organized for DNA replication as multiple replicons. In each replicon, replication is initiated at an origin. With the exception of those in budding yeast, replication origins, including the only one to be isolated so far from a plant, do not appear to embody a specific sequence; rather, they are AT-rich, with short tracts of locally bent DNA. The proteins involved in initiation are remarkably similar across the range of eukaryotes. Nevertheless, their activity may be modified by plant-specific mechanisms, including regulation by plant hormones. The molecular features of initiation are seen in a much simpler form in the Archaea. In particular, where eukaryotes possess a number of closely related proteins that form 'hetero-complexes' (such as the origin recognition complex and the MCM complex), archaeans typically possess one type of protein (e.g. one MCM) that forms a homo-complex. This suggests that several eukaryotic initiation proteins have evolved from archaeal ancestors by gene duplication and divergence.

Quantitative proteomic analysis of chromatin reveals that Ctf18 acts in the DNA replication checkpoint

Yeast cells lacking Ctf18, the major subunit of an alternative Replication Factor C complex, have multiple problems with genome stability. To understand the in vivo function of the Ctf18 complex, we analyzed chromatin composition in a ctf18Δ mutant using the quantitative proteomic technique of stable isotope labeling by amino acids in cell culture. Three hundred and seven of the 491 reported chromosomal proteins were quantitated. The most marked abnormalities occurred when cells were challenged with the replication inhibitor hydroxyurea. Compared with wild type, hydroxyurea-treated ctf18Δ cells exhibited increased chromatin association of replisome progression complex components including Cdc45, Ctf4, and GINS complex subunits, the polymerase processivity clamp PCNA and the single-stranded DNA-binding complex RPA. Chromatin composition abnormalities observed in ctf18Δ cells were very similar to those of an mrc1Δ mutant, which is defective in the activating the Rad53 checkpoint kinase in response to DNA replication stress. We found that ctf18Δ cells are also defective in Rad53 activation, revealing that the Ctf18 complex is required for engagement of the DNA replication checkpoint. Inappropriate initiation of replication at late origins, because of loss of the checkpoint, probably causes the elevated level of chromatin-bound replisome proteins in the ctf18Δ mutant. The role of Ctf18 in checkpoint activation is not shared by all Replication Factor C-like complexes, because proteomic analysis revealed that cells lacking Elg1 (the major subunit of a different Replication Factor C-like complex) display a different spectrum of chromatin abnormalities. Identification of Ctf18 as a checkpoint protein highlights the usefulness of chromatin proteomic analysis for understanding the in vivo function of proteins that mediate chromatin transactions.

Analysis of replication profiles reveals key role of RFC-Ctf18 in yeast replication stress response

Maintenance of genome integrity relies on surveillance mechanisms that detect and signal arrested replication forks. Although evidence from budding yeast indicates that the DNA replication checkpoint (DRC) is primarily activated by single-stranded DNA (ssDNA), studies in higher eukaryotes have implicated primer ends in this process. To identify factors that signal primed ssDNA in Saccharomyces cerevisiae, we have screened a collection of checkpoint mutants for their ability to activate the DRC, using the repression of late origins as readout for checkpoint activity. This quantitative analysis reveals that neither RFC(Rad24) and the 9-1-1 clamp nor the alternative clamp loader RFC(Elg1) is required to signal paused forks. In contrast, we found that RFC(Ctf18) is essential for the Mrc1-dependent activation of Rad53 and for the maintenance of paused forks. These data identify RFC(Ctf18) as a key DRC mediator, potentially bridging Mrc1 and primed ssDNA to signal paused forks.

Evaluating genome-scale approaches to eukaryotic DNA replication

Mechanisms regulating where and when eukaryotic DNA replication initiates remain a mystery. Recently, genome-scale methods have been brought to bear on this problem. The identification of replication origins and their associated proteins in yeasts is a well-integrated investigative tool, but corresponding data sets from multicellular organisms are scarce. By contrast, standardized protocols for evaluating replication timing have generated informative data sets for most eukaryotic systems. Here, I summarize the genome-scale methods that are most frequently used to analyse replication in eukaryotes, the kinds of questions each method can address and the technical hurdles that must be overcome to gain a complete understanding of the nature of eukaryotic replication origins.