DNA Replication Timing Enters the Single-Cell Era

Embryonic genome instability upon DNA replication timing program emergence

Faithful DNA replication is essential for genome integrity1-4. Under-replicated DNA leads to defects in chromosome segregation, which are common during embryogenesis5-8. However, the regulation of DNA replication remains poorly understood in early mammalian embryos. Here we constructed a single-cell genome-wide DNA replication atlas of pre-implantation mouse embryos and identified an abrupt replication program switch accompanied by a transient period of genomic instability. In 1- and 2-cell embryos, we observed the complete absence of a replication timing program, and the entire genome replicated gradually and uniformly using extremely slow-moving replication forks. In 4-cell embryos, a somatic-cell-like replication timing program commenced abruptly. However, the fork speed was still slow, S phase was extended, and markers of replication stress, DNA damage and repair increased. This was followed by an increase in break-type chromosome segregation errors specifically during the 4-to-8-cell division with breakpoints enriched in late-replicating regions. These errors were rescued by nucleoside supplementation, which accelerated fork speed and reduced the replication stress. By the 8-cell stage, forks gained speed, S phase was no longer extended and chromosome aberrations decreased. Thus, a transient period of genomic instability exists during normal mouse development, preceded by an S phase lacking coordination between replisome-level regulation and megabase-scale replication timing regulation, implicating a link between their coordination and genome stability.

The dynamical organization of the core pluripotency transcription factors responds to differentiation cues in early S-phase

DNA replication in stem cells is a major challenge for pluripotency preservation and cell fate decisions. This process involves massive changes in the chromatin architecture and the reorganization of many transcription-related molecules in different spatial and temporal scales. Pluripotency is controlled by the master transcription factors (TFs) OCT4, SOX2 and NANOG that partition into condensates in the nucleus of embryonic stem cells. These condensates are proposed to play relevant roles in the regulation of gene expression and the maintenance of pluripotency. Here, we asked whether the dynamical distribution of the pluripotency TFs changes during the cell cycle, particularly during DNA replication. Since the S phase is considered to be a window of opportunity for cell fate decisions, we explored if differentiation cues in G1 phase trigger changes in the distribution of these TFs during the subsequent S phase. Our results show a spatial redistribution of TFs condensates during DNA replication which was not directly related to chromatin compaction. Additionally, fluorescence fluctuation spectroscopy revealed TF-specific, subtle changes in the landscape of TF-chromatin interactions, consistent with their particularities as key players of the pluripotency network. Moreover, we found that differentiation stimuli in the preceding G1 phase triggered a relatively fast and massive reorganization of pluripotency TFs in early-S phase. Particularly, OCT4 and SOX2 condensates dissolved whereas the lifetimes of TF-chromatin interactions increased suggesting that the reorganization of condensates is accompanied with a change in the landscape of TF-chromatin interactions. Notably, NANOG showed impaired interactions with chromatin in stimulated early-S cells in line with its role as naïve pluripotency TF. Together, these findings provide new insights into the regulation of the core pluripotency TFs during DNA replication of embryonic stem cells and highlight their different roles at early differentiation stages.

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.

Rif1-Dependent Control of Replication Timing

Successful duplication of the genome requires the accurate replication of billions of base pairs of DNA within a relatively short time frame. Failure to accurately replicate the genome results in genomic instability and a host of diseases. To faithfully and rapidly replicate the genome, DNA replication must be tightly regulated and coordinated with many other nuclear processes. These regulations, however, must also be flexible as replication kinetics can change through development and differentiation. Exactly how DNA replication is regulated and how this regulation changes through development is an active field of research. One aspect of genome duplication where much remains to be discovered is replication timing (RT), which dictates when each segment of the genome is replicated during S phase. All organisms display some level of RT, yet the precise mechanisms that govern RT remain are not fully understood. The study of Rif1, a protein that actively regulates RT from yeast to humans, provides a key to unlock the underlying molecular mechanisms controlling RT. The paradigm for Rif1 function is to delay helicase activation within certain regions of the genome, causing these regions to replicate late in S phase. Many questions, however, remain about the intricacies of Rif1 function. Here, we review the current models for the activity of Rif1 with the goal of trying to understand how Rif1 functions to establish the RT program.

Mapping Replication Timing in Single Mammalian Cells

Replication timing (RT) is the temporal order in which genomic DNA is replicated during S phase. Early and late replication correlate with transcriptionally active and inactive chromatin compartments, but mechanistic links between large-scale chromosome structure, transcription, and replication are still enigmatic. A proper RT program is necessary to maintain the global epigenome that defines cell identity, suggesting that RT is critical for epigenome integrity by facilitating the assembly of different types of chromatin at different times during S phase. RT is regulated during development and has been found to be altered in disease. Thus, RT can identify stable epigenetic differences distinguishing cell types, and can be used to help stratify patient outcomes and identify markers of disease. Most methods to profile RT require thousands of S-phase cells. In cases where cells are rare (e.g., early-stage embryos or rare primary cell types) or consist of a heterogeneous mixture of cell states (e.g., differentiation intermediates), or when the interest is in determining the degree of stable epigenetic heterogeneity within a population of cells, single-cell measurements of RT are necessary. We have previously developed single cell Repli-seq, a method to measure replication timing in single cells using DNA copy number quantification. To date, however, single-cell Repli-seq suffers from relatively low throughput and high costs. Here, we describe an improved single-cell Repli-seq protocol that uses degenerate oligonucleotide-primed PCR (DOP-PCR) for uniform whole-genome amplification and uniquely barcoded primers that permit early pooling of single-cell samples into a single library preparation. We also provide a bioinformatics platform for analysis of the data. The improved throughput and decreased costs of this method relative to previously published single-cell Repli-seq protocols should make it considerably more accessible to a broad range of investigators. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Whole Genome Amplification (WGA) of single cells and sequence library construction. Basic Protocol 2: Deriving and displaying single-cell replication timing data from whole genome sequencing.

Applications of Single-Cell DNA Sequencing

Over the past decade, genomic analyses of single cells-the fundamental units of life-have become possible. Single-cell DNA sequencing has shed light on biological questions that were previously inaccessible across diverse fields of research, including somatic mutagenesis, organismal development, genome function, and microbiology. Single-cell DNA sequencing also promises significant future biomedical and clinical impact, spanning oncology, fertility, and beyond. While single-cell approaches that profile RNA and protein have greatly expanded our understanding of cellular diversity, many fundamental questions in biology and important biomedical applications require analysis of the DNA of single cells. Here, we review the applications and biological questions for which single-cell DNA sequencing is uniquely suited or required. We include a discussion of the fields that will be impacted by single-cell DNA sequencing as the technology continues to advance.

Chromatin remodeling in replication-uncoupled maintenance DNA methylation and chromosome stability: Insights from ICF syndrome studies

Immunodeficiency, centromeric instability, and facial anomalies (ICF) syndrome is characterized by frequent appearance of multiradial chromosomes, which are distinctive chromosome fusions that occur at hypomethylated pericentromeric regions comprising repetitive sequences, in activated lymphocytes. The syndrome is caused by mutations in DNMT3B, ZBTB24, CDCA7, or HELLS. De novo DNA methylation is likely defective in patients with ICF syndrome harboring mutations in DNMT3B, whereas accumulating evidence suggests that replication‐uncoupled maintenance DNA methylation of late‐replicating regions is impaired in patients with ICF syndrome harboring mutations in ZBTB24, CDCA7, or HELLS. ZBTB24 is a transcriptional activator of CDCA7, and CDCA7 and HELLS compose a chromatin remodeling complex and are involved in the maintenance DNA methylation through an interaction with UHRF1 in a feed‐forward manner. Furthermore, our recent studies possibly provided the missing link between DNA hypomethylation and the formation of the abnormal chromosomes; it could occur via aberrant transcription from the hypomethylated regions, followed by pathological R‐loop formation. The homologous‐recombination dominant condition caused by a defect in nonhomologous end joining observed in several types of ICF syndrome could facilitate the formation of multiradial chromosomes. Here, the latest knowledge regarding maintenance DNA methylation and chromosome stability provided by those studies is reviewed.

Congenital Diseases of DNA Replication: Clinical Phenotypes and Molecular Mechanisms

Deoxyribonucleic acid (DNA) replication can be divided into three major steps: initiation, elongation and termination. Each time a human cell divides, these steps must be reiteratively carried out. Disruption of DNA replication can lead to genomic instability, with the accumulation of point mutations or larger chromosomal anomalies such as rearrangements. While cancer is the most common class of disease associated with genomic instability, several congenital diseases with dysfunctional DNA replication give rise to similar DNA alterations. In this review, we discuss all congenital diseases that arise from pathogenic variants in essential replication genes across the spectrum of aberrant replisome assembly, origin activation and DNA synthesis. For each of these conditions, we describe their clinical phenotypes as well as molecular studies aimed at determining the functional mechanisms of disease, including the assessment of genomic stability. By comparing and contrasting these diseases, we hope to illuminate how the disruption of DNA replication at distinct steps affects human health in a surprisingly cell-type-specific manner.

Mapping replication timing domains genome wide in single mammalian cells with single-cell DNA replication sequencing

Replication timing (RT) domains are stable units of chromosome structure that are regulated in the context of development and disease. Conventional genome-wide RT mapping methods require many S-phase cells for either the effective enrichment of replicating DNA through bromodeoxyuridine (BrdU) immunoprecipitation or the determination of copy-number differences during S-phase, which precludes their application to non-abundant cell types and single cells. Here, we provide a simple, cost-effective, and robust protocol for single-cell DNA replication sequencing (scRepli-seq). The scRepli-seq methodology relies on whole-genome amplification (WGA) of genomic DNA (gDNA) from single S-phase cells and next-generation sequencing (NGS)-based determination of copy-number differences that arise between replicated and unreplicated DNA. Haplotype-resolved scRepli-seq, which distinguishes pairs of homologous chromosomes within a single cell, is feasible by using single-nucleotide polymorphism (SNP)/indel information. We also provide computational pipelines for quality control, normalization, and binarization of the scRepli-seq data. The experimental portion of this protocol (before sequencing) takes 3 d.

Control of DNA replication timing in the 3D genome

The 3D organization of mammalian chromatin was described more than 30 years ago by visualizing sites of DNA synthesis at different times during the S phase of the cell cycle. These early cytogenetic studies revealed structurally stable chromosome domains organized into subnuclear compartments. Active-gene-rich domains in the nuclear interior replicate early, whereas more condensed chromatin domains that are largely at the nuclear and nucleolar periphery replicate later. During the past decade, this spatiotemporal DNA replication programme has been mapped along the genome and found to correlate with epigenetic marks, transcriptional activity and features of 3D genome architecture such as chromosome compartments and topologically associated domains. But the causal relationship between these features and DNA replication timing and the regulatory mechanisms involved have remained an enigma. The recent identification of cis-acting elements regulating the replication time and 3D architecture of individual replication domains and of long non-coding RNAs that coordinate whole chromosome replication provide insights into such mechanisms.