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.

Epitaxial Kagome Thin Films as a Platform for Topological Flat Bands.

Systems with flat bands are ideal for studying strongly correlated electronic states and related phenomena. Among them, kagome-structured metals such as CoSn have been recognized as promising candidates due to the proximity between the flat bands and the Fermi level. A key next step will be to realize epitaxial kagome thin films with flat bands to enable tuning of the flat bands across the Fermi level via electrostatic gating or strain. Here, we report the band structures of epitaxial CoSn thin films grown directly on the insulating substrates. Flat bands are observed by using synchrotron-based angle-resolved photoemission spectroscopy (ARPES). The band structure is consistent with density functional theory (DFT) calculations, and the transport properties are quantitatively explained by the band structure and semiclassical transport theory. Our work paves the way to realize flat band-induced phenomena through fine-tuning of flat bands in kagome materials.

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.

Synthesis, Magnetic Properties, and Electronic Structure of Magnetic Topological Insulator MnBi2Se4.

The intrinsic magnetic topological insulators MnBi2Te4 and MnBi2Se4 support novel topological states related to symmetry breaking by magnetic order. Unlike MnBi2Te4, the study of MnBi2Se4 has been inhibited by the lack of bulk crystals, as the van der Waals (vdW) crystal is not the thermodynamic equilibrium phase. Here, we report the layer-by-layer synthesis of vdW MnBi2Se4 crystals using nonequilibrium molecular beam epitaxy. Atomic-resolution scanning transmission electron microscopy and scanning tunneling microscopy identify a well-ordered vdW crystal with septuple-layer base units. The magnetic properties agree with the predicted layered antiferromagnetic ordering but disagree with its predicted out-of-plane orientation. Instead, our samples exhibit an easy-plane anisotropy, which is explained by including dipole-dipole interactions. Angle-resolved photoemission spectroscopy reveals the gapless Dirac-like surface state, which demonstrates that MnBi2Se4 is a topological insulator above the magnetic-ordering temperature. These studies show that MnBi2Se4 is a promising candidate for exploring rich topological phases of layered antiferromagnetic topological insulators.

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.

The Lived Experiences of African American Centenarians: Narrative Exploration From a Life Course Theory Perspective.

BACKGROUND AND OBJECTIVES: The purpose of this study was to examine examples of inequality experienced by African American centenarians throughout the life course. RESEARCH DESIGN AND METHODS: Data were derived from 111 oral history transcripts from the Oklahoma Oral History Project. A subsample of all centenarians who identified as African American (n = 12) was utilized. An interpretive phenomenological approach and thematic analysis were used to review and code each narrative. RESULTS: Findings revealed 2 themes: (a) historical time and place with subthemes on family lineage, segregation, discrimination, and racism, and societal changes; (b) resource buffers against inequality with subthemes on multigenerational transmission, support and cultural connectedness, religion and spirituality, and self-reliance. Several centenarians acknowledged having grandparents who were slaves and others reported hidden family lineage due to complicated mixed-race dynamics. Work histories were affected by the historical era, with limited "domestic" positions for women. Experiences with education varied, with one centenarian commenting on his legacy of teaching at historically Black schools whereas others expressed difficulty in keeping up with school and regret for not finishing. DISCUSSION AND IMPLICATIONS: This is one of the first studies to explore life course experiences specific to African American centenarians. Centenarians in this study possessed survival skills and various support systems including community and religious/spiritual support, perhaps buffering against the impact of cumulative racial inequality.

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.