Global reorganization of replication domains during embryonic stem cell differentiation

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.

Cortical Development and Brain Malformations: Insights From the Differential Regulation of Early Events of DNA Replication

During the development of the cortex distinct populations of Neural Stem Cells (NSCs) are defined by differences in their cell cycle duration, self-renewal capacity and transcriptional profile. A key difference across the distinct populations of NSCs is the length of G1 phase, where the licensing of the DNA replication origins takes place by the assembly of a pre-replicative complex. Licensing of DNA replication is a process that is adapted accordingly to the cell cycle length of NSCs to secure the timed duplication of the genome. Moreover, DNA replication should be efficiently coordinated with ongoing transcription for the prevention of conflicts that would impede the progression of both processes, compromising the normal course of development. In the present review we discuss how the differential regulation of the licensing and initiation of DNA replication in different cortical NSCs populations is integrated with the properties of these stem cells populations. Moreover, we examine the implication of the initial steps of DNA replication in the pathogenetic mechanisms of neurodevelopmental defects and Zika virus-related microcephaly, highlighting the significance of the differential regulation of DNA replication during brain development.

Origin Firing Regulations to Control Genome Replication Timing

Complete genome duplication is essential for genetic homeostasis over successive cell generations. Higher eukaryotes possess a complex genome replication program that involves replicating the genome in units of individual chromatin domains with a reproducible order or timing. Two types of replication origin firing regulations ensure complete and well-timed domain-wise genome replication: (1) the timing of origin firing within a domain must be determined and (2) enough origins must fire with appropriate positioning in a short time window to avoid inter-origin gaps too large to be fully copied. Fundamental principles of eukaryotic origin firing are known. We here discuss advances in understanding the regulation of origin firing to control firing time. Work with yeasts suggests that eukaryotes utilise distinct molecular pathways to determine firing time of distinct sets of origins, depending on the specific requirements of the genomic regions to be replicated. Although the exact nature of the timing control processes varies between eukaryotes, conserved aspects exist: (1) the first step of origin firing, pre-initiation complex (pre-IC formation), is the regulated step, (2) many regulation pathways control the firing kinase Dbf4-dependent kinase, (3) Rif1 is a conserved mediator of late origin firing and (4) competition between origins for limiting firing factors contributes to firing timing. Characterization of the molecular timing control pathways will enable us to manipulate them to address the biological role of replication timing, for example, in cell differentiation and genome instability.

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

Here, we report a single-cell DNA replication sequencing method, scRepli-seq, a genome-wide methodology that measures copy number differences between replicated and unreplicated DNA. Using scRepli-seq, we demonstrate that replication-domain organization is conserved among individual mouse embryonic stem cells (mESCs). Differentiated mESCs exhibited distinct profiles, which were also conserved among cells. Haplotype-resolved scRepli-seq revealed similar replication profiles of homologous autosomes, while the inactive X chromosome was clearly replicated later than its active counterpart. However, a small degree of cell-to-cell replication-timing heterogeneity was present, which was smallest at the beginning and the end of S phase. In addition, developmentally regulated domains were found to deviate from others and showed a higher degree of heterogeneity, thus suggesting a link to developmental plasticity. Moreover, allelic expression imbalance was found to strongly associate with replication-timing asynchrony. Our results form a foundation for single-cell-level understanding of DNA replication regulation and provide insights into three-dimensional genome organization.

Identifying cis Elements for Spatiotemporal Control of Mammalian DNA Replication

The 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 early-to-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.

CERENKOV2: improved detection of functional noncoding SNPs using data-space geometric features

BACKGROUND

We previously reported on CERENKOV, an approach for identifying regulatory single nucleotide polymorphisms (rSNPs) that is based on 246 annotation features. CERENKOV uses the xgboost classifier and is designed to be used to find causal noncoding SNPs in loci identified by genome-wide association studies (GWAS). We reported that CERENKOV has state-of-the-art performance (by two traditional measures and a novel GWAS-oriented measure, AVGRANK) in a comparison to nine other tools for identifying functional noncoding SNPs, using a comprehensive reference SNP set (OSU17, 15,331 SNPs). Given that SNPs are grouped within loci in the reference SNP set and given the importance of the data-space manifold geometry for machine-learning model selection, we hypothesized that within-locus inter-SNP distances would have class-based distributional biases that could be exploited to improve rSNP recognition accuracy. We thus defined an intralocus SNP "radius" as the average data-space distance from a SNP to the other intralocus neighbors, and explored radius likelihoods for five distance measures.

RESULTS

We expanded the set of reference SNPs to 39,083 (the OSU18 set) and extracted CERENKOV SNP feature data. We computed radius empirical likelihoods and likelihood densities for rSNPs and control SNPs, and found significant likelihood differences between rSNPs and control SNPs. We fit parametric models of likelihood distributions for five different distance measures to obtain ten log-likelihood features that we combined with the 248-dimensional CERENKOV feature matrix. On the OSU18 SNP set, we measured the classification accuracy of CERENKOV with and without the new distance-based features, and found that the addition of distance-based features significantly improves rSNP recognition performance as measured by AUPVR, AUROC, and AVGRANK. Along with feature data for the OSU18 set, the software code for extracting the base feature matrix, estimating ten distance-based likelihood ratio features, and scoring candidate causal SNPs, are released as open-source software CERENKOV2.

CONCLUSIONS

Accounting for the locus-specific geometry of SNPs in data-space significantly improved the accuracy with which noncoding rSNPs can be computationally identified.

Replication timing and epigenome remodelling are associated with the nature of chromosomal rearrangements in cancer

DNA replication timing is known to facilitate the establishment of the epigenome, however, the intimate connection between replication timing and changes to the genome and epigenome in cancer remain largely uncharacterised. Here, we perform Repli-Seq and integrated epigenome analyses and demonstrate that genomic regions that undergo long-range epigenetic deregulation in prostate cancer also show concordant differences in replication timing. A subset of altered replication timing domains are conserved across cancers from different tissue origins. Notably, late-replicating regions in cancer cells display a loss of DNA methylation, and a switch in heterochromatin features from H3K9me3-marked constitutive to H3K27me3-marked facultative heterochromatin. Finally, analysis of 214 prostate and 35 breast cancer genomes reveal that late-replicating regions are prone to cis and early-replication to trans chromosomal rearrangements. Together, our data suggests that the nature of chromosomal rearrangement in cancer is related to the spatial and temporal positioning and altered epigenetic states of early-replicating compared to late-replicating loci.

DNA fragility in the parallel evolution of pelvic reduction in stickleback fish

Evolution generates a remarkable breadth of living forms, but many traits evolve repeatedly, by mechanisms that are still poorly understood. A classic example of repeated evolution is the loss of pelvic hindfins in stickleback fish (Gasterosteus aculeatus). Repeated pelvic loss maps to recurrent deletions of a pelvic enhancer of the Pitx1 gene. Here, we identify molecular features contributing to these recurrent deletions. Pitx1 enhancer sequences form alternative DNA structures in vitro and increase double-strand breaks and deletions in vivo. Enhancer mutability depends on DNA replication direction and is caused by TG-dinucleotide repeats. Modeling shows that elevated mutation rates can influence evolution under demographic conditions relevant for sticklebacks and humans. DNA fragility may thus help explain why the same loci are often used repeatedly during parallel adaptive evolution.

Linking the organization of DNA replication with genome maintenance

The spatial and temporal organization of genome duplication, also referred to as the replication program, is defined by the distribution and the activities of the sites of replication initiation across the genome. Alterations to the replication profile are associated with cell fate changes during development and in pathologies, but the importance of undergoing S phase with distinct and specific programs remains largely unexplored. We have recently addressed this question, focusing on the interplay between the replication program and genome maintenance. In particular, we demonstrated that when cells encounter challenges to DNA synthesis, the organization of DNA replication drives the response to replication stress that is mediated by the ATR/Rad3 checkpoint pathway, thus shaping the pattern of genome instability along the chromosomes. In this review, we present the major findings of our study and discuss how they may bring new perspectives to our understanding of the biological importance of the replication program.

Common fragile site instability in normal cells: Lessons and perspectives

Mechanisms and events related to common fragile site (CFS) instability are well known in cancer cells. Here, we argue that normal cells remain an important experimental model to address questions related to CFS instability in the absence of alterations in cell cycle and DNA damage repair pathways, which are common features acquired in cancer. Furthermore, a major gap of knowledge concerns the stability of CFSs during gametogenesis. CFS instability in meiotic or postmeiotic stages of the germ cell line could generate chromosome deletions or large rearrangements. This in turn can lead to the functional loss of the several CFS-associated genes with tumor suppressor function. Our hypothesis is that such mutations can potentially result in genetic predisposition to develop cancer. Indirect evidence for CFS instability in human germ cells has been provided by genomic investigations in family pedigrees associated with genetic disease. The issue of CFS instability in the germ cell line should represent one of the future efforts, and may take advantage of the existence of sequence and functional conservation of CFSs between rodents and humans.

DNA mismatch repair preferentially safeguards actively transcribed genes

DNA mismatch repair (MMR) is an evolutionally conserved genome maintenance pathway and is well known for its role in maintaining replication fidelity by correcting biosynthetic errors generated during DNA replication. However, recent studies have shown that MMR preferentially protects actively transcribed genes from mutation during both DNA replication and transcription. This review describes the recent discoveries in this area. Potential mechanisms by which MMR safeguards actively transcribed genes are also discussed.

Chromatin conformation and transcriptional activity are permissive regulators of DNA replication initiation in Drosophila

Chromatin structure has emerged as a key contributor to spatial and temporal control over the initiation of DNA replication. However, despite genome-wide correlations between early replication of gene-rich, accessible euchromatin and late replication of gene-poor, inaccessible heterochromatin, a causal relationship between chromatin structure and replication initiation remains elusive. Here, we combined histone gene engineering and whole-genome sequencing in Drosophila to determine how perturbing chromatin structure affects replication initiation. We found that most pericentric heterochromatin remains late replicating in H3K9R mutants, even though H3K9R pericentric heterochromatin is depleted of HP1a, more accessible, and transcriptionally active. These data indicate that HP1a loss, increased chromatin accessibility, and elevated transcription do not result in early replication of heterochromatin. Nevertheless, a small amount of pericentric heterochromatin with increased accessibility replicates earlier in H3K9R mutants. Transcription is de-repressed in these regions of advanced replication but not in those regions of the H3K9R mutant genome that replicate later, suggesting that transcriptional repression may contribute to late replication. We also explored relationships among chromatin, transcription, and replication in euchromatin by analyzing H4K16R mutants. In Drosophila, the X Chromosome gene expression is up-regulated twofold and replicates earlier in XY males than it does in XX females. We found that H4K16R mutation prevents normal male development and abrogates hyperexpression and earlier replication of the male X, consistent with previously established genome-wide correlations between transcription and early replication. In contrast, H4K16R females are viable and fertile, indicating that H4K16 modification is dispensable for genome replication and gene expression.

Mapping mammalian replication domains using the ion torrent semiconductor sequencing platform

Here, we show that semiconductor-based sequencing technology can be used to map mammalian replication domains, chromosomal units with similar DNA replication timing. Replicating DNA purified from mammalian cells was successfully sequenced by the Ion Torrent platform. The resultant replication domain map of mouse embryonic stem cells is comparable to those obtained by the conventional microarray-based method.

Dynamic changes in ORC localization and replication fork progression during tissue differentiation

Background

Genomic regions repressed for DNA replication, resulting in either delayed replication in S phase or underreplication in polyploid cells, are thought to be controlled by inhibition of replication origin activation. Studies in Drosophila polytene cells, however, raised the possibility that impeding replication fork progression also plays a major role.

Results

We exploited genomic regions underreplicated (URs) with tissue specificity in Drosophila polytene cells to analyze mechanisms of replication repression. By localizing the Origin Recognition Complex (ORC) in the genome of the larval fat body and comparing this to ORC binding in the salivary gland, we found that sites of ORC binding show extensive tissue specificity. In contrast, there are common domains nearly devoid of ORC in the salivary gland and fat body that also have reduced density of ORC binding sites in diploid cells. Strikingly, domains lacking ORC can still be replicated in some polytene tissues, showing absence of ORC and origins is insufficient to repress replication. Analysis of the width and location of the URs with respect to ORC position indicates that whether or not a genomic region lacking ORC is replicated is controlled by whether replication forks formed outside the region are inhibited.

Conclusions

These studies demonstrate that inhibition of replication fork progression can block replication across genomic regions that constitutively lack ORC. Replication fork progression can be inhibited in both tissue-specific and genome region-specific ways. Consequently, when evaluating sources of genome instability it is important to consider altered control of replication forks in response to differentiation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4992-3) contains supplementary material, which is available to authorized users.

Cellular senescence induces replication stress with almost no affect on DNA replication timing

Organismal aging entails a gradual decline of normal physiological functions and a major contributor to this decline is withdrawal of the cell cycle, known as senescence. Senescence can result from telomere diminution leading to a finite number of population doublings, known as replicative senescence (RS), or from oncogene overexpression, as a protective mechanism against cancer. Senescence is associated with large-scale chromatin re-organization and changes in gene expression. Replication stress is a complex phenomenon, defined as the slowing or stalling of replication fork progression and/or DNA synthesis, which has serious implications for genome stability, and consequently in human diseases. Aberrant replication fork structures activate the replication stress response leading to the activation of dormant origins, which is thought to be a safeguard mechanism to complete DNA replication on time. However, the relationship between replicative stress and the changes in the spatiotemporal program of DNA replication in senescence progression remains unclear. Here, we studied the DNA replication program during senescence progression in proliferative and pre-senescent cells from donors of various ages by single DNA fiber combing of replicated DNA, origin mapping by sequencing short nascent strands and genome-wide profiling of replication timing (TRT). We demonstrate that, progression into RS leads to reduced replication fork rates and activation of dormant origins, which are the hallmarks of replication stress. However, with the exception of a delay in RT of the CREB5 gene in all pre-senescent cells, RT was globally unaffected by replication stress during entry into either oncogene-induced or RS. Consequently, we conclude that RT alterations associated with physiological and accelerated aging, do not result from senescence progression. Our results clarify the interplay between senescence, aging and replication programs and demonstrate that RT is largely resistant to replication stress.

MCM2 promotes symmetric inheritance of modified histones during DNA replication

During genome replication, parental histones are recycled to newly replicated DNA with their posttranslational modifications (PTMs). Whether sister chromatids inherit modified histones evenly remains unknown. We measured histone PTM partition to sister chromatids in embryonic stem cells. We found that parental histones H3-H4 segregate to both daughter DNA strands with a weak leading-strand bias, skewing partition at topologically associating domain (TAD) borders and enhancers proximal to replication initiation zones. Segregation of parental histones to the leading strand increased markedly in cells with histone-binding mutations in MCM2, part of the replicative helicase, exacerbating histone PTM sister chromatid asymmetry. This work reveals how histones are inherited to sister chromatids and identifies a mechanism by which the replication machinery ensures symmetric cell division.

Modeling three-dimensional genomic organization in evolution and pathogenesis

The regulation of gene expression is mediated via the complex three-dimensional (3D) conformation of the genetic material and its interactions with various intracellular factors. Various experimental and computational approaches have been developed in recent years for understating the relation between the 3D conformation of the genome and the phenotypes of cells in normal condition and diseases. In this review, we will discuss novel approaches for analyzing and modeling the 3D genomic conformation, focusing on deciphering disease-causing mutations that affect gene expression. We conclude that as this is a very challenging mission, an important direction should involve the comparative analysis of various 3D models from various organisms or cells.

The organization of genome duplication is a critical determinant of the landscape of genome maintenance

Genome duplication is essential for cell proliferation, and the mechanisms regulating its execution are highly conserved. These processes give rise to a spatiotemporal organization of replication initiation across the genome, referred to as the replication program. Despite the identification of such programs in diverse eukaryotic organisms, their biological importance for cellular physiology remains largely unexplored. We address this fundamental question in the context of genome maintenance, taking advantage of the inappropriate origin firing that occurs when fission yeast cells lacking the Rad3/ATR checkpoint kinase are subjected to replication stress. Using this model, we demonstrate that the replication program quantitatively dictates the extent of origin de-regulation and the clustered localization of these events. Furthermore, our results uncover an accumulation of abnormal levels of single-stranded DNA (ssDNA) and the Rad52 repair protein at de-regulated origins. We show that these loci constitute a defining source of the overall ssDNA and Rad52 hotspots in the genome, generating a signature pattern of instability along the chromosomes. We then induce a genome-wide reprogramming of origin usage and evaluate its consequences in our experimental system. This leads to a complete redistribution of the sites of both inappropriate initiation and associated Rad52 recruitment. We therefore conclude that the organization of genome duplication governs the checkpoint control of origin-associated hotspots of instability and plays an integral role in shaping the landscape of genome maintenance.

Regulation of the program of DNA replication by CDK: new findings and perspectives

Progression through the cell cycle is driven by the activities of the cyclin-dependent kinase (CDK) family of enzymes, which establish an ordered passage through the cell cycle phases. CDK activity is crucial for the cellular transitions from G1 to S and G2 to M, which are highly controlled to promote the faithful duplication of the genetic material and the transmission of the genome into daughter cells, respectively. While oscillations in CDK activity are essential for cell division, how its specific dynamics may shape cellular processes remains an open question. Recently, we have investigated the potential role of CDK in establishing the profile of replication initiation along the chromosomes, also referred to as the replication program. Our results demonstrated that the timing and level of CDK activity at G1/S provide two critical and independent inputs that modulate the pattern of origin usage. In this review, we will present the conclusions of our study and discuss the implications of our findings for cellular function and physiology.

Replication timing and nuclear structure

DNA replication proceeds along spatially and temporally coordinated patterns within the nucleus, thus protecting the genome during the synthesis of new genetic material. While we have been able to visualize replication patterns on DNA fibers for 50 years, recent developments and discoveries have provided a greater insight into how DNA replication is controlled. In this review, we highlight many of these discoveries. Of great interest are the physiological role of the replication timing program, cis and trans-acting factors that modulate replication timing and the effects of chromatin structure on the replication timing program. We also discuss future directions in the study of replication timing.

Rif1 prolongs the embryonic S phase at the Drosophila mid-blastula transition

In preparation for dramatic morphogenetic events of gastrulation, rapid embryonic cell cycles slow at the mid-blastula transition (MBT). In Drosophila melanogaster embryos, down-regulation of cyclin-dependent kinase 1 (Cdk1) activity initiates this slowing by delaying replication of heterochromatic satellite sequences and extending S phase. We found that Cdk1 activity inhibited the chromatin association of Rap1 interacting factor 1 (Rif1), a candidate repressor of replication. Furthermore, Rif1 bound selectively to satellite sequences following Cdk1 down-regulation at the MBT. In the next S phase, Rif1 dissociated from different satellites in an orderly schedule that anticipated their replication. Rif1 lacking potential phosphorylation sites failed to dissociate and dominantly prevented completion of replication. Loss of Rif1 in mutant embryos shortened the post-MBT S phase and rescued embryonic cell cycles disrupted by depletion of the S phase–promoting kinase, cell division cycle 7 (Cdc7). Our work shows that Rif1 and S phase kinases compose a replication timer controlling first the developmental onset of late replication and then the precise schedule of replication within S phase. In addition, we describe how onset of late replication fits into the progressive maturation of heterochromatin during development.

Cells divide rapidly in the early embryos of most animals. However, during a conserved period of development known as the mid-blastula transition (MBT), the cell cycle slows down dramatically. In Drosophila embryos, genome duplication abruptly slows to initiate this cell cycle prolongation. This is achieved through the onset of late replication, a well-recognized phenomenon in which specific sequences of the genome await replication until long after other sequences have finished. Even though this temporal program of replication is a major determinant of the duration of S phase, the factors involved in this process remain unknown. Here, we use genetics and real-time microscopy to visualize replication in developing fly embryos and show that the protein Rap1 interacting factor 1 (Rif1) mediates the introduction of late replication at the MBT. We find that at this stage, Rif1 binds to and selectively delays the replication of large blocks of repetitive DNA known as satellite sequences. During the rapid cell cycles before the MBT, we show that the cyclin-dependent kinase 1 (Cdk1) prevents Rif1 from slowing down DNA replication by driving its removal from the chromatin. The developmental down-regulation of Cdk1 at the MBT allows Rif1 to associate with the satellite sequences and initiate cell cycle slowing. Our work provides new insights into the temporal programming of S phase and into the embryonic origin of late replication.

Profiling DNA Replication Timing Using Zebrafish as an In Vivo Model System

DNA replication timing is an important cellular characteristic, exhibiting significant relationships with chromatin structure, transcription, and DNA mutation rates. Changes in replication timing occur during development and in cancer, but the role replication timing plays in development and disease is not known. Zebrafish were recently established as an in vivo model system to study replication timing. Here is detailed the protocols for using the zebrafish to determine DNA replication timing. After sorting cells from embryos and adult zebrafish, high-resolution genome-wide DNA replication timing patterns can be constructed by determining changes in DNA copy number through analysis of next generation sequencing data. The zebrafish model system allows for evaluation of the replication timing changes that occur in vivo throughout development, and can also be used to assess changes in individual cell types, disease models, or mutant lines. These methods will enable studies investigating the mechanisms and determinants of replication timing establishment and maintenance during development, the role replication timing plays in mutations and tumorigenesis, and the effects of perturbing replication timing on development and disease.

Chromatin conformation regulates the coordination between DNA replication and transcription

Chromatin is the template for the basic processes of replication and transcription, making the maintenance of chromosomal integrity critical for cell viability. To elucidate how dividing cells respond to alterations in chromatin structure, here we analyse the replication programme of primary cells with altered chromatin configuration caused by the genetic ablation of the HMGB1 gene, or three histone H1 genes. We find that loss of chromatin compaction in H1-depleted cells triggers the accumulation of stalled forks and DNA damage as a consequence of transcription–replication conflicts. In contrast, reductions in nucleosome occupancy due to the lack of HMGB1 cause faster fork progression without impacting the initiation landscape or fork stability. Thus, perturbations in chromatin integrity elicit a range of responses in the dynamics of DNA replication and transcription, with different consequences on replicative stress. These findings have broad implications for our understanding of how defects in chromatin structure contribute to genomic instability.

The maintenance of chromatin integrity during replication is critical for cell viability. Here the authors study how dividing cells respond to alterations in chromatin structure and find that these elicit a range of responses in the dynamics of DNA replication and consequences on replicative stress.

Srf destabilizes cellular identity by suppressing cell-type-specific gene expression programs

Multicellular organisms consist of multiple cell types. The identity of these cells is primarily maintained by cell-type-specific gene expression programs; however, mechanisms that suppress these programs are poorly defined. Here we show that serum response factor (Srf), a transcription factor that is activated by various extracellular stimuli, can repress cell-type-specific genes and promote cellular reprogramming to pluripotency. Manipulations that decrease β-actin monomer quantity result in the nuclear accumulation of Mkl1 and the activation of Srf, which downregulate cell-type-specific genes and alter the epigenetics of regulatory regions and chromatin organization. Mice overexpressing Srf exhibit various pathologies including an ulcerative colitis-like symptom and a metaplasia-like phenotype in the pancreas. Our results demonstrate an unexpected function of Srf via a mechanism by which extracellular stimuli actively destabilize cell identity and suggest Srf involvement in a wide range of diseases.

The transcription factor Srf is a central regulator of immediate-early and actin cytoskeletal genes. Here the authors show that Srf is activated by repression of β-actin, promoting iPSC reprogramming of neural progenitor cells and hepatoblasts by repressing cell-type specific genes.

H3K36me3-mediated mismatch repair preferentially protects actively transcribed genes from mutation

Histone H3 trimethylation at lysine 36 (H3K36me3) is an important histone mark involved in both transcription elongation and DNA mismatch repair (MMR). It is known that H3K36me3 recruits the mismatch-recognition protein MutSα to replicating chromatin via its physical interaction with MutSα's PWWP domain, but the exact role of H3K36me3 in transcription is undefined. Using ChIP combined with whole-genome DNA sequencing analysis, we demonstrate here that H3K36me3, together with MutSα, is involved in protecting against mutation, preferentially in actively transcribed genomic regions. We found that H3K36me3 and MutSα are much more co-enriched in exons and actively transcribed regions than in introns and nontranscribed regions. The H3K36me3-MutSα co-enrichment correlated with a much lower mutation frequency in exons and actively transcribed regions than in introns and nontranscribed regions. Correspondingly, depleting H3K36me3 or disrupting the H3K36me3-MutSα interaction elevated the spontaneous mutation frequency in actively transcribed genes, but it had little influence on the mutation frequency in nontranscribed or transcriptionally inactive regions. Similarly, H₂O₂-induced mutations, which mainly cause base oxidations, preferentially occurred in actively transcribed genes in MMR-deficient cells. The data presented here suggest that H3K36me3-mediated MMR preferentially safeguards actively transcribed genes not only during replication by efficiently correcting mispairs in early replicating chromatin but also during transcription by directly or indirectly removing DNA lesions associated with a persistently open chromatin structure.

Genome-wide analysis of replication timing by next-generation sequencing with E/L Repli-seq

This protocol is an extension to: Nat. Protoc. 6, 870-895 (2014); doi:10.1038/nprot.2011.328; published online 02 June 2011Cycling cells duplicate their DNA content during S phase, following a defined program called replication timing (RT). Early- and late-replicating regions differ in terms of mutation rates, transcriptional activity, chromatin marks and subnuclear position. Moreover, RT is regulated during development and is altered in diseases. Here, we describe E/L Repli-seq, an extension of our Repli-chip protocol. E/L Repli-seq is a rapid, robust and relatively inexpensive protocol for analyzing RT by next-generation sequencing (NGS), allowing genome-wide assessment of how cellular processes are linked to RT. Briefly, cells are pulse-labeled with BrdU, and early and late S-phase fractions are sorted by flow cytometry. Labeled nascent DNA is immunoprecipitated from both fractions and sequenced. Data processing leads to a single bedGraph file containing the ratio of nascent DNA from early versus late S-phase fractions. The results are comparable to those of Repli-chip, with the additional benefits of genome-wide sequence information and an increased dynamic range. We also provide computational pipelines for downstream analyses, for parsing phased genomes using single-nucleotide polymorphisms (SNPs) to analyze RT allelic asynchrony, and for direct comparison to Repli-chip data. This protocol can be performed in up to 3 d before sequencing, and requires basic cellular and molecular biology skills, as well as a basic understanding of Unix and R.

Genome-Wide Analysis of the Arabidopsis Replication Timing Program

Eukaryotes use a temporally regulated process, known as the replication timing program, to ensure that their genomes are fully and accurately duplicated during S phase. Replication timing programs are predictive of genomic features and activity and are considered to be functional readouts of chromatin organization. Although replication timing programs have been described for yeast and animal systems, much less is known about the temporal regulation of plant DNA replication or its relationship to genome sequence and chromatin structure. We used the thymidine analog, 5-ethynyl-2'-deoxyuridine, in combination with flow sorting and Repli-Seq to describe, at high-resolution, the genome-wide replication timing program for Arabidopsis (Arabidopsis thaliana) Col-0 suspension cells. We identified genomic regions that replicate predominantly during early, mid, and late S phase, and correlated these regions with genomic features and with data for chromatin state, accessibility, and long-distance interaction. Arabidopsis chromosome arms tend to replicate early while pericentromeric regions replicate late. Early and mid-replicating regions are gene-rich and predominantly euchromatic, while late regions are rich in transposable elements and primarily heterochromatic. However, the distribution of chromatin states across the different times is complex, with each replication time corresponding to a mixture of states. Early and mid-replicating sequences interact with each other and not with late sequences, but early regions are more accessible than mid regions. The replication timing program in Arabidopsis reflects a bipartite genomic organization with early/mid-replicating regions and late regions forming separate, noninteracting compartments. The temporal order of DNA replication within the early/mid compartment may be modulated largely by chromatin accessibility.

CDK activity provides temporal and quantitative cues for organizing genome duplication

In eukaryotes, the spatial and temporal organization of genome duplication gives rise to distinctive profiles of replication origin usage along the chromosomes. While it has become increasingly clear that these programs are important for cellular physiology, the mechanisms by which they are determined and modulated remain elusive. Replication initiation requires the function of cyclin-dependent kinases (CDKs), which associate with various cyclin partners to drive cell proliferation. Surprisingly, although we possess detailed knowledge of the CDK regulators and targets that are crucial for origin activation, little is known about whether CDKs play a critical role in establishing the genome-wide pattern of origin selection. We have addressed this question in the fission yeast, taking advantage of a simplified cell cycle network in which cell proliferation is driven by a single cyclin-CDK module. This system allows us to precisely control CDK activity in vivo using chemical genetics. First, in contrast to previous reports, our results clearly show that distinct cyclin-CDK pairs are not essential for regulating specific subsets of origins and for establishing a normal replication program. Importantly, we then demonstrate that the timing at which CDK activity reaches the S phase threshold is critical for the organization of replication in distinct efficiency domains, while the level of CDK activity at the onset of S phase is a dose-dependent modulator of overall origin efficiencies. Our study therefore implicates these different aspects of CDK regulation as versatile mechanisms for shaping the architecture of DNA replication across the genome.

The duplication of the genetic material is a highly conserved and tightly regulated process that is essential for cell proliferation. DNA synthesis initiates at sites called origins that are distributed throughout the genome. Replication origins are not all used equivalently, and their patterns of activation along the chromosomes give rise to specific profiles, or programs, of DNA replication. These programs change during development and in response to external stimuli, and we have previously shown that they have important consequences for cellular function. However, we still do not understand the mechanisms by which cells establish different replication patterns. Here we investigate the role of the cyclin-dependent kinase (CDK) family of proteins, whose activities are critical for cell cycle progression, in regulating the organization of genome duplication. Taking advantage of a system that allows us to precisely modulate CDK activity levels in living cells, we demonstrate that both the temporal and quantitative controls of CDK function are crucial for determining distinct programs of DNA replication. Our work therefore uncovers a fundamental link between CDK activity, a central input in a variety of cellular and developmental processes, and the architecture of genome duplication.

Three classes of recurrent DNA break clusters in brain progenitors identified by 3D proximity-based break joining assay

We recently discovered 27 recurrent DNA double-strand break (DSB) clusters (RDCs) in mouse neural stem/progenitor cells (NSPCs). Most RDCs occurred across long, late-replicating RDC genes and were found only after mild inhibition of DNA replication. RDC genes share intriguing characteristics, including encoding surface proteins that organize brain architecture and neuronal junctions, and are genetically implicated in neuropsychiatric disorders and/or cancers. RDC identification relies on high-throughput genome-wide translocation sequencing (HTGTS), which maps recurrent DSBs based on their translocation to "bait" DSBs in specific chromosomal locations. Cellular heterogeneity in 3D genome organization allowed unequivocal identification of RDCs on 14 different chromosomes using HTGTS baits on three mouse chromosomes. Additional candidate RDCs were also implicated, however, suggesting that some RDCs were missed. To more completely identify RDCs, we exploited our finding that joining of two DSBs occurs more frequently if they lie on the same cis chromosome. Thus, we used CRISPR/Cas9 to introduce specific DSBs into each mouse chromosome in NSPCs that were used as bait for HTGTS libraries. This analysis confirmed all 27 previously identified RDCs and identified many new ones. NSPC RDCs fall into three groups based on length, organization, transcription level, and replication timing of genes within them. While mostly less robust, the largest group of newly defined RDCs share many intriguing characteristics with the original 27. Our findings also revealed RDCs in NSPCs in the absence of induced replication stress, and support the idea that the latter treatment augments an already active endogenous process.

Single-cell replication profiling to measure stochastic variation in mammalian replication timing

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.

[G-quadruplex: key controllers of human genome duplication]

The correct duplication of the human genome is under the control of a spatiotemporal program that determines where and when replication forks start. This regulation thus mainly operates on replication start sites named replication origins. During the S-phase, about 50 000 origins fire in one human cell. However, the normal or perturbed progression of replication forks also strongly impacts on replication. Recently, several studies have put forward the role of a noncanonical DNA structure, the G-quadruplex, in the control of genome duplication. In this review, we describe the major impact of this structure on starting points and on the progression of replication forks.

Binding of high mobility group A proteins to the mammalian genome occurs as a function of AT-content

Genomic location can inform on potential function and recruitment signals for chromatin-associated proteins. High mobility group (Hmg) proteins are of similar size as histones with Hmga1 and Hmga2 being particularly abundant in replicating normal tissues and in cancerous cells. While several roles for Hmga proteins have been proposed we lack a comprehensive description of their genomic location as a function of chromatin, DNA sequence and functional domains. Here we report such a characterization in mouse embryonic stem cells in which we introduce biotin-tagged constructs of wild-type and DNA-binding domain mutants. Comparative analysis of the genome-wide distribution of Hmga proteins reveals pervasive binding, a feature that critically depends on a functional DNA-binding domain and which is shared by both Hmga proteins. Assessment of the underlying queues instructive for this binding modality identifies AT richness, defined as high frequency of A or T bases, as the major criterion for local binding. Additionally, we show that other chromatin states such as those linked to cis-regulatory regions have little impact on Hmga binding both in stem and differentiated cells. As a consequence, Hmga proteins are preferentially found at AT-rich regions such as constitutively heterochromatic regions but are absent from enhancers and promoters arguing for a limited role in regulating individual genes. In line with this model, we show that genetic deletion of Hmga proteins in stem cells causes limited transcriptional effects and that binding is conserved in neuronal progenitors. Overall our comparative study describing the in vivo binding modality of Hmga1 and Hmga2 identifies the proteins’ preference for AT-rich DNA genome-wide and argues against a suggested function of Hmga at regulatory regions. Instead we discover pervasive binding with enrichment at regions of higher AT content irrespective of local variation in chromatin modifications.

We investigated the chromosomal location of a group of highly abundant nuclear proteins. Our genome-wide results for Hmga1 and Hmga2 reveal a unique binding modality indicating preference for DNA rich in A or T bases in vivo. Importantly this preferential binding to AT-rich sequences occurs throughout the genome irrespectively of other local chromatin features. Genomic location and loss of function experiments challenge the view that Hmga proteins act as local modulators of transcriptional regulation but rather argue for a role as structural components of chromatin.

The Three-Dimensional Organization of Mammalian Genomes

Animal development depends on not only the linear genome sequence that embeds millions of cis-regulatory elements, but also the three-dimensional (3D) chromatin architecture that orchestrates the interplay between cis-regulatory elements and their target genes. Compared to our knowledge of the cis-regulatory sequences, the understanding of the 3D genome organization in human and other eukaryotes is still limited. Recent advances in technologies to map the 3D genome architecture have greatly accelerated the pace of discovery. Here, we review emerging concepts of chromatin organization in mammalian cells, discuss the dynamics of chromatin conformation during development, and highlight important roles for chromatin organization in cancer and other human diseases.

The 4D nucleome project

The 4D Nucleome Network aims to develop and apply approaches to map the structure and dynamics of the human and mouse genomes in space and time with the goal of gaining deeper mechanistic insights into how the nucleus is organized and functions. The project will develop and benchmark experimental and computational approaches for measuring genome conformation and nuclear organization, and investigate how these contribute to gene regulation and other genome functions. Validated experimental technologies will be combined with biophysical approaches to generate quantitative models of spatial genome organization in different biological states, both in cell populations and in single cells.

Genomic Analysis of the DNA Replication Timing Program during Mitotic S Phase in Maize (Zea mays) Root Tips

All plants and animals must replicate their DNA, using a regulated process to ensure that their genomes are completely and accurately replicated. DNA replication timing programs have been extensively studied in yeast and animal systems, but much less is known about the replication programs of plants. We report a novel adaptation of the "Repli-seq" assay for use in intact root tips of maize (Zea mays) that includes several different cell lineages and present whole-genome replication timing profiles from cells in early, mid, and late S phase of the mitotic cell cycle. Maize root tips have a complex replication timing program, including regions of distinct early, mid, and late S replication that each constitute between 20 and 24% of the genome, as well as other loci corresponding to ∼32% of the genome that exhibit replication activity in two different time windows. Analyses of genomic, transcriptional, and chromatin features of the euchromatic portion of the maize genome provide evidence for a gradient of early replicating, open chromatin that transitions gradually to less open and less transcriptionally active chromatin replicating in mid S phase. Our genomic level analysis also demonstrated that the centromere core replicates in mid S, before heavily compacted classical heterochromatin, including pericentromeres and knobs, which replicate during late S phase.

Repliscan: a tool for classifying replication timing regions

BACKGROUND

Replication timing experiments that use label incorporation and high throughput sequencing produce peaked data similar to ChIP-Seq experiments. However, the differences in experimental design, coverage density, and possible results make traditional ChIP-Seq analysis methods inappropriate for use with replication timing.

RESULTS

To accurately detect and classify regions of replication across the genome, we present Repliscan. Repliscan robustly normalizes, automatically removes outlying and uninformative data points, and classifies Repli-seq signals into discrete combinations of replication signatures. The quality control steps and self-fitting methods make Repliscan generally applicable and more robust than previous methods that classify regions based on thresholds.

CONCLUSIONS

Repliscan is simple and effective to use on organisms with different genome sizes. Even with analysis window sizes as small as 1 kilobase, reliable profiles can be generated with as little as 2.4x coverage.

Dormant origins as a built-in safeguard in eukaryotic DNA replication against genome instability and disease development

DNA replication is a prerequisite for cell proliferation, yet it can be increasingly challenging for a eukaryotic cell to faithfully duplicate its genome as its size and complexity expands. Dormant origins now emerge as a key component for cells to successfully accomplish such a demanding but essential task. In this perspective, we will first provide an overview of the fundamental processes eukaryotic cells have developed to regulate origin licensing and firing. With a special focus on mammalian systems, we will then highlight the role of dormant origins in preventing replication-associated genome instability and their functional interplay with proteins involved in the DNA damage repair response for tumor suppression. Lastly, deficiencies in the origin licensing machinery will be discussed in relation to their influence on stem cell maintenance and human diseases.

Single-cell multi-omics sequencing of mouse early embryos and embryonic stem cells

Single-cell epigenome sequencing techniques have recently been developed. However, the combination of different layers of epigenome sequencing in an individual cell has not yet been achieved. Here, we developed a single-cell multi-omics sequencing technology (single-cell COOL-seq) that can analyze the chromatin state/nucleosome positioning, DNA methylation, copy number variation and ploidy simultaneously from the same individual mammalian cell. We used this method to analyze the reprogramming of the chromatin state and DNA methylation in mouse preimplantation embryos. We found that within < 12 h of fertilization, each individual cell undergoes global genome demethylation together with the rapid and global reprogramming of both maternal and paternal genomes to a highly opened chromatin state. This was followed by decreased openness after the late zygote stage. Furthermore, from the late zygote to the 4-cell stage, the residual DNA methylation is preferentially preserved on intergenic regions of the paternal alleles and intragenic regions of maternal alleles in each individual blastomere. However, chromatin accessibility is similar between paternal and maternal alleles in each individual cell from the late zygote to the blastocyst stage. The binding motifs of several pluripotency regulators are enriched at distal nucleosome depleted regions from as early as the 2-cell stage. This indicates that the cis-regulatory elements of such target genes have been primed to an open state from the 2-cell stage onward, long before pluripotency is eventually established in the ICM of the blastocyst. Genes may be classified into homogeneously open, homogeneously closed and divergent states based on the chromatin accessibility of their promoter regions among individual cells. This can be traced to step-wise transitions during preimplantation development. Our study offers the first single-cell and parental allele-specific analysis of the genome-scale chromatin state and DNA methylation dynamics at single-base resolution in early mouse embryos and provides new insights into the heterogeneous yet highly ordered features of epigenomic reprogramming during this process.

A stitch in time: Replicate early and escape dosage compensation to express more

María Gómez previews work by Müller and Nieduszynski that describes a biological role for conserved replication timing differences in eukaryotic cells.

The biological significance of conserved replication timing patterns in eukaryotic genomes remains a mystery. In this issue, Müller and Nieduszynski (2017. J. Cell Biol.https://doi.org/10.1083/jcb.201701061) find that early replication is a requirement for the highest expression levels of certain genes.