Dynamic inosinome profiles reveal novel patient stratification and gender-specific differences in glioblastoma

Is Metal Stabilization of the Leaving Group Required or Can Lysine Facilitate Phosphodiester Bond Cleavage in Nucleic Acids? A Computational Study of EndoV

Endonuclease V (EndoV) is a single-metal-dependent enzyme that repairs deaminated DNA nucleobases in cells by cleaving the phosphodiester bond, and this enzyme has proven to be a powerful tool in biotechnology and medicine. The catalytic mechanism used by EndoV must be understood to design new disease detection and therapeutic solutions and further exploit the enzyme in interdisciplinary applications. This study has used a mixed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) approach to compare eight distinct catalytic pathways and provides the first proposed mechanism for bacterial EndoV. The calculations demonstrate that mechanisms involving either direct or indirect metal coordination to the leaving group of the substrate previously proposed for other nucleases are unlikely for EndoV, regardless of the general base (histidine, aspartate, and substrate phosphate moiety). Instead, distinct catalytic pathways are characterized for EndoV that involve K139 stabilizing the leaving group, a metal-coordinated water stabilizing the transition structure, and either H214 or a substrate phosphate group activating the water nucleophile. In silico K139A and H214A mutational results support the newly proposed roles of these residues. Although this is a previously unseen combination of general base, general acid, and metal-binding architecture for a one-metal-dependent endonuclease, our proposed catalytic mechanisms are fully consistent with experimental kinetic, structural, and mutational data. In addition to substantiating a growing body of literature, suggesting that one metal is enough to catalyze P-O bond cleavage in nucleic acids, this new fundamental understanding of the catalytic function will promote the exploration of new and improved applications of EndoV.

Genome-wide perturbations of A-to-I RNA editing dysregulated circular RNAs promoting the development of cervical cancer

Cervical cancer, the second most common female malignant tumor, seriously threatens women's health and lives. Despite the availability of the HPV vaccine, effective treatment options for cervical cancer are still lacking. New research perspectives now clarify that RNA editing dysregulation and changes in circRNA expression are jointly involved in disease pathogenesis, so molecular changes associated with circRNA and RNA editing may provide clues for the development of new therapeutic strategies for cervical cancer. In this study, we designed a series of pipelines to identify and analyze dysregulated RNA editing events in circRNAs. Our findings indicate a decrease in A-to-I RNA editing levels in cervical cancer compared to normal tissues, and editing may influence the back-splicing process of circRNAs through structural modifications of Alu elements. Moreover, our research reveals that RNA editing could modulate circRNA biogenesis by influencing RNA binding protein (RBP) binding on a transcriptome-wide scale, as well as influence the expression and coding potential of circRNAs. Importantly, we identified three RNA editing sites that could serve as potential biomarkers. In summary, our study presents a comprehensive landscape of RNA editing perturbations in circRNAs, providing new insights into the complex relationship between RNA editing and circRNA dysregulation in cervical cancer.

Contribution of A-to-I RNA editing, M6A RNA Methylation, and Alternative Splicing to physiological brain aging and neurodegenerative diseases

Aging is a physiological and progressive phenomenon in all organisms' life cycle, characterized by the accumulation of degenerative processes triggered by several alterations within molecular pathways. These changes compromise cell fate, resulting in the loss of functions in tissues throughout the body, including the brain. Physiological brain aging has been linked to structural and functional alterations, as well as to an increased risk of neurodegenerative diseases. Post-transcriptional RNA modifications modulate mRNA coding properties, stability, translatability, expanding the coding capacity of the genome, and are involved in all cellular processes. Among mRNA post-transcriptional modifications, the A-to-I RNA editing, m6A RNA Methylation and Alternative Splicing play a critical role in all the phases of a neuronal cell life cycle and alterations in their mechanisms of action significantly contribute to aging and neurodegeneration. Here we review our current understanding of the contribution of A-to-I RNA editing, m6A RNA Methylation, and Alternative Splicing to physiological brain aging process and neurodegenerative diseases.

Extensive Recoding of the Neural Proteome in Cephalopods by RNA Editing

The coleoid cephalopods have the largest brains, and display the most complex behaviors, of all invertebrates. The molecular and cellular mechanisms that underlie these remarkable advancements remain largely unexplored. Early molecular cloning studies of squid ion channel transcripts uncovered an unusually large number of A→I RNA editing sites that recoded codons. Further cloning of other neural transcripts showed a similar pattern. The advent of deep-sequencing technologies and the associated bioinformatics allowed the mapping of RNA editing events across the entire neural transcriptomes of various cephalopods. The results were remarkable: They contained orders of magnitude more recoding editing sites than any other taxon. Although RNA editing sites are abundant in most multicellular metazoans, they rarely recode. In cephalopods, the majority of neural transcripts are recoded. Recent studies have focused on whether these events are adaptive, as well as other noncanonical aspects of cephalopod RNA editing.

The RNA editing landscape in acute myeloid leukemia reveals associations with disease mutations and clinical outcome

Several studies have documented aberrant RNA editing patterns across multiple tumors across large patient cohorts from The Cancer Genome Atlas (TCGA). However, studies on understanding the role of RNA editing in acute myeloid leukemia (AML) have been limited to smaller sample sizes. Using high throughput transcriptomic data from the TCGA, we demonstrated higher levels of editing as a predictor of poor outcome within the AML patient samples. Moreover, differential editing patterns were observed across individual AML genotypes. We also could demonstrate a negative association between the degree of editing and mRNA abundance for some transcripts, identifying the potential regulatory potential of RNA-editing in altering gene expression in AML. Further edQTL analysis suggests potential cis-regulatory mechanisms in RNA editing variation. Our work suggests a functional and regulatory role of RNA editing in the pathogenesis of AML and we extended our analysis to gain insight into the factors influencing altered levels of editing.

Profiling RNA Editing in Single Cells

RNA editing is a widespread molecular phenomenon occurring in a variety of organisms. In humans, it mainly involves the deamination of adenosine to inosine (A-to-I) in double-stranded RNAs by ADAR enzymes. A-to-I RNA editing has been investigated in different tissues as well as in diverse experimental and pathological conditions. By contrast, its biological role in single cells has not been explored in depth. Recent methodologies for cell sorting in combination with deep sequencing technologies have enabled the study of eukaryotic transcriptomes at single cell resolution, paving the way to the profiling of their epitranscriptomic dynamics.Here we describe a step-by-step protocol to detect and characterize A-to-I events occurring in publicly available single-cell RNAseq experiments from human alpha and beta pancreatic cells.

ADAR2 Protein Is Associated with Overall Survival in GBM Patients and Its Decrease Triggers the Anchorage-Independent Cell Growth Signature

Background: Epitranscriptomic mechanisms, such as A-to-I RNA editing mediated by ADAR deaminases, contribute to cancer heterogeneity and patients’ stratification. ADAR enzymes can change the sequence, structure, and expression of several RNAs, affecting cancer cell behavior. In glioblastoma, an overall decrease in ADAR2 RNA level/activity has been reported. However, no data on ADAR2 protein levels in GBM patient tissues are available; and most data are based on ADARs overexpression experiments. Methods: We performed IHC analysis on GBM tissues and correlated ADAR2 levels and patients’ overall survival. We silenced ADAR2 in GBM cells, studied cell behavior, and performed a gene expression/editing analysis. Results: GBM tissues do not all show a low/no ADAR2 level, as expected by previous studies. Although, different amounts of ADAR2 protein were observed in different patients, with a low level correlating with a poor patient outcome. Indeed, reducing the endogenous ADAR2 protein in GBM cells promotes cell proliferation and migration and changes the cell’s program to an anchorage-independent growth mode. In addition, deep-seq data and bioinformatics analysis indicated multiple RNAs are differently expressed/edited upon siADAR2. Conclusion: ADAR2 protein is an important deaminase in GBM and its amount correlates with patient prognosis.

ADAR3 activates NF-κB signaling and promotes glioblastoma cell resistance to temozolomide

The RNA binding protein ADAR3 is expressed exclusively in the brain and reported to have elevated expression in tumors of patients suffering from glioblastoma compared to adjacent brain tissue. Yet, other studies have indicated that glioblastoma tumors exhibit hemizygous deletions of the genomic region encompassing ADAR3 (10p15.3). As the molecular and cellular consequences of altered ADAR3 expression are largely unknown, here we directly examined the impacts of elevated ADAR3 in a glioblastoma cell line model. Transcriptome-wide sequencing revealed 641 differentially expressed genes between control and ADAR3-expressing U87-MG glioblastoma cells. A vast majority of these genes belong to pathways involved in glioblastoma progression and are regulated by NF-κB signaling. Biochemical and molecular analysis indicated that ADAR3-expressing U87-MG cells exhibit increased NF-κB activation, and treatment with an NF-κB inhibitor abrogated the impacts of ADAR3 on gene expression. Similarly, we found that increased cell survival of ADAR3-expressing cells to temozolomide, the preferred chemotherapeutic for glioblastoma, was due to increased NF-κB activity. Aberrant constitutive NF-κB activation is a common event in glioblastoma and can impact both tumor progression and resistance to treatment. Our results suggest that elevated ADAR3 promotes NF-κB activation and a gene expression program that provides a growth advantage to glioblastoma cells.

RNA binding by ADAR3 inhibits adenosine-to-inosine editing and promotes expression of immune response protein MAVS

Members of the ADAR family of double-stranded RNA–binding proteins regulate one of the most abundant RNA modifications in humans, the deamination of adenosine to inosine. Several transcriptome-wide studies have been carried out to identify RNA targets of the active deaminases ADAR1 and ADAR2. However, our understanding of ADAR3, the brain-specific deaminase-deficient ADAR family member, is limited to a few transcripts. In this study, we identified over 3300 transcripts bound by ADAR3 and observed that binding of ADAR3 correlated with reduced editing of over 400 sites in the glioblastoma transcriptome. We further investigated the impact of ADAR3 on gene regulation of the transcript that encodes MAVS, an essential protein in the innate immune response pathway. We observed reduced editing in the MAVS 3′ UTR in cells expressing increased ADAR3 or reduced ADAR1 suggesting ADAR3 acts as a negative regulator of ADAR1-mediated editing. While neither ADAR1 knockdown or ADAR3 overexpression affected MAVS mRNA expression, we demonstrate increased ADAR3 expression resulted in upregulation of MAVS protein expression. In addition, we created a novel genetic mutant of ADAR3 that exhibited enhanced RNA binding and MAVS upregulation compared with wildtype ADAR3. Interestingly, this ADAR3 mutant no longer repressed RNA editing, suggesting ADAR3 has a unique regulatory role beyond altering editing levels. Altogether, this study provides the first global view of ADAR3-bound RNAs in glioblastoma cells and identifies both a role for ADAR3 in repressing ADAR1-mediated editing and an RNA-binding dependent function of ADAR3 in regulating MAVS expression.

RNA Editing in Glioma as a Sexually Dimorphic Prognostic Factor That Affects mRNA Abundance in Fatty Acid Metabolism and Inflammation Pathways

RNA editing alters the nucleotide sequence and has been associated with cancer progression. However, little is known about its prognostic and regulatory roles in glioma, one of the most common types of primary brain tumors. We characterized and analyzed RNA editomes of glioblastoma and isocitrate dehydrogenase mutated (IDH-MUT) gliomas from The Cancer Genome Atlas and the Chinese Glioma Genome Atlas (CGGA). We showed that editing change during glioma progression was another layer of molecular alterations and that editing profiles predicted the prognosis of glioblastoma and IDH-MUT gliomas in a sex-dependent manner. Hyper-editing was associated with poor survival in females but better survival in males. Moreover, noncoding editing events impacted mRNA abundance of the host genes. Genes associated with inflammatory response (e.g., EIF2AK2, a key mediator of innate immunity) and fatty acid oxidation (e.g., acyl-CoA oxidase 1, the rate-limiting enzyme in fatty acid β-oxidation) were editing-regulated and associated with glioma progression. The above findings were further validated in CGGA samples. Establishment of the prognostic and regulatory roles of RNA editing in glioma holds promise for developing editing-based therapeutic strategies against glioma progression. Furthermore, sexual dimorphism at the epitranscriptional level highlights the importance of developing sex-specific treatments for glioma.

A Systematic Review of Common and Brain-Disease-Specific RNA Editing Alterations Providing Novel Insights into Neurological and Neurodegenerative Disease Manifestations

RNA editing contributes to transcriptome diversification through RNA modifications in relation to genome-encoded information (RNA–DNA differences, RDDs). The deamination of Adenosine (A) to Inosine (I) or Cytidine (C) to Uridine (U) is the most common type of mammalian RNA editing. It occurs as a nuclear co- and/or post-transcriptional event catalyzed by ADARs (Adenosine deaminases acting on RNA) and APOBECs (apolipoprotein B mRNA editing enzyme catalytic polypeptide-like genes). RNA editing may modify the structure, stability, and processing of a transcript. This review focuses on RNA editing in psychiatric, neurological, neurodegenerative (NDs), and autoimmune brain disorders in humans and rodent models. We discuss targeted studies that focus on RNA editing in specific neuron-enriched transcripts with well-established functions in neuronal activity, and transcriptome-wide studies, enabled by recent technological advances. We provide comparative editome analyses between human disease and corresponding animal models. Data suggest RNA editing to be an emerging mechanism in disease development, displaying common and disease-specific patterns. Commonly edited RNAs represent potential disease-associated targets for therapeutic and diagnostic values. Currently available data are primarily descriptive, calling for additional research to expand global editing profiles and to provide disease mechanistic insights. The potential use of RNA editing events as disease biomarkers and available tools for RNA editing identification, classification, ranking, and functional characterization that are being developed will enable comprehensive analyses for a better understanding of disease(s) pathogenesis and potential cures.

Landscape of adenosine-to-inosine RNA recoding across human tissues

RNA editing by adenosine deaminases changes the information encoded in the mRNA from its genomic blueprint. Editing of protein-coding sequences can introduce novel, functionally distinct, protein isoforms and diversify the proteome. The functional importance of a few recoding sites has been appreciated for decades. However, systematic methods to uncover these sites perform poorly, and the full repertoire of recoding in human and other mammals is unknown. Here we present a new detection approach, and analyze 9125 GTEx RNA-seq samples, to produce a highly-accurate atlas of 1517 editing sites within the coding region and their editing levels across human tissues. Single-cell RNA-seq data shows protein recoding contributes to the variability across cell subpopulations. Most highly edited sites are evolutionary conserved in non-primate mammals, attesting for adaptation. This comprehensive set can facilitate understanding of the role of recoding in human physiology and diseases.

Gabay et al. provide a highly-accurate atlas of recoding by A-to-I RNA editing in human, profiled across tissues and cell subpopulations. Most highly edited sites are evolutionary conserved in non-primate mammals, attesting for adaptation.

ADAR1-mediated RNA editing links ganglioside catabolism to glioblastoma stem cell maintenance

Glioblastoma (GBM) is the most common and lethal primary malignant brain tumor, containing GBM stem cells (GSCs) that contribute to therapeutic resistance and relapse. Exposing potential GSC vulnerabilities may provide therapeutic strategies against GBM. Here, we interrogated the role of Adenosine-to-Inosine (A-to-I) RNA editing mediated by ADAR1 (adenosine deaminase acting on RNA 1) in GSCs and found that both ADAR1 and global RNA editomes were elevated in GSCs compared to normal neural stem cells (NSCs). ADAR1 inactivation or blocking the upstream JAK/STAT pathway through TYK2 inhibition impaired GSC self-renewal and stemness. Downstream of ADAR1, RNA editing of the 3'UTR of GM2A, a key ganglioside catabolism activator, proved to be critical, as interfering with ganglioside catabolism showed similar functional impact on GSCs as ADAR1 disruption. These findings reveal RNA editing links ganglioside catabolism to GSC self-renewal and stemness, exposing a potential vulnerability of GBM for therapeutic intervention.

Potential Effect of Post-Transcriptional Substitutions of Tyrosine for Cysteine Residues on Transformation of Amyloidogenic Proteins

The review considers the reasons and consequences of post-transcriptional tyrosine substitutions for cysteine residues. Main attention is paid to the Tyr/Cys substitutions that arise during gene expression in bacterial systems at the stage of protein translation as a result of misrecognition of the similar mRNA codons. Notably, translation errors generally occur relatively rarely - from 10^-4 to 10^-3 errors per codon for E. coli cells, but in some cases the error rate increases significantly. For example, this is typical for certain pairs of codons, when the culture conditions change or in the presence of antibiotics. Thus, with overproduction of the recombinant human alpha-synuclein in E. coli cells, the content of the mutant form with the replacement of Tyr136 (UAC codon) with a cysteine residue (UGC codon) can reach 50%. Possible reasons for the increased production of alpha-synuclein with the Tyr136Cys substitution are considered, as well as consequences of the presence of mutant forms in preparations of amyloidogenic proteins when studying their pathological transformation in vitro. A separate section is devoted to the Tyr/Cys substitutions occurring due to mRNA editing by adenosine deaminases, which is typical for eukaryotic organisms, and the possible role of this process in the amyloid transformation of proteins associated with neurodegenerative diseases.

Epitranscriptomic Modifications Modulate Normal and Pathological Functions in CNS

RNA is more than just a combination of four genetically encoded nucleobases as it carries extra information in the form of epitranscriptomic modifications. Diverse chemical groups attach covalently to RNA to enhance the plasticity of cellular transcriptome. The reversible and dynamic nature of epitranscriptomic modifications allows RNAs to achieve rapid and context-specific gene regulation. Dedicated cellular machinery comprising of writers, erasers, and readers drives the epitranscriptomic signaling. Epitranscriptomic modifications control crucial steps of mRNA metabolism such as splicing, export, localization, stability, degradation, and translation. The majority of the epitranscriptomic modifications are highly abundant in the brain and contribute to activity-dependent gene expression. Thus, they regulate the vital physiological processes of the brain, such as synaptic plasticity, neurogenesis, and stress response. Furthermore, epitranscriptomic alterations influence the progression of several neurologic disorders. This review discussed the molecular mechanisms of epitranscriptomic regulation in neurodevelopmental and neuropathological conditions with the goal to identify novel therapeutic targets.

Genome-wide perturbations of Alu expression and Alu-associated post-transcriptional regulations distinguish oligodendroglioma from other gliomas

Alu is a primate-specific repeat element in the human genome and has been increasingly appreciated as a regulatory element in many biological processes. But the appreciation of Alu has been limited in tumorigenesis, especially for brain tumor. To investigate the relevance of Alu to the gliomagenesis, we studied Alu element-associated post-transcriptional processes and the RNA expression of the element by performing RNA-seq for a total of 41 pairs of neurotypical and diverse glioma brain tissues. We find that A-to-I editing and circular RNA levels, as well as Alu RNA expression, are decreased overall in gliomas, compared to normal tissue. Interestingly, grade 2 oligodendrogliomas are least affected in A-to-I editing and circular RNA levels among gliomas, whereas they have a higher proportion of down-regulated Alu subfamilies, compared to the other gliomas. These findings collectively imply a unique pattern of Alu-associated transcriptomes in grade 2 oligodendroglioma, providing an insight to gliomagenesis from the perspective of an evolutionary genetic element.

The A to I editing landscape in melanoma and its relation to clinical outcome

RNA editing refers to non-transient RNA modifications that occur after transcription and prior to translation by the ribosomes. RNA editing is more widespread in cancer cells than in non-transformed cells and is associated with tumorigenesis of various cancer tissues. However, RNA editing can also generate neo-antigens that expose tumour cells to host immunosurveillance. Global RNA editing in melanoma and its relevance to clinical outcome currently remain poorly characterized. The present study compared RNA editing as well as gene expression in tumour cell lines from melanoma patients of short or long metastasis-free survival, patients relapsing or not after immuno- and targeted therapy and tumours harbouring BRAF or NRAS mutations. Overall, our results showed that NTRK gene expression can be a marker of resistance to BRAF and MEK inhibition and gives some insights of candidate genes as potential biomarkers. In addition, this study revealed an increase in Adenosine-to-Inosine editing in Alu regions and in non-repetitive regions, including the hyperediting of the MOK and DZIP3 genes in relapsed tumour samples during targeted therapy and of the ZBTB11 gene in NRAS mutated melanoma cells. Therefore, RNA editing could be a promising tool for identifying predictive markers, tumour neoantigens and targetable pathways that could help in preventing relapses during immuno- or targeted therapies.

A survey of RNA editing at single-cell resolution links interneurons to schizophrenia and autism

Conversion of adenosine to inosine in RNA by ADAR enzymes, termed "RNA editing," is essential for healthy brain development. Editing is dysregulated in neuropsychiatric diseases, but has not yet been investigated at scale at the level of individual neurons. We quantified RNA editing sites in nuclear transcriptomes of 3055 neurons from six cortical regions of a neurotypical female donor, and found 41,930 sites present in at least ten nuclei. Most sites were located within Alu repeats in introns or 3' UTRs, and approximately 80% were cataloged in public RNA editing databases. We identified 9285 putative novel editing sites, 29% of which were also detectable in unrelated donors. Intersection with results from bulk RNA-seq studies provided cell-type and spatial context for 1730 sites that are differentially edited in schizophrenic brain donors, and 910 such sites in autistic donors. Autism-related genes were also enriched with editing sites predicted to modify RNA structure. Inhibitory neurons showed higher overall transcriptome editing than excitatory neurons, and the highest editing rates were observed in the frontal cortex. We used generalized linear models to identify differentially edited sites and genes between cell types. Twenty nine genes were preferentially edited in excitatory neurons, and 43 genes were edited more heavily in inhibitory neurons, including RBFOX1, its target genes, and genes in the autism-associated Prader-Willi locus (15q11). The abundance of SNORD115/116 genes from locus 15q11 was positively associated with editing activity across the transcriptome. We contend that insufficient editing of autism-related genes in inhibitory neurons may contribute to the specific perturbation of those cells in autism.

Differential and Common Signatures of miRNA Expression and Methylation in Childhood Central Nervous System Malignancies: An Experimental and Computational Approach

Epigenetic modifications are considered of utmost significance for tumor ontogenesis and progression. Especially, it has been found that miRNA expression, as well as DNA methylation plays a significant role in central nervous system tumors during childhood. A total of 49 resected brain tumors from children were used for further analysis. DNA methylation was identified with methylation-specific MLPA and, in particular, for the tumor suppressor genes CASP8, RASSF1, MGMT, MSH6, GATA5, ATM1, TP53, and CADM1. miRNAs were identified with microarray screening, as well as selected samples, were tested for their mRNA expression levels. CASP8, RASSF1 were the most frequently methylated genes in all tumor samples. Simultaneous methylation of genes manifested significant results with respect to tumor staging, tumor type, and the differentiation of tumor and control samples. There was no significant dependence observed with the methylation of one gene promoter, rather with the simultaneous presence of all detected methylated genes' promoters. miRNA expression was found to be correlated to gene methylation. Epigenetic regulation appears to be of major importance in tumor progression and pathophysiology, making it an imperative field of study.

ADAR-mediated RNA editing of DNA:RNA hybrids is required for DNA double strand break repair

The maintenance of genomic stability requires the coordination of multiple cellular tasks upon the appearance of DNA lesions. RNA editing, the post-transcriptional sequence alteration of RNA, has a profound effect on cell homeostasis, but its implication in the response to DNA damage was not previously explored. Here we show that, in response to DNA breaks, an overall change of the Adenosine-to-Inosine RNA editing is observed, a phenomenon we call the RNA Editing DAmage Response (REDAR). REDAR relies on the checkpoint kinase ATR and the recombination factor CtIP. Moreover, depletion of the RNA editing enzyme ADAR2 renders cells hypersensitive to genotoxic agents, increases genomic instability and hampers homologous recombination by impairing DNA resection. Such a role of ADAR2 in DNA repair goes beyond the recoding of specific transcripts, but depends on ADAR2 editing DNA:RNA hybrids to ease their dissolution.

Direct Immunodetection of Global A-to-I RNA Editing Activity with a Chemiluminescent Bioassay

Adenosine-to-inosine (A-to-I) editing is a conserved eukaryotic RNA modification that contributes to development, immune response, and overall cellular function. Here, we utilize Endonuclease V (EndoV), which binds specifically to inosine in RNA, to develop an EndoV-linked immunosorbency assay (EndoVLISA) as a rapid, plate-based chemiluminescent method for measuring global A-to-I editing signatures in cellular RNA. We first optimize and validate our assay with chemically synthesized oligonucleotides. We then demonstrate rapid detection of inosine content in treated cell lines, demonstrating equivalent performance against current standard RNA-seq approaches. Lastly, we deploy our EndoVLISA for profiling differential A-to-I RNA editing signatures in normal and diseased human tissue, illustrating the utility of our platform as a diagnostic bioassay. Together, the EndoVLISA method is cost-effective, straightforward, and utilizes common laboratory equipment, offering a highly accessible new approach for studying A-to-I editing. Moreover, the multi-well plate format makes this the first assay amenable for direct high-throughput quantification of A-to-I editing for applications in disease detection and drug development.

Identification of A-to-I RNA editing profiles and their clinical relevance in lung adenocarcinoma

Adenosine-to-inosine (A-to-I) RNA editing is a widespread posttranscriptional modification that has been shown to play an important role in tumorigenesis. Here, we evaluated a total of 19,316 RNA editing sites in the tissues of 80 lung adenocarcinoma (LUAD) patients from our Nanjing Lung Cancer Cohort (NJLCC) and 486 LUAD patients from the TCGA database. The global RNA editing level was significantly increased in tumor tissues and was highly heterogeneous across patients. The high RNA editing level in tumors was attributed to both RNA (ADAR1 expression) and DNA alterations (mutation load). Consensus clustering on RNA editing sites revealed a new molecular subtype (EC3) that was associated with the poorest prognosis of LUAD patients. Importantly, the new classification was independent of classic molecular subtypes based on gene expression or DNA methylation. We further proposed a simplified model including eight RNA editing sites to accurately distinguish the EC3 subtype in our patients. The model was further validated in the TCGA dataset and had an area under the curve (AUC) of the receiver operating characteristic curve of 0.93 (95%CI: 0.91-0.95). In addition, we found that LUAD cell lines with the EC3 subtype were sensitive to four chemotherapy drugs. These findings highlighted the importance of RNA editing events in the tumorigenesis of LUAD and provided insight into the application of RNA editing in the molecular subtyping and clinical treatment of cancer.

Regulation of RNA editing by intracellular acidification

The hydrolytic deamination of adenosine-to-inosine (A-to-I) by RNA editing is a widespread post-transcriptional modification catalyzed by the adenosine deaminase acting on RNA (ADAR) family of proteins. ADAR-mediated RNA editing modulates cellular pathways involved in innate immunity, RNA splicing, RNA interference, and protein recoding, and has been investigated as a strategy for therapeutic intervention of genetic disorders. Despite advances in basic and translational research, the mechanisms regulating RNA editing are poorly understood. Though several trans-acting regulators of editing have been shown to modulate ADAR protein expression, previous studies have not identified factors that modulate ADAR catalytic activity. Here, we show that RNA editing increases upon intracellular acidification, and that these effects are predominantly explained by both enhanced ADAR base-flipping and deamination rate at acidic pH. We also show that the extent of RNA editing increases with the reduction in pH associated with conditions of cellular hypoxia.

A mark of disease: how mRNA modifications shape genetic and acquired pathologies

RNA modifications have recently emerged as a widespread and complex facet of gene expression regulation. Counting more than 170 distinct chemical modifications with far-reaching implications for RNA fate, they are collectively referred to as the epitranscriptome. These modifications can occur in all RNA species, including messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs). In mRNAs the deposition, removal, and recognition of chemical marks by writers, erasers and readers influence their structure, localization, stability, and translation. In turn, this modulates key molecular and cellular processes such as RNA metabolism, cell cycle, apoptosis, and others. Unsurprisingly, given their relevance for cellular and organismal functions, alterations of epitranscriptomic marks have been observed in a broad range of human diseases, including cancer, neurological and metabolic disorders. Here, we will review the major types of mRNA modifications and editing processes in conjunction with the enzymes involved in their metabolism and describe their impact on human diseases. We present the current knowledge in an updated catalog. We will also discuss the emerging evidence on the crosstalk of epitranscriptomic marks and what this interplay could imply for the dynamics of mRNA modifications. Understanding how this complex regulatory layer can affect the course of human pathologies will ultimately lead to its exploitation toward novel epitranscriptomic therapeutic strategies.

REDIportal: millions of novel A-to-I RNA editing events from thousands of RNAseq experiments

RNA editing is a relevant epitranscriptome phenomenon able to increase the transcriptome and proteome diversity of eukaryotic organisms. ADAR mediated RNA editing is widespread in humans in which millions of A-to-I changes modify thousands of primary transcripts. RNA editing has pivotal roles in the regulation of gene expression or modulation of the innate immune response or functioning of several neurotransmitter receptors. Massive transcriptome sequencing has fostered the research in this field. Nonetheless, different aspects of the RNA editing biology are still unknown and need to be elucidated. To support the study of A-to-I RNA editing we have updated our REDIportal catalogue raising its content to about 16 millions of events detected in 9642 human RNAseq samples from the GTEx project by using a dedicated pipeline based on the HPC version of the REDItools software. REDIportal now allows searches at sample level, provides overviews of RNA editing profiles per each RNAseq experiment, implements a Gene View module to look at individual events in their genic context and hosts the CLAIRE database. Starting from this novel version, REDIportal will start collecting non-human RNA editing changes for comparative genomics investigations. The database is freely available at http://srv00.recas.ba.infn.it/atlas/index.html.

ADAR1 is a new target of METTL3 and plays a pro-oncogenic role in glioblastoma by an editing-independent mechanism

Background

N⁶-methyladenosine (m6A) and adenosine-to-inosine (A-to-I) RNA editing are two of the most abundant RNA modification events affecting adenosines in mammals. Both these RNA modifications determine mRNA fate and play a pivotal role in tumor development and progression.

Results

Here, we show that METTL3, upregulated in glioblastoma, methylates ADAR1 mRNA and increases its protein level leading to a pro-tumorigenic mechanism connecting METTL3, YTHDF1, and ADAR1. We show that ADAR1 plays a cancer-promoting role independently of its deaminase activity by binding CDK2 mRNA, underlining the importance of ADARs as essential RNA-binding proteins for cell homeostasis as well as cancer progression. Additionally, we show that ADAR1 knockdown is sufficient to strongly inhibit glioblastoma growth in vivo.

Conclusions

Hence, our findings underscore METTL3/ADAR1 axis as a novel crucial pathway in cancer progression that connects m6A and A-to-I editing post-transcriptional events.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02271-9.

Gender and race interact to influence survival disparities in head and neck cancer

Gender and race disparities in head and neck squamous cell carcinoma (HSNCC) survival are independently well documented, but no prior studies have examined the joint effect of these factors on HSNCC outcomes. We aim to comprehensively estimate the effect of gender and race on overall survival in HNSCC. We constructed a retrospective cohort from the National Cancer Database for primary HNSCC of the larynx, hypopharynx, oral cavity, and oropharynx from 2010 to 2015. We used Kaplan-Meier curves and Cox proportional hazards regressions to calculate hazard ratios adjusting for treatment type, age, insurance, staging classifications, and comorbidities. Oral cavity cancer was significantly more common among Hispanic and White females compared to other sites. Female non-oropharyngeal HNSCC cases had better five-year overall survival than males (56.3% versus 54.4%, respectively), though Black females (52.8%) had poorer survival than both White (56.2%) and Hispanic (57.9%) males. There were significant differences in oropharyngeal cancer by HPV status. Notably, Black females with HPV-positive oropharyngeal OPSCC had far worse survival than any other race and gender group. These results persisted even when adjusting for potential mediating factors. Clearly gender is a significant prognosticator for HNSCC and has meaningful interactions with race. The distinct site distributions across gender and race reveal important insights into HNSCC among females. Taking into account these gender disparities while considering race is essential to providing appropriate care to head and neck patients and accurately counselling these individuals on prognosis and outcomes.

De Novo A-to-I RNA Editing Discovery in lncRNA

Simple Summary

Long non-coding RNAs are emerging as key regulators of gene expression at both transcriptional and translational levels, and their alterations (in expression or sequence) are linked to tumorigenesis and tumor progression. RNA editing has the unique ability to change the RNA sequence without altering the integrity or sequence of genomic DNA, with adenosine to inosine (A-to-I) RNA editing being the most common event in humans. With the ability to change the genetic information after transcription, RNA editing is an essential player in the transcriptome and proteome enrichment; however, when deregulated, it can contribute to cell transformation. In this article, we performed the first deep de novo editing survey in lncRNA, demonstrating that RNA editing is a pervasive phenomenon involving lncRNAs important in the brain and brain cancer. Our study will open a new field of research in which the interplay between lncRNA and RNA editing can add novel insights into cancer.

Abstract

Background: Adenosine to inosine (A-to-I) RNA editing is the most frequent editing event in humans. It converts adenosine to inosine in double-stranded RNA regions (in coding and non-coding RNAs) through the action of the adenosine deaminase acting on RNA (ADAR) enzymes. Long non-coding RNAs, particularly abundant in the brain, account for a large fraction of the human transcriptome, and their important regulatory role is becoming progressively evident in both normal and transformed cells. Results: Herein, we present a bioinformatic analysis to generate a comprehensive inosinome picture in long non-coding RNAs (lncRNAs), using an ad hoc index and searching for de novo editing events in the normal brain cortex as well as in glioblastoma, a highly aggressive human brain cancer. We discovered >10,000 new sites and 335 novel lncRNAs that undergo editing, never reported before. We found a generalized downregulation of editing at multiple lncRNA sites in glioblastoma samples when compared to the normal brain cortex. Conclusion: Overall, our study discloses a novel layer of complexity that controls lncRNAs in the brain and brain cancer.

Statistical inference of differential RNA-editing sites from RNA-sequencing data by hierarchical modeling

MOTIVATION

RNA-sequencing (RNA-seq) enables global identification of RNA-editing sites in biological systems and disease. A salient step in many studies is to identify editing sites that statistically associate with treatment (e.g. case versus control) or covary with biological factors, such as age. However, RNA-seq has technical features that incumbent tests (e.g. t-test and linear regression) do not consider, which can lead to false positives and false negatives.

RESULTS

In this study, we demonstrate the limitations of currently used tests and introduce the method, RNA-editing tests (REDITs), a suite of tests that employ beta-binomial models to identify differential RNA editing. The tests in REDITs have higher sensitivity than other tests, while also maintaining the type I error (false positive) rate at the nominal level. Applied to the GTEx dataset, we unveil RNA-editing changes associated with age and gender, and differential recoding profiles between brain regions.

AVAILABILITY AND IMPLEMENTATION

REDITs are implemented as functions in R and freely available for download at https://github.com/gxiaolab/REDITs. The repository also provides a code example for leveraging parallelization using multiple cores.

Global RNA editome landscape discovers reduced RNA editing in glioma: loss of editing of gamma-amino butyric acid receptor alpha subunit 3 (GABRA3) favors glioma migration and invasion

Background

Gliomas are the most common and lethal type of intracranial tumors. With the current treatment regime, the median survival of patients with grade IV glioma (glioblastoma/GBM) remains at 14-16 months. RNA editing modifies the function and regulation of transcripts. The development of glial tumors may be caused by altered RNA editing events.

Methods

In this study, we uncover the global RNA editome landscape of glioma patients from RNA-seq data of control, lower grade glioma (LGG) and GBM samples (n = 1,083).

Results

A-to-I editing events were found to comprise 80% of the total editing events of which 96% were located in the Alu regions. The total RNA editing events were found to be reduced in glioma compared to control samples. More specifically, we found Gamma-aminobutyric acid type A receptor alpha3 (GABRA3) to be edited (c.1026 A-to-G; pI343M) in 73% (editing ratio 0.8) of control samples compared to LGG (28.96%; 0.47) and GBM (5.2%; 0.53) samples. GABRA3 transcript level was found to be downregulated in glioma compared to control in a grade-specific manner with GBMs having the lowest level of the transcript. Further, GABRA3 transcripts were observed to be higher in edited compared to unedited glioma samples. The transcript and protein levels of exogenously expressed gene were found to be higher for edited compared to unedited GABRA3 in glioma cells. Further, exogenously expressed edited GABRA3 inhibited migration and invasion of glioma cells efficiently but not the unedited GABRA3.

Conclusion

Collectively, our study discovered a reduction in RNA editing during glioma development. We further demonstrate that elevated RNA editing maintains a high level of GABRA3 RNA and protein in normal glial cells which provides a less migratory environment for the normal functioning of the brain. In contrast, the reduction in GABRA3 protein levels, due to lower stability of unedited RNA, results in the loss of function which confers an aggressive phenotype to GBM tumor.

Altered Regulation of adipomiR Editing with Aging

Adipose dysfunction with aging increases risk to insulin resistance and other chronic metabolic diseases. We previously showed functional changes in microRNAs involved in pre-adipocyte differentiation with aging resulting in adipose dysfunction. However, the mechanisms leading to this dysfunction in microRNAs in adipose tissue (adipomiRs) during aging are not well understood. We determined the longitudinal changes in expression of adipomiRs and studied their regulatory mechanisms, such as miRNA biogenesis and editing, in an aging rodent model, with Fischer344 × Brown-Norway hybrid rats at ages ranging from 3 to 30 months (male/females, n > 8). Expression of adipomiRs and their edited forms were determined by small-RNA sequencing. RT-qPCR was used to measure the mRNA expression of biogenesis and editing enzymes. Sanger sequencing was used to validate editing with aging. Differential expression of adipomiRs involved in adipocyte differentiation and insulin signaling was altered with aging. Sex- and age-specific changes in edited adipomiRs were observed. An increase in miRNA biogenesis and editing enzymes (ADARs and their splice variants) were observed with increasing age, more so in female than male rats. The adipose dysfunction observed with age is attributed to differences in editing of adipomiRs, suggesting a novel regulatory pathway in aging.

RNA modifications in brain tumorigenesis

RNA modifications are emerging as critical regulators in cancer biology, thanks to their ability to influence gene expression and the predominant protein isoforms expressed during cell proliferation, migration, and other pro-oncogenic properties. The reversibility and dynamic nature of post-transcriptional RNA modifications allow cells to quickly adapt to microenvironmental changes. Recent literature has revealed that the deregulation of RNA modifications can promote a plethora of developmental diseases, including tumorigenesis. In this review, we will focus on four key post-transcriptional RNA modifications which have been identified as contributors to the pathogenesis of brain tumors: m6A, alternative polyadenylation, alternative splicing and adenosine to inosine modifications. In addition to the role of RNA modifications in brain tumor progression, we will also discuss potential opportunities to target these processes to improve the dismal prognosis for brain tumors.

Sex differences in cancer mechanisms

We now know that cancer is many different diseases, with great variation even within a single histological subtype. With the current emphasis on developing personalized approaches to cancer treatment, it is astonishing that we have not yet systematically incorporated the biology of sex differences into our paradigms for laboratory and clinical cancer research. While some sex differences in cancer arise through the actions of circulating sex hormones, other sex differences are independent of estrogen, testosterone, or progesterone levels. Instead, these differences are the result of sexual differentiation, a process that involves genetic and epigenetic mechanisms, in addition to acute sex hormone actions. Sexual differentiation begins with fertilization and continues beyond menopause. It affects virtually every body system, resulting in marked sex differences in such areas as growth, lifespan, metabolism, and immunity, all of which can impact on cancer progression, treatment response, and survival. These organismal level differences have correlates at the cellular level, and thus, males and females can fundamentally differ in their protections and vulnerabilities to cancer, from cellular transformation through all stages of progression, spread, and response to treatment. Our goal in this review is to cover some of the robust sex differences that exist in core cancer pathways and to make the case for inclusion of sex as a biological variable in all laboratory and clinical cancer research. We finish with a discussion of lab- and clinic-based experimental design that should be used when testing whether sex matters and the appropriate statistical models to apply in data analysis for rigorous evaluations of potential sex effects. It is our goal to facilitate the evaluation of sex differences in cancer in order to improve outcomes for all patients.

Quantifying RNA Editing in Deep Transcriptome Datasets

Massive transcriptome sequencing through the RNAseq technology has enabled quantitative transcriptome-wide investigation of co-/post-transcriptional mechanisms such as alternative splicing and RNA editing. The latter is abundant in human transcriptomes in which million adenosines are deaminated into inosines by the ADAR enzymes. RNA editing modulates the innate immune response and its deregulation has been associated with different human diseases including autoimmune and inflammatory pathologies, neurodegenerative and psychiatric disorders, and tumors. Accurate profiling of RNA editing using deep transcriptome data is still a challenge, and the results depend strongly on processing and alignment steps taken. Accurate calling of the inosinome repertoire, however, is required to reliably quantify RNA editing and, in turn, investigate its biological and functional role across multiple samples. Using real RNAseq data, we demonstrate the impact of different bioinformatics steps on RNA editing detection and describe the main metrics to quantify its level of activity.

Selective Enrichment of A-to-I Edited Transcripts from Cellular RNA Using Endonuclease V

Creating accurate maps of A-to-I RNA editing activity is vital to improving our understanding of the biological role of this process and harnessing it as a signal for disease diagnosis. Current RNA sequencing techniques are susceptible to random sampling limitations due to the complexity of the transcriptome and require large amounts of RNA material, specialized instrumentation, and high read counts to accurately interrogate A-to-I editing sites. To address these challenges, we show that Escherichia coli Endonuclease V (eEndoV), an inosine-cleaving enzyme, can be repurposed to bind and isolate A-to-I edited transcripts from cellular RNA. While Mg²⁺ enables eEndoV to catalyze RNA cleavage, we show that similar levels of Ca²⁺ instead promote binding of inosine without cleavage and thus enable high affinity capture of inosine in RNA. We leverage this capability to demonstrate EndoVIPER-seq (Endonuclease V inosine precipitation enrichment sequencing) as a facile and effective method to enrich A-to-I edited transcripts prior to RNA-seq, producing significant increases in the coverage and detection of identified editing sites. We envision the use of this approach as a straightforward and cost-effective strategy to improve the epitranscriptomic informational density of RNA samples, facilitating a deeper understanding of the functional roles of A-to-I editing.

AKT-Dependent Phosphorylation of ADAR1p110 and ADAR2 Represents a New and Important Link Between Cell Signaling and RNA Editing

RNA editing is a process by which nascent RNA transcripts are covalently modified, thus enhancing the complexity of the transcriptome. The most common modifications are deaminations of adenosine to inosine at sites of complex RNA secondary structure, a process that is carried out by the adenosine deaminase acting on double-strand RNA (ADAR) family of RNA editases. Although much has been learned about the ADAR family members since their discovery, very little information on their post-transcriptional regulation has been reported. Similar to most proteins, the ADAR family members are post-translationally modified at multiple sites. We recently reported that members of the AKT kinase family directly phosphorylate ADAR1p110 and ADAR2 on a conserved threonine within the catalytic domain of the protein. Phosphorylation was observed to differentially inhibit the enzymatic activity of the ADAR proteins toward known RNA substrates. The direct downstream involvement of the AKT kinases in multiple major signaling pathways associated with cell survival, growth, glucose metabolism (insulin signaling), and differentiation is well established; thus, the AKT kinases represent a link between ADAR-dependent A-to-I editing and major signal transduction pathways that are necessary for cell maintenance and development.

Investigating RNA editing in deep transcriptome datasets with REDItools and REDIportal

RNA editing is a widespread post-transcriptional mechanism able to modify transcripts through insertions/deletions or base substitutions. It is prominent in mammals, in which millions of adenosines are deaminated to inosines by members of the ADAR family of enzymes. A-to-I RNA editing has a plethora of biological functions, but its detection in large-scale transcriptome datasets is still an unsolved computational task. To this aim, we developed REDItools, the first software package devoted to the RNA editing profiling in RNA-sequencing (RNAseq) data. It has been successfully used in human transcriptomes, proving the tissue and cell type specificity of RNA editing as well as its pervasive nature. Outcomes from large-scale REDItools analyses on human RNAseq data have been collected in our specialized REDIportal database, containing more than 4.5 million events. Here we describe in detail two bioinformatic procedures based on our computational resources, REDItools and REDIportal. In the first procedure, we outline a workflow to detect RNA editing in the human cell line NA12878, for which transcriptome and whole genome data are available. In the second procedure, we show how to identify dysregulated editing at specific recoding sites in post-mortem brain samples of Huntington disease donors. On a 64-bit computer running Linux with ≥32 GB of random-access memory (RAM), both procedures should take ~76 h, using 4 to 24 cores. Our protocols have been designed to investigate RNA editing in different organisms with available transcriptomic and/or genomic reads. Scripts to complete both procedures and a docker image are available at https://github.com/BioinfoUNIBA/REDItools.

ADAR, the carcinogenesis mechanisms of ADAR and related clinical applications

Adenosine deaminases acting on RNA (ADARs) catalyze the conversion of adenosine (A) to inosine (I) in double-stranded RNA, which can change the codons after transcription. Abnormal ADAR editing is present in a variety of cancers. However, the study of the biological effects of ADARs in cancer is not very deep. Here, we review current important ADAR-mediated editing events, related carcinogenic mechanisms and applications in clinical medicine. Further exploration in ADARs can provide a new direction for cancer treatment.

Genome-wide quantification of ADAR adenosine-to-inosine RNA editing activity

Adenosine-to-inosine (A-to-I) RNA editing by the adenosine deaminase that acts on RNA (ADAR) enzymes is a common RNA modification, preventing false activation of the innate immune system by endogenous double-stranded RNAs. Methods for quantification of ADAR activity are sought after, due to an increasing interest in the role of ADARs in cancer and autoimmune disorders, as well as attempts to harness the ADAR enzymes for RNA engineering. Here, we present the Alu editing index (AEI), a robust and simple-to-use computational tool devised for this purpose. We describe its properties and demonstrate its superiority to current quantification methods of ADAR activity. The AEI is used to map global editing across a large dataset of healthy human samples and identify putative regulators of ADAR, as well as previously unknown factors affecting the observed Alu editing levels. These should be taken into account in future comparative studies of ADAR activity. The AEI tool is available at https://github.com/a2iEditing/RNAEditingIndexer.

Radiosensitization and a Less Aggressive Phenotype of Human Malignant Glioma Cells Expressing Isocitrate Dehydrogenase 1 (IDH1) Mutant Protein: Dissecting the Mechanisms

The presence of an isocitrate dehydrogenase 1 (IDH1) mutation is associated with a less aggressive phenotype, increased sensitivity to radiation, and increased overall survival in patients with diffuse glioma. Based on in vitro experimentations in malignant glioma cell lines, the consequences on cellular processes of IDH1^R132H expression were analyzed. The results revealed that IDH1^R132H expression enhanced the radiation induced accumulation of residual γH2AX foci and decreased the amount of glutathione (GSH) independent of the oxygen status. In addition, expression of the mutant IDH1 caused a significant increase of cell stiffness and induced an altered organization of the cytoskeleton, which has been shown to reinforce cell stiffness. Furthermore, IDH1^R132H expression decreased the expression of vimentin, an important component of the cytoskeleton and regulator of the cell stiffness. The results emphasize the important role of mutant IDH1 in treatment of patients with diffuse gliomas especially in response to radiation. Hence, detection of the genetic status of IDH1 before therapy massively expands the utility of immunohistochemistry to accurately distinguish patients with a less aggressive and radiosensitive IDH1-mutant diffuse glioma suitable for radiotherapy from those with a more aggressive IDH1-wildtype diffuse glioma who might benefit from an individually intensified therapy comprising radiotherapy and alternative medical treatments.