Global reorganization of replication domains during embryonic stem cell differentiation

Chromatin-centric insights into DNA replication

DNA replication ensures the precise transmission of genetic information from parent to daughter cells. In eukaryotes, this process involves the replication of every base pair within a highly complex chromatin environment, encompassing multiple levels of chromatin structure and various chromatin metabolic processes. Recent evidence has demonstrated that DNA replication is strictly regulated in both temporal and spatial dimensions by factors such as 3D genome structure and transcription, which is crucial for maintaining genomic stability in each cell cycle. In this review, we discuss the diverse mechanisms that govern eukaryotic DNA replication, emphasizing the roles of chromatin architecture and transcriptional activity within the mammalian chromatin landscape. These insights provide a foundation for future investigations in this field.

RIF1 controls replication timing in early mouse embryos independently of lamina-associated nuclear organization

Cells must duplicate their genome before they divide to ensure equal transmission of genetic information. The genome is replicated with a defined temporal order, replication timing (RT), which is cell-type specific and linked to 3D-genome organization. During mammalian development, RT is initially not well defined and becomes progressively consolidated from the 4-cell stage. However, the molecular regulators are unknown. Here, by combining loss-of-function analysis with genome-wide investigation of RT in mouse embryos, we identify Rap1 interacting factor 1 (RIF1) as a regulator of the progressive consolidation of RT. Embryos without RIF1 show DNA replication features of an early, more totipotent state. RIF1 regulates the progressive stratification of RT values and its depletion leads to global RT changes and a more heterogeneous RT program. Developmental RT changes are disentangled from changes in transcription and nuclear organization, specifically nuclear lamina association. Our work provides molecular understanding of replication and genome organization at the beginning of mammalian development.

Cell-cycle dependence of bursty gene expression: insights from fitting mechanistic models to single-cell RNA-seq data

Bursty gene expression is characterized by two intuitive parameters, burst frequency and burst size, the cell-cycle dependence of which has not been extensively profiled at the transcriptome level. In this study, we estimate the burst parameters per allele in the G1 and G2/M cell-cycle phases for thousands of mouse genes by fitting mechanistic models of gene expression to messenger RNA count data, obtained by sequencing of single cells whose cell-cycle position has been inferred using a deep-learning method. We find that upon DNA replication, the median burst frequency approximately halves, while the burst size remains mostly unchanged. Genome-wide distributions of the burst parameter ratios between the G2/M and G1 phases are broad, indicating substantial heterogeneity in transcriptional regulation. We also observe a significant negative correlation between the burst frequency and size ratios, suggesting that regulatory processes do not independently control the burst parameters. We show that to accurately estimate the burst parameter ratios, mechanistic models must explicitly account for gene copy number variation and extrinsic noise due to the coupling of transcription to cell age across the cell cycle, but corrections for technical noise due to imperfect capture of RNA molecules in sequencing experiments are less critical.

Roles for the 3D genome in the cell cycle, DNA replication, and double strand break repair

Large eukaryotic genomes are packaged into the restricted area of the nucleus to protect the genetic code and provide a dedicated environment to read, copy and repair DNA. The physical organisation of the genome into chromatin loops and self-interacting domains provides the basic structural units of genome architecture. These structural arrangements are complex, multi-layered, and highly dynamic and influence how different regions of the genome interact. The role of chromatin structures during transcription via enhancer-promoter interactions is well established. Less understood is how nuclear architecture influences the plethora of chromatin transactions during DNA replication and repair. In this review, we discuss how genome architecture is regulated during the cell cycle to influence the positioning of replication origins and the coordination of DNA double strand break repair. The role of genome architecture in these cellular processes highlights its critical involvement in preserving genome integrity and cancer prevention.

KMT2C/KMT2D-dependent H3K4me1 mediates changes in DNA replication timing and origin activity during a cell fate transition

Mammalian genomes replicate in a cell-type-specific order during the S phase, correlated to transcriptional activity, histone modifications, and chromatin structure. The causal relationships between these features and DNA replication timing (RT), especially during cell fate changes, are largely unknown. Using machine learning, we quantify 21 chromatin features predicting local RT and RT changes during differentiation in embryonic stem cells (ESCs). About one-third of the genome shows RT changes during differentiation. Chromatin features accurately predict both steady-state RT and RT changes. Histone H3 lysine 4 monomethylation (H3K4me1), catalyzed by KMT2C and KMT2D (KMT2C/D), emerges as a top predictor. Loss of KMT2C/D or their enzymatic activities impairs RT changes during differentiation. This correlates with local H3K4me1 loss and reduced replication origin firing, while transcription remains largely unaffected. Our findings reveal KMT2C/D-dependent H3K4me1 as a key regulator of RT and replication initiation, a role that likely impacts diseases associated with KMT2C/D mutations.

Master transcription factor binding sites constitute the core of early replication control elements

Eukaryotic genomes replicate in a defined temporal order called the replication timing (RT) program. RT is developmentally regulated with potential to drive cell fate transitions, but mechanisms controlling RT remain elusive. We previously identified "Early Replication Control Elements" (ERCEs) necessary for early RT, domain-wide transcription, 3D chromatin architecture and compartmentalization in mouse embryonic stem cells (mESCs) but, deletions identifying ERCEs were large and encompassed many putative regulatory elements. Here, we show that ERCEs are compound elements whose RT activity can largely be accounted for by multiple sites of diverse master transcription factor binding (subERCEs), distinguished from other such sites by their long-range interactions. While deletion of subERCEs had large effects on both transcription and RT, deleting transcription start sites eliminated nearly all transcription with moderate effects on RT. Our results suggest a model in which subERCEs respond to diverse master transcription factors by functioning both as transcription enhancers and as elements that organize chromatin domains structurally and support early RT, potentially providing a feed-forward loop to drive robust epigenomic change during cell fate transitions.

R-loops acted on by RNase H1 influence DNA replication timing and genome stability in Leishmania

Genomes in eukaryotes normally undergo DNA replication in a choreographed temporal order, resulting in early and late replicating chromosome compartments. Leishmania, a human protozoan parasite, displays an unconventional DNA replication program in which the timing of DNA replication completion is chromosome size-dependent: larger chromosomes complete replication later then smaller ones. Here we show that both R-loops and RNase H1, a ribonuclease that resolves RNA-DNA hybrids, accumulate in Leishmania major chromosomes in a pattern that reflects their replication timing. Furthermore, we demonstrate that such differential organisation of R-loops, RNase H1 and DNA replication timing across the parasite's chromosomes correlates with size-dependent differences in chromatin accessibility, G quadruplex distribution and sequence content. Using conditional gene excision, we show that loss of RNase H1 leads to transient growth perturbation and permanently abrogates the differences in DNA replication timing across chromosomes, as well as altering levels of aneuploidy and increasing chromosome instability in a size-dependent manner. This work provides a link between R-loop homeostasis and DNA replication timing in a eukaryotic parasite and demonstrates that orchestration of DNA replication dictates levels of genome plasticity in Leishmania.

Transcription can be sufficient, but is not necessary, to advance replication timing

DNA replication timing (RT) is correlated with transcription during cell fate changes but there are many exceptions and our understanding of this relationship suffers from a paucity of reductionist approaches. Here, we manipulated length and strength of transcription in hybrid-genome mouse embryonic stem cells (mESCs) at a single locus upstream of the silent, late replicating, Pleiotrophin (Ptn) gene, directly comparing RT to nascent transcription rates at engineered vs. wild-type alleles. First, we inserted four reporter genes that differ only in their promoter. Two promoters transcribed the reporter gene at high rates and advanced RT. The other two transcribed at lower rates and did not advance RT. Since these promoters may prove useful in applications where effects on RT are undesirable, we confirmed the inability of one of them to advance RT at numerous ectopic sites. We next juxtaposed these same four promoters upstream of the Ptn transcription start site where they all transcribed the 96kb Ptn gene and advanced RT to different extents correlated with transcription rates. Indeed, a doxycycline-responsive promoter, which could not advance RT when induced as a small reporter gene, elicited a rapid and reversible RT advance proportional to the rate of transcription, providing direct evidence that transcription itself can advance RT. However, deletion of the Ptn promoter and enhancer, followed by directed differentiation to neural precursors, eliminated induction of transcription throughout the entire Ptn replication domain, without preventing the switch to early replication. Our results provide a solid empirical base with which to re-evaluate many decades of literature, demonstrating that length and strength of transcription is sufficient but not necessary to advance RT. Our results also provide a robust system in which to rapidly effect an RT change, permitting mechanistic studies of the role of transcription in RT and the consequences of RT changes to epigenomic remodeling.

Genome organization and stability in mammalian pre-implantation development

A largely stable genome is required for normal development, even as genetic change is an integral aspect of reproduction, genetic adaptation and evolution. Recent studies highlight a critical window of mammalian development with intrinsic DNA replication stress and genome instability in the first cell divisions after fertilization. Patterns of DNA replication and genome stability are established very early in mammals, alongside patterns of nuclear organization, and before the emergence of gene expression patterns, and prior to cell specification and germline formation. The study of DNA replication and genome stability in the mammalian embryo provides a unique cellular system due to the resetting of the epigenome to a totipotent state, and the de novo establishment of the patterns of nuclear organization, gene expression, DNA methylation, histone modifications and DNA replication. Studies on DNA replication and genome stability in the early mammalian embryo is relevant for understanding both normal and disease-causing genetic variation, and to uncover basic principles of genome regulation.

Structure and dynamics of nuclear A/B compartments and subcompartments

Mammalian chromosomes form a hierarchical structure within the cell nucleus, from chromatin loops, megabase (Mb)-sized topologically associating domains (TADs) to larger-scale A/B compartments. The molecular basis of the structures of loops and TADs has been actively studied. However, the A and B compartments, which correspond to early-replicating euchromatin and late-replicating heterochromatin, respectively, are still relatively unexplored. In this review, we focus on the A/B compartments, discuss their close relationship to DNA replication timing (RT), and introduce recent findings on the features of subcompartments revealed by detailed classification of the A/B compartments. In doing so, we speculate on the structure, potential function, and developmental dynamics of A/B compartments and subcompartments in mammalian cells.

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.

DNA replication in early mammalian embryos is patterned, predisposing lamina-associated regions to fragility

DNA replication in differentiated cells follows a defined program, but when and how it is established during mammalian development is not known. Here we show using single-cell sequencing, that late replicating regions are established in association with the B compartment and the nuclear lamina from the first cell cycle after fertilization on both maternal and paternal genomes. Late replicating regions contain a relative paucity of active origins and few but long genes and low G/C content. In both bovine and mouse embryos, replication timing patterns are established prior to embryonic genome activation. Chromosome breaks, which form spontaneously in bovine embryos at sites concordant with human embryos, preferentially locate to late replicating regions. In mice, late replicating regions show enhanced fragility due to a sparsity of dormant origins that can be activated under conditions of replication stress. This pattern predisposes regions with long neuronal genes to fragility and genetic change prior to separation of soma and germ cell lineages. Our studies show that the segregation of early and late replicating regions is among the first layers of genome organization established after fertilization.

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.

Developmental Changes in Genome Replication Progression in Pluripotent versus Differentiated Human Cells

DNA replication is a fundamental process ensuring the maintenance of the genome each time cells divide. This is particularly relevant early in development when cells divide profusely, later giving rise to entire organs. Here, we analyze and compare the genome replication progression in human embryonic stem cells, induced pluripotent stem cells, and differentiated cells. Using single-cell microscopic approaches, we map the spatio-temporal genome replication as a function of chromatin marks/compaction level. Furthermore, we mapped the replication timing of subchromosomal tandem repeat regions and interspersed repeat sequence elements. Albeit the majority of these genomic repeats did not change their replication timing from pluripotent to differentiated cells, we found developmental changes in the replication timing of rDNA repeats. Comparing single-cell super-resolution microscopic data with data from genome-wide sequencing approaches showed comparable numbers of replicons and large overlap in origins numbers and genomic location among developmental states with a generally higher origin variability in pluripotent cells. Using ratiometric analysis of incorporated nucleotides normalized per replisome in single cells, we uncovered differences in fork speed throughout the S phase in pluripotent cells but not in somatic cells. Altogether, our data define similarities and differences on the replication program and characteristics in human cells at different developmental states.

CDCA7-associated global aberrant DNA hypomethylation translates to localized, tissue-specific transcriptional responses

Disruption of cell division cycle associated 7 (CDCA7) has been linked to aberrant DNA hypomethylation, but the impact of DNA methylation loss on transcription has not been investigated. Here, we show that CDCA7 is critical for maintaining global DNA methylation levels across multiple tissues in vivo. A pathogenic Cdca7 missense variant leads to the formation of large, aberrantly hypomethylated domains overlapping with the B genomic compartment but without affecting the deposition of H3K9 trimethylation (H3K9me3). CDCA7-associated aberrant DNA hypomethylation translated to localized, tissue-specific transcriptional dysregulation that affected large gene clusters. In the brain, we identify CDCA7 as a transcriptional repressor and epigenetic regulator of clustered protocadherin isoform choice. Increased protocadherin isoform expression frequency is accompanied by DNA methylation loss, gain of H3K4 trimethylation (H3K4me3), and increased binding of the transcriptional regulator CCCTC-binding factor (CTCF). Overall, our in vivo work identifies a key role for CDCA7 in safeguarding tissue-specific expression of gene clusters via the DNA methylation pathway.

An emerging paradigm in epigenetic marking: coordination of transcription and replication

DNA replication and RNA transcription both utilize DNA as a template and therefore need to coordinate their activities. The predominant theory in the field is that in order for the replication fork to proceed, transcription machinery has to be evicted from DNA until replication is complete. If that does not occur, these machineries collide, and these collisions elicit various repair mechanisms which require displacement of one of the enzymes, often RNA polymerase, in order for replication to proceed. This model is also at the heart of the epigenetic bookmarking theory, which implies that displacement of RNA polymerase during replication requires gradual re-building of chromatin structure, which guides recruitment of transcriptional proteins and resumption of transcription. We discuss these theories but also bring to light newer data that suggest that these two processes may not be as detrimental to one another as previously thought. This includes findings suggesting that these processes can occur without fork collapse and that RNA polymerase may only be transiently displaced during DNA replication. We discuss potential mechanisms by which RNA polymerase may be retained at the replication fork and quickly rebind to DNA post-replication. These discoveries are important, not only as new evidence as to how these two processes are able to occur harmoniously but also because they have implications on how transcriptional programs are maintained through DNA replication. To this end, we also discuss the coordination of replication and transcription in light of revising the current epigenetic bookmarking theory of how the active gene status can be transmitted through S phase.

Incomplete reprogramming of DNA replication timing in induced pluripotent stem cells

Induced pluripotent stem cells (iPSCs) are the foundation of cell therapy. Differences in gene expression, DNA methylation, and chromatin conformation, which could affect differentiation capacity, have been identified between iPSCs and embryonic stem cells (ESCs). Less is known about whether DNA replication timing, a process linked to both genome regulation and genome stability, is efficiently reprogrammed to the embryonic state. To answer this, we compare genome-wide replication timing between ESCs, iPSCs, and cells reprogrammed by somatic cell nuclear transfer (NT-ESCs). While NT-ESCs replicate their DNA in a manner indistinguishable from ESCs, a subset of iPSCs exhibits delayed replication at heterochromatic regions containing genes downregulated in iPSCs with incompletely reprogrammed DNA methylation. DNA replication delays are not the result of gene expression or DNA methylation aberrations and persist after cells differentiate to neuronal precursors. Thus, DNA replication timing can be resistant to reprogramming and influence the quality of iPSCs.

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.

MLL3/MLL4 enzymatic activity shapes DNA replication timing

Mammalian genomes are replicated in a precise order during S phase, which is cell-type-specific1-3 and correlates with local transcriptional activity2,4-8, chromatin modifications9 and chromatin architecture1,10,11,12. However, the causal relationships between these features and the key regulators of DNA replication timing (RT) are largely unknown. Here, machine learning was applied to quantify chromatin features, including epigenetic marks, histone variants and chromatin architectural factors, best predicting local RT under steady-state and RT changes during early embryonic stem (ES) cell differentiation. About one-third of genome exhibited RT changes during the differentiation. Combined, chromatin features predicted steady-state RT and RT changes with high accuracy. Of these features, histone H3 lysine 4 monomethylation (H3K4me1) catalyzed by MLL3/4 (also known as KMT2C/D) emerged as a top predictor. Loss of Mll3/4 (but not Mll3 alone) or their enzymatic activity resulted in erasure of genome-wide RT dynamics during ES cell differentiation. Sites that normally gain H3K4me1 in a MLL3/4-dependent fashion during the transition failed to transition towards earlier RT, often with transcriptional activation unaffected. Further analysis revealed a requirement for MLL3/4 in promoting DNA replication initiation zones through MCM2 recruitment, providing a direct link for its role in regulating RT. Our results uncover MLL3/4-dependent H3K4me1 as a functional regulator of RT and highlight a causal relationship between the epigenome and RT that is largely uncoupled from transcription. These findings uncover a previously unknown role for MLL3/4-dependent chromatin functions which is likely relevant to the numerous diseases associated with MLL3/4 mutations.

The nucleolar protein GNL3 prevents resection of stalled replication forks

Faithful DNA replication requires specific proteins that protect replication forks and so prevent the formation of DNA lesions that may damage the genome. Identification of new proteins involved in this process is essential to understand how DNA lesions accumulate in cancer cells and how they tolerate them. Here, we show that human GNL3/nucleostemin, a GTP-binding protein localized mostly in the nucleolus and highly expressed in cancer cells, prevents nuclease-dependent resection of nascent DNA in response to replication stress. We demonstrate that inhibiting origin firing reduces resection. This suggests that the heightened replication origin activation observed upon GNL3 depletion largely drives the observed DNA resection probably due to the exhaustion of the available RPA pool. We show that GNL3 and DNA replication initiation factor ORC2 interact in the nucleolus and that the concentration of GNL3 in the nucleolus is required to limit DNA resection. We propose that the control of origin firing by GNL3 through the sequestration of ORC2 in the nucleolus is critical to prevent nascent DNA resection in response to replication stress.

Where and when to start: Regulating DNA replication origin activity in eukaryotic genomes

In eukaryotic genomes, hundreds to thousands of potential start sites of DNA replication named origins are dispersed across each of the linear chromosomes. During S-phase, only a subset of origins is selected in a stochastic manner to assemble bidirectional replication forks and initiate DNA synthesis. Despite substantial progress in our understanding of this complex process, a comprehensive 'identity code' that defines origins based on specific nucleotide sequences, DNA structural features, the local chromatin environment, or 3D genome architecture is still missing. In this article, we review the genetic and epigenetic features of replication origins in yeast and metazoan chromosomes and highlight recent insights into how this flexibility in origin usage contributes to nuclear organization, cell growth, differentiation, and genome stability.

Repli-seq Sample Preparation using Cell Sorting with Cell-Permeant Dyes

Replication timing is significantly correlated with gene expression and chromatin organization, changes dynamically during cell differentiation, and is altered in diseased states. Genome-wide analysis of replication timing is performed in actively replicating cells by Repli-seq. Current methods for Repli-seq require cells to be fixed in large numbers. This is a barrier for sample types that are sensitive to fixation or are in very limited numbers. In this article, we outline optimized methods to process live cells and intact nuclei for Repli-seq. Our protocol enables the processing of a smaller number of cells per sample and reduces processing time and sample loss while obtaining high-quality data. Further, for samples that tend to form clumps and are difficult to dissociate into a single-cell suspension, we also outline methods for isolation, staining, and processing of nuclei for Repli-seq. The Repli-seq data obtained from live cells and intact nuclei are comparable to those obtained from the standard protocols. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Live cell isolation and staining Alternate Protocol: Nuclei isolation and staining.

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.

Replication stress causes delayed mitotic entry and chromosome 12 fragility at the ANKS1B large neuronal gene in human induced pluripotent stem cells

Substantial background level of replication stress is a feature of embryonic and induced pluripotent stem cells (iPSCs), which can predispose to numerical and structural chromosomal instability, including recurrent aberrations of chromosome 12. In differentiated cells, replication stress-sensitive genomic regions, including common fragile sites, are widely mapped through mitotic chromosome break induction by mild aphidicolin treatment, an inhibitor of replicative polymerases. IPSCs exhibit lower apoptotic threshold and higher repair capacity hindering fragile site mapping. Caffeine potentiates genotoxic effects and abrogates G2/M checkpoint delay induced by chemical and physical mutagens. Using 5-ethynyl-2'-deoxyuridine (EdU) for replication labeling, we characterized the mitotic entry dynamics of asynchronous iPSCs exposed to aphidicolin and/or caffeine. Under the adjusted timing of replication stress exposure accounting revealed cell cycle delay, higher metaphase chromosome breakage rate was observed in iPSCs compared to primary lymphocytes. Using differential chromosome staining and subsequent locus-specific fluorescent in situ hybridization, we mapped the FRA12L fragile site spanning the large neuronal ANKS1B gene at 12q23.1, which may contribute to recurrent chromosome 12 missegregation and rearrangements in iPSCs. Publicly available data on the ANKS1B genetic alterations and their possible functional impact are reviewed. Our study provides the first evidence of common fragile site induction in iPSCs and reveals potential somatic instability of a clinically relevant gene during early human development and in vitro cell expansion.

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

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

Replication dynamics identifies the folding principles of the inactive X chromosome

Chromosome-wide late replication is an enigmatic hallmark of the inactive X chromosome (Xi). How it is established and what it represents remains obscure. By single-cell DNA replication sequencing, here we show that the entire Xi is reorganized to replicate rapidly and uniformly in late S-phase during X-chromosome inactivation (XCI), reflecting its relatively uniform structure revealed by 4C-seq. Despite this uniformity, only a subset of the Xi became earlier replicating in SmcHD1-mutant cells. In the mutant, these domains protruded out of the Xi core, contacted each other and became transcriptionally reactivated. 4C-seq suggested that they constituted the outermost layer of the Xi even before XCI and were rich in escape genes. We propose that this default positioning forms the basis for their inherent heterochromatin instability in cells lacking the Xi-binding protein SmcHD1 or exhibiting XCI escape. These observations underscore the importance of 3D genome organization for heterochromatin stability and gene regulation.

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.

Reorganization of the DNA replication landscape during adipogenesis is closely linked with adipogenic gene expression

The temporal order of DNA replication along the chromosomes is thought to reflect the transcriptional competence of the genome. During differentiation of mouse 3T3-L1 cells into adipocytes, cells undergo one or two rounds of cell division called mitotic clonal expansion (MCE). MCE is an essential step for adipogenesis; however, little is known about the regulation of DNA replication during this period. Here, we performed genome-wide mapping of replication timing (RT) in mouse 3T3-L1 cells before and during MCE, and identified a number of chromosomal regions shifting toward either earlier or later replication through two rounds of replication. These RT changes were confirmed in individual cells by single-cell DNA-replication sequencing. Coordinate changes between a shift toward earlier replication and transcriptional activation of adipogenesis-associated genes were observed. RT changes occurred before the full expression of these genes, indicating that RT reorganization might contribute to the mature adipocyte phenotype. To support this, cells undergoing two rounds of DNA replication during MCE had a higher potential to differentiate into lipid droplet-accumulating adipocytes, compared with cells undergoing a single round of DNA replication and non-replicating cells.

Quantifying propagation of DNA methylation and hydroxymethylation with iDEMS

DNA methylation is a critical epigenetic mark in mammalian cells. Many aspects of DNA methylation maintenance have been characterized; however, the exact kinetics of post-replicative methylation maintenance remain a subject of debate. Here we develop isolation of DNA by 5-ethynyl-deoxyuridine labelling for mass spectrometry (iDEMS), a highly sensitive, quantitative mass spectrometry-based method for measuring DNA modifications on metabolically labelled DNA. iDEMS reveals an unexpectedly hemi-methylated landscape on nascent DNA. Combining iDEMS with metabolic labelling reveals that methylation maintenance is outpaced by cell division in mouse embryonic stem cells. Our approach shows that hydroxymethylation is perpetually asymmetric between sister strands in favour of the parental, template strand. iDEMS can be coupled with immunoprecipitation of chromatin proteins, revealing features of DNA methylation-histone modification crosstalk and suggesting a model for interplay between methylation and nucleosome assembly. iDEMS therefore elucidates long-standing questions about DNA modification propagation and provides an important orthogonal technology to understanding this process in dynamic cellular contexts.

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.

The multiple de novo copy number variant (MdnCNV) phenomenon presents with peri-zygotic DNA mutational signatures and multilocus pathogenic variation

BACKGROUND

The multiple de novo copy number variant (MdnCNV) phenotype is described by having four or more constitutional de novo CNVs (dnCNVs) arising independently throughout the human genome within one generation. It is a rare peri-zygotic mutational event, previously reported to be seen once in every 12,000 individuals referred for genome-wide chromosomal microarray analysis due to congenital abnormalities. These rare families provide a unique opportunity to understand the genetic factors of peri-zygotic genome instability and the impact of dnCNV on human diseases.

METHODS

Chromosomal microarray analysis (CMA), array-based comparative genomic hybridization, short- and long-read genome sequencing (GS) were performed on the newly identified MdnCNV family to identify de novo mutations including dnCNVs, de novo single-nucleotide variants (dnSNVs), and indels. Short-read GS was performed on four previously published MdnCNV families for dnSNV analysis. Trio-based rare variant analysis was performed on the newly identified individual and four previously published MdnCNV families to identify potential genetic etiologies contributing to the peri-zygotic genomic instability. Lin semantic similarity scores informed quantitative human phenotype ontology analysis on three MdnCNV families to identify gene(s) driving or contributing to the clinical phenotype.

RESULTS

In the newly identified MdnCNV case, we revealed eight de novo tandem duplications, each ~ 1 Mb, with microhomology at 6/8 breakpoint junctions. Enrichment of de novo single-nucleotide variants (SNV; 6/79) and de novo indels (1/12) was found within 4 Mb of the dnCNV genomic regions. An elevated post-zygotic SNV mutation rate was observed in MdnCNV families. Maternal rare variant analyses identified three genes in distinct families that may contribute to the MdnCNV phenomenon. Phenotype analysis suggests that gene(s) within dnCNV regions contribute to the observed proband phenotype in 3/3 cases. CNVs in two cases, a contiguous gene duplication encompassing PMP22 and RAI1 and another duplication affecting NSD1 and SMARCC2, contribute to the clinically observed phenotypic manifestations.

CONCLUSIONS

Characteristic features of dnCNVs reported here are consistent with a microhomology-mediated break-induced replication (MMBIR)-driven mechanism during the peri-zygotic period. Maternal genetic variants in DNA repair genes potentially contribute to peri-zygotic genomic instability. Variable phenotypic features were observed across a cohort of three MdnCNV probands, and computational quantitative phenotyping revealed that two out of three had evidence for the contribution of more than one genetic locus to the proband's phenotype supporting the hypothesis of de novo multilocus pathogenic variation (MPV) in those families.

Dynamic chromosomal interactions and control of heterochromatin positioning by Ki-67

Ki-67 is a chromatin-associated protein with a dynamic distribution pattern throughout the cell cycle and is thought to be involved in chromatin organization. The lack of genomic interaction maps has hampered a detailed understanding of its roles, particularly during interphase. By pA-DamID mapping in human cell lines, we find that Ki-67 associates with large genomic domains that overlap mostly with late-replicating regions. Early in interphase, when Ki-67 is present in pre-nucleolar bodies, it interacts with these domains on all chromosomes. However, later in interphase, when Ki-67 is confined to nucleoli, it shows a striking shift toward small chromosomes. Nucleolar perturbations indicate that these cell cycle dynamics correspond to nucleolar maturation during interphase, and suggest that nucleolar sequestration of Ki-67 limits its interactions with larger chromosomes. Furthermore, we demonstrate that Ki-67 does not detectably control chromatin-chromatin interactions during interphase, but it competes with the nuclear lamina for interaction with late-replicating DNA, and it controls replication timing of (peri)centromeric regions. Together, these results reveal a highly dynamic choreography of genome interactions and roles for Ki-67 in heterochromatin organization.

Relationships between genome-wide R-loop distribution and classes of recurrent DNA breaks in neural stem/progenitor cells

Recent studies revealed classes of recurrent DNA double-strand breaks (DSBs) in neural stem/progenitor cells, including transcription-associated, promoter-proximal breaks and recurrent DSB clusters in late-replicating, long neural genes that may give rise to somatic brain mosaicism. The mechanistic factors promoting these different classes of DSBs in neural stem/progenitor cells are not understood. Here, we elucidated the genome-wide landscape of RNA:DNA hybrid structures called “R-loops” in primary neural stem/progenitor cells undergoing aphidicolin-induced, mild replication stress to assess the potential contribution of R-loops to the different, recurrent classes of DNA break “hotspots”. We find that R-loops in neural stem/progenitor cells undergoing mild replication stress are present primarily in early-replicating, transcribed regions and in genes with promoter GC skew that are associated with cell lineage-specific processes. Surprisingly, most long, neural genes that form recurrent DSB clusters do not show R-loop formation under conditions of mild replication stress. Our findings are consistent with a role of R-loop-associated processes in promoter-proximal DNA break formation in highly transcribed, early replicating regions but suggest that R-loops do not drive replication stress-induced, recurrent DSB cluster formation in most long, neural genes.

Impact of Chromosomal Context on Origin Selection and the Replication Program

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

Nuclear Compartments: An Incomplete Primer to Nuclear Compartments, Bodies, and Genome Organization Relative to Nuclear Architecture

This work reviews nuclear compartments, defined broadly to include distinct nuclear structures, bodies, and chromosome domains. It first summarizes original cytological observations before comparing concepts of nuclear compartments emerging from microscopy versus genomic approaches and then introducing new multiplexed imaging approaches that promise in the future to meld both approaches. I discuss how previous models of radial distribution of chromosomes or the binary division of the genome into A and B compartments are now being refined by the recognition of more complex nuclear compartmentalization. The poorly understood question of how these nuclear compartments are established and maintained is then discussed, including through the modern perspective of phase separation, before moving on to address possible functions of nuclear compartments, using the possible role of nuclear speckles in modulating gene expression as an example. Finally, the review concludes with a discussion of future questions for this field.

Intrinsic neural stem cell properties define brain hypersensitivity to genotoxic stress

Impaired replication has been previously linked to growth retardation and microcephaly; however, why the brain is critically affected compared with other organs remains elusive. Here, we report the differential response between early neural progenitors (neuroepithelial cells [NECs]) and fate-committed neural progenitors (NPs) to replication licensing defects. Our results show that, while NPs can tolerate altered expression of licensing factors, NECs undergo excessive replication stress, identified by impaired replication, increased DNA damage, and defective cell-cycle progression, leading eventually to NEC attrition and microcephaly. NECs that possess a short G1 phase license and activate more origins than NPs, by acquiring higher levels of DNA-bound MCMs. In vivo G1 shortening in NPs induces DNA damage upon impaired licensing, suggesting that G1 length correlates with replication stress hypersensitivity. Our findings propose that NECs possess distinct cell-cycle characteristics to ensure fast proliferation, although these inherent features render them susceptible to genotoxic stress.

Telomere-to-telomere human DNA replication timing profiles

The spatiotemporal organization of DNA replication produces a highly robust and reproducible replication timing profile. Sequencing-based methods for assaying replication timing genome-wide have become commonplace, but regions of high repeat content in the human genome have remained refractory to analysis. Here, we report the first nearly-gapless telomere-to-telomere replication timing profiles in human, using the T2T-CHM13 genome assembly and sequencing data for five cell lines. We find that replication timing can be successfully assayed in centromeres and large blocks of heterochromatin. Centromeric regions replicate in mid-to-late S-phase and contain replication-timing peaks at a similar density to other genomic regions, while distinct families of heterochromatic satellite DNA differ in their bias for replicating in late S-phase. The high degree of consistency in centromeric replication timing across chromosomes within each cell line prompts further investigation into the mechanisms dictating that some cell lines replicate their centromeres earlier than others, and what the consequences of this variation are.

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.

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.

HIRA-dependent boundaries between H3 variants shape early replication in mammals

The lack of a consensus DNA sequence defining replication origins in mammals has led researchers to consider chromatin as a means to specify these regions. However, to date, there is no mechanistic understanding of how this could be achieved and maintained given that nucleosome disruption occurs with each fork passage and with transcription. Here, by genome-wide mapping of the de novo deposition of the histone variants H3.1 and H3.3 in human cells during S phase, we identified how their dual deposition mode ensures a stable marking with H3.3 flanked on both sides by H3.1. These H3.1/H3.3 boundaries correspond to the initiation zones of early origins. Loss of the H3.3 chaperone HIRA leads to the concomitant disruption of H3.1/H3.3 boundaries and initiation zones. We propose that the HIRA-dependent deposition of H3.3 preserves H3.1/H3.3 boundaries by protecting them from H3.1 invasion linked to fork progression, contributing to a chromatin-based definition of early replication zones.

DNA replication fork speed underlies cell fate changes and promotes reprogramming

Totipotency emerges in early embryogenesis, but its molecular underpinnings remain poorly characterized. In the present study, we employed DNA fiber analysis to investigate how pluripotent stem cells are reprogrammed into totipotent-like 2-cell-like cells (2CLCs). We show that totipotent cells of the early mouse embryo have slow DNA replication fork speed and that 2CLCs recapitulate this feature, suggesting that fork speed underlies the transition to a totipotent-like state. 2CLCs emerge concomitant with DNA replication and display changes in replication timing (RT), particularly during the early S-phase. RT changes occur prior to 2CLC emergence, suggesting that RT may predispose to gene expression changes and consequent reprogramming of cell fate. Slowing down replication fork speed experimentally induces 2CLCs. In vivo, slowing fork speed improves the reprogramming efficiency of somatic cell nuclear transfer. Our data suggest that fork speed regulates cellular plasticity and that remodeling of replication features leads to changes in cell fate and reprogramming.

Totipotent cells in mouse embryos and 2-cell-like cells have slow DNA replication fork speed. Perturbations that slow replication fork speed promote 2-cell-like cell emergence and improve somatic cell nuclear transfer reprogramming and formation of induced pluripotent stem cell colonies.

Human spinal cord in vitro differentiation pace is initially maintained in heterologous embryonic environments

Species-specific differentiation pace in vitro indicates that some aspects of neural differentiation are governed by cell intrinsic properties. Here we describe a novel in vitro human neural-rosette assay that recapitulates dorsal spinal cord differentiation but proceeds more rapidly than in the human embryo, suggesting that it lacks endogenous signalling dynamics. To test whether in vitro conditions represent an intrinsic differentiation pace, human iPSC-derived neural rosettes were challenged by grafting into the faster differentiating chicken embryonic neural tube iso-chronically, or hetero-chronically into older embryos. In both contexts in vitro differentiation pace was initially unchanged, while long-term analysis revealed iso-chronic slowed and hetero-chronic conditions promoted human neural differentiation. Moreover, hetero-chronic conditions did not alter the human neural differentiation programme, which progressed to neurogenesis, while the host embryo advanced into gliogenesis. This study demonstrates that intrinsic properties limit human differentiation pace, and that timely extrinsic signals are required for progression through an intrinsic human neural differentiation programme.

Chemically Induced Chromosomal Interaction (CICI) method to study chromosome dynamics and its biological roles

Numerous intra- and inter-chromosomal contacts have been mapped in eukaryotic genomes, but it remains challenging to link these 3D structures to their regulatory functions. To establish the causal relationships between chromosome conformation and genome functions, we  develop a method, Chemically Induced Chromosomal Interaction (CICI), to selectively perturb the chromosome conformation at targeted loci. Using this method, long-distance chromosomal interactions can be induced dynamically between intra- or inter-chromosomal loci pairs, including the ones with very low Hi-C contact frequencies. Measurement of CICI formation time allows us to probe chromosome encounter dynamics between different loci pairs across the cell cycle. We also conduct two functional tests of CICI. We perturb the chromosome conformation near a DNA double-strand break and observe altered donor preference in homologous recombination; we force interactions between early and late-firing DNA replication origins and find no significant changes in replication timing. These results suggest that chromosome conformation plays a deterministic role in homology-directed DNA repair, but not in the establishment of replication timing. Overall, our study demonstrates that CICI is a powerful tool to study chromosome dynamics and 3D genome function.

Methods to selectively manipulate specific long-distance chromosomal interactions are limited. Here the authors develop a method called Chemically Induced Chromosomal Interaction (CICI) to engineer interactions and demonstrate that 3D conformation plays a causal role in establishing donor DNA preference during DNA repair.

Chromatin Velocity reveals epigenetic dynamics by single-cell profiling of heterochromatin and euchromatin

Recent efforts have succeeded in surveying open chromatin at the single-cell level, but high-throughput, single-cell assessment of heterochromatin and its underlying genomic determinants remains challenging. We engineered a hybrid transposase including the chromodomain (CD) of the heterochromatin protein-1α (HP-1α), which is involved in heterochromatin assembly and maintenance through its binding to trimethylation of the lysine 9 on histone 3 (H3K9me3), and developed a single-cell method, single-cell genome and epigenome by transposases sequencing (scGET-seq), that, unlike single-cell assay for transposase-accessible chromatin with sequencing (scATAC-seq), comprehensively probes both open and closed chromatin and concomitantly records the underlying genomic sequences. We tested scGET-seq in cancer-derived organoids and human-derived xenograft (PDX) models and identified genetic events and plasticity-driven mechanisms contributing to cancer drug resistance. Next, building upon the differential enrichment of closed and open chromatin, we devised a method, Chromatin Velocity, that identifies the trajectories of epigenetic modifications at the single-cell level. Chromatin Velocity uncovered paths of epigenetic reorganization during stem cell reprogramming and identified key transcription factors driving these developmental processes. scGET-seq reveals the dynamics of genomic and epigenetic landscapes underlying any cellular processes.

SAF-A promotes origin licensing and replication fork progression to ensure robust DNA replication

The organisation of chromatin is closely intertwined with biological activities of chromosome domains, including transcription and DNA replication status. Scaffold-attachment factor A (SAF-A), also known as heterogeneous nuclear ribonucleoprotein U (HNRNPU), contributes to the formation of open chromatin structure. Here, we demonstrate that SAF-A promotes the normal progression of DNA replication and enables resumption of replication after inhibition. We report that cells depleted of SAF-A show reduced origin licensing in G1 phase and, consequently, reduced origin activation frequency in S phase. Replication forks also progress less consistently in cells depleted of SAF-A, contributing to reduced DNA synthesis rate. Single-cell replication timing analysis revealed two distinct effects of SAF-A depletion: first, the boundaries between early- and late-replicating domains become more blurred; and second, SAF-A depletion causes replication timing changes that tend to bring regions of discordant domain compartmentalisation and replication timing into concordance. Associated with these defects, SAF-A-depleted cells show elevated formation of phosphorylated histone H2AX (γ-H2AX) and tend to enter quiescence. Overall, we find that SAF-A protein promotes robust DNA replication to ensure continuing cell proliferation.

Chromatin and Nuclear Dynamics in the Maintenance of Replication Fork Integrity

Replication of the eukaryotic genome is a highly regulated process and stringent control is required to maintain genome integrity. In this review, we will discuss the many aspects of the chromatin and nuclear environment that play key roles in the regulation of both unperturbed and stressed replication. Firstly, the higher order organisation of the genome into A and B compartments, topologically associated domains (TADs) and sub-nuclear compartments has major implications in the control of replication timing. In addition, the local chromatin environment defined by non-canonical histone variants, histone post-translational modifications (PTMs) and enrichment of factors such as heterochromatin protein 1 (HP1) plays multiple roles in normal S phase progression and during the repair of replicative damage. Lastly, we will cover how the spatial organisation of stalled replication forks facilitates the resolution of replication stress.

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.