Replication Timing Patterns Research Articles

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.

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.

Collective regulation of chromatin modifications predicts replication timing during cell cycle.

SUMMARYReplication timing (RT) associates with genome architecture, while having a mixed relationship to histone marks. By profiling replication at high resolution and assessing broad histone marks across the cell cycle at the resolution of RT with and without genetic perturbation, we address the causal relationship between histone marks and RT. Four primary chromatin states, including an uncharacterized H3K36me2 state, emerge and define 97% of the mappable genome. RT and local replication patterns (e.g., initiation zones) quantitatively associate with chromatin states, histone mark dynamics, and spatial chromatin structure. Manipulation of broad histone marks and enhancer elements by overexpressing the histone H3 lysine 9/36 tri-demethylase KDM4A impacts RT across 11% of the genome. Broad histone modification changes were strong predictors of the observed RT alterations. Lastly, replication within H3K36me2-enriched neighborhoods is sensitive to KDM4A overexpression and is controlled at a megabase scale. These studies establish a role for collective chromatin mark regulation in modulating RT.

Promoters are key organizers of the duplication of vertebrate genomes.

In vertebrates, single cell analyses of replication timing patterns brought to light a very well controlled program suggesting a tight regulation on initiation sites. Mapping of replication origins with different methods has revealed discrete preferential sites, enriched in promoters and potential G-quadruplex motifs, which can aggregate into initiation zones spanning several tens of kilobases (kb). Another characteristic of replication origins is a nucleosome-free region (NFR). A modified yeast strain containing a humanized origin recognition complex (ORC) fires new origins at NFRs revealing their regulatory role. In cooperation with NFRs, the histone variant H2A.Z facilitates ORC loading through di-methylation of lysine 20 of histone H4. Recent studies using genome editing methods show that efficient initiation sites associated with transcriptional activity can synergize over several tens of kb by establishing physical contacts and lead to the formation of early domains of DNA replication demonstrating a co-regulation between replication initiation and transcription.

The capacity of origins to load MCM establishes replication timing patterns.

Loading of the MCM replicative helicase at origins of replication is a highly regulated process that precedes DNA replication in all eukaryotes. The stoichiometry of MCM loaded at origins has been proposed to be a key determinant of when those origins initiate replication during S phase. Nevertheless, the genome-wide regulation of MCM loading stoichiometry and its direct effect on replication timing remain unclear. In order to investigate why some origins load more MCM than others, we perturbed MCM levels in budding yeast cells and, for the first time, directly measured MCM levels and replication timing in the same experiment. Reduction of MCM levels through degradation of Mcm4, one of the six obligate components of the MCM complex, slowed progression through S phase and increased sensitivity to replication stress. Reduction of MCM levels also led to differential loading at origins during G1, revealing origins that are sensitive to reductions in MCM and others that are not. Sensitive origins loaded less MCM under normal conditions and correlated with a weak ability to recruit the origin recognition complex (ORC). Moreover, reduction of MCM loading at specific origins of replication led to a delay in their replication during S phase. In contrast, overexpression of MCM had no effects on cell cycle progression, relative MCM levels at origins, or replication timing, suggesting that, under optimal growth conditions, cellular MCM levels are not limiting for MCM loading. Our results support a model in which the loading capacity of origins is the primary determinant of MCM stoichiometry in wild-type cells, but that stoichiometry is controlled by origins' ability to recruit ORC and compete for MCM when MCM becomes limiting.

Characteristics of mutational signatures of unknown etiology.

Although not all somatic mutations are cancer drivers, their mutational signatures, i.e. the patterns of genomic alterations at a genome-wide scale, provide insights into past exposure to mutagens, DNA damage and repair processes. Computational deconvolution of somatic mutation patterns and expert curation pan-cancer studies have identified a number of mutational signatures associated with point mutations, dinucleotide substitutions, insertions and deletions, and rearrangements, and have established etiologies for a subset of these signatures. However, the mechanisms underlying nearly one-third of all mutational signatures are not yet understood. The signatures with established etiology and those with hitherto unknown origin appear to have some differences in strand bias, GC content and nucleotide context diversity. It is possible that some of the hitherto ‘unknown’ signatures predominantly occur outside gene regions. While nucleotide contexts might be adequate to establish etiologies of some mutational signatures, in other cases additional features, such as broader (epi)genomic contexts, including chromatin, replication timing, processivity and local mutational patterns, may help fully understand the underlying DNA damage and repair processes. Nonetheless, remarkable progress in characterization of mutational signatures has provided fundamental insights into the biology of cancer, informed disease etiology and opened up new opportunities for cancer prevention, risk management, and therapeutic decision making.

Local rewiring of genome–nuclear lamina interactions by transcription

Transcriptionally inactive genes are often positioned at the nuclear lamina (NL), as part of large lamina‐associated domains (LADs). Activation of such genes is often accompanied by repositioning toward the nuclear interior. How this process works and how it impacts flanking chromosomal regions are poorly understood. We addressed these questions by systematic activation or inactivation of individual genes, followed by detailed genome‐wide analysis of NL interactions, replication timing, and transcription patterns. Gene activation inside LADs typically causes NL detachment of the entire transcription unit, but rarely more than 50–100 kb of flanking DNA, even when multiple neighboring genes are activated. The degree of detachment depends on the expression level and the length of the activated gene. Loss of NL interactions coincides with a switch from late to early replication timing, but the latter can involve longer stretches of DNA. Inactivation of active genes can lead to increased NL contacts. These extensive datasets are a resource for the analysis of LAD rewiring by transcription and reveal a remarkable flexibility of interphase chromosomes.

Abstract 4611: Identifying drivers of mammary tumorigenesis & elucidating the mechanisms of cancer initiation in DNA replication defective Chaos3 mice

Abstract Distinguishing commonly deleted genes as drivers or passengers of human cancers are key to delineating mechanisms of tumor initiation. The “Chaos3” mouse model of spontaneous breast cancer carries a missense mutation in a gene called Mcm4 (minichromosome maintenance 4), causing destabilization of the MCM2-7 replicative helicase, which in turn causes elevated genomic instability. Approximately 80% of female mice homozygous for this mutation in the C3HeB/FeJ strain background develop mammary tumors with an average latency of 12 months. Genomic analysis of these tumors revealed a recurrent subset of genes with copy number alterations such as Arid1a and Nf1. This project aims to: 1) test whether the genes Arid1a and Nf1 act as mammary tumor drivers in the Chaos3 breast cancer model; 2) understand the genome-wide effects of Arid1a loss in mammary tumorigenesis; 3) identifying potential therapeutic targets for Arid1a-deficient mammary tumors; and 4) unravel the mechanism by which destabilization of MCM2-7 replicative helicase leads to cancer onset. To achieve these aims, several experiments such as conditional mutagenesis, CUTandRUN, RNA-seq, and replication timing profiling were conducted to address my aims. Preliminary data indicates that heterozygous loss of Arid1ain Chaos3 mammary tumors produces a distinct transcriptional profile compared to mammary tumors without Arid1a deletion. Additionally, Chaos3 primary cells have markedly altered replication timing patterns in certain regions of the genome that might explain increased susceptibility to mutations. My experiments can validate new candidates like Arid1a as a driver of sporadic breast cancer in humans. In addition, these studies will not only lend insight into how key tumor suppressors are deleted in the Chaos3 model, but also potentially reveal features of genomic regions that are particularly susceptible to genetic or environmental perturbations to DNA replication. The Chaos3 sporadic breast cancer mouse model mimics sporadic human luminal breast cancer and can be used as a tool to better understand the biology and genetics of a common disease that plagues women worldwide. Citation Format: Marquita Winters, John Schimenti. Identifying drivers of mammary tumorigenesis & elucidating the mechanisms of cancer initiation in DNA replication defective Chaos3 mice [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4611.

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.

Profiling DNA Replication Timing Using Zebrafish as an In Vivo Model System.

DNA replication timing is an important cellular characteristic, exhibiting significant relationships with chromatin structure, transcription, and DNA mutation rates. Changes in replication timing occur during development and in cancer, but the role replication timing plays in development and disease is not known. Zebrafish were recently established as an in vivo model system to study replication timing. Here is detailed the protocols for using the zebrafish to determine DNA replication timing. After sorting cells from embryos and adult zebrafish, high-resolution genome-wide DNA replication timing patterns can be constructed by determining changes in DNA copy number through analysis of next generation sequencing data. The zebrafish model system allows for evaluation of the replication timing changes that occur in vivo throughout development, and can also be used to assess changes in individual cell types, disease models, or mutant lines. These methods will enable studies investigating the mechanisms and determinants of replication timing establishment and maintenance during development, the role replication timing plays in mutations and tumorigenesis, and the effects of perturbing replication timing on development and disease.

Bacterial artificial chromosomes establish replication timing and sub-nuclear compartment de novo as extra-chromosomal vectors.

The role of DNA sequence in determining replication timing (RT) and chromatin higher order organization remains elusive. To address this question, we have developed an extra-chromosomal replication system (E-BACs) consisting of ∼200 kb human bacterial artificial chromosomes (BACs) modified with Epstein-Barr virus (EBV) stable segregation elements. E-BACs were stably maintained as autonomous mini-chromosomes in EBNA1-expressing HeLa or human induced pluripotent stem cells (hiPSCs) and established distinct RT patterns. An E-BAC harboring an early replicating chromosomal region replicated early during S phase, while E-BACs derived from RT transition regions (TTRs) and late replicating regions replicated in mid to late S phase. Analysis of E-BAC interactions with cellular chromatin (4C-seq) revealed that the early replicating E-BAC interacted broadly throughout the genome and preferentially with the early replicating compartment of the nucleus. In contrast, mid- to late-replicating E-BACs interacted with more specific late replicating chromosomal segments, some of which were shared between different E-BACs. Together, we describe a versatile system in which to study the structure and function of chromosomal segments that are stably maintained separately from the influence of cellular chromosome context.

A stitch in time: Replicate early and escape dosage compensation to express more.

The biological significance of conserved replication timing patterns in eukaryotic genomes remains a mystery. In this issue, Müller and Nieduszynski (2017. J. Cell Biol. https://doi.org/10.1083/jcb.201701061) find that early replication is a requirement for the highest expression levels of certain genes.

Genetic Variation Meets Replication Origins

Genome replication programs are highly orchestrated. In this issue, Koren and colleagues leverage whole-genome sequencing data to discover that DNA replication timing patterns differ between individuals and are associated with genetic variants. These so-called rtQTLs represent a new form of quantitative trait loci that may influence gene dosage and mutational frequency and provide mechanistic insights into the regulation of DNA replication.

Cell line differences in replication timing of human glutamate receptor genes and other large genes associated with neural disease

There is considerable current interest in the function of epigenetic mechanisms in neuroplasticity with regard to learning and memory formation and to a range of neural diseases. Previously, we described replication timing on human chromosome 21q in the THP-1 human cell line (2n = 46, XY) and showed that several genes associated with neural diseases, such as the neuronal glutamate receptor subunit GluR-5 (GRIK1) and amyloid precursor protein (APP), were located in regions where replication timing transitioned from early to late S phase. Here, we compared replication timing of all known human glutamate receptor genes (26 genes in total) and APP in 6 different human cell lines including human neuron-related cell lines. Replication timings were obtained by integrating our previously reported data with new data generated here and information from the online database ReplicationDomain. We found that many of the glutamate receptor genes were clearly located in replication timing transition zones in neural precursor cells, but this relationship was less clear in embryonic stem cells before neural differentiation; in the latter, the genes were often located in later replication timing zones that displayed DNA hypermethylation. Analysis of selected large glutamate receptor genes (>200 kb), and of APP, showed that their precise replication timing patterns differed among the cell lines. We propose that the transition zones of DNA replication timing are altered by epigenetic mechanisms, and that these changes may affect the neuroplasticity that is important to memory and learning, and may also have a role in the development of neural diseases.

A mechanistic role for DNA methylation in endothelial cell (EC)-enriched gene expression: relationship with DNA replication timing

AbstractProximal promoter DNA methylation has been shown to be important for regulating gene expression. However, its relative contribution to the cell-specific expression of endothelial cell (EC)-enriched genes has not been defined. We used methyl-DNA immunoprecipitation and bisulfite conversion to analyze the DNA methylation profile of EC-enriched genes in ECs vs nonexpressing cell types, both in vitro and in vivo. We show that prototypic EC-enriched genes exhibit functional differential patterns of DNA methylation in proximal promoter regions of most (eg, CD31, von Willebrand factor [vWF], VE-cadherin, and intercellular adhesion molecule-2), but not all (eg, VEGFR-1 and VEGFR-2), EC-enriched genes. Comparable findings were evident in cultured ECs, human blood origin ECs, and murine aortic ECs. Promoter-reporter episomal transfection assays for endothelial nitric oxide synthase, VE-cadherin, and vWF indicated functional promoter activity in cell types where the native gene was not active. Inhibition of DNA methyltransferase activity indicated important functional relevance. Importantly, profiling DNA replication timing patterns indicated that EC-enriched gene promoters with differentially methylated regions replicate early in S-phase in both expressing and nonexpressing cell types. Collectively, these studies highlight the functional importance of promoter DNA methylation in controlling vascular EC gene expression.

Abstract 5150: What is in control of DNA replication timing.

Abstract The initiation of DNA replication at specific sites in the genome is known to be far from perfect, where known initiation sites may or may not fire in a given single cell during specific cell cycle. At the same time, from larger scale perspective, one finds that the timing of DNA replication is highly orchestrated, where specific regions in the genome are copied at specific, reproducible time points throughout S-phase. In an attempt to reconcile these two somewhat discordant observations, many researchers have assumed that a yet unknown, active regulatory mechanism in the cell is responsible for coordinating replication timing. Here we present a mechanistic explanation for genome-wide replication timing patterns. Our model is purely stochastic, does not contain any active control components and reproduces existing replication timing data with extreme fidelity (correlation > 0.9, about the same as the correlation between technical replicates of the experiment). The main findings of our work are • Existence of reproducible genome-wide replication patterns do not logically necessitate existence of a complex, active control system in the cell • Integrative analysis of recently released ENCODE data identifies semi-local chromatin structure as the main determinant of replication timing • Plasticity of replication timing across different cell types are a direct consequence of changes in the chromatin structure between the different cell differentiation states • Replication timing is a “systems” phenomenon emerging from the collective action of many uncoordinated, stochastic replication initiation sites. Replication timing cannot be understood by investigating individual initiation sites: deleting up to 30% of high probability initiation sites does not significantly change replication timing Citation Format: Yevgeniy Gindin, Paul S. Meltzer, Sven Bilke. What is in control of DNA replication timing. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 5150. doi:10.1158/1538-7445.AM2013-5150

Clonal allelic predetermination of immunoglobulin-κ rearrangement

Although most genes are expressed biallelically, a number of key genomic sites--including immune and olfactory receptor regions--are controlled monoallelically in a stochastic manner, with some cells expressing the maternal allele and others the paternal allele in the target tissue. Very little is known about how this phenomenon is regulated and programmed during development. Here, using mouse immunoglobulin-κ (Igκ) as a model system, we demonstrate that although individual haematopoietic stem cells are characterized by allelic plasticity, early lymphoid lineage cells become committed to the choice of a single allele, and this decision is then stably maintained in a clonal manner that predetermines monoallelic rearrangement in B cells. This is accompanied at the molecular level by underlying allelic changes in asynchronous replication timing patterns at the κ locus. These experiments may serve to define a new concept of stem cell plasticity.

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

DNA replication in eukaryotes initiates from discrete genomic regions, termed origins, according to a strict and often tissue-specific temporal program. However, the genetic program that controls activation of replication origins has still not been fully elucidated in mammalian cells. Previously, we measured replication timing at the sequence level along human chromosomes 11q and 21q. In the present study, we sought to obtain a greater understanding of the relationship between replication timing programs and human chromosomes by analysis of the timing of replication of a single human chromosome 11 that had been transferred into the Chinese hamster ovary (CHO) cell line by chromosome engineering. Timing of replication was compared for three 11q chromosomal regions in the transformed CHO cell line (CHO(h11)) and the original human fibroblast cell line, namely, the R/G-band boundary at 11q13.5/q14.1, the centromere and the distal telomere. We found that the pattern of replication timing in and around the R/G band boundary at 11q13.5/q14.1 was similar in CHO(h11) cells and fibroblasts. The 11q centromeric region, which replicates late in human fibroblasts, replicated in the second half of S phase in CHO(h11) cells. By contrast, however, the telomeric region at 11q25, which is late replicating in fibroblasts (and in several other human cell lines), replicated in the first half of S phase or in very early S phase in CHO(h11) cells. Our observations suggest that the replication timing programs of the R/G-band boundary and the centromeric region of human chromosome 11q are maintained in CHO(h11) cells, whereas that for the telomeric region is altered. The replication timing program of telomeric regions on human chromosomes might be regulated by specific mechanisms that differ from those for other chromosomal regions.

Reprogramming of DNA Replication Timing

Replication timing is an important developmentally regulated regional property that is correlated with chromosome structure and gene expression, but little is known about the establishment and maintenance of these patterns. Here we followed the fate of replication timing patterns in cells that undergo reprogramming either through somatic-cell nuclear transplantation or by the generation of induced pluripotential stem cells. We have investigated three different paradigms, stage-specific replication timing, parental allele-specific asynchrony (imprinted regions), and random allelic asynchronous replication. In all cases, somatic replication timing patterns were reset exactly at the appropriate stage in early development and could be properly established upon re-differentiation. Taken together, these results suggest that, unlike DNA methylation, the molecular mechanisms governing replication timing are not only stable but can also be easily reprogrammed.

Sequencing newly replicated DNA reveals widespread plasticity in human replication timing

Faithful transmission of genetic material to daughter cells involves a characteristic temporal order of DNA replication, which may play a significant role in the inheritance of epigenetic states. We developed a genome-scale approach--Repli Seq--to map temporally ordered replicating DNA using massively parallel sequencing and applied it to study regional variation in human DNA replication time across multiple human cell types. The method requires as few as 8,000 cytometry-fractionated cells for a single analysis, and provides high-resolution DNA replication patterns with respect to both cell-cycle time and genomic position. We find that different cell types exhibit characteristic replication signatures that reveal striking plasticity in regional replication time patterns covering at least 50% of the human genome. We also identified autosomal regions with marked biphasic replication timing that include known regions of monoallelic expression as well as many previously uncharacterized domains. Comparison with high-resolution genome-wide profiles of DNaseI sensitivity revealed that DNA replication typically initiates within foci of accessible chromatin comprising clustered DNaseI hypersensitive sites, and that replication time is better correlated with chromatin accessibility than with gene expression. The data collectively provide a unique, genome-wide picture of the epigenetic compartmentalization of the human genome and suggest that cell-lineage specification involves extensive reprogramming of replication timing patterns.