Global reorganization of replication domains during embryonic stem cell differentiation

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.

G-quadruplex binding protein Rif1, a key regulator of replication timing

DNA replication is spatially and temporally regulated during S phase to execute efficient and coordinated duplication of entire genome. Various epigenomic mechanisms operate to regulate the timing and locations of replication. Among them, Rif1 plays a major role to shape the 'replication domains' that dictate which segments of the genome are replicated when and where in the nuclei. Rif1 achieves this task by generating higher-order chromatin architecture near nuclear membrane and by recruiting a protein phosphatase. Rif1 is a G4 binding protein, and G4 binding activity of Rif1 is essential for replication timing regulation in fission yeast. In this article, we first summarize strategies by which cells regulate their replication timing and then describe how Rif1 and its interaction with G4 contribute to regulation of chromatin architecture and replication timing.

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.

A unified framework for inferring the multi-scale organization of chromatin domains from Hi-C

Chromosomes are giant chain molecules organized into an ensemble of three-dimensional structures characterized with its genomic state and the corresponding biological functions. Despite the strong cell-to-cell heterogeneity, the cell-type specific pattern demonstrated in high-throughput chromosome conformation capture (Hi-C) data hints at a valuable link between structure and function, which makes inference of chromatin domains (CDs) from the pattern of Hi-C a central problem in genome research. Here we present a unified method for analyzing Hi-C data to determine spatial organization of CDs over multiple genomic scales. By applying statistical physics-based clustering analysis to a polymer physics model of the chromosome, our method identifies the CDs that best represent the global pattern of correlation manifested in Hi-C. The multi-scale intra-chromosomal structures compared across different cell types uncover the principles underlying the multi-scale organization of chromatin chain: (i) Sub-TADs, TADs, and meta-TADs constitute a robust hierarchical structure. (ii) The assemblies of compartments and TAD-based domains are governed by different organizational principles. (iii) Sub-TADs are the common building blocks of chromosome architecture. Our physically principled interpretation and analysis of Hi-C not only offer an accurate and quantitative view of multi-scale chromatin organization but also help decipher its connections with genome function.

An array of square blocks and checkerboard patterns characteristic to Hi-C data reflects the multi-scale organization of the chromatin chain. Deciphering the structures of chromatin domains from Hi-C and associating them with genome function are open problems of great importance in genome research. However, most existing methods are specialized in finding domains at different scales, making it difficult to integrate the solutions. Here we develop a unified framework for modeling and inferring domain structures over multiple scales, based on a physical model of the chromosome that reflects its nature as a three-dimensional object. Our method efficiently explores the space of domain solutions at different genomic scales, and systematically infers the chromatin domains over multiple scales from Hi-C data by employing a single tuning parameter. Our principled interpretation of Hi-C not only offers a quantitative view of multi-scale chromatin organization but also helps understand its connections with genome function.

Engineering 3D genome organization

Cancers and developmental disorders are associated with alterations in the 3D genome architecture in space and time (the fourth dimension). Mammalian 3D genome organization is complex and dynamic and plays an essential role in regulating gene expression and cellular function. To study the causal relationship between genome function and its spatio-temporal organization in the nucleus, new technologies for engineering and manipulating the 3D organization of the genome have been developed. In particular, CRISPR-Cas technologies allow programmable manipulation at specific genomic loci, enabling unparalleled opportunities in this emerging field of 3D genome engineering. We review advances in mammalian 3D genome engineering with a focus on recent manipulative technologies using CRISPR-Cas and related technologies.

The Temporal Order of DNA Replication Shaped by Mammalian DNA Methyltransferases

Multiple epigenetic pathways underlie the temporal order of DNA replication (replication timing) in the contexts of development and disease. DNA methylation by DNA methyltransferases (Dnmts) and downstream chromatin reorganization and transcriptional changes are thought to impact DNA replication, yet this remains to be comprehensively tested. Using cell-based and genome-wide approaches to measure replication timing, we identified a number of genomic regions undergoing subtle but reproducible replication timing changes in various Dnmt-mutant mouse embryonic stem (ES) cell lines that included a cell line with a drug-inducible Dnmt3a2 expression system. Replication timing within pericentromeric heterochromatin (PH) was shown to be correlated with redistribution of H3K27me3 induced by DNA hypomethylation: Later replicating PH coincided with H3K27me3-enriched regions. In contrast, this relationship with H3K27me3 was not evident within chromosomal arm regions undergoing either early-to-late (EtoL) or late-to-early (LtoE) switching of replication timing upon loss of the Dnmts. Interestingly, Dnmt-sensitive transcriptional up- and downregulation frequently coincided with earlier and later shifts in replication timing of the chromosomal arm regions, respectively. Our study revealed the previously unrecognized complex and diverse effects of the Dnmts loss on the mammalian DNA replication landscape.

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.

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.

Sequential peripheral enrichment of H2A.Zac and H3K9me2 during trophoblast differentiation in human embryonic stem cells

During the transition from pluripotency to a lineage-committed state, chromatin undergoes large-scale changes in structure, involving covalent modification of histone tails, use of histone variants and gene position changes with respect to the nuclear periphery. Here, using high-resolution microscopy and quantitative image analysis, we surveyed a panel of histone modifications for changes in nuclear peripheral enrichment during differentiation of human embryonic stem cells to a trophoblast-like lineage. We found two dynamic modifications at the nuclear periphery, acetylation of histone H2A.Z (H2A.Zac), and dimethylation of histone H3 at lysine 9 (H3K9me2). We demonstrate successive peripheral enrichment of these markers, with H2A.Zac followed by H3K9me2, over the course of 4 days. We find that H3K9me2 increases concomitantly with, but independently of, expression of lamin A, since deletion of lamin A did not affect H3K9me2 enrichment. We further show that inhibition of histone deacetylases causes persistent and increased H2A.Z acetylation at the periphery, delayed H3K9me2 enrichment and failure to differentiate. Our results show a concerted change in the nature of peripheral chromatin occurs upon differentiation into the trophoblast state.

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.

Epigenetic regulation of replication origin assembly: A role for histone H1 and chromatin remodeling factors

During early embryonic development in several metazoans, accurate DNA replication is ensured by high number of replication origins. This guarantees rapid genome duplication coordinated with fast cell divisions. In Xenopus laevis embryos this program switches to one with a lower number of origins at a developmental stage known as mid-blastula transition (MBT) when cell cycle length increases and gene transcription starts. Consistent with this regulation, somatic nuclei replicate poorly when transferred to eggs, suggesting the existence of an epigenetic memory suppressing replication assembly origins at all available sites. Recently, it was shown that histone H1 imposes a non-permissive chromatin configuration preventing replication origin assembly on somatic nuclei. This somatic state can be erased by SSRP1, a subunit of the FACT complex. Here, we further develop the hypothesis that this novel form of epigenetic memory might impact on different areas of vertebrate biology going from nuclear reprogramming to cancer development.

The Function and Evolution of Motile DNA Replication Systems in Ciliates

DNA replication is a ubiquitous and conserved cellular process. However, regulation of DNA replication is only understood in a small fraction of organisms that poorly represent the diversity of genetic systems in nature. Here we used computational and experimental approaches to examine the function and evolution of one such system, the replication band (RB) in spirotrich ciliates, which is a localized, motile hub that traverses the macronucleus while replicating DNA. We show that the RB can take unique forms in different species, from polar bands to a "replication envelope," where replication initiates at the nuclear periphery before advancing inward. Furthermore, we identify genes involved in cellular transport, including calcium transporters and cytoskeletal regulators, that are associated with the RB and may be involved in its function and translocation. These findings highlight the evolution and diversity of DNA replication systems and provide insights into the regulation of nuclear organization and processes.

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.

Detection of base analogs incorporated during DNA replication by nanopore sequencing

DNA synthesis is a fundamental requirement for cell proliferation and DNA repair, but no single method can identify the location, direction and speed of replication forks with high resolution. Mammalian cells have the ability to incorporate thymidine analogs along with the natural A, T, G and C bases during DNA synthesis, which allows for labeling of replicating or repaired DNA. Here, we demonstrate the use of the Oxford Nanopore Technologies MinION to detect 11 different thymidine analogs including CldU, BrdU, IdU as well as EdU alone or coupled to Biotin and other bulky adducts in synthetic DNA templates. We also show that the large adduct Biotin can be distinguished from the smaller analog IdU, which opens the possibility of using analog combinations to identify the location and direction of DNA synthesis. Furthermore, we detect IdU label on single DNA molecules in the genome of mouse pluripotent stem cells and using CRISPR/Cas9-mediated enrichment, determine replication rates using newly synthesized DNA strands in human mitochondrial DNA. We conclude that this novel method, termed Replipore sequencing, has the potential for on target examination of DNA replication in a wide range of biological contexts.

Cell cycle dynamics of lamina-associated DNA

In mammalian interphase nuclei, more than one thousand large genomic regions are positioned at the nuclear lamina (NL). These lamina‐associated domains (LADs) are involved in gene regulation and may provide a backbone for the folding of interphase chromosomes. Little is known about the dynamics of LADs during interphase, in particular at the onset of G1 phase and during DNA replication. We developed an antibody‐based variant of the DamID technology (named pA‐DamID) that allows us to map and visualize genome–NL interactions with high temporal resolution. Application of pA‐DamID combined with synchronization and cell sorting experiments reveals that LAD–NL contacts are generally rapidly established early in G1 phase. However, LADs on the distal ~25 Mb of most chromosomes tend to contact the NL first and then gradually detach, while centromere‐proximal LADs accumulate gradually at the NL. Furthermore, our data indicate that S‐phase chromatin shows transiently increased lamin interactions. These findings highlight a dynamic choreography of LAD–NL contacts during interphase progression and illustrate the usefulness of pA‐DamID to study the dynamics of genome compartmentalization.

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.

3D genome organization contributes to genome instability at fragile sites

Common fragile sites (CFSs) are regions susceptible to replication stress and are hotspots for chromosomal instability in cancer. Several features were suggested to underlie CFS instability, however, these features are prevalent across the genome. Therefore, the molecular mechanisms underlying CFS instability remain unclear. Here, we explore the transcriptional profile and DNA replication timing (RT) under mild replication stress in the context of the 3D genome organization. The results reveal a fragility signature, comprised of a TAD boundary overlapping a highly transcribed large gene with APH-induced RT-delay. This signature enables precise mapping of core fragility regions in known CFSs and identification of novel fragile sites. CFS stability may be compromised by incomplete DNA replication and repair in TAD boundaries core fragility regions leading to genomic instability. The identified fragility signature will allow for a more comprehensive mapping of CFSs and pave the way for investigating mechanisms promoting genomic instability in cancer.

Common fragile sites are regions susceptible to replication stress and are prone to chromosomal instability. Here, the authors, by analyzing the contribution of 3D chromatin organization, identify and characterize a fragility signature and precisely map these fragility regions.

Efficient, quick and easy-to-use DNA replication timing analysis with START-R suite

DNA replication must be faithful and follow a well-defined spatiotemporal program closely linked to transcriptional activity, epigenomic marks, intranuclear structures, mutation rate and cell fate determination. Among the readouts of the spatiotemporal program of DNA replication, replication timing analyses require not only complex and time-consuming experimental procedures, but also skills in bioinformatics. We developed a dedicated Shiny interactive web application, the START-R (Simple Tool for the Analysis of the Replication Timing based on R) suite, which analyzes DNA replication timing in a given organism with high-throughput data. It reduces the time required for generating and analyzing simultaneously data from several samples. It automatically detects different types of timing regions and identifies significant differences between two experimental conditions in ∼15 min. In conclusion, START-R suite allows quick, efficient and easier analyses of DNA replication timing for all organisms. This novel approach can be used by every biologist. It is now simpler to use this method in order to understand, for example, whether 'a favorite gene or protein' has an impact on replication process or, indirectly, on genomic organization (as Hi-C experiments), by comparing the replication timing profiles between wild-type and mutant cell lines.

Superresolution imaging reveals spatiotemporal propagation of human replication foci mediated by CTCF-organized chromatin structures

Mammalian DNA replication is initiated at numerous replication origins, which are clustered into thousands of replication domains (RDs) across the genome. However, it remains unclear whether the replication origins within each RD are activated stochastically or preferentially near certain chromatin features. To understand how DNA replication in single human cells is regulated at the sub-RD level, we directly visualized and quantitatively characterized the spatiotemporal organization, morphology, and in situ epigenetic signatures of individual replication foci (RFi) across S-phase at superresolution using stochastic optical reconstruction microscopy. Importantly, we revealed a hierarchical radial pattern of RFi propagation dynamics that reverses directionality from early to late S-phase and is diminished upon caffeine treatment or CTCF knockdown. Together with simulation and bioinformatic analyses, our findings point to a "CTCF-organized REplication Propagation" (CoREP) model, which suggests a nonrandom selection mechanism for replication activation at the sub-RD level during early S-phase, mediated by CTCF-organized chromatin structures. Collectively, these findings offer critical insights into the key involvement of local epigenetic environment in coordinating DNA replication across the genome and have broad implications for our conceptualization of the role of multiscale chromatin architecture in regulating diverse cell nuclear dynamics in space and time.

Rif1 Functions in a Tissue-Specific Manner To Control Replication Timing Through Its PP1-Binding Motif

Replication initiation in eukaryotic cells occurs asynchronously throughout S phase, yielding early- and late-replicating regions of the genome, a process known as replication timing (RT). RT changes during development to ensure accurate genome duplication and maintain genome stability. To understand the relative contributions that cell lineage, cell cycle, and replication initiation regulators have on RT, we utilized the powerful developmental systems available in Drosophila melanogaster We generated and compared RT profiles from mitotic cells of different tissues and from mitotic and endocycling cells of the same tissue. Our results demonstrate that cell lineage has the largest effect on RT, whereas switching from a mitotic to an endoreplicative cell cycle has little to no effect on RT. Additionally, we demonstrate that the RT differences we observed in all cases are largely independent of transcriptional differences. We also employed a genetic approach in these same cell types to understand the relative contribution the eukaryotic RT control factor, Rif1, has on RT control. Our results demonstrate that Rif1 can function in a tissue-specific manner to control RT. Importantly, the Protein Phosphatase 1 (PP1) binding motif of Rif1 is essential for Rif1 to regulate RT. Together, our data support a model in which the RT program is primarily driven by cell lineage and is further refined by Rif1/PP1 to ultimately generate tissue-specific RT programs.

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.

Local rewiring of genome-nuclear lamina interactions by transcription

Transcriptionally inactive genes are often positioned at the nuclear lamina (NL), as part of large lamina‐associated domains (LADs). Activation of such genes is often accompanied by repositioning toward the nuclear interior. How this process works and how it impacts flanking chromosomal regions are poorly understood. We addressed these questions by systematic activation or inactivation of individual genes, followed by detailed genome‐wide analysis of NL interactions, replication timing, and transcription patterns. Gene activation inside LADs typically causes NL detachment of the entire transcription unit, but rarely more than 50–100 kb of flanking DNA, even when multiple neighboring genes are activated. The degree of detachment depends on the expression level and the length of the activated gene. Loss of NL interactions coincides with a switch from late to early replication timing, but the latter can involve longer stretches of DNA. Inactivation of active genes can lead to increased NL contacts. These extensive datasets are a resource for the analysis of LAD rewiring by transcription and reveal a remarkable flexibility of interphase chromosomes.

Genome anchoring to nuclear landmarks drives functional compartmentalization of the nuclear space

The spatial organization of the genome contributes to essential functions such as transcription and chromosome integrity maintenance. The principles governing nuclear compartmentalization have been the focus of considerable research over the last decade. In these studies, the genome–nuclear structure interactions emerged as a main driver of this particular 3D genome organization. In this review, we describe the interactions between the genome and four major landmarks of the nucleus: the nuclear lamina, the nuclear pores, the pericentromeric heterochromatin and the nucleolus. We present the recent studies that identify sequences bound to these different locations and address the tethering mechanisms. We give an overview of the relevance of this organization in development and disease. Finally, we discuss the dynamic aspects and self-organizing properties that allow this complex architecture to be inherited.

Ubiquitin Dynamics in Stem Cell Biology: Current Challenges and Perspectives

Ubiquitination plays a central role in the regulation of stem cell self-renewal, propagation, and differentiation. In this review, the functions of ubiquitin dynamics in a myriad of cellular processes, acting along side the pluripotency network, to regulate embryonic stem cell identity are highlighted. The implication of deubiquitinases (DUBs) and E3 Ubiquitin (Ub) ligases in cellular functions beyond protein degradation is reported, including key functions in the regulation of mRNA stability, protein translation, and intra-cellular trafficking; and how it affects cell metabolism, the micro-environment, and chromatin organization is discussed. Finally, unsolved issues in the field are emphasized and will need to be tackled in order to fully understand the contribution of ubiquitin dynamics to stem cell self-renewal and differentiation.

Recent advances in the spatial organization of the mammalian genome

The mammalian genome is complex and presents a dynamic structural organization that reflects function. Organization of the genome inside the mammalian nucleus impacts all nuclear processes including but not limited to transcription, replication and repair, and in many biological contexts such as early development, differentiation and physiological adaptations. However, there is limited understating of how 3D organization of the mammalian genome regulates different nuclear processes. Recent advances in microscopy and a myriad of genomics methods -- ropelled by next-generation sequencing -- have advanced our knowledge of genome organization to a great extent. In this review, we discuss nuclear compartments in general and recent advances in the understanding of how mammalian genome is organized in these compartments with an emphasis on dynamics at the nuclear periphery.

Methods for mapping 3D chromosome architecture

Determining how chromosomes are positioned and folded within the nucleus is critical to understanding the role of chromatin topology in gene regulation. Several methods are available for studying chromosome architecture, each with different strengths and limitations. Established imaging approaches and proximity ligation-based chromosome conformation capture (3C) techniques (such as DNA-FISH and Hi-C, respectively) have revealed the existence of chromosome territories, functional nuclear landmarks (such as splicing speckles and the nuclear lamina) and topologically associating domains. Improvements to these methods and the recent development of ligation-free approaches, including GAM, SPRITE and ChIA-Drop, are now helping to uncover new aspects of 3D genome topology that confirm the nucleus to be a complex, highly organized organelle.

Evolution of replication origins in vertebrate genomes: rapid turnover despite selective constraints

The replication program of vertebrate genomes is driven by the chromosomal distribution and timing of activation of tens of thousands of replication origins. Genome-wide studies have shown the association of origins with promoters and CpG islands, and their enrichment in G-quadruplex motifs (G4). However, the genetic determinants driving their activity remain poorly understood. To gain insight on the constraints operating on origins, we conducted the first evolutionary comparison of origins across vertebrates. We generated a genome-wide map of chicken origins (the first of a bird genome), and performed a comparison with human and mouse maps. The analysis of intra-species polymorphism revealed a strong depletion of genetic diversity at the core of replication initiation loci. This depletion is not linked to the presence of G4 motifs, promoters or CpG islands. In contrast, we show that origins experienced a rapid turnover during vertebrate evolution, since pairwise comparisons of origin maps revealed that <24% of them are conserved among vertebrates. This study unravels the existence of a novel determinant of origins, the precise functional role of which remains to be determined. Despite the importance of replication initiation for the fitness of organisms, the distribution of origins along vertebrate chromosomes is highly flexible.

Synchronized replication of genes encoding the same protein complex in fast-proliferating cells

DNA replication perturbs the dosage balance among genes; at mid-S phase, early-replicating genes have doubled their copies while late-replicating ones have not. Dosage imbalance among genes, especially within members of a protein complex, is toxic to cells. However, the molecular mechanisms that cells use to deal with such imbalance remain not fully understood. Here, we validate at the genomic scale that the dosage between early- and late-replicating genes is imbalanced in HeLa cells. We propose the synchronized replication hypothesis that genes sensitive to stoichiometric relationships will be replicated simultaneously to maintain stoichiometry. In support of this hypothesis, we observe that genes encoding the same protein complex have similar replication timing but mainly in fast-proliferating cells such as embryonic stem cells and cancer cells. We find that the synchronized replication observed in cancer cells, but not in slow-proliferating differentiated cells, is due to convergent evolution during tumorigenesis that restores synchronized replication timing within protein complexes. Taken together, our study reveals that the demand for dosage balance during S phase plays an important role in the optimization of the replication-timing program; this selection is relaxed during differentiation as the cell cycle prolongs and is restored during tumorigenesis as the cell cycle shortens.

Assessment of the Utility of Gene Positioning Biomarkers in the Stratification of Prostate Cancers

There is a pressing need for additional clinical biomarkers to predict the aggressiveness of individual cancers. Here, we examine the potential usefulness of spatial genome organization as a prognostic tool for prostate cancer. Using fluorescence in situ hybridization on formalin-fixed, paraffin embedded human prostate tissue specimens, we compared the nuclear positions of four genes between clinically relevant subgroups of prostate tissues. We find that directional repositioning of SP100 and TGFB3 gene loci stratifies prostate cancers of differing Gleason scores. A more peripheral position of SP100 and TGFB3 in the nucleus, compared to benign tissues, is associated with low Gleason score cancers, whereas more internal positioning correlates with higher Gleason scores. Conversely, LMNA is more internally positioned in many non-metastatic prostate cancers, while its position is indistinguishable from benign tissue in metastatic cancer. The false positive rates were relatively low, whereas, the false negative rates of single or combinations of genes were high, limiting the clinical utility of this assay in its current form. Nevertheless, our findings of subtype-specific gene positioning patterns in prostate cancer provides proof-of-concept for the potential usefulness of spatial gene positioning for prognostic applications, and encourage further exploration of spatial gene positioning patterns to identify novel clinically relevant molecular biomarkers, which may aid treatment decisions for cancer patients.

Replication timing networks reveal a link between transcription regulatory circuits and replication timing control

DNA replication occurs in a defined temporal order known as the replication timing (RT) program and is regulated during development, coordinated with 3D genome organization and transcriptional activity. However, transcription and RT are not sufficiently coordinated to predict each other, suggesting an indirect relationship. Here, we exploit genome-wide RT profiles from 15 human cell types and intermediate differentiation stages derived from human embryonic stem cells to construct different types of RT regulatory networks. First, we constructed networks based on the coordinated RT changes during cell fate commitment to create highly complex RT networks composed of thousands of interactions that form specific functional subnetwork communities. We also constructed directional regulatory networks based on the order of RT changes within cell lineages, and identified master regulators of differentiation pathways. Finally, we explored relationships between RT networks and transcriptional regulatory networks (TRNs) by combining them into more complex circuitries of composite and bipartite networks. Results identified novel trans interactions linking transcription factors that are core to the regulatory circuitry of each cell type to RT changes occurring in those cell types. These core transcription factors were found to bind cooperatively to sites in the affected replication domains, providing provocative evidence that they constitute biologically significant directional interactions. Our findings suggest a regulatory link between the establishment of cell-type-specific TRNs and RT control during lineage specification.

Germline DNA replication timing shapes mammalian genome composition

Mammalian DNA replication is a highly organized and regulated process. Large, Mb-sized regions are replicated at defined times along S-phase. Replication Timing (RT) is thought to play a role in shaping the mammalian genome by affecting mutation rates. Previous analyses relied on somatic RT profiles. However, only germline mutations are passed on to offspring and affect genomic composition. Therefore, germ cell RT information is necessary to evaluate the influences of RT on the mammalian genome. We adapted the RT mapping technique for limited amounts of cells, and measured RT from two stages in the mouse germline - primordial germ cells (PGCs) and spermatogonial stem cells (SSCs). RT in germline cells exhibited stronger correlations to both mutation rate and recombination hotspots density than those of RT in somatic tissues, emphasizing the importance of using correct tissues-of-origin for RT profiling. Germline RT maps exhibited stronger correlations to additional genetic features including GC-content, transposable elements (SINEs and LINEs), and gene density. GC content stratification and multiple regression analysis revealed independent contributions of RT to SINE, gene, mutation, and recombination hotspot densities. Together, our results establish a central role for RT in shaping multiple levels of mammalian genome composition.

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.

Paradoxical association of TET loss of function with genome-wide DNA hypomethylation

Cancer genomes are characterized by focal increases in DNA methylation, co-occurring with widespread hypomethylation. Here, we show that TET loss of function results in a similar genomic footprint. Both 5hmC in wild-type (WT) genomes and DNA hypermethylation in TET-deficient genomes are largely confined to the active euchromatic compartment, consistent with the known functions of TET proteins in DNA demethylation and the known distribution of 5hmC at transcribed genes and active enhancers. In contrast, an unexpected DNA hypomethylation noted in multiple TET-deficient genomes is primarily observed in the heterochromatin compartment. In a mouse model of T cell lymphoma driven by TET deficiency (Tet2/3 DKO T cells), genomic analysis of malignant T cells revealed DNA hypomethylation in the heterochromatic genomic compartment, as well as reactivation of repeat elements and enrichment for single-nucleotide alterations, primarily in heterochromatic regions of the genome. Moreover, hematopoietic stem/precursor cells (HSPCs) doubly deficient for Tet2 and Dnmt3a displayed greater losses of DNA methylation than HSPCs singly deficient for Tet2 or Dnmt3a alone, potentially explaining the unexpected synergy between DNMT3A and TET2 mutations in myeloid and lymphoid malignancies. Tet1-deficient cells showed decreased localization of DNMT3A in the heterochromatin compartment compared with WT cells, pointing to a functional interaction between TET and DNMT proteins and providing a potential explanation for the hypomethylation observed in TET-deficient genomes. Our data suggest that TET loss of function may at least partially underlie the characteristic pattern of global hypomethylation coupled to regional hypermethylation observed in diverse cancer genomes, and highlight the potential contribution of heterochromatin hypomethylation to oncogenesis.

RT States: systematic annotation of the human genome using cell type-specific replication timing programs

MOTIVATION

The replication timing (RT) program has been linked to many key biological processes including cell fate commitment, 3D chromatin organization and transcription regulation. Significant technology progress now allows to characterize the RT program in the entire human genome in a high-throughput and high-resolution fashion. These experiments suggest that RT changes dynamically during development in coordination with gene activity. Since RT is such a fundamental biological process, we believe that an effective quantitative profile of the local RT program from a diverse set of cell types in various developmental stages and lineages can provide crucial biological insights for a genomic locus.

RESULTS

In this study, we explored recurrent and spatially coherent combinatorial profiles from 42 RT programs collected from multiple lineages at diverse differentiation states. We found that a Hidden Markov Model with 15 hidden states provide a good model to describe these genome-wide RT profiling data. Each of the hidden state represents a unique combination of RT profiles across different cell types which we refer to as 'RT states'. To understand the biological properties of these RT states, we inspected their relationship with chromatin states, gene expression, functional annotation and 3D chromosomal organization. We found that the newly defined RT states possess interesting genome-wide functional properties that add complementary information to the existing annotation of the human genome.

AVAILABILITY AND IMPLEMENTATION

R scripts for inferring HMM models and Perl scripts for further analysis are available https://github.com/PouletAxel/script_HMM_Replication_timing.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome wide decrease of DNA replication eye density at the midblastula transition of Xenopus laevis

During the first rapid divisions of early development in many species, the DNA:cytoplasm ratio increases until the midblastula transition (MBT) when transcription resumes and cell cycles lengthen. S phase is very rapid in early embryos, about 20-30 times faster than in differentiated cells. Using a combination of DNA fiber studies and a Xenopus laevis embryonic in vitro replication system, we show that S phase slows down shortly after the MBT owing to a genome wide decrease of replication eye density. Increasing the dNTP pool did not accelerate S phase or increase replication eye density implying that dNTPs are not rate limiting for DNA replication at the Xenopus MBT. Increasing the ratio of DNA:cytoplasm in egg extracts faithfully recapitulates changes in the spatial replication program in embryos, supporting the hypothesis that titration of soluble limiting factors could explain the observed changes in the DNA replication program at the MBT in Xenopus laevis.

Rapid Irreversible Transcriptional Reprogramming in Human Stem Cells Accompanied by Discordance between Replication Timing and Chromatin Compartment

The temporal order of DNA replication is regulated during development and is highly correlated with gene expression, histone modifications and 3D genome architecture. We tracked changes in replication timing, gene expression, and chromatin conformation capture (Hi-C) A/B compartments over the first two cell cycles during differentiation of human embryonic stem cells to definitive endoderm. Remarkably, transcriptional programs were irreversibly reprogrammed within the first cell cycle and were largely but not universally coordinated with replication timing changes. Moreover, changes in A/B compartment and several histone modifications that normally correlate strongly with replication timing showed weak correlation during the early cell cycles of differentiation but showed increased alignment in later differentiation stages and in terminally differentiated cell lines. Thus, epigenetic cell fate transitions during early differentiation can occur despite dynamic and discordant changes in otherwise highly correlated genomic properties.

Two contrasting classes of nucleolus-associated domains in mouse fibroblast heterochromatin

In interphase eukaryotic cells, almost all heterochromatin is located adjacent to the nucleolus or to the nuclear lamina, thus defining nucleolus-associated domains (NADs) and lamina-associated domains (LADs), respectively. Here, we determined the first genome-scale map of murine NADs in mouse embryonic fibroblasts (MEFs) via deep sequencing of chromatin associated with purified nucleoli. We developed a Bioconductor package called NADfinder and demonstrated that it identifies NADs more accurately than other peak-calling tools, owing to its critical feature of chromosome-level local baseline correction. We detected two distinct classes of NADs. Type I NADs associate frequently with both the nucleolar periphery and the nuclear lamina, and generally display characteristics of constitutive heterochromatin, including late DNA replication, enrichment of H3K9me3, and little gene expression. In contrast, Type II NADs associate with nucleoli but do not overlap with LADs. Type II NADs tend to replicate earlier, display greater gene expression, and are more often enriched in H3K27me3 than Type I NADs. The nucleolar associations of both classes of NADs were confirmed via DNA-FISH, which also detected Type I but not Type II probes enriched at the nuclear lamina. Type II NADs are enriched in distinct gene classes, including factors important for differentiation and development. In keeping with this, we observed that a Type II NAD is developmentally regulated, and present in MEFs but not in undifferentiated embryonic stem (ES) cells.

Gene expression dysregulation domains are not a specific feature of Down syndrome

Down syndrome (DS), trisomy of human chromosome 21 (Hsa21), results in a broad range of phenotypes. A recent study reported that DS cells show genome-wide transcriptional changes in which up- or down-regulated genes are clustered in gene expression dysregulation domains (GEDDs). GEDDs were also reported in fibroblasts derived from a DS mouse model duplicated for some Hsa21-orthologous genes, indicating cross-species conservation of this phenomenon. Here we investigate GEDDs using the Dp1Tyb mouse model of DS, which is duplicated for the entire Hsa21-orthologous region of mouse chromosome 16. Our statistical analysis shows that GEDDs are present both in DS cells and in Dp1Tyb mouse fibroblasts and hippocampus. However, we find that GEDDs do not depend on the DS genotype but occur whenever gene expression changes. We conclude that GEDDs are not a specific feature of DS but instead result from the clustering of co-regulated genes, a function of mammalian genome organisation.

Gene expression dysregulation domains (GEDDs) have been reported in Down syndrome (DS) cells, where changes in gene expression are clustered. Here the authors find that, while GEDDs are present in DS cells and in the Dp1Tyb mouse model of DS, GEDDs do not depend on the DS genotype and occur whenever gene expression changes, suggesting they result from the clustering of co-regulated genes as a function of mammalian genome organisation.

Sir2 suppresses transcription-mediated displacement of Mcm2-7 replicative helicases at the ribosomal DNA repeats

Repetitive DNA sequences within eukaryotic heterochromatin are poorly transcribed and replicate late in S-phase. In Saccharomyces cerevisiae, the histone deacetylase Sir2 is required for both transcriptional silencing and late replication at the repetitive ribosomal DNA arrays (rDNA). Despite the widespread association between transcription and replication timing, it remains unclear how transcription might impinge on replication, or vice versa. Here we show that, when silencing of an RNA polymerase II (RNA Pol II)-transcribed non-coding RNA at the rDNA is disrupted by SIR2 deletion, RNA polymerase pushes and thereby relocalizes replicative Mcm2-7 helicases away from their loading sites to an adjacent region with low nucleosome occupancy, and this relocalization is associated with increased rDNA origin efficiency. Our results suggest a model in which two of the major defining features of heterochromatin, transcriptional silencing and late replication, are mechanistically linked through suppression of polymerase-mediated displacement of replication initiation complexes.

Eukaryotic genomes typically contain large regions of repetitive DNA, referred to as heterochromatin, that are both transcriptionally silent and late replicating. We provide a possible explanation for the association between transcriptional silencing and late replication. Budding yeast contains a histone deacetylase called SIR2 that was originally identified as a transcriptional repressor, but was later also found to ensure late replication of repetitive ribosomal DNA (rDNA) sequences. We show that the transcription that occurs in the absence of SIR2 directly displaces the helicase required for replication initiation at the rDNA. This work represents an important advance in understanding the interplay between transcription and replication at repetitive sequences by directly linking transcription with replication machinery in heterochromatin.

The Landscape of L1 Retrotransposons in the Human Genome Is Shaped by Pre-insertion Sequence Biases and Post-insertion Selection

L1 retrotransposons are transposable elements and major contributors of genetic variation in humans. Where L1 integrates into the genome can directly impact human evolution and disease. Here, we experimentally induced L1 retrotransposition in cells and mapped integration sites at nucleotide resolution. At local scales, L1 integration is mostly restricted by genome sequence biases and the specificity of the L1 machinery. At regional scales, L1 shows a broad capacity for integration into all chromatin states, in contrast to other known mobile genetic elements. However, integration is influenced by the replication timing of target regions, suggesting a link to host DNA replication. The distribution of new L1 integrations differs from those of preexisting L1 copies, which are significantly reshaped by natural selection. Our findings reveal that the L1 machinery has evolved to efficiently target all genomic regions and underline a predominant role for post-integrative processes on the distribution of endogenous L1 elements.

Next-Generation Sequencing Enables Spatiotemporal Resolution of Human Centromere Replication Timing

Centromeres serve a critical function in preserving genome integrity across sequential cell divisions, by mediating symmetric chromosome segregation. The repetitive, heterochromatic nature of centromeres is thought to be inhibitory to DNA replication, but has also led to their underrepresentation in human reference genome assemblies. Consequently, centromeres have been excluded from genomic replication timing analyses, leaving their time of replication unresolved. However, the most recent human reference genome, hg38, included models of centromere sequences. To establish the experimental requirements for achieving replication timing profiles for centromeres, we sequenced G₁- and S-phase cells from five human cell lines, and aligned the sequence reads to hg38. We were able to infer DNA replication timing profiles for the centromeres in each of the five cell lines, which showed that centromere replication occurs in mid-to-late S phase. Furthermore, we found that replication timing was more variable between cell lines in the centromere regions than expected, given the distribution of variation in replication timing genome-wide. These results suggest the potential of these, and future, sequence models to enable high-resolution studies of replication in centromeres and other heterochromatic regions.

Replication timing in Drosophila and its peculiarities in polytene chromosomes

Drosophila melanogaster is one of the popular model organisms in DNA replication studies. Since the 1960s, DNA replication of polytene chromosomes has been extensively studied by cytological methods. In the recent two decades, the progress in our understanding of DNA replication was associated with new techniques. Use of fluorescent dyes increased the resolution of cytological methods significantly. High-throughput methods allowed analysis of DNA replication on a genome scale, as well as its correlation with chromatin structure and gene activi ty. Precise mapping of the cytological structures of polytene chromosomes to the genome assembly allowed comparison of replication between polytene chromosomes and chromosomes of diploid cells. New features of replication characteristic for D. melanogaster were described for both diploid and polytene chromosomes. Comparison of genomic replication profiles revealed a significant similarity between Drosophila and other well-studi ed eukaryotic species, such as human. Early replication is often confined to intensely transcribed gene-dense regions characterized by multiple replication initiation sites. Features of DNA replication in Drosophila might be explained by a compact genome. The organization of replication in polytene chromosomes has much in common with the organization of replication in chromosomes in diploid cells. The most important feature of replication in polytene chromosomes is its low rate and the dependence of S-phase duration on many factors: external and internal, local and global. The speed of replication forks in D. melanogaster polytene chromosomes is affected by SUUR and Rif1 proteins. It is not known yet how universal the mechanisms associated with these factors are, but their study is very promising.

Insights into the Link between the Organization of DNA Replication and the Mutational Landscape

The generation of a complete and accurate copy of the genetic material during each cell cycle is integral to cell growth and proliferation. However, genetic diversity is essential for adaptation and evolution, and the process of DNA replication is a fundamental source of mutations. Genome alterations do not accumulate randomly, with variations in the types and frequencies of mutations that arise in different genomic regions. Intriguingly, recent studies revealed a striking link between the mutational landscape of a genome and the spatial and temporal organization of DNA replication, referred to as the replication program. In our review, we discuss how this program may contribute to shaping the profile and spectrum of genetic alterations, with implications for genome dynamics and organismal evolution in natural and pathological contexts.

DNA Replication Timing Enters the Single-Cell Era

In mammalian cells, DNA replication timing is controlled at the level of megabase (Mb)-sized chromosomal domains and correlates well with transcription, chromatin structure, and three-dimensional (3D) genome organization. Because of these properties, DNA replication timing is an excellent entry point to explore genome regulation at various levels and a variety of studies have been carried out over the years. However, DNA replication timing studies traditionally required at least tens of thousands of cells, and it was unclear whether the replication domains detected by cell population analyses were preserved at the single-cell level. Recently, single-cell DNA replication profiling methods became available, which revealed that the Mb-sized replication domains detected by cell population analyses were actually well preserved in individual cells. In this article, we provide a brief overview of our current knowledge on DNA replication timing regulation in mammals based on cell population studies, outline the findings from single-cell DNA replication profiling, and discuss future directions and challenges.

Inactivation of CDK12 Delineates a Distinct Immunogenic Class of Advanced Prostate Cancer

Using integrative genomic analysis of 360 metastatic castration-resistant prostate cancer (mCRPC) samples, we identified a novel subtype of prostate cancer typified by biallelic loss of CDK12 that is mutually exclusive with tumors driven by DNA repair deficiency, ETS fusions, and SPOP mutations. CDK12 loss is enriched in mCRPC relative to clinically localized disease and characterized by focal tandem duplications (FTDs) that lead to increased gene fusions and marked differential gene expression. FTDs associated with CDK12 loss result in highly recurrent gains at loci of genes involved in the cell cycle and DNA replication. CDK12 mutant cases are baseline diploid and do not exhibit DNA mutational signatures linked to defects in homologous recombination. CDK12 mutant cases are associated with elevated neoantigen burden ensuing from fusion-induced chimeric open reading frames and increased tumor T cell infiltration/clonal expansion. CDK12 inactivation thereby defines a distinct class of mCRPC that may benefit from immune checkpoint immunotherapy.