Dynamic changes in replication timing and gene expression during lineage specification of human pluripotent stem cells

Inferring replication timing and proliferation dynamics from single-cell DNA sequencing data

Dysregulated DNA replication is a cause and a consequence of aneuploidy in cancer, yet the interplay between copy number alterations (CNAs), replication timing (RT) and cell cycle dynamics remain understudied in aneuploid tumors. We developed a probabilistic method, PERT, for simultaneous inference of cell-specific replication and copy number states from single-cell whole genome sequencing (scWGS) data. We used PERT to investigate clone-specific RT and proliferation dynamics in  >50,000 cells obtained from aneuploid and clonally heterogeneous cell lines, xenografts and primary cancers. We observed bidirectional relationships between RT and CNAs, with CNAs affecting X-inactivation producing the largest RT shifts. Additionally, we found that clone-specific S-phase enrichment positively correlated with ground-truth proliferation rates in genomically stable but not unstable cells. Together, these results demonstrate robust computational identification of S-phase cells from scWGS data, and highlight the importance of RT and cell cycle properties in studying the genomic evolution of aneuploid tumors.

Sex-specific DNA-replication in the early mammalian embryo

The timing of DNA replication in mammals is crucial for minimizing errors and influenced by genome usage and chromatin states. Replication timing in the newly formed mammalian embryo remains poorly understood. Here, we have investigated replication timing in mouse zygotes and 2-cell embryos, revealing that zygotes lack a conventional replication timing program, which then emerges in 2-cell embryos. This program differs from embryonic stem cells and generally correlates with transcription and genome compartmentalization of both parental genomes. However, consistent and systematic differences existed between the replication timing of the two parental genomes, including considerably later replication of maternal pericentromeric regions compared to paternal counterparts. Moreover, maternal chromatin modified by Polycomb Repressive Complexes in the oocyte, undergoes early replication, despite belonging to the typically late-replicating B-compartment of the genome. This atypical and asynchronous replication of the two parental genomes may advance our understanding of replication stress in early human embryos and trigger strategies to reduce errors and aneuploidies.

Emergence of enhancers at late DNA replicating regions

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

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

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

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.

Three-dimensional nuclear organisation and the DNA replication timing program

In eukaryotic cells, genome duplication is temporally organised according to a program referred to as the replication-timing (RT) program. The RT of individual genomic domains strikingly parallels the three-dimensional architecture of their chromatin contacts and subnuclear distribution. However, it is unclear whether this correspondence is coincidental or whether it indicates a causal and regulatory relationship. In either case, the nature of the molecular mechanisms ensuring this spatio-temporal coordination is still unknown. Here, we review recent evidence that begins to uncover the existence of a shared molecular machinery at the core of the spatio-temporal co-regulation of DNA replication and genome architecture. Finally, we discuss the outstanding, key question of the biological role of their coordination.

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.

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.

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.

Both cell autonomous and non-autonomous processes modulate the association between replication timing and mutation rate

Cancer somatic mutations are the product of multiple mutational and repair processes, some of which are tightly associated with DNA replication. Mutation rates (MR) are known to be higher in late replication timing (RT) regions, but different processes can affect this association. Systematic analysis of the mutational landscape of 2787 tumors from 32 tumor types revealed that approximately one third of the tumor samples show weak association between replication timing and mutation rate. Further analyses revealed that those samples have unique mutational signatures and are enriched with mutations in genes involved in DNA replication, DNA repair and chromatin structure. Surprisingly, analysis of differentially expressed genes between weak and strong RT-MR association groups revealed that tumors with weak association are enriched with genes associated with cell-cell communication and the immune system, suggesting a non-autonomous response to DNA damage.

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.

The landscape of somatic mutations in lymphoblastoid cell lines

Somatic mutations have important biological ramifications while exerting substantial rate, type, and genomic location heterogeneity. Yet, their sporadic occurrence makes them difficult to study at scale and across individuals. Lymphoblastoid cell lines (LCLs), a model system for human population and functional genomics, harbor large numbers of somatic mutations and have been extensively genotyped. By comparing 1,662 LCLs, we report that the mutational landscape of the genome varies across individuals in terms of the number of mutations, their genomic locations, and their spectra; this variation may itself be modulated by somatic trans-acting mutations. Mutations attributed to the translesion DNA polymerase η follow two different modes of formation, with one mode accounting for the hypermutability of the inactive X chromosome. Nonetheless, the distribution of mutations along the inactive X chromosome appears to follow an epigenetic memory of the active form.

Cell-type specificity of the human mutation landscape with respect to DNA replication dynamics

The patterns of genomic mutations are associated with various genomic features, most notably late replication timing, yet it remains contested which mutation types and signatures relate to DNA replication dynamics and to what extent. Here, we perform high-resolution comparisons of mutational landscapes between lymphoblastoid cell lines, chronic lymphocytic leukemia tumors, and three colon adenocarcinoma cell lines, including two with mismatch repair deficiency. Using cell-type-matched replication timing profiles, we demonstrate that mutation rates exhibit heterogeneous replication timing associations among cell types. This cell-type heterogeneity extends to the underlying mutational pathways, as mutational signatures show inconsistent replication timing bias between cell types. Moreover, replicative strand asymmetries exhibit similar cell-type specificity, albeit with different relationships to replication timing than mutation rates. Overall, we reveal an underappreciated complexity and cell-type specificity of mutational pathways and their relationship to replication timing.

Mutational signatures association with replication timing in normal cells reveals similarities and differences with matched cancer tissues

Mutational signatures' association with replication timing (RT) has been studied in cancer samples, but the RT distribution of somatic mutations in non-cancerous cells was only minimally explored. Here, we performed comprehensive analyses of mutational signatures in 2.9 million somatic mutations across multiple non-cancerous tissues, stratified by early and late RT regions. We found that many mutational processes are active mainly or solely in early RT, such as SBS16 in hepatocytes and SBS88 in the colon, or in late RT, such as SBS4 in lung and hepatocytes, and SBS18 across many tissues. The two ubiquitous signatures, SBS1 and SBS5, showed late and early bias, respectively, across multiple tissues and in mutations representing germ cells. We also performed a direct comparison with cancer samples in 4 matched tissue-cancer types. Unexpectedly, while for most signatures the RT bias was consistent in normal tissue and in cancer, we found that SBS1's late RT bias is lost in cancer.

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.

The evolution of the human DNA replication timing program

DNA is replicated according to a defined spatiotemporal program that is linked to both gene regulation and genome stability. The evolutionary forces that have shaped replication timing programs in eukaryotic species are largely unknown. Here, we studied the molecular causes and consequences of replication timing evolution across 94 humans, 95 chimpanzees, and 23 rhesus macaques. Replication timing differences recapitulated the species' phylogenetic tree, suggesting continuous evolution of the DNA replication timing program in primates. Hundreds of genomic regions had significant replication timing variation between humans and chimpanzees, of which 66 showed advances in replication origin firing in humans, while 57 were delayed. Genes overlapping these regions displayed correlated changes in expression levels and chromatin structure. Many human-chimpanzee variants also exhibited interindividual replication timing variation, pointing to ongoing evolution of replication timing at these loci. Association of replication timing variation with genetic variation revealed that DNA sequence evolution can explain replication timing variation between species. Taken together, DNA replication timing shows substantial and ongoing evolution in the human lineage that is driven by sequence alterations and could impact regulatory evolution at specific genomic sites.

Strand Asymmetries Across Genomic Processes

Across biological systems, a number of genomic processes, including transcription, replication, DNA repair, and transcription factor binding, display intrinsic directionalities. These directionalities are reflected in the asymmetric distribution of nucleotides, motifs, genes, transposon integration sites, and other functional elements across the two complementary strands. Strand asymmetries, including GC skews and mutational biases, have shaped the nucleotide composition of diverse organisms. The investigation of strand asymmetries often serves as a method to understand underlying biological mechanisms, including protein binding preferences, transcription factor interactions, retrotransposition, DNA damage and repair preferences, transcription-replication collisions, and mutagenesis mechanisms. Research into this subject also enables the identification of functional genomic sites, such as replication origins and transcription start sites. Improvements in our ability to detect and quantify DNA strand asymmetries will provide insights into diverse functionalities of the genome, the contribution of different mutational mechanisms in germline and somatic mutagenesis, and our knowledge of genome instability and evolution, which all have significant clinical implications in human disease, including cancer. In this review, we describe key developments that have been made across the field of genomic strand asymmetries, as well as the discovery of associated mechanisms.

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.

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

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

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

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

Global early replication disrupts gene expression and chromatin conformation in a single cell cycle

BACKGROUND

The early embryonic divisions of many organisms, including fish, flies, and frogs, are characterized by a very rapid S-phase caused by high rates of replication initiation. In somatic cells, S-phase is much longer due to both a reduction in the total number of initiation events and the imposition of a temporal order of origin activation. The physiological importance of changes in the rate and timing of replication initiation in S-phase remains unclear.

RESULTS

Here we assess the importance of the temporal control of replication initiation using a conditional system in budding yeast to drive the early replication of the majority of origins in a single cell cycle. We show that global early replication disrupts the expression of over a quarter of all genes. By deleting individual origins, we show that delaying replication is sufficient to restore normal gene expression, directly implicating origin firing control in this regulation. Global early replication disrupts nucleosome positioning and transcription factor binding during S-phase, suggesting that the rate of S-phase is important to regulate the chromatin landscape.

CONCLUSIONS

Together, these data provide new insight into the role of the temporal control of origin firing during S-phase for coordinating replication, gene expression, and chromatin establishment as occurs in the early embryo.

Genome architecture plasticity underlies DNA replication timing dynamics in cell differentiation

During the S-phase of eukaryotic cell cycle, DNA is replicated in a dedicatedly regulated temporal order, with regions containing active and inactive genes replicated early and late, respectively. Recent advances in sequencing technology allow us to explore the connection between replication timing (RT), histone modifications, and three-dimensional (3D) chromatin structure in diverse cell types. To characterize the dynamics during cell differentiation, corresponding sequencing data for human embryonic stem cells and four differentiated cell types were collected. By comparing RT and its extent of conservation before and after germ layer specification, the human genome was partitioned into distinct categories. Each category is then subject to comparisons on genomic, epigenetic, and chromatin 3D structural features. As expected, while constitutive early and late replication regions showed active and inactive features, respectively, dynamic regions with switched RT showed intermediate features. Surprisingly, although early-to-late replication and late-to-early replication regions showed similar histone modification patterns in hESCs, their structural preferences were opposite. Specifically, in hESCs, early-to-late replication regions tended to appear in the B compartment and large topologically associated domains, while late-to-early replication regions showed the opposite. Our results uncover the coordinated regulation of RT and 3D genome structure that underlies the loss of pluripotency and lineage commitment and indicate the importance and potential roles of genome architecture in biological processes.

Comprehensive analysis of DNA replication timing across 184 cell lines suggests a role for MCM10 in replication timing regulation

Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. Cells from a single patient carrying MCM10 mutations demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in the mutated MCM10 cells were predominantly comprised of replication delays and initiation site gains and losses. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel candidate modulator of DNA replication timing.

Determination of human DNA replication origin position and efficiency reveals principles of initiation zone organisation

Replication of the human genome initiates within broad zones of ∼150 kb. The extent to which firing of individual DNA replication origins within initiation zones is spatially stochastic or localised at defined sites remains a matter of debate. A thorough characterisation of the dynamic activation of origins within initiation zones is hampered by the lack of a high-resolution map of both their position and efficiency. To address this shortcoming, we describe a modification of initiation site sequencing (ini-seq), based on density substitution. Newly replicated DNA is rendered 'heavy-light' (HL) by incorporation of BrdUTP while unreplicated DNA remains 'light-light' (LL). Replicated HL-DNA is separated from unreplicated LL-DNA by equilibrium density gradient centrifugation, then both fractions are subjected to massive parallel sequencing. This allows precise mapping of 23,905 replication origins simultaneously with an assignment of a replication initiation efficiency score to each. We show that origin firing within early initiation zones is not randomly distributed. Rather, origins are arranged hierarchically with a set of very highly efficient origins marking zone boundaries. We propose that these origins explain much of the early firing activity arising within initiation zones, helping to unify the concept of replication initiation zones with the identification of discrete replication origin sites.

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.

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.

Strand asymmetry influences mismatch resolution during a single-strand annealing

Background

Biases of DNA repair can shape the nucleotide landscape of genomes at evolutionary timescales. The molecular mechanisms of those biases are still poorly understood because it is difficult to isolate the contributions of DNA repair from those of DNA damage.

Results

Here, we develop a genome-wide assay whereby the same DNA lesion is repaired in different genomic contexts. We insert thousands of barcoded transposons carrying a reporter of DNA mismatch repair in the genome of mouse embryonic stem cells. Upon inducing a double-strand break between tandem repeats, a mismatch is generated if the break is repaired through single-strand annealing. The resolution of the mismatch showed a 60–80% bias in favor of the strand with the longest 3′ flap. The location of the lesion in the genome and the type of mismatch had little influence on the bias. Instead, we observe a complete reversal of the bias when the longest 3′ flap is moved to the opposite strand by changing the position of the double-strand break in the reporter.

Conclusions

These results suggest that the processing of the double-strand break has a major influence on the repair of mismatches during a single-strand annealing.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02665-3.

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.

Cancer Mutational Processes Vary in Their Association with Replication Timing and Chromatin Accessibility

Cancer somatic mutations are the product of multiple mutational and repair processes, both of which are tightly associated with DNA replication. Distinctive patterns of somatic mutation accumulation, termed mutational signatures, are indicative of processes sustained within tumors. However, the association of various mutational processes with replication timing (RT) remains an open question. In this study, we systematically analyzed the mutational landscape of 2,787 tumors from 32 tumor types separately for early and late replicating regions using sequence context normalization and chromatin data to account for sequence and chromatin accessibility differences. To account for sequence differences between various genomic regions, an artificial genome-based approach was developed to expand the signature analyses to doublet base substitutions and small insertions and deletions. The association of mutational processes and RT was signature specific: Some signatures were associated with early or late replication (such as SBS7b and SBS7a, respectively), and others had no association. Most associations existed even after normalizing for genome accessibility. A focused mutational signature identification approach was also developed that uses RT information to improve signature identification; this approach found that SBS16, which is biased toward early replication, is strongly associated with better survival rates in liver cancer. Overall, this novel and comprehensive approach provides a better understanding of the etiology of mutational signatures, which may lead to improved cancer prevention, diagnosis, and treatment. SIGNIFICANCE: Many mutational processes associate with early or late replication timing regions independently of chromatin accessibility, enabling development of a focused identification approach to improve mutational signature detection.

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.

Collective regulation of chromatin modifications predicts replication timing during cell cycle

Replication timing (RT) associates with genome architecture, while having a mixed relationship to histone marks. By profiling replication at high resolution and assessing broad histone marks across the cell cycle at the resolution of RT with and without genetic perturbation, we address the causal relationship between histone marks and RT. Four primary chromatin states, including an uncharacterized H3K36me2 state, emerge and define 97% of the mappable genome. RT and local replication patterns (e.g., initiation zones) quantitatively associate with chromatin states, histone mark dynamics, and spatial chromatin structure. Manipulation of broad histone marks and enhancer elements by overexpressing the histone H3 lysine 9/36 tri-demethylase KDM4A impacts RT across 11% of the genome. Broad histone modification changes were strong predictors of the observed RT alterations. Lastly, replication within H3K36me2-enriched neighborhoods is sensitive to KDM4A overexpression and is controlled at a megabase scale. These studies establish a role for collective chromatin mark regulation in modulating RT.

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.

Development of a three-layer consecutive gene delivery system for enhanced bone regeneration

This study developed a three-layer consecutive gene delivery system (T-CGDS) for timely gene delivery into human mesenchymal stem cells (hMSCs). The timing of transcription factor expression is important to effectively induce bone differentiation. Therefore, a three-layered nanocomposite was fabricated using differently sized gold nanoparticles to promote bone regeneration and osteogenic differentiation. The core layer comprised 80 nm gold nanoparticles coupled with ATF4 pDNA. Following coating with heparin-conjugated Pluronic F-127 (HP-F127), 50 nm gold nanoparticles coupled with SP7 pDNA were added to fabricate a bi-layer system. After further coating with HP-F127, 20 nm gold nanoparticles combined with RUNX2 pDNA were added. Consequently, a T-CGDS measuring 350-450 nm was fabricated. Genes were released for more than 8 days, while the size of the T-CGDS decreased over time. When the T-CGDS was applied to hMSCs, the gene in the outer layer (RUNX2) was expressed first, followed by those in the middle (SP7) and core (ATF4) layers. The T-CGDS effectively induced bone differentiation and regeneration in vitro and in vivo. Timely delivery of the ATF4 gene to stem cells via the T-CGDS can greatly assist osteogenic differentiation involved in bone regeneration.

Chromosomal Mcm2-7 distribution and the genome replication program in species from yeast to humans

The spatio-temporal program of genome replication across eukaryotes is thought to be driven both by the uneven loading of pre-replication complexes (pre-RCs) across the genome at the onset of S-phase, and by differences in the timing of activation of these complexes during S phase. To determine the degree to which distribution of pre-RC loading alone could account for chromosomal replication patterns, we mapped the binding sites of the Mcm2-7 helicase complex (MCM) in budding yeast, fission yeast, mouse and humans. We observed similar individual MCM double-hexamer (DH) footprints across the species, but notable differences in their distribution: Footprints in budding yeast were more sharply focused compared to the other three organisms, consistent with the relative sequence specificity of replication origins in S. cerevisiae. Nonetheless, with some clear exceptions, most notably the inactive X-chromosome, much of the fluctuation in replication timing along the chromosomes in all four organisms reflected uneven chromosomal distribution of pre-replication complexes.

Gene-rich regions of the genome tend to replicate earlier in S phase than do repetitive and other non-genic regions. This may be an evolutionary consequence of the fact that replication later in S phase is associated with higher frequencies of mutation and genome rearrangement. Replication timing along the chromosome is determined by 1) events prior to S-phase that specify the locations where DNA replication can be initiated, referred to as origin licensing; and 2) the timing of activation of these licensed origins during S-phase, referred to as origin firing. To determine the relative importance of these two mechanisms, here we identify both the binding sites and the abundance of a key component of the origin licensing machinery in budding yeast, fission yeast, mice, and humans, namely the replicative helicase complex. We discovered that, with a few notable exceptions, which include the inactive X chromosome in mammals, the program of replication timing can be largely explained simply on the basis of origin licensing. Our results support a model for replication timing that emphasizes stochastic firing of origins that have been licensed before S phase begins.

Fused in sarcoma regulates DNA replication timing and kinetics

Fused in sarcoma (FUS) encodes an RNA-binding protein with diverse roles in transcriptional activation and RNA splicing. While oncogenic fusions of FUS and transcription factor DNA-binding domains are associated with soft tissue sarcomas, dominant mutations in FUS can cause amyotrophic lateral sclerosis. FUS has also been implicated in genome maintenance. However, the underlying mechanisms of its actions in genome stability are unknown. Here, we applied gene editing, functional reconstitution, and integrated proteomics and transcriptomics to illuminate roles for FUS in DNA replication and repair. Consistent with a supportive role in DNA double-strand break repair, FUS-deficient cells exhibited subtle alterations in the recruitment and retention of double-strand break-associated factors, including 53BP1 and BRCA1. FUS-/- cells also exhibited reduced proliferative potential that correlated with reduced speed of replication fork progression, diminished loading of prereplication complexes, enhanced micronucleus formation, and attenuated expression and splicing of S-phase-associated genes. Finally, FUS-deficient cells exhibited genome-wide alterations in DNA replication timing that were reversed upon re-expression of FUS complementary DNA. We also showed that FUS-dependent replication domains were enriched in transcriptionally active chromatin and that FUS was required for the timely replication of transcriptionally active DNA. These findings suggest that alterations in DNA replication kinetics and programming contribute to genome instability and functional defects in FUS-deficient cells.

Unraveling Ewing sarcoma tumorigenesis originating from patient-derived Mesenchymal Stem Cells

Ewing sarcoma (EwS) is characterized by pathognomonic translocations, most frequently fusing EWSR1 with FLI1. An estimated 30% of EwS tumors also display genetic alterations in STAG2, TP53, or CDKN2A (SPC). Numerous attempts to develop relevant EwS models from primary human cells have been unsuccessful in faithfully recapitulating the phenotypic, transcriptomic and epigenetic features of EwS. In this study, by engineering the t(11;22)(q24;q12) translocation together with a combination of SPC mutations, we generated a wide collection of immortalized cells (EWIma cells) tolerating EWSR1-FLI1 expression from primary mesenchymal stem cells (MSC) derived from an EwS patient. Within this model, SPC alterations strongly favored EwS oncogenicity. Xenograft experiments with independent EWIma cells induced tumors and metastases in mice, which displayed bona fide features of EwS. EWIma cells presented balanced but also more complex translocation profiles mimicking chromoplexy, which is frequently observed in EwS and other cancers. Collectively, these results demonstrate that bone marrow-derived MSCs are a source of origin for EwS and also provide original experimental models to investigate Ewing sarcomagenesis.

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.

Mammalian DNA Replication Timing

Immediately following the discovery of the structure of DNA and the semi-conservative replication of the parental DNA sequence into two new DNA strands, it became apparent that DNA replication is organized in a temporal and spatial fashion during the S phase of the cell cycle, correlated with the large-scale organization of chromatin in the nucleus. After many decades of limited progress, technological advances in genomics, genome engineering, and imaging have finally positioned the field to tackle mechanisms underpinning the temporal and spatial regulation of DNA replication and the causal relationships between DNA replication and other features of large-scale chromosome structure and function. In this review, we discuss these major recent discoveries as well as expectations for the coming decade.

Human L1 Transposition Dynamics Unraveled with Functional Data Analysis

Long INterspersed Elements-1 (L1s) constitute >17% of the human genome and still actively transpose in it. Characterizing L1 transposition across the genome is critical for understanding genome evolution and somatic mutations. However, to date, L1 insertion and fixation patterns have not been studied comprehensively. To fill this gap, we investigated three genome-wide data sets of L1s that integrated at different evolutionary times: 17,037 de novo L1s (from an L1 insertion cell-line experiment conducted in-house), and 1,212 polymorphic and 1,205 human-specific L1s (from public databases). We characterized 49 genomic features-proxying chromatin accessibility, transcriptional activity, replication, recombination, etc.-in the ±50 kb flanks of these elements. These features were contrasted between the three L1 data sets and L1-free regions using state-of-the-art Functional Data Analysis statistical methods, which treat high-resolution data as mathematical functions. Our results indicate that de novo, polymorphic, and human-specific L1s are surrounded by different genomic features acting at specific locations and scales. This led to an integrative model of L1 transposition, according to which L1s preferentially integrate into open-chromatin regions enriched in non-B DNA motifs, whereas they are fixed in regions largely free of purifying selection-depleted of genes and noncoding most conserved elements. Intriguingly, our results suggest that L1 insertions modify local genomic landscape by extending CpG methylation and increasing mononucleotide microsatellite density. Altogether, our findings substantially facilitate understanding of L1 integration and fixation preferences, pave the way for uncovering their role in aging and cancer, and inform their use as mutagenesis tools in genetic studies.

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.

Formation of a multi-layered 3-dimensional structure of the heterochromatin compartment during early mammalian development

During embryogenesis in mammals, the 3-dimensional (3D) genome organization changes globally in parallel with transcription changes in a cell-type specific manner. This involves the progressive formation of heterochromatin, the best example of which is the inactive X chromosome (Xi) in females, originally discovered as a compact 3D structure at the nuclear periphery known as the Barr body. The heterochromatin formation on the autosomes and the Xi is tightly associated with the differentiation state and the developmental potential of cells, making it an ideal readout of the cellular epigenetic state. At a glance, the heterochromatin appears to be uniform. However, recent studies are beginning to reveal a more complex picture, with multiple hierarchical levels co-existing within the heterochromatin compartment. Such hierarchical levels appear to exist in the heterochromatin compartment on autosomes as well as on the Xi. Here, we review recent progress in our understanding of the 3D genome organization changes during the period of differentiation surrounding pluripotency in vivo and in vitro, with a focus on the heterochromatin compartment. We first look at the whole genome, then focus on the Xi, and discuss their differences and similarities. Finally, we present a unified view of how the heterochromatin compartment is formed and regulated during early development. In particular, we emphasize that there are multiple layers within the heterochromatic compartment on both the autosomes and the Xi, with regulatory mechanisms common and specific to each layer.

Congenital Diseases of DNA Replication: Clinical Phenotypes and Molecular Mechanisms

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

Developmental differences in genome replication program and origin activation

To ensure error-free duplication of all (epi)genetic information once per cell cycle, DNA replication follows a cell type and developmental stage specific spatio-temporal program. Here, we analyze the spatio-temporal DNA replication progression in (un)differentiated mouse embryonic stem (mES) cells. Whereas telomeres replicate throughout S-phase, we observe mid S-phase replication of (peri)centromeric heterochromatin in mES cells, which switches to late S-phase replication upon differentiation. This replication timing reversal correlates with and depends on an increase in condensation and a decrease in acetylation of chromatin. We further find synchronous duplication of the Y chromosome, marking the end of S-phase, irrespectively of the pluripotency state. Using a combination of single-molecule and super-resolution microscopy, we measure molecular properties of the mES cell replicon, the number of replication foci active in parallel and their spatial clustering. We conclude that each replication nanofocus in mES cells corresponds to an individual replicon, with up to one quarter representing unidirectional forks. Furthermore, with molecular combing and genome-wide origin mapping analyses, we find that mES cells activate twice as many origins spaced at half the distance than somatic cells. Altogether, our results highlight fundamental developmental differences on progression of genome replication and origin activation in pluripotent cells.

Comparing DNA replication programs reveals large timing shifts at centromeres of endocycling cells in maize roots

Plant cells undergo two types of cell cycles–the mitotic cycle in which DNA replication is coupled to mitosis, and the endocycle in which DNA replication occurs in the absence of cell division. To investigate DNA replication programs in these two types of cell cycles, we pulse labeled intact root tips of maize (Zea mays) with 5-ethynyl-2’-deoxyuridine (EdU) and used flow sorting of nuclei to examine DNA replication timing (RT) during the transition from a mitotic cycle to an endocycle. Comparison of the sequence-based RT profiles showed that most regions of the maize genome replicate at the same time during S phase in mitotic and endocycling cells, despite the need to replicate twice as much DNA in the endocycle and the fact that endocycling is typically associated with cell differentiation. However, regions collectively corresponding to 2% of the genome displayed significant changes in timing between the two types of cell cycles. The majority of these regions are small with a median size of 135 kb, shift to a later RT in the endocycle, and are enriched for genes expressed in the root tip. We found larger regions that shifted RT in centromeres of seven of the ten maize chromosomes. These regions covered the majority of the previously defined functional centromere, which ranged between 1 and 2 Mb in size in the reference genome. They replicate mainly during mid S phase in mitotic cells but primarily in late S phase of the endocycle. In contrast, the immediately adjacent pericentromere sequences are primarily late replicating in both cell cycles. Analysis of CENH3 enrichment levels in 8C vs 2C nuclei suggested that there is only a partial replacement of CENH3 nucleosomes after endocycle replication is complete. The shift to later replication of centromeres and possible reduction in CENH3 enrichment after endocycle replication is consistent with a hypothesis that centromeres are inactivated when their function is no longer needed.

In traditional cell division, or mitosis, a cell’s genetic material is duplicated and then split between two daughter cells. In contrast, in some specialized cell types, the DNA is duplicated a second time without an intervening division step, resulting in cells that carry twice as much DNA. This phenomenon, which is called the endocycle, is common during plant development. At each step, DNA replication follows an ordered program in which highly compacted DNA is unraveled and replicated in sections at different times during the synthesis (S) phase. In plants, it is unclear whether traditional and endocycle programs are the same, especially since endocycling cells are typically in the process of differentiation. Using root tips of maize, we found that in comparison to replication in the mitotic cell cycle, there is a small portion of the genome whose replication in the endocycle is shifted in time, usually to later in S phase. Some of these regions are scattered around the genome and mostly coincide with active genes. However, the most prominent shifts occur in centromeres. The shift to later replication in centromeres is noteworthy because they orchestrate the process of separating duplicated chromosomes into daughter cells, a function that is not needed in the endocycle.

Programming pluripotent stem cells: Can't teach an old cell new DNA replication tricks

Pluripotent stem cells differentiate with varying efficiencies depending on the method of reprogramming that created them. In this issue, Paniza et al. (2020. J. Cell Biol.https://doi.org/10.1083/jcb.201909163) demonstrate that cells with lower differentiation potential retain some features of somatic DNA replication origin utilization and suffer more frequent DNA damage.

Replication timing alterations in leukemia affect clinically relevant chromosome domains

Human B-cell precursor acute lymphoid leukemias (BCP-ALLs) comprise a group of genetically and clinically distinct disease entities with features of differentiation arrest at known stages of normal B-lineage differentiation. We previously showed that BCP-ALL cells display unique and clonally heritable, stable DNA replication timing (RT) programs (ie, programs describing the variable order of replication and subnuclear 3D architecture of megabase-scale chromosomal units of DNA in different cell types). To determine the extent to which BCP-ALL RT programs mirror or deviate from specific stages of normal human B-cell differentiation, we transplanted immunodeficient mice with quiescent normal human CD34+ cord blood cells and obtained RT signatures of the regenerating B-lineage populations. We then compared these with RT signatures for leukemic cells from a large cohort of BCP-ALL patients with varied genetic subtypes and outcomes. The results identify BCP-ALL subtype-specific features that resemble specific stages of B-cell differentiation and features that seem to be associated with relapse. These results suggest that the genesis of BCP-ALL involves alterations in RT that reflect biologically significant and potentially clinically relevant leukemia-specific epigenetic changes.