Processing of eukaryotic Okazaki fragments by redundant nucleases can be uncoupled from ongoing DNA replication in vivo.

Prior to ligation, each Okazaki fragment synthesized on the lagging strand in eukaryotes must be nucleolytically processed. Nuclease cleavage takes place in the context of 5' flap structures generated via strand-displacement synthesis by DNA polymerase delta. At least three DNA nucleases: Rad27 (Fen1), Dna2 and Exo1, have been implicated in processing Okazaki fragment flaps. However, neither the contributions of individual nucleases to lagging-strand synthesis nor the structure of the DNA intermediates formed in their absence have been fully defined in vivo. By conditionally depleting lagging-strand nucleases and directly analyzing Okazaki fragments synthesized in vivo in Saccharomyces cerevisiae, we conduct a systematic evaluation of the impact of Rad27, Dna2 and Exo1 on lagging-strand synthesis. We find that Rad27 processes the majority of lagging-strand flaps, with a significant additional contribution from Exo1 but not from Dna2. When nuclease cleavage is impaired, we observe a reduction in strand-displacement synthesis as opposed to the widespread generation of long Okazaki fragment 5' flaps, as predicted by some models. Further, using cell cycle-restricted constructs, we demonstrate that both the nucleolytic processing and the ligation of Okazaki fragments can be uncoupled from DNA replication and delayed until after synthesis of the majority of the genome is complete.

Developmental and cancer-associated plasticity of DNA replication preferentially targets GC-poor, lowly expressed and late-replicating regions.

The spatiotemporal program of metazoan DNA replication is regulated during development and altered in cancers. We have generated novel OK-seq, Repli-seq and RNA-seq data to compare the DNA replication and gene expression programs of twelve cancer and non-cancer human cell types. Changes in replication fork directionality (RFD) determined by OK-seq are widespread but more frequent within GC-poor isochores and largely disconnected from transcription changes. Cancer cell RFD profiles cluster with non-cancer cells of similar developmental origin but not with different cancer types. Importantly, recurrent RFD changes are detected in specific tumour progression pathways. Using a model for establishment and early progression of chronic myeloid leukemia (CML), we identify 1027 replication initiation zones (IZs) that progressively change efficiency during long-term expression of the BCR-ABL1 oncogene, being twice more often downregulated than upregulated. Prolonged expression of BCR-ABL1 results in targeting of new IZs and accentuation of previous efficiency changes. Targeted IZs are predominantly located in GC-poor, late replicating gene deserts and frequently silenced in late CML. Prolonged expression of BCR-ABL1 results in massive deletion of GC-poor, late replicating DNA sequences enriched in origin silencing events. We conclude that BCR-ABL1 expression progressively affects replication and stability of GC-poor, late-replicating regions during CML progression.

Zac1 regulates an imprinted gene network critically involved in the control of embryonic growth.

Genomic imprinting is an epigenetic mechanism of regulation that restrains the expression of a small subset of mammalian genes to one parental allele. The reason for the targeting of these approximately 80 genes by imprinting remains uncertain. We show that inactivation of the maternally repressed Zac1 transcription factor results in intrauterine growth restriction, altered bone formation, and neonatal lethality. A meta-analysis of microarray data reveals that Zac1 is a member of a network of coregulated genes comprising other imprinted genes involved in the control of embryonic growth. Zac1 alters the expression of several of these imprinted genes, including Igf2, H19, Cdkn1c, and Dlk1, and it directly regulates the Igf2/H19 locus through binding to a shared enhancer. Accordingly, these data identify a network of imprinted genes, including Zac1, which controls embryonic growth and which may be the basis for the implementation of a common mechanism of gene regulation during mammalian evolution.

Monitoring genome-wide replication fork directionality by Okazaki fragment sequencing in mammalian cells.

The ability to monitor DNA replication fork directionality at the genome-wide scale is paramount for a greater understanding of how genetic and environmental perturbations can impact replication dynamics in human cells. Here we describe a detailed protocol for isolating and sequencing Okazaki fragments from asynchronously growing mammalian cells, termed Okazaki fragment sequencing (Ok-seq), for the purpose of quantitatively determining replication initiation and termination frequencies around specific genomic loci by meta-analyses. Briefly, cells are pulsed with 5-ethynyl-2'-deoxyuridine (EdU) to label newly synthesized DNA, and collected for DNA extraction. After size fractionation on a sucrose gradient, Okazaki fragments are concentrated and purified before click chemistry is used to tag the EdU label with a biotin conjugate that is cleavable under mild conditions. Biotinylated Okazaki fragments are then captured on streptavidin beads and ligated to Illumina adapters before library preparation for Illumina sequencing. The use of Ok-seq to interrogate genome-wide replication fork initiation and termination efficiencies can be applied to all unperturbed, asynchronously growing mammalian cells or under conditions of replication stress, and the assay can be performed in less than 2 weeks.

Transcription shapes DNA replication initiation and termination in human cells.

Although DNA replication is a fundamental aspect of biology, it is not known what determines where DNA replication starts and stops in the human genome. We directly identified and quantitatively compared sites of replication initiation and termination in untransformed human cells. We found that replication preferentially initiates at the transcription start site of genes occupied by high levels of RNA polymerase II, and terminates at their polyadenylation sites, thereby ensuring global co-directionality of transcription and replication, particularly at gene 5' ends. During replication stress, replication initiation is stimulated downstream of genes and termination is redistributed to gene bodies; this globally reorients replication relative to transcription around gene 3' ends. These data suggest that replication initiation and termination are coupled to transcription in human cells, and propose a model for the impact of replication stress on genome integrity.

Replication landscape of the human genome.

Despite intense investigation, human replication origins and termini remain elusive. Existing data have shown strong discrepancies. Here we sequenced highly purified Okazaki fragments from two cell types and, for the first time, quantitated replication fork directionality and delineated initiation and termination zones genome-wide. Replication initiates stochastically, primarily within non-transcribed, broad (up to 150 kb) zones that often abut transcribed genes, and terminates dispersively between them. Replication fork progression is significantly co-oriented with the transcription. Initiation and termination zones are frequently contiguous, sometimes separated by regions of unidirectional replication. Initiation zones are enriched in open chromatin and enhancer marks, even when not flanked by genes, and often border 'topologically associating domains' (TADs). Initiation zones are enriched in origin recognition complex (ORC)-binding sites and better align to origins previously mapped using bubble-trap than λ-exonuclease. This novel panorama of replication reveals how chromatin and transcription modulate the initiation process to create cell-type-specific replication programs.

Best practices for mapping replication origins in eukaryotic chromosomes.

Understanding the regulatory principles ensuring complete DNA replication in each cell division is critical for deciphering the mechanisms that maintain genomic stability. Recent advances in genome sequencing technology facilitated complete mapping of DNA replication sites and helped move the field from observing replication patterns at a handful of single loci to analyzing replication patterns genome-wide. These advances address issues, such as the relationship between replication initiation events, transcription, and chromatin modifications, and identify potential replication origin consensus sequences. This unit summarizes the technological and fundamental aspects of replication profiling and briefly discusses novel insights emerging from mining large datasets, published in the last 3 years, and also describes DNA replication dynamics on a whole-genome scale.

Mapping replication origin sequences in eukaryotic chromosomes.

Recent advances in genome-sequencing technology have led to the complete mapping of DNA replication initiation sites in the human genome. This thorough origin mapping facilitates understanding of the relationship between replication initiation events, transcription, and chromatin modifications, and allows the characterization of consensus sequences of potential replication origins. This unit provides a detailed protocol for isolation and sequence analysis of nascent DNA strands. Two variations of the protocol based on non-overlapping assumptions are described below, addressing potential bias issues for whole-genome analyses.

SUMO2/3 modification of cyclin E contributes to the control of replication origin firing.

The small ubiquitin-like modifier (SUMO) pathway is essential for the maintenance of genome stability. We investigated its possible involvement in the control of DNA replication during S phase by using the Xenopus cell-free system. Here we show that the SUMO pathway is critical to limit the number and, thus, the density of replication origins that are activated in early S phase. We identified cyclin E, which regulates cyclin-dependent kinase 2 (Cdk2) to trigger origin firing, as an S-phase substrate of this pathway. We show that cyclin E is dynamically and highly conjugated to SUMO2/3 on chromatin, independently of Cdk2 activity and origin activation. Moreover, cyclin E is the predominant SUMO2/3 target on chromatin in early S phase, as cyclin E depletion abolishes, while its readdition restores, the SUMO2/3 signal. Together, our data indicate that cyclin E SUMOylation is important for controlling origin firing once the cyclin E-Cdk2 complex is recruited onto replication origins.

From simple bacterial and archaeal replicons to replication N/U-domains.

The Replicon Theory proposed 50 years ago has proven to apply for replicons of the three domains of life. Here, we review our knowledge of genome organization into single and multiple replicons in bacteria, archaea and eukarya. Bacterial and archaeal replicator/initiator systems are quite specific and efficient, whereas eukaryotic replicons show degenerate specificity and efficiency, allowing for complex regulation of origin firing time. We expand on recent evidence that ~50% of the human genome is organized as ~1,500 megabase-sized replication domains with a characteristic parabolic (U-shaped) replication timing profile and linear (N-shaped) gradient of replication fork polarity. These N/U-domains correspond to self-interacting segments of the chromatin fiber bordered by open chromatin zones and replicate by cascades of origin firing initiating at their borders and propagating to their center, possibly by fork-stimulated initiation. The conserved occurrence of this replication pattern in the germline of mammals has resulted over evolutionary times in the formation of megabase-sized domains with an N-shaped nucleotide compositional skew profile due to replication-associated mutational asymmetries. Overall, these results reveal an evolutionarily conserved but developmentally plastic organization of replication that is driving mammalian genome evolution.