Characterization of the replication timing program of 6 human model cell lines

The location and development of Replicon Cluster Domains in early replicating DNA

Background: It has been known for many years that in metazoan cells, replication origins are organised into clusters where origins within each cluster fire near-synchronously. Despite clusters being a fundamental organising principle of metazoan DNA replication, the genomic location of origin clusters has not been documented. Methods: We synchronised human U2OS by thymidine block and release followed by L-mimosine block and release to create a population of cells progressing into S phase with a high degree of synchrony. At different times after release into S phase, cells were pulsed with EdU; the EdU-labelled DNA was then pulled down, sequenced and mapped onto the human genome. Results: The early replicating DNA showed features at a range of scales. Wavelet analysis showed that the major feature of the early replicating DNA was at a size of 500 kb, consistent with clusters of replication origins. Over the first two hours of S phase, these Replicon Cluster Domains broadened in width, consistent with their being enlarged by the progression of replication forks at their outer boundaries. The total replication signal associated with each Replicon Cluster Domain varied considerably, and this variation was reproducible and conserved over time. We provide evidence that this variability in replication signal was at least in part caused by Replicon Cluster Domains being activated at different times in different cells in the population. We also provide evidence that adjacent clusters had a statistical preference for being activated in sequence across a group, consistent with the 'domino' model of replication focus activation order observed by microscopy. Conclusions: We show that early replicating DNA is organised into Replicon Cluster Domains that behave as expected of replicon clusters observed by DNA fibre analysis. The coordinated activation of different Replicon Cluster Domains can generate the replication timing programme by which the genome is duplicated.

Omics Analyses: How to Navigate Through a Constant Data Deluge

Omics data are very valuable for researchers in biology, but the work required to develop a solid expertise in their analysis contrasts with the rapidity with which the omics technologies evolve. Data accumulate in public databases, and despite significant advances in bioinformatics softwares to integrate them, data analysis remains a burden for those who perform experiments. Beyond the issue of dealing with a very large number of results, we believe that working with omics data requires a change in the way scientific problems are solved. In this chapter, we explain pitfalls and tips we found during our functional genomics projects in yeasts. Our main lesson is that, if applying a protocol does not guarantee a successful project, following simple rules can help to become strategic and intentional, thus avoiding an endless drift into an ocean of possibilities.

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.

Low Replicative Stress Triggers Cell-Type Specific Inheritable Advanced Replication Timing

DNA replication timing (RT), reflecting the temporal order of origin activation, is known as a robust and conserved cell-type specific process. Upon low replication stress, the slowing of replication forks induces well-documented RT delays associated to genetic instability, but it can also generate RT advances that are still uncharacterized. In order to characterize these advanced initiation events, we monitored the whole genome RT from six independent human cell lines treated with low doses of aphidicolin. We report that RT advances are cell-type-specific and involve large heterochromatin domains. Importantly, we found that some major late to early RT advances can be inherited by the unstressed next-cellular generation, which is a unique process that correlates with enhanced chromatin accessibility, as well as modified replication origin landscape and gene expression in daughter cells. Collectively, this work highlights how low replication stress may impact cellular identity by RT advances events at a subset of chromosomal domains.

Replication Stress, Genomic Instability, and Replication Timing: A Complex Relationship

The replication-timing program constitutes a key element of the organization and coordination of numerous nuclear processes in eukaryotes. This program is established at a crucial moment in the cell cycle and occurs simultaneously with the organization of the genome, thus indicating the vital significance of this process. With recent technological achievements of high-throughput approaches, a very strong link has been confirmed between replication timing, transcriptional activity, the epigenetic and mutational landscape, and the 3D organization of the genome. There is also a clear relationship between replication stress, replication timing, and genomic instability, but the extent to which they are mutually linked to each other is unclear. Recent evidence has shown that replication timing is affected in cancer cells, although the cause and consequence of this effect remain unknown. However, in-depth studies remain to be performed to characterize the molecular mechanisms of replication-timing regulation and clearly identify different cis- and trans-acting factors. The results of these studies will potentially facilitate the discovery of new therapeutic pathways, particularly for personalized medicine, or new biomarkers. This review focuses on the complex relationship between replication timing, replication stress, and genomic instability.

Cohesin-Mediated Genome Architecture Does Not Define DNA Replication Timing Domains

3D genome organization is strongly predictive of DNA replication timing in mammalian cells. This work tested the extent to which loop-based genome architecture acts as a regulatory unit of replication timing by using an auxin-inducible system for acute cohesin ablation. Cohesin ablation in a population of cells in asynchronous culture was shown not to disrupt patterns of replication timing as assayed by replication sequencing (RepliSeq) or BrdU-focus microscopy. Furthermore, cohesin ablation prior to S phase entry in synchronized cells was similarly shown to not impact replication timing patterns. These results suggest that cohesin-mediated genome architecture is not required for the execution of replication timing patterns in S phase, nor for the establishment of replication timing domains in G1.

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.

Temporal association of ORCA/LRWD1 to late-firing origins during G1 dictates heterochromatin replication and organization

DNA replication requires the recruitment of a pre-replication complex facilitated by Origin Recognition Complex (ORC) onto the chromatin during G1 phase of the cell cycle. The ORC-associated protein (ORCA/LRWD1) stabilizes ORC on chromatin. Here, we evaluated the genome-wide distribution of ORCA using ChIP-seq during specific time points of G1. ORCA binding sites on the G1 chromatin are dynamic and temporally regulated. ORCA association to specific genomic sites decreases as the cells progressed towards S-phase. The majority of the ORCA-bound sites represent replication origins that also associate with the repressive chromatin marks H3K9me3 and methylated-CpGs, consistent with ORCA-bound origins initiating DNA replication late in S-phase. Further, ORCA directly associates with the repressive marks and interacts with the enzymes that catalyze these marks. Regions that associate with both ORCA and H3K9me3, exhibit diminished H3K9 methylation in ORCA-depleted cells, suggesting a role for ORCA in recruiting the H3K9me3 mark at certain genomic loci. Similarly, DNA methylation is altered at ORCA-occupied sites in cells lacking ORCA. Furthermore, repressive chromatin marks influence ORCA's binding on chromatin. We propose that ORCA coordinates with the histone and DNA methylation machinery to establish a repressive chromatin environment at a subset of origins, which primes them for late replication.