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

Chromatin as a Coevolutionary Graph: Modeling the Interplay of Replication with Chromatin Dynamics

Modeling DNA replication poses significant challenges due to the intricate interplay of biophysical processes and the need for precise parameter optimization. In this study, we explore the interactions among three key biophysical factors that influence chromatin folding: replication, loop extrusion, and compartmentalization. Replication forks, known to act as barriers to the motion of loop extrusion factors, also correlate with the phase separation of chromatin into A and B compartments. Our approach integrates three components: (1) a numerical model that takes into advantage single-cell replication timing data to simulate replication fork propagation; (2) a stochastic Monte Carlo simulation that captures the interplay between the biophysical factors, with loop extrusion factors binding, unbinding, and extruding dynamically, while CTCF barriers and replication forks act as static and moving barriers, and a Potts Hamiltonian governs the spreading of epigenetic states driving chromatin compartmentalization; and (3) a 3D OpenMM simulation that reconstructs the chromatin's 3D structure based on the states generated by the stochastic model. To our knowledge, this is the first framework to dynamically integrate and simulate these three biophysical factors, enabling insights into chromatin behavior during replication. Furthermore, we investigate how replication stress alters these dynamics and affects chromatin structure.

DNA replication and replication stress response in the context of nuclear architecture

The DNA replication process needs to be coordinated with other DNA metabolism transactions and must eventually extend to the full genome, regardless of chromatin status, gene expression, secondary structures and DNA lesions. Completeness and accuracy of DNA replication are crucial to maintain genome integrity, limiting transformation in normal cells and offering targeting opportunities for proliferating cancer cells. DNA replication is thus tightly coordinated with chromatin dynamics and 3D genome architecture, and we are only beginning to understand the underlying molecular mechanisms. While much has recently been discovered on how DNA replication initiation is organised and modulated in different genomic regions and nuclear territories-the so-called "DNA replication program"-we know much less on how the elongation of ongoing replication forks and particularly the response to replication obstacles is affected by the local nuclear organisation. Also, it is still elusive how specific components of nuclear architecture participate in the replication stress response. Here, we review known mechanisms and factors orchestrating replication initiation, and replication fork progression upon stress, focusing on recent evidence linking genome organisation and nuclear architecture with the cellular responses to replication interference, and highlighting open questions and future challenges to explore this exciting new avenue of research.

Unveiling human origins of replication using deep learning: accurate prediction and comprehensive analysis

Accurate identification of replication origins (ORIs) is crucial for a comprehensive investigation into the progression of human cell growth and cancer therapy. Here, we proposed a computational approach Ori-FinderH, which can efficiently and precisely predict the human ORIs of various lengths by combining the Z-curve method with deep learning approach. Compared with existing methods, Ori-FinderH exhibits superior performance, achieving an area under the receiver operating characteristic curve (AUC) of 0.9616 for K562 cell line in 10-fold cross-validation. In addition, we also established a cross-cell-line predictive model, which yielded a further improved AUC of 0.9706. The model was subsequently employed as a fitness function to support genetic algorithm for generating artificial ORIs. Sequence analysis through iORI-Euk revealed that a vast majority of the created sequences, specifically 98% or more, incorporate at least one ORI for three cell lines (Hela, MCF7 and K562). This innovative approach could provide more efficient, accurate and comprehensive information for experimental investigation, thereby further advancing the development of this field.

Replisome loading reduces chromatin motion independent of DNA synthesis

Chromatin has been shown to undergo diffusional motion, which is affected during gene transcription by RNA polymerase activity. However, the relationship between chromatin mobility and other genomic processes remains unclear. Hence, we set out to label the DNA directly in a sequence unbiased manner and followed labeled chromatin dynamics in interphase human cells expressing GFP-tagged proliferating cell nuclear antigen (PCNA), a cell cycle marker and core component of the DNA replication machinery. We detected decreased chromatin mobility during the S-phase compared to G1 and G2 phases in tumor as well as normal diploid cells using automated particle tracking. To gain insight into the dynamical organization of the genome during DNA replication, we determined labeled chromatin domain sizes and analyzed their motion in replicating cells. By correlating chromatin mobility proximal to the active sites of DNA synthesis, we showed that chromatin motion was locally constrained at the sites of DNA replication. Furthermore, inhibiting DNA synthesis led to increased loading of DNA polymerases. This was accompanied by accumulation of the single-stranded DNA binding protein on the chromatin and activation of DNA helicases further restricting local chromatin motion. We, therefore, propose that it is the loading of replisomes but not their catalytic activity that reduces the dynamics of replicating chromatin segments in the S-phase as well as their accessibility and probability of interactions with other genomic regions.

The location and development of Replicon Cluster Domains in early replicating DNA

Background: It has been known for many years that in metazoan cells, replication origins are organised into clusters where origins within each cluster fire near-synchronously. Despite clusters being a fundamental organising principle of metazoan DNA replication, the genomic location of origin clusters has not been documented. Methods: We synchronised human U2OS by thymidine block and release followed by L-mimosine block and release to create a population of cells progressing into S phase with a high degree of synchrony. At different times after release into S phase, cells were pulsed with EdU; the EdU-labelled DNA was then pulled down, sequenced and mapped onto the human genome. Results: The early replicating DNA showed features at a range of scales. Wavelet analysis showed that the major feature of the early replicating DNA was at a size of 500 kb, consistent with clusters of replication origins. Over the first two hours of S phase, these Replicon Cluster Domains broadened in width, consistent with their being enlarged by the progression of replication forks at their outer boundaries. The total replication signal associated with each Replicon Cluster Domain varied considerably, and this variation was reproducible and conserved over time. We provide evidence that this variability in replication signal was at least in part caused by Replicon Cluster Domains being activated at different times in different cells in the population. We also provide evidence that adjacent clusters had a statistical preference for being activated in sequence across a group, consistent with the 'domino' model of replication focus activation order observed by microscopy. Conclusions: We show that early replicating DNA is organised into Replicon Cluster Domains that behave as expected of replicon clusters observed by DNA fibre analysis. The coordinated activation of different Replicon Cluster Domains can generate the replication timing programme by which the genome is duplicated.

Neural network and kinetic modelling of human genome replication reveal replication origin locations and strengths

In human and other metazoans, the determinants of replication origin location and strength are still elusive. Origins are licensed in G1 phase and fired in S phase of the cell cycle, respectively. It is debated which of these two temporally separate steps determines origin efficiency. Experiments can independently profile mean replication timing (MRT) and replication fork directionality (RFD) genome-wide. Such profiles contain information on multiple origins' properties and on fork speed. Due to possible origin inactivation by passive replication, however, observed and intrinsic origin efficiencies can markedly differ. Thus, there is a need for methods to infer intrinsic from observed origin efficiency, which is context-dependent. Here, we show that MRT and RFD data are highly consistent with each other but contain information at different spatial scales. Using neural networks, we infer an origin licensing landscape that, when inserted in an appropriate simulation framework, jointly predicts MRT and RFD data with unprecedented precision and underlies the importance of dispersive origin firing. We furthermore uncover an analytical formula that predicts intrinsic from observed origin efficiency combined with MRT data. Comparison of inferred intrinsic origin efficiencies with experimental profiles of licensed origins (ORC, MCM) and actual initiation events (Bubble-seq, SNS-seq, OK-seq, ORM) show that intrinsic origin efficiency is not solely determined by licensing efficiency. Thus, human replication origin efficiency is set at both the origin licensing and firing steps.

Isoform-specific and ubiquitination dependent recruitment of Tet1 to replicating heterochromatin modulates methylcytosine oxidation

Oxidation of the epigenetic DNA mark 5-methylcytosine by Tet dioxygenases is an established route to diversify the epigenetic information, modulate gene expression and overall cellular (patho-)physiology. Here, we demonstrate that Tet1 and its short isoform Tet1s exhibit distinct nuclear localization during DNA replication resulting in aberrant cytosine modification levels in human and mouse cells. We show that Tet1 is tethered away from heterochromatin via its zinc finger domain, which is missing in Tet1s allowing its targeting to these regions. We find that Tet1s interacts with and is ubiquitinated by CRL4(VprBP). The ubiquitinated Tet1s is then recognized by Uhrf1 and recruited to late replicating heterochromatin. This leads to spreading of 5-methylcytosine oxidation to heterochromatin regions, LINE 1 activation and chromatin decondensation. In summary, we elucidate a dual regulation mechanism of Tet1, contributing to the understanding of how epigenetic information can be diversified by spatio-temporal directed Tet1 catalytic activity.

Convergence of SIRT1 and ATR signaling to modulate replication origin dormancy

During routine genome duplication, many potential replication origins remain inactive or 'dormant'. Such origin dormancy is achieved, in part, by an interaction with the metabolic sensor SIRT1 deacetylase. We report here that dormant origins are a group of consistent, pre-determined genomic sequences that are distinguished from baseline (i.e. ordinarily active) origins by their preferential association with two phospho-isoforms of the helicase component MCM2. During normal unperturbed cell growth, baseline origins, but not dormant origins, associate with a form of MCM2 that is phosphorylated by DBF4-dependent kinase (DDK) on serine 139 (pS139-MCM2). This association facilitates the initiation of DNA replication from baseline origins. Concomitantly, SIRT1 inhibits Ataxia Telangiectasia and Rad3-related (ATR)-kinase-mediated phosphorylation of MCM2 on serine 108 (pS108-MCM2) by deacetylating the ATR-interacting protein DNA topoisomerase II binding protein 1 (TOPBP1), thereby preventing ATR recruitment to chromatin. In cells devoid of SIRT1 activity, or challenged by replication stress, this inhibition is circumvented, enabling ATR-mediated S108-MCM2 phosphorylation. In turn, pS108-MCM2 enables DDK-mediated phosphorylation on S139-MCM2 and facilitates replication initiation at dormant origins. These observations suggest that replication origin dormancy and activation are regulated by distinct post-translational MCM modifications that reflect a balance between SIRT1 activity and ATR signaling.

Nucleases and Co-Factors in DNA Replication Stress Responses

DNA replication stress is a constant threat that cells must manage to proliferate and maintain genome integrity. DNA replication stress responses, a subset of the broader DNA damage response (DDR), operate when the DNA replication machinery (replisome) is blocked or replication forks collapse during S phase. There are many sources of replication stress, such as DNA lesions caused by endogenous and exogenous agents including commonly used cancer therapeutics, and difficult-to-replicate DNA sequences comprising fragile sites, G-quadraplex DNA, hairpins at trinucleotide repeats, and telomeres. Replication stress is also a consequence of conflicts between opposing transcription and replication, and oncogenic stress which dysregulates replication origin firing and fork progression. Cells initially respond to replication stress by protecting blocked replisomes, but if the offending problem (e.g., DNA damage) is not bypassed or resolved in a timely manner, forks may be cleaved by nucleases, inducing a DNA double-strand break (DSB) and providing a means to accurately restart stalled forks via homologous recombination. However, DSBs pose their own risks to genome stability if left unrepaired or misrepaired. Here we focus on replication stress response systems, comprising DDR signaling, fork protection, and fork processing by nucleases that promote fork repair and restart. Replication stress nucleases include MUS81, EEPD1, Metnase, CtIP, MRE11, EXO1, DNA2-BLM, SLX1-SLX4, XPF-ERCC1-SLX4, Artemis, XPG, and FEN1. Replication stress factors are important in cancer etiology as suppressors of genome instability associated with oncogenic mutations, and as potential cancer therapy targets to enhance the efficacy of chemo- and radiotherapeutics.

Preventing excess replication origin activation to ensure genome stability

Cells activate distinctive regulatory pathways that prevent excessive initiation of DNA replication to achieve timely and accurate genome duplication. Excess DNA synthesis is constrained by protein-DNA interactions that inhibit initiation at dormant origins. In parallel, specific modifications of pre-replication complexes prohibit post-replicative origin relicensing. Replication stress ensues when the controls that prevent excess replication are missing in cancer cells, which often harbor extrachromosomal DNA that can be further amplified by recombination-mediated processes to generate chromosomal translocations. The genomic instability that accompanies excess replication origin activation can provide a promising target for therapeutic intervention. Here we review molecular pathways that modulate replication origin dormancy, prevent excess origin activation, and detect, encapsulate, and eliminate persistent excess DNA.

Super-resolution microscopy reveals stochastic initiation of replication in Drosophila polytene chromosomes

Studying the probability distribution of replication initiation along a chromosome is a huge challenge. Drosophila polytene chromosomes in combination with super-resolution microscopy provide a unique opportunity for analyzing the probabilistic nature of replication initiation at the ultrastructural level. Here, we developed a method for synchronizing S-phase induction among salivary gland cells. An analysis of the replication label distribution in the first minutes of S phase and in the following hours after the induction revealed the dynamics of replication initiation. Spatial super-resolution structured illumination microscopy allowed identifying multiple discrete replication signals and to investigate the behavior of replication signals in the first minutes of the S phase at the ultrastructural level. We identified replication initiation zones where initiation occurs stochastically. These zones differ significantly in the probability of replication initiation per time unit. There are zones in which initiation occurs on most strands of the polytene chromosome in a few minutes. In other zones, the initiation on all strands takes several hours. Compact bands are free of replication initiation events, and the replication runs from outer edges to the middle, where band shapes may alter.

Polo-like kinase 1 (Plk1) regulates DNA replication origin firing and interacts with Rif1 in Xenopus

The activation of eukaryotic DNA replication origins needs to be strictly controlled at multiple steps in order to faithfully duplicate the genome and to maintain its stability. How the checkpoint recovery and adaptation protein Polo-like kinase 1 (Plk1) regulates the firing of replication origins during non-challenged S phase remained an open question. Using DNA fiber analysis, we show that immunodepletion of Plk1 in the Xenopus in vitro system decreases replication fork density and initiation frequency. Numerical analyses suggest that Plk1 reduces the overall probability and synchrony of origin firing. We used quantitative chromatin proteomics and co-immunoprecipitations to demonstrate that Plk1 interacts with firing factors MTBP/Treslin/TopBP1 as well as with Rif1, a known regulator of replication timing. Phosphopeptide analysis by LC/MS/MS shows that the C-terminal domain of Rif1, which is necessary for its repressive action on origins through protein phosphatase 1 (PP1), can be phosphorylated in vitro by Plk1 on S2058 in its PP1 binding site. The phosphomimetic S2058D mutant interrupts the Rif1-PP1 interaction and modulates DNA replication. Collectively, our study provides molecular insights into how Plk1 regulates the spatio-temporal replication program and suggests that Plk1 controls origin activation at the level of large chromatin domains in vertebrates.

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

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

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

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

Transcription-coupled structural dynamics of topologically associating domains regulate replication origin efficiency

BACKGROUND

Metazoan cells only utilize a small subset of the potential DNA replication origins to duplicate the whole genome in each cell cycle. Origin choice is linked to cell growth, differentiation, and replication stress. Although various genetic and epigenetic signatures have been linked to the replication efficiency of origins, there is no consensus on how the selection of origins is determined.

RESULTS

We apply dual-color stochastic optical reconstruction microscopy (STORM) super-resolution imaging to map the spatial distribution of origins within individual topologically associating domains (TADs). We find that multiple replication origins initiate separately at the spatial boundary of a TAD at the beginning of the S phase. Intriguingly, while both high-efficiency and low-efficiency origins are distributed homogeneously in the TAD during the G1 phase, high-efficiency origins relocate to the TAD periphery before the S phase. Origin relocalization is dependent on both transcription and CTCF-mediated chromatin structure. Further, we observe that the replication machinery protein PCNA forms immobile clusters around TADs at the G1/S transition, explaining why origins at the TAD periphery are preferentially fired.

CONCLUSION

Our work reveals a new origin selection mechanism that the replication efficiency of origins is determined by their physical distribution in the chromatin domain, which undergoes a transcription-dependent structural re-organization process. Our model explains the complex links between replication origin efficiency and many genetic and epigenetic signatures that mark active transcription. The coordination between DNA replication, transcription, and chromatin organization inside individual TADs also provides new insights into the biological functions of sub-domain chromatin structural dynamics.

Mammalian DNA Replication Timing

Immediately following the discovery of the structure of DNA and the semi-conservative replication of the parental DNA sequence into two new DNA strands, it became apparent that DNA replication is organized in a temporal and spatial fashion during the S phase of the cell cycle, correlated with the large-scale organization of chromatin in the nucleus. After many decades of limited progress, technological advances in genomics, genome engineering, and imaging have finally positioned the field to tackle mechanisms underpinning the temporal and spatial regulation of DNA replication and the causal relationships between DNA replication and other features of large-scale chromosome structure and function. In this review, we discuss these major recent discoveries as well as expectations for the coming decade.

Super-resolved live-cell imaging using random illumination microscopy

Current super-resolution microscopy (SRM) methods suffer from an intrinsic complexity that might curtail their routine use in cell biology. We describe here random illumination microscopy (RIM) for live-cell imaging at super-resolutions matching that of 3D structured illumination microscopy, in a robust fashion. Based on speckled illumination and statistical image reconstruction, easy to implement and user-friendly, RIM is unaffected by optical aberrations on the excitation side, linear to brightness, and compatible with multicolor live-cell imaging over extended periods of time. We illustrate the potential of RIM on diverse biological applications, from the mobility of proliferating cell nuclear antigen (PCNA) in U2OS cells and kinetochore dynamics in mitotic S. pombe cells to the 3D motion of myosin minifilaments deep inside Drosophila tissues. RIM's inherent simplicity and extended biological applicability, particularly for imaging at increased depths, could help make SRM accessible to biology laboratories.

Developmental differences in genome replication program and origin activation

To ensure error-free duplication of all (epi)genetic information once per cell cycle, DNA replication follows a cell type and developmental stage specific spatio-temporal program. Here, we analyze the spatio-temporal DNA replication progression in (un)differentiated mouse embryonic stem (mES) cells. Whereas telomeres replicate throughout S-phase, we observe mid S-phase replication of (peri)centromeric heterochromatin in mES cells, which switches to late S-phase replication upon differentiation. This replication timing reversal correlates with and depends on an increase in condensation and a decrease in acetylation of chromatin. We further find synchronous duplication of the Y chromosome, marking the end of S-phase, irrespectively of the pluripotency state. Using a combination of single-molecule and super-resolution microscopy, we measure molecular properties of the mES cell replicon, the number of replication foci active in parallel and their spatial clustering. We conclude that each replication nanofocus in mES cells corresponds to an individual replicon, with up to one quarter representing unidirectional forks. Furthermore, with molecular combing and genome-wide origin mapping analyses, we find that mES cells activate twice as many origins spaced at half the distance than somatic cells. Altogether, our results highlight fundamental developmental differences on progression of genome replication and origin activation in pluripotent cells.

Are the processes of DNA replication and DNA repair reading a common structural chromatin unit?

Decades of investigation on genomic DNA have brought us deeper insights into its organization within the nucleus and its metabolic mechanisms. This was fueled by the parallel development of experimental techniques and has stimulated model building to simulate genome conformation in agreement with the experimental data. Here, we will discuss our recent discoveries on the chromatin units of DNA replication and DNA damage response. We will highlight their remarkable structural similarities and how both revealed themselves as clusters of nanofocal structures each on the hundred thousand base pair size range corresponding well with chromatin loop sizes. We propose that the function of these two global genomic processes is determined by the loop level organization of chromatin structure with structure dictating function.

Abbreviations: 3D-SIM: 3D-structured illumination microscopy; 3C: chromosome conformation capture; DDR: DNA damage response; FISH: fluorescent in situ hybridization; Hi-C: high conformation capture; HiP-HoP: highly predictive heteromorphic polymer model; IOD: inter-origin distance; LAD: lamina associated domain; STED: stimulated emission depletion microscopy; STORM: stochastic optical reconstruction microscopy; SBS: strings and binders switch model; TAD: topologically associated domain

Chromatin and Nuclear Architecture: Shaping DNA Replication in 3D

In eukaryotes, DNA replication progresses through a finely orchestrated temporal and spatial program. The 3D genome structure and nuclear architecture have recently emerged as fundamental determinants of the replication program. Factors with established roles in replication have been recognized as genome organization regulators. Exploiting paradigms from yeasts and mammals, we discuss how DNA replication is regulated in time and space through DNA-associated trans-acting factors, diffusible limiting replication initiation factors, higher-order chromatin folding, dynamic origin localization, and specific nuclear microenvironments. We present an integrated model for the regulation of DNA replication in 3D and highlight the importance of accurate spatio-temporal regulation of DNA replication in physiology and disease.

Processive DNA synthesis is associated with localized decompaction of constitutive heterochromatin at the sites of DNA replication and repair

Constitutive heterochromatin is considered as a functionally inert genome compartment, important for its architecture and stability. How such stable structure is maintained is not well understood. Here, we apply four different visualization schemes to label it and investigate its dynamics during DNA replication and repair. We show that replisomes assemble over the heterochromatin in a temporally ordered manner. Furthermore, heterochromatin undergoes transient decompaction locally at the active sites of DNA synthesis. Using selective laser microirradiation conditions that lead to damage repaired via processive DNA synthesis, we measured similarly local decompaction of heterochromatin. In both cases, we could not observe large-scale movement of heterochromatin to the domain surface. Instead, the processive DNA synthesis machinery assembled at the replication/repair sites. Altogether, our data are compatible with a progression of DNA replication/repair along the chromatin in a dynamic mode with localized and transient decompaction that does not globally remodels the whole heterochromatin compartment.

Deciphering DNA methylation signatures of pancreatic cancer and pancreatitis

BACKGROUND

Chronic pancreatitis presents a high risk of inflammation-related progression to pancreatic cancer. Pancreatic cancer is the fourth leading cause of cancer-related death worldwide. The high mortality rate is directly related to the difficulty in promptly diagnosing the disease, which often presents as overt and advanced. Hence, early diagnosis for pancreatic cancer becomes crucial, propelling research into the molecular and epigenetic landscape of the disease.

MAIN BODY

Recent studies have shown that cell-free DNA methylation profiles from inflammatory diseases or cancer can vary, thus opening a new venue for the development of biomarkers for early diagnosis. In particular, cell-free DNA methylation could be employed in the identification of pre-neoplastic signatures in individuals with suspected pancreatic conditions, representing a specific and non-invasive method of early diagnosis of pancreatic cancer. In this review, we describe the molecular determinants of pancreatic cancer and how these are related to chronic pancreatitis. We will then present an overview of differential methylated genes in the two conditions, highlighting their diagnostic or prognostic potential.

CONCLUSION

Exploiting the relation between abnormally methylated cell-free DNA and pre-neoplastic lesions or chronic pancreatitis may become a game-changing approach for the development of tools for the early diagnosis of pancreatic cancer.

Peripheral re-localization of constitutive heterochromatin advances its replication timing and impairs maintenance of silencing marks

The replication of the genome is a highly organized process, both spatially and temporally. Although a lot is known on the composition of the basic replication machinery, how its activity is regulated is mostly unknown. Several chromatin properties have been proposed as regulators, but a potential role of the nuclear DNA position remains unclear. We made use of the prominent structure and well-defined heterochromatic landscape of mouse pericentric chromosome domains as a well-studied example of late replicating constitutive heterochromatin. We established a method to manipulate its nuclear position and evaluated the effect on replication timing, DNA compaction and epigenetic composition. Using time-lapse microscopy, we observed that constitutive heterochromatin, known to replicate during late S-phase, was replicated in mid S-phase when repositioned to the nuclear periphery. Out-of-schedule replication resulted in deficient post-replicative maintenance of chromatin modifications, namely silencing marks. We propose that repositioned constitutive heterochromatin was activated in trans according to the domino model of origin firing by nearby (mid S) firing origins. In summary, our data provide, on the one hand, a novel approach to manipulate nuclear DNA position and, on the other hand, establish nuclear DNA position as a novel mechanism regulating DNA replication timing and epigenetic maintenance.

Dormant origin signaling during unperturbed replication

Mechanisms that limit origin firing are essential as the ˜50,000 origins that replicate the human genome in unperturbed cells are chosen from an excess of ˜500,000 licensed origins. Computational models of the spatiotemporal pattern of replication foci assume that origins fire stochastically with a domino-like progression that places later firing origins near recent fired origins. These stochastic models of origin firing require dormant origin signaling that inhibits origin firing and suppresses licensed origins for passive replication at a distance of ∼7-120 kbp around replication forks. ATR and CHK1 kinase inhibitors increase origin firing and increase origin density in unperturbed cells. Thus, basal ATR and CHK1 kinase-dependent dormant origin signaling inhibits origin firing and there appear to be two thresholds of ATR kinase signaling. A minority of ATR molecules are activated for ATR and CHK1 kinase-dependent dormant origin signaling and this is essential for DNA replication in unperturbed cells. A majority of ATR molecules are activated for ATR and CHK1 kinase-dependent checkpoint signaling in cells treated with DNA damaging agents that target replication forks. Since ATR and CHK1 kinase inhibitors increase origin firing and this is associated with fork stalling and extensive regions of single-stranded DNA, they are DNA damaging agents. Accordingly, the sequence of administration of ATR and CHK1 kinase inhibitors and DNA damaging agents may impact the DNA damage induced by the combination and the efficacy of cell killing by the combination.

An ATR and CHK1 kinase signaling mechanism that limits origin firing during unperturbed DNA replication

The 50,000 origins that replicate the human genome are selected from an excess of licensed origins. Firing licensed origins that would otherwise be passively replicated is a simple mechanism to recover DNA replication between stalled replication forks. This plasticity in origin use promotes genome stability if an unknown mechanism prevents a subset of origins from firing during unperturbed DNA replication. We describe ATR and CHK1 kinase signaling that suppresses a CDK1 kinase-dependent phosphorylation on the chromatin protein RIF1. The CDK1 kinase-dependent phosphorylation of RIF1 disrupts its interaction with PP1 phosphatase. Thus, ATR and CHK1 stabilize an interaction between RIF1 and PP1 that counteracts CDC7 and CDK2 kinase signaling at licensed origins. This mechanism limits origin firing during unperturbed DNA replication.

DNA damage-induced signaling by ATR and CHK1 inhibits DNA replication, stabilizes stalled and collapsed replication forks, and mediates the repair of multiple classes of DNA lesions. We and others have shown that ATR kinase inhibitors, three of which are currently undergoing clinical trials, induce excessive origin firing during unperturbed DNA replication, indicating that ATR kinase activity limits replication initiation in the absence of damage. However, the origins impacted and the underlying mechanism(s) have not been described. Here, we show that unperturbed DNA replication is associated with a low level of ATR and CHK1 kinase signaling and that inhibition of this signaling induces dormant origin firing at sites of ongoing replication throughout the S phase. We show that ATR and CHK1 kinase inhibitors induce RIF1 Ser2205 phosphorylation in a CDK1-dependent manner, which disrupts an interaction between RIF1 and PP1 phosphatase. Thus, ATR and CHK1 signaling suppresses CDK1 kinase activity throughout the S phase and stabilizes an interaction between RIF1 and PP1 in replicating cells. PP1 dephosphorylates key CDC7 and CDK2 kinase substrates to inhibit the assembly and activation of the replicative helicase. This mechanism limits origin firing during unperturbed DNA replication in human cells.

DNA Replication Timing Enters the Single-Cell Era

In mammalian cells, DNA replication timing is controlled at the level of megabase (Mb)-sized chromosomal domains and correlates well with transcription, chromatin structure, and three-dimensional (3D) genome organization. Because of these properties, DNA replication timing is an excellent entry point to explore genome regulation at various levels and a variety of studies have been carried out over the years. However, DNA replication timing studies traditionally required at least tens of thousands of cells, and it was unclear whether the replication domains detected by cell population analyses were preserved at the single-cell level. Recently, single-cell DNA replication profiling methods became available, which revealed that the Mb-sized replication domains detected by cell population analyses were actually well preserved in individual cells. In this article, we provide a brief overview of our current knowledge on DNA replication timing regulation in mammals based on cell population studies, outline the findings from single-cell DNA replication profiling, and discuss future directions and challenges.

DNA replication dynamics of vole genome and its epigenetic regulation

Background

The genome of some vole rodents exhibit large blocks of heterochromatin coupled to their sex chromosomes. The DNA composition and transcriptional activity of these heterochromatin blocks have been studied, but little is known about their DNA replication dynamics and epigenetic composition.

Results

Here, we show prominent epigenetic marks of the heterochromatic blocks in the giant sex chromosomes of female Microtus cabrerae cells. While the X chromosomes are hypoacetylated and cytosine hypomethylated, they are either enriched for macroH2A and H3K27me3 typical for facultative heterochromatin or for H3K9me3 and HP1 beta typical for constitutive heterochromatin. Using pulse-chase replication labeling and time-lapse microscopy, we found that the heterochromatic block enriched for macroH2A/H3K27me3 of the X chromosome is replicated during mid-S-phase, prior to the heterochromatic block enriched for H3K9me3/HP1 beta, which is replicated during late S-phase. To test whether histone acetylation level regulates its replication dynamics, we induced either global hyperacetylation by pharmacological inhibition or by targeting a histone acetyltransferase to the heterochromatic region of the X chromosomes. Our data reveal that histone acetylation level affects DNA replication dynamics of the sex chromosomes’ heterochromatin and leads to a global reduction in replication fork rate genome wide.

Conclusions

In conclusion, we mapped major epigenetic modifications controlling the structure of the sex chromosome-associated heterochromatin and demonstrated the occurrence of differences in the molecular mechanisms controlling the replication timing of the heterochromatic blocks at the sex chromosomes in female Microtus cabrerae cells. Furthermore, we highlighted a conserved role of histone acetylation level on replication dynamics across mammalian species.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0262-0) contains supplementary material, which is available to authorized users.

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.

Snapshots of archaeal DNA replication and repair in living cells using super-resolution imaging

Using the haloarchaeon Haloferax volcanii as a model, we developed nascent DNA labeling and the functional GFP-labeled single-stranded binding protein RPA2 as novel tools to gain new insight into DNA replication and repair in live haloarchaeal cells. Our quantitative fluorescence microscopy data revealed that RPA2 forms distinct replication structures that dynamically responded to replication stress and DNA damaging agents. The number of the RPA2 foci per cell followed a probabilistic Poisson distribution, implying hitherto unnoticed stochastic cell-to-cell variation in haloarchaeal DNA replication and repair processes. The size range of haloarchaeal replication structures is very similar to those observed earlier in eukaryotic cells. The improved lateral resolution of 3D-SIM fluorescence microscopy allowed proposing that inhibition of DNA synthesis results in localized replication foci clustering and facilitated observation of RPA2 complexes brought about by chemical agents creating DNA double-strand breaks. Altogether our in vivo observations are compatible with earlier in vitro studies on archaeal single-stranded DNA binding proteins. Our work thus underlines the great potential of live cell imaging for unraveling the dynamic nature of transient molecular interactions that underpin fundamental molecular processes in the Third domain of life.

The Protective Role of Dormant Origins in Response to Replicative Stress

Genome stability requires tight regulation of DNA replication to ensure that the entire genome of the cell is duplicated once and only once per cell cycle. In mammalian cells, origin activation is controlled in space and time by a cell-specific and robust program called replication timing. About 100,000 potential replication origins form on the chromatin in the gap 1 (G1) phase but only 20⁻30% of them are active during the DNA replication of a given cell in the synthesis (S) phase. When the progress of replication forks is slowed by exogenous or endogenous impediments, the cell must activate some of the inactive or "dormant" origins to complete replication on time. Thus, the many origins that may be activated are probably key to protect the genome against replication stress. This review aims to discuss the role of these dormant origins as safeguards of the human genome during replicative stress.

Regulation of the initiation of DNA replication in human cells

The origin of species would not have been possible without high fidelity DNA replication and complex genomes evolved with mechanisms that control the initiation of DNA replication at multiple origins on multiple chromosomes such that the genome is duplicated once and only once. The mechanisms that control the assembly and activation of the replicative helicase and the initiation of DNA replication in yeast and Xenopus egg extract systems have been identified and reviewed [1,2]. The goal of this review is to organize currently available data on the mechanisms that control the initiation of DNA replication in human cells.

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

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

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

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

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

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

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

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

Replication timing and nuclear structure

DNA replication proceeds along spatially and temporally coordinated patterns within the nucleus, thus protecting the genome during the synthesis of new genetic material. While we have been able to visualize replication patterns on DNA fibers for 50 years, recent developments and discoveries have provided a greater insight into how DNA replication is controlled. In this review, we highlight many of these discoveries. Of great interest are the physiological role of the replication timing program, cis and trans-acting factors that modulate replication timing and the effects of chromatin structure on the replication timing program. We also discuss future directions in the study of replication timing.

A journey through the microscopic ages of DNA replication

Scientific discoveries and technological advancements are inseparable but not always take place in a coherent chronological manner. In the next, we will provide a seemingly unconnected and serendipitous series of scientific facts that, in the whole, converged to unveil DNA and its duplication. We will not cover here the many and fundamental contributions from microbial genetics and in vitro biochemistry. Rather, in this journey, we will emphasize the interplay between microscopy development culminating on super resolution fluorescence microscopy (i.e., nanoscopy) and digital image analysis and its impact on our understanding of DNA duplication. We will interlace the journey with landmark concepts and experiments that have brought the cellular DNA replication field to its present state.

Anatomy of Mammalian Replication Domains

Genetic information is faithfully copied by DNA replication through many rounds of cell division. In mammals, DNA is replicated in Mb-sized chromosomal units called “replication domains.” While genome-wide maps in multiple cell types and disease states have uncovered both dynamic and static properties of replication domains, we are still in the process of understanding the mechanisms that give rise to these properties. A better understanding of the molecular basis of replication domain regulation will bring new insights into chromosome structure and function.

Order from clutter: selective interactions at mammalian replication origins

Mammalian chromosome duplication progresses in a precise order and is subject to constraints that are often relaxed in developmental disorders and malignancies. Molecular information about the regulation of DNA replication at the chromatin level is lacking because protein complexes that initiate replication seem to bind chromatin indiscriminately. High-throughput sequencing and mathematical modelling have yielded detailed genome-wide replication initiation maps. Combining these maps and models with functional genetic analyses suggests that distinct DNA-protein interactions at subgroups of replication initiation sites (replication origins) modulate the ubiquitous replication machinery and supports an emerging model that delineates how indiscriminate DNA-binding patterns translate into a consistent, organized replication programme.

Regulation of Replication Origins

In eukaryotes, genome duplication starts concomitantly at many replication initiation sites termed replication origins. The replication initiation program is spatially and temporally coordinated to ensure accurate, efficient DNA synthesis that duplicates the entire genome while maintaining other chromatin-dependent functions. Unlike in prokaryotes, not all potential replication origins in eukaryotes are needed for complete genome duplication during each cell cycle. Instead, eukaryotic cells vary the use of initiation sites so that only a fraction of potential replication origins initiate replication each cell cycle. Flexibility in origin choice allows each eukaryotic cell type to utilize different initiation sites, corresponding to unique nuclear DNA packaging patterns. These patterns coordinate replication with gene expression and chromatin condensation. Budding yeast replication origins share a consensus sequence that marks potential initiation sites. Metazoan origins, on the other hand, lack a consensus sequence. Rather, they are associated with a collection of structural features, chromatin packaging features, histone modifications, transcription, and DNA-DNA/DNA-protein interactions. These features confer cell type-specific replication and expression and play an essential role in maintaining genomic stability.

4D Visualization of replication foci in mammalian cells corresponding to individual replicons

Since the pioneering proposal of the replicon model of DNA replication 50 years ago, the predicted replicons have not been identified and quantified at the cellular level. Here, we combine conventional and super-resolution microscopy of replication sites in live and fixed cells with computational image analysis. We complement these data with genome size measurements, comprehensive analysis of S-phase dynamics and quantification of replication fork speed and replicon size in human and mouse cells. These multidimensional analyses demonstrate that replication foci (RFi) in three-dimensional (3D) preserved somatic mammalian cells can be optically resolved down to single replicons throughout S-phase. This challenges the conventional interpretation of nuclear RFi as replication factories, that is, the complex entities that process multiple clustered replicons. Accordingly, 3D genome organization and duplication can be now followed within the chromatin context at the level of individual replicons.