Genome-wide studies highlight indirect links between human replication origins and gene regulation

Integrative analysis of DNA replication origins and ORC-/MCM-binding sites in human cells reveals a lack of overlap

Based on experimentally determined average inter-origin distances of ~100 kb, DNA replication initiates from ~50,000 origins on human chromosomes in each cell cycle. The origins are believed to be specified by binding of factors like the origin recognition complex (ORC) or CTCF or other features like G-quadruplexes. We have performed an integrative analysis of 113 genome-wide human origin profiles (from five different techniques) and five ORC-binding profiles to critically evaluate whether the most reproducible origins are specified by these features. Out of ~7.5 million union origins identified by all datasets, only 0.27% (20,250 shared origins) were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques, suggesting extensive variability in origin usage and identification. Also, 21% of the shared origins overlap with transcriptional promoters, posing a conundrum. Although the shared origins overlap more than union origins with constitutive CTCF-binding sites, G-quadruplex sites, and activating histone marks, these overlaps are comparable or less than that of known transcription start sites, so that these features could be enriched in origins because of the overlap of origins with epigenetically open, promoter-like sequences. Only 6.4% of the 20,250 shared origins were within 1 kb from any of the ~13,000 reproducible ORC-binding sites in human cancer cells, and only 4.5% were within 1 kb of the ~11,000 union MCM2-7-binding sites in contrast to the nearly 100% overlap in the two comparisons in the yeast, Saccharomyces cerevisiae. Thus, in human cancer cell lines, replication origins appear to be specified by highly variable stochastic events dependent on the high epigenetic accessibility around promoters, without extensive overlap between the most reproducible origins and currently known ORC- or MCM-binding sites.

A Cluster of Evolutionarily Recent KRAB Zinc Finger Proteins Protects Cancer Cells from Replicative Stress-Induced Inflammation

Heterochromatin loss and genetic instability enhance cancer progression by favoring clonal diversity, yet uncontrolled replicative stress leads to mitotic catastrophe and inflammatory responses that promote immune rejection. KRAB domain-containing zinc finger proteins (KZFP) contribute to heterochromatin maintenance at transposable elements (TE). Here, we identified an association of upregulation of a cluster of primate-specific KZFPs with poor prognosis, increased copy-number alterations, and changes in the tumor microenvironment in diffuse large B-cell lymphoma (DLBCL). Depleting two of these KZFPs targeting evolutionarily recent TEs, ZNF587 and ZNF417, impaired the proliferation of cells derived from DLBCL and several other tumor types. ZNF587 and ZNF417 depletion led to heterochromatin redistribution, replicative stress, and cGAS-STING-mediated induction of an interferon/inflammatory response, which enhanced susceptibility to macrophage-mediated phagocytosis and increased surface expression of HLA-I, together with presentation of a neoimmunopeptidome. Thus, cancer cells can exploit KZFPs to dampen TE-originating surveillance mechanisms, which likely facilitates clonal expansion, diversification, and immune evasion.

SIGNIFICANCE

Upregulation of a cluster of primate-specific KRAB zinc finger proteins in cancer cells prevents replicative stress and inflammation by regulating heterochromatin maintenance, which could facilitate the development of improved biomarkers and treatments.

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.

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.

Dimeric G-quadruplex motifs-induced NFRs determine strong replication origins in vertebrates

Replication of vertebrate genomes is tightly regulated to ensure accurate duplication, but our understanding of the interplay between genetic and epigenetic factors in this regulation remains incomplete. Here, we investigated the involvement of three elements enriched at gene promoters and replication origins: guanine-rich motifs potentially forming G-quadruplexes (pG4s), nucleosome-free regions (NFRs), and the histone variant H2A.Z, in the firing of origins of replication in vertebrates. We show that two pG4s on the same DNA strand (dimeric pG4s) are sufficient to induce the assembly of an efficient minimal replication origin without inducing transcription in avian DT40 cells. Dimeric pG4s in replication origins are associated with formation of an NFR next to precisely-positioned nucleosomes enriched in H2A.Z on this minimal origin and genome-wide. Thus, our data suggest that dimeric pG4s are important for the organization and duplication of vertebrate genomes. It supports the hypothesis that a nucleosome close to an NFR is a shared signal for the formation of replication origins in eukaryotes.

Loss of Ezh2 function remodels the DNA replication initiation landscape

In metazoan cells, DNA replication initiates from thousands of genomic loci scattered throughout the genome called DNA replication origins. Origins are strongly associated with euchromatin, particularly open genomic regions such as promoters and enhancers. However, over a third of transcriptionally silent genes are associated with DNA replication initiation. Most of these genes are bound and repressed by the Polycomb repressive complex-2 (PRC2) through the repressive H3K27me3 mark. This is the strongest overlap observed for a chromatin regulator with replication origin activity. Here, we asked whether Polycomb-mediated gene repression is functionally involved in recruiting DNA replication origins to transcriptionally silent genes. We show that the absence of EZH2, the catalytic subunit of PRC2, results in increased DNA replication initiation, specifically in the vicinity of EZH2 binding sites. The increase in DNA replication initiation does not correlate with transcriptional de-repression or the acquisition of activating histone marks but does correlate with loss of H3K27me3 from bivalent promoters.

Detection and characterization of constitutive replication origins defined by DNA polymerase epsilon

BACKGROUND

Despite the process of DNA replication being mechanistically highly conserved, the location of origins of replication (ORI) may vary from one tissue to the next, or between rounds of replication in eukaryotes, suggesting flexibility in the choice of locations to initiate replication. Lists of human ORI therefore vary widely in number and location, and there are currently no methods available to compare them. Here, we propose a method of detection of ORI based on somatic mutation patterns generated by the mutator phenotype of damaged DNA polymerase epsilon (POLE).

RESULTS

We report the genome-wide localization of constitutive ORI in POLE-mutated human tumors using whole genome sequencing data. Mutations accumulated after many rounds of replication of unsynchronized dividing cell populations in tumors allow to identify constitutive origins, which we show are shared with high fidelity between individuals and tumor types. Using a Smith-Waterman-like dynamic programming approach, we compared replication origin positions obtained from multiple different methods. The comparison allowed us to define a consensus set of replication origins, identified consistently by multiple ORI detection methods. Many DNA features co-localized with the consensus set of ORI, including chromatin loop anchors, G-quadruplexes, S/MARs, and CpGs. Among all features, the H2A.Z histone exhibited the most significant association.

CONCLUSIONS

Our results show that mutation-based detection of replication origins is a viable approach to determining their location and associated sequence features.

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.

Genome-wide measurement of DNA replication fork directionality and quantification of DNA replication initiation and termination with Okazaki fragment sequencing

Studying the dynamics of genome replication in mammalian cells has been historically challenging. To reveal the location of replication initiation and termination in the human genome, we developed Okazaki fragment sequencing (OK-seq), a quantitative approach based on the isolation and strand-specific sequencing of Okazaki fragments, the lagging strand replication intermediates. OK-seq quantitates the proportion of leftward- and rightward-oriented forks at every genomic locus and reveals the location and efficiency of replication initiation and termination events. Here we provide the detailed experimental procedures for performing OK-seq in unperturbed cultured human cells and budding yeast and the bioinformatics pipelines for data processing and computation of replication fork directionality. Furthermore, we present the analytical approach based on a hidden Markov model, which allows automated detection of ascending, descending and flat replication fork directionality segments revealing the zones of replication initiation, termination and unidirectional fork movement across the entire genome. These tools are essential for the accurate interpretation of human and yeast replication programs. The experiments and the data processing can be accomplished within six days. Besides revealing the genome replication program in fine detail, OK-seq has been instrumental in numerous studies unravelling mechanisms of genome stability, epigenome maintenance and genome evolution.

Anomalies in dye-terminator DNA sequencing caused by a natural G-quadruplex

A G-rich DNA sequence from yeast that can form a non-canonical G-quadruplex structure was cloned into a plasmid vector and subjected to Sanger sequencing using dye-labeled dideoxynucleotides. Two different effects were observed. In one, presence of the G4 sequence on the template strand led to incorrect incorporation of an A residue at an internal position in the G4 sequence. In the other, the nascent strand caused attenuation of the readout coincident with synthesis of the G-rich DNA. The two effects are novel examples of disruption in DNA synthesis caused by a G4 sequence. These results provide a new example of a DNA structure that could influence genomic stability in human cells.

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 directly regulates the frequency of oncogenic chromosomal translocations

Chromosomal translocations result from the joining of DNA double-strand breaks (DSBs) and frequently cause cancer. However, the steps linking DSB formation to DSB ligation remain undeciphered. We report that DNA replication timing (RT) directly regulates lymphomagenic Myc translocations during antibody maturation in B cells downstream of DSBs and independently of DSB frequency. Depletion of minichromosome maintenance complexes alters replication origin activity, decreases translocations, and deregulates global RT. Ablating a single origin at Myc causes an early-to-late RT switch, loss of translocations, and reduced proximity with the immunoglobulin heavy chain ( Igh ) gene, its major translocation partner. These phenotypes were reversed by restoring early RT. Disruption of early RT also reduced tumorigenic translocations in human leukemic cells. Thus, RT constitutes a general mechanism in translocation biogenesis linking DSB formation to DSB ligation.

Determination of human DNA replication origin position and efficiency reveals principles of initiation zone organisation

Replication of the human genome initiates within broad zones of ∼150 kb. The extent to which firing of individual DNA replication origins within initiation zones is spatially stochastic or localised at defined sites remains a matter of debate. A thorough characterisation of the dynamic activation of origins within initiation zones is hampered by the lack of a high-resolution map of both their position and efficiency. To address this shortcoming, we describe a modification of initiation site sequencing (ini-seq), based on density substitution. Newly replicated DNA is rendered 'heavy-light' (HL) by incorporation of BrdUTP while unreplicated DNA remains 'light-light' (LL). Replicated HL-DNA is separated from unreplicated LL-DNA by equilibrium density gradient centrifugation, then both fractions are subjected to massive parallel sequencing. This allows precise mapping of 23,905 replication origins simultaneously with an assignment of a replication initiation efficiency score to each. We show that origin firing within early initiation zones is not randomly distributed. Rather, origins are arranged hierarchically with a set of very highly efficient origins marking zone boundaries. We propose that these origins explain much of the early firing activity arising within initiation zones, helping to unify the concept of replication initiation zones with the identification of discrete replication origin sites.

Impact of Chromosomal Context on Origin Selection and the Replication Program

Eukaryotic DNA replication is regulated by conserved mechanisms that bring about a spatial and temporal organization in which distinct genomic domains are copied at characteristic times during S phase. Although this replication program has been closely linked with genome architecture, we still do not understand key aspects of how chromosomal context modulates the activity of replication origins. To address this question, we have exploited models that combine engineered genomic rearrangements with the unique replication programs of post-quiescence and pre-meiotic S phases. Our results demonstrate that large-scale inversions surprisingly do not affect cell proliferation and meiotic progression, despite inducing a restructuring of replication domains on each rearranged chromosome. Remarkably, these alterations in the organization of DNA replication are entirely due to changes in the positions of existing origins along the chromosome, as their efficiencies remain virtually unaffected genome wide. However, we identified striking alterations in origin firing proximal to the fusion points of each inversion, suggesting that the immediate chromosomal neighborhood of an origin is a crucial determinant of its activity. Interestingly, the impact of genome reorganization on replication initiation is highly comparable in the post-quiescent and pre-meiotic S phases, despite the differences in DNA metabolism in these two physiological states. Our findings therefore shed new light on how origin selection and the replication program are governed by chromosomal architecture.

Cohesin-mediated loop anchors confine the locations of human replication origins

DNA replication occurs through an intricately regulated series of molecular events and is fundamental for genome stability1,2. At present, it is unknown how the locations of replication origins are determined in the human genome. Here we dissect the role of topologically associating domains (TADs)3-6, subTADs7 and loops8 in the positioning of replication initiation zones (IZs). We stratify TADs and subTADs by the presence of corner-dots indicative of loops and the orientation of CTCF motifs. We find that high-efficiency, early replicating IZs localize to boundaries between adjacent corner-dot TADs anchored by high-density arrays of divergently and convergently oriented CTCF motifs. By contrast, low-efficiency IZs localize to weaker dotless boundaries. Following ablation of cohesin-mediated loop extrusion during G1, high-efficiency IZs become diffuse and delocalized at boundaries with complex CTCF motif orientations. Moreover, G1 knockdown of the cohesin unloading factor WAPL results in gained long-range loops and narrowed localization of IZs at the same boundaries. Finally, targeted deletion or insertion of specific boundaries causes local replication timing shifts consistent with IZ loss or gain, respectively. Our data support a model in which cohesin-mediated loop extrusion and stalling at a subset of genetically encoded TAD and subTAD boundaries is an essential determinant of the locations of replication origins in human S phase.

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.

Promoters are key organizers of the duplication of vertebrate genomes

In vertebrates, single cell analyses of replication timing patterns brought to light a very well controlled program suggesting a tight regulation on initiation sites. Mapping of replication origins with different methods has revealed discrete preferential sites, enriched in promoters and potential G-quadruplex motifs, which can aggregate into initiation zones spanning several tens of kilobases (kb). Another characteristic of replication origins is a nucleosome-free region (NFR). A modified yeast strain containing a humanized origin recognition complex (ORC) fires new origins at NFRs revealing their regulatory role. In cooperation with NFRs, the histone variant H2A.Z facilitates ORC loading through di-methylation of lysine 20 of histone H4. Recent studies using genome editing methods show that efficient initiation sites associated with transcriptional activity can synergize over several tens of kb by establishing physical contacts and lead to the formation of early domains of DNA replication demonstrating a co-regulation between replication initiation and transcription.

MTA2 sensitizes gastric cancer cells to PARP inhibition by induction of DNA replication stress

Poly (ADP-ribose) polymerase (PARP) inhibitor olaparib selectively kills cancer cells with BRCA-deficiency and is approved for BRCA-mutated breast, ovarian and pancreatic cancers by FDA. However, phase III study of olaparib failed to show a significant improvement in overall survival in patients with gastric cancer (GC). To discover an effective biomarker for GC patient-selection in olaparib treatment, we analyzed proteomic profiling of 12 GC cell lines. MTA2 was identified to confer sensitivity to olaparib by aggravating olaparib-induced replication stress in cancer cells. Mechanistically, we applied Cleavage Under Targets and Tagmentation assay to find that MTA2 proteins preferentially bind regions of replication origin-associated DNA sequences, which could be enhanced by olaparib treatment. Furthermore, MTA2 was validated here to render cancer cells susceptible to combination of olaparib with ATR inhibitor AZD6738. In general, our study identified MTA2 as a potential biomarker for olaparib sensitivity by aggravating olaparib-induced replication stress.

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.

Replication initiation: Implications in genome integrity

In every cell cycle, billions of nucleotides need to be duplicated within hours, with extraordinary precision and accuracy. The molecular mechanism by which cells regulate the replication event is very complicated, and the entire process begins way before the onset of S phase. During the G1 phase of the cell cycle, cells prepare by assembling essential replication factors to establish the pre-replicative complex at origins, sites that dictate where replication would initiate during S phase. During S phase, the replication process is tightly coupled with the DNA repair system to ensure the fidelity of replication. Defects in replication and any error must be recognized by DNA damage response and checkpoint signaling pathways in order to halt the cell cycle before cells are allowed to divide. The coordination of these processes throughout the cell cycle is therefore critical to achieve genomic integrity and prevent diseases. In this review, we focus on the current understanding of how the replication initiation events are regulated to achieve genome stability.

Binding of the Treslin-MTBP Complex to Specific Regions of the Human Genome Promotes the Initiation of DNA Replication

The processes that control where higher eukaryotic cells initiate DNA replication throughout the genome are not understood clearly. In metazoans, the Treslin-MTBP complex mediates critical final steps in formation of the activated replicative helicase prior to initiation of replication. Here, we map the genome-wide distribution of the MTBP subunit of this complex in human cells. Our results indicate that MTBP binds to at least 30,000 sites in the genome. A majority of these sites reside in regions of open chromatin that contain transcriptional-regulatory elements (e.g., promoters, enhancers, and super-enhancers), which are known to be preferred areas for initiation of replication. Furthermore, many binding sites encompass two genomic features: a nucleosome-free DNA sequence (e.g., G-quadruplex DNA or AP-1 motif) and a nucleosome bearing histone marks characteristic of open chromatin, such as H3K4me2. Taken together, these findings indicate that Treslin-MTBP associates coordinately with multiple genomic signals to promote initiation of replication.

Kumagai and Dunphy show that Treslin-MTBP, activator of the replicative helicase, binds to at least 30,000 sites in the human genome. Many sites contain a nucleosome with active chromatin marks and nucleosome-free DNA (G-quadruplex or AP-1 site). Thus, Treslin-MTBP associates with multiple genomic elements to promote initiation of DNA replication.

Genome-wide analysis of DNA replication and DNA double-strand breaks using TrAEL-seq

Faithful replication of the entire genome requires replication forks to copy large contiguous tracts of DNA, and sites of persistent replication fork stalling present a major threat to genome stability. Understanding the distribution of sites at which replication forks stall, and the ensuing fork processing events, requires genome-wide methods that profile replication fork position and the formation of recombinogenic DNA ends. Here, we describe Transferase-Activated End Ligation sequencing (TrAEL-seq), a method that captures single-stranded DNA 3' ends genome-wide and with base pair resolution. TrAEL-seq labels both DNA breaks and replication forks, providing genome-wide maps of replication fork progression and fork stalling sites in yeast and mammalian cells. Replication maps are similar to those obtained by Okazaki fragment sequencing; however, TrAEL-seq is performed on asynchronous populations of wild-type cells without incorporation of labels, cell sorting, or biochemical purification of replication intermediates, rendering TrAEL-seq far simpler and more widely applicable than existing replication fork direction profiling methods. The specificity of TrAEL-seq for DNA 3' ends also allows accurate detection of double-strand break sites after the initiation of DNA end resection, which we demonstrate by genome-wide mapping of meiotic double-strand break hotspots in a dmc1Δ mutant that is competent for end resection but not strand invasion. Overall, TrAEL-seq provides a flexible and robust methodology with high sensitivity and resolution for studying DNA replication and repair, which will be of significant use in determining mechanisms of genome instability.

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.

A predictable conserved DNA base composition signature defines human core DNA replication origins

DNA replication initiates from multiple genomic locations called replication origins. In metazoa, DNA sequence elements involved in origin specification remain elusive. Here, we examine pluripotent, primary, differentiating, and immortalized human cells, and demonstrate that a class of origins, termed core origins, is shared by different cell types and host ~80% of all DNA replication initiation events in any cell population. We detect a shared G-rich DNA sequence signature that coincides with most core origins in both human and mouse genomes. Transcription and G-rich elements can independently associate with replication origin activity. Computational algorithms show that core origins can be predicted, based solely on DNA sequence patterns but not on consensus motifs. Our results demonstrate that, despite an attributed stochasticity, core origins are chosen from a limited pool of genomic regions. Immortalization through oncogenic gene expression, but not normal cellular differentiation, results in increased stochastic firing from heterochromatin and decreased origin density at TAD borders.

In metazoan the DNA sequence elements characterizing origin specification are unknown. By generating and analysing 19 SNS-seq datasets from different human cell types, the authors reveal a class and features of Core origins of replication which can be predicted by an algorithm.

Arabidopsis DNA Replication Initiates in Intergenic, AT-Rich Open Chromatin

The selection and firing of DNA replication origins play key roles in ensuring that eukaryotes accurately replicate their genomes. This process is not well documented in plants due in large measure to difficulties in working with plant systems. We developed a new functional assay to label and map very early replicating loci that must, by definition, include at least a subset of replication origins. Arabidopsis (Arabidopsis thaliana) cells were briefly labeled with 5-ethynyl-2'-deoxy-uridine, and nuclei were subjected to two-parameter flow sorting. We identified more than 5500 loci as initiation regions (IRs), the first regions to replicate in very early S phase. These were classified as strong or weak IRs based on the strength of their replication signals. Strong initiation regions were evenly spaced along chromosomal arms and depleted in centromeres, while weak initiation regions were enriched in centromeric regions. IRs are AT-rich sequences flanked by more GC-rich regions and located predominantly in intergenic regions. Nuclease sensitivity assays indicated that IRs are associated with accessible chromatin. Based on these observations, initiation of plant DNA replication shows some similarity to, but is also distinct from, initiation in other well-studied eukaryotic systems.

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.

Flying High-Muscle-Specific Underreplication in Drosophila

Drosophila underreplicate the DNA of thoracic nuclei, stalling during S phase at a point that is proportional to the total genome size in each species. In polytene tissues, such as the Drosophila salivary glands, all of the nuclei initiate multiple rounds of DNA synthesis and underreplicate. Yet, only half of the nuclei isolated from the thorax stall; the other half do not initiate S phase. Our question was, why half? To address this question, we use flow cytometry to compare underreplication phenotypes between thoracic tissues. When individual thoracic tissues are dissected and the proportion of stalled DNA synthesis is scored in each tissue type, we find that underreplication occurs in the indirect flight muscle, with the majority of underreplicated nuclei in the dorsal longitudinal muscles (DLM). Half of the DNA in the DLM nuclei stall at S phase between the unreplicated G0 and fully replicated G1. The dorsal ventral flight muscle provides the other source of underreplication, and yet, there, the replication stall point is earlier (less DNA replicated), and the endocycle is initiated. The differences in underreplication and ploidy in the indirect flight muscles provide a new tool to study heterochromatin, underreplication and endocycle control.

Evolution of replication origins in vertebrate genomes: rapid turnover despite selective constraints

The replication program of vertebrate genomes is driven by the chromosomal distribution and timing of activation of tens of thousands of replication origins. Genome-wide studies have shown the association of origins with promoters and CpG islands, and their enrichment in G-quadruplex motifs (G4). However, the genetic determinants driving their activity remain poorly understood. To gain insight on the constraints operating on origins, we conducted the first evolutionary comparison of origins across vertebrates. We generated a genome-wide map of chicken origins (the first of a bird genome), and performed a comparison with human and mouse maps. The analysis of intra-species polymorphism revealed a strong depletion of genetic diversity at the core of replication initiation loci. This depletion is not linked to the presence of G4 motifs, promoters or CpG islands. In contrast, we show that origins experienced a rapid turnover during vertebrate evolution, since pairwise comparisons of origin maps revealed that <24% of them are conserved among vertebrates. This study unravels the existence of a novel determinant of origins, the precise functional role of which remains to be determined. Despite the importance of replication initiation for the fitness of organisms, the distribution of origins along vertebrate chromosomes is highly flexible.

RecQL4 tethering on the pre-replicative complex induces unscheduled origin activation and replication stress in human cells

Sequential activation of DNA replication origins is precisely programmed and critical to maintaining genome stability. RecQL4, a member of the conserved RecQ family of helicases, plays an essential role in the initiation of DNA replication in mammalian cells. Here, we showed that RecQL4 protein tethered on the pre-replicative complex (pre-RC) induces early activation of late replicating origins during S phase. Tethering of RecQL4 or its N terminus on pre-RCs via fusion with Orc4 protein resulted in the recruitment of essential initiation factors, such as Mcm10, And-1, Cdc45, and GINS, increasing nascent DNA synthesis in late replicating origins during early S phase. In this origin activation process, tethered RecQL4 was able to recruit Cdc45 even in the absence of cyclin-dependent kinase (CDK) activity, whereas CDK phosphorylation of RecQL4 N terminus was required for interaction with and origin recruitment of And-1 and GINS. In addition, forced activation of replication origins by RecQL4 tethering resulted in increased replication stress and the accumulation of ssDNAs, which can be recovered by transcription inhibition. Collectively, these results suggest that recruitment of RecQL4 to replication origins is an important step for temporal activation of replication origins during S phase. Further, perturbation of replication timing control by unscheduled origin activation significantly induces replication stress, which is mostly caused by transcription-replication conflicts.

Origins of DNA replication

In all kingdoms of life, DNA is used to encode hereditary information. Propagation of the genetic material between generations requires timely and accurate duplication of DNA by semiconservative replication prior to cell division to ensure each daughter cell receives the full complement of chromosomes. DNA synthesis of daughter strands starts at discrete sites, termed replication origins, and proceeds in a bidirectional manner until all genomic DNA is replicated. Despite the fundamental nature of these events, organisms have evolved surprisingly divergent strategies that control replication onset. Here, we discuss commonalities and differences in replication origin organization and recognition in the three domains of life.

5-hydroxymethylcytosine Marks Mammalian Origins Acting as a Barrier to Replication

In most mammalian cells, DNA replication occurs once, and only once between cell divisions. Replication initiation is a highly regulated process with redundant mechanisms that prevent errant initiation events. In lower eukaryotes, replication is initiated from a defined consensus sequence, whereas a consensus sequence delineating mammalian origin of replication has not been identified. Here we show that 5-hydroxymethylcytosine (5hmC) is present at mammalian replication origins. Our data support the hypothesis that 5hmC has a role in cell cycle regulation. We show that 5hmC level is inversely proportional to proliferation; indeed, 5hmC negatively influences cell division by increasing the time a cell resides in G1. Our data suggest that 5hmC recruits replication-licensing factors, then is removed prior to or during origin firing. Later we propose that TET2, the enzyme catalyzing 5mC to 5hmC conversion, acts as barrier to rereplication. In a broader context, our results significantly advance the understating of 5hmC involvement in cell proliferation and disease states.

Involvement of G-quadruplex regions in mammalian replication origin activity

Genome-wide studies of DNA replication origins revealed that origins preferentially associate with an Origin G-rich Repeated Element (OGRE), potentially forming G-quadruplexes (G4). Here, we functionally address their requirements for DNA replication initiation in a series of independent approaches. Deletion of the OGRE/G4 sequence strongly decreased the corresponding origin activity. Conversely, the insertion of an OGRE/G4 element created a new replication origin. This element also promoted replication of episomal EBV vectors lacking the viral origin, but not if the OGRE/G4 sequence was deleted. A potent G4 ligand, PhenDC3, stabilized G4s but did not alter the global origin activity. However, a set of new, G4-associated origins was created, whereas suppressed origins were largely G4-free. In vitro Xenopus laevis replication systems showed that OGRE/G4 sequences are involved in the activation of DNA replication, but not in the pre-replication complex formation. Altogether, these results converge to the functional importance of OGRE/G4 elements in DNA replication initiation.

Genome-wide mapping and profiling of γH2AX binding hotspots in response to different replication stress inducers

BACKGROUND

Replication stress (RS) gives rise to DNA damage that threatens genome stability. RS can originate from different sources that stall replication by diverse mechanisms. However, the mechanism underlying how different types of RS contribute to genome instability is unclear, in part due to the poor understanding of the distribution and characteristics of damage sites induced by different RS mechanisms.

RESULTS

We use ChIP-seq to map γH2AX binding sites genome-wide caused by aphidicolin (APH), hydroxyurea (HU), and methyl methanesulfonate (MMS) treatments in human lymphocyte cells. Mapping of γH2AX ChIP-seq reveals that APH, HU, and MMS treatments induce non-random γH2AX chromatin binding at discrete regions, suggesting that there are γH2AX binding hotspots in the genome. Characterization of the distribution and sequence/epigenetic features of γH2AX binding sites reveals that the three treatments induce γH2AX binding at largely non-overlapping regions, suggesting that RS may cause damage at specific genomic loci in a manner dependent on the fork stalling mechanism. Nonetheless, γH2AX binding sites induced by the three treatments share common features including compact chromatin, coinciding with larger-than-average genes, and depletion of CpG islands and transcription start sites. Moreover, we observe significant enrichment of SINEs in γH2AX sites in all treatments, indicating that SINEs may be a common barrier for replication polymerases.

CONCLUSIONS

Our results identify the location and common features of genome instability hotspots induced by different types of RS, and help in deciphering the mechanisms underlying RS-induced genetic diseases and carcinogenesis.

Replication dynamics of individual loci in single living cells reveal changes in the degree of replication stochasticity through S phase

Eukaryotic genomes are replicated under the control of a highly sophisticated program during the restricted time period corresponding to S phase. The most widely used replication timing assays, which are performed on populations of millions of cells, suggest that most of the genome is synchronously replicated on homologous chromosomes. We investigated the stochastic nature of this temporal program, by comparing the precise replication times of allelic loci within single vertebrate cells progressing through S phase at six loci replicated from very early to very late. We show that replication timing is strictly controlled for the three loci replicated in the first half of S phase. Out of the three loci replicated in the second part of S phase, two present a significantly more stochastic pattern. Surprisingly, we find that the locus replicated at the very end of S phase, presents stochasticity similar to those replicated in early S phase. We suggest that the richness of loci in efficient origins of replication, which decreases from early- to late-replicating regions, and the strength of interaction with the nuclear lamina may underlie the variation of timing control during S phase.

Replication timing in Drosophila and its peculiarities in polytene chromosomes

Drosophila melanogaster is one of the popular model organisms in DNA replication studies. Since the 1960s, DNA replication of polytene chromosomes has been extensively studied by cytological methods. In the recent two decades, the progress in our understanding of DNA replication was associated with new techniques. Use of fluorescent dyes increased the resolution of cytological methods significantly. High-throughput methods allowed analysis of DNA replication on a genome scale, as well as its correlation with chromatin structure and gene activi ty. Precise mapping of the cytological structures of polytene chromosomes to the genome assembly allowed comparison of replication between polytene chromosomes and chromosomes of diploid cells. New features of replication characteristic for D. melanogaster were described for both diploid and polytene chromosomes. Comparison of genomic replication profiles revealed a significant similarity between Drosophila and other well-studi ed eukaryotic species, such as human. Early replication is often confined to intensely transcribed gene-dense regions characterized by multiple replication initiation sites. Features of DNA replication in Drosophila might be explained by a compact genome. The organization of replication in polytene chromosomes has much in common with the organization of replication in chromosomes in diploid cells. The most important feature of replication in polytene chromosomes is its low rate and the dependence of S-phase duration on many factors: external and internal, local and global. The speed of replication forks in D. melanogaster polytene chromosomes is affected by SUUR and Rif1 proteins. It is not known yet how universal the mechanisms associated with these factors are, but their study is very promising.

Developing Novel G-Quadruplex Ligands: from Interaction with Nucleic Acids to Interfering with Nucleic Acid⁻Protein Interaction

G-quadruplex is a special secondary structure of nucleic acids in guanine-rich sequences of genome. G-quadruplexes have been proved to be involved in the regulation of replication, DNA damage repair, and transcription and translation of oncogenes or other cancer-related genes. Therefore, targeting G-quadruplexes has become a novel promising anti-tumor strategy. Different kinds of small molecules targeting the G-quadruplexes have been designed, synthesized, and identified as potential anti-tumor agents, including molecules directly bind to the G-quadruplex and molecules interfering with the binding between the G-quadruplex structures and related binding proteins. This review will explore the feasibility of G-quadruplex ligands acting as anti-tumor drugs, from basis to application. Meanwhile, since helicase is the most well-defined G-quadruplex-related protein, the most extensive research on the relationship between helicase and G-quadruplexes, and its meaning in drug design, is emphasized.

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.

DNA Topoisomerase I differentially modulates R-loops across the human genome

BACKGROUND

Co-transcriptional R-loops are abundant non-B DNA structures in mammalian genomes. DNA Topoisomerase I (Top1) is often thought to regulate R-loop formation owing to its ability to resolve both positive and negative supercoils. How Top1 regulates R-loop structures at a global level is unknown.

RESULTS

Here, we perform high-resolution strand-specific R-loop mapping in human cells depleted for Top1 and find that Top1 depletion results in both R-loop gains and losses at thousands of transcribed loci, delineating two distinct gene classes. R-loop gains are characteristic for long, highly transcribed, genes located in gene-poor regions anchored to Lamin B1 domains and in proximity to H3K9me3-marked heterochromatic patches. R-loop losses, by contrast, occur in gene-rich regions overlapping H3K27me3-marked active replication initiation regions. Interestingly, Top1 depletion coincides with a block of the cell cycle in G0/G1 phase and a trend towards replication delay.

CONCLUSIONS

Our findings reveal new properties of Top1 in regulating R-loop homeostasis in a context-dependent manner and suggest a potential role for Top1 in modulating the replication process via R-loop formation.

Eukaryotic DNA replication: Orchestrated action of multi-subunit protein complexes

Genome duplication is an essential process to preserve genetic information between generations. The eukaryotic cell cycle is composed of functionally distinct phases: G1, S, G2, and M. One of the key replicative proteins that participate at every stage of DNA replication is the Mcm2-7 complex, a replicative helicase. In the G1 phase, inactive Mcm2-7 complexes are loaded on the replication origins by replication-initiator proteins, ORC and Cdc6. Two kinases, S-CDK and DDK, convert the inactive origin-loaded Mcm2-7 complex to an active helicase, the CMG complex in the S phase. The activated CMG complex begins DNA unwinding and recruits enzymes essential for DNA synthesis to assemble a replisome at the replication fork. After completion of DNA synthesis, the inactive CMG complex on the replicated DNA is removed from chromatin to terminate DNA replication. In this review, we will discuss the structure, function, and regulation of the molecular machines involved in each step of DNA replication.

Conservation and Variation in Strategies for DNA Replication of Kinetoplastid Nuclear Genomes

Introduction:

Understanding how the nuclear genome of kinetoplastid parasites is replicated received experimental stimulus from sequencing of the Leishmania major, Trypanosoma brucei and Trypanosoma cruzi genomes around 10 years ago. Gene annotations suggested key players in DNA replication initiation could not be found in these organisms, despite considerable conservation amongst characterised eukaryotes. Initial studies that indicated trypanosomatids might possess an archaeal-like Origin Recognition Complex (ORC), composed of only a single factor termed ORC1/CDC6, have been supplanted by the more recent identification of an ORC in T. brucei. However, the constituent subunits of T. brucei ORC are highly diverged relative to other eukaryotic ORCs and the activity of the complex appears subject to novel, positive regulation. The availability of whole genome sequences has also allowed the deployment of genome-wide strategies to map DNA replication dynamics, to date in T. brucei and Leishmania. ORC1/CDC6 binding and function in T. brucei displays pronounced overlap with the unconventional organisation of gene expression in the genome. Moreover, mapping of sites of replication initiation suggests pronounced differences in replication dynamics in Leishmania relative to T. brucei.

Conclusion:

Here we discuss what implications these emerging data may have for parasite and eukaryotic biology of DNA replication.

[G-quadruplex: key controllers of human genome duplication]

The correct duplication of the human genome is under the control of a spatiotemporal program that determines where and when replication forks start. This regulation thus mainly operates on replication start sites named replication origins. During the S-phase, about 50 000 origins fire in one human cell. However, the normal or perturbed progression of replication forks also strongly impacts on replication. Recently, several studies have put forward the role of a noncanonical DNA structure, the G-quadruplex, in the control of genome duplication. In this review, we describe the major impact of this structure on starting points and on the progression of replication forks.

DNA Replication Origins in Immunoglobulin Switch Regions Regulate Class Switch Recombination in an R-Loop-Dependent Manner

Class switch recombination (CSR) at the immunoglobulin heavy chain (IgH) locus generates antibody isotypes. CSR depends on double-strand breaks (DSBs) induced by activation-induced cytidine deaminase (AID). Although DSB formation and repair machineries are active in G1 phase, efficient CSR is dependent on cell proliferation and S phase entry; however, the underlying mechanisms are obscure. Here, we show that efficient CSR requires the replicative helicase, the Mcm complex. Mcm proteins are enriched at IgH switch regions during CSR, leading to assembly of facultative replication origins that require Mcm helicase function for productive CSR. Assembly of CSR-associated origins is facilitated by R loops and promotes the physical proximity (synapsis) of recombining switch regions, which is reduced by R loop inhibition or Mcm complex depletion. Thus, R loops contribute to replication origin specification that promotes DSB resolution in CSR. This suggests a mechanism for the dependence of CSR on S phase and cell division.

R Loops in the Regulation of Antibody Gene Diversification

For nearly three decades, R loops have been closely linked with class switch recombination (CSR), the process that generates antibody isotypes and that occurs via a complex cascade initiated by transcription-coupled mutagenesis in switch recombination sequences. R loops form during transcription of switch recombination sequences in vitro and in vivo, and there is solid evidence that R loops are required for efficient class switching. The classical model of R loops posits that they boost mutation rates by generating stable and long tracts of single-stranded DNA that serve as the substrate for activation induced deaminase (AID), the enzyme that initiates the CSR reaction cascade by co-transcriptionally mutating ssDNA in switch recombination sequences. Though logical and compelling, this model has not been supported by in vivo evidence. Indeed, several reports suggest that R loops may not be involved in recruiting AID activity to switch regions, meaning that R loops probably serve other unanticipated roles in CSR. Here, I review the key findings in this field to date and propose hypotheses that could help towards elucidating the precise function of R loops in CSR.

DNA replication origins-where do we begin?

In this review, Prioleau and MacAlpine summarize recent advances in our understanding of how primary sequence, chromatin environment, and nuclear architecture contribute to the dynamic selection and activation of replication origins across diverse cell types and developmental stages.

For more than three decades, investigators have sought to identify the precise locations where DNA replication initiates in mammalian genomes. The development of molecular and biochemical approaches to identify start sites of DNA replication (origins) based on the presence of defining and characteristic replication intermediates at specific loci led to the identification of only a handful of mammalian replication origins. The limited number of identified origins prevented a comprehensive and exhaustive search for conserved genomic features that were capable of specifying origins of DNA replication. More recently, the adaptation of origin-mapping assays to genome-wide approaches has led to the identification of tens of thousands of replication origins throughout mammalian genomes, providing an unprecedented opportunity to identify both genetic and epigenetic features that define and regulate their distribution and utilization. Here we summarize recent advances in our understanding of how primary sequence, chromatin environment, and nuclear architecture contribute to the dynamic selection and activation of replication origins across diverse cell types and developmental stages.

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.

Mechanisms and regulation of DNA replication initiation in eukaryotes

Cellular DNA replication is initiated through the action of multiprotein complexes that recognize replication start sites in the chromosome (termed origins) and facilitate duplex DNA melting within these regions. In a typical cell cycle, initiation occurs only once per origin and each round of replication is tightly coupled to cell division. To avoid aberrant origin firing and re-replication, eukaryotes tightly regulate two events in the initiation process: loading of the replicative helicase, MCM2-7, onto chromatin by the origin recognition complex (ORC), and subsequent activation of the helicase by its incorporation into a complex known as the CMG. Recent work has begun to reveal the details of an orchestrated and sequential exchange of initiation factors on DNA that give rise to a replication-competent complex, the replisome. Here, we review the molecular mechanisms that underpin eukaryotic DNA replication initiation - from selecting replication start sites to replicative helicase loading and activation - and describe how these events are often distinctly regulated across different eukaryotic model organisms.

Molecular Mechanism for Chromatin Regulation During MCM Loading in Mammalian Cells

DNA replication is a fundamental process required for the accurate and timely duplication of chromosomes. During late mitosis to G1 phase, the MCM2-7 complex is loaded onto chromatin in a manner dependent on ORC, CDC6, and Cdt1, and chromatin becomes licensed for replication. Although every eukaryotic organism shares common features in replication control, there are also some differences among species. For example, in higher eukaryotic cells including human cells, no strict sequence specificity has been observed for replication origins, unlike budding yeast or bacterial replication origins. Therefore, elements other than beyond DNA sequences are important for regulating replication. For example, the stability and precise positioning of nucleosomes affects replication control. However, little is known about how nucleosome structure is regulated when replication licensing occurs. During the last decade, histone acetylation enzyme HBO1, chromatin remodeler SNF2H, and histone chaperone GRWD1 have been identified as chromatin-handling factors involved in the promotion of replication licensing. In this review, we discuss how the rearrangement of nucleosome formation by these factors affects replication licensing.