High-resolution, strand-specific R-loop mapping via S9.6-based DNA-RNA immunoprecipitation and high-throughput sequencing

Cut&tag: a powerful epigenetic tool for chromatin profiling

Analysis of transcription factors and chromatin modifications at the genome-wide level provides insights into gene regulatory processes, such as transcription, cell differentiation and cellular response. Chromatin immunoprecipitation is the most popular and powerful approach for mapping chromatin, and other enzyme-tethering techniques have recently become available for living cells. Among these, Cleavage Under Targets and Tagmentation (CUT&Tag) is a relatively novel chromatin profiling method that has rapidly gained popularity in the field of epigenetics since 2019. It has also been widely adapted to map chromatin modifications and TFs in different species, illustrating the association of these chromatin epitopes with various physiological and pathological processes. Scalable single-cell CUT&Tag can be combined with distinct platforms to distinguish cellular identity, epigenetic features and even spatial chromatin profiling. In addition, CUT&Tag has been developed as a strategy for joint profiling of the epigenome, transcriptome or proteome on the same sample. In this review, we will mainly consolidate the applications of CUT&Tag and its derivatives on different platforms, give a detailed explanation of the pros and cons of this technique as well as the potential development trends and applications in the future.

R-loop resolution by ARIP4 helicase promotes androgen-mediated transcription induction

Pausing of RNA polymerase II (Pol II) at transcription start sites (TSSs) primes target genes for productive elongation. Coincidentally, DNA double-strand breaks (DSBs) enrich at highly transcribed and Pol II-paused genes, although their interplay remains undefined. Using androgen receptor (AR) signaling as a model, we have uncovered AR-interacting protein 4 (ARIP4) helicase as a driver of androgen-dependent transcription induction. Chromatin immunoprecipitation sequencing analysis revealed that ARIP4 preferentially co-occupies TSSs with paused Pol II. Moreover, we found that ARIP4 complexes with topoisomerase II beta and mediates transient DSB formation upon hormone stimulation. Accordingly, ARIP4 deficiency compromised release of paused Pol II and resulted in R-loop accumulation at a panel of highly transcribed AR target genes. Last, we showed that ARIP4 binds and unwinds R-loops in vitro and that its expression positively correlates with prostate cancer progression. We propose that androgen stimulation triggers ARIP4-mediated unwinding of R-loops at TSSs, enforcing Pol II pause release to effectively drive an androgen-dependent expression program.

DHX9 SUMOylation is required for the suppression of R-loop-associated genome instability

RNA helicase DHX9 is essential for genome stability by resolving aberrant R-loops. However, its regulatory mechanisms remain unclear. Here we show that SUMOylation at lysine 120 (K120) is crucial for DHX9 function. Preventing SUMOylation at K120 leads to R-loop dysregulation, increased DNA damage, and cell death. Cells expressing DHX9 K120R mutant which cannot be SUMOylated are more sensitive to genotoxic agents and this sensitivity is mitigated by RNase H overexpression. Unlike the mutant, wild-type DHX9 interacts with R-loop-associated proteins such as PARP1 and DDX21 via SUMO-interacting motifs. Fusion of SUMO2 to the DHX9 K120R mutant enhances its association with these proteins, reduces R-loop accumulation, and alleviates survival defects of DHX9 K120R. Our findings highlight the critical role of DHX9 SUMOylation in maintaining genome stability by regulating protein interactions necessary for R-loop balance.

Combined HDAC and eIF4A inhibition: A novel epigenetic therapy for pancreatic adenocarcinoma

Pancreatic ductal adenocarcinoma-(PDAC) needs innovative approaches due to its 12% 5-year survival despite current therapies. We show marked sensitivity of pancreatic cancer cells to the combination of a novel eIF4A inhibitor, des-methyl pateamine A (DMPatA), and a histone deacetylase inhibitor, romidepsin, inducing epigenetic reprogramming as an innovative therapeutic strategy. Exploring the mechanistic activity of this combination showed that with a short duration of romidepsin at low doses, robust acetylation persisted up to 48h with the combination, while histone acetylation rapidly faded with monotherapy. This represents an unexpected mechanism of action against PDAC cells that triggers transcriptional overload, metabolic stress, and augmented DNA damage. Structurally different class I HDAC inhibitors exhibit the same hyperacetylation patterns when co-administered with DMPatA, suggesting a class effect. We show efficacy of this combination regimen against tumor growth in a MIA PaCa-2 xenograft model of PDAC with persistent hyperacetylation confirmed in tumor samples.

STATEMENT OF SIGNIFICANCE

Pancreatic ductal adenocarcinoma, a significant clinical challenge, could benefit from the latent potential of epigenetic therapies like HDAC inhibitors-(HDIs), typically limited to hematological malignancies. Our study shows that a synergistic low dose combination of HDIs with an eIF4A-inhibitor in pancreatic cancer models results in marked pre-clinical efficacy, offering a promising new treatment strategy.

HELLS regulates transcription in T-cell lymphomas by reducing unscheduled R-loops and by facilitating RNAPII progression

Chromatin modifiers are emerging as major determinants of many types of cancers, including Anaplastic Large Cell Lymphomas (ALCL), a family of highly heterogeneous T-cell lymphomas for which therapeutic options are still limited. HELLS is a multifunctional chromatin remodeling protein that affects genomic instability by participating in the DNA damage response. Although the transcriptional function of HELLS has been suggested, no clues on how HELLS controls transcription are currently available. In this study, by integrating different multi-omics and functional approaches, we characterized the transcriptional landscape of HELLS in ALCL. We explored the clinical impact of its transcriptional program in a large cohort of 44 patients with ALCL. We demonstrated that HELLS, loaded at the level of intronic regions of target promoters, facilitates RNA Polymerase II (RNAPII) progression along the gene bodies by reducing the persistence of co-transcriptional R-loops and promoting DNA damage resolution. Importantly, selective knockdown of HELLS sensitizes ALCL cells to different chemotherapeutic agents, showing a synergistic effect. Collectively, our work unveils the role of HELLS in acting as a gatekeeper of ALCL genome stability providing a rationale for drug design.

Molecular insights into regulatory RNAs in the cellular machinery

It is apparent that various functional units within the cellular machinery are derived from RNAs. The evolution of sequencing techniques has resulted in significant insights into approaches for transcriptome studies. Organisms utilize RNA to govern cellular systems, and a heterogeneous class of RNAs is involved in regulatory functions. In particular, regulatory RNAs are increasingly recognized to participate in intricately functioning machinery across almost all levels of biological systems. These systems include those mediating chromatin arrangement, transcription, suborganelle stabilization, and posttranscriptional modifications. Any class of RNA exhibiting regulatory activity can be termed a class of regulatory RNA and is typically represented by noncoding RNAs, which constitute a substantial portion of the genome. These RNAs function based on the principle of structural changes through cis and/or trans regulation to facilitate mutual RNA‒RNA, RNA‒DNA, and RNA‒protein interactions. It has not been clearly elucidated whether regulatory RNAs identified through deep sequencing actually function in the anticipated mechanisms. This review addresses the dominant properties of regulatory RNAs at various layers of the cellular machinery and covers regulatory activities, structural dynamics, modifications, associated molecules, and further challenges related to therapeutics and deep learning.

Dormant origin firing promotes head-on transcription-replication conflicts at transcription termination sites in response to BRCA2 deficiency

BRCA2 is a tumor suppressor protein responsible for safeguarding the cellular genome from replication stress and genotoxicity, but the specific mechanism(s) by which this is achieved to prevent early oncogenesis remains unclear. Here, we provide evidence that BRCA2 acts as a critical suppressor of head-on transcription-replication conflicts (HO-TRCs). Using Okazaki-fragment sequencing (Ok-seq) and computational analysis, we identified origins (dormant origins) that are activated near the transcription termination sites (TTS) of highly expressed, long genes in response to replication stress. Dormant origins are a source for HO-TRCs, and drug treatments that inhibit dormant origin firing led to a reduction in HO-TRCs, R-loop formation, and DNA damage. Using super-resolution microscopy, we showed that HO-TRC events track with elongating RNA polymerase II, but not with transcription initiation. Importantly, RNase H2 is recruited to sites of HO-TRCs in a BRCA2-dependent manner to help alleviate toxic R-loops associated with HO-TRCs. Collectively, our results provide a mechanistic basis for how BRCA2 shields against genomic instability by preventing HO-TRCs through both direct and indirect means occurring at predetermined genomic sites based on the pre-cancer transcriptome.

Maintenance of R-loop structures by phosphorylated hTERT preserves genome integrity

As aberrant accumulation of RNA-DNA hybrids (R-loops) causes DNA damage and genome instability, cells express regulators of R-loop structures. Here we report that RNA-dependent RNA polymerase (RdRP) activity of human telomerase reverse transcriptase (hTERT) regulates R-loop formation. We found that the phosphorylated form of hTERT (p-hTERT) exhibits RdRP activity in nuclear speckles both in telomerase-positive cells and telomerase-negative cells with alternative lengthening of telomeres (ALT) activity. The p-hTERT did not associate with telomerase RNA component in nuclear speckles but, instead, with TERRA RNAs to resolve R-loops. Targeting of the TERT gene in ALT cells ablated RdRP activity and impaired tumour growth. Using a genome-scale CRISPR loss-of-function screen, we identified Fanconi anaemia/BRCA genes as synthetic lethal partners of hTERT RdRP. Inactivation of RdRP and Fanconi anaemia/BRCA genes caused accumulation of R-loop structures and DNA damage. These findings indicate that RdRP activity of p-hTERT guards against genome instability by removing R-loop structures.

AAV-mediated genome editing is influenced by the formation of R-loops

Recombinant adeno-associated viral vectors (rAAV) hold an intrinsic ability to stimulate homologous recombination (AAV-HR) and are the most used in clinical settings for in vivo gene therapy. However, rAAVs also integrate throughout the genome. Here, we describe DNA-RNA immunoprecipitation sequencing (DRIP-seq) in murine HEPA1-6 hepatoma cells and whole murine liver to establish the similarities and differences in genomic R-loop formation in a transformed cell line and intact tissue. We show enhanced AAV-HR in mice upon genetic and pharmacological upregulation of R-loops. Selecting the highly expressed Albumin gene as a model locus for genome editing in both in vitro and in vivo experiments showed that the R-loop prone, 3' end of Albumin was efficiently edited by AAV-HR, whereas the upstream R-loop-deficient region did not result in detectable vector integration. In addition, we found a positive correlation between previously reported off-target rAAV integration sites and R-loop enriched genomic regions. Thus, we conclude that high levels of R-loops, present in highly transcribed genes, promote rAAV vector genome integration. These findings may shed light on potential mechanisms for improving the safety and efficacy of genome editing by modulating R-loops and may enhance our ability to predict regions most susceptible to off-target insertional mutagenesis by rAAV vectors.

Linear interaction between replication and transcription shapes DNA break dynamics at recurrent DNA break Clusters

Recurrent DNA break clusters (RDCs) are replication-transcription collision hotspots; many are unique to neural progenitor cells. Through high-resolution replication sequencing and a capture-ligation assay in mouse neural progenitor cells experiencing replication stress, we unravel the replication features dictating RDC location and orientation. Most RDCs occur at the replication forks traversing timing transition regions (TTRs), where sparse replication origins connect unidirectional forks. Leftward-moving forks generate telomere-connected DNA double-strand breaks (DSBs), while rightward-moving forks lead to centromere-connected DSBs. Strand-specific mapping for DNA-bound RNA reveals co-transcriptional dual-strand DNA:RNA hybrids present at a higher density in RDC than in other actively transcribed long genes. In addition, mapping RNA polymerase activity uncovers that head-to-head interactions between replication and transcription machinery result in 60% DSB contribution to the head-on compared to 40% for co-directional. Taken together we reveal TTR as a fragile class and show how the linear interaction between transcription and replication impacts genome stability.

Co-transcriptional R-loops-mediated epigenetic regulation drives growth retardation and docetaxel chemosensitivity enhancement in advanced prostate cancer

R-loops are prevalent three-stranded nucleic acid structures, comprising a DNA-RNA hybrid and a displaced single-stranded DNA, that frequently form during transcription and may be attributed to genomic stability and gene expression regulation. It was recently discovered that RNA modification contributes to maintain the stability of R-loops such as N6-methyladenosine (m6A). Yet, m6A-modified R-loops in regulating gene transcription remains poorly understood. Here, we demonstrated that insulin-like growth factor 2 mRNA-binding proteins (IGF2BPs) recognize R-loops in an m6A-dependent way. Consequently, IGF2BPs overexpression leads to increased overall R-loop levels, cell migration inhibition, and cell growth retardation in prostate cancer (PCa) via precluding the binding of DNA methyltransferase 1(DNMT1) to semaphorin 3 F (SEMA3F) promoters. Moreover, the K homology (KH) domains of IGF2BPs are required for their recognition of m6A-containing R-loops and are required for tumor suppressor functions. Overexpression of SEMA3F markedly enhanced docetaxel chemosensitivity in prostate cancer via regulating Hippo pathway. Our findings point to a distinct R-loop resolution pathway mediated by IGF2BPs, emphasizing the functional importance of IGF2BPs as epigenetic R-loop readers in transcriptional genetic regulation and cancer biology.The manuscript summarizes the new role of N6-methyladenosine in epigenetic regulation, we introduce the distinct R-loop resolution mediated by IGF2BP proteins in an m6A-dependent way, which probably lead to the growth retardation and docetaxel chemotherapy resistance in prostate cancer. Moreover, our findings first emphasized the functional importance of IGF2BPs as epigenetic R-loop readers in transcriptional genetic regulation and cancer biology. In addition, our research provides a novel RBM15/IGF2BPs/DNMT1 trans-omics regulation m6A axis, indicating the new crosstalk between RNA m6A methylation and DNA methylation in prostate cancer.

A circular RNA-gawky-chromatin regulatory axis modulates stress-induced transcription

In response to heavy metal stress, the RNA-binding protein (RBP) gawky translocates into the nucleus and acts as a chromatin-interacting factor to activate the transcription of many stress-responsive genes. However, the upstream regulators of gawky-mediated transcription and their mechanistic details remain unknown. Here, we identified a class of metal-responsive element-containing circRNAs (MRE circRNAs) which specifically interact with gawky during copper stress. Using classic stress-responsive genes as a readout (Drosophila MT), we found that overexpression of MRE circRNAs led to a significant repression in stress-induced transcription. Mechanistically, MRE circRNAs promote the dissociation of gawky from chromatin and increase its aberrant cytoplasmic accumulation, which ultimately impedes the loading of RNA polymerase II to the active gene loci. The MRE motif serves as an important RNA regulon for maintaining the circRNA-gawky interaction, loss of which impaired the inhibitory effects of MRE circRNAs on gawky. Through RNA-seq analyses, we then identified over 500 additional stress-responsive genes whose induced transcription was attenuated upon MRE circRNA overexpression. Finally, we uncovered the physiological relevance of MRE circRNA-mediated regulation in cellular defense against copper overloading. Taken together, this study proposes that the circRNA-RBP-chromatin axis may represent a fundamental regulatory network for gene expression in eukaryotic cells.

Transcription elongation defects link oncogenic SF3B1 mutations to targetable alterations in chromatin landscape

Transcription and splicing of pre-messenger RNA are closely coordinated, but how this functional coupling is disrupted in human diseases remains unexplored. Using isogenic cell lines, patient samples, and a mutant mouse model, we investigated how cancer-associated mutations in SF3B1 alter transcription. We found that these mutations reduce the elongation rate of RNA polymerase II (RNAPII) along gene bodies and its density at promoters. The elongation defect results from disrupted pre-spliceosome assembly due to impaired protein-protein interactions of mutant SF3B1. The decreased promoter-proximal RNAPII density reduces both chromatin accessibility and H3K4me3 marks at promoters. Through an unbiased screen, we identified epigenetic factors in the Sin3/HDAC/H3K4me pathway, which, when modulated, reverse both transcription and chromatin changes. Our findings reveal how splicing factor mutant states behave functionally as epigenetic disorders through impaired transcription-related changes to the chromatin landscape. We also present a rationale for targeting the Sin3/HDAC complex as a therapeutic strategy.

Accelerated DNA replication fork speed due to loss of R-loops in myelodysplastic syndromes with SF3B1 mutation

Myelodysplastic syndromes (MDS) with mutated SF3B1 gene present features including a favourable outcome distinct from MDS with mutations in other splicing factor genes SRSF2 or U2AF1. Molecular bases of these divergences are poorly understood. Here we find that SF3B1-mutated MDS show reduced R-loop formation predominating in gene bodies associated with intron retention reduction, not found in U2AF1- or SRSF2-mutated MDS. Compared to erythroblasts from SRSF2- or U2AF1-mutated patients, SF3B1-mutated erythroblasts exhibit augmented DNA synthesis, accelerated replication forks, and single-stranded DNA exposure upon differentiation. Importantly, histone deacetylase inhibition using vorinostat restores R-loop formation, slows down DNA replication forks and improves SF3B1-mutated erythroblast differentiation. In conclusion, loss of R-loops with associated DNA replication stress represents a hallmark of SF3B1-mutated MDS ineffective erythropoiesis, which could be used as a therapeutic target.

METTL3 drives telomere targeting of TERRA lncRNA through m6A-dependent R-loop formation: a therapeutic target for ALT-positive neuroblastoma

Telomerase-negative tumors maintain telomere length by alternative lengthening of telomeres (ALT), but the underlying mechanism behind ALT remains poorly understood. A proportion of aggressive neuroblastoma (NB), particularly relapsed tumors, are positive for ALT (ALT+), suggesting that a better dissection of the ALT mechanism could lead to novel therapeutic opportunities. TERRA, a long non-coding RNA (lncRNA) derived from telomere ends, localizes to telomeres in a R-loop-dependent manner and plays a crucial role in telomere maintenance. Here we present evidence that RNA modification at the N6 position of internal adenosine (m6A) in TERRA by the methyltransferase METTL3 is essential for telomere maintenance in ALT+ cells, and the loss of TERRA m6A/METTL3 results in telomere damage. We observed that m6A modification is abundant in R-loop enriched TERRA, and the m6A-mediated recruitment of hnRNPA2B1 to TERRA is critical for R-loop formation. Our findings suggest that m6A drives telomere targeting of TERRA via R-loops, and this m6A-mediated R-loop formation could be a widespread mechanism employed by other chromatin-interacting lncRNAs. Furthermore, treatment of ALT+ NB cells with a METTL3 inhibitor resulted in compromised telomere targeting of TERRA and accumulation of DNA damage at telomeres, indicating that METTL3 inhibition may represent a therapeutic approach for ALT+ NB.

Preferential formation of Z-RNA over intercalated motifs in long noncoding RNA

Secondary structure is a principal determinant of lncRNA function, predominantly regarding scaffold formation and interfaces with target molecules. Noncanonical secondary structures that form in nucleic acids have known roles in regulating gene expression and include G-quadruplexes (G4s), intercalated motifs (iMs), and R-loops (RLs). In this paper, we used the computational tools G4-iM Grinder and QmRLFS-finder to predict the formation of each of these structures throughout the lncRNA transcriptome in comparison to protein-coding transcripts. The importance of the predicted structures in lncRNAs in biological contexts was assessed by combining our results with publicly available lncRNA tissue expression data followed by pathway analysis. The formation of predicted G4 (pG4) and iM (piM) structures in select lncRNA sequences was confirmed in vitro using biophysical experiments under near-physiological conditions. We find that the majority of the tested pG4s form highly stable G4 structures, and identify many previously unreported G4s in biologically important lncRNAs. In contrast, none of the piM sequences are able to form iM structures, consistent with the idea that RNA is unable to form stable iMs. Unexpectedly, these C-rich sequences instead form Z-RNA structures, which have not been previously observed in regions containing cytosine repeats and represent an interesting and underexplored target for protein-RNA interactions. Our results highlight the prevalence and potential structure-associated functions of noncanonical secondary structures in lncRNAs, and show G4 and Z-RNA structure formation in many lncRNA sequences for the first time, furthering the understanding of the structure-function relationship in lncRNAs.

Protocol for detection of in vitro R-loop formation using dot blots

Dot-blot analysis is a technique that allows for fast and convenient detection and identification of nucleic acids and proteins. Here, we provide a guide for nucleic acid isolation from eukaryotic cells and sample processing to detect RNA/DNA hybrids. We then provide detailed steps to quantify dot signal intensity. This protocol can be adapted for screening conditions that result in the accumulation of R-loops. For complete details on the use and execution of this protocol, please refer to Smith et al.1.

AKT1 interacts with DHX9 to Mitigate R Loop-Induced Replication Stress in Ovarian Cancer

PARP inhibitor (PARPi)-resistant BRCA-mutant (BRCAm) high-grade serous ovarian cancer (HGSOC) represents a new clinical challenge with unmet therapeutic needs. Here, we performed a quantitative high-throughput drug combination screen that identified the combination of an ATR inhibitor (ATRi) and an AKT inhibitor (AKTi) as an effective treatment strategy for both PARPi-sensitive and PARPi-resistant BRCAm HGSOC. The ATRi and AKTi combination induced DNA damage and R loop-mediated replication stress (RS). Mechanistically, the kinase domain of AKT1 directly interacted with DHX9 and facilitated recruitment of DHX9 to R loops. AKTi increased ATRi-induced R loop-mediated RS by mitigating recruitment of DHX9 to R loops. Moreover, DHX9 was upregulated in tumors from patients with PARPi-resistant BRCAm HGSOC, and high coexpression of DHX9 and AKT1 correlated with worse survival. Together, this study reveals an interaction between AKT1 and DHX9 that facilitates R loop resolution and identifies combining ATRi and AKTi as a rational treatment strategy for BRCAm HGSOC irrespective of PARPi resistance status.

SIGNIFICANCE

Inhibition of the AKT and ATR pathways cooperatively induces R loop-associated replication stress in high-grade serous ovarian cancer, providing rationale to support the clinical development of AKT and ATR inhibitor combinations. See related commentary by Ramanarayanan and Oberdoerffer, p. 793.

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.

Targeting DHX9 Triggers Tumor-Intrinsic Interferon Response and Replication Stress in Small Cell Lung Cancer

Activating innate immunity in cancer cells through cytoplasmic nucleic acid sensing pathways, a phenomenon known as "viral mimicry," has emerged as an effective strategy to convert immunologically "cold" tumors into "hot." Through a curated CRISPR-based screen of RNA helicases, we identified DExD/H-box helicase 9 (DHX9) as a potent repressor of double-stranded RNA (dsRNA) in small cell lung cancers (SCLC). Depletion of DHX9 induced accumulation of cytoplasmic dsRNA and triggered tumor-intrinsic innate immunity. Intriguingly, ablating DHX9 also induced aberrant accumulation of R-loops, which resulted in an increase of DNA damage-derived cytoplasmic DNA and replication stress in SCLCs. In vivo, DHX9 deletion promoted a decrease in tumor growth while inducing a more immunogenic tumor microenvironment, invigorating responsiveness to immune-checkpoint blockade. These findings suggest that DHX9 is a crucial repressor of tumor-intrinsic innate immunity and replication stress, representing a promising target for SCLC and other "cold" tumors in which genomic instability contributes to pathology.

SIGNIFICANCE

One promising strategy to trigger an immune response within tumors and enhance immunotherapy efficacy is by inducing endogenous "virus-mimetic" nucleic acid accumulation. Here, we identify DHX9 as a viral-mimicry-inducing factor involved in the suppression of double-stranded RNAs and R-loops and propose DHX9 as a novel target to enhance antitumor immunity. See related commentary by Chiappinelli, p. 389. This article is featured in Selected Articles from This Issue, p. 384.

Linear Interaction Between Replication and Transcription Shapes DNA Break Dynamics at Recurrent DNA Break Clusters

Recurrent DNA break clusters (RDCs) are replication-transcription collision hotspots; many are unique to neural progenitor cells. Through high-resolution replication sequencing and a capture-ligation assay in mouse neural progenitor cells experiencing replication stress, we unraveled the replication features dictating RDC location and orientation. Most RDCs occur at the replication forks traversing timing transition regions (TTRs), where sparse replication origins connect unidirectional forks. Leftward-moving forks generate telomere-connected DNA double-strand breaks (DSBs), while rightward-moving forks lead to centromere-connected DSBs. Strand-specific mapping for DNA-bound RNA revealed co-transcriptional dual-strand DNA:RNA hybrids present at a higher density in RDC than in other actively transcribed long genes. In addition, mapping RNA polymerase activity revealed that head-to-head interactions between replication and transcription machinery resulted in 60% DSB contribution to the head-on compared to 40% for co-directional. Our findings revealed TTR as a novel fragile class and highlighted how the linear interaction between transcription and replication impacts genome stability.

The putative RNA helicase DDX1 associates with the nuclear RNA exosome and modulates RNA/DNA hybrids (R-loops)

The RNA exosome is a ribonuclease complex that mediates both RNA processing and degradation. This complex is evolutionarily conserved, ubiquitously expressed, and required for fundamental cellular functions, including rRNA processing. The RNA exosome plays roles in regulating gene expression and protecting the genome, including modulating the accumulation of RNA-DNA hybrids (R-loops). The function of the RNA exosome is facilitated by cofactors, such as the RNA helicase MTR4, which binds/remodels RNAs. Recently, missense mutations in RNA exosome subunit genes have been linked to neurological diseases. One possibility to explain why missense mutations in genes encoding RNA exosome subunits lead to neurological diseases is that the complex may interact with cell- or tissue-specific cofactors that are impacted by these changes. To begin addressing this question, we performed immunoprecipitation of the RNA exosome subunit, EXOSC3, in a neuronal cell line (N2A), followed by proteomic analyses to identify novel interactors. We identified the putative RNA helicase, DDX1, as an interactor. DDX1 plays roles in double-strand break repair, rRNA processing, and R-loop modulation. To explore the functional connections between EXOSC3 and DDX1, we examined the interaction following double-strand breaks and analyzed changes in R-loops in N2A cells depleted for EXOSC3 or DDX1 by DNA/RNA immunoprecipitation followed by sequencing. We find that EXOSC3 interaction with DDX1 is decreased in the presence of DNA damage and that loss of EXOSC3 or DDX1 alters R-loops. These results suggest EXOSC3 and DDX1 interact during events of cellular homeostasis and potentially suppress unscrupulous expression of genes promoting neuronal projection.

TDP-43 chronic deficiency leads to dysregulation of transposable elements and gene expression by affecting R-loop and 5hmC crosstalk

TDP-43 is an RNA/DNA-binding protein that forms aggregates in various brain disorders. TDP-43 engages in many aspects of RNA metabolism, but its molecular roles in regulating genes and transposable elements (TEs) have not been extensively explored. Chronic TDP-43 knockdown impairs cell proliferation and cellular responses to DNA damage. At the molecular level, TDP-43 chronic deficiency affects gene expression either locally or distally by concomitantly altering the crosstalk between R-loops and 5-hydroxymethylcytosine (5hmC) in gene bodies and long-range enhancer/promoter interactions. Furthermore, TDP-43 knockdown induces substantial disease-relevant TE activation by influencing their R-loop and 5hmC homeostasis in a locus-specific manner. Together, our findings highlight the genomic roles of TDP-43 in modulating R-loop-5hmC coordination in coding genes, distal regulatory elements, and TEs, presenting a general and broad molecular mechanism underlying the contributions of proteinopathies to the etiology of neurodegenerative disorders.

In vivo detection of DNA secondary structures using permanganate/S1 footprinting with direct adapter ligation and sequencing (PDAL-Seq)

DNA secondary structures are essential elements of the genomic landscape, playing a critical role in regulating various cellular processes. These structures refer to G-quadruplexes, cruciforms, Z-DNA or H-DNA structures, amongst others (collectively called 'non-B DNA'), which DNA molecules can adopt beyond the B conformation. DNA secondary structures have significant biological roles, and their landscape is dynamic and can rearrange due to various factors, including changes in cellular conditions, temperature, and DNA-binding proteins. Understanding this dynamic nature is crucial for unraveling their functions in cellular processes. Detecting DNA secondary structures remains a challenge. Conventional methods, such as gel electrophoresis and chemical probing, have limitations in terms of sensitivity and specificity. Emerging techniques, including next-generation sequencing and single-molecule approaches, offer promise but face challenges since these techniques are mostly limited to only one type of secondary structure. Here we describe an updated version of a technique permanganate/S1 nuclease footprinting, which uses potassium permanganate to trap single-stranded DNA regions as found in many non-B structures, in combination with S1 nuclease digest and adapter ligation to detect genome-wide non-B formation. To overcome technical hurdles, we combined this method with direct adapter ligation and sequencing (PDAL-Seq). Furthermore, we established a user-friendly pipeline available on Galaxy to standardize PDAL-Seq data analysis. This optimized method allows the analysis of many types of DNA secondary structures that form in a living cell and will advance our knowledge of their roles in health and disease.

ATR phosphorylates DHX9 at serine 321 to suppress R-loop accumulation upon genotoxic stress

Aberrant DNA/RNA hybrids (R-loops) formed during transcription and replication disturbances pose threats to genome stability. DHX9 is an RNA helicase involved in R-loop resolution, but how DHX9 is regulated in response to genotoxic stress remains unclear. Here we report that DHX9 is phosphorylated at S321 and S688, with S321 phosphorylation primarily induced by ATR after DNA damage. Phosphorylation of DHX9 at S321 promotes its interaction with γH2AX, BRCA1 and RPA, and is required for its association with R-loops under genotoxic stress. Inhibition of ATR or expression of the non-phosphorylatable DHX9S321A prevents DHX9 from interacting with RPA and R-loops, leading to the accumulation of stress-induced R-loops. Furthermore, depletion of RPA reduces the association between DHX9 and γH2AX, and in vitro binding analysis confirms a direct interaction between DHX9 and RPA. Notably, cells with the non-phosphorylatable DHX9S321A variant exhibit hypersensitivity to genotoxic stress, while those expressing the phosphomimetic DHX9S321D variant prevent R-loop accumulation and display resistance to DNA damage agents. In summary, we uncover a new mechanism by which ATR directly regulates DHX9 through phosphorylation to eliminate stress-induced R-loops.

Aberrant R-loop-mediated immune evasion, cellular communication, and metabolic reprogramming affect cancer progression: a single-cell analysis

Dysregulation of R-loop homeostasis is closely related to various human diseases, including cancer. However, the causality of aberrant R-loops in tumor progression remains unclear. In this study, using single-cell RNA-sequencing datasets from lung adenocarcinoma (LUAD), we constructed an R-loop scoring model to characterize the R-loop state according to the identified R-loop regulators related to EGFR mutations, tissue origins, and TNM stage. We then evaluated the relationships of the R-loop score with the tumor microenvironment (TME) and treatment response. Furthermore, the potential roles of FANCI-mediated R-loops in LUAD were explored using a series of in vitro experiments. Results showed that malignant cells with low R-loop scores displayed glycolysis and epithelial-mesenchymal transition pathway activation and immune escape promotion, thereby hampering the antitumor therapeutic effects. Cell communication analysis suggested that low R-loop scores contributed to T cell exhaustion. We subsequently validated the prognostic value of R-loop scores by using bulk transcriptome datasets across 33 tumor types. The R-loop scoring model well predicted patients' therapeutic response to targeted therapy, chemotherapy, or immunotherapy in 32 independent cohorts. Remarkably, changes in R-loop distribution mediated by FANCI deficiency blocked the activity of Ras signaling pathway, suppressing tumor-cell proliferation and dissemination. In conclusion, this study reveals the underlying molecular mechanism of metabolic reprogramming and T cell exhaustion under R-loop score patterns, and the changes in R-loops mediated by R-loop regulators resulting in tumor progression. Therefore, incorporating anticancer methods based on R-loop or R-loop regulators into the treatment schemes of precision medicine may be beneficial.

Looping out of control: R-loops in transcription-replication conflict

Transcription-replication conflict is a major cause of replication stress that arises when replication forks collide with the transcription machinery. Replication fork stalling at sites of transcription compromises chromosome replication fidelity and can induce DNA damage with potentially deleterious consequences for genome stability and organismal health. The block to DNA replication by the transcription machinery is complex and can involve stalled or elongating RNA polymerases, promoter-bound transcription factor complexes, or DNA topology constraints. In addition, studies over the past two decades have identified co-transcriptional R-loops as a major source for impairment of DNA replication forks at active genes. However, how R-loops impede DNA replication at the molecular level is incompletely understood. Current evidence suggests that RNA:DNA hybrids, DNA secondary structures, stalled RNA polymerases, and condensed chromatin states associated with R-loops contribute to the of fork progression. Moreover, since both R-loops and replication forks are intrinsically asymmetric structures, the outcome of R-loop-replisome collisions is influenced by collision orientation. Collectively, the data suggest that the impact of R-loops on DNA replication is highly dependent on their specific structural composition. Here, we will summarize our current understanding of the molecular basis for R-loop-induced replication fork progression defects.

Complex interplay between FMRP and DHX9 during DNA replication stress

Mutations in, or deficiency of, fragile X messenger ribonucleoprotein (FMRP) is responsible for the Fragile X syndrome (FXS), the most common cause for inherited intellectual disability. FMRP is a nucleocytoplasmic protein, primarily characterized as a translation repressor with poorly understood nuclear function(s). We recently reported that FXS patient cells lacking FMRP sustain higher level of DNA double-strand breaks (DSBs) than normal cells, specifically at sequences prone to forming R-loops, a phenotype further exacerbated by DNA replication stress. Moreover, expression of FMRP, and not an FMRPI304N mutant known to cause FXS, reduced R-loop-associated DSBs. We subsequently reported that recombinant FMRP directly binds R-loops, primarily through the carboxyl terminal intrinsically disordered region. Here, we show that FMRP directly interacts with an RNA helicase, DHX9. This interaction, which is mediated by the amino terminal structured domain of FMRP, is reduced with FMRPI304N. We also show that FMRP inhibits DHX9 helicase activity on RNA:DNA hybrids and the inhibition is also dependent on the amino terminus. Furthermore, the FMRPI304N mutation causes both FMRP and DHX9 to persist on the chromatin in replication stress. These results suggest an antagonistic relationship between FMRP and DHX9 at the chromatin, where their proper interaction leads to dissociation of both proteins from the fully resolved R-loop. We propose that the absence or the loss of function of FMRP leads to persistent presence of DHX9 or both proteins, respectively, on the unresolved R-loop, ultimately leading to DSBs. Our study sheds new light on our understanding of the genome functions of FMRP.

S9.6-based hybrid capture immunoassay for pathogen detection

The detection of pathogens is critical for clinical diagnosis and public health surveillance. Detection is usually done with nucleic acid-based tests (NATs) and rapid antigen tests (e.g., lateral flow assays [LFAs]). Although NATs are more sensitive and specific, their use is often limited in resource-poor settings due to specialized requirements. To address this limitation, we developed a rapid DNA-RNA Hybrid Capture immunoassay (HC) that specifically detects RNA from pathogens. This assay utilizes a unique monoclonal antibody, S9.6, which binds DNA-RNA hybrids. Biotinylated single-stranded DNA probes are hybridized to target RNAs, followed by hybrid capture on streptavidin and detection with S9.6. The HC-ELISA assay can detect as few as 104 RNA molecules that are 2.2 kb in length. We also adapted this assay into a LFA format, where captured Bacillus anthracis rpoB RNA of 3.5 kb length was detectable from a bacterial load equivalent to 107 CFU per 100 mg of mouse tissue using either HC-ELISA or HC-LFA. Importantly, we also demonstrated the versatility of HC by detecting other pathogens, including SARS-CoV-2 and Toxoplasma gondii, showing its potential for broad pathogen detection. Notably, HC does not require amplification of the target nucleic acid and utilizes economical formats like ELISA and LFA, making it suitable for use in sentinel labs for pathogen detection or as a molecular tool in basic research laboratories. Our study highlights the potential of HC as a sensitive and versatile method for RNA-based pathogen detection.

HELZ promotes R loop resolution to facilitate DNA double-strand break repair by homologous recombination

R loops are RNA-DNA hybrid containing structures involved in diverse cellular processes, including DNA double-strand break (DSB) repair. R loop homeostasis involving the formation and resolution of R loops is critical for DSB repair, and its dysregulation leads to genome instability. Here we show that the HELZ helicase promotes R loop resolution to facilitate DSB repair by homologous recombination (HR). HELZ depletion causes hypersensitivity to DSB-inducing agents, and HELZ localizes and binds to DSBs. HELZ depletion further leads to genomic instability in a R loop dependent manner and the accumulation of R loops globally and at DSBs. HELZ binds to R loops in response to DSBs and promotes their resolution, thereby facilitating HR to promote genome integrity. Our findings thus define a role for HELZ in promoting the resolution of R loops critical for DSB repair by HR.

TATDN2 resolution of R-loops is required for survival of BRCA1-mutant cancer cells

BRCA1-deficient cells have increased IRE1 RNase, which degrades multiple microRNAs. Reconstituting expression of one of these, miR-4638-5p, resulted in synthetic lethality in BRCA1-deficient cancer cells. We found that miR-4638-5p represses expression of TATDN2, a poorly characterized member of the TATD nuclease family. We discovered that human TATDN2 has RNA 3' exonuclease and endonuclease activity on double-stranded hairpin RNA structures. Given the cleavage of hairpin RNA by TATDN2, and that BRCA1-deficient cells have difficulty resolving R-loops, we tested whether TATDN2 could resolve R-loops. Using in vitro biochemical reconstitution assays, we found TATDN2 bound to R-loops and degraded the RNA strand but not DNA of multiple forms of R-loops in vitro in a Mg2+-dependent manner. Mutations in amino acids E593 and E705 predicted by Alphafold-2 to chelate an essential Mg2+ cation completely abrogated this R-loop resolution activity. Depleting TATDN2 increased cellular R-loops, DNA damage and chromosomal instability. Loss of TATDN2 resulted in poor replication fork progression in the presence of increased R-loops. Significantly, we found that TATDN2 is essential for survival of BRCA1-deficient cancer cells, but much less so for cognate BRCA1-repleted cancer cells. Thus, we propose that TATDN2 is a novel target for therapy of BRCA1-deficient cancers.

Noncoding function of super enhancer derived mRNA in modulating neighboring gene expression and TAD interaction

Super enhancers are important regulators of gene expression that often overlap with protein-coding genes. However, it is unclear whether the overlapping protein-coding genes and the mRNA derived from them contribute to enhancer activity. Using an erythroid-specific super enhancer that overlaps the Cpox gene as a model, we found that Cpox mRNA has a non-coding function in regulating neighboring protein-coding genes, eRNA expression and TAD interactions. Depletion of Cpox mRNA leads to accumulation of H3K27me3 and release of p300 from the Cpox locus, activating an intra-TAD enhancer and gene expression. Additionally, we identified a head-to-tail interaction between the TAD boundary genes Cpox and Dcbld2 that is facilitated by a novel type of repressive loop anchored by p300 and PRC2/H3K27me3. Our results uncover a regulatory role for mRNA transcribed within a super enhancer context and provide insight into head-to-tail inter-gene interaction in the regulation of gene expression and oncogene activation.

Transcription-Replication Conflicts as a Source of Genome Instability

Transcription and replication both require large macromolecular complexes to act on a DNA template, yet these machineries cannot simultaneously act on the same DNA sequence. Conflicts between the replication and transcription machineries (transcription-replication conflicts, or TRCs) are widespread in both prokaryotes and eukaryotes and have the capacity to both cause DNA damage and compromise complete, faithful replication of the genome. This review will highlight recent studies investigating the genomic locations of TRCs and the mechanisms by which they may be prevented, mitigated, or resolved. We address work from both model organisms and mammalian systems but predominantly focus on multicellular eukaryotes owing to the additional complexities inherent in the coordination of replication and transcription in the context of cell type-specific gene expression and higher-order chromatin organization.

The TDRD3-USP9X complex and MIB1 regulate TOP3B homeostasis and prevent deleterious TOP3B cleavage complexes

TOP3B is stabilized by TDRD3. Hypothesizing that TDRD3 recruits a deubiquitinase, we find that TOP3B interacts with USP9X via TDRD3. Inactivation of USP9X destabilizes TOP3B, and depletion of both TDRD3 and USP9X does not promote further TOP3B ubiquitylation. Additionally, we observe that MIB1 mediates the ubiquitylation and proteasomal degradation of TOP3B by directly interacting with TOP3B independently of TDRD3. Combined depletion of USP9X, TDRD3 and MIB1 causes no additional increase in TOP3B levels compared to MIB1 knockdown alone indicating that the TDRD3-USP9X complex works downstream of MIB1. To comprehend why cells degrade TOP3B in the absence of TDRD3, we measured TOP3Bccs. Lack of TDRD3 increases TOP3Bccs in DNA and RNA, and induced R-loops, γH2AX and growth defect. Biochemical experiments confirm that TDRD3 increases the turnover of TOP3B. Our work provides molecular insights into the mechanisms by which TDRD3 protect cells from deleterious TOP3Bccs which are otherwise removed by TRIM41.

In-Cell Stability Prediction of RNA/DNA Hybrid Duplexes for Designing Oligonucleotides Aimed at Therapeutics

In cells, the formation of RNA/DNA hybrid duplexes regulates gene expression and modification. The environment inside cellular organelles is heterogeneously crowded with high concentrations of biomolecules that affect the structure and stability of RNA/DNA hybrid duplexes. However, the detailed environmental effects remain unclear. Therefore, the mechanistic details of the effect of such molecular crowding were investigated at the molecular level by using thermodynamic and nuclear magnetic resonance analyses, revealing structure-dependent destabilization of the duplexes under crowded conditions. The transition from B- to A-like hybrid duplexes due to a change in conformation of the DNA strand guided by purine-pyrimidine asymmetry significantly increased the hydration number, which resulted in greater destabilization by the addition of cosolutes. By quantifying the individual contributions of environmental factors and the bulk structure of the duplex, we developed a set of parameters that predict the stability of hybrid duplexes with conformational dissimilarities under diverse crowding conditions. A comparison of the effects of environmental conditions in living cells and in vitro crowded solutions on hybrid duplex formation using the Förster resonance energy transfer technique established the applicability of our parameters to living cells. Moreover, our derived parameters can be used to estimate the efficiency of transcriptional inhibition, genome editing, and silencing techniques in cells. This supports the usefulness of our parameters for the visualization of cellular mechanisms of gene expression and the development of nucleic acid-based therapeutics targeting different cells.

Regulation and function of R-loops at repetitive elements

R-loops are atypical, three-stranded nucleic acid structures that contain a stretch of RNA:DNA hybrids and an unpaired, single stranded DNA loop. R-loops are physiological relevant and can act as regulators of gene expression, chromatin structure, DNA damage repair and DNA replication. However, unscheduled and persistent R-loops are mutagenic and can mediate replication-transcription conflicts, leading to DNA damage and genome instability if left unchecked. Detailed transcriptome analysis unveiled that 85% of the human genome, including repetitive regions, hold transcriptional activity. This anticipates that R-loops management plays a central role for the regulation and integrity of genomes. This function is expected to have a particular relevance for repetitive sequences that make up to 75% of the human genome. Here, we review the impact of R-loops on the function and stability of repetitive regions such as centromeres, telomeres, rDNA arrays, transposable elements and triplet repeat expansions and discuss their relevance for associated pathological conditions.

TERRA R-loops connect and protect sister telomeres in mitosis

Resolution of cohesion between sister telomeres in human cells depends on TRF1-mediated recruitment of the polyADP-ribosyltransferase tankyrase to telomeres. In human aged cells, due to insufficient recruitment of TRF1/tankyrase to shortened telomeres, sisters remain cohered in mitosis. This persistent cohesion plays a protective role, but the mechanism by which sisters remain cohered is not well understood. Here we show that telomere repeat-containing RNA (TERRA) holds sister telomeres together through RNA-DNA hybrid (R-loop) structures. We show that a tankyrase-interacting partner, the RNA-binding protein C19orf43, is required for repression of TERRA R-loops. Persistent telomere cohesion in C19orf43-depleted cells is counteracted by RNaseH1, confirming that RNA-DNA hybrids hold sisters together. Consistent with a protective role for persistent telomere cohesion, depletion of C19orf43 in aged cells reduces DNA damage and delays replicative senescence. We propose that the inherent inability of shortened telomeres to recruit R-loop-repressing machinery permits a controlled onset of senescence.

Small molecule targeting of transcription-replication conflict for selective chemotherapy

Targeting transcription replication conflicts, a major source of endogenous DNA double-stranded breaks and genomic instability could have important anticancer therapeutic implications. Proliferating cell nuclear antigen (PCNA) is critical to DNA replication and repair processes. Through a rational drug design approach, we identified a small molecule PCNA inhibitor, AOH1996, which selectively kills cancer cells. AOH1996 enhances the interaction between PCNA and the largest subunit of RNA polymerase II, RPB1, and dissociates PCNA from actively transcribed chromatin regions, while inducing DNA double-stranded breaks in a transcription-dependent manner. Attenuation of RPB1 interaction with PCNA, by a point mutation in RPB1's PCNA-binding region, confers resistance to AOH1996. Orally administrable and metabolically stable, AOH1996 suppresses tumor growth as a monotherapy or as a combination treatment but causes no discernable side effects. Inhibitors of transcription replication conflict resolution may provide a new and unique therapeutic avenue for exploiting this cancer-selective vulnerability.

Mouse HP1γ regulates TRF1 expression and telomere stability

AIMS

Telomeric repeat-containing RNAs are long non-coding RNAs generated from the telomeres. TERRAs are essential for the establishment of heterochromatin marks at telomeres, which serve for the binding of members of the heterochromatin protein 1 (HP1) protein family of epigenetic modifiers involved with chromatin compaction and gene silencing. While HP1γ is enriched on gene bodies of actively transcribed human and mouse genes, it is unclear if its transcriptional role is important for HP1γ function in telomere cohesion and telomere maintenance. We aimed to study the effect of mouse HP1γ on the transcription of telomere factors and molecules that can affect telomere maintenance.

MAIN METHODS

We investigated the telomere function of HP1γ by using HP1γ deficient mouse embryonic fibroblasts (MEFs). We used gene expression analysis of HP1γ deficient MEFs and validated the molecular and mechanistic consequences of HP1γ loss by telomere FISH, immunofluorescence, RT-qPCR and DNA-RNA immunoprecipitation (DRIP).

KEY FINDINGS

Loss of HP1γ in primary MEFs led to a downregulation of various telomere and telomere-accessory transcripts, including the shelterin protein TRF1. Its downregulation is associated with increased telomere replication stress and DNA damage (γH2AX), effects more profound in females. We suggest that the source for the impaired telomere maintenance is a consequence of increased telomeric DNA-RNA hybrids and TERRAs arising at and from mouse chromosomes 18 and X.

SIGNIFICANCE

Our results suggest an important transcriptional control by mouse HP1γ of various telomere factors including TRF1 protein and TERRAs that has profound consequences on telomere stability, with a potential sexually dimorphic nature.

R-Loops in Genome Instability and Cancer

R-loops are unique, three-stranded nucleic acid structures that primarily form when an RNA molecule displaces one DNA strand and anneals to the complementary DNA strand in a double-stranded DNA molecule. R-loop formation can occur during natural processes, such as transcription, in which the nascent RNA molecule remains hybridized with the template DNA strand, while the non-template DNA strand is displaced. However, R-loops can also arise due to many non-natural processes, including DNA damage, dysregulation of RNA degradation pathways, and defects in RNA processing. Despite their prevalence throughout the whole genome, R-loops are predominantly found in actively transcribed gene regions, enabling R-loops to serve seemingly controversial roles. On one hand, the pathological accumulation of R-loops contributes to genome instability, a hallmark of cancer development that plays a role in tumorigenesis, cancer progression, and therapeutic resistance. On the other hand, R-loops play critical roles in regulating essential processes, such as gene expression, chromatin organization, class-switch recombination, mitochondrial DNA replication, and DNA repair. In this review, we summarize discoveries related to the formation, suppression, and removal of R-loops and their influence on genome instability, DNA repair, and oncogenic events. We have also discussed therapeutical opportunities by targeting pathological R-loops.

Promoter R-Loops Recruit U2AF1 to Modulate Its Phase Separation and RNA Splicing

R-loops and guanine quadruplexes (G4s) are secondary structures of nucleic acids that are ubiquitously present in cells and are enriched in promoter regions of genes. By employing a bioinformatic approach based on overlap analysis of transcription factor chromatin immunoprecipitation sequencing (ChIP-seq) data sets, we found that many splicing factors, including U2AF1 whose recognition of the 3' splicing site is crucial for pre-mRNA splicing, exhibit pronounced enrichment at endogenous R-loop- and DNA G4-structure loci in promoter regions of human genes. We also revealed that U2AF1 binds directly to R-loops and DNA G4 structures at a low-nM binding affinity. Additionally, we showed the ability of U2AF1 to undergo phase separation, which could be stimulated by binding with R-loops, but not duplex DNA, RNA/DNA hybrid, DNA G4, or single-stranded RNA. We also demonstrated that U2AF1 binds to promoter R-loops in human cells, and this binding competes with U2AF1's interaction with 3' splicing site and leads to augmented distribution of RNA polymerase II (RNAPII) to promoters over gene bodies, thereby modulating cotranscriptional pre-mRNA splicing. Together, we uncovered a group of candidate proteins that can bind to both R-loops and DNA G4s, revealed the direct and strong interactions of U2AF1 with these nucleic acid structures, and established a biochemical rationale for U2AF1's occupancy in gene promoters. We also unveiled that interaction with R-loops promotes U2AF1's phase separation, and our work suggests that U2AF1 modulates pre-mRNA splicing by regulating RNAPII's partition in transcription initiation versus elongation.

Long noncoding RNA TUG1 regulates smooth muscle cell differentiation via KLF4-myocardin axis

Abdominal aortic aneurysms (AAAs) are asymptomatic vascular diseases that have life-threatening outcomes. Smooth muscle cell (SMC) dysfunction plays an important role in AAA development. The contribution of non-coding genome, specifically the role of long non-coding RNAs (lncRNAs) in SMC dysfunction, is relatively unexplored. We investigated the role of lncRNA TUG1 in SMC dysfunction. To identify potential lncRNAs relevant to SMC functionality, lncRNA profiling was performed in angiotensin-II-treated SMCs. AAA was induced by angiotensin-II treatment in mice. Transcriptional regulation of TUG1 was studied using promoter luciferase and chromatin-immuno-precipitation experiments. Gain-or-loss-of-function experiments were performed in vitro to investigate TUG1-mediated regulation of SMC function. Immunoprecipitation experiments were conducted to elucidate the mechanism underlying TUG1-mediated SMC dysfunction. TUG1 was upregulated in SMCs following angiotensin-II treatment. Similarly, TUG1 levels were elevated in abdominal aorta in a mouse model of angiotensin-II-induced AAA. Further investigations showed that angiotensin-II-induced TUG1 expression could be suppressed by inhibiting Notch-signaling pathway, both in vitro and in mouse AAA model and that TUG1 is a direct transcriptional target of the Notch pathway. In aneurysmal tissues, TUG1 expression was inversely correlated with the expression of SMC contractile genes. Overexpression of TUG1 repressed SMC differentiation in vitro, whereas siRNA/shRNA-mediated TUG1 knockdown showed an opposite effect. Mechanistically, TUG1 interacts with transcriptional repressor KLF4 and facilitates its recruitment to myocardin promoter ultimately leading to the repression of SMC differentiation. In summary, our study uncovers a novel role for the lncRNA TUG1 wherein it modulates SMC differentiation via the KLF4-myocardin axis, which may have potential implications in AAA development.NEW & NOTEWORTHY TUG1 is an angiotensin-II-induced long noncoding RNA that mediates smooth muscle cell (SMC) dysfunction through interaction with transcriptional repressor KLF4.

Binding of the nuclear ribonucleoprotein family member FUS to RNA prevents R-loop RNA:DNA hybrid structures

The protein FUS (FUSed in sarcoma) is a metazoan RNA-binding protein that influences RNA production by all three nuclear polymerases. FUS also binds nascent transcripts, RNA processing factors, RNA polymerases, and transcription machinery. Here, we explored the role of FUS binding interactions for activity during transcription. In vitro run-off transcription assays revealed FUS-enhanced RNA produced by a non-eukaryote polymerase. The activity also reduced the formation of R-loops between RNA products and their DNA template. Analysis by domain mutation and deletion indicated RNA-binding was required for activity. We interpret that FUS binds and sequesters nascent transcripts to prevent R-loops from forming with nearby DNA. DRIP-seq analysis showed that a knockdown of FUS increased R-loop enrichment near expressed genes. Prevention of R-loops by FUS binding to nascent transcripts has the potential to affect transcription by any RNA polymerase, highlighting the broad impact FUS can have on RNA metabolism in cells and disease.

The RNA m5C modification in R-loops as an off switch of Alt-NHEJ

The roles of R-loops and RNA modifications in homologous recombination (HR) and other DNA double-stranded break (DSB) repair pathways remain poorly understood. Here, we find that DNA damage-induced RNA methyl-5-cytosine (m5C) modification in R-loops plays a crucial role to regulate PARP1-mediated poly ADP-ribosylation (PARylation) and the choice of DSB repair pathways at sites of R-loops. Through bisulfite sequencing, we discover that the methyltransferase TRDMT1 preferentially generates m5C after DNA damage in R-loops across the genome. In the absence of m5C, R-loops activate PARP1-mediated PARylation both in vitro and in cells. Concurrently, m5C promotes transcription-coupled HR (TC-HR) while suppressing PARP1-dependent alternative non-homologous end joining (Alt-NHEJ), favoring TC-HR over Alt-NHEJ in transcribed regions as the preferred repair pathway. Importantly, simultaneous disruption of both TC-HR and Alt-NHEJ with TRDMT1 and PARP or Polymerase θ inhibitors prevents alternative DSB repair and exhibits synergistic cytotoxic effects on cancer cells, suggesting an effective strategy to exploit genomic instability in cancer therapy.

KAS-Analyzer: a novel computational framework for exploring KAS-seq data

Motivation

Kethoxal-assisted ssDNA sequencing (KAS-seq) is rapidly gaining popularity as a robust and effective approach to study the nascent dynamics of transcriptionally engaged RNA polymerases through profiling of genome-wide single-stranded DNA (ssDNA). Its latest variant, spKAS-seq, a strand-specific version of KAS-seq, has been developed to map genome-wide R-loop structures by detecting imbalances of ssDNA on two strands. However, user-friendly, open-source computational tools tailored for KAS-seq data are still lacking.

Results

Here, we introduce KAS-Analyzer, the first comprehensive computational framework aimed at streamlining and enhancing the analysis and interpretation of KAS-seq and spKAS-seq data. In addition to standard analyses, KAS-Analyzer offers many novel tools specifically designed for KAS-seq data, including, but not limited to: calculation of transcription-related metrics, identification of single-stranded transcribing (SST) enhancers, high-resolution mapping of R-loops, and differential RNA polymerase activity analysis. We provided a detailed overview of KAS-seq data and its diverse applications through the implementation of KAS-Analyzer. Using the example time-course KAS-seq datasets, we further showcase the robust capabilities of KAS-Analyzer for investigating dynamic transcriptional regulatory programs in response to UVB radiation.

Availability and implementation

KAS-Analyzer is available at https://github.com/Ruitulyu/KAS-Analyzer.

p53-dependent R-loop formation and HPV pathogenesis

R-loops are trimeric RNA: DNA hybrids that are important physiological regulators of transcription; however, their aberrant formation or turnover leads to genomic instability and DNA breaks. High-risk human papillomaviruses (HPV) are the causative agents of genital as well as oropharyngeal cancers and exhibit enhanced amounts of DNA breaks. The levels of R-loops were found to be increased up to 50-fold in cells that maintain high-risk HPV genomes and were readily detected in squamous cell cervical carcinomas in vivo but not in normal cells. The high levels of R-loops in HPV-positive cells were present on both viral and cellular sites together with RNase H1, an enzyme that controls their resolution. Depletion of RNase H1 in HPV-positive cells further increased R-loop levels, resulting in impaired viral transcription and replication along with reduced expression of the DNA repair genes such as FANCD2 and ATR, both of which are necessary for viral functions. Overexpression of RNase H1 decreased total R-loop levels, resulting in a reduction of DNA breaks by over 50%. Furthermore, increased RNase H1 expression blocked viral transcription and replication while enhancing the expression of factors in the innate immune regulatory pathway. This suggests that maintaining elevated R-loop levels is important for the HPV life cycle. The E6 viral oncoprotein was found to be responsible for inducing high levels of R-loops by inhibiting p53's transcriptional activity. Our studies indicate that high R-loop levels are critical for HPV pathogenesis and that this depends on suppressing the p53 pathway.

MEN1 is a regulator of alternative splicing and prevents R-loop-induced genome instability through suppression of RNA polymerase II elongation

The fidelity of alternative splicing (AS) patterns is essential for growth development and cell fate determination. However, the scope of the molecular switches that regulate AS remains largely unexplored. Here we show that MEN1 is a previously unknown splicing regulatory factor. MEN1 deletion resulted in reprogramming of AS patterns in mouse lung tissue and human lung cancer cells, suggesting that MEN1 has a general function in regulating alternative precursor mRNA splicing. MEN1 altered exon skipping and the abundance of mRNA splicing isoforms of certain genes with suboptimal splice sites. Chromatin immunoprecipitation and chromosome walking assays revealed that MEN1 favored the accumulation of RNA polymerase II (Pol II) in regions encoding variant exons. Our data suggest that MEN1 regulates AS by slowing the Pol II elongation rate and that defects in these processes trigger R-loop formation, DNA damage accumulation and genome instability. Furthermore, we identified 28 MEN1-regulated exon-skipping events in lung cancer cells that were closely correlated with survival in patients with lung adenocarcinoma, and MEN1 deficiency sensitized lung cancer cells to splicing inhibitors. Collectively, these findings led to the identification of a novel biological role for menin in maintaining AS homeostasis and link this role to the regulation of cancer cell behavior.

Identification and characterization of topoisomerase III beta poisons

We designed and carried out a high-throughput screen for compounds that trap topoisomerase III beta (TOP3B poisons) by developing a Comparative Cellular Cytotoxicity Screen. We found a bisacridine compound NSC690634 and a thiacyanine compound NSC96932 that preferentially sensitize cell lines expressing TOP3B, indicating that they target TOP3B. These compounds trap TOP3B cleavage complex (TOP3Bcc) in cells and in vitro and predominately act on RNA, leading to high levels of RNA-TOP3Bccs. NSC690634 also leads to enhanced R-loops in a TOP3B-dependent manner. Preliminary structural activity studies show that the lengths of linkers between the two aromatic moieties in each compound are critical; altering the linker length completely abolishes the trapping of TOP3Bccs. Both of our lead compounds share a similar structural motif, which can serve as a base for further modification. They may also serve in anticancer, antiviral, and/or basic research applications.

Integrative modeling of lncRNA-chromatin interaction maps reveals diverse mechanisms of nuclear retention

BACKGROUND

Many long non-coding RNAs, known to be involved in transcriptional regulation, are enriched in the nucleus and interact with chromatin. However, their mechanisms of chromatin interaction and the served cellular functions are poorly understood. We sought to characterize the mechanisms of lncRNA nuclear retention by systematically mapping the sequence and chromatin features that distinguish lncRNA-interacting genomic segments.

RESULTS

We found DNA 5-mer frequencies to be predictive of chromatin interactions for all lncRNAs, suggesting sequence-specificity as a global theme in the interactome. Sequence features representing protein-DNA and protein-RNA binding motifs revealed potential mechanisms for specific lncRNAs. Complementary to these global themes, transcription-related features and DNA-RNA triplex formation potential were noted to be highly predictive for two mutually exclusive sets of lncRNAs. DNA methylation was also noted to be a significant predictor, but only when combined with other epigenomic features.

CONCLUSIONS

Taken together, our statistical findings suggest that a group of lncRNAs interacts with transcriptionally inactive chromatin through triplex formation, whereas another group interacts with transcriptionally active regions and is involved in DNA Damage Response (DDR) through formation of R-loops. Curiously, we observed a strong pattern of enrichment of 5-mers in four potentially interacting entities: lncRNA-bound DNA tiles, lncRNAs, miRNA seed sequences, and repeat elements. This finding points to a broad sequence-based network of interactions that may underlie regulation of fundamental cellular functions. Overall, this study reveals diverse sequence and chromatin features related to lncRNA-chromatin interactions, suggesting potential mechanisms of nuclear retention and regulatory function.

The exoribonuclease XRN2 mediates degradation of the long non-coding telomeric RNA TERRA

The telomeric repeat-containing RNA, TERRA, associates with both telomeric DNA and telomeric proteins, often forming RNA:DNA hybrids (R-loops). TERRA is most abundant in cancer cells utilizing the alternative lengthening of telomeres (ALT) pathway for telomere maintenance, suggesting that persistent TERRA R-loops may contribute to activation of the ALT mechanism. Therefore, we sought to identify the enzyme(s) that regulate TERRA metabolism in mammalian cells. Here, we identify that the 5'-3' exoribonuclease XRN2 regulates the stability of TERRA RNA. Moreover, while stabilization of TERRA alone was insufficient to drive ALT, depletion of XRN2 in ALT-positive cells led to a significant increase in TERRA R-loops and exacerbated ALT activity. Together, our findings highlight XRN2 as a key determinant of TERRA metabolism and telomere stability in cancer cells that rely on the ALT pathway.