Analysis Of Replication Timing Research Articles

Repli-seq Sample Preparation using Cell Sorting with Cell-Permeant Dyes.

Replication timing is significantly correlated with gene expression and chromatin organization, changes dynamically during cell differentiation, and is altered in diseased states. Genome-wide analysis of replication timing is performed in actively replicating cells by Repli-seq. Current methods for Repli-seq require cells to be fixed in large numbers. This is a barrier for sample types that are sensitive to fixation or are in very limited numbers. In this article, we outline optimized methods to process live cells and intact nuclei for Repli-seq. Our protocol enables the processing of a smaller number of cells per sample and reduces processing time and sample loss while obtaining high-quality data. Further, for samples that tend to form clumps and are difficult to dissociate into a single-cell suspension, we also outline methods for isolation, staining, and processing of nuclei for Repli-seq. The Repli-seq data obtained from live cells and intact nuclei are comparable to those obtained from the standard protocols. © 2023 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Live cell isolation and staining Alternate Protocol: Nuclei isolation and staining.

Epigenomic signatures associated with spontaneous and replication stress-induced DNA double strand breaks.

Common fragile sites (CFSs) are specific regions of all individuals' genome that are predisposed to DNA double strand breaks (DSBs) and undergo subsequent rearrangements. CFS formation can be induced in vitro by mild level of DNA replication stress, such as DNA polymerase inhibition or nucleotide pool disturbance. The mechanisms of CFS formation have been linked to DNA replication timing control, transcription activities, as well as chromatin organization. However, it is unclear what specific cis- or trans-factors regulate the interplay between replication and transcription that determine CFS formation. We recently reported genome-wide mapping of DNA DSBs under replication stress induced by aphidicolin in human lymphoblastoids for the first time. Here, we systematically compared these DSBs with regards to nearby epigenomic features mapped in the same cell line from published studies. We demonstrate that aphidicolin-induced DSBs are strongly correlated with histone 3 lysine 36 trimethylation, a marker for active transcription. We further demonstrate that this DSB signature is a composite effect by the dual treatment of aphidicolin and its solvent, dimethylsulfoxide, the latter of which potently induces transcription on its own. We also present complementing evidence for the association between DSBs and 3D chromosome architectural domains with high density gene cluster and active transcription. Additionally, we show that while DSBs were detected at all but one of the fourteen finely mapped CFSs, they were not enriched in the CFS core sequences and rather demarcated the CFS core region. Related to this point, DSB density was not higher in large genes of greater than 300kb, contrary to reported enrichment of CFS sites at these large genes. Finally, replication timing analyses demonstrate that the CFS core region contain initiation events, suggesting that altered replication dynamics are responsible for CFS formation in relatively higher level of replication stress.

Replication stress generates distinctive landscapes of DNA copy number alterations and chromosome scale losses

BackgroundA major driver of cancer chromosomal instability is replication stress, the slowing or stalling of DNA replication. How replication stress and genomic instability are connected is not known. Aphidicolin-induced replication stress induces breakages at common fragile sites, but the exact causes of fragility are debated, and acute genomic consequences of replication stress are not fully explored.ResultsWe characterize DNA copy number alterations (CNAs) in single, diploid non-transformed cells, caused by one cell cycle in the presence of either aphidicolin or hydroxyurea. Multiple types of CNAs are generated, associated with different genomic regions and features, and observed copy number landscapes are distinct between aphidicolin and hydroxyurea-induced replication stress. Coupling cell type-specific analysis of CNAs to gene expression and single-cell replication timing analyses pinpointed the causative large genes of the most recurrent chromosome-scale CNAs in aphidicolin. These are clustered on chromosome 7 in RPE1 epithelial cells but chromosome 1 in BJ fibroblasts. Chromosome arm level CNAs also generate acentric lagging chromatin and micronuclei containing these chromosomes.ConclusionsChromosomal instability driven by replication stress occurs via focal CNAs and chromosome arm scale changes, with the latter confined to a very small subset of chromosome regions, potentially heavily skewing cancer genome evolution. Different inducers of replication stress lead to distinctive CNA landscapes providing the opportunity to derive copy number signatures of specific replication stress mechanisms. Single-cell CNA analysis thus reveals the impact of replication stress on the genome, providing insights into the molecular mechanisms which fuel chromosomal instability in cancer.

Optimized Repli-seq: improved DNA replication timing analysis by next-generation sequencing.

The human genome is divided into functional units that replicate at specific times during S-phase. This temporal program is known as replication timing (RT) and is coordinated with the spatial organization of the genome and transcriptional activity. RT is also cell type-specific, dynamically regulated during development, and alterations in RT are observed in multiple diseases. Thus, the precise measure of RT is critical to understand the role of RT in gene function regulation. Distinct methods for assaying the RT program exist; however, conventional methods require thousands of cells as input, prohibiting its applicability to samples with limited cell numbers such as those from disease patients or from early developing embryos. Although single-cell RT analyses have been developed, these methods are low throughput, require generation of numerous libraries, increased sequencing costs, and produce low resolution data. Here, we developed an improved method to measure RT genome-wide that enables high-resolution analysis of low input samples. This method incorporates direct cell sorting into lysis buffer, as well as DNA fragmentation and library preparation in a single tube, resulting in higher yields, increased quality, and reproducibility with decreased costs. We also performed a systematic data processing analysis to provide standardized parameters for RT measurement. This optimized method facilitates RT analysis and will enable its application to a broad range of studies investigating the role of RT in gene expression, nuclear architecture, and disease.

Kronos scRT: a uniform framework for single-cell replication timing analysis

Mammalian genomes are replicated in a cell type-specific order and in coordination with transcription and chromatin organization. Currently, single-cell replication studies require individual processing of sorted cells, yielding a limited number (<100) of cells. Here, we develop Kronos scRT, a software for single-cell Replication Timing (scRT) analysis. Kronos scRT does not require a specific platform or cell sorting, which allows investigating large datasets obtained from asynchronous cells. By applying our tool to published data as well as droplet-based single-cell whole-genome sequencing data generated in this study, we exploit scRT from thousands of cells for different mouse and human cell lines. Our results demonstrate that although genomic regions are frequently replicated around their population average RT, replication can occur stochastically throughout S phase. Altogether, Kronos scRT allows fast and comprehensive investigations of the RT programme at the single-cell resolution for both homogeneous and heterogeneous cell populations.

SAF-A promotes origin licensing and replication fork progression to ensure robust DNA replication.

The organisation of chromatin is closely intertwined with biological activities of chromosome domains, including transcription and DNA replication status. Scaffold-attachment factor A (SAF-A), also known as heterogeneous nuclear ribonucleoprotein U (HNRNPU), contributes to the formation of open chromatin structure. Here, we demonstrate that SAF-A promotes the normal progression of DNA replication and enables resumption of replication after inhibition. We report that cells depleted of SAF-A show reduced origin licensing in G1 phase and, consequently, reduced origin activation frequency in S phase. Replication forks also progress less consistently in cells depleted of SAF-A, contributing to reduced DNA synthesis rate. Single-cell replication timing analysis revealed two distinct effects of SAF-A depletion: first, the boundaries between early- and late-replicating domains become more blurred; and second, SAF-A depletion causes replication timing changes that tend to bring regions of discordant domain compartmentalisation and replication timing into concordance. Associated with these defects, SAF-A-depleted cells show elevated formation of phosphorylated histone H2AX (γ-H2AX) and tend to enter quiescence. Overall, we find that SAF-A protein promotes robust DNA replication to ensure continuing cell proliferation.

Replication timing analysis in polyploid cells reveals Rif1 uses multiple mechanisms to promote underreplication in Drosophila.

Regulation of DNA replication and copy number is necessary to promote genome stability and maintain cell and tissue function. DNA replication is regulated temporally in a process known as replication timing (RT). Rap1-interacting factor 1 (Rif1) is a key regulator of RT and has a critical function in copy number control in polyploid cells. Previously, we demonstrated that Rif1 functions with SUUR to inhibit replication fork progression and promote underreplication (UR) of specific genomic regions. How Rif1-dependent control of RT factors into its ability to promote UR is unknown. By applying a computational approach to measure RT in Drosophila polyploid cells, we show that SUUR and Rif1 have differential roles in controlling UR and RT. Our findings reveal that Rif1 acts to promote late replication, which is necessary for SUUR-dependent underreplication. Our work provides new insight into the process of UR and its links to RT.

Orc4 spatiotemporally stabilizes centromeric chromatin.

The establishment of centromeric chromatin and its propagation by the centromere-specific histone CENPA is mediated by epigenetic mechanisms in most eukaryotes. DNA replication origins, origin binding proteins, and replication timing of centromere DNA are important determinants of centromere function. The epigenetically regulated regional centromeres in the budding yeast Candida albicans have unique DNA sequences that replicate earliest in every chromosome and are clustered throughout the cell cycle. In this study, the genome-wide occupancy of the replication initiation protein Orc4 reveals its abundance at all centromeres in C. albicans. Orc4 is associated with four different DNA sequence motifs, one of which coincides with tRNA genes (tDNA) that replicate early and cluster together in space. Hi-C combined with genome-wide replication timing analyses identify that early replicating Orc4-bound regions interact with themselves stronger than with late replicating Orc4-bound regions. We simulate a polymer model of chromosomes of C. albicans and propose that the early replicating and highly enriched Orc4-bound sites preferentially localize around the clustered kinetochores. We also observe that Orc4 is constitutively localized to centromeres, and both Orc4 and the helicase Mcm2 are essential for cell viability and CENPA stability in C. albicans. Finally, we show that new molecules of CENPA are recruited to centromeres during late anaphase/telophase, which coincides with the stage at which the CENPA-specific chaperone Scm3 localizes to the kinetochore. We propose that the spatiotemporal localization of Orc4 within the nucleus, in collaboration with Mcm2 and Scm3, maintains centromeric chromatin stability and CENPA recruitment in C. albicans.

Efficient, quick and easy-to-use DNA replication timing analysis with START-R suite.

DNA replication must be faithful and follow a well-defined spatiotemporal program closely linked to transcriptional activity, epigenomic marks, intranuclear structures, mutation rate and cell fate determination. Among the readouts of the spatiotemporal program of DNA replication, replication timing analyses require not only complex and time-consuming experimental procedures, but also skills in bioinformatics. We developed a dedicated Shiny interactive web application, the START-R (Simple Tool for the Analysis of the Replication Timing based on R) suite, which analyzes DNA replication timing in a given organism with high-throughput data. It reduces the time required for generating and analyzing simultaneously data from several samples. It automatically detects different types of timing regions and identifies significant differences between two experimental conditions in ∼15 min. In conclusion, START-R suite allows quick, efficient and easier analyses of DNA replication timing for all organisms. This novel approach can be used by every biologist. It is now simpler to use this method in order to understand, for example, whether ‘a favorite gene or protein’ has an impact on replication process or, indirectly, on genomic organization (as Hi-C experiments), by comparing the replication timing profiles between wild-type and mutant cell lines.

How do CD8+ T cells complete DNA replication in less than 4 hours?

Abstract We demonstrate that activated CD8+ T cells from Pmel mice replicate their genome within 4 hours which is significantly faster than most cancer cells. Many genotoxic chemotherapeutic agents target the rapid proliferation of cancer cells. However, activated T cells also expand rapidly and this can lead to lymphopenia in patients receiving these chemotherapeutic agents. Understanding how CD8+ T cells complete DNA replication in less than 4 hours could lead to the development of treatment schedules that spare CD8+ T cells. Three possible mechanisms may underlie the rapid DNA replication in CD8+ T cells: 1) faster replication fork velocity, 2) higher density of origin firing, or 3) modified replication timing program, allowing for simultaneous replication of early- and late-replicating regions of the genome. Using DNA combing, we show that the replication fork velocity and inter-origin distance is similar in unperturbed CD8+ T cells and B16 cancer cells. Repli-seq analysis, a NGS genome-wide replication timing analysis, suggest that regions of the genome replicated late in S phase in fibroblasts are replicated earlier in CD8+ T cells, suggesting that the rapid DNA replication in CD8+ T cells is a consequence of a modified, truncated replication timing program. We also show for the first time that actively replicating CD8+ T cells have dormant origins suggesting that this mechanism of genome stability is maintained in activated CD8+ T cells.

Transcription-mediated organization of the replication initiation program across large genes sets common fragile sites genome-wide

Common fragile sites (CFSs) are chromosome regions prone to breakage upon replication stress known to drive chromosome rearrangements during oncogenesis. Most CFSs nest in large expressed genes, suggesting that transcription could elicit their instability; however, the underlying mechanisms remain elusive. Genome-wide replication timing analyses here show that stress-induced delayed/under-replication is the hallmark of CFSs. Extensive genome-wide analyses of nascent transcripts, replication origin positioning and fork directionality reveal that 80% of CFSs nest in large transcribed domains poor in initiation events, replicated by long-travelling forks. Forks that travel long in late S phase explains CFS replication features, whereas formation of sequence-dependent fork barriers or head-on transcription–replication conflicts do not. We further show that transcription inhibition during S phase, which suppresses transcription–replication encounters and prevents origin resetting, could not rescue CFS stability. Altogether, our results show that transcription-dependent suppression of initiation events delays replication of large gene bodies, committing them to instability.

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 G1- 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.

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.

PREP1 tumor suppressor protects the late-replicating DNA by controlling its replication timing and symmetry

The synthesis of middle-to-late-replicating DNA can be affected independently of the rest of the genome by down-regulating the tumor suppressor PREP1 (PKNOX1). Indeed, DNA combing shows that PREP1 down-regulation affects DNA replication rate, increases the number of simultaneously firing origins and the asymmetry of DNA replication, leading to DNA damage. Genome-wide analysis of replication timing by Repli-seq shows that, upon PREP1 down-regulation, 25% of the genome is replicated earlier in the S-phase. The targeted DNA sequences correspond to Lamin-Associated Domains (LADs), and include late-replicating (LRRs) and temporal transition regions (TTRs). Notably, the distribution of PREP1 DNA binding sites and of its target genes indicates that DNA replication defects are independent of the overall PREP1 transcriptional activity. Finally, PREP1 down-regulation causes a substantial decrease in Lamin B1 levels. This suggests that DNA is released from the nuclear lamina earlier than in the control cells and is available for replication, thus explaining timing defects and DNA damage.This is the first evidence that the replication timing of a specific fraction of the human genome is affected by PREP1 tumor suppressor. This previously unknown function might significantly contribute to the genomic instability observed in human tumors.

Analysis of Single-Locus Replication Timing in Asynchronous Cycling Cells.

In higher eucaryotes, not all the genome is replicated simultaneously: there are parts of the genome that replicate at the beginning of S-phase (early S-phase), others that are replicated later. In each cell, early replicating genomic regions are alternated to late-replicating regions. In general, eucaryotic genomes are organized into structural domains where genes showing the same epigenetic state replicate at the same time.Here, we will describe the protocol that we routinely used for the analysis of replication timing of specific loci in Drosophila embryonic cell lines (S2 and S3) based on BrdU labeling and FACS sorting of different S-phase fractions (early, mid, late) of asynchronous cycling cells.

At short telomeres Tel1 directs early replication and phosphorylates Rif1.

The replication time of Saccharomyces cerevisiae telomeres responds to TG1–3 repeat length, with telomeres of normal length replicating late during S phase and short telomeres replicating early. Here we show that Tel1 kinase, which is recruited to short telomeres, specifies their early replication, because we find a tel1Δ mutant has short telomeres that nonetheless replicate late. Consistent with a role for Tel1 in driving early telomere replication, initiation at a replication origin close to an induced short telomere was reduced in tel1Δ cells, in an S phase blocked by hydroxyurea. The telomeric chromatin component Rif1 mediates late replication of normal telomeres and is a potential substrate of Tel1 phosphorylation, so we tested whether Tel1 directs early replication of short telomeres by inactivating Rif1. A strain lacking both Rif1 and Tel1 behaves like a rif1Δ mutant by replicating its telomeres early, implying that Tel1 can counteract the delaying effect of Rif1 to control telomere replication time. Proteomic analyses reveals that in yku70Δ cells that have short telomeres, Rif1 is phosphorylated at Tel1 consensus sequences (S/TQ sites), with phosphorylation of Serine-1308 being completely dependent on Tel1. Replication timing analysis of a strain mutated at these phosphorylation sites, however, suggested that Tel1-mediated phosphorylation of Rif1 is not the sole mechanism of replication timing control at telomeres. Overall, our results reveal two new functions of Tel1 at shortened telomeres: phosphorylation of Rif1, and specification of early replication by counteracting the Rif1-mediated delay in initiation at nearby replication origins.

Analysis of replication timing using synchronized budding yeast cultures.

Eukaryotic DNA replication exhibits at once extraordinary fidelity and substantial plasticity. The importance of the apparent presence of a replication temporal program on a population level has been the subject of intense debate of late. Such debate has been, to a great extent, facilitated by methods that permit the description and analysis of replication dynamics in various model organisms, both globally and at a single-molecule level. Each of these methods provides a unique view of the replication process, and also presents challenges and questions in the interpretation of experimental observations. Thus, wider applications of these methods in different genetic backgrounds and in different organisms would doubtless enable us to better understand the execution and regulation of chromosomal DNA synthesis as well as its impact on genome maintenance.

PcG-Mediated Higher-Order Chromatin Structures Modulate Replication Programs at the Drosophila BX-C

Polycomb group proteins (PcG) exert conserved epigenetic functions that convey maintenance of repressed transcriptional states, via post-translational histone modifications and high order structure formation. During S-phase, in order to preserve cell identity, in addition to DNA information, PcG-chromatin-mediated epigenetic signatures need to be duplicated requiring a tight coordination between PcG proteins and replication programs. However, the interconnection between replication timing control and PcG functions remains unknown. Using Drosophila embryonic cell lines, we find that, while presence of specific PcG complexes and underlying transcription state are not the sole determinants of cellular replication timing, PcG-mediated higher-order structures appear to dictate the timing of replication and maintenance of the silenced state. Using published datasets we show that PRC1, PRC2, and PhoRC complexes differently correlate with replication timing of their targets. In the fully repressed BX-C, loss of function experiments revealed a synergistic role for PcG proteins in the maintenance of replication programs through the mediation of higher-order structures. Accordingly, replication timing analysis performed on two Drosophila cell lines differing for BX-C gene expression states, PcG distribution, and chromatin domain conformation revealed a cell-type-specific replication program that mirrors lineage-specific BX-C higher-order structures. Our work suggests that PcG complexes, by regulating higher-order chromatin structure at their target sites, contribute to the definition and the maintenance of genomic structural domains where genes showing the same epigenetic state replicate at the same time.

Systematic Determination of Replication Activity Type Highlights Interconnections between Replication, Chromatin Structure and Nuclear Localization

DNA replication is a highly regulated process, with each genomic locus replicating at a distinct time of replication (ToR). Advances in ToR measurement technology enabled several genome-wide profiling studies that revealed tight associations between ToR and general genomic features and a remarkable ToR conservation in mammals. Genome wide studies further showed that at the hundreds kb-to-megabase scale the genome can be divided into constant ToR regions (CTRs) in which the replication process propagates at a faster pace due to the activation of multiple origins and temporal transition regions (TTRs) in which the replication process propagates at a slower pace. We developed a computational tool that assigns a ToR to every measured locus and determines its replication activity type (CTR versus TTR). Our algorithm, ARTO (Analysis of Replication Timing and Organization), uses signal processing methods to fit a constant piece-wise linear curve to the measured raw data. We tested our algorithm and provide performance and usability results. A Matlab implementation of ARTO is available at http://bioinfo.cs.technion.ac.il/people/zohar/ARTO/. Applying our algorithm to ToR data measured in multiple mouse and human samples allowed precise genome-wide ToR determination and replication activity type characterization. Analysis of the results highlighted the plasticity of the replication program. For example, we observed significant ToR differences in 10–25% of the genome when comparing different tissue types. Our analyses also provide evidence for activity type differences in up to 30% of the probes. Integration of the ToR data with multiple aspects of chromosome organization characteristics suggests that ToR plays a role in shaping the regional chromatin structure. Namely, repressive chromatin marks, are associated with late ToR both in TTRs and CTRs. Finally, characterization of the differences between TTRs and CTRs, with matching ToR, revealed that TTRs are associated with compact chromatin and are located significantly closer to the nuclear envelope. Supplementary material is available. Raw and processed data were deposited in Geo (GSE17236).

Genome-wide analysis of replication timing in mammalian cells: Troubleshooting problems encountered when comparing different cell types

DNA is replicated in a defined temporal order that is developmentally regulated and constitutes a unique and stable fingerprint of a given cell type. Recently, we developed a robust assay to profile replication timing genome wide that can be applied to essentially any proliferating cell population. Asynchronously cycling cells are pulse labeled with the nucleotide analog 5-bromo-2-deoxyuridine (BrdU). The cells are sorted into S-phase fractions on the basis of DNA content using flow cytometry. BrdU-labeled DNA from each fraction is immunoprecipitated (BrdU IP), amplified, differentially labeled and co-hybridized to a whole-genome comparative genomic hybridization microarray (or sequenced). Since the basic steps of this protocol have been detailed elsewhere, here we focus on problems encountered when adapting this protocol to different cell types or tissue sources and modifications that have been successfully applied to troubleshoot these problems. There is an increasing demand for such studies to address how replication is regulated during development, its relationship to chromatin architecture and other chromosome functions, and the relevance of cell culture models to regulation in the native organismal niche.