The replication domain model: regulating replicon firing in the context of large-scale chromosome architecture

RIF1 Links Replication Timing with Fork Reactivation and DNA Double-Strand Break Repair

Replication timing (RT) is a cellular program to coordinate initiation of DNA replication in all origins within the genome. RIF1 (replication timing regulatory factor 1) is a master regulator of RT in human cells. This role of RIF1 is associated with binding G4-quadruplexes and changes in 3D chromatin that may suppress origin activation over a long distance. Many effects of RIF1 in fork reactivation and DNA double-strand (DSB) repair (DSBR) are underlined by its interaction with TP53BP1 (tumor protein p53 binding protein). In G1, RIF1 acts antagonistically to BRCA1 (BRCA1 DNA repair associated), suppressing end resection and homologous recombination repair (HRR) and promoting non-homologous end joining (NHEJ), contributing to DSBR pathway choice. RIF1 is an important element of intra-S-checkpoints to recover damaged replication fork with the involvement of HRR. High-resolution microscopic studies show that RIF1 cooperates with TP53BP1 to preserve 3D structure and epigenetic markers of genomic loci disrupted by DSBs. Apart from TP53BP1, RIF1 interact with many other proteins, including proteins involved in DNA damage response, cell cycle regulation, and chromatin remodeling. As impaired RT, DSBR and fork reactivation are associated with genomic instability, a hallmark of malignant transformation, RIF1 has a diagnostic, prognostic, and therapeutic potential in cancer. Further studies may reveal other aspects of common regulation of RT, DSBR, and fork reactivation by RIF1.

The Temporal Order of DNA Replication Shaped by Mammalian DNA Methyltransferases

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

DNA replication and chromosome positioning throughout the interphase in three-dimensional space of plant nuclei

Despite much recent progress, our understanding of the principles of plant genome organization and its dynamics in three-dimensional space of interphase nuclei remains surprisingly limited. Notably, it is not clear how these processes could be affected by the size of a plant's nuclear genome. In this study, DNA replication timing and interphase chromosome positioning were analyzed in seven Poaceae species that differ in their genome size. To provide a comprehensive picture, a suite of advanced, complementary methods was used: labeling of newly replicated DNA by ethynyl-2'-deoxyuridine, isolation of nuclei at particular cell cycle phases by flow cytometric sorting, three-dimensional immunofluorescence in situ hybridization, and confocal microscopy. Our results revealed conserved dynamics of DNA replication in all species, and a similar replication timing order for telomeres and centromeres, as well as for euchromatin and heterochromatin regions, irrespective of genome size. Moreover, stable chromosome positioning was observed while transitioning through different stages of interphase. These findings expand upon earlier studies in suggesting that a more complex interplay exists between genome size, organization of repetitive DNA sequences along chromosomes, and higher order chromatin structure and its maintenance in interphase, albeit controlled by currently unknown factors.

Genomic methods for measuring DNA replication dynamics

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density, and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analog incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

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

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

Peripheral re-localization of constitutive heterochromatin advances its replication timing and impairs maintenance of silencing marks

The replication of the genome is a highly organized process, both spatially and temporally. Although a lot is known on the composition of the basic replication machinery, how its activity is regulated is mostly unknown. Several chromatin properties have been proposed as regulators, but a potential role of the nuclear DNA position remains unclear. We made use of the prominent structure and well-defined heterochromatic landscape of mouse pericentric chromosome domains as a well-studied example of late replicating constitutive heterochromatin. We established a method to manipulate its nuclear position and evaluated the effect on replication timing, DNA compaction and epigenetic composition. Using time-lapse microscopy, we observed that constitutive heterochromatin, known to replicate during late S-phase, was replicated in mid S-phase when repositioned to the nuclear periphery. Out-of-schedule replication resulted in deficient post-replicative maintenance of chromatin modifications, namely silencing marks. We propose that repositioned constitutive heterochromatin was activated in trans according to the domino model of origin firing by nearby (mid S) firing origins. In summary, our data provide, on the one hand, a novel approach to manipulate nuclear DNA position and, on the other hand, establish nuclear DNA position as a novel mechanism regulating DNA replication timing and epigenetic maintenance.

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

MOTIVATION

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

RESULTS

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

AVAILABILITY AND IMPLEMENTATION

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

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Bayesian inference of origin firing time distributions, origin interference and licencing probabilities from Next Generation Sequencing data

DNA replication is a stochastic process with replication forks emanating from multiple replication origins. The origins must be licenced in G1, and the replisome activated at licenced origins in order to generate bi-directional replication forks in S-phase. Differential firing times lead to origin interference, where a replication fork from an origin can replicate through and inactivate neighbouring origins (origin obscuring). We developed a Bayesian algorithm to characterize origin firing statistics from Okazaki fragment (OF) sequencing data. Our algorithm infers the distributions of firing times and the licencing probabilities for three consecutive origins. We demonstrate that our algorithm can distinguish partial origin licencing and origin obscuring in OF sequencing data from Saccharomyces cerevisiae and human cell types. We used our method to analyse the decreased origin efficiency under loss of Rat1 activity in S. cerevisiae, demonstrating that both reduced licencing and increased obscuring contribute. Moreover, we show that robust analysis is possible using only local data (across three neighbouring origins), and analysis of the whole chromosome is not required. Our algorithm utilizes an approximate likelihood and a reversible jump sampling technique, a methodology that can be extended to analysis of other mechanistic processes measurable through Next Generation Sequencing data.

ChECing out Rif1 action in freely cycling cells

In buddying yeast, like all eukaryotes examined so far, DNA replication is under temporal control, such that some origins fire early and some late during S phase. This replication timing program is established in G1 phase, where chromatin states are thought to prevent binding of key-limiting initiation factors at late-firing origins. Although many factors are involved in replication initiation, a new player, Rif1, has recently entered the scene, with a spate of papers revealing a global role for the protein in the control of replication initiation timing from yeasts to humans. Since budding yeast Rif1 was known to bind only to telomeric and silent mating loci regions, it remained controversial whether Rif1 acts directly at replication origins or instead influences origin activity indirectly. In this perspective, we discuss our recent finding that Rif1 binds directly to the replication origins that it controls. In this study, we also found that Rif1's regulatory activity at origins is best revealed by an assay (sort-seq) that measures replication in unperturbed, freely cycling cultures, as opposed to commonly used protocols in which cells are first blocked in the G1 phase of the cell cycle by mating pheromone, then released into a synchronous S phase. Finally, we discuss how the sequestration of Rif1 at telomeres, through an interaction with the arrays of Rap1 molecules bound there, plays an important role in limiting Rif1's action primarily to telomere-proximal replication origins.

DNA Replication: From Radioisotopes to Click Chemistry

The replication of nuclear and mitochondrial DNA are basic processes assuring the doubling of the genetic information of eukaryotic cells. In research of the basic principles of DNA replication, and also in the studies focused on the cell cycle, an important role is played by artificially-prepared nucleoside and nucleotide analogues that serve as markers of newly synthesized DNA. These analogues are incorporated into the DNA during DNA replication, and are subsequently visualized. Several methods are used for their detection, including the highly popular click chemistry. This review aims to provide the readers with basic information about the various possibilities of the detection of replication activity using nucleoside and nucleotide analogues, and to show the strengths and weaknesses of those different detection systems, including click chemistry for microscopic studies.

Histone H4K20 tri-methylation at late-firing origins ensures timely heterochromatin replication

Among other targets, the protein lysine methyltransferase PR-Set7 induces histone H4 lysine 20 monomethylation (H4K20me1), which is the substrate for further methylation by the Suv4-20h methyltransferase. Although these enzymes have been implicated in control of replication origins, the specific contribution of H4K20 methylation to DNA replication remains unclear. Here, we show that H4K20 mutation in mammalian cells, unlike in Drosophila, partially impairs S-phase progression and protects from DNA re-replication induced by stabilization of PR-Set7. Using Epstein-Barr virus-derived episomes, we further demonstrate that conversion of H4K20me1 to higher H4K20me2/3 states by Suv4-20h is not sufficient to define an efficient origin per se, but rather serves as an enhancer for MCM2-7 helicase loading and replication activation at defined origins. Consistent with this, we find that Suv4-20h-mediated H4K20 tri-methylation (H4K20me3) is required to sustain the licensing and activity of a subset of ORCA/LRWD1-associated origins, which ensure proper replication timing of late-replicating heterochromatin domains. Altogether, these results reveal Suv4-20h-mediated H4K20 tri-methylation as a critical determinant in the selection of active replication initiation sites in heterochromatin regions of mammalian genomes.

A journey through the microscopic ages of DNA replication

Scientific discoveries and technological advancements are inseparable but not always take place in a coherent chronological manner. In the next, we will provide a seemingly unconnected and serendipitous series of scientific facts that, in the whole, converged to unveil DNA and its duplication. We will not cover here the many and fundamental contributions from microbial genetics and in vitro biochemistry. Rather, in this journey, we will emphasize the interplay between microscopy development culminating on super resolution fluorescence microscopy (i.e., nanoscopy) and digital image analysis and its impact on our understanding of DNA duplication. We will interlace the journey with landmark concepts and experiments that have brought the cellular DNA replication field to its present state.

Defining multiple, distinct, and shared spatiotemporal patterns of DNA replication and endoreduplication from 3D image analysis of developing maize (Zea mays L.) root tip nuclei

Spatiotemporal patterns of DNA replication have been described for yeast and many types of cultured animal cells, frequently after cell cycle arrest to aid in synchronization. However, patterns of DNA replication in nuclei from plants or naturally developing organs remain largely uncharacterized. Here we report findings from 3D quantitative analysis of DNA replication and endoreduplication in nuclei from pulse-labeled developing maize root tips. In both early and middle S phase nuclei, flow-sorted on the basis of DNA content, replicative labeling was widely distributed across euchromatic regions of the nucleoplasm. We did not observe the perinuclear or perinucleolar replicative labeling patterns characteristic of middle S phase in mammals. Instead, the early versus middle S phase patterns in maize could be distinguished cytologically by correlating two quantitative, continuous variables, replicative labeling and DAPI staining. Early S nuclei exhibited widely distributed euchromatic labeling preferentially localized to regions with weak DAPI signals. Middle S nuclei also exhibited widely distributed euchromatic labeling, but the label was preferentially localized to regions with strong DAPI signals. Highly condensed heterochromatin, including knobs, replicated during late S phase as previously reported. Similar spatiotemporal replication patterns were observed for both mitotic and endocycling maize nuclei. These results revealed that maize euchromatin exists as an intermingled mixture of two components distinguished by their condensation state and replication timing. These different patterns might reflect a previously described genome organization pattern, with "gene islands" mostly replicating during early S phase followed by most of the intergenic repetitive regions replicating during middle S phase.

Genome architecture: from linear organisation of chromatin to the 3D assembly in the nucleus

The genetic information is stored in the eukaryotic nucleus in the form of chromatin. This is a macromolecular entity that includes genomic DNA and histone proteins that form nucleosomes, plus a large variety of chromatin-associated non-histone proteins. Chromatin is structurally and functionally organised at various levels. One reveals the linear topography of DNA, histones and their post-translational modifications and non-histone proteins along each chromosome. This level provides regulatory information about the association of genomic elements with particular signatures that have been used to define chromatin states. Importantly, these chromatin states correlate with structural and functional genomic features. Another regulatory layer is established at the level of the 3D organisation of chromatin within the nucleus, which has been revealed clearly as non-random. Instead, a variety of intra- and inter-chromosomal genomic domains with specific epigenetic and functional properties has been identified. In this review, we discuss how the recent advances in genomic approaches have contributed to our understanding of these two levels of genome architecture. We have emphasised our analysis with the aim of integrating information available for yeast, Arabidopsis, Drosophila, and mammalian cells. We consider that this comparative study helps define common and unique features in each system, providing a basis to better understand the complexity of genome organisation.

Predicting chromatin organization using histone marks

Genome-wide mapping of three dimensional chromatin organization is an important yet technically challenging task. To aid experimental effort and to understand the determinants of long-range chromatin interactions, we have developed a computational model integrating Hi-C and histone mark ChIP-seq data to predict two important features of chromatin organization: chromatin interaction hubs and topologically associated domain (TAD) boundaries. Our model accurately and robustly predicts these features across datasets and cell types. Cell-type specific histone mark information is required for prediction of chromatin interaction hubs but not for TAD boundaries. Our predictions provide a useful guide for the exploration of chromatin organization.

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

Duplication of the genome in mammalian cells occurs in a defined temporal order referred to as its replication-timing (RT) program. RT changes dynamically during development, regulated in units of 400-800 kb referred to as replication domains (RDs). Changes in RT are generally coordinated with transcriptional competence and changes in subnuclear position. We generated genome-wide RT profiles for 26 distinct human cell types, including embryonic stem cell (hESC)-derived, primary cells and established cell lines representing intermediate stages of endoderm, mesoderm, ectoderm, and neural crest (NC) development. We identified clusters of RDs that replicate at unique times in each stage (RT signatures) and confirmed global consolidation of the genome into larger synchronously replicating segments during differentiation. Surprisingly, transcriptome data revealed that the well-accepted correlation between early replication and transcriptional activity was restricted to RT-constitutive genes, whereas two-thirds of the genes that switched RT during differentiation were strongly expressed when late replicating in one or more cell types. Closer inspection revealed that transcription of this class of genes was frequently restricted to the lineage in which the RT switch occurred, but was induced prior to a late-to-early RT switch and/or down-regulated after an early-to-late RT switch. Analysis of transcriptional regulatory networks showed that this class of genes contains strong regulators of genes that were only expressed when early replicating. These results provide intriguing new insight into the complex relationship between transcription and RT regulation during human development.

Structural organization of human replication timing domains

Recent analysis of genome-wide epigenetic modification data, mean replication timing (MRT) profiles and chromosome conformation data in mammals have provided increasing evidence that flexibility in replication origin usage is regulated locally by the epigenetic landscape and over larger genomic distances by the 3D chromatin architecture. Here, we review the recent results establishing some link between replication domains and chromatin structural domains in pluripotent and various differentiated cell types in human. We reconcile the originally proposed dichotomic picture of early and late constant timing regions that replicate by multiple rather synchronous origins in separated nuclear compartments of open and closed chromatins, with the U-shaped MRT domains bordered by "master" replication origins specified by a localized (∼200-300 kb) zone of open and transcriptionally active chromatin from which a replication wave likely initiates and propagates toward the domain center via a cascade of origin firing. We discuss the relationships between these MRT domains, topologically associated domains and lamina-associated domains. This review sheds a new light on the epigenetically regulated global chromatin reorganization that underlies the loss of pluripotency and the determination of differentiation properties.

Depletion of the chromatin looping proteins CTCF and cohesin causes chromatin compaction: insight into chromatin folding by polymer modelling

Folding of the chromosomal fibre in interphase nuclei is an important element in the regulation of gene expression. For instance, physical contacts between promoters and enhancers are a key element in cell-type-specific transcription. We know remarkably little about the principles that control chromosome folding. Here we explore the view that intrachromosomal interactions, forming a complex pattern of loops, are a key element in chromosome folding. CTCF and cohesin are two abundant looping proteins of interphase chromosomes of higher eukaryotes. To investigate the role of looping in large-scale (supra Mb) folding of human chromosomes, we knocked down the gene that codes for CTCF and the one coding for Rad21, an essential subunit of cohesin. We measured the effect on chromosome folding using systematic 3D fluorescent in situ hybridization (FISH). Results show that chromatin becomes more compact after reducing the concentration of these two looping proteins. The molecular basis for this counter-intuitive behaviour is explored by polymer modelling usingy the Dynamic Loop model (Bohn M, Heermann DW (2010) Diffusion-driven looping provides a consistent framework for chromatin organization. PLoS ONE 5: e12218.). We show that compaction can be explained by selectively decreasing the number of short-range loops, leaving long-range looping unchanged. In support of this model prediction it has recently been shown by others that CTCF and cohesin indeed are responsible primarily for short-range looping. Our results suggest that the local and the overall changes in of chromosome structure are controlled by a delicate balance between short-range and long-range loops, allowing easy switching between, for instance, open and more compact chromatin states.

Replication initiation and genome instability: a crossroads for DNA and RNA synthesis

Nuclear DNA replication requires the concerted action of hundreds of proteins to efficiently unwind and duplicate the entire genome while also retaining epigenetic regulatory information. Initiation of DNA replication is tightly regulated, rapidly firing thousands of origins once the conditions to promote rapid and faithful replication are in place, and defects in replication initiation lead to proliferation defects, genome instability, and a range of developmental abnormalities. Interestingly, DNA replication in metazoans initiates in actively transcribed DNA, meaning that replication initiation occurs in DNA that is co-occupied with tens of thousands of poised and active RNA polymerase complexes. Active transcription can induce genome instability, particularly during DNA replication, as RNA polymerases can induce torsional stress, formation of secondary structures, and act as a physical barrier to other enzymes involved in DNA metabolism. Here we discuss the challenges facing mammalian DNA replication, their impact on genome instability, and the development of cancer.

Chromatin structure and replication origins: determinants of chromosome replication and nuclear organization

The DNA replication program is, in part, determined by the epigenetic landscape that governs local chromosome architecture and directs chromosome duplication. Replication must coordinate with other biochemical processes occurring concomitantly on chromatin, such as transcription and remodeling, to insure accurate duplication of both genetic and epigenetic features and to preserve genomic stability. The importance of genome architecture and chromatin looping in coordinating cellular processes on chromatin is illustrated by two recent sets of discoveries. First, chromatin-associated proteins that are not part of the core replication machinery were shown to affect the timing of DNA replication. These chromatin-associated proteins could be working in concert, or perhaps in competition, with the transcriptional machinery and with chromatin modifiers to determine the spatial and temporal organization of replication initiation events. Second, epigenetic interactions are mediated by DNA sequences that determine chromosomal replication. In this review, we summarize recent findings and current models linking spatial and temporal regulation of the replication program with epigenetic signaling. We discuss these issues in the context of the genome's three-dimensional structure with an emphasis on events occurring during the initiation of DNA replication.

Allele-specific genome-wide profiling in human primary erythroblasts reveal replication program organization

We have developed a new approach to characterize allele-specific timing of DNA replication genome-wide in human primary basophilic erythroblasts. We show that the two chromosome homologs replicate at the same time in about 88% of the genome and that large structural variants are preferentially associated with asynchronous replication. We identified about 600 megabase-sized asynchronously replicated domains in two tested individuals. The longest asynchronously replicated domains are enriched in imprinted genes suggesting that structural variants and parental imprinting are two causes of replication asynchrony in the human genome. Biased chromosome X inactivation in one of the two individuals tested was another source of detectable replication asynchrony. Analysis of high-resolution TimEX profiles revealed small variations termed timing ripples, which were undetected in previous, lower resolution analyses. Timing ripples reflect highly reproducible, variations of the timing of replication in the 100 kb-range that exist within the well-characterized megabase-sized replication timing domains. These ripples correspond to clusters of origins of replication that we detected using novel nascent strands DNA profiling methods. Analysis of the distribution of replication origins revealed dramatic differences in initiation of replication frequencies during S phase and a strong association, in both synchronous and asynchronous regions, between origins of replication and three genomic features: G-quadruplexes, CpG Islands and transcription start sites. The frequency of initiation in asynchronous regions was similar in the two homologs. Asynchronous regions were richer in origins of replication than synchronous regions.

DNA replication in mammalian cells proceeds according to a distinct order. Genes that are expressed tend to replicate before genes that are not expressed. We report here that we have developed a method to measure the timing of replication of the maternal and paternal chromosomes separately. We found that the paternal and maternal chromosomes replicate at exactly the same time in the large majority of the genome and that the 12% of the genome that replicated asynchronously was enriched in imprinted genes and in structural variants. Previous experiments have shown that chromosomes could be divided into replication timing domains that are a few hundred thousand to a few megabases in size. We show here that these domains can be divided into sub-domains defined by ripples in the timing profile. These ripples corresponded to clusters of origins of replication. Finally, we show that the frequency of initiation in asynchronous regions was similar in the two homologs.