Genome-wide detection of DNA double-strand breaks by in-suspension BLISS

CDK12 controls transcription at damaged genes and prevents MYC-induced transcription-replication conflicts

The identification of genes involved in replicative stress is key to understanding cancer evolution and to identify therapeutic targets. Here, we show that CDK12 prevents transcription-replication conflicts (TRCs) and the activation of cytotoxic replicative stress upon deregulation of the MYC oncogene. CDK12 was recruited at damaged genes by PARP-dependent DDR-signaling and elongation-competent RNAPII, to repress transcription. Either loss or chemical inhibition of CDK12 led to DDR-resistant transcription of damaged genes. Loss of CDK12 exacerbated TRCs in MYC-overexpressing cells and led to the accumulation of double-strand DNA breaks, occurring between co-directional early-replicating regions and transcribed genes. Overall, our data demonstrate that CDK12 protects genome integrity by repressing transcription of damaged genes, which is required for proper resolution of DSBs at oncogene-induced TRCs. This provides a rationale that explains both how CDK12 deficiency can promote tandem duplications of early-replicated regions during tumor evolution, and how CDK12 targeting can exacerbate replicative-stress in tumors.

Adaptive and Maladaptive DNA Breaks in Neuronal Physiology and Alzheimer's Disease

DNA strand breaks excessively accumulate in the brains of patients with Alzheimer's disease (AD). While traditionally considered random, deleterious events, neuron activity itself induces DNA breaks, and these "adaptive" breaks help mediate synaptic plasticity and memory formation. Recent studies mapping the brain DNA break landscape reveal that despite a net increase in DNA breaks in ectopic genomic hotspots, adaptive DNA breaks around synaptic genes are lost in AD brains, and this is associated with transcriptomic dysregulation. Additionally, relationships exist between mitochondrial dysfunction, a hallmark of AD, and DNA damage, such that mitochondrial dysfunction may perturb adaptive DNA break formation, while DNA breaks may conversely impair mitochondrial function. A failure of DNA break physiology could, therefore, potentially contribute to AD pathogenesis.

DSB profiles in human spermatozoa highlight the role of TMEJ in the male germline

The male mammalian germline is characterized by substantial chromatin remodeling associated with the transition from histones to protamines during spermatogenesis, followed by the reversal to nucleohistones in the male pronucleus preceding the zygotic genome activation. Both transitions are associated with the extensive formation of DNA double-strand breaks (DSBs), requiring an estimated 5 to 10 million transient DSBs per spermatozoa. Additionally, the high transcription rate in early stages of spermatogenesis leads to transcription-coupled damage preceding meiotic homologous recombination, potentially further contributing to the DSB landscape in mature spermatozoa. Once meiosis is completed, spermatozoa remain haploid and therefore cannot rely on error-free homologous recombination, but instead depend on error-prone classical non-homologous end joining (cNHEJ). This DNA damage/repair-scenario is proposed to be one of the main causes of the observed paternal mutation propensity in human evolution. Recent studies have shown that DSBs in the male pronucleus are repaired by maternally provided Polθ in Caenorhabditis elegans through Polθ-mediated end joining (TMEJ). Additionally, population genetic datasets have revealed a preponderance of TMEJ signatures associated with human variation. Since these signatures are the result of the combined effect of TMEJ and DSB formation in spermatozoa and male pronuclei, we used a BLISS-based protocol to analyze recurrent DSBs in mature human sperm heads as a proxy of the male pronucleus before zygotic chromatin remodeling. The DSBs were found to be enriched in (YR)n short tandem repeats and in evolutionarily young SINEs, reminiscent to patterns observed in murine spermatids, indicating evolutionary hotspots of recurrent DSB formation in mammalian spermatozoa. Additionally, we detected a similar DSB pattern in diploid human IMR90 cells when cNHEJ was selectively inhibited, indicating the significant impact of absent cNHEJ on the sperm DSB landscape. Strikingly, regions associated with most retained histones, and therefore less condensed chromatin, were not strongly enriched with recurrent DSBs. In contrast, the fraction of retained H3K27me3 in the mature spermatozoa displayed a strong association with recurrent DSBs. DSBs in H3K27me3 are associated with a preference for TMEJ over cNHEJ during repair. We hypothesize that the retained H3K27me3 may trigger transgenerational DNA repair by priming maternal Polθ to these regions.

Methods and applications of genome-wide profiling of DNA damage and rare mutations

DNA damage is a threat to genome integrity and can be a cause of many human diseases, owing to either changes in the chemical structure of DNA or conversion of the damage into a mutation, that is, a permanent change in DNA sequence. Determining the exact positions of DNA damage and ensuing mutations in the genome are important for identifying mechanisms of disease aetiology when characteristic mutations are prevalent and probably causative in a particular disease. However, this approach is challenging particularly when levels of DNA damage are low, for example, as a result of chronic exposure to environmental agents or certain endogenous processes, such as the generation of reactive oxygen species. Over the past few years, a comprehensive toolbox of genome-wide methods has been developed for the detection of DNA damage and rare mutations at single-nucleotide resolution in mammalian cells. Here, we review and compare these methods, describe their current applications and discuss future research questions that can now be addressed.

Protocol for mapping physiological DSBs using in-suspension break labeling in situ and sequencing

Physiological double-stranded breaks (DSBs) are a major source of genomic instability. Here, we present a protocol for mapping physiological DSBs by in-suspension break labeling in situ and sequencing (sBLISS) in a single-nucleotide resolution. We describe steps for cell fixation, labeling of DSBs, DNA isolation followed by in vitro transcription (IVT), reverse transcription, and library preparation. sBLISS provides a map of DSBs over the genome and can be used to study the role of different factors in DSB formation. For complete details on the use and execution of this protocol, please refer to Hidmi et al.1.

Uncovering the dynamics of precise repair at CRISPR/Cas9-induced double-strand breaks

CRISPR/Cas9 is widely used for precise mutagenesis through targeted DNA double-strand breaks (DSBs) induction followed by error-prone repair. A better understanding of this process requires measuring the rates of cutting, error-prone, and precise repair, which have remained elusive so far. Here, we present a molecular and computational toolkit for multiplexed quantification of DSB intermediates and repair products by single-molecule sequencing. Using this approach, we characterize the dynamics of DSB induction, processing and repair at endogenous loci along a 72 h time-course in tomato protoplasts. Combining this data with kinetic modeling reveals that indel accumulation is determined by the combined effect of the rates of DSB induction processing of broken ends, and precise versus error repair. In this study, 64-88% of the molecules were cleaved in the three targets analyzed, while indels ranged between 15-41%. Precise repair accounts for most of the gap between cleavage and error repair, representing up to 70% of all repair events. Altogether, this system exposes flux in the DSB repair process, decoupling induction and repair dynamics, and suggesting an essential role of high-fidelity repair in limiting the efficiency of CRISPR-mediated mutagenesis.

ARID1A regulates DNA repair through chromatin organization and its deficiency triggers DNA damage-mediated anti-tumor immune response

AT-rich interaction domain protein 1A (ARID1A), a SWI/SNF chromatin remodeling complex subunit, is frequently mutated across various cancer entities. Loss of ARID1A leads to DNA repair defects. Here, we show that ARID1A plays epigenetic roles to promote both DNA double-strand breaks (DSBs) repair pathways, non-homologous end-joining (NHEJ) and homologous recombination (HR). ARID1A is accumulated at DSBs after DNA damage and regulates chromatin loops formation by recruiting RAD21 and CTCF to DSBs. Simultaneously, ARID1A facilitates transcription silencing at DSBs in transcriptionally active chromatin by recruiting HDAC1 and RSF1 to control the distribution of activating histone marks, chromatin accessibility, and eviction of RNAPII. ARID1A depletion resulted in enhanced accumulation of micronuclei, activation of cGAS-STING pathway, and an increased expression of immunomodulatory cytokines upon ionizing radiation. Furthermore, low ARID1A expression in cancer patients receiving radiotherapy was associated with higher infiltration of several immune cells. The high mutation rate of ARID1A in various cancer types highlights its clinical relevance as a promising biomarker that correlates with the level of immune regulatory cytokines and estimates the levels of tumor-infiltrating immune cells, which can predict the response to the combination of radio- and immunotherapy.

Genomic origin, fragmentomics, and transcriptional properties of long cell-free DNA molecules in human plasma

Recent studies have revealed an unexplored population of long cell-free DNA (cfDNA) molecules in human plasma using long-read sequencing technologies. However, the biological properties of long cfDNA molecules (>500 bp) remain largely unknown. To this end, we have investigated the origins of long cfDNA molecules from different genomic elements. Analysis of plasma cfDNA using long-read sequencing reveals an uneven distribution of long molecules from across the genome. Long cfDNA molecules show overrepresentation in euchromatic regions of the genome, in sharp contrast to short DNA molecules. We observe a stronger relationship between the abundance of long molecules and mRNA gene expression levels, compared with short molecules (Pearson's r = 0.71 vs. -0.14). Moreover, long and short molecules show distinct fragmentation patterns surrounding CpG sites. Leveraging the cleavage preferences surrounding CpG sites, the combined cleavage ratios of long and short molecules can differentiate patients with hepatocellular carcinoma (HCC) from non-HCC subjects (AUC = 0.87). We also investigated knockout mice in which selected nuclease genes had been inactivated in comparison with wild-type mice. The proportion of long molecules originating from transcription start sites are lower in Dffb-deficient mice but higher in Dnase1l3-deficient mice compared with that of wild-type mice. This work thus provides new insights into the biological properties and potential clinical applications of long cfDNA molecules.

TOP1 and R-loops facilitate transcriptional DSBs at hypertranscribed cancer driver genes

DNA double-stranded breaks (DSBs) pose a significant threat to genomic integrity, and their generation during essential cellular processes like transcription remains poorly understood. In this study, we employ several techniques to map DSBs, R-loops, and topoisomerase 1 cleavage complex (TOP1cc) to comprehensively investigate the interplay between transcription, DSBs, topoisomerase 1 (TOP1), and R-loops. Our findings reveal the presence of DSBs at highly expressed genes enriched with TOP1 and R-loops. Remarkably, transcription-associated DSBs at these loci are significantly reduced upon depletion of R-loops and TOP1, uncovering the pivotal roles of TOP1 and R-loops in transcriptional DSB formation. By elucidating the intricate interplay between TOP1cc trapping, R-loops, and DSBs, our study provides insights into the mechanisms underlying transcription-associated genomic instability. Moreover, we establish a link between transcriptional DSBs and early molecular changes driving cancer development, highlighting the distinct etiology and molecular characteristics of driver mutations compared to passenger mutations.

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.

RAD51-mediated R-loop formation acts to repair transcription-associated DNA breaks driving antigenic variation in Trypanosoma brucei

RNA-DNA hybrids are epigenetic features of all genomes that intersect with many processes, including transcription, telomere homeostasis, and centromere function. Increasing evidence suggests that RNA-DNA hybrids can provide two conflicting roles in the maintenance and transmission of genomes: They can be the triggers of DNA damage, leading to genome change, or can aid the DNA repair processes needed to respond to DNA lesions. Evasion of host immunity by African trypanosomes, such as Trypanosoma brucei, relies on targeted recombination of silent Variant Surface Glycoprotein (VSG) genes into a specialized telomeric locus that directs transcription of just one VSG from thousands. How such VSG recombination is targeted and initiated is unclear. Here, we show that a key enzyme of T. brucei homologous recombination, RAD51, interacts with RNA-DNA hybrids. In addition, we show that RNA-DNA hybrids display a genome-wide colocalization with DNA breaks and that this relationship is impaired by mutation of RAD51. Finally, we show that RAD51 acts to repair highly abundant, localised DNA breaks at the single transcribed VSG and that mutation of RAD51 alters RNA-DNA hybrid abundance at 70 bp repeats both around the transcribed VSG and across the silent VSG archive. This work reveals a widespread, generalised role for RNA-DNA hybrids in directing RAD51 activity during recombination and uncovers a specialised application of this interplay during targeted DNA break repair needed for the critical T. brucei immune evasion reaction of antigenic variation.

Let’s not take DNA breaks for granted. The importance of direct detection of DNA breaks for the successful development of DDR inhibitors

Successful development of a drug candidate requires availability of robust methods that enable precise and quantitative assessment of the biological effects exerted by the molecule of interest. In case of DNA Damage Response inhibitors, the most proximal readout of their efficiency is the level of induced DNA damage, usually - DNA breaks. Here we review the methods that are currently used for the assessment of the level of DNA damage, with special attention to their specificity and sensitivity. We also discuss the most common problems and challenges related to the classic IF or IHC methods that indirectly report on the activation of DNA repair mechanisms as the downstream effects of occurrence of the DNA lesions. Finally, we highlight the advent of new tools, such as STRIDE, which have the potential to transform the landscape of DDR functional biomarkers.

A NPAS4-NuA4 complex couples synaptic activity to DNA repair

Neuronal activity is crucial for adaptive circuit remodelling but poses an inherent risk to the stability of the genome across the long lifespan of postmitotic neurons1-5. Whether neurons have acquired specialized genome protection mechanisms that enable them to withstand decades of potentially damaging stimuli during periods of heightened activity is unknown. Here we identify an activity-dependent DNA repair mechanism in which a new form of the NuA4-TIP60 chromatin modifier assembles in activated neurons around the inducible, neuronal-specific transcription factor NPAS4. We purify this complex from the brain and demonstrate its functions in eliciting activity-dependent changes to neuronal transcriptomes and circuitry. By characterizing the landscape of activity-induced DNA double-strand breaks in the brain, we show that NPAS4-NuA4 binds to recurrently damaged regulatory elements and recruits additional DNA repair machinery to stimulate their repair. Gene regulatory elements bound by NPAS4-NuA4 are partially protected against age-dependent accumulation of somatic mutations. Impaired NPAS4-NuA4 signalling leads to a cascade of cellular defects, including dysregulated activity-dependent transcriptional responses, loss of control over neuronal inhibition and genome instability, which all culminate to reduce organismal lifespan. In addition, mutations in several components of the NuA4 complex are reported to lead to neurodevelopmental and autism spectrum disorders. Together, these findings identify a neuronal-specific complex that couples neuronal activity directly to genome preservation, the disruption of which may contribute to developmental disorders, neurodegeneration and ageing.

DNA fragility at the KMT2A/MLL locus: insights from old and new technologies

The Mixed-Lineage Leukaemia (MLL/KMT2A) gene is frequently rearranged in childhood and adult acute leukaemia (AL) and in secondary leukaemias occurring after therapy with DNA topoisomerase targeting anti-cancer agents such as etoposide (t-AL). MLL/KMT2A chromosome translocation break sites in AL patients fall within an 8 kb breakpoint cluster region (BCR). Furthermore, MLL/KMT2A break sites in t-AL frequently occur in a much smaller region, or hotspot, towards the 3' end of the BCR, close to the intron 11/exon 12 boundary. These findings have prompted considerable effort to uncover mechanisms behind the apparent fragility of the BCR and particularly the t-AL hotspot. Recent genome-wide analyses have demonstrated etoposide-induced DNA cleavage within the BCR, and it is tempting to conclude that this cleavage explains the distribution of translocation break sites in t-AL. However, the t-AL hotspot and the centre of the observed preferential DNA cleavage are offset by over 250 nucleotides, suggesting additional factors contribute to the distribution of t-AL break sites. We review these recent genomic datasets along with older experimental results, analysis of TOP2 DNA cleavage site preferences and DNA secondary structure features that may lead to break site selection in t-AL MLL/KMT2A translocations.

Quantitative evaluation of DNA double-strand breaks (DSBs) through single-molecule observation

By adapting the method of single molecular observation for individual DNAs, it will be shown that reliable analysis of double-strand breaks, DSBs, becomes possible for various kinds of damage sources. Single DNA above the size of several-tens kilo base-pairs exhibits the length scale above several μm, indicating that their whole conformation is visible with fluorescence microscopy by adding suitable fluoresce dye to the solution. Various examples of the quantitative evaluation on DSBs are described, together with the evaluation of the protective effects of anti-oxidants.

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.

An atlas of endogenous DNA double-strand breaks arising during human neural cell fate determination

Endogenous DNA double-strand breaks (DSBs) occurring in neural cells have been implicated in the pathogenesis of neurodevelopmental disorders (NDDs). Currently, a genomic map of endogenous DSBs arising during human neurogenesis is missing. Here, we applied in-suspension Breaks Labeling In Situ and Sequencing (sBLISS), RNA-Seq, and Hi-C to chart the genomic landscape of DSBs and relate it to gene expression and genome architecture in 2D cultures of human neuroepithelial stem cells (NES), neural progenitor cells (NPC), and post-mitotic neural cells (NEU). Endogenous DSBs were enriched at the promoter and along the gene body of transcriptionally active genes, at the borders of topologically associating domains (TADs), and around chromatin loop anchors. NDD risk genes harbored significantly more DSBs in comparison to other protein-coding genes, especially in NEU cells. We provide sBLISS, RNA-Seq, and Hi-C datasets for each differentiation stage, and all the scripts needed to reproduce our analyses. Our datasets and tools represent a unique resource that can be harnessed to investigate the role of genome fragility in the pathogenesis of NDDs.

R-loop proximity proteomics identifies a role of DDX41 in transcription-associated genomic instability

Transcription poses a threat to genomic stability through the formation of R-loops that can obstruct progression of replication forks. R-loops are three-stranded nucleic acid structures formed by an RNA-DNA hybrid with a displaced non-template DNA strand. We developed RNA-DNA Proximity Proteomics to map the R-loop proximal proteome of human cells using quantitative mass spectrometry. We implicate different cellular proteins in R-loop regulation and identify a role of the tumor suppressor DDX41 in opposing R-loop and double strand DNA break accumulation in promoters. DDX41 is enriched in promoter regions in vivo, and can unwind RNA-DNA hybrids in vitro. R-loop accumulation upon loss of DDX41 is accompanied with replication stress, an increase in the formation of double strand DNA breaks and transcriptome changes associated with the inflammatory response. Germline loss-of-function mutations in DDX41 lead to predisposition to acute myeloid leukemia in adulthood. We propose that R-loop accumulation and genomic instability-associated inflammatory response may contribute to the development of familial AML with mutated DDX41.

A recurrent chromosomal inversion suffices for driving escape from oncogene-induced senescence via subTAD reorganization

Oncogene-induced senescence (OIS) is an inherent and important tumor suppressor mechanism. However, if not removed timely via immune surveillance, senescent cells also have detrimental effects. Although this has mostly been attributed to the senescence-associated secretory phenotype (SASP) of these cells, we recently proposed that "escape" from the senescent state is another unfavorable outcome. The mechanism underlying this phenomenon remains elusive. Here, we exploit genomic and functional data from a prototypical human epithelial cell model carrying an inducible CDC6 oncogene to identify an early-acquired recurrent chromosomal inversion that harbors a locus encoding the circadian transcription factor BHLHE40. This inversion alone suffices for BHLHE40 activation upon CDC6 induction and driving cell cycle re-entry of senescent cells, and malignant transformation. Ectopic overexpression of BHLHE40 prevented induction of CDC6-triggered senescence. We provide strong evidence in support of replication stress-induced genomic instability being a causative factor underlying "escape" from oncogene-induced senescence.

Somatic Copy Number Alterations in Human Cancers: An Analysis of Publicly Available Data From The Cancer Genome Atlas

Somatic copy number alterations (SCNAs) are a pervasive trait of human cancers that contributes to tumorigenesis by affecting the dosage of multiple genes at the same time. In the past decade, The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) initiatives have generated and made publicly available SCNA genomic profiles from thousands of tumor samples across multiple cancer types. Here, we present a comprehensive analysis of 853,218 SCNAs across 10,729 tumor samples belonging to 32 cancer types using TCGA data. We then discuss current models for how SCNAs likely arise during carcinogenesis and how genomic SCNA profiles can inform clinical practice. Lastly, we highlight open questions in the field of cancer-associated SCNAs.

Systematic overview on the most widespread techniques for inducing and visualizing the DNA double-strand breaks

DNA double-strand breaks (DSBs) are one of the most frequent causes of initiating cancerous malformations, therefore, to reduce the risk, cells have developed sophisticated DNA repair mechanisms. These pathways ensure proper cellular function and genome integrity. However, any alteration or malfunction during DNA repair can influence cellular homeostasis, as improper recognition of the DNA damage or dysregulation of the repair process can lead to genome instability. Several powerful methods have been established to extend our current knowledge in the field of DNA repair. For this reason, in this review, we focus on the methods used to study DSB repair, and we summarize the advantages and disadvantages of the most commonly used techniques currently available for the site-specific induction of DSBs and the subsequent tracking of the repair processes in human cells. We highlight methods that are suitable for site-specific DSB induction (by restriction endonucleases, CRISPR-mediated DSB induction and laser microirradiation) as well as approaches [e.g., fluorescence-, confocal- and super-resolution microscopy, chromatin immunoprecipitation (ChIP), DSB-labeling and sequencing techniques] to visualize and follow the kinetics of DSB repair.

Emerging Technologies for Genome-Wide Profiling of DNA Breakage

Genome instability is associated with myriad human diseases and is a well-known feature of both cancer and neurodegenerative disease. Until recently, the ability to assess DNA damage-the principal driver of genome instability-was limited to relatively imprecise methods or restricted to studying predefined genomic regions. Recently, new techniques for detecting DNA double strand breaks (DSBs) and single strand breaks (SSBs) with next-generation sequencing on a genome-wide scale with single nucleotide resolution have emerged. With these new tools, efforts are underway to define the "breakome" in normal aging and disease. Here, we compare the relative strengths and weaknesses of these technologies and their potential application to studying neurodegenerative diseases.