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.

RNA helicase DDX3 regulates RAD51 localization and DNA damage repair in Ewing sarcoma.

We previously demonstrated that RNA helicase DDX3X (DDX3) can be a therapeutic target in Ewing sarcoma (EWS), but its role in EWS biology remains unclear. The present work demonstrates that DDX3 plays a unique role in DNA damage repair (DDR). We show that DDX3 interacts with several proteins involved in homologous recombination, including RAD51, RECQL1, RPA32, and XRCC2. In particular, DDX3 colocalizes with RAD51 and RNA:DNA hybrid structures in the cytoplasm of EWS cells. Inhibition of DDX3 RNA helicase activity increases cytoplasmic RNA:DNA hybrids, sequestering RAD51 in the cytoplasm, which impairs nuclear translocation of RAD51 to sites of double-stranded DNA breaks, thus increasing sensitivity of EWS to radiation treatment, both in vitro and in vivo. This discovery lays the foundation for exploring new therapeutic approaches directed at manipulating DDR protein localization in solid tumors.

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.

RNA Helicase DDX3 Regulates RAD51 Localization and DNA Damage Repair in Ewing Sarcoma.

We previously demonstrated that RNA helicase DDX3X (DDX3) can be a therapeutic target in Ewing sarcoma (EWS), but its role in EWS biology remains unclear. The present work demonstrates that DDX3 plays a unique role in DNA damage repair (DDR). We show that DDX3 interacts with several proteins involved in homologous recombination, including RAD51, RECQL1, RPA32, and XRCC2. In particular, DDX3 colocalizes with RAD51 and RNA:DNA hybrid structures in the cytoplasm of EWS cells. Inhibition of DDX3 RNA helicase activity increases cytoplasmic RNA:DNA hybrids, sequestering RAD51 in the cytoplasm, which impairs nuclear translocation of RAD51 to sites of double-stranded DNA breaks thus increasing sensitivity of EWS to radiation treatment, both in vitro and in vivo. This discovery lays the foundation for exploring new therapeutic approaches directed at manipulating DDR protein localization in solid tumors.

Epigenetic Control of Translation Checkpoint and Tumor Progression via RUVBL1-EEF1A1 Axis.

Epigenetic dysregulation is reported in multiple cancers including Ewing sarcoma (EwS). However, the epigenetic networks underlying the maintenance of oncogenic signaling and therapeutic response remain unclear. Using a series of epigenetics- and complex-focused CRISPR screens, RUVBL1, the ATPase component of NuA4 histone acetyltransferase complex, is identified to be essential for EwS tumor progression. Suppression of RUVBL1 leads to attenuated tumor growth, loss of histone H4 acetylation, and ablated MYC signaling. Mechanistically, RUVBL1 controls MYC chromatin binding and modulates the MYC-driven EEF1A1 expression and thus protein synthesis. High-density CRISPR gene body scan pinpoints the critical MYC interacting residue in RUVBL1. Finally, this study reveals the synergism between RUVBL1 suppression and pharmacological inhibition of MYC in EwS xenografts and patient-derived samples. These results indicate that the dynamic interplay between chromatin remodelers, oncogenic transcription factors, and protein translation machinery can provide novel opportunities for combination cancer therapy.

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/.

An international working group consensus report for the prioritization of molecular biomarkers for Ewing sarcoma.

The advent of dose intensified interval compressed therapy has improved event-free survival for patients with localized Ewing sarcoma (EwS) to 78% at 5 years. However, nearly a quarter of patients with localized tumors and 60-80% of patients with metastatic tumors suffer relapse and die of disease. In addition, those who survive are often left with debilitating late effects. Clinical features aside from stage have proven inadequate to meaningfully classify patients for risk-stratified therapy. Therefore, there is a critical need to develop approaches to risk stratify patients with EwS based on molecular features. Over the past decade, new technology has enabled the study of multiple molecular biomarkers in EwS. Preliminary evidence requiring validation supports copy number changes, and loss of function mutations in tumor suppressor genes as biomarkers of outcome in EwS. Initial studies of circulating tumor DNA demonstrated that diagnostic ctDNA burden and ctDNA clearance during induction are also associated with outcome. In addition, fusion partner should be a pre-requisite for enrollment on EwS clinical trials, and the fusion type and structure require further study to determine prognostic impact. These emerging biomarkers represent a new horizon in our understanding of disease risk and will enable future efforts to develop risk-adapted treatment.

Systematic Analysis of Long Non-Coding RNA Genes in Nonalcoholic Fatty Liver Disease.

The largest solid organ in humans, the liver, performs a variety of functions to sustain life. When damaged, cells in the liver can regenerate themselves to maintain normal liver physiology. However, some damage is beyond repair, which necessitates liver transplantation. Increasing rates of obesity, Western diets (i.e., rich in processed carbohydrates and saturated fats), and cardiometabolic diseases are interlinked to liver diseases, including non-alcoholic fatty liver disease (NAFLD), which is a collective term to describe the excess accumulation of fat in the liver of people who drink little to no alcohol. Alarmingly, the prevalence of NAFLD extends to 25% of the world population, which calls for the urgent need to understand the disease mechanism of NAFLD. Here, we performed secondary analyses of published RNA sequencing (RNA-seq) data of NAFLD patients compared to healthy and obese individuals to identify long non-coding RNAs (lncRNAs) that may underly the disease mechanism of NAFLD. Similar to protein-coding genes, many lncRNAs are dysregulated in NAFLD patients compared to healthy and obese individuals, suggesting that understanding the functions of dysregulated lncRNAs may shed light on the pathology of NAFLD. To demonstrate the functional importance of lncRNAs in the liver, loss-of-function experiments were performed for one NAFLD-related lncRNA, LINC01639, which showed that it is involved in the regulation of genes related to apoptosis, TNF/TGF, cytokine signaling, and growth factors as well as genes upregulated in NAFLD. Since there is no lncRNA database focused on the liver, especially NAFLD, we built a web database, LiverDB, to further facilitate functional and mechanistic studies of hepatic lncRNAs.

Development and characterization of a mass cytometry panel for detecting the effect of acute doxorubicin exposure on murine cardiac nonmyocytes.

Childhood cancer survivors (CCSs) face lifelong side effects related to their treatment with chemotherapy. Anthracycline agents, such as doxorubicin (DOX), are important in the treatment of childhood cancers but are associated with cardiotoxicity. Cardiac toxicities represent a significant source of chronic disability that cancer survivors face; despite this, the chronic cardiotoxicity phenotype and how it relates to acute toxicity remains poorly defined. To address this critical knowledge gap, we studied the acute effect of DOX on murine cardiac nonmyocytes in vivo. Determination of the acute cellular effects of DOX on nonmyocytes, a cell pool with finite replicative capacity, provides a basis for understanding the pathogenesis of the chronic heart disease that CCSs face. To investigate the acute cellular effects of DOX, we present single-cell RNA sequencing (scRNAseq) data from homeostatic cardiac nonmyocytes and compare it with preexisting datasets, as well as a novel CyTOF datasets. SCANPY, a python-based single-cell analysis, was used to assess the heterogeneity of cells detected in scRNAseq and CyTOF. To further assist in CyTOF data annotation, joint analyses of scRNAseq and CyTOF data using an artificial neural network known as sparse autoencoder for clustering, imputation, and embedding (SAUCIE) are performed. Lastly, the panel is tested on a mouse model of acute DOX exposure at two time points (24 and 72 h) after the last dose of doxorubicin and examined with joint clustering. In sum, we report the first ever CyTOF study of cardiac nonmyocytes and characterize the effect of acute DOX exposure with scRNAseq and CyTOF.NEW & NOTEWORTHY We describe the first mass cytometry studies of murine cardiac nonmyocytes. The mass cytometry panel is compared with single-cell RNA sequencing data. Homeostatic cardiac nonmyocytes are characterized by mass cytometry to identify and quantify four major cell populations: endothelial cells, fibroblasts, leukocytes, and pericytes. The single-cell acute nonmyocyte response to doxorubicin is studied at 24 and 72 h after doxorubicin exposure given daily for 5 days at a dose of 4 mg/kg/day.

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.

Current Status of Epitranscriptomic Marks Affecting lncRNA Structures and Functions.

Long non-coding RNAs (lncRNAs) belong to a class of non-protein-coding RNAs with their lengths longer than 200 nucleotides. Most of the mammalian genome is transcribed as RNA, yet only a small percent of the transcribed RNA corresponds to exons of protein-coding genes. Thus, the number of lncRNAs is predicted to be several times higher than that of protein-coding genes. Because of sheer number of lncRNAs, it is often difficult to elucidate the functions of all lncRNAs, especially those arising from their relationship to their binding partners, such as DNA, RNA, and proteins. Due to their binding to other macromolecules, it has become evident that the structures of lncRNAs influence their functions. In this regard, the recent development of epitranscriptomics (the field of study to investigate RNA modifications) has become important to further elucidate the structures and functions of lncRNAs. In this review, the current status of lncRNA structures and functions influenced by epitranscriptomic marks is discussed.

FibroDB: Expression Analysis of Protein-Coding and Long Non-Coding RNA Genes in Fibrosis.

Most long non-coding RNAs (lncRNAs) are expressed at lower levels than protein-coding genes and their expression is often restricted to specific cell types, certain time points during development, and various stress and disease conditions, respectively. To revisit this long-held concept, we focused on fibroblasts, a common cell type in various organs and tissues. Using fibroblasts and changes in their expression profiles during fibrosis as a model system, we show that the overall expression level of lncRNA genes is significantly lower than that of protein-coding genes. Furthermore, we identified lncRNA genes whose expression is upregulated during fibrosis. Using dermal fibroblasts as a model, we performed loss-of-function experiments and show that the knockdown of the lncRNAs LINC00622 and LINC01711 result in gene expression changes associated with cellular and inflammatory responses, respectively. Since there are no lncRNA databases focused on fibroblasts and fibrosis, we built a web application, FibroDB, to further promote functional and mechanistic studies of fibrotic lncRNAs.

BRCA2 Promotes Spontaneous Homologous Recombination In Vivo.

BACKGROUND: BRCA2 is known to be a tumor suppressor involved in homologous recombination repair and presumed to prevent genome instability in normal tissues prior to the development of tumors. Typical assessment of BRCA2 deficiency on the genome involves cell-based models using cancer cells with mixed genetic contexts, but the role in normal tissue in vivo has not been clearly demonstrated. METHODS: Using conditional deletion of Brca2 exon 11, the region containing all eight BRC repeats, in the retinal pigment epithelium and the pink-eyed unstable mouse model, we evaluate the frequency of DNA deletion events. RESULTS: In the current study, we show that conditional loss of Brca2 exon 11 results in a decreased frequency of spontaneous homologous recombination compared to wild-type mice. Of note, we observe no apparent concomitant increase in events that indicate single-strand annealing by the pink-eyed unstable mouse model. CONCLUSIONS: Therefore, our results demonstrate that BRCA2, as expected, is required for high-fidelity homologous recombination DNA repair in normal tissues, here in a tissue undergoing normal proliferation through normal development.

Correlation AnalyzeR: functional predictions from gene co-expression correlations.

BACKGROUND: Co-expression correlations provide the ability to predict gene functionality within specific biological contexts, such as different tissue and disease conditions. However, current gene co-expression databases generally do not consider biological context. In addition, these tools often implement a limited range of unsophisticated analysis approaches, diminishing their utility for exploring gene functionality and gene relationships. Furthermore, they typically do not provide the summary visualizations necessary to communicate these results, posing a significant barrier to their utilization by biologists without computational skills. RESULTS: We present Correlation AnalyzeR, a user-friendly web interface for exploring co-expression correlations and predicting gene functions, gene-gene relationships, and gene set topology. Correlation AnalyzeR provides flexible access to its database of tissue and disease-specific (cancer vs normal) genome-wide co-expression correlations, and it also implements a suite of sophisticated computational tools for generating functional predictions with user-friendly visualizations. In the usage example provided here, we explore the role of BRCA1-NRF2 interplay in the context of bone cancer, demonstrating how Correlation AnalyzeR can be effectively implemented to generate and support novel hypotheses. CONCLUSIONS: Correlation AnalyzeR facilitates the exploration of poorly characterized genes and gene relationships to reveal novel biological insights. The database and all analysis methods can be accessed as a web application at https://gccri.bishop-lab.uthscsa.edu/correlation-analyzer/ and as a standalone R package at https://github.com/Bishop-Laboratory/correlationAnalyzeR .

Bloom syndrome in a Mexican American family with rhabdomyosarcoma: evidence of a Mexican founder mutation.

Bloom syndrome is a rare autosomal recessive disorder with less than 300 cases reported in the literature. Bloom syndrome is characterized by chromosome instability, physical stigmata, growth deficiency, immunodeficiency, and a predisposition to cancer, most commonly leukemias, although solid tumors are reported as well. Bloom syndrome occurs in multiple ethnic groups with a higher incidence in persons of Ashkenazi Jewish origin. Few patients of Hispanic ethnicity have been reported. We report here a Mexican American family with a BLM pathogenic variant, c.2506_2507delAG, previously reported in a single patient from Mexico. In this family of four siblings, three have phenotypic features of Bloom syndrome, and BLM gene mutation was homozygous in these affected individuals. Our proband developed a rhabdomyosarcoma. Analysis of surrounding markers in the germline DNA revealed a common haplotype, suggesting a previously unrecognized founder mutation in the Hispanic population of Mexican origin.

Nucleolar RNA polymerase II drives ribosome biogenesis.

Proteins are manufactured by ribosomes-macromolecular complexes of protein and RNA molecules that are assembled within major nuclear compartments called nucleoli1,2. Existing models suggest that RNA polymerases I and III (Pol I and Pol III) are the only enzymes that directly mediate the expression of the ribosomal RNA (rRNA) components of ribosomes. Here we show, however, that RNA polymerase II (Pol II) inside human nucleoli operates near genes encoding rRNAs to drive their expression. Pol II, assisted by the neurodegeneration-associated enzyme senataxin, generates a shield comprising triplex nucleic acid structures known as R-loops at intergenic spacers flanking nucleolar rRNA genes. The shield prevents Pol I from producing sense intergenic noncoding RNAs (sincRNAs) that can disrupt nucleolar organization and rRNA expression. These disruptive sincRNAs can be unleashed by Pol II inhibition, senataxin loss, Ewing sarcoma or locus-associated R-loop repression through an experimental system involving the proteins RNaseH1, eGFP and dCas9 (which we refer to as 'red laser'). We reveal a nucleolar Pol-II-dependent mechanism that drives ribosome biogenesis, identify disease-associated disruption of nucleoli by noncoding RNAs, and establish locus-targeted R-loop modulation. Our findings revise theories of labour division between the major RNA polymerases, and identify nucleolar Pol II as a major factor in protein synthesis and nuclear organization, with potential implications for health and disease.

Cohesin SA1 and SA2 are RNA binding proteins that localize to RNA containing regions on DNA.

Cohesin SA1 (STAG1) and SA2 (STAG2) are key components of the cohesin complex. Previous studies have highlighted the unique contributions by SA1 and SA2 to 3D chromatin organization, DNA replication fork progression, and DNA double-strand break (DSB) repair. Recently, we discovered that cohesin SA1 and SA2 are DNA binding proteins. Given the recently discovered link between SA2 and RNA-mediated biological pathways, we investigated whether or not SA1 and SA2 directly bind to RNA using a combination of bulk biochemical assays and single-molecule techniques, including atomic force microscopy (AFM) and the DNA tightrope assay. We discovered that both SA1 and SA2 bind to various RNA containing substrates, including ssRNA, dsRNA, RNA:DNA hybrids, and R-loops. Importantly, both SA1 and SA2 localize to regions on dsDNA that contain RNA. We directly compared the SA1/SA2 binding and R-loops sites extracted from Chromatin Immunoprecipitation sequencing (ChIP-seq) and DNA-RNA Immunoprecipitation sequencing (DRIP-Seq) data sets, respectively. This analysis revealed that SA1 and SA2 binding sites overlap significantly with R-loops. The majority of R-loop-localized SA1 and SA2 are also sites where other subunits of the cohesin complex bind. These results provide a new direction for future investigation of the diverse biological functions of SA1 and SA2.

Reconstruction of Ewing Sarcoma Developmental Context from Mass-Scale Transcriptomics Reveals Characteristics of EWSR1-FLI1 Permissibility.

Ewing sarcoma is an aggressive pediatric cancer of enigmatic cellular origins typically resulting from a single translocation event t (11; 22) (q24; q12). The resulting fusion gene, EWSR1-FLI1, is toxic or unstable in most primary tissues. Consequently, attempts to model Ewing sarcomagenesis have proven unsuccessful thus far, highlighting the need to identify the cellular features which permit stable EWSR1-FLI1 expression. By re-analyzing publicly available RNA-Sequencing data with manifold learning techniques, we uncovered a group of Ewing-like tissues belonging to a developmental trajectory between pluripotent, neuroectodermal, and mesodermal cell states. Furthermore, we demonstrated that EWSR1-FLI1 expression levels control the activation of these developmental trajectories within Ewing sarcoma cells. Subsequent analysis and experimental validation demonstrated that the capability to resolve R-loops and mitigate replication stress are probable prerequisites for stable EWSR1-FLI1 expression in primary tissues. Taken together, our results demonstrate how EWSR1-FLI1 hijacks developmental gene programs and advances our understanding of Ewing sarcomagenesis.

Phase 1/2 trial of talazoparib in combination with temozolomide in children and adolescents with refractory/recurrent solid tumors including Ewing sarcoma: A Children's Oncology Group Phase 1 Consortium study (ADVL1411).

PURPOSE: We conducted a phase 1/2 trial of the poly(ADP-ribose) polymerase 1/2 inhibitor talazoparib in combination with low-dose temozolomide (TMZ) to determine the dose-limiting toxicities (DLTs), recommended phase 2 dose (RP2D), and pharmacokinetics of this combination in children with recurrent/refractory solid tumors; and to explore clinical activity in Ewing sarcoma (EWS) (NCT02116777). METHODS: Talazoparib (400-600 µg/m2 /dose, maximum daily dose 800-1000 µg) was administered q.d. or b.i.d. orally on day 1 followed by q.d. dosing concomitant with q.d. dosing of oral TMZ (20-55 mg/m2 /day) on days 2 to 6, every 28 days. RESULTS: Forty patients, aged 4 to 25 years, were enrolled. Talazoparib was increased to 600 µg/m2 /dose b.i.d. on day 1, and q.d. thereafter, with 20 mg/m2 /day of TMZ, without DLTs. TMZ was subsequently increased, during which dose-limiting neutropenia and thrombocytopenia occurred in two of three subjects at 55 mg/m2 /day, two of six subjects at 40 mg/m2 /day, and one of six subjects at 30 mg/m2 /day. During dose-finding, two of five EWS and four of 25 non-EWS subjects experienced prolonged stable disease (SD), and one subject with malignant glioma experienced a partial response. In phase 2, 0 of 10 EWS subjects experienced an objective response; two experienced prolonged SD. CONCLUSIONS: Talazoparib and low-dose TMZ are tolerated in children with recurrent/refractory solid tumors. Reversible neutropenia and thrombocytopenia were dose limiting. The RP2D is talazoparib 600 µg/m2 b.i.d. on day 1 followed by 600 µg/m2 q.d. on days 2 to 6 (daily maximum 1000 µg) in combination with temozolomide 30 mg/m2 /day on days 2 to 6. Antitumor activity was not observed in EWS, and limited antitumor activity was observed in central nervous system tumors.