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

Control of DNA replication in vitro using a reversible replication barrier

A major obstacle to studying DNA replication is that it involves asynchronous and highly delocalized events. A reversible replication barrier overcomes this limitation and allows replication fork movement to be synchronized and localized, facilitating the study of replication fork function and replication coupled repair. Here we provide details on establishing a reversible replication barrier in vitro and using it to monitor different aspects of DNA replication. DNA template containing an array of lac operator (lacO) sequences is first bound to purified lac repressor (LacR). This substrate is then replicated in vitro using a biochemical replication system, which results in replication forks stalled on either side of the LacR array regardless of when or where they arise. Once replication forks are synchronized at the barrier, isopropyl-β-D-thiogalactopyranoside can be added to disrupt LacR binding so that replication forks synchronously resume synthesis. We describe how this approach can be employed to control replication fork elongation, termination, stalling and uncoupling, as well as assays that can be used to monitor these processes. We also explain how this approach can be adapted to control whether replication forks encounter a DNA lesion on the leading or lagging strand template and whether a converging fork is present. The required reagents can be prepared in 1-2 weeks and experiments using this approach are typically performed over 1-3 d. The main requirements for utilizing the LacR replication barrier are basic biochemical expertise and access to an in vitro system to study DNA replication. Investigators should also be trained in working with radioactive materials.

Single molecule MATAC-seq reveals key determinants of DNA replication origin efficiency

Stochastic origin activation gives rise to significant cell-to-cell variability in the pattern of genome replication. The molecular basis for heterogeneity in efficiency and timing of individual origins is a long-standing question. Here, we developed Methylation Accessibility of TArgeted Chromatin domain Sequencing (MATAC-Seq) to determine single-molecule chromatin accessibility of four specific genomic loci. MATAC-Seq relies on preferential modification of accessible DNA by methyltransferases combined with Nanopore-Sequencing for direct readout of methylated DNA-bases. Applying MATAC-Seq to selected early-efficient and late-inefficient yeast replication origins revealed large heterogeneity of chromatin states. Disruption of INO80 or ISW2 chromatin remodeling complexes leads to changes at individual nucleosomal positions that correlate with changes in their replication efficiency. We found a chromatin state with an accessible nucleosome-free region in combination with well-positioned +1 and +2 nucleosomes as a strong predictor for efficient origin activation. Thus, MATAC-Seq identifies the large spectrum of alternative chromatin states that co-exist on a given locus previously masked in population-based experiments and provides a mechanistic basis for origin activation heterogeneity during eukaryotic DNA replication. Consequently, our single-molecule chromatin accessibility assay will be ideal to define single-molecule heterogeneity across many fundamental biological processes such as transcription, replication, or DNA repair in vitro and ex vivo.

Where and when to start: Regulating DNA replication origin activity in eukaryotic genomes

In eukaryotic genomes, hundreds to thousands of potential start sites of DNA replication named origins are dispersed across each of the linear chromosomes. During S-phase, only a subset of origins is selected in a stochastic manner to assemble bidirectional replication forks and initiate DNA synthesis. Despite substantial progress in our understanding of this complex process, a comprehensive 'identity code' that defines origins based on specific nucleotide sequences, DNA structural features, the local chromatin environment, or 3D genome architecture is still missing. In this article, we review the genetic and epigenetic features of replication origins in yeast and metazoan chromosomes and highlight recent insights into how this flexibility in origin usage contributes to nuclear organization, cell growth, differentiation, and genome stability.

Origins of DNA replication in eukaryotes

Errors occurring during DNA replication can result in inaccurate replication, incomplete replication, or re-replication, resulting in genome instability that can lead to diseases such as cancer or disorders such as autism. A great deal of progress has been made toward understanding the entire process of DNA replication in eukaryotes, including the mechanism of initiation and its control. This review focuses on the current understanding of how the origin recognition complex (ORC) contributes to determining the location of replication initiation in the multiple chromosomes within eukaryotic cells, as well as methods for mapping the location and temporal patterning of DNA replication. Origin specification and configuration vary substantially between eukaryotic species and in some cases co-evolved with gene-silencing mechanisms. We discuss the possibility that centromeres and origins of DNA replication were originally derived from a common element and later separated during evolution.

3D chromatin connectivity underlies replication origin efficiency in mouse embryonic stem cells

In mammalian cells, chromosomal replication starts at thousands of origins at which replisomes are assembled. Replicative stress triggers additional initiation events from 'dormant' origins whose genomic distribution and regulation are not well understood. In this study, we have analyzed origin activity in mouse embryonic stem cells in the absence or presence of mild replicative stress induced by aphidicolin, a DNA polymerase inhibitor, or by deregulation of origin licensing factor CDC6. In both cases, we observe that the majority of stress-responsive origins are also active in a small fraction of the cell population in a normal S phase, and stress increases their frequency of activation. In a search for the molecular determinants of origin efficiency, we compared the genetic and epigenetic features of origins displaying different levels of activation, and integrated their genomic positions in three-dimensional chromatin interaction networks derived from high-depth Hi-C and promoter-capture Hi-C data. We report that origin efficiency is directly proportional to the proximity to transcriptional start sites and to the number of contacts established between origin-containing chromatin fragments, supporting the organization of origins in higher-level DNA replication factories.

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.

High-throughput analysis of single human cells reveals the complex nature of DNA replication timing control

DNA replication initiates from replication origins firing throughout S phase. Debate remains about whether origins are a fixed set of loci, or a loose agglomeration of potential sites used stochastically in individual cells, and about how consistent their firing time is. We develop an approach to profile DNA replication from whole-genome sequencing of thousands of single cells, which includes in silico flow cytometry, a method for discriminating replicating and non-replicating cells. Using two microfluidic platforms, we analyze up to 2437 replicating cells from a single sample. The resolution and scale of the data allow focused analysis of replication initiation sites, demonstrating that most occur in confined genomic regions. While initiation order is remarkably similar across cells, we unexpectedly identify several subtypes of initiation regions in late-replicating regions. Taken together, high throughput, high resolution sequencing of individual cells reveals previously underappreciated variability in replication initiation and progression.

Mapping replication forks, one replicon at a time

Claussin et al. (2022) present an elegant approach to replication fork mapping that combines single-molecule resolution with genome-wide coverage to provide unprecedented insight into the robust nature of DNA replication.

f = m*a: A Framework for Investigating the Regulation of Replication Timing

Stochastic models of replication timing posit that origin firing timing is regulated by origin firing probability, with early-firing origins having a high probability of firing and late-firing origins having a lower probability. However, they offer no insight into why one origin should have a higher firing probability than another. Here, a simple framework is suggested for how to approach the question by noting that the firing probability (f) must be the product of the stoichiometry of the MCM replicative helicase loaded at the origin (m) and the probability with which that MCM is activated (a). This framework emphasizes that mechanistic understanding of replication timing must focus on MCM loading and activation and can be simplified to the equation f = m*a.

Chromosomal Mcm2-7 distribution and the genome replication program in species from yeast to humans

The spatio-temporal program of genome replication across eukaryotes is thought to be driven both by the uneven loading of pre-replication complexes (pre-RCs) across the genome at the onset of S-phase, and by differences in the timing of activation of these complexes during S phase. To determine the degree to which distribution of pre-RC loading alone could account for chromosomal replication patterns, we mapped the binding sites of the Mcm2-7 helicase complex (MCM) in budding yeast, fission yeast, mouse and humans. We observed similar individual MCM double-hexamer (DH) footprints across the species, but notable differences in their distribution: Footprints in budding yeast were more sharply focused compared to the other three organisms, consistent with the relative sequence specificity of replication origins in S. cerevisiae. Nonetheless, with some clear exceptions, most notably the inactive X-chromosome, much of the fluctuation in replication timing along the chromosomes in all four organisms reflected uneven chromosomal distribution of pre-replication complexes.

Gene-rich regions of the genome tend to replicate earlier in S phase than do repetitive and other non-genic regions. This may be an evolutionary consequence of the fact that replication later in S phase is associated with higher frequencies of mutation and genome rearrangement. Replication timing along the chromosome is determined by 1) events prior to S-phase that specify the locations where DNA replication can be initiated, referred to as origin licensing; and 2) the timing of activation of these licensed origins during S-phase, referred to as origin firing. To determine the relative importance of these two mechanisms, here we identify both the binding sites and the abundance of a key component of the origin licensing machinery in budding yeast, fission yeast, mice, and humans, namely the replicative helicase complex. We discovered that, with a few notable exceptions, which include the inactive X chromosome in mammals, the program of replication timing can be largely explained simply on the basis of origin licensing. Our results support a model for replication timing that emphasizes stochastic firing of origins that have been licensed before S phase begins.

DNA molecular combing-based replication fork directionality profiling

The replication strategy of metazoan genomes is still unclear, mainly because definitive maps of replication origins are missing. High-throughput methods are based on population average and thus may exclusively identify efficient initiation sites, whereas inefficient origins go undetected. Single-molecule analyses of specific loci can detect both common and rare initiation events along the targeted regions. However, these usually concentrate on positioning individual events, which only gives an overview of the replication dynamics. Here, we computed the replication fork directionality (RFD) profiles of two large genes in different transcriptional states in chicken DT40 cells, namely untranscribed and transcribed DMD and CCSER1 expressed at WT levels or overexpressed, by aggregating hundreds of oriented replication tracks detected on individual DNA fibres stretched by molecular combing. These profiles reconstituted RFD domains composed of zones of initiation flanking a zone of termination originally observed in mammalian genomes and were highly consistent with independent population-averaging profiles generated by Okazaki fragment sequencing. Importantly, we demonstrate that inefficient origins do not appear as detectable RFD shifts, explaining why dispersed initiation has remained invisible to population-based assays. Our method can both generate quantitative profiles and identify discrete events, thereby constituting a comprehensive approach to study metazoan genome replication.

Single-Molecule Techniques to Study Chromatin

Besides the basic organization in nucleosome core particles (NCPs), eukaryotic chromatin is further packed through interactions with numerous protein complexes including transcription factors, chromatin remodeling and modifying enzymes. This nucleoprotein complex provides the template for many important biological processes, such as DNA replication, transcription, and DNA repair. Thus, to understand the molecular basis of these DNA transactions, it is critical to define individual changes of the chromatin structure at precise genomic regions where these machineries assemble and drive biological reactions. Single-molecule approaches provide the only possible solution to overcome the heterogenous nature of chromatin and monitor the behavior of individual chromatin transactions in real-time. In this review, we will give an overview of currently available single-molecule methods to obtain mechanistic insights into nucleosome positioning, histone modifications and DNA replication and transcription analysis-previously unattainable with population-based assays.

Consequences and Resolution of Transcription-Replication Conflicts

Transcription–replication conflicts occur when the two critical cellular machineries responsible for gene expression and genome duplication collide with each other on the same genomic location. Although both prokaryotic and eukaryotic cells have evolved multiple mechanisms to coordinate these processes on individual chromosomes, it is now clear that conflicts can arise due to aberrant transcription regulation and premature proliferation, leading to DNA replication stress and genomic instability. As both are considered hallmarks of aging and human diseases such as cancer, understanding the cellular consequences of conflicts is of paramount importance. In this article, we summarize our current knowledge on where and when collisions occur and how these encounters affect the genome and chromatin landscape of cells. Finally, we conclude with the different cellular pathways and multiple mechanisms that cells have put in place at conflict sites to ensure the resolution of conflicts and accurate genome duplication.

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.

Positive and Negative Regulation of DNA Replication Initiation

Genomic DNA is replicated every cell cycle by the programmed activation of replication origins at specific times and chromosomal locations. The factors that define the locations of replication origins and their typical activation times in eukaryotic cells are poorly understood. Previous studies highlighted the role of activating factors and epigenetic modifications in regulating replication initiation. Here, we review the role that repressive pathways - and their alleviation - play in establishing the genomic landscape of replication initiation. Several factors mediate this repression, in particular, factors associated with inactive chromatin. Repression can support organized, yet stochastic, replication initiation, and its absence could explain instances of rapid and random replication or re-replication.

FORK-seq: replication landscape of the Saccharomyces cerevisiae genome by nanopore sequencing

Genome replication mapping methods profile cell populations, masking cell-to-cell heterogeneity. Here, we describe FORK-seq, a nanopore sequencing method to map replication of single DNA molecules at 200-nucleotide resolution. By quantifying BrdU incorporation along pulse-chased replication intermediates from Saccharomyces cerevisiae, we orient 58,651 replication tracks reproducing population-based replication directionality profiles and map 4964 and 4485 individual initiation and termination events, respectively. Although most events cluster at known origins and fork merging zones, 9% and 18% of initiation and termination events, respectively, occur at many locations previously missed. Thus, FORK-seq reveals the full extent of cell-to-cell heterogeneity in DNA replication.

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.

Dynamics of DNA replication in a eukaryotic cell

Each genomic locus in a eukaryotic cell has a distinct average time of replication during S phase that depends on the spatial and temporal pattern of replication initiation events. Replication timing can affect genomic integrity because late replication is associated with an increased mutation rate. For most eukaryotes, the features of the genome that specify the location and timing of initiation events are unknown. To investigate these features for the fission yeast, Schizosaccharomyces pombe, we developed an integrative model to analyze large single-molecule and global genomic datasets. The model provides an accurate description of the complex dynamics of S. pombe DNA replication at high resolution. We present evidence that there are many more potential initiation sites in the S. pombe genome than previously identified and that the distribution of these sites is primarily determined by two factors: the sequence preferences of the origin recognition complex (ORC), and the interference of transcription with the assembly or stability of prereplication complexes (pre-RCs). We suggest that in addition to directly interfering with initiation, transcription has driven the evolution of the binding properties of ORC in S. pombe and other eukaryotic species to target pre-RC assembly to regions of the genome that are less likely to be transcribed.

Shining a Spotlight on DNA: Single-Molecule Methods to Visualise DNA

The ability to watch single molecules of DNA has revolutionised how we study biological transactions concerning nucleic acids. Many strategies have been developed to manipulate DNA molecules to investigate mechanical properties, dynamics and protein–DNA interactions. Imaging methods using small molecules and protein-based probes to visualise DNA have propelled our understanding of complex biochemical reactions involving DNA. This review focuses on summarising some of the methodological developments made to visualise individual DNA molecules and discusses how these probes have been used in single-molecule biophysical assays.

Noncoding transcription influences the replication initiation program through chromatin regulation

In eukaryotic organisms, replication initiation follows a temporal program. Among the parameters that regulate this program in Saccharomyces cerevisiae, chromatin structure has been at the center of attention without considering the contribution of transcription. Here, we revisit the replication initiation program in the light of widespread genomic noncoding transcription. We find that noncoding RNA transcription termination in the vicinity of autonomously replicating sequences (ARSs) shields replication initiation from transcriptional readthrough. Consistently, high natural nascent transcription correlates with low ARS efficiency and late replication timing. High readthrough transcription is also linked to increased nucleosome occupancy and high levels of H3K36me3. Moreover, forcing ARS readthrough transcription promotes these chromatin features. Finally, replication initiation defects induced by increased transcriptional readthrough are partially rescued in the absence of H3K36 methylation. Altogether, these observations indicate that natural noncoding transcription into ARSs influences replication initiation through chromatin regulation.

The organization of genome duplication is a critical determinant of the landscape of genome maintenance

Genome duplication is essential for cell proliferation, and the mechanisms regulating its execution are highly conserved. These processes give rise to a spatiotemporal organization of replication initiation across the genome, referred to as the replication program. Despite the identification of such programs in diverse eukaryotic organisms, their biological importance for cellular physiology remains largely unexplored. We address this fundamental question in the context of genome maintenance, taking advantage of the inappropriate origin firing that occurs when fission yeast cells lacking the Rad3/ATR checkpoint kinase are subjected to replication stress. Using this model, we demonstrate that the replication program quantitatively dictates the extent of origin de-regulation and the clustered localization of these events. Furthermore, our results uncover an accumulation of abnormal levels of single-stranded DNA (ssDNA) and the Rad52 repair protein at de-regulated origins. We show that these loci constitute a defining source of the overall ssDNA and Rad52 hotspots in the genome, generating a signature pattern of instability along the chromosomes. We then induce a genome-wide reprogramming of origin usage and evaluate its consequences in our experimental system. This leads to a complete redistribution of the sites of both inappropriate initiation and associated Rad52 recruitment. We therefore conclude that the organization of genome duplication governs the checkpoint control of origin-associated hotspots of instability and plays an integral role in shaping the landscape of genome maintenance.

The evolution of the temporal program of genome replication

Genome replication is highly regulated in time and space, but the rules governing the remodeling of these programs during evolution remain largely unknown. We generated genome-wide replication timing profiles for ten Lachancea yeasts, covering a continuous evolutionary range from closely related to more divergent species. We show that replication programs primarily evolve through a highly dynamic evolutionary renewal of the cohort of active replication origins. We found that gained origins appear with low activity yet become more efficient and fire earlier as they evolutionarily age. By contrast, origins that are lost comprise the complete range of firing strength. Additionally, they preferentially occur in close vicinity to strong origins. Interestingly, despite high evolutionary turnover, active replication origins remain regularly spaced along chromosomes in all species, suggesting that origin distribution is optimized to limit large inter-origin intervals. We propose a model on the evolutionary birth, death, and conservation of active replication origins.

Temporal programs of genome replication show different levels of conservation between closely or distantly related species. Here, the authors generate genome-wide replication timing profiles for ten yeast species, and analyze their evolutionary dynamics.

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.

Analysis of DNA Metabolism in Fission Yeast

The fission yeast Schizosaccharomyces pombe is an excellent model organism to study DNA metabolism, in which the DNA replication and repair mechanisms are evolutionarily conserved. In this introduction we describe a range of methods commonly used to study aspects of DNA metabolism in fission yeast, focusing on approaches used for the analysis of genome stability, DNA replication, and DNA repair. We describe the use of a minichromosome, Ch¹⁶, for monitoring different aspects of genome stability. We introduce two-dimensional gel electrophoresis and immunofluorescent visualization of combed DNA molecules for the analysis of DNA replication. Further, we introduce a pulsed field gel electrophoresis (PFGE) assay to physically monitor chromosome integrity, which can be used in conjunction with a DNA double-strand break (DSB) repair assay to genetically quantitate different DSB repair and misrepair outcomes, including gross chromosomal rearrangements, in fission yeast.

Analysis of Fission Yeast Single DNA Molecules on the Megabase Scale Using DNA Combing

DNA combing enables the quantitative analysis of DNA replication, DNA recombination, DNA-protein interaction, and DNA methylation along genomic single DNA molecules at 1 kb resolution. However, DNA combing has been restricted to short 200-500 kb long DNA fragments, which introduces significant bias in data analysis. An improved DNA combing methodology that allows to routinely image Mb-scale single DNA molecules and occasionally up to full-length fission yeast chromosomes is presented in this chapter. DNA combing of Mb-scale single DNA molecules can be applied to accurately measure the dynamic properties of DNA replication such as the rate of origin firing, replication fork velocity, fork directionality and the frequency of fork blockage. In addition, Mb-scale single DNA molecules enable the quantitative analysis of complex genomic rearrangements including gross chromosomal translocations, repetitive DNA sequences, large deletions, and duplications, which are difficult to investigate with deep sequencing strategies.

Replicating Large Genomes: Divide and Conquer

Complete duplication of large metazoan chromosomes requires thousands of potential initiation sites, only a small fraction of which are selected in each cell cycle. Assembly of the replication machinery is highly conserved and tightly regulated during the cell cycle, but the sites of initiation are highly flexible, and their temporal order of firing is regulated at the level of large-scale multi-replicon domains. Importantly, the number of replication forks must be quickly adjusted in response to replication stress to prevent genome instability. Here we argue that large genomes are divided into domains for exactly this reason. Once established, domain structure abrogates the need for precise initiation sites and creates a scaffold for the evolution of other chromosome functions.

Chromosome Duplication in Saccharomyces cerevisiae

The accurate and complete replication of genomic DNA is essential for all life. In eukaryotic cells, the assembly of the multi-enzyme replisomes that perform replication is divided into stages that occur at distinct phases of the cell cycle. Replicative DNA helicases are loaded around origins of DNA replication exclusively during G₁ phase. The loaded helicases are then activated during S phase and associate with the replicative DNA polymerases and other accessory proteins. The function of the resulting replisomes is monitored by checkpoint proteins that protect arrested replisomes and inhibit new initiation when replication is inhibited. The replisome also coordinates nucleosome disassembly, assembly, and the establishment of sister chromatid cohesion. Finally, when two replisomes converge they are disassembled. Studies in Saccharomyces cerevisiae have led the way in our understanding of these processes. Here, we review our increasingly molecular understanding of these events and their regulation.

Nucleosome occupancy as a novel chromatin parameter for replication origin functions

Eukaryotic DNA replication initiates from multiple discrete sites in the genome, termed origins of replication (origins). Prior to S phase, multiple origins are poised to initiate replication by recruitment of the pre-replicative complex (pre-RC). For proper replication to occur, origin activation must be tightly regulated. At the population level, each origin has a distinct firing time and frequency of activation within S phase. Many studies have shown that chromatin can strongly influence initiation of DNA replication. However, the chromatin parameters that affect properties of origins have not been thoroughly established. We found that nucleosome occupancy in G1 varies greatly around origins across the S. cerevisiae genome, and nucleosome occupancy around origins significantly correlates with the activation time and efficiency of origins, as well as pre-RC formation. We further demonstrate that nucleosome occupancy around origins in G1 is established during transition from G2/M to G1 in a pre-RC-dependent manner. Importantly, the diminished cell-cycle changes in nucleosome occupancy around origins in the orc1-161 mutant are associated with an abnormal global origin usage profile, suggesting that proper establishment of nucleosome occupancy around origins is a critical step for regulation of global origin activities. Our work thus establishes nucleosome occupancy as a novel and key chromatin parameter for proper origin regulation.

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.

Single-Molecule Analysis of Replicating Yeast Chromosomes

The faithful replication of eukaryotic chromosomal DNA occurs during S phase once per cell cycle. Replication is highly regulated and is initiated at special structures, termed origins, from which replication forks move out bidirectionally. A wide variety of techniques have been developed to study the features and kinetics of replication. Many of these, such as those based on flow cytometry and two-dimensional and pulsed-field gel electrophoresis, give a population-level view of replication. However, an alternative approach, DNA fiber analysis, which was originally developed more than 50 years ago, has the advantage of revealing features of replication at the level of individual DNA fibers. Initially based on autoradiography, this technique has been superseded by immunofluorescence-based detection of incorporated halogenated thymidine analogs. Furthermore, derivations of this technique have been developed to distribute and stretch the labeled DNA fibers uniformly on optically clear surfaces. As described here, one such technique-DNA combing, in which DNA is combed onto silanized coverslips-has been used successfully to monitor replication fork progression and origin usage in budding yeast.

Molecular Combing of Single DNA Molecules on the 10 Megabase Scale

DNA combing allows the investigation of DNA replication on genomic single DNA molecules, but the lengths that can be analysed have been restricted to molecules of 200-500 kb. We have improved the DNA combing procedure so that DNA molecules can be analysed up to the length of entire chromosomes in fission yeast and up to 12 Mb fragments in human cells. Combing multi-Mb-scale DNA molecules revealed previously undetected origin clusters in fission yeast and shows that in human cells replication origins fire stochastically forming clusters of fired origins with an average size of 370 kb. We estimate that a single human cell forms around 3200 clusters at mid S-phase and fires approximately 100,000 origins to complete genome duplication. The procedure presented here will be adaptable to other organisms and experimental conditions.

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.

Replication timing is regulated by the number of MCMs loaded at origins

Replication timing is a crucial aspect of genome regulation that is strongly correlated with chromatin structure, gene expression, DNA repair, and genome evolution. Replication timing is determined by the timing of replication origin firing, which involves activation of MCM helicase complexes loaded at replication origins. Nonetheless, how the timing of such origin firing is regulated remains mysterious. Here, we show that the number of MCMs loaded at origins regulates replication timing. We show for the first time in vivo that multiple MCMs are loaded at origins. Because early origins have more MCMs loaded, they are, on average, more likely to fire early in S phase. Our results provide a mechanistic explanation for the observed heterogeneity in origin firing and help to explain how defined replication timing profiles emerge from stochastic origin firing. These results establish a framework in which further mechanistic studies on replication timing, such as the strong effect of heterochromatin, can be pursued.

The mechanism of DNA replication termination in vertebrates

Eukaryotic DNA replication terminates when replisomes from adjacent replication origins converge. Termination involves local completion of DNA synthesis, decatenation of daughter molecules, and replisome disassembly. Termination has been difficult to study because termination events are generally asynchronous and sequence non-specific. To overcome these challenges, we paused converging replisomes with a site-specific barrier in Xenopus egg extracts. Upon removal of the barrier, forks underwent synchronous and site-specific termination, allowing mechanistic dissection of this process. We show that DNA synthesis does not slow detectably as forks approach each other and that leading strands pass each other unhindered before undergoing ligation to downstream lagging strands. Dissociation of CMG helicases occurs only after the final ligation step, and is not required for completion of DNA synthesis, strongly suggesting that converging CMGs pass one another and dissociate from double-stranded DNA. This termination mechanism allows rapid completion of DNA synthesis while avoiding premature replisome disassembly

Replication dynamics in fission and budding yeasts through DNA polymerase tracking

The dynamics of eukaryotic DNA polymerases has been difficult to establish because of the difficulty of tracking them along the chromosomes during DNA replication. Recent work has addressed this problem in the yeasts Schizosaccharomyces pombe and Saccharomyces cerevisiae through the engineering of replicative polymerases to render them prone to incorporating ribonucleotides at high rates. Their use as tracers of the passage of each polymerase has provided a picture of unprecedented resolution of the organization of replicons and replication origins in the two yeasts and has uncovered important differences between them. Additional studies have found an overlapping distribution of DNA polymorphisms and the junctions of Okazaki fragments along mononucleosomal DNA. This sequence instability is caused by the premature release of polymerase δ and the retention of non proof‐read DNA tracts replicated by polymerase α. The possible implementation of these new experimental approaches in multicellular organisms opens the door to the analysis of replication dynamics under a broad range of genetic backgrounds and physiological or pathological conditions.

Behavior of replication origins in Eukaryota - spatio-temporal dynamics of licensing and firing

Although every organism shares some common features of replication, this process varies greatly among eukaryotic species. Current data show that mathematical models of the organization of origins based on possibility theory may be applied (and remain accurate) in every model organism i.e. from yeast to humans. The major differences lie within the dynamics of origin firing and the regulation mechanisms that have evolved to meet new challenges throughout the evolution of the organism. This article elaborates on the relations between chromatin structure, organization of origins, their firing times and the impact that these features can have on genome stability, showing both differences and parallels inside the eukaryotic domain.

Volatility of Mutator Phenotypes at Single Cell Resolution

Mutator phenotypes accelerate the evolutionary process of neoplastic transformation. Historically, the measurement of mutation rates has relied on scoring the occurrence of rare mutations in target genes in large populations of cells. Averaging mutation rates over large cell populations assumes that new mutations arise at a constant rate during each cell division. If the mutation rate is not constant, an expanding mutator population may contain subclones with widely divergent rates of evolution. Here, we report mutation rate measurements of individual cell divisions of mutator yeast deficient in DNA polymerase ε proofreading and base-base mismatch repair. Our data are best fit by a model in which cells can assume one of two distinct mutator states, with mutation rates that differ by an order of magnitude. In error-prone cell divisions, mutations occurred on the same chromosome more frequently than expected by chance, often in DNA with similar predicted replication timing, consistent with a spatiotemporal dimension to the hypermutator state. Mapping of mutations onto predicted replicons revealed that mutations were enriched in the first half of the replicon as well as near termination zones. Taken together, our findings show that individual genome replication events exhibit an unexpected volatility that may deepen our understanding of the evolution of mutator-driven malignancies.

Mutations fuel microbial evolution and cancer. Cells with an increased rate of mutation are said to have a “mutator phenotype” and adapt more rapidly than non-mutator cells. Our study utilizes a novel way of measuring mutation rates of individual cell divisions to show that mutator cells can adopt one of two mutation rates that differ tenfold in magnitude. This mutator volatility suggests that the rates of mutation accumulation may vary widely within the same clone of mutator cells. Understanding how to modulate the mutator state may provide an avenue to treat certain cancers.

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.

Inferring the spatiotemporal DNA replication program from noisy data

We generalize a stochastic model of DNA replication to the case where replication-origin-initiation rates vary locally along the genome and with time. Using this generalized model, we address the inverse problem of inferring initiation rates from experimental data concerning replication in cell populations. Previous work based on curve fitting depended on arbitrarily chosen functional forms for the initiation rate, with free parameters that were constrained by the data. We introduce a nonparametric method of inference that is based on Gaussian process regression. The method replaces specific assumptions about the functional form of the initiation rate with more general prior expectations about the smoothness of variation of this rate, along the genome and in time. Using this inference method, we recover, with high precision, simulated replication schemes from noisy data that are typical of current experiments.

Elongated unique DNA strand deposition on microstructured substrate by receding meniscus assembly and capillary force

Ordered deposition of elongated DNA molecules was achieved by the forced dewetting of a DNA solution droplet over a microstructured substrate. This technique allows trapping, uncoiling, and deposition of DNA fragments without the need of a physicochemical anchoring of the molecule and results in the combing of double stranded DNA from the edge of microwells on a polydimethylsiloxane (PDMS) substrate. The technique involves scanning a droplet of DNA solution caught between a movable blade and a PDMS substrate containing an array of microwells. The deposition and elongation appears when the receding meniscus dewets microwells, the latter acting here as a perturbation in the dewetting line forcing the water film to break locally. Thus, DNA molecules can be deposited in an ordered manner and elongated conformation based solely on a physical phenomenon, allowing uncoiled DNA molecules to be observed in all their length. However, the exact mechanism that governs the deposition of DNA strands is not well understood. This paper is an analysis of the physical phenomenon occurring in the deposition process and is based on observations made with the use of high frame/second rate video microscopy.

Epigenetic landscape for initiation of DNA replication

The key genetic process of DNA replication is initiated at specific sites referred to as replication origins. In eukaryotes, origins of DNA replication are not specified by a defined nucleotide sequence. Recent studies have shown that the structural context and topology of DNA sequence, chromatin features, and its transcriptional activity play an important role in origin choice. During differentiation and development, significant changes in chromatin organization and transcription occur, influencing origin activity and choice. In the last few years, a number of different genome-wide studies have broadened the understanding of replication origin regulation. In this review, we discuss the epigenetic factors and mechanisms that modulate origin choice and firing.

A personal reflection on the replicon theory: from R1 plasmid to replication timing regulation in human cells

Fifty years after the Replicon Theory was originally presented, detailed mechanistic insight into prokaryotic replicons has been obtained and rapid progress is being made to elucidate the more complex regulatory mechanisms of replicon regulation in eukaryotic cells. Here, I present my personal perspectives on how studies of model replicons have contributed to our understanding of the basic mechanisms of DNA replication as well as the evolution of replication regulation in human cells. I will also discuss how replication regulation contributes to the stable maintenance of the genome and how disruption of replication regulation leads to human diseases.

Bubble-seq analysis of the human genome reveals distinct chromatin-mediated mechanisms for regulating early- and late-firing origins

We have devised a method for isolating virtually pure and comprehensive libraries of restriction fragments that contained replication initiation sites (bubbles) in vivo. We have now sequenced and mapped the bubble-containing fragments from GM06990, a near-normal EBV-transformed lymphoblastoid cell line, and have compared origin distributions with a comprehensive replication timing study recently published for this cell line. We find that early-firing origins, which represent ∼32% of all origins, overwhelmingly represent zones, associate only marginally with active transcription units, are localized within large domains of open chromatin, and are significantly associated with DNase I hypersensitivity. Origin "density" falls from early- to mid-S-phase, but rises again in late S-phase to levels only 17% lower than in early S-phase. Unexpectedly, late origin density calculated on the 1-Mb scale increases as a function of increasing chromatin compaction. Furthermore, the median efficiency of origins in late-replicating, heterochromatic domains is only 25% lower than in early-replicating euchromatic loci. Thus, the activation of early- and late-firing origins must be regulated by quintessentially different mechanisms. The aggregate data can be unified into a model in which initiation site selection is driven almost entirely by epigenetic factors that fashion both the long-range and local chromatin environments, with underlying DNA sequence and local transcriptional activity playing only minor roles. Importantly, the comprehensive origin map we have prepared for GM06990 overlaps moderately well with origin maps recently reported for the genomes of four different human cell lines based on the distributions of small nascent strands.

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.

Stochastic association of neighboring replicons creates replication factories in budding yeast

Inside the nucleus, DNA replication is organized at discrete sites called replication factories, consisting of DNA polymerases and other replication proteins. Replication factories play important roles in coordinating replication and in responding to replication stress. However, it remains unknown how replicons are organized for processing at each replication factory. Here we address this question using budding yeast. We analyze how individual replicons dynamically organized a replication factory using live-cell imaging and investigate how replication factories were structured using super-resolution microscopy. Surprisingly, we show that the grouping of replicons within factories is highly variable from cell to cell. Once associated, however, replicons stay together relatively stably to maintain replication factories. We derive a coherent genome-wide mathematical model showing how neighboring replicons became associated stochastically to form replication factories, which was validated by independent microscopy-based analyses. This study not only reveals the fundamental principles promoting replication factory organization in budding yeast, but also provides insight into general mechanisms by which chromosomes organize sub-nuclear structures.

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.

Next-generation sequencing of Okazaki fragments extracted from Saccharomyces cerevisiae

Genome-wide Okazaki fragment distribution can differentiate the discontinuous from the semi-discontinuous DNA replication model. Here, we investigated the genome-wide Okazaki fragment distribution in Saccharomyces cerevisiae S288C. We improved the method based upon lambda exonuclease digestion to purify Okazaki fragments from S288C yeast cells, followed by Illumina sequencing. The distribution of Okazaki fragments around confirmed replication origins, including two highly efficient replication origins, supported the discontinuous DNA replication model. In S. cerevisiae mitochondria, Okazaki fragments were overrepresented in the transcribed regions, indicating the interplay between transcription and DNA replication.