Chromatin structure fully determines replication timing program in human cells

Eukaryotic DNA replication follows a strict temporal program where genomic loci are replicated at precise times during the S phase of the cell cycle. Yet, the mechanism in control of the timing program in metazoan cells is poorly understood. In a recent publication, the authors proposed an intuitive stochastic model of DNA replication and showed that it predicts replication timing with an accuracy approaching the level of experimental biological repeats. Here, we discuss an extended software implementation of the mechanistic model: Replicon. This package allows interested researchers to predict the global replication timing program in human cells from chromatin data.

Recent evidence suggests that the timing of DNA replication is coordinated across megabase-scale domains in metazoan genomes, yet the importance of this aspect of genome organization is unclear. Here we show that replication timing is remarkably conserved between human and mouse, uncovering large regions that may have been governed by similar replication dynamics since these species have diverged. This conservation is both tissue-specific and independent of the genomic G+C content conservation. Moreover, we show that time of replication is globally conserved despite numerous large-scale genome rearrangements. We systematically identify rearrangement fusion points and demonstrate that replication time can be locally diverged at these loci. Conversely, rearrangements are shown to be correlated with early replication and physical chromosomal proximity. These results suggest that large chromosomal domains of coordinated replication are shuffled by evolution while conserving the large-scale nuclear architecture of the genome.

The transmission of genetic information relies on a coordinated network of cell cycle controls. Abnormalities in this network can result in genomic instability and lead to the transformation of normal cells into cancer cells. Chromosomal DNA replication is not only central to cellular division but also plays a crucial role in the maintenance of genomic integrity. DNA replication errors increase genetic instability, and may be a causative factor in diseases such as cancer and neuronal disorders. Replication in eukaryotes initiates from discrete genomic regions, termed origins, according to a strict, often tissue-specific, temporal program. The genetic program that controls activation of replication origins in mammalian cells has still not been elucidated. There is evidence that specification of replication sites and timing of replication are dynamic processes that are regulated by tissue-specific and developmental cues and that are responsive to epigenetic modifications. Here, we focus on the spatiotemporal regulation of DNA replication in the human genome. There is growing evidence that chromosome band patterns and epigenetic transformation of chromatin influence the timing of replication. On the basis of this evidence, we propose that the chromatin regions showing switches in replication timing from early to late in S phase are correlated with chromosome band boundaries. These chromatin regions generally display transitions in GC contents and include more non-B-form DNA structures than other genomic regions. We also examine here the effect of changes in replication timing on genomic stability and the possible role of replication timing in the etiology of diseases such as cancer. Replication timing assays are one of many promising techniques under investigation that may in future allow much earlier cancer detection than is possible today.

Genomes in eukaryotes normally undergo DNA replication in a choreographed temporal order, resulting in early and late replicating chromosome compartments. Leishmania, a human protozoan parasite, displays an unconventional DNA replication program in which the timing of DNA replication completion is chromosome size-dependent: larger chromosomes complete replication later then smaller ones. Here we show that both R-loops and RNase H1, a ribonuclease that resolves RNA-DNA hybrids, accumulate in Leishmania major chromosomes in a pattern that reflects their replication timing. Furthermore, we demonstrate that such differential organisation of R-loops, RNase H1 and DNA replication timing across the parasite’s chromosomes correlates with size-dependent differences in chromatin accessibility, G quadruplex distribution and sequence content. Using conditional gene excision, we show that loss of RNase H1 leads to transient growth perturbation and permanently abrogates the differences in DNA replication timing across chromosomes, as well as altering levels of aneuploidy and increasing chromosome instability in a size-dependent manner. This work provides a link between R-loop homeostasis and DNA replication timing in a eukaryotic parasite and demonstrates that orchestration of DNA replication dictates levels of genome plasticity in Leishmania.

Replication timing in a single human chromosome 11 transferred into the Chinese hamster ovary (CHO) cell line

The organization of mammalian DNA replication is poorly understood. We have produced high-resolution dynamic maps of the timing of replication in human erythroid, mesenchymal, and embryonic stem (ES) cells using TimEX, a method that relies on gaussian convolution of massive, highly redundant determinations of DNA copy-number variations during S phase to produce replication timing profiles. We first obtained timing maps of 3% of the genome using high-density oligonucleotide tiling arrays and then extended the TimEX method genome-wide using massively parallel sequencing. We show that in untransformed human cells, timing of replication is highly regulated and highly synchronous, and that many genomic segments are replicated in temporal transition regions devoid of initiation, where replication forks progress unidirectionally from origins that can be hundreds of kilobases away. Absence of initiation in one transition region is shown at the molecular level by single molecule analysis of replicated DNA (SMARD). Comparison of ES and erythroid cells replication patterns revealed that these cells replicate about 20% of their genome in different quarters of S phase. Importantly, we detected a strong inverse relationship between timing of replication and distance to the closest expressed gene. This relationship can be used to predict tissue-specific timing of replication profiles from expression data and genomic annotations. We also provide evidence that early origins of replication are preferentially located near highly expressed genes, that mid-firing origins are located near moderately expressed genes, and that late-firing origins are located far from genes.

Dynamic Alterations of Replication Timing in Mammalian Cells

SUMMARYMammalian genomes replicate in a cell-type-specific order during the S phase, correlated to transcriptional activity, histone modifications, and chromatin structure. The causal relationships between these features and DNA replication timing (RT), especially during cell fate changes, are largely unknown. Using machine learning, we quantify 21 chromatin features predicting local RT and RT changes during differentiation in embryonic stem cells (ESCs). About one-third of the genome shows RT changes during differentiation. Chromatin features accurately predict both steady-state RT and RT changes. Histone H3 lysine 4 monomethylation (H3K4me1), catalyzed by KMT2C and KMT2D (KMT2C/D), emerges as a top predictor. Loss of KMT2C/D or their enzymatic activities impairs RT changes during differentiation. This correlates with local H3K4me1 loss and reduced replication origin firing, while transcription remains largely unaffected. Our findings reveal KMT2C/D-dependent H3K4me1 as a key regulator of RT and replication initiation, a role that likely impacts diseases associated with KMT2C/D mutations.

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.

Replication timing alterations in leukemia affect clinically relevant chromosome domains

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.

During the last decades it has been shown that the replication timing program in metazoans is related to chromosome structure, the nuclear positioning and AT/GC content of chromosomal loci, their patterns of histone modifications, and their transcriptional regulation. Here, the current state of knowledge concerning these relationships is reviewed. An integrated view on structure-function relationships in the nucleus is provided and the determination and functional role of the replication timing program is discussed in this context. A corresponding comprehensive model is developed and a key aspect of this model is the suggestion that mammalian chromosomes are organized into stable units equivalent to replicon clusters. It is proposed that the nuclear positions of these units would depend on their histone modifications and determine the replication timing of the whole unit. It is furthermore predicted that replication timing is only indirectly linked to transcriptional regulation and contributes to the maintenance of gene expression patterns. These clear predictions, and the fact that the tools are at hand now to further test them, open an avenue towards solving the long standing problem on how replication timing is determined in metazoan cells.

Three-dimensional nuclear organisation and the DNA replication timing program

The temporal order of DNA replication [replication timing (RT)] is correlated with chromatin modifications and three-dimensional genome architecture; however, causal links have not been established, largely because of an inability to manipulate the global RT program. We show that loss of RIF1 causes near-complete elimination of the RT program by increasing heterogeneity between individual cells. RT changes are coupled with widespread alterations in chromatin modifications and genome compartmentalization. Conditional depletion of RIF1 causes replication-dependent disruption of histone modifications and alterations in genome architecture. These effects were magnified with successive cycles of altered RT. These results support models in which the timing of chromatin replication and thus assembly plays a key role in maintaining the global epigenetic state.

Studies on DNA replication during synchronized vegetative growth and gametic differentiation in Chlamydomonas reinhardtii