Genome-wide analysis of DNA replication and DNA double-strand breaks using TrAEL-seq

Human DNA replication initiation sites are specified epigenetically by oxidation of 5-methyl-deoxycytidine

DNA replication initiates at tens of thousands of sites on the human genome during each S phase. However, no consensus DNA sequence has been found that specifies the locations of these replication origins. Here, we investigate modifications of human genomic DNA by density equilibrium centrifugation and DNA sequencing. We identified short discrete sites with increased density during quiescence and G1 phase that overlap with DNA replication origins before their activation in S phase. The increased density is due to the oxidation of 5-methyl-deoxycytidines by ten-eleven-translocation DNA dioxygenase (TET) enzymes at GC-rich domains. Reversible inhibition of de novo methylation and of subsequent oxidation of deoxycytidines results in a reversible inhibition of DNA replication and of cell proliferation. Our findings suggest a mechanism for the epigenetic specification and semiconservative inheritance of DNA replication origin sites in human cells that also provides a stable integral DNA replication licence to support once-per-cell cycle control of origin activation.

Keeping it safe: control of meiotic chromosome breakage

Meiotic cells introduce numerous programmed DNA double-strand breaks (DSBs) into their genome to stimulate crossover recombination. DSB numbers must be high enough to ensure each homologous chromosome pair receives the obligate crossover required for accurate meiotic chromosome segregation. However, every DSB also increases the risk of aberrant or incomplete DNA repair, and thus genome instability. To mitigate these risks, meiotic cells have evolved an intricate network of controls that modulates the timing, levels, and genomic location of meiotic DSBs. This Review summarizes our current understanding of these controls with a particular focus on the mechanisms that prevent meiotic DSB formation at the wrong time or place, thereby guarding the genome from potentially catastrophic meiotic errors.

Distinct chromatin regulators downmodulate meiotic axis formation and DNA break induction at chromosome ends

In many organisms, meiotic crossover recombination is suppressed near the extreme ends of chromosomes. Here, we identified two chromatin modifiers, the histone methyltransferase Dot1 and the Sir silencing complex, as regulators of this process in Saccharomyces cerevisiae. We show that the recombination-promoting axis proteins Red1 and Hop1, but not the axis-associated cohesin Rec8, are significantly reduced within 20 kb of telomeres compared to the chromosome interior. Dot1, which preferentially methylates histones in the chromosome interior, is required for this pattern by directing Red1 binding toward the chromosome interior. In parallel, the Sir complex suppresses the induction of meiotic DNA double-strand breaks (DSBs) at chromosome ends. Sir-dependent DSB suppression is independent of axis deposition and occurs in a chromosome end-specific manner that mirrors the spreading and transcriptional silencing activity of the complex, suggesting that the Sir complex suppresses DSB formation by limiting the openness of promoters, the preferred sites of meiotic DSB formation. We conclude that multiple chromatin-based mechanisms collaborate to achieve a robust reduction of meiotic recombination near chromosome ends.

The double life of mammalian DNA replication origins

Mammalian DNA replication origins have been historically difficult to identify and their determinants are still unresolved. Here, we first review methods developed over the last decades to map replication initiation sites either directly via initiation intermediates or indirectly via determining replication fork directionality profiles. We also discuss the factors that may specify these sites as replication initiation sites. Second, we address the controversy that has emerged from these results over whether origins are narrowly defined and localized to specific sites or are more dispersed and organized into broad zones. Ample evidence in favor of both scenarios currently creates an impression of unresolved confusion in the field. We attempt to formulate a synthesis of both models and to reconcile discrepant findings. It is evident that not only one approach is sufficient in isolation but that the combination of several is instrumental toward understanding initiation sites in mammalian genomes. We argue that an aggregation of several individual and often inefficient initiation sites into larger initiation zones and the existence of efficient unidirectional initiation sites and fork stalling at the borders of initiation zones can reconcile the different observations.

A tale of two strands: Decoding chromatin replication through strand-specific sequencing

DNA replication, a fundamental process in all living organisms, proceeds with continuous synthesis of the leading strand by DNA polymerase ε (Pol ε) and discontinuous synthesis of the lagging strand by polymerase δ (Pol δ). This inherent asymmetry at each replication fork necessitates the development of methods to distinguish between these two nascent strands in vivo. Over the past decade, strand-specific sequencing strategies, such as enrichment and sequencing of protein-associated nascent DNA (eSPAN) and Okazaki fragment sequencing (OK-seq), have become essential tools for studying chromatin replication in eukaryotic cells. In this review, we outline the foundational principles underlying these methodologies and summarize key mechanistic insights into DNA replication, parental histone transfer, epigenetic inheritance, and beyond, gained through their applications. Finally, we discuss the limitations and challenges of current techniques, highlighting the need for further technological innovations to better understand the dynamics and regulation of chromatin replication in eukaryotic cells.

Two-ended recombination at a Flp-nickase-broken replication fork

Replication fork collision with a DNA nick can generate a one-ended break, fostering genomic instability. The opposing fork's collision with the nick could form a second DNA end, enabling conservative repair by homologous recombination (HR). To study mechanisms of nickase-induced HR, we developed the Flp recombinase "step arrest" nickase in mammalian cells. A Flp-nick induces two-ended, BRCA2/RAD51-dependent short tract gene conversion (STGC), BRCA2/RAD51-independent long tract gene conversion, and discoordinated two-ended invasions. HR pathways induced by a replication-independent break and the Flp-nickase differ in their dependence on BRCA1, MRE11, and CtIP. To determine the origin of the second DNA end during Flp-nickase-induced STGC, we blocked the opposing fork using a Tus/Ter replication fork barrier (RFB). Flp-nickase-induced STGC remained robust and two ended. Thus, a single replication fork's collision with a Flp-nick triggers two-ended HR, possibly reflecting replicative bypass of lagging strand nicks. This response may limit genomic instability during replication of nicked DNA.

Precise Mapping of Physiological DSBs Using In-Suspension Break Labeling In Situ and Sequencing (sBLISS)

DNA double-strand breaks (DSBs) are a major source of genomic instability. Physiological DSBs are naturally occurring breaks that happen during normal cellular processes. Unlike DNA breaks resulting from DNA damage due to external factors like radiation or chemicals, physiological DSBs play critical roles in various normal biological functions. Some key processes involving physiological DSBs include V(D)J recombination, transcription, and replication. These breaks are typically tightly controlled and are part of the cellular machinery designed to maintain and enhance genomic integrity and diversity. However, if these breaks are misrepaired or left unrepaired, they can contribute to genomic instability, potentially leading to senescence and diseases such as cancer. Here, we outline various methods commonly employed to detect physiological DSBs and introduce a detailed, step-by-step protocol for mapping these breaks using the in-suspension break labeling in situ and sequencing (sBLISS) technique. sBLISS offers single nucleotide resolution and is versatile enough to be applied to any cell type amenable to single-cell suspension. This comprehensive approach not only enhances our understanding of DSBs but also aids in the exploration of their roles in genomic instability.

Mutations in the DNA processivity factor POL30 predispose the FLO11 locus to epigenetic instability in S. cerevisiae

The FLO genes in Saccharomyces cerevisiae are repressed by heterochromatin formation, involving histone deacetylases, transcription factors and non-coding RNAs. Here, we report that mutations in the processivity factor POL30 (PCNA) that show transient derepression at the subtelomeres and the mating-type loci do not derepress FLO loci. However, deletions of the replisome stability factors RRM3 and TOF1 along with pol30 mutations induced flocculation phenotypes. The phenotypes correlated with increased expression of reporter proteins driven by the FLO11 promoter, the frequency of silent to active conversions of FLO11, and reduced expression of the regulatory long non-coding RNAs ICR1 and PWR1. Alterations in the local replication landscape of FLO11 indicate a link between defects in the fork protection complex and the stability of gene silencing. Analyses of these mutants at the subtelomeres and the HMLα locus showed a similar derepression phenotype and suggest transient instability of both active and silent states of FLO11. We conclude that RRM3 and TOF1 interact differentially with the pol30 mutations to promote transient derepression or complete epigenetic conversions of FLO11. We suggest that the interaction between POL30, RRM3 and TOF1 is essential to maintain epigenetic stability at the studied loci.

The fork protection complex generates DNA topological stress-induced DNA damage while ensuring full and faithful genome duplication

The fork protection complex (FPC), composed of Mrc1, Tof1, and Csm3, supports rapid and stable DNA replication. Here, we show that FPC activity also introduces DNA damage by increasing DNA topological stress during replication. Mrc1 action increases DNA topological stress during plasmid replication, while Mrc1 or Tof1 activity causes replication stress and DNA damage within topologically constrained regions. We show that the recruitment of Top1 to the fork by Tof1 suppresses the DNA damage generated in these loci. While FPC activity introduces some DNA damage due to increased topological stress, the FPC is also necessary to prevent DNA damage in long replicons across the genome, indicating that the FPC is required for complete and faithful genome duplication. We conclude that FPC regulation must balance ensuring full genome duplication through rapid replication with minimizing the consequential DNA topological stress-induced DNA damage caused by rapid replication through constrained regions.

Replication initiation sites and zones in the mammalian genome: Where are they located and how are they defined?

Eukaryotic DNA replication is a tightly controlled process that occurs in two main steps, i.e., licensing and firing, which take place in the G1 and S phases of the cell cycle, respectively. In Saccharomyces cerevisiae, the budding yeast, replication origins contain consensus sequences that are recognized and bound by the licensing factor Orc1-6, which then recruits the replicative Mcm2-7 helicase. By contrast, mammalian initiation sites lack such consensus sequences, and the mammalian ORC does not exhibit sequence specificity. Studies performed over the past decades have identified replication initiation sites in the mammalian genome using sequencing-based assays, raising the question of whether replication initiation occurs at confined sites or in broad zones across the genome. Although recent reports have shown that the licensed MCMs in mammalian cells are broadly distributed, suggesting that ORC-dependent licensing may not determine the initiation sites/zones, they are predominantly located upstream of actively transcribed genes. This review compares the mechanism of replication initiation in yeast and mammalian cells, summarizes the sequencing-based technologies used for the identification of initiation sites/zones, and proposes a possible mechanism of initiation-site/zone selection in mammalian cells. Future directions and challenges in this field are also discussed.

Structure and repair of replication-coupled DNA breaks

Using CRISPR-Cas9 nicking enzymes, we examined the interaction between the replication machinery and single-strand breaks, one of the most common forms of endogenous DNA damage. We show that replication fork collapse at leading-strand nicks generates resected single-ended double-strand breaks (seDSBs) that are repaired by homologous recombination (HR). If these seDSBs are not promptly repaired, arrival of adjacent forks creates double-ended DSBs (deDSBs), which could drive genomic scarring in HR-deficient cancers. deDSBs can also be generated directly when the replication fork bypasses lagging-strand nicks. Unlike deDSBs produced independently of replication, end resection at nick-induced seDSBs and deDSBs is BRCA1-independent. Nevertheless, BRCA1 antagonizes 53BP1 suppression of RAD51 filament formation. These results highlight distinctive mechanisms that maintain replication fork stability.

Novel WRN Helicase Inhibitors Selectively Target Microsatellite-Unstable Cancer Cells

Microsatellite-unstable (MSI) cancers require WRN helicase to resolve replication stress due to expanded DNA (TA)n dinucleotide repeats. WRN is a promising synthetic lethal target for MSI tumors, and WRN inhibitors are in development. In this study, we used CRISPR-Cas9 base editing to map WRN residues critical for MSI cells, validating the helicase domain as the primary drug target. Fragment-based screening led to the development of potent and highly selective WRN helicase covalent inhibitors. These compounds selectively suppressed MSI model growth in vitro and in vivo by mimicking WRN loss, inducing DNA double-strand breaks at expanded TA repeats and DNA damage. Assessment of biomarkers in preclinical models linked TA-repeat expansions and mismatch repair alterations to compound activity. Efficacy was confirmed in immunotherapy-resistant organoids and patient-derived xenograft models. The discovery of potent, selective covalent WRN inhibitors provides proof of concept for synthetic lethal targeting of WRN in MSI cancer and tools to dissect WRN biology. Significance: We report the discovery and characterization of potent, selective WRN helicase inhibitors for MSI cancer treatment, with biomarker analysis and evaluation of efficacy in vivo and in immunotherapy-refractory preclinical models. These findings pave the way to translate WRN inhibition into MSI cancer therapies and provide tools to investigate WRN biology. See related commentary by Wainberg, p. 1369.

Monitoring and quantifying replication fork dynamics with high-throughput methods

Before each cell division, eukaryotic cells must replicate their chromosomes to ensure the accurate transmission of genetic information. Chromosome replication involves more than just DNA duplication; it also includes chromatin assembly, inheritance of epigenetic marks, and faithful resumption of all genomic functions after replication. Recent progress in quantitative technologies has revolutionized our understanding of the complexity and dynamics of DNA replication forks at both molecular and genomic scales. Here, we highlight the pivotal role of these novel methods in uncovering the principles and mechanisms of chromosome replication. These technologies have illuminated the regulation of genome replication programs, quantified the impact of DNA replication on genomic mutations and evolutionary processes, and elucidated the mechanisms of replication-coupled chromatin assembly and epigenome maintenance.

S-phase checkpoint prevents leading strand degradation from strand-associated nicks at stalled replication forks

The S-phase checkpoint is involved in coupling DNA unwinding with nascent strand synthesis and is critical to maintain replication fork stability in conditions of replicative stress. However, its role in the specific regulation of leading and lagging strands at stalled forks is unclear. By conditionally depleting RNaseH2 and analyzing polymerase usage genome-wide, we examine the enzymology of DNA replication during a single S-phase in the presence of replicative stress and show that there is a differential regulation of lagging and leading strands. In checkpoint proficient cells, lagging strand replication is down-regulated through an Elg1-dependent mechanism. Nevertheless, when checkpoint function is impaired we observe a defect specifically at the leading strand, which was partially dependent on Exo1 activity. Further, our genome-wide mapping of DNA single-strand breaks reveals that strand discontinuities highly accumulate at the leading strand in HU-treated cells, whose dynamics are affected by checkpoint function and Exo1 activity. Our data reveal an unexpected role of Exo1 at the leading strand and support a model of fork stabilization through prevention of unrestrained Exo1-dependent resection of leading strand-associated nicks after fork stalling.

Two-ended recombination at a Flp-nickase-broken replication fork

Collision of a replication fork with a DNA nick is thought to generate a one-ended break, fostering genomic instability. Collision of the opposing converging fork with the nick could, in principle, form a second DNA end, enabling conservative repair by homologous recombination (HR). To study mechanisms of nickase-induced HR, we developed the Flp recombinase "step arrest" nickase in mammalian cells. Flp-nickase-induced HR entails two-ended, BRCA2/RAD51-dependent short tract gene conversion (STGC), BRCA2/RAD51-independent long tract gene conversion, and discoordinated two-ended invasions. HR induced by a replication-independent break and by the Flp-nickase differ in their dependence on BRCA1 . To determine the origin of the second DNA end during Flp-nickase-induced STGC, we blocked the opposing fork using a site-specific Tus/ Ter replication fork barrier. Flp-nickase-induced STGC remained robust and two-ended. Thus, collision of a single replication fork with a Flp-nick can trigger two-ended HR, possibly reflecting replicative bypass of lagging strand nicks. This response may limit genomic instability during replication of a nicked DNA template.

Transcription as source of genetic heterogeneity in budding yeast

Transcription presents challenges to genome stability both directly, by altering genome topology and exposing single-stranded DNA to chemical insults and nucleases, and indirectly by introducing obstacles to the DNA replication machinery. Such obstacles include the RNA polymerase holoenzyme itself, DNA-bound regulatory factors, G-quadruplexes and RNA-DNA hybrid structures known as R-loops. Here, we review the detrimental impacts of transcription on genome stability in budding yeast, as well as the mitigating effects of transcription-coupled nucleotide excision repair and of systems that maintain DNA replication fork processivity and integrity. Interactions between DNA replication and transcription have particular potential to induce mutation and structural variation, but we conclude that such interactions must have only minor effects on DNA replication by the replisome with little if any direct mutagenic outcome. However, transcription can significantly impair the fidelity of replication fork rescue mechanisms, particularly Break Induced Replication, which is used to restart collapsed replication forks when other means fail. This leads to de novo mutations, structural variation and extrachromosomal circular DNA formation that contribute to genetic heterogeneity, but only under particular conditions and in particular genetic contexts, ensuring that the bulk of the genome remains extremely stable despite the seemingly frequent interactions between transcription and DNA replication.

A Cluster of Evolutionarily Recent KRAB Zinc Finger Proteins Protects Cancer Cells from Replicative Stress-Induced Inflammation

Heterochromatin loss and genetic instability enhance cancer progression by favoring clonal diversity, yet uncontrolled replicative stress leads to mitotic catastrophe and inflammatory responses that promote immune rejection. KRAB domain-containing zinc finger proteins (KZFP) contribute to heterochromatin maintenance at transposable elements (TE). Here, we identified an association of upregulation of a cluster of primate-specific KZFPs with poor prognosis, increased copy-number alterations, and changes in the tumor microenvironment in diffuse large B-cell lymphoma (DLBCL). Depleting two of these KZFPs targeting evolutionarily recent TEs, ZNF587 and ZNF417, impaired the proliferation of cells derived from DLBCL and several other tumor types. ZNF587 and ZNF417 depletion led to heterochromatin redistribution, replicative stress, and cGAS-STING-mediated induction of an interferon/inflammatory response, which enhanced susceptibility to macrophage-mediated phagocytosis and increased surface expression of HLA-I, together with presentation of a neoimmunopeptidome. Thus, cancer cells can exploit KZFPs to dampen TE-originating surveillance mechanisms, which likely facilitates clonal expansion, diversification, and immune evasion.

SIGNIFICANCE

Upregulation of a cluster of primate-specific KRAB zinc finger proteins in cancer cells prevents replicative stress and inflammation by regulating heterochromatin maintenance, which could facilitate the development of improved biomarkers and treatments.

SSBlazer: a genome-wide nucleotide-resolution model for predicting single-strand break sites

Single-strand breaks are the major DNA damage in the genome and serve a crucial role in various biological processes. To reveal the significance of single-strand breaks, multiple sequencing-based single-strand break detection methods have been developed, which are costly and unfeasible for large-scale analysis. Hence, we propose SSBlazer, an explainable and scalable deep learning framework for single-strand break site prediction at the nucleotide level. SSBlazer is a lightweight model with robust generalization capabilities across various species and is capable of numerous unexplored SSB-related applications.

DNA Damage Atlas: an atlas of DNA damage and repair

DNA damage and its improper repair are the major source of genomic alterations responsible for many human diseases, particularly cancer. To aid researchers in understanding the underlying mechanisms of genome instability, a number of genome-wide profiling approaches have been developed to monitor DNA damage and repair events. The rapid accumulation of published datasets underscores the critical necessity of a comprehensive database to curate sequencing data on DNA damage and repair intermediates. Here, we present DNA Damage Atlas (DDA, http://www.bioinformaticspa.com/DDA/), the first large-scale repository of DNA damage and repair information. Currently, DDA comprises 6,030 samples from 262 datasets by 59 technologies, covering 16 species, 10 types of damage and 135 treatments. Data collected in DDA was processed through a standardized workflow, including quality checks, hotspots identification and a series of feature characterization for the hotspots. Notably, DDA encompasses analyses of highly repetitive regions, ribosomal DNA and telomere. DDA offers a user-friendly interface that facilitates browsing, searching, genome browser visualization, hotspots comparison and data downloading, enabling convenient and thorough exploration for datasets of interest. In summary, DDA will stand as a valuable resource for research in genome instability and its association with diseases.

Checkpoint kinase interaction with DNA polymerase alpha regulates replication progression during stress

Background: In eukaryotes, replication stress activates a checkpoint response, which facilitates genome duplication by stabilising the replisome. How the checkpoint kinases regulate the replisome remains poorly understood. The aim of this study is to identify new targets of checkpoint kinases within the replisome during replication stress. Methods: Here we use an unbiased biotin proximity-ligation approach in Saccharomyces cerevisiae to identify new interactors and substrates of the checkpoint kinase Rad53 in vivo. Results: From this screen, we identified the replication initiation factor Sld7 as a Rad53 substrate, and Pol1, the catalytic subunit of polymerase a, as a Rad53-interactor. We showed that CDK phosphorylation of Pol1 mediates its interaction with Rad53. Combined with other interactions between Rad53 and the replisome, this Rad53-Pol1 interaction is important for viability and replisome progression during replication stress. Conclusions: Together, we explain how the interactions of Rad53 with the replisome are controlled by both replication stress and the cell cycle, and why these interactions might be important for coordinating the stabilisation of both the leading and lagging strand machineries.

DEtail-seq is an ultra-efficient and convenient method for meiotic DNA break profiling in multiple organisms

Programmed DNA double-strand break (DSB) formation is a crucial step in meiotic recombination, yet techniques for high-efficiency and precise mapping of the 3' ends of DSBs are still in their infancy. Here, we report a novel technique, named DNA End tailing and sequencing (DEtail-seq), which can directly and ultra-efficiently characterize the 3' ends of meiotic DSBs with near single-nucleotide resolution in a variety of species, including yeast, mouse, and human. We find that the 3' ends of meiotic DSBs are stable without significant resection in budding yeast. Meiotic DSBs are strongly enriched in de novo H3K4me3 peaks in the mouse genome at leptotene stage. We also profile meiotic DSBs in human and find DSB hotspots are enriched near the common fragile sites during human meiosis, especially at CCCTC-binding factor (CTCF)-associated enhancers. Therefore, DEtail-seq provides a powerful method to detect DSB ends in various species, and our results provide new insights into the distribution and regulation of meiotic DSB hotspots.

OKseqHMM: a genome-wide replication fork directionality analysis toolkit

During each cell division, tens of thousands of DNA replication origins are co-ordinately activated to ensure the complete duplication of the human genome. However, replication fork progression can be challenged by many factors, including co-directional and head-on transcription-replication conflicts (TRC). Head-on TRCs are more dangerous for genome integrity. To study the direction of replication fork movement and TRCs, we developed a bioinformatics toolkit called OKseqHMM (https://github.com/CL-CHEN-Lab/OK-Seq, https://doi.org/10.5281/zenodo.7428883). Then, we used OKseqHMM to analyse a large number of datasets obtained by Okazaki fragment sequencing to directly measure the genome-wide replication fork directionality (RFD) and to accurately predict replication initiation and termination at a fine resolution in organisms including yeast, mouse and human. We also successfully applied our analysis to other genome-wide sequencing techniques that also contain RFD information (e.g. eSPAN, TrAEL-seq). Our toolkit can be used to predict replication initiation and fork progression direction genome-wide in a wide range of cell models and growth conditions. Comparing the replication and transcription directions allows identifying loci at risk of TRCs, particularly head-on TRCs, and investigating their role in genome instability by checking DNA damage data, which is of prime importance for human health.

Genome-wide measurement of DNA replication fork directionality and quantification of DNA replication initiation and termination with Okazaki fragment sequencing

Studying the dynamics of genome replication in mammalian cells has been historically challenging. To reveal the location of replication initiation and termination in the human genome, we developed Okazaki fragment sequencing (OK-seq), a quantitative approach based on the isolation and strand-specific sequencing of Okazaki fragments, the lagging strand replication intermediates. OK-seq quantitates the proportion of leftward- and rightward-oriented forks at every genomic locus and reveals the location and efficiency of replication initiation and termination events. Here we provide the detailed experimental procedures for performing OK-seq in unperturbed cultured human cells and budding yeast and the bioinformatics pipelines for data processing and computation of replication fork directionality. Furthermore, we present the analytical approach based on a hidden Markov model, which allows automated detection of ascending, descending and flat replication fork directionality segments revealing the zones of replication initiation, termination and unidirectional fork movement across the entire genome. These tools are essential for the accurate interpretation of human and yeast replication programs. The experiments and the data processing can be accomplished within six days. Besides revealing the genome replication program in fine detail, OK-seq has been instrumental in numerous studies unravelling mechanisms of genome stability, epigenome maintenance and genome evolution.

From fluorescent foci to sequence: Illuminating DNA double strand break repair by high-throughput sequencing technologies

Technologies to study DNA double-strand break (DSB) repair have traditionally mostly relied on fluorescence read-outs, either by microscopy or flow cytometry. The advent of high throughput sequencing (HTS) has created fundamentally new opportunities to study the mechanisms underlying DSB repair. Here, we review the suite of HTS-based assays that are used to study three different aspects of DNA repair: detection of broken ends, protein recruitment and pathway usage. We highlight new opportunities that HTS technology offers towards a better understanding of the DSB repair process.

RRM2 enhances MYCN-driven neuroblastoma formation and acts as a synergistic target with CHK1 inhibition

High-risk neuroblastoma, a pediatric tumor originating from the sympathetic nervous system, has a low mutation load but highly recurrent somatic DNA copy number variants. Previously, segmental gains and/or amplifications allowed identification of drivers for neuroblastoma development. Using this approach, combined with gene dosage impact on expression and survival, we identified ribonucleotide reductase subunit M2 (RRM2) as a candidate dependency factor further supported by growth inhibition upon in vitro knockdown and accelerated tumor formation in a neuroblastoma zebrafish model coexpressing human RRM2 with MYCN. Forced RRM2 induction alleviates excessive replicative stress induced by CHK1 inhibition, while high RRM2 expression in human neuroblastomas correlates with high CHK1 activity. MYCN-driven zebrafish tumors with RRM2 co-overexpression exhibit differentially expressed DNA repair genes in keeping with enhanced ATR-CHK1 signaling activity. In vitro, RRM2 inhibition enhances intrinsic replication stress checkpoint addiction. Last, combinatorial RRM2-CHK1 inhibition acts synergistic in high-risk neuroblastoma cell lines and patient-derived xenograft models, illustrating the therapeutic potential.

Stimulation of adaptive gene amplification by origin firing under replication fork constraint

Adaptive mutations can cause drug resistance in cancers and pathogens, and increase the tolerance of agricultural pests and diseases to chemical treatment. When and how adaptive mutations form is often hard to discern, but we have shown that adaptive copy number amplification of the copper resistance gene CUP1 occurs in response to environmental copper due to CUP1 transcriptional activation. Here we dissect the mechanism by which CUP1 transcription in budding yeast stimulates copy number variation (CNV). We show that transcriptionally stimulated CNV requires TREX-2 and Mediator, such that cells lacking TREX-2 or Mediator respond normally to copper but cannot acquire increased resistance. Mediator and TREX-2 can cause replication stress by tethering transcribed loci to nuclear pores, a process known as gene gating, and transcription at the CUP1 locus causes a TREX-2-dependent accumulation of replication forks indicative of replication fork stalling. TREX-2-dependent CUP1 gene amplification occurs by a Rad52 and Rad51-mediated homologous recombination mechanism that is enhanced by histone H3K56 acetylation and repressed by Pol32 and Pif1. CUP1 amplification is also critically dependent on late-firing replication origins present in the CUP1 repeats, and mutations that remove or inactivate these origins strongly suppress the acquisition of copper resistance. We propose that replicative stress imposed by nuclear pore association causes replication bubbles from these origins to collapse soon after activation, leaving a tract of H3K56-acetylated chromatin that promotes secondary recombination events during elongation after replication fork re-start events. The capacity for inefficient replication origins to promote copy number variation renders certain genomic regions more fragile than others, and therefore more likely to undergo adaptive evolution through de novo gene amplification.

Top1p targeting by Fob1p at the ribosomal Replication Fork Barrier does not account for camptothecin sensitivity in Saccharomyces cerevisiae cells

Camptothecin (CPT) is a specific inhibitor of the DNA topoisomerase I (Top1p), currently used in cancer therapy, which induces DNA damage and cell death. Top1p is highly active at the repeated ribosomal DNA locus (rDNA) to relax DNA supercoiling caused by elevated transcription and replication occurring in opposite directions. Fob1p interacts with, and stabilizes, Top1p at the rDNA Replication Fork Barrier (rRFB), where replication and transcription converge. Here, we have investigated if the absence of Fob1p and the consequent loss of Top1p specific targeting to the rRFB impact the sensitivity and the cell cycle progression of wild-type cells to CPT. We have also investigated the consequences of the absence of Fob1p in rad52∆ mutants, which are affected in the repair of CPT-induced DNA damage by homologous recombination. The results show that CPT sensitivity and the global cell cycle progression in cells exposed to CPT is not changed in the absence of Fob1p. Moreover, we have observed in fob1∆ cells treated with CPT that the homologous recombination factor Rad52p still congregates in the shape of foci in the nucleolus, which hosts the rDNA. This suggests that, in the absence of Fob1p, Top1p is still recruited to the rDNA, presumably at sequences other than the rRFB, and its inhibition by CPT leads to recombination events.

Genome-wide mapping of genomic DNA damage: methods and implications

Exposures from the external and internal environments lead to the modification of genomic DNA, which is implicated in the cause of numerous diseases, including cancer, cardiovascular, pulmonary and neurodegenerative diseases, together with ageing. However, the precise mechanism(s) linking the presence of damage, to impact upon cellular function and pathogenesis, is far from clear. Genomic location of specific forms of damage is likely to be highly informative in understanding this process, as the impact of downstream events (e.g. mutation, microsatellite instability, altered methylation and gene expression) on cellular function will be positional—events at key locations will have the greatest impact. However, until recently, methods for assessing DNA damage determined the totality of damage in the genomic location, with no positional information. The technique of “mapping DNA adductomics” describes the molecular approaches that map a variety of forms of DNA damage, to specific locations across the nuclear and mitochondrial genomes. We propose that integrated comparison of this information with other genome-wide data, such as mutational hotspots for specific genotoxins, tumour-specific mutation patterns and chromatin organisation and transcriptional activity in non-cancerous lesions (such as nevi), pre-cancerous conditions (such as polyps) and tumours, will improve our understanding of how environmental toxins lead to cancer. Adopting an analogous approach for non-cancer diseases, including the development of genome-wide assays for other cellular outcomes of DNA damage, will improve our understanding of the role of DNA damage in pathogenesis more generally.

DNA molecular combing-based replication fork directionality profiling

The replication strategy of metazoan genomes is still unclear, mainly because definitive maps of replication origins are missing. High-throughput methods are based on population average and thus may exclusively identify efficient initiation sites, whereas inefficient origins go undetected. Single-molecule analyses of specific loci can detect both common and rare initiation events along the targeted regions. However, these usually concentrate on positioning individual events, which only gives an overview of the replication dynamics. Here, we computed the replication fork directionality (RFD) profiles of two large genes in different transcriptional states in chicken DT40 cells, namely untranscribed and transcribed DMD and CCSER1 expressed at WT levels or overexpressed, by aggregating hundreds of oriented replication tracks detected on individual DNA fibres stretched by molecular combing. These profiles reconstituted RFD domains composed of zones of initiation flanking a zone of termination originally observed in mammalian genomes and were highly consistent with independent population-averaging profiles generated by Okazaki fragment sequencing. Importantly, we demonstrate that inefficient origins do not appear as detectable RFD shifts, explaining why dispersed initiation has remained invisible to population-based assays. Our method can both generate quantitative profiles and identify discrete events, thereby constituting a comprehensive approach to study metazoan genome replication.