High-resolution Repli-Seq defines the temporal choreography of initiation, elongation and termination of replication in mammalian cells

Parent-of-origin-specific DNA replication timing is confined to large imprinted regions

Genomic imprinting involves differential DNA methylation and gene expression between homologous paternal and maternal loci. It remains unclear, however, whether DNA replication also shows parent-of-origin-specific patterns at imprinted or other genomic regions. Here, we investigate genome-wide asynchronous DNA replication utilizing uniparental human embryonic stem cells containing either maternal-only (parthenogenetic) or paternal-only (androgenetic) DNA. Four clusters of imprinted genes exhibited differential replication timing based on parent of origin, while the remainder of the genome, 99.82%, showed no significant replication asynchrony between parental origins. Active alleles in imprinted gene clusters replicated earlier than their inactive counterparts. At the Prader-Willi syndrome locus, replication asynchrony spanned virtually the entirety of S phase. Replication asynchrony was carried through differentiation to neuronal precursor cells in a manner consistent with gene expression. This study establishes asynchronous DNA replication as a hallmark of large imprinted gene clusters.

In vivo DNA replication dynamics unveil aging-dependent replication stress

The genome duplication program is affected by multiple factors in vivo, including developmental cues, genotoxic stress, and aging. Here, we monitored DNA replication initiation dynamics in regenerating livers of young and old mice after partial hepatectomy to investigate the impact of aging. In young mice, the origin firing sites were well defined; the majority were located 10-50 kb upstream or downstream of expressed genes, and their position on the genome was conserved in human cells. Old mice displayed the same replication initiation sites, but origin firing was inefficient and accompanied by a replication stress response. Inhibitors of the ATR checkpoint kinase fully restored origin firing efficiency in the old mice but at the expense of an inflammatory response and without significantly enhancing the fraction of hepatocytes entering the cell cycle. These findings unveil aging-dependent replication stress and a crucial role of ATR in mitigating the stress-associated inflammation, a hallmark of aging.

Replication initiation sites and zones in the mammalian genome: Where are they located and how are they defined?

Eukaryotic DNA replication is a tightly controlled process that occurs in two main steps, i.e., licensing and firing, which take place in the G1 and S phases of the cell cycle, respectively. In Saccharomyces cerevisiae, the budding yeast, replication origins contain consensus sequences that are recognized and bound by the licensing factor Orc1-6, which then recruits the replicative Mcm2-7 helicase. By contrast, mammalian initiation sites lack such consensus sequences, and the mammalian ORC does not exhibit sequence specificity. Studies performed over the past decades have identified replication initiation sites in the mammalian genome using sequencing-based assays, raising the question of whether replication initiation occurs at confined sites or in broad zones across the genome. Although recent reports have shown that the licensed MCMs in mammalian cells are broadly distributed, suggesting that ORC-dependent licensing may not determine the initiation sites/zones, they are predominantly located upstream of actively transcribed genes. This review compares the mechanism of replication initiation in yeast and mammalian cells, summarizes the sequencing-based technologies used for the identification of initiation sites/zones, and proposes a possible mechanism of initiation-site/zone selection in mammalian cells. Future directions and challenges in this field are also discussed.

Bulk synthesis and beyond: The roles of eukaryotic replicative DNA polymerases

An organism's genomic DNA must be accurately duplicated during each cell cycle. DNA synthesis is catalysed by DNA polymerase enzymes, which extend nucleotide polymers in a 5' to 3' direction. This inherent directionality necessitates that one strand is synthesised forwards (leading), while the other is synthesised backwards discontinuously (lagging) to couple synthesis to the unwinding of duplex DNA. Eukaryotic cells possess many diverse polymerases that coordinate to replicate DNA, with the three main replicative polymerases being Pol α, Pol δ and Pol ε. Studies conducted in yeasts and human cells utilising mutant polymerases that incorporate molecular signatures into nascent DNA implicate Pol ε in leading strand synthesis and Pol α and Pol δ in lagging strand replication. Recent structural insights have revealed how the spatial organization of these enzymes around the core helicase facilitates their strand-specific roles. However, various challenging situations during replication require flexibility in the usage of these enzymes, such as during replication initiation or encounters with replication-blocking adducts. This review summarises the roles of the replicative polymerases in bulk DNA replication and explores their flexible and dynamic deployment to complete genome replication. We also examine how polymerase usage patterns can inform our understanding of global replication dynamics by revealing replication fork directionality to identify regions of replication initiation and termination.

Beyond A and B Compartments: how major nuclear locales define nuclear genome organization and function

Models of nuclear genome organization often propose a binary division into active versus inactive compartments, yet they overlook nuclear bodies. Here we integrated analysis of sequencing and image-based data to compare genome organization in four human cell types relative to three different nuclear locales: the nuclear lamina, nuclear speckles, and nucleoli. Whereas gene expression correlates mostly with nuclear speckle proximity, DNA replication timing correlates with proximity to multiple nuclear locales. Speckle attachment regions emerge as DNA replication initiation zones whose replication timing and gene composition vary with their attachment frequency. Most facultative LADs retain a partially repressed state as iLADs, despite their positioning in the nuclear interior. Knock out of two lamina proteins, Lamin A and LBR, causes a shift of H3K9me3-enriched LADs from lamina to nucleolus, and a reciprocal relocation of H3K27me3-enriched partially repressed iLADs from nucleolus to lamina. Thus, these partially repressed iLADs appear to compete with LADs for nuclear lamina attachment with consequences for replication timing. The nuclear organization in adherent cells is polarized with nuclear bodies and genomic regions segregating both radially and relative to the equatorial plane. Together, our results underscore the importance of considering genome organization relative to nuclear locales for a more complete understanding of the spatial and functional organization of the human genome.

Genome-wide identification of replication fork stalling/pausing sites and the interplay between RNA Pol II transcription and DNA replication progression

BACKGROUND

DNA replication progression can be affected by the presence of physical barriers like the RNA polymerases, leading to replication stress and DNA damage. Nonetheless, we do not know how transcription influences overall DNA replication progression.

RESULTS

To characterize sites where DNA replication forks stall and pause, we establish a genome-wide approach to identify them. This approach uses multiple timepoints during S-phase to identify replication fork/stalling hotspots as replication progresses through the genome. These sites are typically associated with increased DNA damage, overlapped with fragile sites and with breakpoints of rearrangements identified in cancers but do not overlap with replication origins. Overlaying these sites with a genome-wide analysis of RNA polymerase II transcription, we find that replication fork stalling/pausing sites inside genes are directly related to transcription progression and activity. Indeed, we find that slowing down transcription elongation slows down directly replication progression through genes. This indicates that transcription and replication can coexist over the same regions. Importantly, rearrangements found in cancers overlapping transcription-replication collision sites are detected in non-transformed cells and increase following treatment with ATM and ATR inhibitors. At the same time, we find instances where transcription activity favors replication progression because it reduces histone density.

CONCLUSIONS

Altogether, our findings highlight how transcription and replication overlap during S-phase, with both positive and negative consequences for replication fork progression and genome stability by the coexistence of these two processes.

Linear interaction between replication and transcription shapes DNA break dynamics at recurrent DNA break Clusters

Recurrent DNA break clusters (RDCs) are replication-transcription collision hotspots; many are unique to neural progenitor cells. Through high-resolution replication sequencing and a capture-ligation assay in mouse neural progenitor cells experiencing replication stress, we unravel the replication features dictating RDC location and orientation. Most RDCs occur at the replication forks traversing timing transition regions (TTRs), where sparse replication origins connect unidirectional forks. Leftward-moving forks generate telomere-connected DNA double-strand breaks (DSBs), while rightward-moving forks lead to centromere-connected DSBs. Strand-specific mapping for DNA-bound RNA reveals co-transcriptional dual-strand DNA:RNA hybrids present at a higher density in RDC than in other actively transcribed long genes. In addition, mapping RNA polymerase activity uncovers that head-to-head interactions between replication and transcription machinery result in 60% DSB contribution to the head-on compared to 40% for co-directional. Taken together we reveal TTR as a fragile class and show how the linear interaction between transcription and replication impacts genome stability.

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Nucleation and spreading maintain Polycomb domains every cell cycle

Gene repression by the Polycomb pathway is essential for metazoan development. Polycomb domains, characterized by trimethylation of histone H3 lysine 27 (H3K27me3), carry the memory of repression and hence need to be maintained to counter the dilution of parental H3K27me3 with unmodified H3 during replication. Yet, how locus-specific H3K27me3 is maintained through replication is unclear. To understand H3K27me3 recovery post-replication, we first define nucleation sites within each Polycomb domain in mouse embryonic stem cells. To map dynamics of H3K27me3 domains across the cell cycle, we develop CUT&Flow (coupling cleavage under target and tagmentation with flow cytometry). We show that post-replication recovery of Polycomb domains occurs by nucleation and spreading, using the same nucleation sites used during de novo domain formation. By using Polycomb repressive complex 2 (PRC2) subunit-specific inhibitors, we find that PRC2 targets nucleation sites post-replication independent of pre-existing H3K27me3. Thus, competition between H3K27me3 deposition and nucleosome turnover drives both de novo domain formation and maintenance during every cell cycle.

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 Decade of Discovery-Eukaryotic Replisome Disassembly at Replication Termination

The eukaryotic replicative helicase (CMG complex) is assembled during DNA replication initiation in a highly regulated manner, which is described in depth by other manuscripts in this Issue. During DNA replication, the replicative helicase moves through the chromatin, unwinding DNA and facilitating nascent DNA synthesis by polymerases. Once the duplication of a replicon is complete, the CMG helicase and the remaining components of the replisome need to be removed from the chromatin. Research carried out over the last ten years has produced a breakthrough in our understanding, revealing that replication termination, and more specifically replisome disassembly, is indeed a highly regulated process. This review brings together our current understanding of these processes and highlights elements of the mechanism that are conserved or have undergone divergence throughout evolution. Finally, we discuss events beyond the classic termination of DNA replication in S-phase and go over the known mechanisms of replicative helicase removal from chromatin in these particular situations.

Select EZH2 inhibitors enhance viral mimicry effects of DNMT inhibition through a mechanism involving NFAT:AP-1 signaling

DNA methyltransferase inhibitor (DNMTi) efficacy in solid tumors is limited. Colon cancer cells exposed to DNMTi accumulate lysine-27 trimethylation on histone H3 (H3K27me3). We propose this Enhancer of Zeste Homolog 2 (EZH2)-dependent repressive modification limits DNMTi efficacy. Here, we show that low-dose DNMTi treatment sensitizes colon cancer cells to select EZH2 inhibitors (EZH2is). Integrative epigenomic analysis reveals that DNMTi-induced H3K27me3 accumulates at genomic regions poised with EZH2. Notably, combined EZH2i and DNMTi alters the epigenomic landscape to transcriptionally up-regulate the calcium-induced nuclear factor of activated T cells (NFAT):activating protein 1 (AP-1) signaling pathway. Blocking this pathway limits transcriptional activating effects of these drugs, including transposable element and innate immune response gene expression involved in viral defense. Analysis of primary human colon cancer specimens reveals positive correlations between DNMTi-, innate immune response-, and calcium signaling-associated transcription profiles. Collectively, we show that compensatory EZH2 activity limits DNMTi efficacy in colon cancer and link NFAT:AP-1 signaling to epigenetic therapy-induced viral mimicry.

Linear Interaction Between Replication and Transcription Shapes DNA Break Dynamics at Recurrent DNA Break Clusters

Recurrent DNA break clusters (RDCs) are replication-transcription collision hotspots; many are unique to neural progenitor cells. Through high-resolution replication sequencing and a capture-ligation assay in mouse neural progenitor cells experiencing replication stress, we unraveled the replication features dictating RDC location and orientation. Most RDCs occur at the replication forks traversing timing transition regions (TTRs), where sparse replication origins connect unidirectional forks. Leftward-moving forks generate telomere-connected DNA double-strand breaks (DSBs), while rightward-moving forks lead to centromere-connected DSBs. Strand-specific mapping for DNA-bound RNA revealed co-transcriptional dual-strand DNA:RNA hybrids present at a higher density in RDC than in other actively transcribed long genes. In addition, mapping RNA polymerase activity revealed that head-to-head interactions between replication and transcription machinery resulted in 60% DSB contribution to the head-on compared to 40% for co-directional. Our findings revealed TTR as a novel fragile class and highlighted how the linear interaction between transcription and replication impacts genome stability.

UHRF1 ubiquitin ligase activity supports the maintenance of low-density CpG methylation

The RING E3 ubiquitin ligase UHRF1 is an established cofactor for DNA methylation inheritance. Nucleosomal engagement through histone and DNA interactions directs UHRF1 ubiquitin ligase activity toward lysines on histone H3 tails, creating binding sites for DNMT1 through ubiquitin interacting motifs (UIM1 and UIM2). Here, we profile contributions of UHRF1 and DNMT1 to genome-wide DNA methylation inheritance and dissect specific roles for ubiquitin signaling in this process. We reveal DNA methylation maintenance at low-density CpGs is vulnerable to disruption of UHRF1 ubiquitin ligase activity and DNMT1 ubiquitin reading activity through UIM1. Hypomethylation of low-density CpGs in this manner induces formation of partially methylated domains (PMD), a methylation signature observed across human cancers. Furthermore, disrupting DNMT1 UIM2 function abolishes DNA methylation maintenance. Collectively, we show DNMT1-dependent DNA methylation inheritance is a ubiquitin-regulated process and suggest a disrupted UHRF1-DNMT1 ubiquitin signaling axis contributes to the development of PMDs in human cancers.

Cell Cycle-Dependent TICRR/TRESLIN and MTBP Chromatin Binding Mechanisms and Patterns

The selection of replication origins is a defining characteristic of DNA replication in eukaryotes, yet its mechanism in humans has not been well-defined. In this study, we use Cut&Run to examine genomic binding locations for TICRR/TRESLIN and MTBP, the human orthologs for the yeast DNA replication initiation factors Sld3 and Sld7. We mapped TRESLIN and MTBP binding in HCT116 colorectal cancer cells using asynchronous and G1 synchronized populations. Our data show that TRESLIN and MTBP binding patterns are more defined in a G1 synchronized population compared to asynchronously cycling cells. We also examined whether TRESLIN and MTBP are dependent on one another for binding. Our data suggest MTBP is dependent on TRESLIN for proper association with chromatin during G1 but not S phase. Finally, we asked whether TRESLIN and MTBP binding to chromatin requires licensed origins. Using cell lines with a non-degradable inducible Geminin to inhibit licensing, we show TRESLIN and MTBP binding does not require loaded MCMs. Altogether, our Cut&Run data provides evidence for a chromatin binding mechanism of TRESLIN-MTBP during G1 that is dependent on TRESLIN and does not require interactions with licensed origins.

Cell cycle gene alterations associate with a redistribution of mutation risk across chromosomal domains in human cancers

Mutations in human cells exhibit increased burden in heterochromatic, late DNA replication time (RT) chromosomal domains, with variation in mutation rates between tissues mirroring variation in heterochromatin and RT. We observed that regional mutation risk further varies between individual tumors in a manner independent of cell type, identifying three signatures of domain-scale mutagenesis in >4,000 tumor genomes. The major signature reflects remodeling of heterochromatin and of the RT program domains seen across tumors, tissues and cultured cells, and is robustly linked with higher expression of cell proliferation genes. Regional mutagenesis is associated with loss of activity of the tumor-suppressor genes RB1 and TP53, consistent with their roles in cell cycle control, with distinct mutational patterns generated by the two genes. Loss of regional heterogeneity in mutagenesis is associated with deficiencies in various DNA repair pathways. These mutation risk redistribution processes modify the mutation supply towards important genes, diverting the course of somatic evolution.

The chromatin-associated lncREST ensures effective replication stress response by promoting the assembly of fork signaling factors

Besides the well-characterized protein network involved in the replication stress response, several regulatory RNAs have been shown to play a role in this critical process. However, it has remained elusive whether they act locally at the stressed forks. Here, by investigating the RNAs localizing on chromatin upon replication stress induced by hydroxyurea, we identified a set of lncRNAs upregulated in S-phase and controlled by stress transcription factors. Among them, we demonstrate that the previously uncharacterized lncRNA lncREST (long non-coding RNA REplication STress) is transcriptionally controlled by p53 and localizes at stressed replication forks. LncREST-depleted cells experience sustained replication fork progression and accumulate un-signaled DNA damage. Under replication stress, lncREST interacts with the protein NCL and assists in engaging its interaction with RPA. The loss of lncREST is associated with a reduced NCL-RPA interaction and decreased RPA on chromatin, leading to defective replication stress signaling and accumulation of mitotic defects, resulting in apoptosis and a reduction in tumorigenic potential of cancer cells. These findings uncover the function of a lncRNA in favoring the recruitment of replication proteins to sites of DNA replication.

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

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

Emergence of replication timing during early mammalian development

DNA replication enables genetic inheritance across the kingdoms of life. Replication occurs with a defined temporal order known as the replication timing (RT) programme, leading to organization of the genome into early- or late-replicating regions. RT is cell-type specific, is tightly linked to the three-dimensional nuclear organization of the genome1,2 and is considered an epigenetic fingerprint3. In spite of its importance in maintaining the epigenome4, the developmental regulation of RT in mammals in vivo has not been explored. Here, using single-cell Repli-seq5, we generated genome-wide RT maps of mouse embryos from the zygote to the blastocyst stage. Our data show that RT is initially not well defined but becomes defined progressively from the 4-cell stage, coinciding with strengthening of the A and B compartments. We show that transcription contributes to the precision of the RT programme and that the difference in RT between the A and B compartments depends on RNA polymerase II at zygotic genome activation. Our data indicate that the establishment of nuclear organization precedes the acquisition of defined RT features and primes the partitioning of the genome into early- and late-replicating domains. Our work sheds light on the establishment of the epigenome at the beginning of mammalian development and reveals the organizing principles of genome organization.

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 5 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% were reproducibly obtained in at least 20 independent SNS-seq datasets and contained in initiation zones identified by each of three other techniques (20,250 shared origins), suggesting extensive variability in origin usage and identification. 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, S. 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.

Transcription-Replication Conflicts as a Source of Genome Instability

Transcription and replication both require large macromolecular complexes to act on a DNA template, yet these machineries cannot simultaneously act on the same DNA sequence. Conflicts between the replication and transcription machineries (transcription-replication conflicts, or TRCs) are widespread in both prokaryotes and eukaryotes and have the capacity to both cause DNA damage and compromise complete, faithful replication of the genome. This review will highlight recent studies investigating the genomic locations of TRCs and the mechanisms by which they may be prevented, mitigated, or resolved. We address work from both model organisms and mammalian systems but predominantly focus on multicellular eukaryotes owing to the additional complexities inherent in the coordination of replication and transcription in the context of cell type-specific gene expression and higher-order chromatin organization.

Asymmetric distribution of parental H3K9me3 in S phase silences L1 elements

In eukaryotes, repetitive DNA sequences are transcriptionally silenced through histone H3 lysine 9 trimethylation (H3K9me3). Loss of silencing of the repeat elements leads to genome instability and human diseases, including cancer and ageing1-3. Although the role of H3K9me3 in the establishment and maintenance of heterochromatin silencing has been extensively studied4-6, the pattern and mechanism that underlie the partitioning of parental H3K9me3 at replicating DNA strands are unknown. Here we report that H3K9me3 is preferentially transferred onto the leading strands of replication forks, which occurs predominantly at long interspersed nuclear element (LINE) retrotransposons (also known as LINE-1s or L1s) that are theoretically transcribed in the head-on direction with replication fork movement. Mechanistically, the human silencing hub (HUSH) complex interacts with the leading-strand DNA polymerase Pol ε and contributes to the asymmetric segregation of H3K9me3. Cells deficient in Pol ε subunits (POLE3 and POLE4) or the HUSH complex (MPP8 and TASOR) show compromised H3K9me3 asymmetry and increased LINE expression. Similar results were obtained in cells expressing a MPP8 mutant defective in H3K9me3 binding and in TASOR mutants with reduced interactions with Pol ε. These results reveal an unexpected mechanism whereby the HUSH complex functions with Pol ε to promote asymmetric H3K9me3 distribution at head-on LINEs to suppress their expression in S phase.

Replication stress in mammalian embryo development, differentiation, and reprogramming

Duplicating a genome of 3 billion nucleotides is challenged by a variety of obstacles that can cause replication stress and affect the integrity of the genome. Recent studies show that replication fork slowing and stalling is prevalent in early mammalian development, resulting in genome instability and aneuploidy, and constituting a barrier to development in human reproduction. Genome instability resulting from DNA replication stress is a barrier to the cloning of animals and to the reprogramming of differentiated cells to induced pluripotent stem cells, as well as a barrier to cell transformation. Remarkably, the regions most impacted by replication stress are shared in these different cellular contexts, affecting long genes and flanking intergenic areas. In this review we integrate our knowledge of DNA replication stress in mammalian embryos, in programming, and in reprogramming, and we discuss a potential role for fragile sites in sensing replication stress and restricting cell cycle progression in health and disease.

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

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

TRAIP resolves DNA replication-transcription conflicts during the S-phase of unperturbed cells

Cell division is the basis for the propagation of life and requires accurate duplication of all genetic information. DNA damage created during replication (replication stress) is a major cause of cancer, premature aging and a spectrum of other human disorders. Over the years, TRAIP E3 ubiquitin ligase has been shown to play a role in various cellular processes that govern genome integrity and faultless segregation. TRAIP is essential for cell viability, and mutations in TRAIP ubiquitin ligase activity lead to primordial dwarfism in patients. Here, we have determined the mechanism of inhibition of cell proliferation in TRAIP-depleted cells. We have taken advantage of the auxin induced degron system to rapidly degrade TRAIP within cells and to dissect the importance of various functions of TRAIP in different stages of the cell cycle. We conclude that upon rapid TRAIP degradation, specifically in S-phase, cells cease to proliferate, arrest in G2 stage of the cell cycle and undergo senescence. Our findings reveal that TRAIP works in S-phase to prevent DNA damage at transcription start sites, caused by replication-transcription conflicts.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

Characterization of Unidirectional Replication Forks in the Mouse Genome

Origins of replication are genomic regions in which replication initiates in a bidirectional manner. Recently, a new methodology (origin-derived single-stranded DNA sequencing; ori-SSDS) was developed that allows the detection of replication initiation in a strand-specific manner. Reanalysis of the strand-specific data revealed that 18-33% of the peaks are non-symmetrical, suggesting a single direction of replication. Analysis of replication fork direction data revealed that these are origins of replication in which the replication is paused in one of the directions, probably due to the existence of a replication fork barrier. Analysis of the unidirectional origins revealed a preference of G4 quadruplexes for the blocked leading strand. Taken together, our analysis identified hundreds of genomic locations in which the replication initiates only in one direction, and suggests that G4 quadruplexes may serve as replication fork barriers in such places.

An autonomous mathematical model for the mammalian cell cycle

A mathematical model for the mammalian cell cycle is developed as a system of 13 coupled nonlinear ordinary differential equations. The variables and interactions included in the model are based on detailed consideration of available experimental data. A novel feature of the model is inclusion of cycle tasks such as origin licensing and initiation, nuclear envelope breakdown and kinetochore attachment, and their interactions with controllers (molecular complexes involved in cycle control). Other key features are that the model is autonomous, except for a dependence on external growth factors; the variables are continuous in time, without instantaneous resets at phase boundaries; mechanisms to prevent rereplication are included; and cycle progression is independent of cell size. Eight variables represent cell cycle controllers: the Cyclin D1-Cdk4/6 complex, APC^Cdh1, SCF^βTrCP, Cdc25A, MPF, NuMA, the securin-separase complex, and separase. Five variables represent task completion, with four for the status of origins and one for kinetochore attachment. The model predicts distinct behaviors corresponding to the main phases of the cell cycle, showing that the principal features of the mammalian cell cycle, including restriction point behavior, can be accounted for in a quantitative mechanistic way based on known interactions among cycle controllers and their coupling to tasks. The model is robust to parameter changes, in that cycling is maintained over at least a five-fold range of each parameter when varied individually. The model is suitable for exploring how extracellular factors affect cell cycle progression, including responses to metabolic conditions and to anti-cancer therapies.

Mistimed origin licensing and activation stabilize common fragile sites under tight DNA-replication checkpoint activation

Genome integrity requires replication to be completed before chromosome segregation. The DNA-replication checkpoint (DRC) contributes to this coordination by inhibiting CDK1, which delays mitotic onset. Under-replication of common fragile sites (CFSs), however, escapes surveillance, resulting in mitotic chromosome breaks. Here we asked whether loose DRC activation induced by modest stresses commonly used to destabilize CFSs could explain this leakage. We found that tightening DRC activation or CDK1 inhibition stabilizes CFSs in human cells. Repli-Seq and molecular combing analyses showed a burst of replication initiations implemented in mid S-phase across a subset of late-replicating sequences, including CFSs, while the bulk genome was unaffected. CFS rescue and extra-initiations required CDC6 and CDT1 availability in S-phase, implying that CDK1 inhibition permits mistimed origin licensing and firing. In addition to delaying mitotic onset, tight DRC activation therefore supports replication completion of late origin-poor domains at risk of under-replication, two complementary roles preserving genome stability.

Replication timing and transcriptional control: beyond cause and effect - part IV

Decades of work on the spatiotemporal organization of mammalian DNA replication timing (RT) continues to unveil novel correlations with aspects of transcription and chromatin organization but, until recently, mechanisms regulating RT and the biological significance of the RT program had been indistinct. We now know that the RT program is both influenced by and necessary to maintain chromatin structure, forming an epigenetic positive feedback loop. Moreover, the discovery of specific cis-acting elements regulating mammalian RT at both the domain and the whole-chromosome level has revealed multiple cell-type-specific and developmentally regulated mechanisms of RT control. We review recent evidence for diverse mechanisms employed by different cell types to regulate their RT programs and the biological significance of RT regulation during development.

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.

Elucidating the structure and function of the nucleus-The NIH Common Fund 4D Nucleome program

Genomic architecture appears to play crucial roles in health and a variety of diseases. How nuclear structures reorganize over different timescales is elusive, partly because the tools needed to probe and perturb them are not as advanced as needed by the field. To fill this gap, the National Institutes of Health Common Fund started a program in 2015, called the 4D Nucleome (4DN), with the goal of developing and ultimately applying technologies to interrogate the structure and function of nuclear organization in space and time.

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.

Diverse silent chromatin states modulate genome compartmentalization and loop extrusion barriers

The relationships between chromosomal compartmentalization, chromatin state and function are poorly understood. Here by profiling long-range contact frequencies in HCT116 colon cancer cells, we distinguish three silent chromatin states, comprising two types of heterochromatin and a state enriched for H3K9me2 and H2A.Z that exhibits neutral three-dimensional interaction preferences and which, to our knowledge, has not previously been characterized. We find that heterochromatin marked by H3K9me3, HP1α and HP1β correlates with strong compartmentalization. We demonstrate that disruption of DNA methyltransferase activity greatly remodels genome compartmentalization whereby domains lose H3K9me3-HP1α/β binding and acquire the neutrally interacting state while retaining late replication timing. Furthermore, we show that H3K9me3-HP1α/β heterochromatin is permissive to loop extrusion by cohesin but refractory to CTCF binding. Together, our work reveals a dynamic structural and organizational diversity of the silent portion of the genome and establishes connections between the regulation of chromatin state and chromosome organization, including an interplay between DNA methylation, compartmentalization and loop extrusion.

Computational workflow for integrative analyses of DNA replication timing, epigenomic, and transcriptomic data

Temporal profiling of DNA replication timing (RT) in combination with chromatin modifications, chromatin accessibility, and gene expression provides new insights into the causal relationships between chromatin and RT during cell cycle. Here, we describe a protocol for in-depth integrative computational analyses of Repli-seq, ATAC-seq, RNA-seq, and ChIP-seq or CUT&RUN data for multiple marks at various time points across cell cycle and changes in their interrelationships upon an experimental perturbation (e.g., knockdown or overexpression of a regulatory protein). For complete details on the use and execution of this protocol, please refer to Van Rechem et al. (2021).

Optimized Repli-seq: improved DNA replication timing analysis by next-generation sequencing

The human genome is divided into functional units that replicate at specific times during S-phase. This temporal program is known as replication timing (RT) and is coordinated with the spatial organization of the genome and transcriptional activity. RT is also cell type-specific, dynamically regulated during development, and alterations in RT are observed in multiple diseases. Thus, the precise measure of RT is critical to understand the role of RT in gene function regulation. Distinct methods for assaying the RT program exist; however, conventional methods require thousands of cells as input, prohibiting its applicability to samples with limited cell numbers such as those from disease patients or from early developing embryos. Although single-cell RT analyses have been developed, these methods are low throughput, require generation of numerous libraries, increased sequencing costs, and produce low resolution data. Here, we developed an improved method to measure RT genome-wide that enables high-resolution analysis of low input samples. This method incorporates direct cell sorting into lysis buffer, as well as DNA fragmentation and library preparation in a single tube, resulting in higher yields, increased quality, and reproducibility with decreased costs. We also performed a systematic data processing analysis to provide standardized parameters for RT measurement. This optimized method facilitates RT analysis and will enable its application to a broad range of studies investigating the role of RT in gene expression, nuclear architecture, and disease.

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.

Global landscape of replicative DNA polymerase usage in the human genome

The division of labour among DNA polymerase underlies the accuracy and efficiency of replication. However, the roles of replicative polymerases have not been directly established in human cells. We developed polymerase usage sequencing (Pu-seq) in HCT116 cells and mapped Polε and Polα usage genome wide. The polymerase usage profiles show Polε synthesises the leading strand and Polα contributes mainly to lagging strand synthesis. Combining the Polε and Polα profiles, we accurately predict the genome-wide pattern of fork directionality plus zones of replication initiation and termination. We confirm that transcriptional activity contributes to the pattern of initiation and termination and, by separately analysing the effect of transcription on co-directional and converging forks, demonstrate that coupled DNA synthesis of leading and lagging strands is compromised by transcription in both co-directional and convergent forks. Polymerase uncoupling is particularly evident in the vicinity of large genes, including the two most unstable common fragile sites, FRA3B and FRA3D, thus linking transcription-induced polymerase uncoupling to chromosomal instability. Together, our result demonstrated that Pu-seq in human cells provides a powerful and straightforward methodology to explore DNA polymerase usage and replication fork dynamics.

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.

Epigenetic control of chromosome-associated lncRNA genes essential for replication and stability

ASARs are long noncoding RNA genes that control replication timing of entire human chromosomes in cis. The three known ASAR genes are located on human chromosomes 6 and 15, and are essential for chromosome integrity. To identify ASARs on all human chromosomes we utilize a set of distinctive ASAR characteristics that allow for the identification of hundreds of autosomal loci with epigenetically controlled, allele-restricted behavior in expression and replication timing of coding and noncoding genes, and is distinct from genomic imprinting. Disruption of noncoding RNA genes at five of five tested loci result in chromosome-wide delayed replication and chromosomal instability, validating their ASAR activity. In addition to the three known essential cis-acting chromosomal loci, origins, centromeres, and telomeres, we propose that all mammalian chromosomes also contain "Inactivation/Stability Centers" that display allele-restricted epigenetic regulation of protein coding and noncoding ASAR genes that are essential for replication and stability of each chromosome.

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.

RRM2 enhances MYCN-driven neuroblastoma formation and acts as a synergistic target with CHK1 inhibition

High-risk neuroblastoma, a pediatric tumor originating from the sympathetic nervous system, has a low mutation load but highly recurrent somatic DNA copy number variants. Previously, segmental gains and/or amplifications allowed identification of drivers for neuroblastoma development. Using this approach, combined with gene dosage impact on expression and survival, we identified ribonucleotide reductase subunit M2 (RRM2) as a candidate dependency factor further supported by growth inhibition upon in vitro knockdown and accelerated tumor formation in a neuroblastoma zebrafish model coexpressing human RRM2 with MYCN. Forced RRM2 induction alleviates excessive replicative stress induced by CHK1 inhibition, while high RRM2 expression in human neuroblastomas correlates with high CHK1 activity. MYCN-driven zebrafish tumors with RRM2 co-overexpression exhibit differentially expressed DNA repair genes in keeping with enhanced ATR-CHK1 signaling activity. In vitro, RRM2 inhibition enhances intrinsic replication stress checkpoint addiction. Last, combinatorial RRM2-CHK1 inhibition acts synergistic in high-risk neuroblastoma cell lines and patient-derived xenograft models, illustrating the therapeutic potential.

Single-Cell Mononucleotide Microsatellite Analysis Reveals Differential Insertion-Deletion Dynamics in Mouse T Cells

Microsatellite sequences are particularly prone to slippage during DNA replication, forming insertion-deletion loops that, if left unrepaired, result in de novo mutations (expansions or contractions of the repeat array). Mismatch repair (MMR) is a critical DNA repair mechanism that corrects these insertion-deletion loops, thereby maintaining microsatellite stability. MMR deficiency gives rise to the molecular phenotype known as microsatellite instability (MSI). By sequencing MMR-proficient and -deficient (Mlh1+/+ and Mlh1−/−) single-cell exomes from mouse T cells, we reveal here several previously unrecognized features of in vivo MSI. Specifically, mutational dynamics of insertions and deletions were different on multiple levels. Factors that associated with propensity of mononucleotide microsatellites to insertions versus deletions were: microsatellite length, nucleotide composition of the mononucleotide tract, gene length and transcriptional status, as well replication timing. Here, we show on a single-cell level that deletions — the predominant MSI type in MMR-deficient cells — are preferentially associated with longer A/T tracts, long or transcribed genes and later-replicating genes.

New Era of Mapping and Understanding Common Fragile Sites: An Updated Review on Origin of Chromosome Fragility

Common fragile sites (CFSs) are specific genomic loci prone to forming gaps or breakages upon replication perturbation, which correlate well with chromosomal rearrangement and copy number variation. CFSs have been actively studied due to their important pathophysiological relevance in different diseases such as cancer and neurological disorders. The genetic locations and sequences of CFSs are crucial to understanding the origin of such unstable sites, which require reliable mapping and characterizing approaches. In this review, we will inspect the evolving techniques for CFSs mapping, especially genome-wide mapping and sequencing of CFSs based on current knowledge of CFSs. We will also revisit the well-established hypotheses on the origin of CFSs fragility, incorporating novel findings from the comprehensive analysis of finely mapped CFSs regarding their locations, sequences, and replication/transcription, etc. This review will present the most up-to-date picture of CFSs and, potentially, a new framework for future research of CFSs.

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

Little is known about replication fork velocity variations along eukaryotic genomes, since reference techniques to determine fork speed either provide no sequence information or suffer from low throughput. Here we present NanoForkSpeed, a nanopore sequencing-based method to map and extract the velocity of individual forks detected as tracks of the thymidine analogue bromodeoxyuridine incorporated during a brief pulse-labelling of asynchronously growing cells. NanoForkSpeed retrieves previous Saccharomyces cerevisiae mean fork speed estimates (≈2 kb/min) in the BT1 strain exhibiting highly efficient bromodeoxyuridine incorporation and wild-type growth, and precisely quantifies speed changes in cells with altered replisome progression or exposed to hydroxyurea. The positioning of >125,000 fork velocities provides a genome-wide map of fork progression based on individual fork rates, showing a uniform fork speed across yeast chromosomes except for a marked slowdown at known pausing sites.

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.

Kronos scRT: a uniform framework for single-cell replication timing analysis

Mammalian genomes are replicated in a cell type-specific order and in coordination with transcription and chromatin organization. Currently, single-cell replication studies require individual processing of sorted cells, yielding a limited number (<100) of cells. Here, we develop Kronos scRT, a software for single-cell Replication Timing (scRT) analysis. Kronos scRT does not require a specific platform or cell sorting, which allows investigating large datasets obtained from asynchronous cells. By applying our tool to published data as well as droplet-based single-cell whole-genome sequencing data generated in this study, we exploit scRT from thousands of cells for different mouse and human cell lines. Our results demonstrate that although genomic regions are frequently replicated around their population average RT, replication can occur stochastically throughout S phase. Altogether, Kronos scRT allows fast and comprehensive investigations of the RT programme at the single-cell resolution for both homogeneous and heterogeneous cell populations.

A scalable approach to explore DNA replication in single cells reveals that although aneuploidy does not have a major impact on the pattern of replication, different cell types and sub-populations display distinguished replication paths.

Fiber-Like Organization as a Basic Principle for Euchromatin Higher-Order Structure

A detailed understanding of the principles of the structural organization of genetic material is of great importance for elucidating the mechanisms of differential regulation of genes in development. Modern ideas about the spatial organization of the genome are based on a microscopic analysis of chromatin structure and molecular data on DNA–DNA contact analysis using Chromatin conformation capture (3C) technology, ranging from the “polymer melt” model to a hierarchical folding concept. Heterogeneity of chromatin structure depending on its functional state and cell cycle progression brings another layer of complexity to the interpretation of structural data and requires selective labeling of various transcriptional states under nondestructive conditions. Here, we use a modified approach for replication timing-based metabolic labeling of transcriptionally active chromatin for ultrastructural analysis. The method allows pre-embedding labeling of optimally structurally preserved chromatin, thus making it compatible with various 3D-TEM techniques including electron tomography. By using variable pulse duration, we demonstrate that euchromatic genomic regions adopt a fiber-like higher-order structure of about 200 nm in diameter (chromonema), thus providing support for a hierarchical folding model of chromatin organization as well as the idea of transcription and replication occurring on a highly structured chromatin template.