qDRIP: a method to quantitatively assess RNA-DNA hybrid formation genome-wide

Base-pair resolution reveals clustered R-loops and DNA damage-susceptible R-loops

R-loops are pervasive triplex nucleic acid structures across diverse organisms, yet their biological functions remain incompletely understood. Here, we develop R-loop identification assisted by nucleases and sequencing (RIAN-seq), a nuclease-assisted, antibody-free sequencing technology, to map R-loops at base-pair resolution. By digesting single-stranded RNA (ssRNA), single-stranded DNA (ssDNA), and double-stranded DNA (dsDNA) with nuclease P1, T5 exonuclease, and lambda exonuclease while preserving RNA:DNA hybrids, RIAN-seq achieves unprecedented precision in identifying the position and size of R-loops, detecting an order of magnitude more R-loops than existing methods. Approximately 50% of RNA:DNA hybrids span between 60 and 130 bp, with many forming previously undetectable clusters. Clustered R-loops at promoters recruit zinc-finger proteins VEZF1 and SP5, enhancing transcription in a number-dependent manner and resisting transcriptional perturbation. Conversely, R-loops featuring the Y(C/T)M(A/C)CAG motif at both ends contribute to DNA damage, a phenomenon conserved from yeast to mammalian cells. Our findings reveal a dual role for R-loops: clustered R-loops promote gene expression, while YMCAG-associated R-loops compromise genome stability.

DNA lesions can frequently precede DNA:RNA hybrid accumulation

While DNA:RNA hybrids contribute to multiple genomic transactions, their unscheduled formation is a recognized source of DNA lesions. Here, through a suite of systematic screens, we rather observed that a wide range of yeast mutant situations primarily triggering DNA damage actually leads to hybrid accumulation. Focusing on Okazaki fragment processing, we establish that genic hybrids can actually form as a consequence of replication-born discontinuities such as unprocessed flaps or unligated Okazaki fragments. Strikingly, such "post-lesion" DNA:RNA hybrids neither detectably contribute to genetic instability, nor disturb gene expression, as opposed to "pre-lesion" hybrids formed upon defective mRNA biogenesis, e.g., in THO complex mutants. Post-lesion hybrids similarly arise in distinct genomic instability situations, triggered by pharmacological or genetic manipulation of DNA-dependent processes, both in yeast and human cells. Altogether, our data establish that the accumulation of transcription-born DNA:RNA hybrids can occur as a consequence of various types of natural or pathological DNA lesions, yet do not necessarily aggravate their genotoxicity.

Transcription-replication conflicts drive R-loop-dependent nucleosome eviction and require DOT1L activity for transcription recovery

Progressing transcription and replication machineries profoundly impact their underlying chromatin template. Consequently, transcription-replication conflict (TRC) sites are vulnerable to chromatin and epigenome alterations, provoking genome instability. Here, we engineered an inducible TRC reporter system using a genome-integrated R-loop-prone sequence and characterized the dynamic changes of the local chromatin structure inflicted by TRCs, leading to reduced nucleosome occupancy and replication fork blockage. Strikingly, inducing a small number of TRCs on the genome results in a measurable global replication stress response. Furthermore, we find a TRC-dependent increase in H3K79 methylation specifically at the R-loop forming TRC site. Accordingly, inhibition of the H3K79 methyltransferase DOT1L leads to reduced transcriptional output and an exacerbated DNA damage response, suggesting that deposition of this mark is required for effective transcription recovery and resolution of TRCs. Our work shows the molecular dynamics and reveals a specific epigenetic modifier bookmarking TRC sites, relevant to cancer and other diseases.

Looping forward: exploring R-loop processing and therapeutic potential

Recently, there has been increasing interest in the complex relationship between transcription and genome stability, with specific attention directed toward the physiological significance of molecular structures known as R-loops. These structures arise when an RNA strand invades into the DNA duplex, and their formation is involved in a wide range of regulatory functions affecting gene expression, DNA repair processes or cell homeostasis. The persistent presence of R-loops, if not effectively removed, contributes to genome instability, underscoring the significance of the factors responsible for their resolution and modification. In this review, we provide a comprehensive overview of how R-loop processing can drive either a beneficial or a harmful outcome. Additionally, we explore the potential for manipulating such structures to devise rationalized therapeutic strategies targeting the aberrant accumulation of R-loops.

Clues to transcription/replication collision-induced DNA damage: it was RNAP, in the chromosome, with the fork

DNA replication and RNA transcription processes compete for the same DNA template and, thus, frequently collide. These transcription-replication collisions are thought to lead to genomic instability, which places a selective pressure on organisms to avoid them. Here, we review the predisposing causes, molecular mechanisms, and downstream consequences of transcription-replication collisions (TRCs) with a strong emphasis on prokaryotic model systems, before contrasting prokaryotic findings with cases in eukaryotic systems. Current research points to genomic structure as the primary determinant of steady-state TRC levels and RNA polymerase regulation as the primary inducer of excess TRCs. We review the proposed mechanisms of TRC-induced DNA damage, attempting to clarify their mechanistic requirements. Finally, we discuss what drives genomes to select against TRCs.

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, may 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.

NAT10 resolves harmful nucleolar R-loops depending on its helicase domain and acetylation of DDX21

BACKGROUND

Aberrant accumulation of R-loops leads to DNA damage, genome instability and even cell death. Therefore, the timely removal of harmful R-loops is essential for the maintenance of genome integrity. Nucleolar R-loops occupy up to 50% of cellular R-loops due to the frequent activation of Pol I transcription. However, the mechanisms involved in the nucleolar R-loop resolution remain elusive. The nucleolar acetyltransferase NAT10 harbors a putative RecD helicase domain (RHD), however, if NAT10 acts in the R-loop resolution is still unknown.

METHODS

NAT10 knockdown cell lines were constructed using CRISPR/Cas9 technology and short hairpin RNA targeting NAT10 mRNA, respectively. The level of R-loops was detected by immunofluorescent staining combined with RNase H treatment. The helicase activity of NAT10 or DDX21 was determined by in vitro helicase experiment. The interaction between NAT10 and DDX21 was verified by co-immunoprecipitation, immunofluorescent staining and GST pull-down experiments. Acetylation sites of DDX21 by NAT10 were analyzed by mass spectrometry. NAT10 knockdown-induced DNA damage was evaluated by immunofluorescent staining and Western blot detecting γH2AX.

RESULTS

Depletion of NAT10 led to the accumulation of nucleolar R-loops. NAT10 resolves R-loops through an RHD in vitro and in cells. However, Flag-NAT10 ∆RHD mutant still partially reduced R-loop levels in the NAT10-depleted cells, suggesting that NAT10 might resolve R-loops through additional pathways. Further, the acetyltransferase activity of NAT10 is required for the nucleolar R-loop resolution. NAT10 acetylates DDX21 at K236 and K573 to enhance the helicase activity of DDX21 to unwind nucleolar R-loops. The helicase activity of DDX21 significantly decreased by Flag-DDX21 2KR and increased by Flag-DDX21 2KQ in cells and in vitro. Consequently, NAT10 depletion-induced nucleolar R-loop accumulation led to DNA damage, which was rescued by co-expression of Flag-DDX21 2KQ and Flag-NAT10 G641E, demonstrating that NAT10 resolves nucleolar R-loops through bipartite pathways.

CONCLUSION

We demonstrate that NAT10 is a novel R-loop resolvase and it resolves nucleolar R-loops depending on its helicase activity and acetylation of DDX21. The cooperation of NAT10 and DDX21 provides comprehensive insights into the nucleolar R-loop resolution for maintaining genome stability.

Regulation of R-Loops in DNA Tumor Viruses

R-loops are triple-stranded nucleic acid structures that occur when newly synthesized single-stranded RNA anneals to duplex DNA upon the collision of replication forks with transcription complexes. These RNA-DNA hybrids facilitate several transcriptional processes in the cell and have been described extensively in the literature. Recently, evidence has emerged that R-loops are key regulators of DNA tumor virus transcription and the replication of their lifecycle. Studies have demonstrated that R-loops on the Human Papillomavirus (HPV) genome must be resolved to maintain genome maintenance and avoid viral integration, a hallmark of HPV cancers. For Epstein-Barr virus (EBV), R-loops are formed at the oriLyt to establish lytic replication. Structural maintenance of chromosome proteins 5/6 (SMC5/6) bind to these viral R-loops to repress EBV lytic replication. Most viruses in the herpesvirales order, such as KSHV, contain R-loop-forming sequences. In this perspective, we will describe the current, although limited, literature demonstrating the importance of RNA-DNA hybrids to regulate DNA virus transcription. We will also detail potential new areas of R-loop research and how these viruses can be used as tools to study the growing field of R-loops.

Simultaneous Visualization of R-Loops/RNA:DNA Hybrids and Replication Forks in a DNA Combing Assay

R-loops, structures that play a crucial role in various biological processes, are integral to gene expression, the maintenance of genome stability, and the formation of epigenomic signatures. When these R-loops are deregulated, they can contribute to the development of serious health conditions, including cancer and neurodegenerative diseases. The detection of R-loops is a complex process that involves several approaches. These include S9.6 antibody- or RNAse H-based immunoprecipitation, non-denaturing bisulfite footprinting, gel electrophoresis, and electron microscopy. Each of these methods offers unique insights into the nature and behavior of R-loops. In our study, we introduce a novel protocol that has been developed based on a single-molecule DNA combing assay. This innovative approach allows for the direct and simultaneous visualization of RNA:DNA hybrids and replication forks, providing a more comprehensive understanding of these structures. Our findings confirm the transcriptional origin of the hybrids, adding to the body of knowledge about their formation. Furthermore, we demonstrate that these hybrids have an inhibitory effect on the progression of replication forks, highlighting their potential impact on DNA replication and cellular function.

Deep learning-enhanced R-loop prediction provides mechanistic implications for repeat expansion diseases

R-loops play diverse functional roles, but controversial genomic localization of R-loops have emerged from experimental approaches, posing significant challenges for R-loop research. The development and application of an accurate computational tool for studying human R-loops remains an unmet need. Here, we introduce DeepER, a deep learning-enhanced R-loop prediction tool. DeepER showcases outstanding performance compared to existing tools, facilitating accurate genome-wide annotation of R-loops and a deeper understanding of the position- and context-dependent effects of nucleotide composition on R-loop formation. DeepER also unveils a strong association between certain tandem repeats and R-loop formation, opening a new avenue for understanding the mechanisms underlying some repeat expansion diseases. To facilitate broader utilization, we have developed a user-friendly web server as an integral component of R-loopBase. We anticipate that DeepER will find extensive applications in the field of R-loop research.

ULI-ssDRIP-seq revealed R-loop dynamics during vertebrate early embryogenesis

R-loop, a chromatin structure containing one RNA:DNA hybrid and one unpaired single-stranded DNA, plays multiple biological roles. However, due to technical limitations, the landscapes and potential functions of R-loops during embryogenesis remain elusive. Here, we developed a quantitative and high-resolution ultra-low input R-loop profiling method, named ULI-ssDRIP-seq, which can map global R-loops with as few as 1000 cells. By using ULI-ssDRIP-seq, we reveal the R-loop dynamics in the zebrafish from gametes to early embryos. In oocytes, the R-loop level is relatively low in most regions of the nuclear genome, except maternal-inherited rDNA and mitochondrial genome. The correlation between R-loop and CG methylation dynamics during early development is relatively weak. Furthermore, either up- or down-regulation of global R-loops by knockdown or overexpression of RNase H1 causes a delay of embryonic development with dramatic expression changes in zygotic and maternal genes. This study provides comprehensive R-loop landscapes during early vertebrate embryogenesis and demonstrates the implication of R-loops in embryonic development.

Unraveling the complexity: Advanced methods in analyzing DNA, RNA, and protein interactions

Exploring the intricate interplay within and between nucleic acids, as well as their interactions with proteins, holds pivotal significance in unraveling the molecular complexities steering cancer initiation and progression. To investigate these interactions, a diverse array of highly specific and sensitive molecular techniques has been developed. The selection of a particular technique depends on the specific nature of the interactions. Typically, researchers employ an amalgamation of these different techniques to obtain a comprehensive and holistic understanding of inter- and intramolecular interactions involving DNA-DNA, RNA-RNA, DNA-RNA, or protein-DNA/RNA. Examining nucleic acid conformation reveals alternative secondary structures beyond conventional ones that have implications for cancer pathways. Mutational hotspots in cancer often lie within sequences prone to adopting these alternative structures, highlighting the importance of investigating intra-genomic and intra-transcriptomic interactions, especially in the context of mutations, to deepen our understanding of oncology. Beyond these intramolecular interactions, the interplay between DNA and RNA leads to formations like DNA:RNA hybrids (known as R-loops) or even DNA:DNA:RNA triplex structures, both influencing biological processes that ultimately impact cancer. Protein-nucleic acid interactions are intrinsic cellular phenomena crucial in both normal and pathological conditions. In particular, genetic mutations or single amino acid variations can alter a protein's structure, function, and binding affinity, thus influencing cancer progression. It is thus, imperative to understand the differences between wild-type (WT) and mutated (MT) genes, transcripts, and proteins. The review aims to summarize the frequently employed methods and techniques for investigating interactions involving nucleic acids and proteins, highlighting recent advancements and diverse adaptations of each technique.

Isolation and molecular characterization of a novel relapsing fever group Borrelia from the white-eared opossum Didelphis albiventris in Brazil

This study aimed to detect, isolate and to characterize by molecular methods a relapsing fever group (RFG) Borrelia in white-eared opossums (Didelphis albiventris) from Brazil. During 2015-2018, when opossums (Didelphis spp.) were captured in six municipalities of the state of São Paulo, Brazil, molecular analyses revealed the presence of a novel RFG Borrelia sp. in the blood of seven opossums (Didelphis albiventris), out of 142 sampled opossums (4.9% infection rate). All seven infected opossums were from a single location (Ribeirão Preto municipality). In a subsequent field study in Ribeirão Preto during 2021, two new opossums (D. albiventris) were captured, of which one contained borrelial DNA in its blood. Macerated tissues from this infected opossum were inoculated into laboratory animals (rodents and rabbits) and two big-eared opossums (Didelphis aurita), which had blood samples examined daily via dark-field microscopy. No spirochetes were visualized in the blood of the laboratory animals. Contrastingly, spirochetes were visualized in the blood of the two D. aurita opossums between 12 and 25 days after inoculation. Blood samples from these opossums were used for a multi-locus sequencing typing (MLST) based on six borrelial loci. Phylogenies inferred from MLST genes positioned the sequenced Borrelia genotype into the RFG borreliae clade basally to borreliae of the Asian-African group, forming a monophyletic group with another Brazilian isolate, "Candidatus B. caatinga". Based on this concatenated phylogenetic analysis, which supports that the new borrelial isolate corresponds to a putative new species, we propose the name "Candidatus Borrelia mimona".

Phylogenetic analysis and divergence time estimation of Lycium species in China based on the chloroplast genomes

BACKGROUND

Lycium is an economically and ecologically important genus of shrubs, consisting of approximately 70 species distributed worldwide, 15 of which are located in China. Despite the economic and ecological importance of Lycium, its phylogeny, interspecific relationships, and evolutionary history remain relatively unknown. In this study, we constructed a phylogeny and estimated divergence time based on the chloroplast genomes (CPGs) of 15 species, including subspecies, of the genus Lycium from China.

RESULTS

We sequenced and annotated 15 CPGs in this study. Comparative analysis of these genomes from these Lycium species revealed a typical quadripartite structure, with a total sequence length ranging from 154,890 to 155,677 base pairs (bp). The CPGs was highly conserved and moderately differentiated. Through annotation, we identified a total of 128-132 genes. Analysis of the boundaries of inverted repeat (IR) regions showed consistent positioning: the junctions of the IRb/LSC region were located in rps19 in all Lycium species, IRb/SSC between the ycf1 and ndhF genes, and SSC/IRa within the ycf1 gene. Sequence variation in the SSC region exceeded that in the IR region. We did not detect major expansions or contractions in the IR region or rearrangements or insertions in the CPGs of the 15 Lycium species. Comparative analyses revealed five hotspot regions in the CPG: trnR(UCU), atpF-atpH, ycf3-trnS(GGA), trnS(GGA), and trnL-UAG, which could potentially serve as molecular markers. In addition, phylogenetic tree construction based on the CPG indicated that the 15 Lycium species formed a monophyletic group and were divided into two typical subbranches and three minor branches. Molecular dating suggested that Lycium diverged from its sister genus approximately 17.7 million years ago (Mya) and species diversification within the Lycium species of China primarily occurred during the recent Pliocene epoch.

CONCLUSION

The divergence time estimation presented in this study will facilitate future research on Lycium, aid in species differentiation, and facilitate diverse investigations into this economically and ecologically important genus.

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.

Approaches for Mapping and Analysis of R-loops

R-loops are nucleic acid structures composed of a DNA:RNA hybrid with a displaced non-template single-stranded DNA. Current approaches to identify and map R-loop formation across the genome employ either an antibody targeted against R-loops (S9.6) or a catalytically inactivated form of RNase H1 (dRNH1), a nuclease that can bind and resolve DNA:RNA hybrids via RNA exonuclease activity. This overview article outlines several ways to map R-loops using either methodology, explaining the differences and similarities among the approaches. Bioinformatic analysis of R-loops involves several layers of quality control and processing before visualizing the data. This article provides resources and tools that can be used to accurately process R-loop mapping data and explains the advantages and disadvantages of the resources as compared to one another. © 2024 Wiley Periodicals LLC.

Resolving the polycistronic aftermath: Essential role of topoisomerase IA in preventing R-loops in Leishmania

Kinetoplastid parasites are "living bridges" in the evolution from prokaryotes to higher eukaryotes. The near-intronless genome of the kinetoplastid Leishmania exhibits polycistronic transcription which can facilitate R-loop formation. Therefore, to prevent such DNA-RNA hybrids, Leishmania has retained prokaryotic-like DNA Topoisomerase IA (LdTOPIA) in the course of evolution. LdTOPIA is an essential enzyme that is expressed ubiquitously and is adapted for the compartmentalized eukaryotic form in harboring functional bipartite nuclear localization signals. Although exhibiting greater homology to mycobacterial TOPIA, LdTOPIA could functionally complement the growth lethality of Escherichia coli TOPIA null GyrB ts strain at non-permissive temperatures. Purified LdTOPIA exhibits Mg2+-dependent relaxation of only negatively supercoiled DNA and preference towards single-stranded DNA substrates. LdTOPIA prevents nuclear R-loops as conditional LdTOPIA downregulated parasites exhibit R-loop formation and thereby parasite killing. The clinically used tricyclic antidepressant, norclomipramine could specifically inhibit LdTOPIA and lead to R-loop formation and parasite elimination. This comprehensive study therefore paves an avenue for drug repurposing against Leishmania.

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.

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.

Compartmentalizing damaged DNA: A double-edged sword

A recent publication in Nature by Arnould et al.1 describes a novel chromatin compartment, termed "damaged" or "D compartment," that facilitates the repair of DNA double-strand breaks but also increases the risk of potentially oncogenic translocation formation.

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.

Live-cell imaging unveils distinct R-loop populations with heterogeneous dynamics

We have developed RHINO, a genetically encoded sensor that selectively binds RNA:DNA hybrids enabling live-cell imaging of cellular R-loops. RHINO comprises a tandem array of three copies of the RNA:DNA hybrid binding domain of human RNase H1 connected by optimized linker segments and fused to a fluorescent protein. This tool allows the measurement of R-loop abundance and dynamics in live cells with high specificity and sensitivity. Using RHINO, we provide a kinetic framework for R-loops at nucleoli, telomeres and protein-coding genes. Our findings demonstrate that R-loop dynamics vary significantly across these regions, potentially reflecting the distinct roles R-loops play in different chromosomal contexts. RHINO is a powerful tool for investigating the role of R-loops in cellular processes and their contribution to disease development and progression.

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.

Chromatin compartmentalization regulates the response to DNA damage

The DNA damage response is essential to safeguard genome integrity. Although the contribution of chromatin in DNA repair has been investigated1,2, the contribution of chromosome folding to these processes remains unclear3. Here we report that, after the production of double-stranded breaks (DSBs) in mammalian cells, ATM drives the formation of a new chromatin compartment (D compartment) through the clustering of damaged topologically associating domains, decorated with γH2AX and 53BP1. This compartment forms by a mechanism that is consistent with polymer-polymer phase separation rather than liquid-liquid phase separation. The D compartment arises mostly in G1 phase, is independent of cohesin and is enhanced after pharmacological inhibition of DNA-dependent protein kinase (DNA-PK) or R-loop accumulation. Importantly, R-loop-enriched DNA-damage-responsive genes physically localize to the D compartment, and this contributes to their optimal activation, providing a function for DSB clustering in the DNA damage response. However, DSB-induced chromosome reorganization comes at the expense of an increased rate of translocations, also observed in cancer genomes. Overall, we characterize how DSB-induced compartmentalization orchestrates the DNA damage response and highlight the critical impact of chromosome architecture in genomic instability.

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.

R-Loops at Chromosome Ends: From Formation, Regulation, and Cellular Consequence

Telomeric repeat containing RNA (TERRA) is transcribed from subtelomeric regions to telomeres. TERRA RNA can invade telomeric dsDNA and form telomeric R-loop structures. A growing body of evidence suggests that TERRA-mediated R-loops are critical players in telomere length homeostasis. Here, we will review current knowledge on the regulation of R-loop levels at telomeres. In particular, we will discuss how the central player TERRA and its binding proteins modulate R-loop levels through various mechanisms. We will further provide an overview of the consequences of TERRA-mediated persistent or unscheduled R-loops at telomeres in human ALT cancers and other organisms, with a focus on telomere length regulation after replication interference-induced damage and DNA homologous recombination-mediated repair.

Direct visualization of transcription-replication conflicts reveals post-replicative DNA:RNA hybrids

Transcription-replication collisions (TRCs) are crucial determinants of genome instability. R-loops were linked to head-on TRCs and proposed to obstruct replication fork progression. The underlying mechanisms, however, remained elusive due to the lack of direct visualization and of non-ambiguous research tools. Here, we ascertained the stability of estrogen-induced R-loops on the human genome, visualized them directly by electron microscopy (EM), and measured R-loop frequency and size at the single-molecule level. Combining EM and immuno-labeling on locus-specific head-on TRCs in bacteria, we observed the frequent accumulation of DNA:RNA hybrids behind replication forks. These post-replicative structures are linked to fork slowing and reversal across conflict regions and are distinct from physiological DNA:RNA hybrids at Okazaki fragments. Comet assays on nascent DNA revealed a marked delay in nascent DNA maturation in multiple conditions previously linked to R-loop accumulation. Altogether, our findings suggest that TRC-associated replication interference entails transactions that follow initial R-loop bypass by the replication fork.

R-loops, type I topoisomerases and cancer

R-loops are abundant and dynamic structures ubiquitously present in human cells both in the nuclear and mitochondrial genomes. They form in cis in the wake of transcription complexes and in trans apart from transcription complexes. In this review, we focus on the relationship between R-loops and topoisomerases, and cancer genomics and therapies. We summarize the topological parameters associated with the formation and resolution of R-loops, which absorb and release high levels of genomic negative supercoiling (Sc-). We review the deleterious consequences of excessive R-loops and rationalize how human type IA (TOP3B) and type IB (TOP1) topoisomerases regulate and resolve R-loops in coordination with helicase and RNase H enzymes. We also review the drugs (topoisomerase inhibitors, splicing inhibitors, G4 stabilizing ligands) and cancer predisposing genes (BRCA1/2, transcription, and splicing genes) known to induce R-loops, and whether stabilizing R-loops and thereby inducing genomic damage can be viewed as a strategy for cancer treatment.

The chromatin-associated RNAs in gene regulation and cancer

Eukaryotic genomes are prevalently transcribed into many types of RNAs that translate into proteins or execute gene regulatory functions. Many RNAs associate with chromatin directly or indirectly and are called chromatin-associated RNAs (caRNAs). To date, caRNAs have been found to be involved in gene and transcriptional regulation through multiple mechanisms and have important roles in different types of cancers. In this review, we first present different categories of caRNAs and the modes of interaction between caRNAs and chromatin. We then detail the mechanisms of chromatin-associated nascent RNAs, chromatin-associated noncoding RNAs and emerging m⁶A on caRNAs in transcription and gene regulation. Finally, we discuss the roles of caRNAs in cancer as well as epigenetic and epitranscriptomic mechanisms contributing to cancer, which could provide insights into the relationship between different caRNAs and cancer, as well as tumor treatment and intervention.

Detecting R-Loop Formation Using a Plasmid-Based In Vitro Transcription Assay

R-loops are three-stranded nucleic acid structures that consist of a DNA-RNA hybrid and a displaced single-stranded DNA. Since it was first reported by Ronald Davis and colleagues over 40 years ago, the study of R-loops has become an increasingly expanded area of research. Numerous factors have been identified to modulate the dynamic formation and resolution of R-loops, which are critical for proper controls of gene expression and genome stability. Along the lines of these discoveries, various biochemical and cellular assays have been developed to detect R-loop changes in vitro and in vivo. In this chapter, we describe a protocol for measuring R-loop formation using a plasmid-based in vitro transcription assay. The R-loop formed is then detected and quantified by using gel mobility, antibody staining, and DNA-RNA immunoprecipitation (DRIP)-qPCR assays. Unlike the helicase assay that uses short R-loop substrates, this assay system introduces DNA topology and active transcription as additional variables that impact R-loop formation, thus, more closely recapitulating in vivo situations. Furthermore, this method can be adopted for investigation of cis-elements and trans-acting factors that influence R-loop formation.

Slot Blot Assay for Detection of R Loops

R loops (DNA-RNA hybrid) are three-stranded nucleic acid structures that comprise of template DNA strand hybridized with the nascent RNA leaving the displaced non-template strand. Although a programmed R loop formation can serve as powerful regulators of gene expression, these structures can also turn into major sources of genomic instability and contribute to the development of diseases. Therefore, understanding how cells prevent the deleterious consequences of R loops yet allow R loop formation to participate in various physiological processes will help to understand how their homeostasis is maintained. Detection and quantitative measurements of R loops are critical that largely relied on S9.6 antibody. Immunofluorescence methods are frequently used to localize and quantify R loops in the cell but they require specialized tools for analysis and relatively expensive; therefore, they are not always useful for initial assessments of R loop accumulation. Here, we describe an improved slot blot protocol to detect and estimate R loops and show its sensitivity and specificity using the S9.6 antibody. Since specific factors protecting cells from harmful R loop accumulation are expanding, this protocol can be used to determine R loop accumulation in research and clinical settings.

R-loop-derived cytoplasmic RNA-DNA hybrids activate an immune response

R-loops are RNA-DNA-hybrid-containing nucleic acids with important cellular roles. Deregulation of R-loop dynamics can lead to DNA damage and genome instability1, which has been linked to the action of endonucleases such as XPG2-4. However, the mechanisms and cellular consequences of such processing have remained unclear. Here we identify a new population of RNA-DNA hybrids in the cytoplasm that are R-loop-processing products. When nuclear R-loops were perturbed by depleting the RNA-DNA helicase senataxin (SETX) or the breast cancer gene BRCA1 (refs. 5-7), we observed XPG- and XPF-dependent cytoplasmic hybrid formation. We identify their source as a subset of stable, overlapping nuclear hybrids with a specific nucleotide signature. Cytoplasmic hybrids bind to the pattern recognition receptors cGAS and TLR3 (ref. 8), activating IRF3 and inducing apoptosis. Excised hybrids and an R-loop-induced innate immune response were also observed in SETX-mutated cells from patients with ataxia oculomotor apraxia type 2 (ref. 9) and in BRCA1-mutated cancer cells10. These findings establish RNA-DNA hybrids as immunogenic species that aberrantly accumulate in the cytoplasm after R-loop processing, linking R-loop accumulation to cell death through the innate immune response. Aberrant R-loop processing and subsequent innate immune activation may contribute to many diseases, such as neurodegeneration and cancer.

spKAS-seq reveals R-loop dynamics using low-input materials by detecting single-stranded DNA with strand specificity

R-loops affect transcription and genome stability. Dysregulation of R-loops is related to human diseases. Genome-wide R-loop mapping typically uses the S9.6 antibody or inactive ribonuclease H, both requiring a large number of cells with varying results observed depending on the approach applied. Here, we present strand-specific kethoxal-assisted single-stranded DNA (ssDNA) sequencing (spKAS-seq) to map R-loops by taking advantage of the presence of a ssDNA in the triplex structure. We show that spKAS-seq detects R-loops and their dynamics at coding sequences, enhancers, and other intergenic regions with as few as 50,000 cells. A joint analysis of R-loops and chromatin-bound RNA binding proteins (RBPs) suggested that R-loops can be RBP binding hotspots on the chromatin.

A common transcriptional mechanism involving R-loop and RNA abasic site regulates an enhancer RNA of APOE

RNA is modified by hundreds of chemical reactions and folds into innumerable shapes. However, the regulatory role of RNA sequence and structure and how dysregulation leads to diseases remain largely unknown. Here, we uncovered a mechanism where RNA abasic sites in R-loops regulate transcription by pausing RNA polymerase II. We found an enhancer RNA, AANCR, that regulates the transcription and expression of apolipoprotein E (APOE). In some human cells such as fibroblasts, AANCR is folded into an R-loop and modified by N-glycosidic cleavage; in this form, AANCR is a partially transcribed nonfunctional enhancer and APOE is not expressed. In contrast, in other cell types including hepatocytes and under stress, AANCR does not form a stable R-loop as its sequence is not modified, so it is transcribed into a full-length enhancer that promotes APOE expression. DNA sequence variants in AANCR are associated significantly with APOE expression and Alzheimer's Disease, thus AANCR is a modifier of Alzheimer's Disease. Besides AANCR, thousands of noncoding RNAs are regulated by abasic sites in R-loops. Together our data reveal the essentiality of the folding and modification of RNA in cellular regulation and demonstrate that dysregulation underlies common complex diseases such as Alzheimer's disease.

Pathophysiological Role and Diagnostic Potential of R-Loops in Cancer and Beyond

R-loops are DNA-RNA hybrids that play multifunctional roles in gene regulation, including replication, transcription, transcription-replication collision, epigenetics, and preserving the integrity of the genome. The aberrant formation and accumulation of unscheduled R-loops can disrupt gene expression and damage DNA, thereby causing genome instability. Recent links between unscheduled R-loop accumulation and the abundance of proteins that modulate R-loop biogenesis have been associated with numerous human diseases, including various cancers. Although R-loops are not necessarily causative for all disease entities described to date, they can perpetuate and even exacerbate the initially disease-eliciting pathophysiology, making them structures of interest for molecular diagnostics. In this review, we discuss the (patho) physiological role of R-loops in health and disease, their surprising diagnostic potential, and state-of-the-art techniques for their detection.

An updated overview of experimental and computational approaches to identify non-canonical DNA/RNA structures with emphasis on G-quadruplexes and R-loops

Multiple types of non-canonical nucleic acid structures play essential roles in DNA recombination and replication, transcription, and genomic instability and have been associated with several human diseases. Thus, an increasing number of experimental and bioinformatics methods have been developed to identify these structures. To date, most reviews have focused on the features of non-canonical DNA/RNA structure formation, experimental approaches to mapping these structures, and the association of these structures with diseases. In addition, two reviews of computational algorithms for the prediction of non-canonical nucleic acid structures have been published. One of these reviews focused only on computational approaches for G4 detection until 2020. The other mainly summarized the computational tools for predicting cruciform, H-DNA and Z-DNA, in which the algorithms discussed were published before 2012. Since then, several experimental and computational methods have been developed. However, a systematic review including the conformation, sequencing mapping methods and computational prediction strategies for these structures has not yet been published. The purpose of this review is to provide an updated overview of conformation, current sequencing technologies and computational identification methods for non-canonical nucleic acid structures, as well as their strengths and weaknesses. We expect that this review will aid in understanding how these structures are characterised and how they contribute to related biological processes and diseases.

TERRA and RAD51AP1 promote alternative lengthening of telomeres through an R- to D-loop switch

Alternative lengthening of telomeres (ALT), a telomerase-independent process maintaining telomeres, is mediated by break-induced replication (BIR). RAD52 promotes ALT by facilitating D-loop formation, but ALT also occurs through a RAD52-independent BIR pathway. Here, we show that the telomere non-coding RNA TERRA forms dynamic telomeric R-loops and contributes to ALT activity in RAD52 knockout cells. TERRA forms R-loops in vitro and at telomeres in a RAD51AP1-dependent manner. The formation of R-loops by TERRA increases G-quadruplexes (G4s) at telomeres. G4 stabilization enhances ALT even when TERRA is depleted, suggesting that G4s act downstream of R-loops to promote BIR. In vitro, the telomeric R-loops assembled by TERRA and RAD51AP1 generate G4s, which persist after R-loop resolution and allow formation of telomeric D-loops without RAD52. Thus, the dynamic telomeric R-loops formed by TERRA and RAD51AP1 enable the RAD52-independent ALT pathway, and G4s orchestrate an R- to D-loop switch at telomeres to stimulate BIR.

Protocol to use RNaseH1-based CRISPR to modulate locus-associated R-loops

Modulating R-loop triplex nucleic acid structures reveals their roles across the genome. However, common approaches cannot ascribe functions to R-loops in a locus-associated manner. This protocol presents the use of a locus-associated R-loop-modulating system (dubbed LasR), which employs an inducible RNaseH1-EGFP-dCas9 chimaera. We detail the in silico design of sgRNAs and their transfection with the chimaera, and outline steps confirming RNaseH1-EGFP-dCas9 expression, localization, locus-targeted association, and R-loop modulation in cis or trans using immunoblotting, microscopy, and chromatin and DNA-RNA immunoprecipitation. For complete details on the use and execution of this protocol, please refer to Abraham et al. (2020).

Mitotic DNA synthesis is caused by transcription-replication conflicts in BRCA2-deficient cells

Aberrant replication causes cells lacking BRCA2 to enter mitosis with under-replicated DNA, which activates a repair mechanism known as mitotic DNA synthesis (MiDAS). Here, we identify genome-wide the sites where MiDAS reactions occur when BRCA2 is abrogated. High-resolution profiling revealed that these sites are different from MiDAS at aphidicolin-induced common fragile sites in that they map to genomic regions replicating in the early S-phase, which are close to early-firing replication origins, are highly transcribed, and display R-loop-forming potential. Both transcription inhibition in early S-phase and RNaseH1 overexpression reduced MiDAS in BRCA2-deficient cells, indicating that transcription-replication conflicts (TRCs) and R-loops are the source of MiDAS. Importantly, the MiDAS sites identified in BRCA2-deficient cells also represent hotspots for genomic rearrangements in BRCA2-mutated breast tumors. Thus, our work provides a mechanism for how tumor-predisposing BRCA2 inactivation links transcription-induced DNA damage with mitotic DNA repair to fuel the genomic instability characteristic of cancer cells.

R-Loop Formation in Meiosis: Roles in Meiotic Transcription-Associated DNA Damage

Meiosis is specialized cell division during gametogenesis that produces genetically unique gametes via homologous recombination. Meiotic homologous recombination entails repairing programmed 200–300 DNA double-strand breaks generated during the early prophase. To avoid interference between meiotic gene transcription and homologous recombination, mammalian meiosis is thought to employ a strategy of exclusively transcribing meiotic or post-meiotic genes before their use. Recent studies have shown that R-loops, three-stranded DNA/RNA hybrid nucleotide structures formed during transcription, play a crucial role in transcription and genome integrity. Although our knowledge about the function of R-loops during meiosis is limited, recent findings in mouse models have suggested that they play crucial roles in meiosis. Given that defective formation of an R-loop can cause abnormal transcription and transcription-coupled DNA damage, the precise regulatory network of R-loops may be essential in vivo for the faithful progression of mammalian meiosis and gametogenesis.

A metagenomic portrait of the microbial community responsible for two decades of bioremediation of poly-contaminated groundwater

Biodegradation of pollutants is a sustainable and cost-effective solution to groundwater pollution. Here, we investigate microbial populations involved in biodegradation of poly-contaminants in a pipeline for heavily contaminated groundwater. Groundwater moves from a polluted park to a treatment plant, where an aerated bioreactor effectively removes the contaminants. While the biomass does not settle in the reactor, sediment is collected afterwards and used to seed the new polluted groundwater via a backwash cycle. The pipeline has successfully operated since 1999, but the biological components in the reactor and the contaminated park groundwater have never been described. We sampled seven points along the pipeline, representing the entire remediation process, and characterized the changing microbial communities using genome-resolved metagenomic analysis. We assembled 297 medium- and high-quality metagenome-assembled genome sequences representing on average 46.3% of the total DNA per sample. We found that the communities cluster into two distinct groups, separating the anaerobic communities in the park groundwater from the aerobic communities inside the plant. In the park, the community is dominated by members of the genus Sulfuricurvum, while the plant is dominated by generalists from the order Burkholderiales. Known aromatic compound biodegradation pathways are four times more abundant in the plant-side communities compared to the park-side. Our findings provide a genome-resolved portrait of the microbial community in a highly effective groundwater treatment system that has treated groundwater with a complex contamination profile for two decades.

At the crossroads of RNA biology, genome integrity and cancer

This article is the synthesis of the scientific presentations that took place during two international courses at Institute Curie, one on post-transcriptional gene regulation and the other on genome instability and human disease, that were joined together in their 2021 edition. This joined course brought together the knowledge on RNA metabolism and the maintenance of genome stability.

Walking a tightrope: The complex balancing act of R-loops in genome stability

Although transcription is an essential cellular process, it is paradoxically also a well-recognized cause of genomic instability. R-loops, non-B DNA structures formed when nascent RNA hybridizes to DNA to displace the non-template strand as single-stranded DNA (ssDNA), are partially responsible for this instability. Yet, recent work has begun to elucidate regulatory roles for R-loops in maintaining the genome. In this review, we discuss the cellular contexts in which R-loops contribute to genomic instability, particularly during DNA replication and double-strand break (DSB) repair. We also summarize the evidence that R-loops participate as an intermediate during repair and may influence pathway choice to preserve genomic integrity. Finally, we discuss the immunogenic potential of R-loops and highlight their links to disease should they become pathogenic.

pncCCND1_B Engages an Inhibitory Protein Network to Downregulate CCND1 Expression upon DNA Damage

Promoter-associated noncoding RNAs (pancRNAs) represent a class of noncoding transcripts driven from the promoter region of protein-coding or non-coding genes that operate as cis-acting elements to regulate the expression of the host gene. PancRNAs act by altering the chromatin structure and recruiting transcription regulators. PncCCND1_B is driven by the promoter region of CCND1 and regulates CCND1 expression in Ewing sarcoma through recruitment of a multi-molecular complex composed of the RNA binding protein Sam68 and the DNA/RNA helicase DHX9. In this study, we investigated the regulation of CCND1 expression in Ewing sarcoma cells upon exposure to chemotherapeutic drugs. Pan-inhibitor screening indicated that etoposide, a drug used for Ewing sarcoma treatment, promotes transcription of pncCCND1_B and repression of CCND1 expression. RNA immunoprecipitation experiments showed increased binding of Sam68 to the pncCCND1_B after treatment, despite the significant reduction in DHX9 protein. This effect was associated with the formation of DNA:RNA duplexes at the CCND1 promoter. Furthermore, Sam68 interacted with HDAC1 in etoposide treated cells, thus contributing to chromatin remodeling and epigenetic changes. Interestingly, inhibition of the ATM signaling pathway by KU 55,933 treatment was sufficient to inhibit etoposide-induced Sam68-HDAC1 interaction without rescuing DHX9 expression. In these conditions, the DNA:RNA hybrids persist, thus contributing to the local chromatin inactivation at the CCND1 promoter region. Altogether, our results show an active role of Sam68 in DNA damage signaling and chromatin remodeling on the CCND1 gene by fine-tuning transitions of epigenetic complexes on the CCND1 promoter.

Genome-Wide Analysis of DNA-RNA Hybrids in Yeast by DRIPc-Seq and DRIP-Seq

DNA-RNA hybrids are required for several natural processes in the cell, such as replication and transcription. However, the misregulation of its metabolism is an important source of genetic instability, a hallmark of diseases including cancer. For this reason, genome-wide detection of DNA-RNA hybrids is becoming essential to identify new factors that play a role in its formation or resolution and to understand the global changes in its dynamics because of genetic alterations or chemical treatments. Here, we describe two different immunoprecipitation-based procedures for the genome-wide profiling of DNA-RNA hybrids in the yeast Saccharomyces cerevisiae: DRIP-seq and DRIPc-seq.

Quantitative DNA-RNA Immunoprecipitation Sequencing with Spike-Ins

R-loops are three-stranded nucleic acid structures, comprising an RNA-DNA hybrid and a displaced strand of ssDNA. R-loops have important physiological roles in cells, but deregulation of R-loop dynamics can also have harmful cellular outcomes. The genome-wide mapping of R-loops offers an unbiased approach to study R-loop biology in a wide range of contexts. Here we present a protocol to sequence RNA-DNA hybrids genome-wide with strand-specificity and high resolution. We also include information on how to prepare and incorporate into the workflow appropriate internal spike-in standards which facilitate accurate normalization of the sequencing signal, thereby providing quantitative insights into R-loop formation between different experimental samples.

DNA:RNA Immunoprecipitation from S. pombe Cells for qPCR and Genome-Wide Sequencing

By temporarily distorting the DNA double helix, the moving RNA polymerases can lead to the formation of non-B DNA structures. One of the most abundant and largest non-B DNA structures in the genome is the R-loop, a three-stranded structure forming when the nascent RNA hybridizes with its DNA template, thereby extruding the non-template DNA strand. Growing evidence suggests that at least a subset of R-loops could induce transcription stress and genome instability, although the direct, primary consequences of R-loop formation on the surrounding chromatin are still unclear.To understand the direct impact of R-loops on transcription and genome stability, accurate and quantitative mapping of R-loops is essential. R-loop mapping is commonly achieved using the antibody-based DNA:RNA Immunoprecipitation (DRIP) strategy. While it is reasonably straightforward to obtain robust DRIP enrichments from human cells, this has proved harder in yeast, where DRIP signals are often relatively weak, with a poor signal-to-noise ratio. Although it is unclear whether such weak signals stem from a technical or a biological reality, they make the accurate quantification of DRIP signals all the more important, especially when deep sequencing is used to monitor and quantify the distribution of R-loops genome-wide. Here we propose a DRIP protocol that has been optimized for the mapping and the quantification of R-loops in Schizosaccharomyces pombe but that can also be used in Saccharomyces cerevisiae. As a result, this protocol can be used to generate calibrated DRIP-seq data, where genomic DNA extracted from S. cerevisiae serves as spike-in reference.

INfrastructure for a PHAge REference Database: Identification of Large-Scale Biases in the Current Collection of Cultured Phage Genomes

Background: With advances in sequencing technology and decreasing costs, the number of phage genomes that have been sequenced has increased markedly in the past decade. Materials and Methods: We developed an automated retrieval and analysis system for phage genomes (https://github.com/RyanCook94/inphared) to produce the INfrastructure for a PHAge REference Database (INPHARED) of phage genomes and associated metadata. Results: As of January 2021, 14,244 complete phage genomes have been sequenced. The INPHARED data set is dominated by phages that infect a small number of bacterial genera, with 75% of phages isolated on only 30 bacterial genera. There is further bias, with significantly more lytic phage genomes (∼70%) than temperate (∼30%) within our database. Collectively, this results in ∼54% of temperate phage genomes originating from just three host genera. With much debate on the carriage of antibiotic resistance genes and their potential safety in phage therapy, we searched for putative antibiotic resistance genes. Frequency of antibiotic resistance gene carriage was found to be higher in temperate phages than in lytic phages and again varied with host. Conclusions: Given the bias of currently sequenced phage genomes, we suggest to fully understand phage diversity, efforts should be made to isolate and sequence a larger number of phages, in particular temperate phages, from a greater diversity of hosts.

R-loopBase: a knowledgebase for genome-wide R-loop formation and regulation

R-loops play versatile roles in many physiological and pathological processes, and are of great interest to scientists in multiple fields. However, controversy about their genomic localization and incomplete understanding of their regulatory network raise great challenges for R-loop research. Here, we present R-loopBase (https://rloopbase.nju.edu.cn) to tackle these pressing issues by systematic integration of genomics and literature data. First, based on 107 high-quality genome-wide R-loop mapping datasets generated by 11 different technologies, we present a reference set of human R-loop zones for high-confidence R-loop localization, and spot conservative genomic features associated with R-loop formation. Second, through literature mining and multi-omics analyses, we curate the most comprehensive list of R-loop regulatory proteins and their targeted R-loops in multiple species to date. These efforts help reveal a global regulatory network of R-loop dynamics and its potential links to the development of cancers and neurological diseases. Finally, we integrate billions of functional genomic annotations, and develop interactive interfaces to search, visualize, download and analyze R-loops and R-loop regulators in a well-annotated genomic context. R-loopBase allows all users, including those with little bioinformatics background to utilize these data for their own research. We anticipate R-loopBase will become a one-stop resource for the R-loop community.

R-loop and its functions at the regulatory interfaces between transcription and (epi)genome

R-loop represents a prevalent and specialized chromatin structure critically involved in a wide range of biological processes. In particular, co-transcriptional R-loops, produced often due to RNA polymerase pausing or RNA biogenesis malfunction, can initiate molecular events to context-dependently regulate local gene transcription and crosstalk with chromatin modifications. Cellular "readers" of R-loops are identified, exerting crucial impacts on R-loop homeostasis and gene regulation. Mounting evidence also supports R-loop deregulation as a frequent, sometimes initiating, event during the development of human pathologies, notably cancer and neurological disorder. The purpose of this review is to cover recent advances in understanding the fundamentals of R-loop biology, which have started to unveil complex interplays of R-loops with factors involved in various biological processes such as transcription, RNA processing and epitranscriptomic modification (such as N6-methyladenosine), DNA damage sensing and repair, and epigenetic regulation.