SummaryEukaryotic chromosomes replicate in a temporal order known as the replication-timing program1. During mammalian development, at least half the genome changes replication timing, primarily in units of 400–800 kb (“replication domains”; RDs), whose positions are preserved in different cell types, conserved between species, and appear to confine long-range effects of chromosome rearrangements2–7. Early and late replication correlate strongly with open and closed chromatin compartments identified by high-resolution chromosome conformation capture (Hi-C), and, to a lesser extent, lamina-associated domains (LADs)4,5,8,9. Recent Hi-C mapping has unveiled a substructure of topologically-associating domains (TADs) that are largely conserved in their positions between cell types and are similar in size to RDs8,10. However, TADs can be further sub-stratified into smaller domains, challenging the significance of structures at any particular scale11,12. Moreover, attempts to reconcile TADs and LADs to replication-timing data have not revealed a common, underlying domain structure8,9,13. Here, we localize boundaries of RDs to the early-replicating border of replication-timing transitions and map their positions in 18 human and 13 mouse cell types. We demonstrate that, collectively, RD boundaries share a near one-to-one correlation with TAD boundaries, whereas within a cell type, adjacent TADs that replicate at similar times obscure RD boundaries, largely accounting for the previously reported lack of alignment. Moreover, cell-type specific replication timing of TADs partitions the genome into two large-scale sub-nuclear compartments revealing that replication-timing transitions are indistinguishable from late-replicating regions in chromatin composition and lamina association and accounting for the reduced correlation of replication timing to LADs and heterochromatin. Our results reconcile cell type specific sub-nuclear compartmentalization with developmentally stable chromosome domains and offer a unified model for large-scale chromosome structure and function.

The eukaryotic genome is replicated according to a specific spatio-temporal programme. However, little is known about both its molecular control and biological significance. Here, we identify mouse Rif1 as a key player in the regulation of DNA replication timing. We show that Rif1 deficiency in primary cells results in an unprecedented global alteration of the temporal order of replication. This effect takes place already in the first S-phase after Rif1 deletion and is neither accompanied by alterations in the transcriptional landscape nor by major changes in the biochemical identity of constitutive heterochromatin. In addition, Rif1 deficiency leads to both defective G1/S transition and chromatin re-organization after DNA replication. Together, these data offer a novel insight into the global regulation and biological significance of the replicationtiming programme in mammalian cells.

Structural variants (SVs) can contribute to oncogenesis through a variety of mechanisms. Despite their importance, the identification of SVs in cancer genomes remains challenging. Here, we present a framework that integrates optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole-genome sequencing to systematically detect SVs in a variety of normal or cancer samples and cell lines. We identify the unique strengths of each method and demonstrate that only integrative approaches can comprehensively identify SVs in the genome. By combining Hi-C and optical mapping, we resolve complex SVs and phase multiple SV events to a single haplotype. Furthermore, we observe widespread structural variation events affecting the functions of noncoding sequences, including the deletion of distal regulatory sequences, alteration of DNA replication timing, and the creation of novel three-dimensional chromatin structural domains. Our results indicate that noncoding SVs may be underappreciated mutational drivers in cancer genomes.

Summary DNA replication is temporally and spatially organized in all eukaryotes, yet the molecular control and biological function of the replication-timing program are poorly understood. A role for three-dimensional chromatin organization has been proposed. Rif1 is required for normal genome-wide regulation of replication timing, but its molecular function is poorly understood. Here we show that in mouse embryonic stem cells Rif1 coats late replicating domains and, together with Lamin B1 identifies the majority of the late replicating genome. Rif1 is an essential determinant of replication timing of non-Lamin B1-bound late domains. We further demonstrate that Rif1 defines and restricts the interactions between replication-timing domains during G1, thereby revealing a novel function of Rif1 as organizer of nuclear architecture. Loss of Rif1 affects both number and replication-timing specificity of the interactions between replication-timing domains. In addition, during S-phase Rif1 ensures temporal coordination of replication of interacting domains. In summary our study identifies Rif1 as the first molecular link between nuclear architecture organization and replication-timing establishment in mammals.

Graphical Abstract Highlights d Early replicating control elements (ERCEs) regulate replication timing d ERCEs regulate A/B compartmentalization and TAD architecture d ERCEs form CTCF-independent loops and have features of enhancer/promoters d ERCEs enable genetic dissection of large-scale chromosome structure and function SUMMARYThe temporal order of DNA replication (replication timing [RT]) is highly coupled with genome architecture, but cis-elements regulating either remain elusive. We created a series of CRISPR-mediated deletions and inversions of a pluripotency-associated topologically associating domain (TAD) in mouse ESCs. CTCF-associated domain boundaries were dispensable for RT. CTCF protein depletion weakened most TAD boundaries but had no effect on RT or A/B compartmentalization genome-wide. By contrast, deletion of three intra-TAD CTCF-independent 3D contact sites caused a domain-wide earlyto-late RT shift, an A-to-B compartment switch, weakening of TAD architecture, and loss of transcription. The dispensability of TAD boundaries and the necessity of these ''early replication control elements'' (ERCEs) was validated by deletions and inversions at additional domains. Our results demonstrate that discrete cis-regulatory elements orchestrate domain-wide RT, A/B compartmentalization, TAD architecture, and transcription, revealing fundamental principles linking genome structure and function.

Mammalian genomes are partitioned into domains that replicate in a defined temporal order. These domains can replicate at similar times in all cell types (constitutive) or at cell type-specific times (developmental). Genome-wide chromatin conformation capture (Hi-C) has revealed sub-megabase topologically associating domains (TADs), which are the structural counterparts of replication domains. Hi-C also segregates inter-TAD contacts into defined 3D spatial compartments that align precisely to genome-wide replication timing profiles. Determinants of the replication-timing program are re-established during early G1 phase of each cell cycle and lost in G2 phase, but it is not known when TAD structure and inter-TAD contacts are re-established after their elimination during mitosis. Here, we use multiplexed 4C-seq to study dynamic changes in chromatin organization during early G1. We find that both establishment of TADs and their compartmentalization occur during early G1, within the same time frame as establishment of the replication-timing program. Once established, this 3D organization is preserved either after withdrawal into quiescence or for the remainder of interphase including G2 phase, implying 3D structure is not sufficient to maintain replication timing. Finally, we find that developmental domains are less well compartmentalized than constitutive domains and display chromatin properties that distinguish them from early and late constitutive domains. Overall, this study uncovers a strong connection between chromatin re-organization during

Mammalian DNA replication is regulated via multi-replicon segments that replicate in a defined temporal order during S-phase. Further, early/late replication of RDs corresponds to active/inactive chromatin interaction compartments. Although replication origins are selected stochastically, variation in replication timing is poorly understood. Here we devise a strategy to measure variation in replication timing using DNA copy number in single mouse embryonic stem cells. We find that borders between replicated and unreplicated DNA are highly conserved between cells, demarcating active and inactive compartments of the nucleus. Fifty percent of replication events deviated from their average replication time by ± 15% of S phase. This degree of variation is similar between cells, between homologs within cells and between all domains genomewide, regardless of their replication timing. These results demonstrate that stochastic variation in replication timing is independent of elements that dictate timing or extrinsic environmental variation.

The epigenome and three-dimensional (3D) –genomic architecture are emerging as key factors in the dynamic regulation of different transcriptional programs required for neuronal functions. Here we utilize an activity-dependent tagging system in mice to determine the epigenetic state, 3D-genome architecture, and transcriptional landscape of engram cells over the lifespan of memory formation and recall. Our findings reveal that memory encoding leads to an epigenetic priming event, marked by increased accessibility of enhancers without corresponding transcriptional changes. Memory consolidation subsequently results in spatial reorganization of large chromatin segments and promoter-enhancer interactions. Finally, with reactivation, engram neurons utilize a subset of de novo long-range interactions, where primed enhancers were brought in contact with their respective promoters to up-regulate genes involved in local protein translation in synaptic compartments. Collectively, our work elucidates the comprehensive transcriptional and epigenomic landscape across the lifespan of memory formation and recall in the hippocampal engram ensemble. The formation and preservation of long-term memories depends on coordinated gene expression and synthesis of synaptic proteins 1 . These molecular processes act within a specific population of neurons, referred to as engram cells 2 – 4 . Recent approaches using activity-dependent expression of reporters, provided a framework for exploring the engram ensemble 5 – 8 , but the molecular mechanisms that govern memory storage and retrieval remain poorly understood. Specifically, epigenetic modifications and 3D -genomic architecture are emerging as a key factors in dynamic regulation of gene expression 9 – 17 , and there is an increasing appreciation of their importance in neuronal function, development and disease 14 , 16 , 18 Here, we utilized the Targeted Recombination in Active Populations (TRAP) mouse model 5 , 6 , in which activated neurons expressing the Activity Regulated Cytoskeleton Associated Protein, ( Arc ) gene, are permanently tagged in an inducible manner. Activated neurons during memory encoding, consolidation and recall were sorted and subjected to nuclear RNA sequencing (nRNA-seq), Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) and chromosome conformation capture (Hi-C). Our data demonstrates that memory encoding leads to a genome-wide increase in chromatin accessibility, without expected changes in gene expression. Furthermore, we demonstrate that late phase of memory consolidation was associated with re-localization of large chromatin segments (sub-compartments) from...