Genome-wide stability of the DNA replication program in single mammalian cells

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.

Quantifying DNA replication speeds in single cells by scEdU-seq

In a human cell, thousands of replication forks simultaneously coordinate duplication of the entire genome. The rate at which this process occurs might depend on the epigenetic state of the genome and vary between, or even within, cell types. To accurately measure DNA replication speeds, we developed single-cell 5-ethynyl-2'-deoxyuridine sequencing to detect nascent replicated DNA. We observed that the DNA replication speed is not constant but increases during S phase of the cell cycle. Using genetic and pharmacological perturbations we were able to alter this acceleration of replication and conclude that DNA damage inflicted by the process of transcription limits the speed of replication during early S phase. In late S phase, during which less-transcribed regions replicate, replication accelerates and approaches its maximum speed.

A rewiring of DNA replication mediated by MRE11 exonuclease underlies primed-to-naive cell de-differentiation

Mouse embryonic stem cells (mESCs) in the primed pluripotency state, which resembles the post-implantation epiblast, can be de-differentiated in culture to a naive state that resembles the pre-implantation inner cell mass. We report that primed-to-naive mESC transition entails a significant slowdown of DNA replication forks and the compensatory activation of dormant origins. Using isolation of proteins on nascent DNA coupled to mass spectrometry, we identify key changes in replisome composition that are responsible for these effects. Naive mESC forks are enriched in MRE11 nuclease and other DNA repair proteins. MRE11 is recruited to newly synthesized DNA in response to transcription-replication conflicts, and its inhibition or genetic downregulation in naive mESCs is sufficient to restore the fork rate of primed cells. Transcriptomic analyses indicate that MRE11 exonuclease activity is required for the complete primed-to-naive mESC transition, demonstrating a direct link between DNA replication dynamics and the mESC de-differentiation process.

scAbsolute: measuring single-cell ploidy and replication status

Cancer cells often exhibit DNA copy number aberrations and can vary widely in their ploidy. Correct estimation of the ploidy of single-cell genomes is paramount for downstream analysis. Based only on single-cell DNA sequencing information, scAbsolute achieves accurate and unbiased measurement of single-cell ploidy and replication status, including whole-genome duplications. We demonstrate scAbsolute's capabilities using experimental cell multiplets, a FUCCI cell cycle expression system, and a benchmark against state-of-the-art methods. scAbsolute provides a robust foundation for single-cell DNA sequencing analysis across different technologies and has the potential to enable improvements in a number of downstream analyses.

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.

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.

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 both bovine and mouse cleavage stage embryos progress through S-phase in a defined pattern. Late replicating regions are associated with the nuclear lamina from the first cell cycle after fertilization, and contain few active origins, and few but long genes. 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 segregation of soma and germ line. Our studies show that the formation of early and late replicating regions is among the first layers of epigenetic regulation established on the mammalian genome after fertilization.

Single-cell DNA replication dynamics in genomically unstable cancers

Dysregulated DNA replication is both a cause and a consequence of aneuploidy, yet the dynamics of DNA replication in aneuploid cell populations remains understudied. We developed a new method, PERT, for inferring cell-specific DNA replication states from single-cell whole genome sequencing, and investigated clone-specific DNA replication dynamics in >50,000 cells obtained from a collection of aneuploid and clonally heterogeneous cell lines, xenografts and primary cancer tissues. Clone replication timing (RT) profiles correlated with future copy number changes in serially passaged cell lines. Cell type was the strongest determinant of RT heterogeneity, while whole genome doubling and mutational process were associated with accumulation of late S-phase cells and weaker RT associations. Copy number changes affecting chromosome X had striking impact on RT, with loss of the inactive X allele shifting replication earlier, and loss of inactive Xq resulting in reactivation of Xp. Finally, analysis of time series xenografts illustrate how cell cycle distributions approximate clone proliferation, recapitulating expected relationships between proliferation and fitness in treatment-naive and chemotherapeutic contexts.

Spontaneous mutagenesis in human cells is controlled by REV1-Polymerase ζ and PRIMPOL

Translesion DNA synthesis (TLS) facilitates replication over damaged or difficult-to-replicate templates by employing specialized DNA polymerases. We investigate the effect on spontaneous mutagenesis of three main TLS control mechanisms: REV1 and PCNA ubiquitylation that recruit TLS polymerases and PRIMPOL that creates post-replicative gaps. Using whole-genome sequencing of cultured human RPE-1 cell clones, we find that REV1 and Polymerase ζ are wholly responsible for one component of base substitution mutagenesis that resembles homologous recombination deficiency, whereas the remaining component that approximates oxidative mutagenesis is reduced in PRIMPOL-/- cells. Small deletions in short repeats appear in REV1-/-PCNAK164R/K164R double mutants, revealing an alternative TLS mechanism. Also, 500-5,000 bp deletions appear in REV1-/- and REV3L-/- mutants, and chromosomal instability is detectable in REV1-/-PRIMPOL-/- cells. Our results indicate that TLS protects the genome from deletions and large rearrangements at the expense of being responsible for the majority of spontaneous base substitutions.

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 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.

Condensed but liquid-like domain organization of active chromatin regions in living human cells

In eukaryotes, higher-order chromatin organization is spatiotemporally regulated as domains, for various cellular functions. However, their physical nature in living cells remains unclear (e.g., condensed domains or extended fiber loops; liquid-like or solid-like). Using novel approaches combining genomics, single-nucleosome imaging, and computational modeling, we investigated the physical organization and behavior of early DNA replicated regions in human cells, which correspond to Hi-C contact domains with active chromatin marks. Motion correlation analysis of two neighbor nucleosomes shows that nucleosomes form physically condensed domains with ~150-nm diameters, even in active chromatin regions. The mean-square displacement analysis between two neighbor nucleosomes demonstrates that nucleosomes behave like a liquid in the condensed domain on the ~150 nm/~0.5 s spatiotemporal scale, which facilitates chromatin accessibility. Beyond the micrometers/minutes scale, chromatin seems solid-like, which may contribute to maintaining genome integrity. Our study reveals the viscoelastic principle of the chromatin polymer; chromatin is locally dynamic and reactive but globally stable.

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.

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.

Polycomb repressive complexes 1 and 2 are each essential for maintenance of X inactivation in extra-embryonic lineages

In female mammals, one of the two X chromosomes becomes inactivated during development by X-chromosome inactivation (XCI). Although Polycomb repressive complex (PRC) 1 and PRC2 have both been implicated in gene silencing, their exact roles in XCI during in vivo development have remained elusive. To this end, we have studied mouse embryos lacking either PRC1 or PRC2. Here we demonstrate that the loss of either PRC has a substantial impact on maintenance of gene silencing on the inactive X chromosome (Xi) in extra-embryonic tissues, with overlapping yet different genes affected, indicating potentially independent roles of the two complexes. Importantly, a lack of PRC1 does not affect PRC2/H3K27me3 accumulation and a lack of PRC2 does not impact PRC1/H2AK119ub1 accumulation on the Xi. Thus PRC1 and PRC2 contribute independently to the maintenance of XCI in early post-implantation extra-embryonic lineages, revealing that both Polycomb complexes can be directly involved and differently deployed in XCI.

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.

Replication stress generates distinctive landscapes of DNA copy number alterations and chromosome scale losses

BACKGROUND

A major driver of cancer chromosomal instability is replication stress, the slowing or stalling of DNA replication. How replication stress and genomic instability are connected is not known. Aphidicolin-induced replication stress induces breakages at common fragile sites, but the exact causes of fragility are debated, and acute genomic consequences of replication stress are not fully explored.

RESULTS

We characterize DNA copy number alterations (CNAs) in single, diploid non-transformed cells, caused by one cell cycle in the presence of either aphidicolin or hydroxyurea. Multiple types of CNAs are generated, associated with different genomic regions and features, and observed copy number landscapes are distinct between aphidicolin and hydroxyurea-induced replication stress. Coupling cell type-specific analysis of CNAs to gene expression and single-cell replication timing analyses pinpointed the causative large genes of the most recurrent chromosome-scale CNAs in aphidicolin. These are clustered on chromosome 7 in RPE1 epithelial cells but chromosome 1 in BJ fibroblasts. Chromosome arm level CNAs also generate acentric lagging chromatin and micronuclei containing these chromosomes.

CONCLUSIONS

Chromosomal instability driven by replication stress occurs via focal CNAs and chromosome arm scale changes, with the latter confined to a very small subset of chromosome regions, potentially heavily skewing cancer genome evolution. Different inducers of replication stress lead to distinctive CNA landscapes providing the opportunity to derive copy number signatures of specific replication stress mechanisms. Single-cell CNA analysis thus reveals the impact of replication stress on the genome, providing insights into the molecular mechanisms which fuel chromosomal instability in cancer.

SmcHD1 underlies the formation of H3K9me3 blocks on the inactive X chromosome in mice

Stable silencing of the inactive X chromosome (Xi) in female mammals is crucial for the development of embryos and their postnatal health. SmcHD1 is essential for stable silencing of the Xi, and its functional deficiency results in derepression of many X-inactivated genes. Although SmcHD1 has been suggested to play an important role in the formation of higher-order chromatin structure of the Xi, the underlying mechanism is largely unknown. Here, we explore the epigenetic state of the Xi in SmcHD1-deficient epiblast stem cells and mouse embryonic fibroblasts in comparison with their wild-type counterparts. The results suggest that SmcHD1 underlies the formation of H3K9me3-enriched blocks on the Xi, which, although the importance of H3K9me3 has been largely overlooked in mice, play a crucial role in the establishment of the stably silenced state. We propose that the H3K9me3 blocks formed on the Xi facilitate robust heterochromatin formation in combination with H3K27me3, and that the substantial loss of H3K9me3 caused by SmcHD1 deficiency leads to aberrant distribution of H3K27me3 on the Xi and derepression of X-inactivated genes.

Large-Scale Chromatin Rearrangements in Cancer

Simple Summary

Cancers have many genetic mutations such as nucleotide changes, deletions, amplifications, and chromosome gains or losses. Some of these genetic alterations directly contribute to the initiation and progression of tumors. In parallel to these genetic changes, cancer cells acquire modifications to their chromatin landscape, i.e., to the marks that are carried by DNA and the histone proteins it is associated with. These “epimutations” have consequences for gene expression and genome stability, and also contribute to tumoral initiation and progression. Some of these chromatin changes are very local, affecting just one or a few genes. In contrast, some chromatin alterations observed in cancer are more widespread and affect a large part of the genome. In this review, we present different types of large-scale chromatin rearrangements in cancer, explain how they may occur, and why they are relevant for cancer diagnosis and treatment.

Abstract

Epigenetic abnormalities are extremely widespread in cancer. Some of them are mere consequences of transformation, but some actively contribute to cancer initiation and progression; they provide powerful new biological markers, as well as new targets for therapies. In this review, we examine the recent literature and focus on one particular aspect of epigenome deregulation: large-scale chromatin changes, causing global changes of DNA methylation or histone modifications. After a brief overview of the one-dimension (1D) and three-dimension (3D) epigenome in healthy cells and of its homeostasis mechanisms, we use selected examples to describe how many different events (mutations, changes in metabolism, and infections) can cause profound changes to the epigenome and fuel cancer. We then present the consequences for therapies and briefly discuss the role of single-cell approaches for the future progress of the field.

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.

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.

Cell cycle dynamics and developmental dynamics of the 3D genome: toward linking the two timescales

In the mammalian cell nucleus, chromosomes are folded differently in interphase and mitosis. Interphase chromosomes are relatively decondensed and display at least two unique layers of higher-order organization: topologically associating domains (TADs) and cell-type-specific A/B compartments, which correlate well with early/late DNA replication timing (RT). In mitosis, these structures rapidly disappear but are gradually reconstructed during G1 phase, coincident with the establishment of the RT program. However, these structures also change dynamically during cell differentiation and reprogramming, and yet we are surprisingly ignorant about the relationship between their cell cycle dynamics and developmental dynamics. In this review, we summarize the recent findings on this topic, discuss how these two processes might be coordinated with each other and its potential significance.

Highly rigid H3.1/H3.2-H3K9me3 domains set a barrier for cell fate reprogramming in trophoblast stem cells

Here, Hada et al. comprehensively analyzed epigenomic features of mouse trophoblast stem cells (TSCs). They used genome-wide, high-throughput analyses to show that the TSC genome contains large-scale (>1-Mb) rigid heterochromatin architectures that have a high degree of histone H3.1/3.2–H3K9me3 accumulation, termed TSC-defined highly heterochromatinized domains (THDs), and are uniquely developed in placental lineage cells that serve to protect them from fate reprogramming to stably maintain placental function.

The placenta is a highly evolved, specialized organ in mammals. It differs from other organs in that it functions only for fetal maintenance during gestation. Therefore, there must be intrinsic mechanisms that guarantee its unique functions. To address this question, we comprehensively analyzed epigenomic features of mouse trophoblast stem cells (TSCs). Our genome-wide, high-throughput analyses revealed that the TSC genome contains large-scale (>1-Mb) rigid heterochromatin architectures with a high degree of histone H3.1/3.2–H3K9me3 accumulation, which we termed TSC-defined highly heterochromatinized domains (THDs). Importantly, depletion of THDs by knockdown of CAF1, an H3.1/3.2 chaperone, resulted in down-regulation of TSC markers, such as Cdx2 and Elf5, and up-regulation of the pluripotent marker Oct3/4, indicating that THDs maintain the trophoblastic nature of TSCs. Furthermore, our nuclear transfer technique revealed that THDs are highly resistant to genomic reprogramming. However, when H3K9me3 was removed, the TSC genome was fully reprogrammed, giving rise to the first TSC cloned offspring. Interestingly, THD-like domains are also present in mouse and human placental cells in vivo, but not in other cell types. Thus, THDs are genomic architectures uniquely developed in placental lineage cells, which serve to protect them from fate reprogramming to stably maintain placental function.

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.

Camptothecin-Induced Replication Stress Affects DNA Replication Profiling by E/L Repli-Seq

E/L Repli-seq is a powerful tool for detecting cell type-specific replication landscapes in mammalian cells, but its potential to monitor DNA replication under replication stress awaits better understanding. Here, we used E/L Repli-seq to examine the temporal order of DNA replication in human retinal pigment epithelium cells treated with the topoisomerase I inhibitor camptothecin. We found that the replication profiles by E/L Repli-seq exhibit characteristic patterns after replication-stress induction, including the loss of specific initiation zones within individual early replication timing domains. We also observed global disappearance of the replication timing domain structures in the profiles, which can be explained by checkpoint-dependent suppression of replication initiation. Thus, our results demonstrate the effectiveness of E/L Repli-seq at identifying cells with replication-stress-induced altered DNA replication programs.

DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

Promoters are key organizers of the duplication of vertebrate genomes

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

Meiotic recombination mirrors patterns of germline replication in mice and humans

Genetic recombination generates novel trait combinations, and understanding how recombination is distributed across the genome is key to modern genetics. The PRDM9 protein defines recombination hotspots; however, megabase-scale recombination patterning is independent of PRDM9. The single round of DNA replication, which precedes recombination in meiosis, may establish these patterns; therefore, we devised an approach to study meiotic replication that includes robust and sensitive mapping of replication origins. We find that meiotic DNA replication is distinct; reduced origin firing slows replication in meiosis, and a distinctive replication pattern in human males underlies the subtelomeric increase in recombination. We detected a robust correlation between replication and both contemporary and historical recombination and found that replication origin density coupled with chromosome size determines the recombination potential of individual chromosomes. Our findings and methods have implications for understanding the mechanisms underlying DNA replication, genetic recombination, and the landscape of mammalian germline variation.

Coexistence of neuronal intranuclear inclusion disease and amyotrophic lateral sclerosis: an autopsy case

Background

Neuronal intranuclear inclusion disease (NIID) is a rare neurodegenerative disease. Pathologically, it is characterized by eosinophilic hyaline intranuclear inclusions in the cells of the visceral organs as well as central, peripheral, and autonomic nervous system cells. Recently, a GGC repeat expansion in the NOTCH2NLC gene has been identified as the etiopathological agent of NIID. Interestingly, this GGC repeat expansion was also reported in some patients with a clinical diagnosis of amyotrophic lateral sclerosis (ALS). However, there are no autopsy-confirmed cases of concurrent NIID and ALS.

Case presentation

A 60-year-old Taiwanese woman reported a four-month history of progressive weakness beginning in the right foot that spread to all four extremities. She was diagnosed with ALS because she met the revised El Escorial diagnostic criteria for definite ALS with upper and lower motor neuron involvement in the cervical, thoracic, and lumbosacral regions. She died of respiratory failure at 22 months from ALS onset, at the age of 62 years. Brain magnetic resonance imaging (MRI) revealed lesions in the medial part of the cerebellar hemisphere, right beside the vermis (paravermal lesions). The subclinical neuropathy, indicated by a nerve conduction study (NCS), prompted a potential diagnosis of NIID. Antemortem skin biopsy and autopsy confirmed the coexistence of pathology consistent with both ALS and NIID. We observed neither eccentric distribution of p62-positive intranuclear inclusions in the areas with abundant large motor neurons nor cytopathological coexistence of ALS and NIID pathology in motor neurons. This finding suggested that ALS and NIID developed independently in this patient.

Conclusions

We describe a case of concurrent NIID and ALS discovered during an autopsy. Abnormal brain MRI findings, including paravermal lesions, could indicate the coexistence of NIID even in patients with ALS showing characteristic clinical phenotypes.

Heterochromatin Networks: Topology, Dynamics, and Function (a Working Hypothesis)

Open systems can only exist by self-organization as pulsing structures exchanging matter and energy with the outer world. This review is an attempt to reveal the organizational principles of the heterochromatin supra-intra-chromosomal network in terms of nonlinear thermodynamics. The accessibility of the linear information of the genetic code is regulated by constitutive heterochromatin (CHR) creating the positional information in a system of coordinates. These features include scale-free splitting-fusing of CHR with the boundary constraints of the nucleolus and nuclear envelope. The analysis of both the literature and our own data suggests a radial-concentric network as the main structural organization principle of CHR regulating transcriptional pulsing. The dynamic CHR network is likely created together with nucleolus-associated chromatin domains, while the alveoli of this network, including springy splicing speckles, are the pulsing transcription hubs. CHR contributes to this regulation due to the silencing position variegation effect, stickiness, and flexible rigidity determined by the positioning of nucleosomes. The whole system acts in concert with the elastic nuclear actomyosin network which also emerges by self-organization during the transcriptional pulsing process. We hypothesize that the the transcriptional pulsing, in turn, adjusts its frequency/amplitudes specified by topologically associating domains to the replication timing code that determines epigenetic differentiation memory.

Derepression of inflammation-related genes link to microglia activation and neural maturation defect in a mouse model of Kleefstra syndrome

Haploinsufficiency of EHMT1, which encodes histone H3 lysine 9 (H3K9) methyltransferase G9a-like protein (GLP), causes Kleefstra syndrome (KS), a complex disorder of developmental delay and intellectual disability. Here, we examined whether postnatal supply of GLP can reverse the neurological phenotypes seen in Ehmt1^Δ/+ mice as a KS model. Ubiquitous GLP supply from the juvenile stage ameliorated behavioral abnormalities in Ehmt1^Δ/+ mice. Postnatal neuron-specific GLP supply was not sufficient for the improvement of abnormal behaviors but still reversed the reduction of H3K9me2 and spine number in Ehmt1^Δ/+ mice. Interestingly, some inflammatory genes, including IL-1β (Il1b), were upregulated and activated microglial cells increased in the Ehmt1^Δ/+ brain, and such phenotypes were also reversed by neuron-specific postnatal GLP supply. Il1b inactivation canceled the microglial and spine number phenotypes in the Ehmt1^Δ/+ mice. Thus, H3K9me2 and some neurological phenotypes are reversible, but behavioral abnormalities are more difficult to improve depending on the timing of GLP supply.

•

Activated microglias increase in a Ehmt1^Δ/+ mouse model of Kleefstra syndrome

•

Diminished H3K9me2 in Ehmt1^Δ/+ mouse neurons is reversed by post-natal GLP supply

•

GLP supply from the juvenile stage can improve abnormal behaviors of Ehmt1^Δ/+ mice

•

Il1b KO cancelles the microglial and spine number phenotypes in the Ehmt1^Δ/+ mice

Low Replicative Stress Triggers Cell-Type Specific Inheritable Advanced Replication Timing

DNA replication timing (RT), reflecting the temporal order of origin activation, is known as a robust and conserved cell-type specific process. Upon low replication stress, the slowing of replication forks induces well-documented RT delays associated to genetic instability, but it can also generate RT advances that are still uncharacterized. In order to characterize these advanced initiation events, we monitored the whole genome RT from six independent human cell lines treated with low doses of aphidicolin. We report that RT advances are cell-type-specific and involve large heterochromatin domains. Importantly, we found that some major late to early RT advances can be inherited by the unstressed next-cellular generation, which is a unique process that correlates with enhanced chromatin accessibility, as well as modified replication origin landscape and gene expression in daughter cells. Collectively, this work highlights how low replication stress may impact cellular identity by RT advances events at a subset of chromosomal domains.

Replication Stress, Genomic Instability, and Replication Timing: A Complex Relationship

The replication-timing program constitutes a key element of the organization and coordination of numerous nuclear processes in eukaryotes. This program is established at a crucial moment in the cell cycle and occurs simultaneously with the organization of the genome, thus indicating the vital significance of this process. With recent technological achievements of high-throughput approaches, a very strong link has been confirmed between replication timing, transcriptional activity, the epigenetic and mutational landscape, and the 3D organization of the genome. There is also a clear relationship between replication stress, replication timing, and genomic instability, but the extent to which they are mutually linked to each other is unclear. Recent evidence has shown that replication timing is affected in cancer cells, although the cause and consequence of this effect remain unknown. However, in-depth studies remain to be performed to characterize the molecular mechanisms of replication-timing regulation and clearly identify different cis- and trans-acting factors. The results of these studies will potentially facilitate the discovery of new therapeutic pathways, particularly for personalized medicine, or new biomarkers. This review focuses on the complex relationship between replication timing, replication stress, and genomic instability.

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

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

Dynamics of transcription-mediated conversion from euchromatin to facultative heterochromatin at the Xist promoter by Tsix

The fine-scale dynamics from euchromatin (EC) to facultative heterochromatin (fHC) has remained largely unclear. Here, we focus on Xist and its silencing initiator Tsix as a paradigm of transcription-mediated conversion from EC to fHC. In mouse epiblast stem cells, induction of Tsix recapitulates the conversion at the Xist promoter. Investigating the dynamics reveals that the conversion proceeds in a stepwise manner. Initially, a transient opened chromatin structure is observed. In the second step, gene silencing is initiated and dependent on Tsix, which is reversible and accompanied by simultaneous changes in multiple histone modifications. At the last step, maintenance of silencing becomes independent of Tsix and irreversible, which correlates with occupation of the -1 position of the transcription start site by a nucleosome and initiation of DNA methylation introduction. This study highlights the hierarchy of multiple chromatin events upon stepwise gene silencing establishment.

Chromosomal coordination and differential structure of asynchronous replicating regions

Stochastic asynchronous replication timing (AS-RT) is a phenomenon in which the time of replication of each allele is different, and the identity of the early allele varies between cells. By taking advantage of stable clonal pre-B cell populations derived from C57BL6/Castaneous mice, we have mapped the genome-wide AS-RT loci, independently of genetic differences. These regions are characterized by differential chromatin accessibility, mono-allelic expression and include new gene families involved in specifying cell identity. By combining population level mapping with single cell FISH, our data reveal the existence of a novel regulatory program that coordinates a fixed relationship between AS-RT regions on any given chromosome, with some loci set to replicate in a parallel and others set in the anti-parallel orientation. Our results show that AS-RT is a highly regulated epigenetic mark established during early embryogenesis that may be used for facilitating the programming of mono-allelic choice throughout development.

Most regions of the mammalian genome replicate both alleles in a synchronous manner, but some loci have been found to replicate asynchronously and the time of replication of each allele is different. Here the authors, by employing clonal mouse cells from a hybrid strain chart replication timing over the entire genome, using polymorphisms to distinguish between the paternal and maternal alleles.

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.

Loss of histone H3.3 results in DNA replication defects and altered origin dynamics in C. elegans

Histone H3.3 is a replication-independent variant of histone H3 with important roles in development, differentiation, and fertility. Here, we show that loss of H3.3 results in replication defects in Caenorhabditis elegans embryos at elevated temperatures. To characterize these defects, we adapt methods to determine replication timing, map replication origins, and examine replication fork progression. Our analysis of the spatiotemporal regulation of DNA replication shows that despite the very rapid embryonic cell cycle, the genome is replicated from early and late firing origins and is partitioned into domains of early and late replication. We find that under temperature stress conditions, additional replication origins become activated. Moreover, loss of H3.3 results in altered replication fork progression around origins, which is particularly evident at stress-activated origins. These replication defects are accompanied by replication checkpoint activation, a delayed cell cycle, and increased lethality in checkpoint-compromised embryos. Our comprehensive analysis of DNA replication in C. elegans reveals the genomic location of replication origins and the dynamics of their firing, and uncovers a role of H3.3 in the regulation of replication origins under stress conditions.

Chromatin and Nuclear Architecture: Shaping DNA Replication in 3D

In eukaryotes, DNA replication progresses through a finely orchestrated temporal and spatial program. The 3D genome structure and nuclear architecture have recently emerged as fundamental determinants of the replication program. Factors with established roles in replication have been recognized as genome organization regulators. Exploiting paradigms from yeasts and mammals, we discuss how DNA replication is regulated in time and space through DNA-associated trans-acting factors, diffusible limiting replication initiation factors, higher-order chromatin folding, dynamic origin localization, and specific nuclear microenvironments. We present an integrated model for the regulation of DNA replication in 3D and highlight the importance of accurate spatio-temporal regulation of DNA replication in physiology and disease.

BAMscale: quantification of next-generation sequencing peaks and generation of scaled coverage tracks

Background

Next-generation sequencing allows genome-wide analysis of changes in chromatin states and gene expression. Data analysis of these increasingly used methods either requires multiple analysis steps, or extensive computational time. We sought to develop a tool for rapid quantification of sequencing peaks from diverse experimental sources and an efficient method to produce coverage tracks for accurate visualization that can be intuitively displayed and interpreted by experimentalists with minimal bioinformatics background. We demonstrate its strength and usability by integrating data from several types of sequencing approaches.

Results

We have developed BAMscale, a one-step tool that processes a wide set of sequencing datasets. To demonstrate the usefulness of BAMscale, we analyzed multiple sequencing datasets from chromatin immunoprecipitation sequencing data (ChIP-seq), chromatin state change data (assay for transposase-accessible chromatin using sequencing: ATAC-seq, DNA double-strand break mapping sequencing: END-seq), DNA replication data (Okazaki fragments sequencing: OK-seq, nascent-strand sequencing: NS-seq, single-cell replication timing sequencing: scRepli-seq) and RNA-seq data. The outputs consist of raw and normalized peak scores (multiple normalizations) in text format and scaled bigWig coverage tracks that are directly accessible to data visualization programs. BAMScale also includes a visualization module facilitating direct, on-demand quantitative peak comparisons that can be used by experimentalists. Our tool can effectively analyze large sequencing datasets (~ 100 Gb size) in minutes, outperforming currently available tools.

Conclusions

BAMscale accurately quantifies and normalizes identified peaks directly from BAM files, and creates coverage tracks for visualization in genome browsers. BAMScale can be implemented for a wide set of methods for calculating coverage tracks, including ChIP-seq and ATAC-seq, as well as methods that currently require specialized, separate tools for analyses, such as splice-aware RNA-seq, END-seq and OK-seq for which no dedicated software is available. BAMscale is freely available on github (https://github.com/ncbi/BAMscale).

Biological phase separation: cell biology meets biophysics

Progress in development of biophysical analytic approaches has recently crossed paths with macromolecule condensates in cells. These cell condensates, typically termed liquid-like droplets, are formed by liquid-liquid phase separation (LLPS). More and more cell biologists now recognize that many of the membrane-less organelles observed in cells are formed by LLPS caused by interactions between proteins and nucleic acids. However, the detailed biophysical processes within the cell that lead to these assemblies remain largely unexplored. In this review, we evaluate recent discoveries related to biological phase separation including stress granule formation, chromatin regulation, and processes in the origin and evolution of life. We also discuss the potential issues and technical advancements required to properly study biological phase separation.

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

BACKGROUND

DNA replication in mammalian cells occurs in a defined temporal order during S phase, known as the replication timing (RT) programme. Replication timing is developmentally regulated and correlated with chromatin conformation and local transcriptional potential. Here, we present RT profiles of unprecedented temporal resolution in two human embryonic stem cell lines, human colon carcinoma line HCT116, and mouse embryonic stem cells and their neural progenitor derivatives.

RESULTS

Fine temporal windows revealed a remarkable degree of cell-to-cell conservation in RT, particularly at the very beginning and ends of S phase, and identified 5 temporal patterns of replication in all cell types, consistent with varying degrees of initiation efficiency. Zones of replication initiation (IZs) were detected throughout S phase and interacted in 3D space preferentially with other IZs of similar firing time. Temporal transition regions were resolved into segments of uni-directional replication punctuated at specific sites by small, inefficient IZs. Sites of convergent replication were divided into sites of termination or large constant timing regions consisting of many synchronous IZs in tandem. Developmental transitions in RT occured mainly by activating or inactivating individual IZs or occasionally by altering IZ firing time, demonstrating that IZs, rather than individual origins, are the units of developmental regulation. Finally, haplotype phasing revealed numerous regions of allele-specific and allele-independent asynchronous replication. Allele-independent asynchronous replication was correlated with the presence of previously mapped common fragile sites.

CONCLUSIONS

Altogether, these data provide a detailed temporal choreography of DNA replication in mammalian cells.

Genomic methods for measuring DNA replication dynamics

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density, and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analog incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

Single-cell DNA replication profiling identifies spatiotemporal developmental dynamics of chromosome organization

In mammalian cells, chromosomes are partitioned into megabase-sized topologically associating domains (TADs). TADs can be in either A (active) or B (inactive) subnuclear compartments, which exhibit early and late replication timing (RT), respectively. Here, we show that A/B compartments change coordinately with RT changes genome wide during mouse embryonic stem cell (mESC) differentiation. While A to B compartment changes and early to late RT changes were temporally inseparable, B to A changes clearly preceded late to early RT changes and transcriptional activation. Compartments changed primarily by boundary shifting, altering the compartmentalization of TADs facing the A/B compartment interface, which was conserved during reprogramming and confirmed in individual cells by single-cell Repli-seq. Differentiating mESCs altered single-cell Repli-seq profiles gradually but uniformly, transiently resembling RT profiles of epiblast-derived stem cells (EpiSCs), suggesting that A/B compartments might also change gradually but uniformly toward a primed pluripotent state. These results provide insights into how megabase-scale chromosome organization changes in individual cells during differentiation.

Metazoan DNA replication origins

DNA replication starts with the opening of DNA at sites called DNA replication origins. From the single sequence-specific DNA replication origin of the small Escherichia coli genome, up to thousands of origins that are necessary to replicate the large human genome, strict sequence specificity has been lost. Nevertheless, genome-wide analyses performed in the recent years, using different mapping methods, demonstrated that there are precise locations along the metazoan genome from which replication initiates. These sites contain relaxed sequence consensus and epigenetic features. There is flexibility in the choice of origins to be used during a given cell cycle, probably imposed by evolution and developmental constraints. Here, we will briefly describe their main features.

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

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

Genome-wide analysis of DNA replication timing in single cells: Yes! We're all individuals

Recent studies have accomplished the extraordinary feat of measuring the exact status of DNA replication in individual cells. We outline how these studies have revealed surprising uniformity in how cells replicate their DNA, and consider the implications of this remarkable technological advance.