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

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.