Insights into common fragile site instability: DNA replication challenges at DNA repeat sequences.

Common fragile sites (CFS) are specific genomic regions prone to chromosomal instability under conditions of DNA replication stress. CFSs manifest as breaks, gaps, and constrictions on metaphase chromosomes under mild replication stress. These replication-sensitive CFS regions are preferentially unstable during cancer development, as reflected by their association with copy number variants (CNVs) frequently arise in most tumor types. Over the years, it became clear that a combination of different characteristics underlies the enhanced sensitivity of CFSs to replication stress. As of today, there is a strong evidence that the core fragility regions along CFSs overlap with actively transcribed large genes with delayed replication timing upon replication stress. Recently, the mechanistic basis for CFS instability was further extended to regions which span topologically associated domain (TAD) boundaries, generating a fragility signature composed of replication, transcription and genome organization. The presence of difficult-to-replicate AT-rich repeats was one of the early features suggested to characterize a subgroup of CFSs. These long stretches of AT-dinucleotide have the potential to fold into stable secondary structures which may impede replication fork progression, leaving the region under-replicated. Here, we focus on the molecular mechanisms underlying repeat instability at CFSs and on the proteins involved in the resolution of secondary structure impediments arising along repetitive sequence elements which are essential for the maintenance of genome stability.

Mutational signatures association with replication timing in normal cells reveals similarities and differences with matched cancer tissues.

Mutational signatures' association with replication timing (RT) has been studied in cancer samples, but the RT distribution of somatic mutations in non-cancerous cells was only minimally explored. Here, we performed comprehensive analyses of mutational signatures in 2.9 million somatic mutations across multiple non-cancerous tissues, stratified by early and late RT regions. We found that many mutational processes are active mainly or solely in early RT, such as SBS16 in hepatocytes and SBS88 in the colon, or in late RT, such as SBS4 in lung and hepatocytes, and SBS18 across many tissues. The two ubiquitous signatures, SBS1 and SBS5, showed late and early bias, respectively, across multiple tissues and in mutations representing germ cells. We also performed a direct comparison with cancer samples in 4 matched tissue-cancer types. Unexpectedly, while for most signatures the RT bias was consistent in normal tissue and in cancer, we found that SBS1's late RT bias is lost in cancer.

Replication timing and genetic instability.

Synchronized activation of DNA replication origins induces genetic instability in lymphoma.

DNA replication timing directly regulates the frequency of oncogenic chromosomal translocations.

Chromosomal translocations result from the joining of DNA double-strand breaks (DSBs) and frequently cause cancer. However, the steps linking DSB formation to DSB ligation remain undeciphered. We report that DNA replication timing (RT) directly regulates lymphomagenic Myc translocations during antibody maturation in B cells downstream of DSBs and independently of DSB frequency. Depletion of minichromosome maintenance complexes alters replication origin activity, decreases translocations, and deregulates global RT. Ablating a single origin at Myc causes an early-to-late RT switch, loss of translocations, and reduced proximity with the immunoglobulin heavy chain (Igh) gene, its major translocation partner. These phenotypes were reversed by restoring early RT. Disruption of early RT also reduced tumorigenic translocations in human leukemic cells. Thus, RT constitutes a general mechanism in translocation biogenesis linking DSB formation to DSB ligation.

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.

RIF1 Links Replication Timing with Fork Reactivation and DNA Double-Strand Break Repair.

Replication timing (RT) is a cellular program to coordinate initiation of DNA replication in all origins within the genome. RIF1 (replication timing regulatory factor 1) is a master regulator of RT in human cells. This role of RIF1 is associated with binding G4-quadruplexes and changes in 3D chromatin that may suppress origin activation over a long distance. Many effects of RIF1 in fork reactivation and DNA double-strand (DSB) repair (DSBR) are underlined by its interaction with TP53BP1 (tumor protein p53 binding protein). In G1, RIF1 acts antagonistically to BRCA1 (BRCA1 DNA repair associated), suppressing end resection and homologous recombination repair (HRR) and promoting non-homologous end joining (NHEJ), contributing to DSBR pathway choice. RIF1 is an important element of intra-S-checkpoints to recover damaged replication fork with the involvement of HRR. High-resolution microscopic studies show that RIF1 cooperates with TP53BP1 to preserve 3D structure and epigenetic markers of genomic loci disrupted by DSBs. Apart from TP53BP1, RIF1 interact with many other proteins, including proteins involved in DNA damage response, cell cycle regulation, and chromatin remodeling. As impaired RT, DSBR and fork reactivation are associated with genomic instability, a hallmark of malignant transformation, RIF1 has a diagnostic, prognostic, and therapeutic potential in cancer. Further studies may reveal other aspects of common regulation of RT, DSBR, and fork reactivation by RIF1.

DNA methylation is required to maintain both DNA replication timing precision and 3D genome organization integrity.

DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

Fused in sarcoma regulates DNA replication timing and kinetics.

Fused in sarcoma (FUS) encodes an RNA-binding protein with diverse roles in transcriptional activation and RNA splicing. While oncogenic fusions of FUS and transcription factor DNA-binding domains are associated with soft tissue sarcomas, dominant mutations in FUS can cause amyotrophic lateral sclerosis. FUS has also been implicated in genome maintenance. However, the underlying mechanisms of its actions in genome stability are unknown. Here, we applied gene editing, functional reconstitution, and integrated proteomics and transcriptomics to illuminate roles for FUS in DNA replication and repair. Consistent with a supportive role in DNA double-strand break repair, FUS-deficient cells exhibited subtle alterations in the recruitment and retention of double-strand break-associated factors, including 53BP1 and BRCA1. FUS-/- cells also exhibited reduced proliferative potential that correlated with reduced speed of replication fork progression, diminished loading of prereplication complexes, enhanced micronucleus formation, and attenuated expression and splicing of S-phase-associated genes. Finally, FUS-deficient cells exhibited genome-wide alterations in DNA replication timing that were reversed upon re-expression of FUS complementary DNA. We also showed that FUS-dependent replication domains were enriched in transcriptionally active chromatin and that FUS was required for the timely replication of transcriptionally active DNA. These findings suggest that alterations in DNA replication kinetics and programming contribute to genome instability and functional defects in FUS-deficient cells.

Heterochromatin Networks: Topology, Dynamics, and Function (a Working Hypothesis).

Open systems can only exist by self-organization as pulsing structures exchanging matter and energy with the outer world. This review is an attempt to reveal the organizational principles of the heterochromatin supra-intra-chromosomal network in terms of nonlinear thermodynamics. The accessibility of the linear information of the genetic code is regulated by constitutive heterochromatin (CHR) creating the positional information in a system of coordinates. These features include scale-free splitting-fusing of CHR with the boundary constraints of the nucleolus and nuclear envelope. The analysis of both the literature and our own data suggests a radial-concentric network as the main structural organization principle of CHR regulating transcriptional pulsing. The dynamic CHR network is likely created together with nucleolus-associated chromatin domains, while the alveoli of this network, including springy splicing speckles, are the pulsing transcription hubs. CHR contributes to this regulation due to the silencing position variegation effect, stickiness, and flexible rigidity determined by the positioning of nucleosomes. The whole system acts in concert with the elastic nuclear actomyosin network which also emerges by self-organization during the transcriptional pulsing process. We hypothesize that the the transcriptional pulsing, in turn, adjusts its frequency/amplitudes specified by topologically associating domains to the replication timing code that determines epigenetic differentiation memory.

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

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

Genome-wide mapping of human DNA replication by optical replication mapping supports a stochastic model of eukaryotic replication.

The heterogeneous nature of eukaryotic replication kinetics and the low efficiency of individual initiation sites make mapping the location and timing of replication initiation in human cells difficult. To address this challenge, we have developed optical replication mapping (ORM), a high-throughput single-molecule approach, and used it to map early-initiation events in human cells. The single-molecule nature of our data and a total of >2,500-fold coverage of the human genome on 27 million fibers averaging ∼300 kb in length allow us to identify initiation sites and their firing probability with high confidence. We find that the distribution of human replication initiation is consistent with inefficient, stochastic activation of heterogeneously distributed potential initiation complexes enriched in accessible chromatin. These observations are consistent with stochastic models of initiation-timing regulation and suggest that stochastic regulation of replication kinetics is a fundamental feature of eukaryotic replication, conserved from yeast to humans.

Chromosomal coordination and differential structure of asynchronous replicating regions.

Stochastic asynchronous replication timing (AS-RT) is a phenomenon in which the time of replication of each allele is different, and the identity of the early allele varies between cells. By taking advantage of stable clonal pre-B cell populations derived from C57BL6/Castaneous mice, we have mapped the genome-wide AS-RT loci, independently of genetic differences. These regions are characterized by differential chromatin accessibility, mono-allelic expression and include new gene families involved in specifying cell identity. By combining population level mapping with single cell FISH, our data reveal the existence of a novel regulatory program that coordinates a fixed relationship between AS-RT regions on any given chromosome, with some loci set to replicate in a parallel and others set in the anti-parallel orientation. Our results show that AS-RT is a highly regulated epigenetic mark established during early embryogenesis that may be used for facilitating the programming of mono-allelic choice throughout development.

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.

High-resolution mapping of mitotic DNA synthesis regions and common fragile sites in the human genome through direct sequencing.

DNA replication stress, a feature of human cancers, often leads to instability at specific genomic loci, such as the common fragile sites (CFSs). Cells experiencing DNA replication stress may also exhibit mitotic DNA synthesis (MiDAS). To understand the physiological function of MiDAS and its relationship to CFSs, we mapped, at high resolution, the genomic sites of MiDAS in cells treated with the DNA polymerase inhibitor aphidicolin. Sites of MiDAS were evident as well-defined peaks that were largely conserved between cell lines and encompassed all known CFSs. The MiDAS peaks mapped within large, transcribed, origin-poor genomic regions. In cells that had been treated with aphidicolin, these regions remained unreplicated even in late S phase; MiDAS then served to complete their replication after the cells entered mitosis. Interestingly, leading and lagging strand synthesis were uncoupled in MiDAS, consistent with MiDAS being a form of break-induced replication, a repair mechanism for collapsed DNA replication forks. Our results provide a better understanding of the mechanisms leading to genomic instability at CFSs and in cancer cells.

Genome-wide high-resolution mapping of mitotic DNA synthesis sites and common fragile sites by direct sequencing.

Common fragile sites (CFSs) are genomic loci prone to the formation of breaks or gaps on metaphase chromosomes. They are hotspots for chromosome rearrangements and structural variations, which have been extensively implicated in carcinogenesis, aging, and other pathological processes. Although many CFSs were identified decades ago, a consensus is still lacking for why they are particularly unstable and sensitive to replication perturbations. This is in part due to the lack of high-resolution mapping data for the vast majority of the CFSs, which has hindered mechanistic interrogations. Here, we seek to map human CFSs with high resolution on a genome-wide scale by sequencing the sites of mitotic DNA synthesis (MiDASeq) that are specific for CFSs. We generated a nucleotide-resolution atlas of MiDAS sites (MDSs) that covered most of the known CFSs, and comprehensively analyzed their sequence characteristics and genomic features. Our data on MDSs tallied well with long-standing hypotheses to explain CFS fragility while highlighting the contributions of late replication timing and large transcription units. Notably, the MDSs also encompassed most of the recurrent double-strand break clusters previously identified in mouse neural stem/progenitor cells, thus bridging evolutionarily conserved break points across species. Moreover, MiDAseq provides an important resource that can stimulate future research on CFSs to further unravel the mechanisms and biological relevance underlying these labile genomic regions.

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.

Using GARDEN-NET and ChAseR to explore human haematopoietic 3D chromatin interaction networks.

We introduce an R package and a web-based visualization tool for the representation, analysis and integration of epigenomic data in the context of 3D chromatin interaction networks. GARDEN-NET allows for the projection of user-submitted genomic features on pre-loaded chromatin interaction networks, exploiting the functionalities of the ChAseR package to explore the features in combination with chromatin network topology properties. We demonstrate the approach using published epigenomic and chromatin structure datasets in haematopoietic cells, including a collection of gene expression, DNA methylation and histone modifications data in primary healthy myeloid cells from hundreds of individuals. These datasets allow us to test the robustness of chromatin assortativity, which highlights which epigenomic features, alone or in combination, are more strongly associated with 3D genome architecture. We find evidence for genomic regions with specific histone modifications, DNA methylation, and gene expression levels to be forming preferential contacts in 3D nuclear space, to a different extent depending on the cell type and lineage. Finally, we examine replication timing data and find it to be the genomic feature most strongly associated with overall 3D chromatin organization at multiple scales, consistent with previous results from the literature.

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.

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.

Synchronized replication of genes encoding the same protein complex in fast-proliferating cells.

DNA replication perturbs the dosage balance among genes; at mid-S phase, early-replicating genes have doubled their copies while late-replicating ones have not. Dosage imbalance among genes, especially within members of a protein complex, is toxic to cells. However, the molecular mechanisms that cells use to deal with such imbalance remain not fully understood. Here, we validate at the genomic scale that the dosage between early- and late-replicating genes is imbalanced in HeLa cells. We propose the synchronized replication hypothesis that genes sensitive to stoichiometric relationships will be replicated simultaneously to maintain stoichiometry. In support of this hypothesis, we observe that genes encoding the same protein complex have similar replication timing but mainly in fast-proliferating cells such as embryonic stem cells and cancer cells. We find that the synchronized replication observed in cancer cells, but not in slow-proliferating differentiated cells, is due to convergent evolution during tumorigenesis that restores synchronized replication timing within protein complexes. Taken together, our study reveals that the demand for dosage balance during S phase plays an important role in the optimization of the replication-timing program; this selection is relaxed during differentiation as the cell cycle prolongs and is restored during tumorigenesis as the cell cycle shortens.