Strand-resolved mutagenicity of DNA damage and repair.

DNA base damage is a major source of oncogenic mutations1. Such damage can produce strand-phased mutation patterns and multiallelic variation through the process of lesion segregation2. Here we exploited these properties to reveal how strand-asymmetric processes, such as replication and transcription, shape DNA damage and repair. Despite distinct mechanisms of leading and lagging strand replication3,4, we observe identical fidelity and damage tolerance for both strands. For small alkylation adducts of DNA, our results support a model in which the same translesion polymerase is recruited on-the-fly to both replication strands, starkly contrasting the strand asymmetric tolerance of bulky UV-induced adducts5. The accumulation of multiple distinct mutations at the site of persistent lesions provides the means to quantify the relative efficiency of repair processes genome wide and at single-base resolution. At multiple scales, we show DNA damage-induced mutations are largely shaped by the influence of DNA accessibility on repair efficiency, rather than gradients of DNA damage. Finally, we reveal specific genomic conditions that can actively drive oncogenic mutagenesis by corrupting the fidelity of nucleotide excision repair. These results provide insight into how strand-asymmetric mechanisms underlie the formation, tolerance and repair of DNA damage, thereby shaping cancer genome evolution.

Cohesin-mediated loop anchors confine the locations of human replication origins.

DNA replication occurs through an intricately regulated series of molecular events and is fundamental for genome stability1,2. At present, it is unknown how the locations of replication origins are determined in the human genome. Here we dissect the role of topologically associating domains (TADs)3-6, subTADs7 and loops8 in the positioning of replication initiation zones (IZs). We stratify TADs and subTADs by the presence of corner-dots indicative of loops and the orientation of CTCF motifs. We find that high-efficiency, early replicating IZs localize to boundaries between adjacent corner-dot TADs anchored by high-density arrays of divergently and convergently oriented CTCF motifs. By contrast, low-efficiency IZs localize to weaker dotless boundaries. Following ablation of cohesin-mediated loop extrusion during G1, high-efficiency IZs become diffuse and delocalized at boundaries with complex CTCF motif orientations. Moreover, G1 knockdown of the cohesin unloading factor WAPL results in gained long-range loops and narrowed localization of IZs at the same boundaries. Finally, targeted deletion or insertion of specific boundaries causes local replication timing shifts consistent with IZ loss or gain, respectively. Our data support a model in which cohesin-mediated loop extrusion and stalling at a subset of genetically encoded TAD and subTAD boundaries is an essential determinant of the locations of replication origins in human S phase.

Allele-specific control of replication timing and genome organization during development.

DNA replication occurs in a defined temporal order known as the replication-timing (RT) program. RT is regulated during development in discrete chromosomal units, coordinated with transcriptional activity and 3D genome organization. Here, we derived distinct cell types from F1 hybrid musculus × castaneus mouse crosses and exploited the high single-nucleotide polymorphism (SNP) density to characterize allelic differences in RT (Repli-seq), genome organization (Hi-C and promoter-capture Hi-C), gene expression (total nuclear RNA-seq), and chromatin accessibility (ATAC-seq). We also present HARP, a new computational tool for sorting SNPs in phased genomes to efficiently measure allele-specific genome-wide data. Analysis of six different hybrid mESC clones with different genomes (C57BL/6, 129/sv, and CAST/Ei), parental configurations, and gender revealed significant RT asynchrony between alleles across ∼12% of the autosomal genome linked to subspecies genomes but not to parental origin, growth conditions, or gender. RT asynchrony in mESCs strongly correlated with changes in Hi-C compartments between alleles but not as strongly with SNP density, gene expression, imprinting, or chromatin accessibility. We then tracked mESC RT asynchronous regions during development by analyzing differentiated cell types, including extraembryonic endoderm stem (XEN) cells, four male and female primary mouse embryonic fibroblasts (MEFs), and neural precursor cells (NPCs) differentiated in vitro from mESCs with opposite parental configurations. We found that RT asynchrony and allelic discordance in Hi-C compartments seen in mESCs were largely lost in all differentiated cell types, accompanied by novel sites of allelic asynchrony at a considerably smaller proportion of the genome, suggesting that genome organization of homologs converges to similar folding patterns during cell fate commitment.

DNA replication timing alterations identify common markers between distinct progeroid diseases.

Progeroid syndromes are rare genetic disorders that phenotypically resemble natural aging. Different causal mutations have been identified, but no molecular alterations have been identified that are in common to these diseases. DNA replication timing (RT) is a robust cell type-specific epigenetic feature highly conserved in the same cell types from different individuals but altered in disease. Here, we characterized DNA RT program alterations in Hutchinson-Gilford progeria syndrome (HGPS) and Rothmund-Thomson syndrome (RTS) patients compared with natural aging and cellular senescence. Our results identified a progeroid-specific RT signature that is common to cells from three HGPS and three RTS patients and distinguishes them from healthy individuals across a wide range of ages. Among the RT abnormalities, we identified the tumor protein p63 gene (TP63) as a gene marker for progeroid syndromes. By using the redifferentiation of four patient-derived induced pluripotent stem cells as a model for the onset of progeroid syndromes, we tracked the progression of RT abnormalities during development, revealing altered RT of the TP63 gene as an early event in disease progression of both HGPS and RTS. Moreover, the RT abnormalities in progeroid patients were associated with altered isoform expression of TP63 Our findings demonstrate the value of RT studies to identify biomarkers not detected by other methods, reveal abnormal TP63 RT as an early event in progeroid disease progression, and suggest TP63 gene regulation as a potential therapeutic target.

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

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

The many faces of the origin recognition complex.

The hetero-hexameric origin recognition complex (ORC) is well known for its separable roles in DNA replication and heterochromatin assembly. However, ORC and its individual subunits have been implicated in diverse cellular activities in both the nucleus and the cytoplasm. Some of ORC's implied functions, such as cell cycle checkpoint control and mitotic chromosome assembly, may be indirectly related to its roles in replication control and/or heterochromatin assembly. Other suggested roles in ribosomal biogenesis and in centrosome and kinetochore function are based on localization/interaction data and are as yet inconclusive. However, recent findings directly link ORC to sister chromatin cohesion, cytokinesis and neural dendritic branching.

Pre-replication complex proteins assemble at regions of low nucleosome occupancy within the Chinese hamster dihydrofolate reductase initiation zone.

Genome-scale mapping of pre-replication complex proteins has not been reported in mammalian cells. Poor enrichment of these proteins at specific sites may be due to dispersed binding, poor epitope availability or cell cycle stage-specific binding. Here, we have mapped sites of biotin-tagged ORC and MCM protein binding in G1-synchronized populations of Chinese hamster cells harboring amplified copies of the dihydrofolate reductase (DHFR) locus, using avidin-affinity purification of biotinylated chromatin followed by high-density microarray analysis across the DHFR locus. We have identified several sites of significant enrichment for both complexes distributed throughout the previously identified initiation zone. Analysis of the frequency of initiations across stretched DNA fibers from the DHFR locus confirmed a broad zone of de-localized initiation activity surrounding the sites of ORC and MCM enrichment. Mapping positions of mononucleosomal DNA empirically and computing nucleosome-positioning information in silico revealed that ORC and MCM map to regions of low measured and predicted nucleosome occupancy. Our results demonstrate that specific sites of ORC and MCM enrichment can be detected within a mammalian initiation zone, and suggest that initiation zones may be regions of generally low nucleosome occupancy where flexible nucleosome positioning permits flexible pre-RC assembly sites.

Brd2 is dispensable for genome compartmentalization and replication timing.

Replication Timing (RT) refers to the temporal order in which the genome is replicated during S phase. Early replicating regions correlate with the transcriptionally active, accessible euchromatin (A) compartment, while late replicating regions correlate with the heterochromatin (B) compartment and repressive histone marks. Previously, widespread A/B genome compartmentalization changes were reported following Brd2 depletion. Since RT and A/B compartmentalization are two of the most highly correlated chromosome properties, we evaluated the effects of Brd2 depletion on RT. We performed E/L Repli-Seq following Brd2 depletion in the previously described Brd2 conditional degron cell line and found no significant alterations in RT after Brd2 KD. This finding prompted us to re-analyze the Micro-C data from the previous publication. We report that we were unable to detect any compartmentalization changes in Brd2 depleted cells compared to DMSO control using the same data. Taken together, our findings demonstrate that Brd2 depletion alone does not affect A/B compartmentalization or RT in mouse embryonic stem cells.

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.

Nucleolus and centromere TSA-Seq reveals variable localization of heterochromatin in different cell types.

Genome differential positioning within interphase nuclei remains poorly explored. We extended and validated 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 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.

Replication timing: a fingerprint for cell identity and pluripotency.

Many types of epigenetic profiling have been used to classify stem cells, stages of cellular differentiation, and cancer subtypes. Existing methods focus on local chromatin features such as DNA methylation and histone modifications that require extensive analysis for genome-wide coverage. Replication timing has emerged as a highly stable cell type-specific epigenetic feature that is regulated at the megabase-level and is easily and comprehensively analyzed genome-wide. Here, we describe a cell classification method using 67 individual replication profiles from 34 mouse and human cell lines and stem cell-derived tissues, including new data for mesendoderm, definitive endoderm, mesoderm and smooth muscle. Using a Monte-Carlo approach for selecting features of replication profiles conserved in each cell type, we identify "replication timing fingerprints" unique to each cell type and apply a k nearest neighbor approach to predict known and unknown cell types. Our method correctly classifies 67/67 independent replication-timing profiles, including those derived from closely related intermediate stages. We also apply this method to derive fingerprints for pluripotency in human and mouse cells. Interestingly, the mouse pluripotency fingerprint overlaps almost completely with previously identified genomic segments that switch from early to late replication as pluripotency is lost. Thereafter, replication timing and transcription within these regions become difficult to reprogram back to pluripotency, suggesting these regions highlight an epigenetic barrier to reprogramming. In addition, the major histone cluster Hist1 consistently becomes later replicating in committed cell types, and several histone H1 genes in this cluster are downregulated during differentiation, suggesting a possible instrument for the chromatin compaction observed during differentiation. Finally, we demonstrate that unknown samples can be classified independently using site-specific PCR against fingerprint regions. In sum, replication fingerprints provide a comprehensive means for cell characterization and are a promising tool for identifying regions with cell type-specific organization.

SPIN reveals genome-wide landscape of nuclear compartmentalization.

We report SPIN, an integrative computational method to reveal genome-wide intranuclear chromosome positioning and nuclear compartmentalization relative to multiple nuclear structures, which are pivotal for modulating genome function. As a proof-of-principle, we use SPIN to integrate nuclear compartment mapping (TSA-seq and DamID) and chromatin interaction data (Hi-C) from K562 cells to identify 10 spatial compartmentalization states genome-wide relative to nuclear speckles, lamina, and putative associations with nucleoli. These SPIN states show novel patterns of genome spatial organization and their relation to other 3D genome features and genome function (transcription and replication timing). SPIN provides critical insights into nuclear spatial and functional compartmentalization.

Replication timing maintains the global epigenetic state in human cells.

The temporal order of DNA replication [replication timing (RT)] is correlated with chromatin modifications and three-dimensional genome architecture; however, causal links have not been established, largely because of an inability to manipulate the global RT program. We show that loss of RIF1 causes near-complete elimination of the RT program by increasing heterogeneity between individual cells. RT changes are coupled with widespread alterations in chromatin modifications and genome compartmentalization. Conditional depletion of RIF1 causes replication-dependent disruption of histone modifications and alterations in genome architecture. These effects were magnified with successive cycles of altered RT. These results support models in which the timing of chromatin replication and thus assembly plays a key role in maintaining the global epigenetic state.

High-resolution Repli-Seq defines the temporal choreography of initiation, elongation and termination of replication in mammalian cells.

BACKGROUND: DNA replication in mammalian cells occurs in a defined temporal order during S phase, known as the replication timing (RT) programme. Replication timing is developmentally regulated and correlated with chromatin conformation and local transcriptional potential. Here, we present RT profiles of unprecedented temporal resolution in two human embryonic stem cell lines, human colon carcinoma line HCT116, and mouse embryonic stem cells and their neural progenitor derivatives. RESULTS: Fine temporal windows revealed a remarkable degree of cell-to-cell conservation in RT, particularly at the very beginning and ends of S phase, and identified 5 temporal patterns of replication in all cell types, consistent with varying degrees of initiation efficiency. Zones of replication initiation (IZs) were detected throughout S phase and interacted in 3D space preferentially with other IZs of similar firing time. Temporal transition regions were resolved into segments of uni-directional replication punctuated at specific sites by small, inefficient IZs. Sites of convergent replication were divided into sites of termination or large constant timing regions consisting of many synchronous IZs in tandem. Developmental transitions in RT occured mainly by activating or inactivating individual IZs or occasionally by altering IZ firing time, demonstrating that IZs, rather than individual origins, are the units of developmental regulation. Finally, haplotype phasing revealed numerous regions of allele-specific and allele-independent asynchronous replication. Allele-independent asynchronous replication was correlated with the presence of previously mapped common fragile sites. CONCLUSIONS: Altogether, these data provide a detailed temporal choreography of DNA replication in mammalian cells.

3D genome organization contributes to genome instability at fragile sites.

Common fragile sites (CFSs) are regions susceptible to replication stress and are hotspots for chromosomal instability in cancer. Several features were suggested to underlie CFS instability, however, these features are prevalent across the genome. Therefore, the molecular mechanisms underlying CFS instability remain unclear. Here, we explore the transcriptional profile and DNA replication timing (RT) under mild replication stress in the context of the 3D genome organization. The results reveal a fragility signature, comprised of a TAD boundary overlapping a highly transcribed large gene with APH-induced RT-delay. This signature enables precise mapping of core fragility regions in known CFSs and identification of novel fragile sites. CFS stability may be compromised by incomplete DNA replication and repair in TAD boundaries core fragility regions leading to genomic instability. The identified fragility signature will allow for a more comprehensive mapping of CFSs and pave the way for investigating mechanisms promoting genomic instability in cancer.

An integrative ENCODE resource for cancer genomics.

ENCODE comprises thousands of functional genomics datasets, and the encyclopedia covers hundreds of cell types, providing a universal annotation for genome interpretation. However, for particular applications, it may be advantageous to use a customized annotation. Here, we develop such a custom annotation by leveraging advanced assays, such as eCLIP, Hi-C, and whole-genome STARR-seq on a number of data-rich ENCODE cell types. A key aspect of this annotation is comprehensive and experimentally derived networks of both transcription factors and RNA-binding proteins (TFs and RBPs). Cancer, a disease of system-wide dysregulation, is an ideal application for such a network-based annotation. Specifically, for cancer-associated cell types, we put regulators into hierarchies and measure their network change (rewiring) during oncogenesis. We also extensively survey TF-RBP crosstalk, highlighting how SUB1, a previously uncharacterized RBP, drives aberrant tumor expression and amplifies the effect of MYC, a well-known oncogenic TF. Furthermore, we show how our annotation allows us to place oncogenic transformations in the context of a broad cell space; here, many normal-to-tumor transitions move towards a stem-like state, while oncogene knockdowns show an opposing trend. Finally, we organize the resource into a coherent workflow to prioritize key elements and variants, in addition to regulators. We showcase the application of this prioritization to somatic burdening, cancer differential expression and GWAS. Targeted validations of the prioritized regulators, elements and variants using siRNA knockdowns, CRISPR-based editing, and luciferase assays demonstrate the value of the ENCODE resource.

The Chinese hamster dihydrofolate reductase replication origin decision point follows activation of transcription and suppresses initiation of replication within transcription units.

Chinese hamster ovary (CHO) cells select specific replication origin sites within the dihydrofolate reductase (DHFR) locus at a discrete point during G1 phase, the origin decision point (ODP). Origin selection is sensitive to transcription but not protein synthesis inhibitors, implicating a pretranslational role for transcription in origin specification. We have constructed a DNA array covering 121 kb surrounding the DHFR locus, to comprehensively investigate replication initiation and transcription in this region. When nuclei isolated within the first 3 h of G1 phase were stimulated to initiate replication in Xenopus egg extracts, replication initiated without any detectable preference for specific sites. At the ODP, initiation became suppressed from within the Msh3, DHFR, and 2BE2121 transcription units. Active transcription was mostly confined to these transcription units, and inhibition of transcription by alpha-amanitin resulted in the initiation of replication within transcription units, indicating that transcription is necessary to limit initiation events to the intergenic region. However, the resumption of DHFR transcription after mitosis took place prior to the ODP and so is not on its own sufficient to suppress initiation of replication. Together, these results demonstrate a remarkable flexibility in sequence selection for initiating replication and implicate transcription as one important component of origin specification at the ODP.

Replication timing alterations in leukemia affect clinically relevant chromosome domains.

Human B-cell precursor acute lymphoid leukemias (BCP-ALLs) comprise a group of genetically and clinically distinct disease entities with features of differentiation arrest at known stages of normal B-lineage differentiation. We previously showed that BCP-ALL cells display unique and clonally heritable, stable DNA replication timing (RT) programs (ie, programs describing the variable order of replication and subnuclear 3D architecture of megabase-scale chromosomal units of DNA in different cell types). To determine the extent to which BCP-ALL RT programs mirror or deviate from specific stages of normal human B-cell differentiation, we transplanted immunodeficient mice with quiescent normal human CD34+ cord blood cells and obtained RT signatures of the regenerating B-lineage populations. We then compared these with RT signatures for leukemic cells from a large cohort of BCP-ALL patients with varied genetic subtypes and outcomes. The results identify BCP-ALL subtype-specific features that resemble specific stages of B-cell differentiation and features that seem to be associated with relapse. These results suggest that the genesis of BCP-ALL involves alterations in RT that reflect biologically significant and potentially clinically relevant leukemia-specific epigenetic changes.

Continuous-Trait Probabilistic Model for Comparing Multi-species Functional Genomic Data.

A large amount of multi-species functional genomic data from high-throughput assays are becoming available to help understand the molecular mechanisms for phenotypic diversity across species. However, continuous-trait probabilistic models, which are key to such comparative analysis, remain under-explored. Here we develop a new model, called phylogenetic hidden Markov Gaussian processes (Phylo-HMGP), to simultaneously infer heterogeneous evolutionary states of functional genomic features in a genome-wide manner. Both simulation studies and real data application demonstrate the effectiveness of Phylo-HMGP. Importantly, we applied Phylo-HMGP to analyze a new cross-species DNA replication timing (RT) dataset from the same cell type in five primate species (human, chimpanzee, orangutan, gibbon, and green monkey). We demonstrate that our Phylo-HMGP model enables discovery of genomic regions with distinct evolutionary patterns of RT. Our method provides a generic framework for comparative analysis of multi-species continuous functional genomic signals to help reveal regions with conserved or lineage-specific regulatory roles.

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.