Genome-wide mapping of individual replication fork velocities using nanopore sequencing

DNA nicks in both leading and lagging strand templates can trigger break-induced replication

Encounters between replication forks and unrepaired DNA single-strand breaks (SSBs) can generate both single-ended and double-ended double-strand breaks (seDSBs and deDSBs). seDSBs can be repaired by break-induced replication (BIR), which is a highly mutagenic pathway that is thought to be responsible for many of the mutations and genome rearrangements that drive cancer development. However, the frequency of BIR's deployment and its ability to be triggered by both leading and lagging template strand SSBs were unclear. Using site- and strand-specific SSBs generated by nicking enzymes, including CRISPR-Cas9 nickase (Cas9n), we demonstrate that leading and lagging template strand SSBs in fission yeast are typically converted into deDSBs that are repaired by homologous recombination. However, both types of SSBs can also trigger BIR, and the frequency of these events increases when fork convergence is delayed and the non-homologous end joining protein Ku70 is deleted.

Telomere-to-telomere DNA replication timing profiling using single-molecule sequencing with Nanotiming

Current temporal studies of DNA replication are either low-resolution or require complex cell synchronisation and/or sorting procedures. Here we introduce Nanotiming, a single-molecule, nanopore sequencing-based method producing high-resolution, telomere-to-telomere replication timing (RT) profiles of eukaryotic genomes by interrogating changes in intracellular dTTP concentration during S phase through competition with its analogue bromodeoxyuridine triphosphate (BrdUTP) for incorporation into replicating DNA. This solely demands the labelling of asynchronously growing cells with an innocuous dose of BrdU during one doubling time followed by BrdU quantification along nanopore reads. We demonstrate in S. cerevisiae model eukaryote that Nanotiming reproduces RT profiles generated by reference methods both in wild-type and mutant cells inactivated for known RT determinants. Nanotiming is simple, accurate, inexpensive, amenable to large-scale analyses, and has the unique ability to access RT of individual telomeres, revealing that Rif1 iconic telomere regulator selectively delays replication of telomeres associated with specific subtelomeric elements.

The fork protection complex generates DNA topological stress-induced DNA damage while ensuring full and faithful genome duplication

The fork protection complex (FPC), composed of Mrc1, Tof1, and Csm3, supports rapid and stable DNA replication. Here, we show that FPC activity also introduces DNA damage by increasing DNA topological stress during replication. Mrc1 action increases DNA topological stress during plasmid replication, while Mrc1 or Tof1 activity causes replication stress and DNA damage within topologically constrained regions. We show that the recruitment of Top1 to the fork by Tof1 suppresses the DNA damage generated in these loci. While FPC activity introduces some DNA damage due to increased topological stress, the FPC is also necessary to prevent DNA damage in long replicons across the genome, indicating that the FPC is required for complete and faithful genome duplication. We conclude that FPC regulation must balance ensuring full genome duplication through rapid replication with minimizing the consequential DNA topological stress-induced DNA damage caused by rapid replication through constrained regions.

Mathematical model for the distribution of DNA replication origins

DNA replication in yeast and in many other organisms starts from well-defined locations on the DNA known as replication origins. The spatial distribution of these origins in the genome is particularly important in ensuring that replication is completed quickly. Cells are more vulnerable to DNA damage and other forms of stress while they are replicating their genome. This raises the possibility that the spatial distribution of origins is under selection pressure. In this paper we investigate the hypothesis that natural selection favors origin distributions leading to shorter replication times. Using a simple mathematical model, we show that this hypothesis leads to two main predictions about the origin distributions: that neighboring origins that are inefficient (less likely to fire) are more likely to be close to each other than efficient origins; and that neighboring origins with larger differences in firing times are more likely to be close to each other than origins with similar firing times. We test these predictions using next-generation sequencing data, and show that they are both supported by the data.

Quantifying DNA replication speeds in single cells by scEdU-seq

In a human cell, thousands of replication forks simultaneously coordinate duplication of the entire genome. The rate at which this process occurs might depend on the epigenetic state of the genome and vary between, or even within, cell types. To accurately measure DNA replication speeds, we developed single-cell 5-ethynyl-2'-deoxyuridine sequencing to detect nascent replicated DNA. We observed that the DNA replication speed is not constant but increases during S phase of the cell cycle. Using genetic and pharmacological perturbations we were able to alter this acceleration of replication and conclude that DNA damage inflicted by the process of transcription limits the speed of replication during early S phase. In late S phase, during which less-transcribed regions replicate, replication accelerates and approaches its maximum speed.

Monitoring and quantifying replication fork dynamics with high-throughput methods

Before each cell division, eukaryotic cells must replicate their chromosomes to ensure the accurate transmission of genetic information. Chromosome replication involves more than just DNA duplication; it also includes chromatin assembly, inheritance of epigenetic marks, and faithful resumption of all genomic functions after replication. Recent progress in quantitative technologies has revolutionized our understanding of the complexity and dynamics of DNA replication forks at both molecular and genomic scales. Here, we highlight the pivotal role of these novel methods in uncovering the principles and mechanisms of chromosome replication. These technologies have illuminated the regulation of genome replication programs, quantified the impact of DNA replication on genomic mutations and evolutionary processes, and elucidated the mechanisms of replication-coupled chromatin assembly and epigenome maintenance.

Identification of 1600 replication origins in S. cerevisiae

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, 'a site that binds Mcm in G1' might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here, we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least three orders of magnitude. Published data quantifying single-stranded DNA (ssDNA) during S phase revealed replication initiation among the most abundant 1600 of these sites, with replication activity decreasing with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5500 sites. Specifically, these sites: (1) appeared in intergenic nucleosome-free regions flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. We conclude that, if sites at which Mcm double hexamers are loaded can function as replication origins, then DNA replication origins are at least threefold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than widely assumed.

The genetic landscape of origins of replication in P. falciparum

Various origin mapping approaches have enabled genome-wide identification of origins of replication (ORI) in model organisms, but only a few studies have focused on divergent organisms. By employing three complementary approaches we provide a high-resolution map of ORIs in Plasmodium falciparum, the deadliest human malaria parasite. We profiled the distribution of origin of recognition complex (ORC) binding sites by ChIP-seq of two PfORC subunits and mapped active ORIs using NFS and SNS-seq. We show that ORIs lack sequence specificity but are not randomly distributed, and group in clusters. Licensing is biased towards regions of higher GC content and associated with G-quadruplex forming sequences (G4FS). While strong transcription likely enhances firing, active origins are depleted from transcription start sites. Instead, most accumulate in transcriptionally active gene bodies. Single molecule analysis of nanopore reads containing multiple initiation events, which could have only come from individual nuclei, showed a relationship between the replication fork pace and the distance to the nearest origin. While some similarities were drawn with the canonic eukaryote model, the distribution of ORIs in P. falciparum is likely shaped by unique genomic features such as extreme AT-richness-a product of evolutionary pressure imposed by the parasitic lifestyle.

MBD1 protects replication fork stability by recruiting PARP1 and controlling transcription-replication conflicts

The replication-stress response is essential to ensure the faithful transmission of genetic information to daughter cells. Although several stress-resolution pathways have been identified to deal with replication stress, the precise regulatory mechanisms for replication fork stability are not fully understood. Our study identified Methyl-CpG Binding Domain 1 (MBD1) as essential for the maintaining genomic stability and protecting stalled replication forks in mammalian cells. Depletion of MBD1 increases DNA lesions and sensitivity to replication stress. Mechanistically, we found that loss of MBD1 leads to the dissociation of Poly(ADP-ribose) polymerase 1 (PARP1) from the replication fork, potentially accelerating fork progression and resulting in higher levels of transcription-replication conflicts (T-R conflicts). Using a proximity ligation assay combined with 5-ethynyl-2'-deoxyuridine, we revealed that the MBD1 and PARP1 proteins were recruited to stalled forks under hydroxyurea (HU) treatment. In addition, our study showed that the level of R-loops also increased in MBD1-delated cells. Without MBD1, stalled replication forks resulting from T-R conflicts were primarily degraded by the DNA2 nuclease. Our findings shed light on a new aspect of MBD1 in maintaining genome stability and providing insights into the mechanisms underlying replication stress response.

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.

Identification of 1600 replication origins in S. cerevisiae

There are approximately 500 known origins of replication in the yeast genome, and the process by which DNA replication initiates at these locations is well understood. In particular, these sites are made competent to initiate replication by loading of the Mcm replicative helicase prior to the start of S phase; thus, "a site to which MCM is bound in G1" might be considered to provide an operational definition of a replication origin. By fusing a subunit of Mcm to micrococcal nuclease, a technique referred to as "Chromatin Endogenous Cleavage", we previously showed that known origins are typically bound by a single Mcm double hexamer, loaded adjacent to the ARS consensus sequence (ACS). Here we extend this analysis from known origins to the entire genome, identifying candidate Mcm binding sites whose signal intensity varies over at least 3 orders of magnitude. Published data quantifying the production of ssDNA during S phase showed clear evidence of replication initiation among the most abundant 1600 of these sites, with replication activity decreasing in concert with Mcm abundance and disappearing at the limit of detection of ssDNA. Three other hallmarks of replication origins were apparent among the most abundant 5,500 sites. Specifically, these sites (1) appeared in intergenic nucleosome-free regions that were flanked on one or both sides by well-positioned nucleosomes; (2) were flanked by ACSs; and (3) exhibited a pattern of GC skew characteristic of replication initiation. Furthermore, the high resolution of this technique allowed us to demonstrate a strong bias for detecting Mcm double-hexamers downstream rather than upstream of the ACS, which is consistent with the directionality of Mcm loading by Orc that has been observed in vitro. We conclude that, if sites at which Mcm double-hexamers are loaded can function as replication origins, then DNA replication origins are at least 3-fold more abundant than previously assumed, and we suggest that replication may occasionally initiate in essentially every intergenic region. These results shed light on recent reports that as many as 15% of replication events initiate outside of known origins, and this broader distribution of replication origins suggest that S phase in yeast may be less distinct from that in humans than is widely assumed.

Replisome-cohesin interactions provided by the Tof1-Csm3 and Mrc1 cohesion establishment factors

The chromosomal cohesin complex establishes sister chromatid cohesion during S phase, which forms the basis for faithful segregation of DNA replication products during cell divisions. Cohesion establishment is defective in the absence of either of three non-essential Saccharomyces cerevisiae replication fork components Tof1-Csm3 and Mrc1. Here, we investigate how these conserved factors contribute to cohesion establishment. Tof1-Csm3 and Mrc1 serve known roles during DNA replication, including replication checkpoint signaling, securing replication fork speed, as well as recruiting topoisomerase I and the histone chaperone FACT. By modulating each of these functions independently, we rule out that one of these known replication roles explains the contribution of Tof1-Csm3 and Mrc1 to cohesion establishment. Instead, using purified components, we reveal direct and multipronged protein interactions of Tof1-Csm3 and Mrc1 with the cohesin complex. Our findings open the possibility that a series of physical interactions between replication fork components and cohesin facilitate successful establishment of sister chromatid cohesion during DNA replication.

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.