GC skew at the 5' and 3' ends of human genes links R-loop formation to epigenetic regulation and transcription termination

Determination of codon pattern and evolutionary forces acting on genes linked to inflammatory bowel disease

Inflammatory bowel disease (IBD) is an inflammatory disorder of the gastrointestinal tract. The present study attempted to understand the codon usage preferences in genes associated with IBD progression. Compositional analysis, codon usage bias (CUB), Relative synonymous codon usage (RSCU), RNA structure, and expression analysis were performed to obtain a comprehensive picture of codon usage in IBD genes. Compositional analysis of 62 IBD-associated genes revealed that G and T are the most and least abundant nucleotides, respectively. ApG, CpA, and TpG dinucleotides were overrepresented or randomly used, while ApC, CpG, GpT, and TpA dinucleotides were either underrepresented or randomly used in genes related to IBD. The codons influencing the codon usage the most in IBD genes were CGC and AGG. A comparison of codon usage between IBD, and pancreatitis (non-IBD inflammatory disease) indicated that only codon CTG codon usage was significantly different between IBD and pancreatitis. At the same time, there were codons ATA, ACA, CGT, CAA, GTA, CCT, ATT, GCT, CGG, TTG, and CAG for whom codon usage was significantly different for IBD and housekeeping gene sets. The results suggest similar codon usage in at least two inflammatory disorders, IBD and pancreatitis. The analysis helps understand the codon biology, factors affecting gene expression of IBD-associated genes, and the evolution of these genes. The study helps reveal the molecular patterns associated with IBD.

M6A-methylated circPOLR2B forms an R-loop and regulates the biological behavior of glioma stem cells through positive feedback loops

Glioma is the most common primary brain tumor, and targeting glioma stem cells (GSCs) has become a key aspect of glioma treatment. In this study, we discovered a molecular network in which circRNA forms an R-loop structure with its parental gene to regulate the biological behavior of GSCs. Genes with abnormal expression in GSCs were screened using RNA-seq and circRNA microarray analyses. The study results showed that high expression of YTHDC1 in GSCs promoted the transportation of N6-methyladenosine (m6A)-modified circPOLR2B from the nucleus to the cytoplasm. Decreased circPOLR2B levels in the nucleus resulted in fewer R-loop structures formed with its parental gene POLR2B. This reduction in R-loop structures relieved the inhibitory effect on POLR2B transcription and upregulated PBX1 expression through alternative polyadenylation (APA) action, thereby promoting the malignant biological behavior of GSCs. Knockdown of YTHDC1, POLR2B, and PBX1 reduced xenograft tumor volume and prolonged the survival of nude mice. The YTHDC1/circPOLR2B/POLR2B/PBX1 axis plays a regulatory role in the biological behavior of GSCs, offering potential targets and novel strategies for the treatment of glioma.

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.

RNA biogenesis and RNA metabolism factors as R-loop suppressors: a hidden role in genome integrity

Genome integrity relies on the accuracy of DNA metabolism, but as appreciated for more than four decades, transcription enhances mutation and recombination frequencies. More recent research provided evidence for a previously unforeseen link between RNA and DNA metabolism, which is often related to the accumulation of DNA-RNA hybrids and R-loops. In addition to physiological roles, R-loops interfere with DNA replication and repair, providing a molecular scenario for the origin of genome instability. Here, we review current knowledge on the multiple RNA factors that prevent or resolve R-loops and consequent transcription-replication conflicts and thus act as modulators of genome dynamics.

Role of senataxin in R-loop-mediated neurodegeneration

Senataxin is an RNA:DNA helicase that plays an important role in the resolution of RNA:DNA hybrids (R-loops) formed during transcription. R-loops are involved in the regulation of biological processes such as immunoglobulin class switching, gene expression and DNA repair. Excessive accumulation of R-loops results in DNA damage and loss of genomic integrity. Senataxin is critical for maintaining optimal levels of R-loops to prevent DNA damage and acts as a genome guardian. Within the nucleus, senataxin interacts with various RNA processing factors and DNA damage response and repair proteins. Senataxin interactors include survival motor neuron and zinc finger protein 1, with whom it co-localizes in sub-nuclear bodies. Despite its ubiquitous expression, mutations in senataxin specifically affect neurons and result in distinct neurodegenerative diseases such as amyotrophic lateral sclerosis type 4 and ataxia with oculomotor apraxia type 2, which are attributed to the gain-of-function and the loss-of-function mutations in senataxin, respectively. In addition, low levels of senataxin (loss-of-function) in spinal muscular atrophy result in the accumulation of R-loops causing DNA damage and motor neuron degeneration. Senataxin may play multiple functions in diverse cellular processes; however, its emerging role in R-loop resolution and maintenance of genomic integrity is gaining attention in the field of neurodegenerative diseases. In this review, we highlight the role of senataxin in R-loop resolution and its potential as a therapeutic target to treat neurodegenerative diseases.

Wiskott Aldrich syndrome protein (WASp)-deficient Th1 cells promote R-loop-driven transcriptional insufficiency and transcription-coupled nucleotide excision repair factor (TC-NER)-driven genome-instability in the pathogenesis of T cell acute lymphoblastic leukemia

BACKGROUND

T-ALL is an aggressive hematological tumor that develops as the result of a multi-step oncogenic process which causes expansion of hematopoietic progenitors that are primed for T cell development to undergo malignant transformation and growth. Even though first-line therapy has a significant response rate, 40% of adult patients and 20% of pediatric patients will relapse. Therefore, there is an unmet need for treatment for relapsed/refractory T-ALL to develop potential targeted therapies.

METHODS

Pediatric T-ALL patient derived T cells were grown under either nonskewingTh0 or Th1-skewing conditions to further process for ChIP-qPCR, RDIP-qPCR and other RT-PCR assays. Endogenous WASp was knocked out using CRISPR-Cas9 and was confirmed using flow cytometry and western blotting. LC-MS/MS was performed to find out proteomic dataset of WASp-interactors generated from Th1-skewed, human primary Th-cells. DNA-damage was assessed by immunofluorescence confocal-imaging and single-cell gel electrophoresis (comet assay). Overexpression of RNaseH1 was also done to restore normal Th1-transcription in WASp-deficient Th1-skewed cells.

RESULTS

We discovered that nuclear-WASp is required for suppressing R-loop production (RNA/DNA-hybrids) at Th1-network genes by ribonucleaseH2 (RNH2) and topoisomerase1. Nuclear-WASp is associated with the factors involved in preventing and dissolving R-loops in Th1 cells. In nuclear- WASp-reduced malignant Th1-cells, R-loops accumulate in vivo and are processed into DNA-breaks by transcription-coupled-nucleotide-excision repair (TC-NER). Several epigenetic modifications were also found to be involved at Th1 gene locus which are responsible for active/repressive marks of particular genes. By demonstrating WASp as a physiologic regulator of programmed versus unprogrammed R-loops, we suggest that the transcriptional role of WASp in vivo extends also to prevent transcription-linked DNA damage during malignancy and through modification of epigenetic dysregulations.

CONCLUSION

Our findings present a provocative possibility of resetting R-loops as a therapeutic intervention to correct both immune deficiency and malignancy in T-cell acute lymphoblastic leukemia patients and a novel role of WASp in the epigenetic regulation of T helper cell differentiation in T-ALL patients, anticipating WASp's requirement for the suppression of T-ALL progression.

Characterizing nascent transcription patterns of PROMPTs, eRNAs, and readthrough transcripts in the ENCODE4 deeply profiled cell lines

Arising as co-products of canonical gene expression, transcription-associated lincRNAs, such as promoter upstream transcripts (PROMPTs), enhancer RNAs (eRNAs), and readthrough (RT) transcripts, are often regarded as byproducts of transcription, although they may be important for the expression of nearby genes. We identified regions of nascent expression of these lincRNA in 16 human cell lines using Bru-seq techniques, and found distinctly regulated patterns of PROMPT, eRNA, and RT transcription using the diverse biochemical approaches in the ENCODE4 deeply profiled cell lines collection. Transcription of these lincRNAs was influenced by sequence-specific features and the local or 3D chromatin landscape. However, these sequence and chromatin features do not describe the full spectrum of lincRNA expression variability we identify, highlighting the complexity of their regulation. This may suggest that transcription-associated lincRNAs are not merely byproducts, but rather that the transcript itself, or the act of its transcription, is important for genomic function.

R-loop and diseases: the cell cycle matters

The cell cycle is a crucial biological process that is involved in cell growth, development, and reproduction. It can be divided into G1, S, G2, and M phases, and each period is closely regulated to ensure the production of two similar daughter cells with the same genetic material. However, many obstacles influence the cell cycle, including the R-loop that is formed throughout this process. R-loop is a triple-stranded structure, composed of an RNA: DNA hybrid and a single DNA strand, which is ubiquitous in organisms from bacteria to mammals. The existence of the R-loop has important significance for the regulation of various physiological processes. However, aberrant accumulation of R-loop due to its limited resolving ability will be detrimental for cells. For example, DNA damage and genomic instability, caused by the R-loop, can activate checkpoints in the cell cycle, which in turn induce cell cycle arrest and cell death. At present, a growing number of factors have been proven to prevent or eliminate the accumulation of R-loop thereby avoiding DNA damage and mutations. Therefore, we need to gain detailed insight into the R-loop resolution factors at different stages of the cell cycle. In this review, we review the current knowledge of factors that play a role in resolving the R-loop at different stages of the cell cycle, as well as how mutations of these factors lead to the onset and progression of diseases.

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.

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.

Telomere-specific regulation of TERRA and its impact on telomere stability

TERRA is a class of telomeric repeat-containing RNAs that are expressed from telomeres in multiple organisms. TERRA transcripts play key roles in telomere maintenance and their physiological levels are essential to maintain the integrity of telomeric DNA. Indeed, deregulated TERRA expression or its altered localization can impact telomere stability by multiple mechanisms including fueling transcription-replication conflicts, promoting resection of chromosome ends, altering the telomeric chromatin, and supporting homologous recombination. Therefore, a fine-tuned control of TERRA is important to maintain the integrity of the genome. Several studies have reported that different cell lines express substantially different levels of TERRA. Most importantly, TERRA levels markedly vary among telomeres of a given cell type, indicating the existence of telomere-specific regulatory mechanisms which may help coordinate TERRA functions. TERRA molecules contain distinct subtelomeric sequences, depending on their telomere of origin, which may instruct specific post-transcriptional modifications or mediate distinct functions. In addition, all TERRA transcripts share a repetitive G-rich sequence at their 3' end which can form DNA:RNA hybrids and fold into G-quadruplex structures. Both structures are involved in TERRA functions and can critically affect telomere stability. In this review, we examine the mechanisms controlling TERRA levels and the impact of their telomere-specific regulation on telomere stability. We compare evidence obtained in different model organisms, discussing recent advances as well as controversies in the field. Furthermore, we discuss the importance of DNA:RNA hybrids and G-quadruplex structures in the context of TERRA biology and telomere maintenance.

Regulation of transcription patterns, poly(ADP-ribose), and RNA-DNA hybrids by the ATM protein kinase

The ataxia telangiectasia mutated (ATM) protein kinase is a master regulator of the DNA damage response and also an important sensor of oxidative stress. Analysis of gene expression in ataxia-telangiectasia (A-T) patient brain tissue shows that large-scale transcriptional changes occur in patient cerebellum that correlate with the expression level and guanine-cytosine (GC) content of transcribed genes. In human neuron-like cells in culture, we map locations of poly(ADP-ribose) and RNA-DNA hybrid accumulation genome-wide with ATM inhibition and find that these marks also coincide with high transcription levels, active transcription histone marks, and high GC content. Antioxidant treatment reverses the accumulation of R-loops in transcribed regions, consistent with the central role of reactive oxygen species in promoting these lesions. Based on these results, we postulate that transcription-associated lesions accumulate in ATM-deficient cells and that the single-strand breaks and PARylation at these sites ultimately generate changes in transcription that compromise cerebellum function and lead to neurodegeneration over time in A-T patients.

Bridging the gap: R-loop mediated genomic instability and its implications in neurological diseases

R-loops, intricate three-stranded structures formed by RNA-DNA hybrids and an exposed non-template DNA strand, are fundamental to various biological phenomena. They carry out essential and contrasting functions within cellular mechanisms, underlining their critical role in maintaining cellular homeostasis. The specific cellular context that dictates R-loop formation determines their function, particularly emphasizing the necessity for their meticulous genomic regulation. Notably, the aberrant formation or misregulation of R-loops is implicated in numerous neurological disorders. This review focuses on the complex interactions between R-loops and double-strand DNA breaks, exploring how R-loop dysregulation potentially contributes to the pathogenesis of various brain disorders, which could provide novel insights into the molecular mechanisms underpinning neurological disease progression and identify potential therapeutic targets by highlighting these aspects.

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.

Human PKD1 sequences form R-loop structures in vitro

Autosomal dominant polycystic kidney disease results from the loss of the PKD1 gene product, polycystin 1. Regulatory mechanisms are unresolved, but an apparent G/C sequence bias in the gene is consistent with co-transcriptional R-loop formation. R-loops regulate gene expression and stability, and they form when newly synthesized RNA extensively pairs with the template DNA to displace the non-template strand. In this study, we tested two human PKD1 sequences for co-transcriptional R-loop formation in vitro. We observed RNase H-sensitive R-loop formation in intron 1 and 22 sequences, but only in one transcriptional orientation. Therefore, R-loops may participate in PKD1 expression or stability.

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.

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.

Regulation of transcription patterns, poly-ADP-ribose, and RNA-DNA hybrids by the ATM protein kinase

The ATM protein kinase is a master regulator of the DNA damage response and also an important sensor of oxidative stress. Analysis of gene expression in Ataxia-telangiectasia patient brain tissue shows that large-scale transcriptional changes occur in patient cerebellum that correlate with expression level and GC content of transcribed genes. In human neuron-like cells in culture we map locations of poly-ADP-ribose and RNA-DNA hybrid accumulation genome-wide with ATM inhibition and find that these marks also coincide with high transcription levels, active transcription histone marks, and high GC content. Antioxidant treatment reverses the accumulation of R-loops in transcribed regions, consistent with the central role of ROS in promoting these lesions. Based on these results we postulate that transcription-associated lesions accumulate in ATM-deficient cells and that the single-strand breaks and PARylation at these sites ultimately generate changes in transcription that compromise cerebellum function and lead to neurodegeneration over time in A-T patients.

Evaluation of potential role of R-loop and G-quadruplex DNA in the fragility of c-MYC during chromosomal translocation associated with Burkitt's lymphoma

t(8;14) translocation is the hallmark of Burkitt's lymphoma and results in c-MYC deregulation. During the translocation, c-MYC gene on chromosome 8 gets juxtaposed to the Ig switch regions on chromosome 14. Although the promoter of c-MYC has been investigated for its mechanism of fragility, little is known about other c-MYC breakpoint regions. We have analyzed the translocation break points at the exon 1/intron 1 of c-MYC locus from patients with Burkitt's lymphoma. Results showed that the breakpoint region, when present on a plasmid, could fold into an R-loop confirmation in a transcription-dependent manner. Sodium bisulfite modification assay revealed significant single-strandedness on chromosomal DNA of Burkitt's lymphoma cell line, Raji, and normal lymphocytes, revealing distinct R-loops covering up to 100 bp region. Besides, ChIP-DRIP analysis reveals that the R-loop antibody can bind to the breakpoint region. Further, we show the formation of stable parallel intramolecular G-quadruplex on non-template strand of the genome. Finally, incubation of purified AID in vitro or overexpression of AID within the cells led to enhanced mutation frequency at the c-MYC breakpoint region. Interestingly, anti-γH2AX can bind to DSBs generated at the c-MYC breakpoint region within the cells. The formation of R-loop and G-quadruplex was found to be mutually exclusive. Therefore, our results suggest that AID can bind to the single-stranded region of the R-loop and G4 DNA, leading to the deamination of cytosines to uracil and induction of DNA breaks in one of the DNA strands, leading to double-strand break, which could culminate in t(8;14) chromosomal translocation.

R-LOOPs on Short Tandem Repeat Expansion Disorders in Neurodegenerative Diseases

Expansions of short tandem repeats (STRs) have been found to be present in more than 50 diseases and have a close connection with neurodegenerative diseases. Transcriptional silencing and R-LOOP formation, RNA-mediated sequestration of RNA-binding proteins (RBPs), gain-of-function (GOF) proteins containing expanded repeats, and repeat-associated non-AUG (RAN) translation of toxic repeat peptides are some potential molecular mechanisms underlying STR expansion disorders. R-LOOP, a byproduct of transcription, is a three-stranded nucleic acid structure with abnormal accumulation that participates in the pathogenesis of STR expansion disorders by inducing DNA damage and genome instability. R-LOOPs can engender a series of DNA damage, such as DNA double-strand breaks (DSBs), single-strand breaks (SSBs), DNA recombination, or mutations in the DNA replication, transcription, or repair processes. In this review, we provide an in-depth discussion of recent advancements in R-LOOP and systematically elaborate on its genetic destabilizing effects in several neurodegenerative diseases. These molecular mechanisms will provide novel targets for drug design and therapeutic upgrading of these devastating diseases.

Cellular Responses to Widespread DNA Replication Stress

Replicative DNA polymerases are blocked by nearly all types of DNA damage. The resulting DNA replication stress threatens genome stability. DNA replication stress is also caused by depletion of nucleotide pools, DNA polymerase inhibitors, and DNA sequences or structures that are difficult to replicate. Replication stress triggers complex cellular responses that include cell cycle arrest, replication fork collapse to one-ended DNA double-strand breaks, induction of DNA repair, and programmed cell death after excessive damage. Replication stress caused by specific structures (e.g., G-rich sequences that form G-quadruplexes) is localized but occurs during the S phase of every cell division. This review focuses on cellular responses to widespread stress such as that caused by random DNA damage, DNA polymerase inhibition/nucleotide pool depletion, and R-loops. Another form of global replication stress is seen in cancer cells and is termed oncogenic stress, reflecting dysregulated replication origin firing and/or replication fork progression. Replication stress responses are often dysregulated in cancer cells, and this too contributes to ongoing genome instability that can drive cancer progression. Nucleases play critical roles in replication stress responses, including MUS81, EEPD1, Metnase, CtIP, MRE11, EXO1, DNA2-BLM, SLX1-SLX4, XPF-ERCC1-SLX4, Artemis, XPG, FEN1, and TATDN2. Several of these nucleases cleave branched DNA structures at stressed replication forks to promote repair and restart of these forks. We recently defined roles for EEPD1 in restarting stressed replication forks after oxidative DNA damage, and for TATDN2 in mitigating replication stress caused by R-loop accumulation in BRCA1-defective cells. We also discuss how insights into biological responses to genome-wide replication stress can inform novel cancer treatment strategies that exploit synthetic lethal relationships among replication stress response factors.

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.

DICER ribonuclease removes harmful R-loops

R-loops, which consist of a DNA-RNA hybrid and a displaced DNA strand, are known to threaten genome integrity. To counteract this, different mechanisms suppress R-loop accumulation by either preventing the hybridization of RNA with the DNA template (RNA biogenesis factors), unwinding the hybrid (DNA-RNA helicases), or degrading the RNA moiety of the R-loop (type H ribonucleases [RNases H]). Thus far, RNases H are the only nucleases known to cleave DNA-RNA hybrids. Now, we show that the RNase DICER also resolves R-loops. Biochemical analysis reveals that DICER acts by specifically cleaving the RNA within R-loops. Importantly, a DICER RNase mutant impaired in R-loop processing causes a strong accumulation of R-loops in cells. Our results thus not only reveal a function of DICER as an R-loop resolvase independent of DROSHA but also provide evidence for the role of multi-functional RNA processing factors in the maintenance of genome integrity in higher eukaryotes.

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.

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.

TERRA stability is regulated by RALY and polyadenylation in a telomere-specific manner

Telomeric repeat-containing RNA (TERRA) is a long non-coding RNA transcribed from telomeres that plays key roles in telomere maintenance. A fraction of TERRA is polyadenylated, and the presence of the poly(A) tail influences TERRA localization and stability. However, the mechanisms of TERRA biogenesis remain mostly elusive. Here, we show that the stability of TERRA transcripts is regulated by the RNA-binding protein associated with lethal yellow mutation (RALY). RALY depletion results in lower TERRA levels, impaired localization of TERRA at telomeres, and ultimately telomere damage. Importantly, we show that TERRA polyadenylation is telomere specific and that RALY preferentially stabilizes non-polyadenylated TERRA transcripts. Finally, we report that TERRA interacts with the poly(A)-binding protein nuclear 1 (PABPN1). Altogether, our results indicate that TERRA stability is regulated by the interplay between RALY and PABPN1, defined by the TERRA polyadenylation state. Our findings also suggest that different telomeres may trigger distinct TERRA-mediated responses.

R-Loops and R-Loop-Binding Proteins in Cancer Progression and Drug Resistance

R-loops are three-stranded DNA/RNA hybrids that form by the annealing of the mRNA transcript to its coding template while displacing the non-coding strand. While R-loop formation regulates physiological genomic and mitochondrial transcription and DNA damage response, imbalanced R-loop formation can be a threat to the genomic integrity of the cell. As such, R-loop formation is a double-edged sword in cancer progression, and perturbed R-loop homeostasis is observed across various malignancies. Here, we discuss the interplay between R-loops and tumor suppressors and oncogenes, with a focus on BRCA1/2 and ATR. R-loop imbalances contribute to cancer propagation and the development of chemotherapy drug resistance. We explore how R-loop formation can cause cancer cell death in response to chemotherapeutics and be used to circumvent drug resistance. As R-loop formation is tightly linked to mRNA transcription, their formation is unavoidable in cancer cells and can thus be explored in novel cancer therapeutics.

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.

PARP1 associates with R-loops to promote their resolution and genome stability

PARP1 is a DNA-dependent ADP-Ribose transferase with ADP-ribosylation activity that is triggered by DNA breaks and non-B DNA structures to mediate their resolution. PARP1 was also recently identified as a component of the R-loop-associated protein-protein interaction network, suggesting a potential role for PARP1 in resolving this structure. R-loops are three-stranded nucleic acid structures that consist of a RNA-DNA hybrid and a displaced non-template DNA strand. R-loops are involved in crucial physiological processes but can also be a source of genome instability if persistently unresolved. In this study, we demonstrate that PARP1 binds R-loops in vitro and associates with R-loop formation sites in cells which activates its ADP-ribosylation activity. Conversely, PARP1 inhibition or genetic depletion causes an accumulation of unresolved R-loops which promotes genomic instability. Our study reveals that PARP1 is a novel sensor for R-loops and highlights that PARP1 is a suppressor of R-loop-associated genomic instability.

Emerging regulatory mechanisms of noncoding RNAs in topologically associating domains

Topologically associating domains (TADs) are integral to spatial genome organization, instructing gene expression, and cell fate. Recently, several advances have uncovered roles for noncoding RNAs (ncRNAs) in the regulation of the form and function of mammalian TADs. Phase separation has also emerged as a potential arbiter of ncRNAs in the regulation of TADs. In this review we discuss the implications of these novel findings in relation to how ncRNAs might structurally and functionally regulate TADs from two perspectives: moderating loop extrusion through interactions with architectural proteins, and facilitating TAD phase separation. Additionally, we propose future studies and directions to investigate these phenomena.

Exploration and analysis of R-loop mapping data with RLBase

R-loops are three-stranded nucleic acid structures formed from the hybridization of RNA and DNA. In 2012, Ginno et al. introduced the first R-loop mapping method. Since that time, dozens of R-loop mapping studies have been conducted, yielding hundreds of publicly available datasets. Current R-loop databases provide only limited access to these data. Moreover, no web tools for analyzing user-supplied R-loop datasets have yet been described. In our recent work, we reprocessed 810 R-loop mapping samples, building the largest R-loop data resource to date. We also defined R-loop consensus regions and developed a framework for R-loop data analysis. Now, we introduce RLBase, a user-friendly database that provides the capability to (i) explore hundreds of public R-loop mapping datasets, (ii) explore R-loop consensus regions, (iii) analyze user-supplied data and (iv) download standardized and reprocessed datasets. RLBase is directly accessible via the following URL: https://gccri.bishop-lab.uthscsa.edu/shiny/rlbase/.

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.

Chromatin regulators in DNA replication and genome stability maintenance during S-phase

The duplication of genetic information is central to life. The replication of genetic information is strictly controlled to ensure that each piece of genomic DNA is copied only once during a cell cycle. Factors that slow or stop replication forks cause replication stress. Replication stress is a major source of genome instability in cancer cells. Multiple control mechanisms facilitate the unimpeded fork progression, prevent fork collapse and coordinate fork repair. Chromatin alterations, caused by histone post-translational modifications and chromatin remodeling, have critical roles in normal replication and in avoiding replication stress and its consequences. This text reviews the chromatin regulators that ensure DNA replication and the proper response to replication stress. We also briefly touch on exploiting replication stress in therapeutic strategies. As chromatin regulators are frequently mutated in cancer, manipulating their activity could provide many possibilities for personalized treatment.

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 landscape in mature human sperm: Regulatory and evolutionary implications

R-loops are three-stranded nucleic acid structures consisting of an RNA:DNA hybrid and a displaced DNA strand. While R-loops pose a potential threat to genome integrity, they constitute 5% of the human genome. The role of R-loops in transcriptional regulation, DNA replication, and chromatin signature is becoming increasingly clear. R-loops are associated with various histone modifications, suggesting that they may modulate chromatin accessibility. To potentially harness transcription-coupled repair mechanisms in the germline, nearly the entire genome is expressed during the early stages of male gametogenesis in mammals, providing ample opportunity for the formation of a transcriptome-dependent R-loop landscape in male germ cells. In this study, our data demonstrated the presence of R-loops in fully mature human and bonobo sperm heads and their partial correspondence to transcribed regions and chromatin structure, which is massively reorganized from mainly histone to mainly protamine-packed chromatin in mature sperm. The sperm R-loop landscape resembles characteristic patterns of somatic cells. Surprisingly, we detected R-loops in both residual histone and protamine-packed chromatin and localize them to still-active retroposons, ALUs and SINE-VNTR-ALUs (SVAs), the latter has recently arisen in hominoid primates. We detected both evolutionarily conserved and species-specific localizations. Comparing our DNA-RNA immunoprecipitation (DRIP) data with published DNA methylation and histone chromatin immunoprecipitation (ChIP) data, we hypothesize that R-loops epigenetically reduce methylation of SVAs. Strikingly, we observe a strong influence of R-loops on the transcriptomes of zygotes from early developmental stages before zygotic genome activation. Overall, these findings suggest that chromatin accessibility influenced by R-loops may represent a system of inherited gene regulation.

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.

Branchpoint translocation by fork remodelers as a general mechanism of R-loop removal

Co-transcriptional R loops arise from stalling of RNA polymerase, leading to the formation of stable DNA:RNA hybrids. Unresolved R loops promote genome instability but are counteracted by helicases and nucleases. Here, we show that branchpoint translocases are a third class of R-loop-displacing enzyme in vitro. In cells, deficiency in the Fanconi-anemia-associated branchpoint translocase FANCM causes R-loop accumulation, particularly after treatment with DNA:RNA-hybrid-stabilizing agents. This correlates with FANCM localization at R-loop-prone regions of the genome. Moreover, other branchpoint translocases associated with human disease, such as SMARCAL1 and ZRANB3, and those from lower organisms can also remove R loops in vitro. Branchpoint translocases are more potent than helicases in resolving R loops, indicating their evolutionary important role in R-loop suppression. In human cells, FANCM, SMARCAL1, and ZRANB3 depletion causes additive effects on R-loop accumulation and DNA damage. Our work reveals a mechanistic basis for R-loop displacement that is linked to genome stability.

Comprehensive computational analysis of epigenetic descriptors affecting CRISPR-Cas9 off-target activity

BACKGROUND

A common issue in CRISPR-Cas9 genome editing is off-target activity, which prevents the widespread use of CRISPR-Cas9 in medical applications. Among other factors, primary chromatin structure and epigenetics may influence off-target activity.

METHODS

In this work, we utilize crisprSQL, an off-target database, to analyze the effect of 19 epigenetic descriptors on CRISPR-Cas9 off-target activity. Termed as 19 epigenetic features/scores, they consist of 6 experimental epigenetic and 13 computed nucleosome organization-related features. In terms of novel features, 15 of the epigenetic scores are newly considered. The 15 newly considered scores consist of 13 freshly computed nucleosome occupancy/positioning scores and 2 experimental features (MNase and DRIP). The other 4 existing scores are experimental features (CTCF, DNase I, H3K4me3, RRBS) commonly used in deep learning models for off-target activity prediction. For data curation, MNase was aggregated from existing experimental nucleosome occupancy data. Based on the sequence context information available in crisprSQL, we also computed nucleosome occupancy/positioning scores for off-target sites.

RESULTS

To investigate the relationship between the 19 epigenetic features and off-target activity, we first conducted Spearman and Pearson correlation analysis. Such analysis shows that some computed scores derived from training-based models and training-free algorithms outperform all experimental epigenetic features. Next, we evaluated the contribution of all epigenetic features in two successful machine/deep learning models which predict off-target activity. We found that some computed scores, unlike all 6 experimental features, significantly contribute to the predictions of both models. As a practical research contribution, we make the off-target dataset containing all 19 epigenetic features available to the research community.

CONCLUSIONS

Our comprehensive computational analysis helps the CRISPR-Cas9 community better understand the relationship between epigenetic features and CRISPR-Cas9 off-target activity.

DDX17 helicase promotes resolution of R-loop-mediated transcription-replication conflicts in human cells

R-loops are three-stranded nucleic acid structures composed of an RNA:DNA hybrid and displaced DNA strand. These structures can halt DNA replication when formed co-transcriptionally in the opposite orientation to replication fork progression. A recent study has shown that replication forks stalled by co-transcriptional R-loops can be restarted by a mechanism involving fork cleavage by MUS81 endonuclease, followed by ELL-dependent reactivation of transcription, and fork religation by the DNA ligase IV (LIG4)/XRCC4 complex. However, how R-loops are eliminated to allow the sequential restart of transcription and replication in this pathway remains elusive. Here, we identified the human DDX17 helicase as a factor that associates with R-loops and counteracts R-loop-mediated replication stress to preserve genome stability. We show that DDX17 unwinds R-loops in vitro and promotes MUS81-dependent restart of R-loop-stalled forks in human cells in a manner dependent on its helicase activity. Loss of DDX17 helicase induces accumulation of R-loops and the formation of R-loop-dependent anaphase bridges and micronuclei. These findings establish DDX17 as a component of the MUS81-LIG4-ELL pathway for resolution of R-loop-mediated transcription-replication conflicts, which may be involved in R-loop unwinding.

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.

Shared properties of gene transfer agent and core genes revealed by comparative genomics of Alphaproteobacteria

Gene transfer agents (GTAs) are phage-like particles that transfer pieces of cellular genomic DNA to other cells. Homologues of the Rhodobacter capsulatus GTA (RcGTA) structural genes are widely distributed in the Alphaproteobacteria and particularly well conserved in the order Rhodobacterales. Possible reasons for their widespread conservation are still being discussed. It has been suggested that these alphaproteobacterial elements originate from a prophage that was present in an ancestral bacterium and subsequently evolved into a GTA that is now widely maintained in extant descendant lineages. Here, we analysed genomic properties that might relate to the conservation of these alphaproteobacterial GTAs. This revealed that the chromosomal locations of the GTA gene clusters are biased. They primarily occur on the leading strand of DNA replication, at large distances from long repetitive elements, and thus are in regions of lower plasticity, and in areas of extreme GC skew, which also accumulate core genes. These extreme GC skew regions arise from the preferential use of codons with an excess of G over C, a distinct phenomenon from the elevated GC content that has previously been found to be associated with GTA genes. The observed properties, along with their high level of conservation, show that GTA genes share multiple features with core genes in the examined lineages of the Alphaproteobacteria.

The pausing zone and control of RNA polymerase II elongation by Spt5: Implications for the pause-release model

The pause-release model of transcription proposes that 40-100 bases from the start site RNA Pol II pauses, followed by release into productive elongation. Pause release is facilitated by the PTEFb phosphorylation of the RNA Pol II elongation factor, Spt5. We mapped paused polymerases by eNET-seq and found frequent pausing in zones that extend ∼0.3-3 kb into genes even when PTEFb is inhibited. The fraction of paused polymerases or pausing propensity declines gradually over several kb and not abruptly as predicted for a discrete pause-release event. Spt5 depletion extends pausing zones, suggesting that it promotes the maturation of elongation complexes to a low-pausing state. The expression of mutants after Spt5 depletion showed that phosphomimetic substitutions in the CTR1 domain diminished pausing throughout genes. By contrast, mutants that prevent the phosphorylation of the Spt5 RNA-binding domain strengthened pausing. Thus, distinct Spt5 phospho-isoforms set the balance between pausing and elongation.

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.

Enhancer RNAs step forward: new insights into enhancer function

Enhancers confer precise spatiotemporal patterns of gene expression in response to developmental and environmental stimuli. Over the last decade, the transcription of enhancer RNAs (eRNAs) – nascent RNAs transcribed from active enhancers – has emerged as a key factor regulating enhancer activity. eRNAs are relatively short-lived RNA species that are transcribed at very high rates but also quickly degraded. Nevertheless, eRNAs are deeply intertwined within enhancer regulatory networks and are implicated in a number of transcriptional control mechanisms. Enhancers show changes in function and sequence over evolutionary time, raising questions about the relationship between enhancer sequences and eRNA function. Moreover, the vast majority of single nucleotide polymorphisms associated with human complex diseases map to the non-coding genome, with causal disease variants enriched within enhancers. In this Primer, we survey the diverse roles played by eRNAs in enhancer-dependent gene expression, evaluating different models for eRNA function. We also explore questions surrounding the genetic conservation of enhancers and how this relates to eRNA function and dysfunction.

Summary: This Primer evaluates the ideas that underpin developing models for eRNA function, exploring cases in which perturbed eRNA function contributes to disease.

Relationships between genome-wide R-loop distribution and classes of recurrent DNA breaks in neural stem/progenitor cells

Recent studies revealed classes of recurrent DNA double-strand breaks (DSBs) in neural stem/progenitor cells, including transcription-associated, promoter-proximal breaks and recurrent DSB clusters in late-replicating, long neural genes that may give rise to somatic brain mosaicism. The mechanistic factors promoting these different classes of DSBs in neural stem/progenitor cells are not understood. Here, we elucidated the genome-wide landscape of RNA:DNA hybrid structures called “R-loops” in primary neural stem/progenitor cells undergoing aphidicolin-induced, mild replication stress to assess the potential contribution of R-loops to the different, recurrent classes of DNA break “hotspots”. We find that R-loops in neural stem/progenitor cells undergoing mild replication stress are present primarily in early-replicating, transcribed regions and in genes with promoter GC skew that are associated with cell lineage-specific processes. Surprisingly, most long, neural genes that form recurrent DSB clusters do not show R-loop formation under conditions of mild replication stress. Our findings are consistent with a role of R-loop-associated processes in promoter-proximal DNA break formation in highly transcribed, early replicating regions but suggest that R-loops do not drive replication stress-induced, recurrent DSB cluster formation in most long, neural genes.

Integrative analysis and prediction of human R-loop binding proteins

In the past decade, there has been a growing appreciation for R-loop structures as important regulators of the epigenome, telomere maintenance, DNA repair, and replication. Given these numerous functions, dozens, or potentially hundreds, of proteins could serve as direct or indirect regulators of R-loop writing, reading, and erasing. In order to understand common properties shared amongst potential R-loop binding proteins, we mined published proteomic studies and distilled 10 features that were enriched in R-loop binding proteins compared with the rest of the proteome. Applying an easy-ensemble machine learning approach, we used these R-loop binding protein-specific features along with their amino acid composition to create random forest classifiers that predict the likelihood of a protein to bind to R-loops. Known R-loop regulating pathways such as splicing, DNA damage repair and chromatin remodeling are highly enriched in our datasets, and we validate 2 new R-loop binding proteins LIG1 and FXR1 in human cells. Together these datasets provide a reference to pursue analyses of novel R-loop regulatory proteins.

Quality-controlled R-loop meta-analysis reveals the characteristics of R-loop consensus regions

R-loops are three-stranded nucleic acid structures formed from the hybridization of RNA and DNA. While the pathological consequences of R-loops have been well-studied to date, the locations, classes, and dynamics of physiological R-loops remain poorly understood. R-loop mapping studies provide insight into R-loop dynamics, but their findings are challenging to generalize. This is due to the narrow biological scope of individual studies, the limitations of each mapping modality, and, in some cases, poor data quality. In this study, we reprocessed 810 R-loop mapping datasets from a wide array of biological conditions and mapping modalities. From this data resource, we developed an accurate R-loop data quality control method, and we reveal the extent of poor-quality data within previously published studies. We then identified a set of high-confidence R-loop mapping samples and used them to define consensus R-loop sites called 'R-loop regions' (RL regions). In the process, we identified a stark divergence between RL regions detected by S9.6 and dRNH-based mapping methods, particularly with respect to R-loop size, location, and colocalization with RNA binding factors. Taken together, this work provides a much-needed method to assess R-loop data quality and offers novel context regarding the differences between dRNH- and S9.6-based R-loop mapping approaches.

Multi-omics to characterize the functional relationships of R-loops with epigenetic modifications, RNAPII transcription and gene expression

Abnormal accumulation of R-loops results in replication stress, genome instability, chromatin alterations and gene silencing. Little research has been done to characterize functional relationships among R-loops, histone marks, RNA polymerase II (RNAPII) transcription and gene regulation. We built extremely randomized trees (ETs) models to predict the genome-wide R-loops using RNAPII and multiple histone modifications chromatin immunoprecipitation (ChIP)-seq, DNase-seq, Global Run-On sequencing (GRO-seq) and R-loop profiling data. We compared the performance of ET models to multiple machine learning approaches, and the proposed ET models achieved the best and extremely robust performances. Epigenetic profiles are highly predictive of R-loops genome-widely and they are strongly associated with R-loop formation. In addition, the presence of R-loops is significantly correlated with RNAPII transcription activity, H3K4me3 and open chromatin around the transcription start site, and H3K9me1 and H3K9me3 around the transcription termination site. RNAPII pausing defects were correlated with 5′R-loops accumulation, and transcriptional termination defects and read-throughs were correlated with 3′R-loops accumulation. Furthermore, we found driver genes with 5′R-loops and RNAPII pausing defects express significantly higher and genes with 3′R-loops and read-through transcription express significantly lower than genes without R-loops. These driver genes are enriched with chromosomal instability, Hippo–Merlin signaling Dysregulation, DNA damage response and TGF-β pathways, indicating R-loops accumulating at the 5′ end of genes play oncogenic roles, whereas at the 3′ end of genes play tumor-suppressive roles in tumorigenesis.

ZFP281-BRCA2 prevents R-loop accumulation during DNA replication

R-loops are prevalent in mammalian genomes and involved in many fundamental cellular processes. Depletion of BRCA2 leads to aberrant R-loop accumulation, contributing to genome instability. Here, we show that ZFP281 cooperates with BRCA2 in preventing R-loop accumulation to facilitate DNA replication in embryonic stem cells. ZFP281 depletion reduces PCNA levels on chromatin and impairs DNA replication. Mechanistically, we demonstrate that ZFP281 can interact with BRCA2, and that BRCA2 is enriched at G/C-rich promoters and requires both ZFP281 and PRC2 for its proper recruitment to the bivalent chromatin at the genome-wide scale. Furthermore, depletion of ZFP281 or BRCA2 leads to accumulation of R-loops over the bivalent regions, and compromises activation of the developmental genes by retinoic acid during stem cell differentiation. In summary, our results reveal that ZFP281 recruits BRCA2 to the bivalent chromatin regions to ensure proper progression of DNA replication through preventing persistent R-loops.

R-loops are prevalent in mammalian genomes and involved in many fundamental cellular processes. Here, Wang et al. report that ZFP281 cooperates with BRCA2 in preventing R-loop accumulation to facilitate DNA replication in embryonic stem cells.