Replication Timing Research Articles

Reducing the replication time for structural estimations: A successful replication of “An Anatomy of International Trade” using GPU computing

AbstractEaton, Kortum, and Kramarz (2011) (EKK) discovered empirical patterns from French manufacturing firms that a baseline firm heterogeneity model could not explain. The authors proposed and estimated a model that closely matches the patterns observed in French data. This paper successfully replicates their findings using author‐provided data, re‐implementing their algorithms in Python and leveraging graphics processing unit computing to significantly boost computational speed. Applying the model to Chinese manufacturing data, despite differences in context, showed that the model explains most observed patterns well.

Inferring replication timing and proliferation dynamics from single-cell DNA sequencing data

Dysregulated DNA replication is a cause and a consequence of aneuploidy in cancer, yet the interplay between copy number alterations (CNAs), replication timing (RT) and cell cycle dynamics remain understudied in aneuploid tumors. We developed a probabilistic method, PERT, for simultaneous inference of cell-specific replication and copy number states from single-cell whole genome sequencing (scWGS) data. We used PERT to investigate clone-specific RT and proliferation dynamics in >50,000 cells obtained from aneuploid and clonally heterogeneous cell lines, xenografts and primary cancers. We observed bidirectional relationships between RT and CNAs, with CNAs affecting X-inactivation producing the largest RT shifts. Additionally, we found that clone-specific S-phase enrichment positively correlated with ground-truth proliferation rates in genomically stable but not unstable cells. Together, these results demonstrate robust computational identification of S-phase cells from scWGS data, and highlight the importance of RT and cell cycle properties in studying the genomic evolution of aneuploid tumors.

Mathematical model for the distribution of DNA replication origins.

DNA replication in yeast and in many other organisms starts from well-defined locations on the DNA known as replication origins. The spatial distribution of these origins in the genome is particularly important in ensuring that replication is completed quickly. Cells are more vulnerable to DNA damage and other forms of stress while they are replicating their genome. This raises the possibility that the spatial distribution of origins is under selection pressure. In this paper we investigate the hypothesis that natural selection favors origin distributions leading to shorter replication times. Using a simple mathematical model, we show that this hypothesis leads to two main predictions about the origin distributions: that neighboring origins that are inefficient (less likely to fire) are more likely to be close to each other than efficient origins; and that neighboring origins with larger differences in firing times are more likely to be close to each other than origins with similar firing times. We test these predictions using next-generation sequencing data, and show that they are both supported by the data.

Nucleolus and centromere Tyramide Signal Amplification-Seq reveals variable localization of heterochromatin in different cell types

Genome differential positioning within interphase nuclei remains poorly explored. We extended and validated Tyramide Signal Amplification (TSA)-seq to map genomic regions near nucleoli and pericentric heterochromatin in four human cell lines. Our study confirmed that smaller chromosomes localize closer to nucleoli but further deconvolved this by revealing a preference for chromosome arms below 36-46 Mbp in length. We identified two lamina associated domain subsets through their differential nuclear lamina versus nucleolar positioning in different cell lines which showed distinctive patterns of DNA replication timing and gene expression across all cell lines. Unexpectedly, active, nuclear speckle-associated genomic regions were found near typically repressive nuclear compartments, which is attributable to the close proximity of nuclear speckles and nucleoli in some cell types, and association of centromeres with nuclear speckles in human embryonic stem cells (hESCs). Our study points to a more complex and variable nuclear genome organization than suggested by current models, as revealed by our TSA-seq methodology.

Centromeric transposable elements and epigenetic status drive karyotypic variation in the eastern hoolock gibbon.

Great apes have maintained a stable karyotype with few large-scale rearrangements; in contrast, gibbons have undergone a high rate of chromosomal rearrangements coincident with rapid centromere turnover. Here we characterize assembled centromeres in the Eastern hoolock gibbon, Hoolock leuconedys (HLE), finding a diverse group of transposable elements (TEs) that differ from the canonical alpha satellites found across centromeres of other apes. We find that HLE centromeres contain a CpG methylation centromere dip region, providing evidence this epigenetic feature is conserved in the absence of satellite arrays; nevertheless, we report a variety of atypical centromeric features, including protein-coding genes and mismatched replication timing. Further, large structural variations define HLE centromeres and distinguish them from other gibbons. Combined with differentially methylated TEs, topologically associated domain boundaries, and segmental duplications at chromosomal breakpoints, we propose that a "perfect storm" of multiple genomic attributes with propensities for chromosome instability shaped gibbon centromere evolution.

Embryonic genome instability upon DNA replication timing program emergence

Faithful DNA replication is essential for genome integrity1–4. Under-replicated DNA leads to defects in chromosome segregation, which are common during embryogenesis5–8. However, the regulation of DNA replication remains poorly understood in early mammalian embryos. Here we constructed a single-cell genome-wide DNA replication atlas of pre-implantation mouse embryos and identified an abrupt replication program switch accompanied by a transient period of genomic instability. In 1- and 2-cell embryos, we observed the complete absence of a replication timing program, and the entire genome replicated gradually and uniformly using extremely slow-moving replication forks. In 4-cell embryos, a somatic-cell-like replication timing program commenced abruptly. However, the fork speed was still slow, S phase was extended, and markers of replication stress, DNA damage and repair increased. This was followed by an increase in break-type chromosome segregation errors specifically during the 4-to-8-cell division with breakpoints enriched in late-replicating regions. These errors were rescued by nucleoside supplementation, which accelerated fork speed and reduced the replication stress. By the 8-cell stage, forks gained speed, S phase was no longer extended and chromosome aberrations decreased. Thus, a transient period of genomic instability exists during normal mouse development, preceded by an S phase lacking coordination between replisome-level regulation and megabase-scale replication timing regulation, implicating a link between their coordination and genome stability.

Sir2 and Fun30 regulate ribosomal DNA replication timing via MCM helicase positioning and nucleosome occupancy.

The association between late replication timing and low transcription rates in eukaryotic heterochromatin is well-known, yet the specific mechanisms underlying this link remain uncertain. In Saccharomyces cerevisiae, the histone deacetylase Sir2 is required for both transcriptional silencing and late replication at the repetitive ribosomal DNA arrays (rDNA). We have previously reported that in the absence of SIR2, a derepressed RNA PolII repositions MCM replicative helicases from their loading site at the ribosomal origin, where they abut well-positioned, high-occupancy nucleosomes, to an adjacent region with lower nucleosome occupancy. By developing a method that can distinguish activation of closely spaced MCM complexes, here we show that the displaced MCMs at rDNA origins have increased firing propensity compared to the nondisplaced MCMs. Furthermore, we found that both, activation of the repositioned MCMs and low occupancy of the adjacent nucleosomes critically depend on the chromatin remodeling activity of FUN30. Our study elucidates the mechanism by which Sir2 delays replication timing, and it demonstrates, for the first time, that activation of a specific replication origin in vivo relies on the nucleosome context shaped by a single chromatin remodeler.

ERCC2 mutations alter the genomic distribution pattern of somatic mutations and are independently prognostic in bladder cancer

Excision repair cross-complementation group 2 (ERCC2) encodes the DNA helicase xeroderma pigmentosum group D, which functions in transcription and nucleotide excision repair. Point mutations in ERCC2 are putative drivers in around 10% of bladder cancers (BLCAs) and a potential positive biomarker for cisplatin therapy response. Nevertheless, the prognostic significance directly attributed to ERCC2 mutations and its pathogenic role in genome instability remain poorly understood. We first demonstrated that mutant ERCC2 is an independent predictor of prognosis in BLCA. We then examined its impact on the somatic mutational landscape using a cohort of ERCC2 wild-type (n= 343) and mutant (n= 39) BLCA whole genomes. The genome-wide distribution of somatic mutations is significantly altered in ERCC2 mutants, including T[C>T]N enrichment, altered replication time correlations, and CTCF-cohesin binding site mutation hotspots. We leverage these alterations to develop a machine learning model for predicting pathogenic ERCC2 mutations, which may be useful to inform treatment of patients with BLCA.

Sex-specific DNA-replication in the early mammalian embryo

The timing of DNA replication in mammals is crucial for minimizing errors and influenced by genome usage and chromatin states. Replication timing in the newly formed mammalian embryo remains poorly understood. Here, we have investigated replication timing in mouse zygotes and 2-cell embryos, revealing that zygotes lack a conventional replication timing program, which then emerges in 2-cell embryos. This program differs from embryonic stem cells and generally correlates with transcription and genome compartmentalization of both parental genomes. However, consistent and systematic differences existed between the replication timing of the two parental genomes, including considerably later replication of maternal pericentromeric regions compared to paternal counterparts. Moreover, maternal chromatin modified by Polycomb Repressive Complexes in the oocyte, undergoes early replication, despite belonging to the typically late-replicating B-compartment of the genome. This atypical and asynchronous replication of the two parental genomes may advance our understanding of replication stress in early human embryos and trigger strategies to reduce errors and aneuploidies.

Resonant noise amplification in a predator-prey model with quasi-discrete generations

Predator-prey models have been shown to exhibit resonance-like behaviour, in which random fluctuations in the number of organisms (demographic noise) are amplified when their frequency is close to the natural oscillatory frequency of the system. This behaviour has been traditionally studied in models with exponentially distributed replication and death times. Here we consider a biologically more realistic model, in which organisms replicate quasi-synchronously such that the distribution of replication times has a narrow maximum at some T>0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$T>0$$\\end{document} corresponding to the mean doubling time. We show that when the frequency of replication f=1/T\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$f=1/T$$\\end{document} is tuned to the natural oscillatory frequency of the predator-prey model, the system exhibits oscillations that are much stronger than in the model with Poissonian (non-synchronous) replication and death. These oscillations lead to population instability and the extinction of one of the species much sooner than in the case of Poissonian replication. The effect can be explained by resonant amplification of coloured noise generated by quasi-synchronous replication events.

Replication timing alterations are associated with mutation acquisition during breast and lung cancer evolution

During each cell cycle, the process of DNA replication timing is tightly regulated to ensure the accurate duplication of the genome. The extent and significance of alterations in this process during malignant transformation have not been extensively explored. Here, we assess the impact of altered replication timing (ART) on cancer evolution by analysing replication-timing sequencing of cancer and normal cell lines and 952 whole-genome sequenced lung and breast tumours. We find that 6%–18% of the cancer genome exhibits ART, with regions with a change from early to late replication displaying an increased mutation rate and distinct mutational signatures. Whereas regions changing from late to early replication contain genes with increased expression and present a preponderance of APOBEC3-mediated mutation clusters and associated driver mutations. We demonstrate that ART occurs relatively early during cancer evolution and that ART may have a stronger correlation with mutation acquisition than alterations in chromatin structure.

Beyond A and B Compartments: how major nuclear locales define nuclear genome organization and function.

Models of nuclear genome organization often propose a binary division into active versus inactive compartments, yet they overlook nuclear bodies. Here we integrated analysis of sequencing and image-based data to compare genome organization in four human cell types relative to three different nuclear locales: the nuclear lamina, nuclear speckles, and nucleoli. Whereas gene expression correlates mostly with nuclear speckle proximity, DNA replication timing correlates with proximity to multiple nuclear locales. Speckle attachment regions emerge as DNA replication initiation zones whose replication timing and gene composition vary with their attachment frequency. Most facultative LADs retain a partially repressed state as iLADs, despite their positioning in the nuclear interior. Knock out of two lamina proteins, Lamin A and LBR, causes a shift of H3K9me3-enriched LADs from lamina to nucleolus, and a reciprocal relocation of H3K27me3-enriched partially repressed iLADs from nucleolus to lamina. Thus, these partially repressed iLADs appear to compete with LADs for nuclear lamina attachment with consequences for replication timing. The nuclear organization in adherent cells is polarized with nuclear bodies and genomic regions segregating both radially and relative to the equatorial plane. Together, our results underscore the importance of considering genome organization relative to nuclear locales for a more complete understanding of the spatial and functional organization of the human genome.

DNA replication in early mammalian embryos is patterned, predisposing lamina-associated regions to fragility

DNA replication in differentiated cells follows a defined program, but when and how it is established during mammalian development is not known. Here we show using single-cell sequencing, that late replicating regions are established in association with the B compartment and the nuclear lamina from the first cell cycle after fertilization on both maternal and paternal genomes. Late replicating regions contain a relative paucity of active origins and few but long genes and low G/C content. In both bovine and mouse embryos, replication timing patterns are established prior to embryonic genome activation. Chromosome breaks, which form spontaneously in bovine embryos at sites concordant with human embryos, preferentially locate to late replicating regions. In mice, late replicating regions show enhanced fragility due to a sparsity of dormant origins that can be activated under conditions of replication stress. This pattern predisposes regions with long neuronal genes to fragility and genetic change prior to separation of soma and germ cell lineages. Our studies show that the segregation of early and late replicating regions is among the first layers of genome organization established after fertilization.

ASAR lncRNAs control DNA replication timing through interactions with multiple hnRNP/RNA binding proteins.

ASARs are a family of very-long noncoding RNAs that control replication timing on individual human autosomes, and are essential for chromosome stability. The eight known ASAR lncRNAs remain closely associated with their parent chromosomes. Analysis of RNA-protein interaction data (from ENCODE) revealed numerous RBPs with significant interactions with multiple ASAR lncRNAs, with several hnRNPs as abundant interactors. An ~7 kb domain within the ASAR6-141 lncRNA shows a striking density of RBP interaction sites. Genetic deletion and ectopic integration assays indicate that this ~7 kb RNA binding protein domain contains functional sequences for controlling replication timing of entire chromosomes in cis. shRNA-mediated depletion of 10 different RNA binding proteins, including HNRNPA1, HNRNPC, HNRNPL, HNRNPM, HNRNPU, or HNRNPUL1, results in dissociation of ASAR lncRNAs from their chromosome territories, and disrupts the synchronous replication that occurs on all autosome pairs, recapitulating the effect of individual ASAR knockouts on a genome-wide scale. Our results further demonstrate the role that ASARs play during the temporal order of genome-wide replication, and we propose that ASARs function as essential RNA scaffolds for the assembly of hnRNP complexes that help maintain the structural integrity of each mammalian chromosome.

ASAR lncRNAs control DNA replication timing through interactions with multiple hnRNP/RNA binding proteins

ASARs are a family of very-long noncoding RNAs that control replication timing on individual human autosomes, and are essential for chromosome stability. The eight known ASAR lncRNAs remain closely associated with their parent chromosomes. Analysis of RNA-protein interaction data (from ENCODE) revealed numerous RBPs with significant interactions with multiple ASAR lncRNAs, with several hnRNPs as abundant interactors. An ~7 kb domain within the ASAR6-141 lncRNA shows a striking density of RBP interaction sites. Genetic deletion and ectopic integration assays indicate that this ~7 kb RNA binding protein domain contains functional sequences for controlling replication timing of entire chromosomes in cis. shRNA-mediated depletion of 10 different RNA binding proteins, including HNRNPA1, HNRNPC, HNRNPL, HNRNPM, HNRNPU, or HNRNPUL1, results in dissociation of ASAR lncRNAs from their chromosome territories, and disrupts the synchronous replication that occurs on all autosome pairs, recapitulating the effect of individual ASAR knockouts on a genome-wide scale. Our results further demonstrate the role that ASARs play during the temporal order of genome-wide replication, and we propose that ASARs function as essential RNA scaffolds for the assembly of hnRNP complexes that help maintain the structural integrity of each mammalian chromosome.

DBF4, not DRF1, is the crucial regulator of CDC7 kinase at replication forks.

CDC7 kinase is crucial for DNA replication initiation and is involved in fork processing and replication stress response. Human CDC7 requires the binding of either DBF4 or DRF1 for its activity. However, it is unclear whether the two regulatory subunits target CDC7 to a specific set of substrates, thus having different biological functions, or if they act redundantly. Using genome editing technology, we generated isogenic cell lines deficient in either DBF4 or DRF1: these cells are viable but present signs of genomic instability, indicating that both can independently support CDC7 for bulk DNA replication. Nonetheless, DBF4-deficient cells show altered replication efficiency, partial deficiency in MCM helicase phosphorylation, and alterations in the replication timing of discrete genomic regions. Notably, we find that CDC7 function at replication forks is entirely dependent on DBF4 and not on DRF1. Thus, DBF4 is the primary regulator of CDC7 activity, mediating most of its functions in unperturbed DNA replication and upon replication interference.

PHF2 regulates genome topology and DNA replication in neural stem cells via cohesin.

Cohesin plays a crucial role in the organization of topologically-associated domains (TADs), which influence gene expression and DNA replication timing. Whether epigenetic regulators may affect TADs via cohesin to mediate DNA replication remains elusive. Here, we discover that the histone demethylase PHF2 associates with RAD21, a core subunit of cohesin, to regulate DNA replication in mouse neural stem cells (NSC). PHF2 loss impairs DNA replication due to the activation of dormant replication origins in NSC. Notably, the PHF2/RAD21 co-bound genomic regions are characterized by CTCF enrichment and epigenomic features that resemble efficient, active replication origins, and can act as boundaries to separate adjacent domains. Accordingly, PHF2 loss weakens TADs and chromatin loops at the co-bound loci due to reduced RAD21 occupancy. The observed topological and DNA replication defects in PHF2 KO NSC support a cohesin-dependent mechanism. Furthermore, we demonstrate that the PHF2/RAD21 complex exerts little effect on gene regulation, and that PHF2's histone-demethylase activity is dispensable for normal DNA replication and proliferation of NSC. We propose that PHF2 may serve as a topological accessory to cohesin for cohesin localization to TADs and chromatin loops, where cohesin represses dormant replication origins directly or indirectly, to sustain DNA replication in NSC.

The epidermal circadian clock integrates and subverts brain signals to guarantee skin homeostasis

In mammals, the circadian clock network drives daily rhythms of tissue-specific homeostasis. To dissect daily inter-tissue communication, we constructed a mouse minimal clock network comprising only two nodes: the peripheral epidermal clock and the central brain clock. By transcriptomic and functional characterization of this isolated connection, we identified a gatekeeping function of the peripheral tissue clock with respect to systemic inputs. The epidermal clock concurrently integrates and subverts brain signals to ensure timely execution of epidermal daily physiology. Timely cell-cycle termination in the epidermal stem cell compartment depends upon incorporation of clock-driven signals originating from the brain. In contrast, the epidermal clock corrects or outcompetes potentially disruptive feeding-related signals to ensure the optimal timing of DNA replication. Together, we present an approach for cataloging the systemic dependencies of daily temporal organization in a tissue and identify an essential gate-keeping function of peripheral circadian clocks that guarantees tissue homeostasis.

Rationally designed chromosome fusion does not prevent rapid growth of Vibrio natriegens

DNA replication is essential for the proliferation of all cells. Bacterial chromosomes are replicated bidirectionally from a single origin of replication, with replication proceeding at about 1000 bp per second. For the model organism, Escherichia coli, this translates into a replication time of about 40 min for its 4.6 Mb chromosome. Nevertheless, E. coli can propagate by overlapping replication cycles with a maximum short doubling time of 20 min. The fastest growing bacterium known, Vibrio natriegens, is able to replicate with a generation time of less than 10 min. It has a bipartite genome with chromosome sizes of 3.2 and 1.9 Mb. Is simultaneous replication from two origins a prerequisite for its rapid growth? We fused the two chromosomes of V. natriegens to create a strain carrying one chromosome with a single origin of replication. Compared to the parental, this strain showed no significant deviation in growth rate. This suggests that the split genome is not a prerequisite for rapid growth.

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.