RADAR: a rigorously annotated database of A-to-I RNA editing

Combined analysis of somatic mutations and gene expression reveals nuclear speckles-associated enhanced stemness in gingivobuccal carcinoma under DNA damage response

Smokeless tobacco chewing habits in India lead to a high prevalence of Gingivobuccal oral squamous cell carcinoma (OSCC-GB). Cancer stem cells (CSCs) are a sub-population of cancer cells within a tumor with stem-like properties and are believed to contribute to tumor initiation, progression, increased resistance to drug therapy, and promote post-therapeutic cancer relapse. An RNA-sequencing data-based combined analysis of somatic mutations and gene expression was performed to explore the role of CSCs in disease progression using the novel Indian-origin OSCC-GB cell line 'IIOC019' from a patient with tobacco-chewing habit. The identified DNA damage-related known mutational signature (1 bp T/(A) nucleotide insertions and C>T mutations) indicates the impact of smokeless tobacco-related carcinogens in the IIOC019 cell line. The differentially expressed somatic variants, functional impact predictions, and survival analysis reveal the role of DNA damage response (DDR)-related genes in OSCC-GB, with the SON gene as a significant player. The study suggests that the loss-of-function in a somatic variant of the SON gene is linked to nuclear speckles-associated enhanced stemness and increased risk of disease progression in OSCC-GB under DDR conditions. The newly identified CSC-associated somatic variants appear to promote cancer spread, local recurrence, and resistance to chemotherapy or radiotherapy, contributing to the high mortality rates among Indian OSCC-GB patients.

Advances in RNA editing in hematopoiesis and associated malignancies

ABSTRACT

Adenosine-to-inosine (A-to-I) RNA editing is a prevalent RNA modification essential for cell survival. The process is catalyzed by the adenosine deaminase acting on RNA (ADAR) enzyme family that converts adenosines in double-stranded RNAs (dsRNAs) into inosines, which are read as guanosines during translation. Deep sequencing has helped to reveal that A-to-I editing occurs across various types of RNAs, affecting their functions. RNA editing detection is now so sophisticated that we can achieve a high level of accuracy and sensitivity to identify low-abundance edited events. Consequently, A-to-I editing has been implicated in various biological processes, including immune and stress responses, cancer progression, and stem cell fate determination. In particular, a crucial role for this process has been recently reported in hematopoietic cell development and hematologic malignancy progression. Results from genetic mouse models have demonstrated the impact of ADARs' catalytic activity on hematopoietic cells, complemented by insights from human cell studies. Meanwhile, clinical studies have implicated ADAR enzymes and RNA editing events in hematologic malignancies and highlighted their potential as prognostic indicators. In this review, we outline the regulatory mechanisms of RNA editing in both normal hematopoiesis and hematologic malignancies. We then speculate on how targeting ADAR expression and site-specific RNA substrates might serve as a therapeutic avenue for affected patients.

Phylogenetic and structural analysis of Hydra ADAR

Adenosine deaminases acting on RNAs (ADARs) perform adenosine-to-inosine (A-to-I) RNA editing for essential biological functions. While studies of editing sites in diverse animals have revealed unique biological roles of ADAR editing including temperature adaptation and reproductive maturation, rigorous biochemical and structural studies of these ADARs are lacking. Here, we present a phylogenetic sequence analysis and AlphaFold computational structure prediction to reveal that medusozoan ADAR2s contain five dsRNA binding domains (dsRBDs) with several RNA binding residues in the dsRBDs and deaminase domain conserved. Additionally, we identified evolutionary divergence between the medusozoan (e.g. Hydra) and anthozoan cnidarian subphyla. The anthozoan ADAR deaminase domains more closely resemble human ADARs with longer 5' RNA binding loops, glutamate base-flipping residues, and a conserved TWDG dimerization motif. Conversely, medusozoan ADAR deaminase domains have short 5' binding loops, glutamine flipping residues, and non-conserved helix dimerization motif. We also report the direct detection of A-to-I RNA editing by an ADAR ortholog from the freshwater cnidarian Hydra vulgaris (hyADAR). We solved the crystal structure of the monomeric deaminase domain of hyADAR (hyADARd) to 2.0 Å resolution, showing conserved active site architecture and the presence of a buried inositol hexakisphosphate known to be required for ADAR activity. In addition, these data demonstrate that medusozoans have evolved novel ADAR structural features, however the physiological consequence of this remains unknown. In addition, these results provide a framework for biochemically and structurally characterizing ADARs from evolutionarily distant organisms to understand the diverse roles of ADAR editing amongst metazoans.

Unveiling RNA Editing by ADAR and APOBEC Protein Gene Families

RNA editing is a crucial post-transcriptional modification that alters the transcriptome and proteome and affects many cellular processes, including splicing, microRNA specificity, stability of RNA molecules, and protein structure. Enzymes from the adenosine deaminase acting on RNA (ADAR) and apolipoprotein B mRNA editing catalytic polypeptide-like (APOBEC) protein families mediate RNA editing and can alter a variety of non-coding and coding RNAs, including all regions of mRNA molecules, leading to tumor development and progression. This review provides novel insights into the potential use of RNA editing parameters, such as editing levels, expression of ADAR and APOBEC genes, and specifically edited genes, as biomarkers for cancer progression, distinguishing it from previous studies that focused on isolated aspects of RNA editing mechanisms. The methodological section offers clues to accelerate high-throughput analysis of RNA or DNA sequencing data for the identification of RNA editing events.

Improved RNA base editing with guide RNAs mimicking highly edited endogenous ADAR substrates

Adenosine deaminase acting on RNA (ADAR)-mediated RNA base editing offers a safer alternative to genome editing for specific clinical applications because of nonpermanent editing of targets. Current guide RNA (gRNA) designs feature a fully complementary specificity domain with an A-C mismatch at the targeted adenosine. However, perfectly matched dsRNA is not the most effective ADAR substrate. Here we introduce MIRROR (mimicking inverted repeats to recruit ADARs using engineered oligoribonucleotides), an approach that implements structural motifs derived from highly edited inverted Alu repeats in human tissues to enable rational gRNA design for ADAR recruitment. We demonstrated that MIRROR is applicable to both short chemically synthesized gRNAs with modifications and long biologically generated gRNAs and surpasses current state-of-the-art approaches in both gRNA forms. It enhances editing efficiency by up to 5.7-fold in multiple human cell types and primary hepatocytes from an alpha-1 antitrypsin deficiency mouse model. Our findings improve programmable RNA editing in vitro and in vivo by rational design through the screening of highly edited natural substrate mimics.

Biochemical profiling and structural basis of ADAR1-mediated RNA editing

ADAR1 regulates RNA-induced immune responses by converting adenosine to inosine in double-stranded RNA. Mutations in ADAR1 are associated with human autoimmune disease, and targeting ADAR1 has been proposed for cancer immunotherapy. However, the molecular mechanisms underlying ADAR1-mediated editing remain unclear. Here, we provide detailed biochemical and structural characterizations of human ADAR1. Our biochemical profiling reveals that ADAR1 editing is both sequence and RNA-duplex-length dependent but can well tolerate mismatches near the editing site. High-resolution ADAR1-RNA complex structures, combined with mutagenesis, elucidate RNA binding, substrate selection, dimerization, and the essential role of RNA-binding domain 3. The ADAR1 structures also help explain the potential defects of disease-associated mutations, where biochemical and RNA sequencing analysis further indicate some of the mutations preferentially impact the editing of RNAs with short duplexes. These findings unveil the molecular basis of ADAR1 editing and provide insights into its immune-regulatory functions and therapeutic potential.

An Ultimate Question for Functional A-to-I mRNA Editing: Why Not a Genomic G?

A-to-I mRNA editing resembles A-to-G mutations. Functional mRNA editing, representing only a corner of total editing events, can be inferred from the experimental removal of editing. However, it is intuitive to ask why evolution chose RNA editing rather than directly (and simply) changing the genomic sequence to G? If G is better than A, then drift or constructive neutral evolution (CNE) theory can explain the emergence of such editing, but it is still unclear why the exemplified conserved editing is perfectly maintained without observing any subsequent A-to-G DNA mutations? Virtually every functional and conserved mRNA editing site faces this ultimate question until one justifies that being editable is better than a hardwired genomic allele. While the advantage of editability has been validated in fungi, this ultimate question has not been answered for any functional editing sites in animals. By providing several conceptual arguments and specific examples, we propose that proving the evolutionary adaptiveness of an editing site is far more difficult than revealing its function.

REDInet: a temporal convolutional network-based classifier for A-to-I RNA editing detection harnessing million known events

A-to-I ribonucleic acid (RNA) editing detection is still a challenging task. Current bioinformatics tools rely on empirical filters and whole genome sequencing or whole exome sequencing data to remove background noise, sequencing errors, and artifacts. Sometimes they make use of cumbersome and time-consuming computational procedures. Here, we present REDInet, a temporal convolutional network-based deep learning algorithm, to profile RNA editing in human RNA sequencing (RNAseq) data. It has been trained on REDIportal RNA editing sites, the largest collection of human A-to-I changes from >8000 RNAseq data of the genotype-tissue expression project. REDInet can classify editing events with high accuracy harnessing RNAseq nucleotide frequencies of 101-base windows without the need for coupled genomic data.

Editome Disease Knowledgebase v2.0: an updated resource of editome-disease associations through literature curation and integrative analysis

Motivation

Editome Disease Knowledgebase (EDK) is a curated resource of knowledge between RNA editome and human diseases. Since its first release in 2018, a number of studies have discovered previously uncharacterized editome-disease associations and generated an abundance of RNA editing datasets. Thus, it is desirable to make significant updates for EDK by incorporating more editome-disease associations as well as their related editing profiles.

Results

Here, we present EDK v2.0, an updated version of editome-disease associations based on both literature curation and integrative analysis. EDK v2.0 incorporates a curated collection of 1097 editome-disease associations involving 115 diseases from 321 publications. Meanwhile, based on a standardized pipeline, EDK v2.0 provides RNA editing profiles from 48 datasets covering 2536 samples across 55 diseases. Through differential analysis on RNA editing, it further identifies a total of 7190 differential edited genes and 86 242 differential editing sites (DESs), leading to 266 339 DES-disease associations. Moreover, a curated list of 28 160 cis-RNA editing QTL associations, 458 187 DES-RNA binding protein associations, and 21 DES-RNA secondary structure associations are annotated and added to EDK v2.0. Additionally, it is equipped with a series of user-friendly tools to facilitate RNA editing online analysis.

Availability and implementation

https://ngdc.cncb.ac.cn/edk/.

M6Allele: a toolkit for detection of allele-specific RNA N6-methyladenosine modifications

BACKGROUND

Allelic gene-specific regulatory events are crucial mechanisms in organisms, pivotal to many fundamental biological processes such as embryonic development and chromosome inactivation. Allelic gene imbalance manifests at both RNA expression and epigenetic levels. Recent research has unveiled allelic-specific regulation of RNA N6-methyladenosine (m6A), emphasizing the need for its precise identification. However, prevailing approaches primarily focus on screening allele-specific genetic variations associated with m6A, but not truly identify allelic m6A events. Therefore, the construction of a novel algorithm dedicated to identifying allele-specific m6A (ASm6A) signals is still necessary for comprehensively understanding the regulatory mechanism of ASm6A.

FINDINGS

To address this limitation, we have developed a meta-analysis approach using hierarchical Bayesian models to accurately detect ASm6A events at the peak level from MeRIP-seq data. For user convenience, we introduce a unified analysis pipeline named M6Allele, streamlining the assessment of significant ASm6A across single and paired samples. Applying M6Allele to MeRIP-seq data analysis of pulmonary fibrosis and lung adenocarcinoma reveals enrichment of ASm6A events in key regulatory genes associated with these diseases, suggesting their potential involvement in disease regulation.

CONCLUSIONS

Our effort provides a method for precisely identifying ASm6A events at the peak level, elucidates the interplay of m6A with human health and disease genetics, and paves a new visual angle for disease research. The M6Allele software is freely available at https://github.com/RenLabBioinformatics/M6Allele under the MIT license.

REDIportal: toward an integrated view of the A-to-I editing

A-to-I RNA editing is the most common non-transient epitranscriptome modification. It plays several roles in human physiology and has been linked to several disorders. Large-scale deep transcriptome sequencing has fostered the characterization of A-to-I editing at the single nucleotide level and the development of dedicated computational resources. REDIportal is a unique and specialized database collecting ∼16 million of putative A-to-I editing sites designed to face the current challenges of epitranscriptomics. Its running version has been enriched with sites from the TCGA project (using data from 31 studies). REDIportal provides an accurate, sustainable and accessible tool enriched with interconnections with widespread ELIXIR core resources such as Ensembl, RNAcentral, UniProt and PRIDE. Additionally, REDIportal now includes information regarding RNA editing in putative double-stranded RNAs, relevant for the immune-related roles of editing, as well as an extended catalog of recoding events. Finally, we report a reliability score per site calculated using a deep learning model trained using a huge collection of positive and negative instances. REDIportal is available at http://srv00.recas.ba.infn.it/atlas/.

ncRNA Editing: Functional Characterization and Computational Resources

Non-coding RNAs (ncRNAs) play crucial roles in gene expression regulation, translation, and disease development, including cancer. They are classified by size in short and long non-coding RNAs. This chapter focuses on the functional implications of adenosine-to-inosine (A-to-I) RNA editing in both short (e.g., miRNAs) and long ncRNAs. RNA editing dynamically alters the sequence and structure of primary transcripts, impacting ncRNA biogenesis and function. Notable findings include the role of miRNA editing in promoting glioblastoma invasiveness, characterizing RNA editing hotspots across cancers, and its implications in thyroid cancer and ischemia. This chapter also highlights bioinformatics resources and next-generation sequencing (NGS) technologies that enable comprehensive ncRNAome studies and genome-wide RNA editing detection. Dysregulation of RNA editing machinery has been linked to various human diseases, emphasizing the potential of RNA editing as a biomarker and therapeutic target. This overview integrates current knowledge and computational tools for studying ncRNA editing, providing insights into its biological significance and clinical applications.

Nanopore sequencing to detect A-to-I editing sites

Adenosine-to-inosine (A-to-I) RNA editing, mediated by the ADAR family of enzymes, is pervasive in metazoans and functions as an important mechanism to diversify the proteome and control gene expression. Over the years, there have been multiple efforts to comprehensively map the editing landscape in different organisms and in different disease states. As inosine (I) is recognized largely as guanosine (G) by cellular machineries including the reverse transcriptase, editing sites can be detected as A-to-G changes during sequencing of complementary DNA (cDNA). However, such an approach is indirect and can be confounded by genomic single nucleotide polymorphisms (SNPs) and DNA mutations. Moreover, past studies rely primarily on the Illumina platform, which generates short sequencing reads that can be challenging to map. Recently, nanopore direct RNA sequencing has emerged as a powerful technology to address the issues. Here, we describe the use of the technology together with deep learning models that we have developed, named Dinopore (Detection of inosine with nanopore sequencing), to interrogate the A-to-I editome of any organism.

Narrowing down the candidates of beneficial A-to-I RNA editing by comparing the recoding sites with uneditable counterparts

Adar-mediated adenosine-to-inosine (A-to-I) RNA editing mainly occurs in nucleus and diversifies the transcriptome in a flexible manner. It has been a challenging task to identify beneficial editing sites from the sea of total editing events. The functional Ser>Gly auto-recoding site in insect Adar gene has uneditable Ser codons in ancestral nodes, indicating the selective advantage to having an editable status. Here, we extended this case study to more metazoan species, and also looked for all Drosophila recoding events with potential uneditable synonymous codons. Interestingly, in D. melanogaster, the abundant nonsynonymous editing is enriched in the codons that have uneditable counterparts, but the Adar Ser>Gly case suggests that the editable orthologous codons in other species are not necessarily edited. The use of editable versus ancestral uneditable codon is a smart way to infer the selective advantage of RNA editing, and priority might be given to these editing sites for functional studies due to the feasibility to construct an uneditable allele. Our study proposes an idea to narrow down the candidates of beneficial recoding sites. Meanwhile, we stress that the matched transcriptomes are needed to verify the conservation of editing events during evolution.

Conserved A-to-I RNA editing with non-conserved recoding expands the candidates of functional editing sites

Adenosine-to-inosine (A-to-I) RNA editing recodes the genome and confers flexibility for the organisms to adapt to the environment. It is believed that RNA recoding sites are well suited for facilitating adaptive evolution by increasing the proteomic diversity in a temporal-spatial manner. The function and essentiality of a few conserved recoding sites are recognized. However, the experimentally discovered functional sites only make up a small corner of the total sites, and there is still the need to expand the repertoire of such functional sites with bioinformatic approaches. In this study, we define a new category of RNA editing sites termed 'conserved editing with non-conserved recoding' and systematically identify such sites in Drosophila editomes, figuring out their selection pressure and signals of adaptation at inter-species and intra-species levels. Surprisingly, conserved editing sites with non-conserved recoding are not suppressed and are even slightly overrepresented in Drosophila. DNA mutations leading to such cases are also favoured during evolution, suggesting that the function of those recoding events in different species might be diverged, specialized, and maintained. Finally, structural prediction suggests that such recoding in potassium channel Shab might increase ion permeability and compensate the effect of low temperature. In conclusion, conserved editing with non-conserved recoding might be functional as well. Our study provides novel aspects in considering the adaptive evolution of RNA editing sites and meanwhile expands the candidates of functional recoding sites for future validation.

Bioinformatic approaches for accurate assessment of A-to-I editing in complete transcriptomes

A-to-I RNA editing is an RNA modification that alters the RNA sequence relative to the its genomic blueprint. It is catalyzed by double-stranded RNA-specific adenosine deaminase (ADAR) enzymes, and contributes to the complexity and diversification of the proteome. Advancement in the study of A-to-I RNA editing has been facilitated by computational approaches for accurate mapping and quantification of A-to-I RNA editing based on sequencing data. In this chapter we review some of the main computational approaches currently used, describe potential hurdles, challenges and pitfalls, and discuss possible ways to mitigate them.

GGNBP2 regulates MDA5 sensing triggered by self double-stranded RNA following loss of ADAR1 editing

Adenosine-to-inosine (A-to-I) editing of double-stranded RNA (dsRNA) by ADAR1 is an essential modifier of the immunogenicity of cellular dsRNA. The role of MDA5 in sensing unedited cellular dsRNA and the downstream activation of type I interferon (IFN) signaling are well established. However, we have an incomplete understanding of pathways that modify the response to unedited dsRNA. We performed a genome-wide CRISPR screen and showed that GGNBP2, CNOT10, and CNOT11 interact and regulate sensing of unedited cellular dsRNA. We found that GGNBP2 acts between dsRNA transcription and its cytoplasmic sensing by MDA5. GGNBP2 loss prevented induction of type I IFN and autoinflammation after the loss of ADAR1 editing activity by modifying the subcellular distribution of endogenous A-to-I editing substrates and reducing cytoplasmic dsRNA load. These findings reveal previously undescribed pathways to modify diseases associated with ADAR mutations and may be determinants of response or resistance to small-molecule ADAR1 inhibitors.

Small and long non-coding RNAs: Past, present, and future

Since the introduction of the central dogma of molecular biology in 1958, various RNA species have been discovered. Messenger RNAs transmit genetic instructions from DNA to make proteins, a process facilitated by housekeeping non-coding RNAs (ncRNAs) such as small nuclear RNAs (snRNAs), ribosomal RNAs (rRNAs), and transfer RNAs (tRNAs). Over the past four decades, a wide array of regulatory ncRNAs have emerged as crucial players in gene regulation. In celebration of Cell's 50th anniversary, this Review explores our current understanding of the most extensively studied regulatory ncRNAs-small RNAs and long non-coding RNAs (lncRNAs)-which have profoundly shaped the field of RNA biology and beyond. While small RNA pathways have been well documented with clearly defined mechanisms, lncRNAs exhibit a greater diversity of mechanisms, many of which remain unknown. This Review covers pivotal events in their discovery, biogenesis pathways, evolutionary traits, action mechanisms, functions, and crosstalks among ncRNAs. We also highlight their roles in pathophysiological contexts and propose future research directions to decipher the unknowns of lncRNAs by leveraging lessons from small RNAs.

Smooth muscle expression of RNA editing enzyme ADAR1 controls activation of RNA sensor MDA5 in atherosclerosis

Mapping the genomic architecture of complex disease has been predicated on the understanding that genetic variants influence disease risk through modifying gene expression. However, recent discoveries have revealed that a significant burden of disease heritability in common autoinflammatory disorders and coronary artery disease is mediated through genetic variation modifying post-transcriptional modification of RNA through adenosine-to-inosine (A-to-I) RNA editing. This common RNA modification is catalyzed by ADAR enzymes, where ADAR1 edits specific immunogenic double stranded RNA (dsRNA) to prevent activation of the double strand RNA (dsRNA) sensor MDA5 ( IFIH1 ) and stimulation of an interferon stimulated gene (ISG) response. Multiple lines of human genetic data indicate impaired RNA editing and increased dsRNA sensing by MDA5 to be an important mechanism of coronary artery disease (CAD) risk. Here, we provide a crucial link between observations in human genetics and mechanistic cell biology leading to progression of CAD. Through analysis of human atherosclerotic plaque, we implicate the vascular smooth muscle cell (SMC) to have a unique requirement for RNA editing, and that ISG induction occurs in SMC phenotypic modulation, implicating MDA5 activation. Through culture of human coronary artery SMCs, generation of a conditional SMC specific Adar1 deletion mouse model on a pro-atherosclerosis background with additional constitutive deletion of MDA5 ( Ifih1 ), and with incorporation of single cell RNA sequencing cellular profiling, we further show that Adar1 controls SMC phenotypic state by regulating Mda5 activation, is required to maintain vascular integrity, and controls progression of atherosclerosis and vascular calcification. Through this work, we describe a fundamental mechanism of CAD, where cell type and context specific RNA editing and sensing of dsRNA mediates disease progression, bridging our understanding of human genetics and disease causality.

One Sentence Summary

Smooth muscle expression of RNA editing enzyme ADAR1 regulates activation of double strand RNA sensor MDA5 in novel mechanism of atherosclerosis.

LoDEI: a robust and sensitive tool to detect transcriptome-wide differential A-to-I editing in RNA-seq data

RNA editing is a highly conserved process. Adenosine deaminase acting on RNA (ADAR) mediated deamination of adenosine (A-to-I editing) is associated with human disease and immune checkpoint control. Functional implications of A-to-I editing are currently of broad interest to academic and industrial research as underscored by the fast-growing number of clinical studies applying base editors as therapeutic tools. Analyzing the dynamics of A-to-I editing, in a biological or therapeutic context, requires the sensitive detection of differential A-to-I editing, a currently unmet need. We introduce the local differential editing index (LoDEI) to detect differential A-to-I editing in RNA-seq datasets using a sliding-window approach coupled with an empirical q value calculation that detects more A-to-I editing sites at the same false-discovery rate compared to existing methods. LoDEI is validated on known and novel datasets revealing that the oncogene MYCN increases and that a specific small non-coding RNA reduces A-to-I editing.

Global A-to-I RNA editing during myogenic differentiation of goat MuSCs

Background

RNA editing, especially A-to-I editing sites, is a common RNA modification critical for stem cell differentiation, muscle development, and disease occurrence. Unveiling comprehensive RNA A-to-I editing events associated with myogenesis of the skeletal muscle satellite cells (MuSCs) is essential for extending our knowledge of the mechanism underpinning muscle development.

Results

A total of 9,632 RNA editing sites (RESs) were screened in the myoblasts (GM), myocytes (DM1), and myotubes (DM5) samples. Among these sites, 4,559 A-to-I edits were classified and further analyzed. There were 3,266 A-to-I sites in the protein-coding region, out of which 113 missense sites recoded protein. Notably, five A-to-I sites in the 3' UTR of four genes (TRAF6, NALF1, SLC38A1, ENSCHIG00000019092) altered their targeted miRNAs. Furthermore, a total of 370 A-to-I sites with different editing levels were detected, including FBN1, MYH10, GSK3B, CSNK1D, and PRKACB genes. These genes were predominantly enriched in the cytoskeleton in muscle cells, the hippo signaling pathway, and the tight junction. Furthermore, we identified 14 hub genes (TUFM, GSK3B, JAK2, RPSA, YARS1, CDH2, PRKACB, RUNX1, NOTCH2, CDC23, VCP, FBN1, RARS1, MEF2C) that potentially related to muscle development. Additionally, 123 stage-specific A-to-I editing sites were identified, with 43 sites in GM, 25 in DM1, and 55 in DM5 samples. These stage-specific edited genes significantly enriched essential biological pathways, including the cell cycle, oocyte meiosis, motor proteins, and hedgehog signaling pathway.

Conclusion

We systematically identified the RNA editing events in proliferating and differentiating goat MuSCs, which was crucial for expanding our understanding of the regulatory mechanisms of muscle development.

DEMINING: A deep learning model embedded framework to distinguish RNA editing from DNA mutations in RNA sequencing data

Precise calling of promiscuous adenosine-to-inosine RNA editing sites from transcriptomic datasets is hindered by DNA mutations and sequencing/mapping errors. Here, we present a stepwise computational framework, called DEMINING, to distinguish RNA editing and DNA mutations directly from RNA sequencing datasets, with an embedded deep learning model named DeepDDR. After transfer learning, DEMINING can also classify RNA editing sites and DNA mutations from non-primate sequencing samples. When applied in samples from acute myeloid leukemia patients, DEMINING uncovers previously underappreciated DNA mutation and RNA editing sites; some associated with the upregulated expression of host genes or the production of neoantigens.

A comprehensive atlas of pig RNA editome across 23 tissues reveals RNA editing affecting interaction mRNA-miRNAs

RNA editing is a co-transcriptional/post-transcriptional modification that is mediated by the ADAR enzyme family. Profiling of RNA editing is very limited in pigs. In this study, we collated 3813 RNA-seq data from the public repositories across 23 tissues and carried out comprehensive profiling of RNA editing in pigs. In total, 127,927 A-to-I RNA-editing sites were detected. Our analysis showed that 98.2% of RNA-editing sites were located within repeat regions, primarily within the pig-specific SINE retrotransposon PRE-1/Pre0_SS elements. Subsequently, we focused on analyzing specific RNA-editing sites (SESs) in skeletal muscle tissues. Functional enrichment analyses suggested that they were enriched in signaling pathways associated with muscle cell differentiation, including DMD, MYOD1, and CAV1 genes. Furthermore, we discovered that RNA editing event in the 3'UTR of CFLAR mRNA influenced miR-708-5p binding in this region. In this study, the panoramic RNA-editing landscape of different tissues of pigs was systematically mapped, and RNA-editing sites and genes involved in muscle cell differentiation were identified. In summary, we identified modifications to pig RNA-editing sites and provided candidate targets for further validation.

An isoform-resolution transcriptomic atlas of colorectal cancer from long-read single-cell sequencing

Colorectal cancer (CRC) ranks as the second leading cause of cancer deaths globally. In recent years, short-read single-cell RNA sequencing (scRNA-seq) has been instrumental in deciphering tumor heterogeneities. However, these studies only enable gene-level quantification but neglect alterations in transcript structures arising from alternative end processing or splicing. In this study, we integrated short- and long-read scRNA-seq of CRC samples to build an isoform-resolution CRC transcriptomic atlas. We identified 394 dysregulated transcript structures in tumor epithelial cells, including 299 resulting from various combinations of splicing events. Second, we characterized genes and isoforms associated with epithelial lineages and subpopulations exhibiting distinct prognoses. Among 31,935 isoforms with novel junctions, 330 were supported by The Cancer Genome Atlas RNA-seq and mass spectrometry data. Finally, we built an algorithm that integrated novel peptides derived from open reading frames of recurrent tumor-specific transcripts with mass spectrometry data and identified recurring neoepitopes that may aid the development of cancer vaccines.

Loss of ADAR1 protein induces changes in small RNA landscape in hepatocytes

In recent years, numerous evidence has been accumulated about the extent of A-to-I editing in human RNAs and the key role ADAR1 plays in the cellular editing machinery. It has been shown that A-to-I editing occurrence and frequency are tissue-specific and essential for some tissue development, such as the liver. To study the effect of ADAR1 function in hepatocytes, we have created Huh7.5 ADAR1 KO cell lines. Upon IFN treatment, the Huh7.5 ADAR1 KO cells show rapid arrest of growth and translation, from which they do not recover. We analyzed translatome changes by using a method based on sequencing of separate polysome profile RNA fractions. We found significant changes in the transcriptome and translatome of the Huh7.5 ADAR1 KO cells. The most prominent changes include negatively affected transcription by RNA polymerase III and the deregulation of snoRNA and Y RNA levels. Furthermore, we observed that ADAR1 KO polysomes are enriched in mRNAs coding for proteins pivotal in a wide range of biological processes such as RNA localization and RNA processing, whereas the unbound fraction is enriched mainly in mRNAs coding for ribosomal proteins and translational factors. This indicates that ADAR1 plays a more relevant role in small RNA metabolism and ribosome biogenesis.

Adenosine-to-Inosine RNA editing in cancer: molecular mechanisms and downstream targets

Adenosine-to-Inosine (A-to-I), one of the most prevalent RNA modifications, has recently garnered significant attention. The A-to-I modification actively contributes to biological and pathological processes by affecting the structure and function of various RNA molecules, including double stranded RNA, transfer RNA, microRNA, and viral RNA. Increasing evidence suggests that A-to-I plays a crucial role in the development of human disease, particularly in cancer, and aberrant A-to-I levels are closely associated with tumorigenesis and progression through regulation of the expression of multiple oncogenes and tumor suppressor genes. Currently, the underlying molecular mechanisms of A-to-I modification in cancer are not comprehensively understood. Here, we review the latest advances regarding the A-to-I editing pathways implicated in cancer, describing their biological functions and their connections to the disease.

RBPs: an RNA editor's choice

RNA-binding proteins (RBPs) play a key role in gene expression and post-transcriptional RNA regulation. As integral components of ribonucleoprotein complexes, RBPs are susceptible to genomic and RNA Editing derived amino acid substitutions, impacting functional interactions. This article explores the prevalent RNA Editing of RBPs, unravelling the complex interplay between RBPs and RNA Editing events. Emphasis is placed on their influence on single amino acid variants (SAAVs) and implications for disease development. The role of Proteogenomics in identifying SAAVs is briefly discussed, offering insights into the RBP landscape. RNA Editing within RBPs emerges as a promising target for precision medicine, reshaping our understanding of genetic and epigenetic variations in health and disease.

Bioinformatics for Inosine: Tools and Approaches to Trace This Elusive RNA Modification

Inosine is a nucleotide resulting from the deamination of adenosine in RNA. This chemical modification process, known as RNA editing, is typically mediated by a family of double-stranded RNA binding proteins named Adenosine Deaminase Acting on dsRNA (ADAR). While the presence of ADAR orthologs has been traced throughout the evolution of metazoans, the existence and extension of RNA editing have been characterized in a more limited number of animals so far. Undoubtedly, ADAR-mediated RNA editing plays a vital role in physiology, organismal development and disease, making the understanding of the evolutionary conservation of this phenomenon pivotal to a deep characterization of relevant biological processes. However, the lack of direct high-throughput methods to reveal RNA modifications at single nucleotide resolution limited an extended investigation of RNA editing. Nowadays, these methods have been developed, and appropriate bioinformatic pipelines are required to fully exploit this data, which can complement existing approaches to detect ADAR editing. Here, we review the current literature on the "bioinformatics for inosine" subject and we discuss future research avenues in the field.

Learning from the Codon Table: Convergent Recoding Provides Novel Understanding on the Evolution of A-to-I RNA Editing

Adenosine-to-inosine (A-to-I) RNA editing recodes the genetic information. Apart from diversifying the proteome, another tempting advantage of RNA recoding is to correct deleterious DNA mutation and restore ancestral allele. Solid evidences for beneficial restorative editing are very rare in animals. By searching for "convergent recoding" under a phylogenetic context, we proposed this term for judging the potential restorative functions of particular editing site. For the well-known mammalian Gln>Arg (Q>R) recoding site, its ancestral state in vertebrate genomes was the pre-editing Gln, and all 470 available mammalian genomes strictly avoid other three equivalent ways to achieve Arg in protein. The absence of convergent recoding from His>Arg, or synonymous mutations on Gln codons, could be attributed to the strong maintenance on editing motif and structure, but the absence of direct A-to-G mutation is extremely unexpected. With similar ideas, we found cases of convergent recoding in Drosophila genus, reducing the possibility of their restorative function. In summary, we defined an interesting scenario of convergent recoding, the occurrence of which could be used as preliminary judgements for whether a recoding site has a sole restorative role. Our work provides novel insights to the natural selection and evolution of RNA editing.

A-to-I RNA co-editing predicts clinical outcomes and is associated with immune cells infiltration in hepatocellular carcinoma

Aberrant RNA editing has emerged as a pivotal factor in the pathogenesis of hepatocellular carcinoma (HCC), but the impact of RNA co-editing within HCC remains underexplored. We used a multi-step algorithm to construct an RNA co-editing network in HCC, and found that HCC-related RNA editings are predominantly centralized within the network. Furthermore, five pairs of risk RNA co-editing events were significantly correlated with the overall survival in HCC. Based on presence of risk RNA co-editings resulted in the categorization of HCC patients into high-risk and low-risk groups. Disparities in immune cell infiltrations were observed between the two groups, with the high-risk group exhibiting a greater abundance of exhausted T cells. Additionally, seven genes associated with risk RNA co-editing pairs were identified, whose expression effectively differentiates HCC tumor samples from normal ones. Our research offers an innovative perspective on the etiology and potential therapeutics for HCC.

ADAR-Mediated A>I(G) RNA Editing in the Genotoxic Drug Response of Breast Cancer

Epitranscriptomics is a field that delves into post-transcriptional changes. Among these modifications, the conversion of adenosine to inosine, traduced as guanosine (A>I(G)), is one of the known RNA-editing mechanisms, catalyzed by ADARs. This type of RNA editing is the most common type of editing in mammals and contributes to biological diversity. Disruption in the A>I(G) RNA-editing balance has been linked to diseases, including several types of cancer. Drug resistance in patients with cancer represents a significant public health concern, contributing to increased mortality rates resulting from therapy non-responsiveness and disease progression, representing the greatest challenge for researchers in this field. The A>I(G) RNA editing is involved in several mechanisms over the immunotherapy and genotoxic drug response and drug resistance. This review investigates the relationship between ADAR1 and specific A>I(G) RNA-edited sites, focusing particularly on breast cancer, and the impact of these sites on DNA damage repair and the immune response over anti-cancer therapy. We address the underlying mechanisms, bioinformatics, and in vitro strategies for the identification and validation of A>I(G) RNA-edited sites. We gathered databases related to A>I(G) RNA editing and cancer and discussed the potential clinical and research implications of understanding A>I(G) RNA-editing patterns. Understanding the intricate role of ADAR1-mediated A>I(G) RNA editing in breast cancer holds significant promise for the development of personalized treatment approaches tailored to individual patients' A>I(G) RNA-editing profiles.

Detecting haplotype-specific transcript variation in long reads with FLAIR2

BACKGROUND

RNA-seq has brought forth significant discoveries regarding aberrations in RNA processing, implicating these RNA variants in a variety of diseases. Aberrant splicing and single nucleotide variants (SNVs) in RNA have been demonstrated to alter transcript stability, localization, and function. In particular, the upregulation of ADAR, an enzyme that mediates adenosine-to-inosine editing, has been previously linked to an increase in the invasiveness of lung adenocarcinoma cells and associated with splicing regulation. Despite the functional importance of studying splicing and SNVs, the use of short-read RNA-seq has limited the community's ability to interrogate both forms of RNA variation simultaneously.

RESULTS

We employ long-read sequencing technology to obtain full-length transcript sequences, elucidating cis-effects of variants on splicing changes at a single molecule level. We develop a computational workflow that augments FLAIR, a tool that calls isoform models expressed in long-read data, to integrate RNA variant calls with the associated isoforms that bear them. We generate nanopore data with high sequence accuracy from H1975 lung adenocarcinoma cells with and without knockdown of ADAR. We apply our workflow to identify key inosine isoform associations to help clarify the prominence of ADAR in tumorigenesis.

CONCLUSIONS

Ultimately, we find that a long-read approach provides valuable insight toward characterizing the relationship between RNA variants and splicing patterns.

A-to-I Editing Is Subtype-Specific in Non-Hodgkin Lymphomas

Cancer is a complex and heterogeneous disease, in which a number of genetic and epigenetic changes occur in tumor onset and progression. Recent studies indicate that changes at the RNA level are also involved in tumorigenesis, such as adenosine-to-inosine (A-to-I) RNA editing. Here, we systematically investigate transcriptome-wide A-to-I editing events in a large number of samples from Non-Hodgkin lymphomas (NHLs). Using a computational pipeline that determines significant differences in editing level between NHL and normal samples at known A-to-I editing sites, we identify a number of differentially edited editing sites between NHL subtypes and normal samples. Most of the differentially edited sites are located in non-coding regions, and many such sites show a strong correlation between gene expression level and editing efficiency, indicating that RNA editing might have direct consequences for the cancer cell's aberrant gene regulation status in these cases. Moreover, we establish a strong link between RNA editing and NHL by demonstrating that NHL and normal samples and even NHL subtypes can be distinguished based on genome-wide RNA editing profiles alone. Our study establishes a strong link between RNA editing, cancer and aberrant gene regulation in NHL.

The implications of APOBEC3-mediated C-to-U RNA editing for human disease

Intra-organism biodiversity is thought to arise from epigenetic modification of constituent genes and post-translational modifications of translated proteins. Here, we show that post-transcriptional modifications, like RNA editing, may also contribute. RNA editing enzymes APOBEC3A and APOBEC3G catalyze the deamination of cytosine to uracil. RNAsee (RNA site editing evaluation) is a computational tool developed to predict the cytosines edited by these enzymes. We find that 4.5% of non-synonymous DNA single nucleotide polymorphisms that result in cytosine to uracil changes in RNA are probable sites for APOBEC3A/G RNA editing; the variant proteins created by such polymorphisms may also result from transient RNA editing. These polymorphisms are associated with over 20% of Medical Subject Headings across ten categories of disease, including nutritional and metabolic, neoplastic, cardiovascular, and nervous system diseases. Because RNA editing is transient and not organism-wide, future work is necessary to confirm the extent and effects of such editing in humans.

RNA modifications in pulmonary diseases

Threatening public health, pulmonary disease (PD) encompasses diverse lung injuries like chronic obstructive PD, pulmonary fibrosis, asthma, pulmonary infections due to pathogen invasion, and fatal lung cancer. The crucial involvement of RNA epigenetic modifications in PD pathogenesis is underscored by robust evidence. These modifications not only shape cell fates but also finely modulate the expression of genes linked to disease progression, suggesting their utility as biomarkers and targets for therapeutic strategies. The critical RNA modifications implicated in PDs are summarized in this review, including N6-methylation of adenosine, N1-methylation of adenosine, 5-methylcytosine, pseudouridine (5-ribosyl uracil), 7-methylguanosine, and adenosine to inosine editing, along with relevant regulatory mechanisms. By shedding light on the pathology of PDs, these summaries could spur the identification of new biomarkers and therapeutic strategies, ultimately paving the way for early PD diagnosis and treatment innovation.

Emerging roles of circular RNAs in nasopharyngeal carcinoma: functions and implications

Nasopharyngeal carcinoma (NPC) is a distinct malignancy primarily prevalent in Southern China and Southeast Asia. Circular RNAs (circRNAs), a class of non-coding RNAs, are evolutionarily conserved and exhibit remarkable stability. Their dysregulation has been observed in various cancers, including NPC. In this review, we investigate the pivotal role of circRNAs in NPC, focusing specifically on their involvement in tumor proliferation, apoptosis, metastasis, angiogenesis, stemness, metabolism, and the tumor microenvironment. We highlight the diagnostic and prognostic potential of circRNAs in NPC, emphasizing their utility as biomarkers for early detection, disease monitoring, and prediction of treatment outcomes. Additionally, we explore the therapeutic implications of circRNAs in NPC, highlighting their potential for targeted therapies.

New comparative genomic evidence supporting the proteomic diversification role of A-to-I RNA editing in insects

Adenosine-to-inosine (A-to-I) RNA editing, resembling A-to-G mutation, confers adaptiveness by increasing proteomic diversity in a temporal-spatial manner. This evolutionary theory named "proteomic diversifying hypothesis" has only partially been tested in very few organisms like Drosophila melanogaster, mainly by observing the positive selection on nonsynonymous editing events. To find additional genome-wide evidences supporting this interesting assumption, we retrieved the genomes of four Drosophila species and collected 20 deep-sequenced transcriptomes of different developmental stages and neuron populations of D. melanogaster. We systematically profiled the RNA editomes in these samples and performed meticulous comparative genomic analyses. Further evidences were found to support the diversifying hypothesis. (1) None of the nonsynonymous editing sites in D. melanogaster had ancestral G-alleles, while the silent editing sites had an unignorable fraction of ancestral G-alleles; (2) Only very few nonsynonymous editing sites in D. melanogaster had corresponding G-alleles derived in the genomes of sibling species, and the fraction of such situation was significantly lower than that of silent editing sites; (3) The few nonsynonymous editing with corresponding G-alleles had significantly more variable editing levels (across samples) than other nonsynonymous editing sites in D. melanogaster. The proteomic diversifying nature of RNA editing in Drosophila excludes the restorative role which favors an ancestral G-allele. The few fixed G-alleles in sibling species might facilitate the adaptation to particular environment and the corresponding nonsynonymous editing in D. melanogaster would introduce stronger advantage of flexible proteomic diversification. With multi-Omics data, our study consolidates the nature of evolutionary significance of A-to-I RNA editing sites in model insects.

Single cell analyses of cancer cells identified two regulatorily and functionally distinct categories in differentially expressed genes among tumor subclones

To explore the feature of cancer cells and tumor subclones, we analyzed 101,065 single-cell transcriptomes from 12 colorectal cancer (CRC) patients and 92 single cell genomes from one of these patients. We found cancer cells, endothelial cells and stromal cells in tumor tissue expressed much more genes and had stronger cell-cell interactions than their counterparts in normal tissue. We identified copy number variations (CNVs) in each cancer cell and found correlation between gene copy number and expression level in cancer cells at single cell resolution. Analysis of tumor subclones inferred by CNVs showed accumulation of mutations in each tumor subclone along lineage trajectories. We found differentially expressed genes (DEGs) between tumor subclones had two populations: DEGCNV and DEGreg. DEGCNV, showing high CNV-expression correlation and whose expression differences depend on the differences of CNV level, enriched in housekeeping genes and cell adhesion associated genes. DEGreg, showing low CNV-expression correlation and mainly in low CNV variation regions and regions without CNVs, enriched in cytokine signaling genes. Furthermore, cell-cell communication analyses showed that DEGCNV tends to involve in cell-cell contact while DEGreg tends to involve in secreted signaling, which further support that DEGCNV and DEGreg are two regulatorily and functionally distinct categories.

Innate immune sensing of macromolecule homeostasis

The innate immune system uses a distinct set of germline-encoded pattern recognition receptors to recognize molecular patterns initially thought to be unique to microbial invaders, named pathogen-associated molecular patterns. The concept was later further developed to include similar molecular patterns originating from host cells during tissue damage, known as damage-associated molecular patterns. However, recent advances in the mechanism of monogenic inflammatory diseases have highlighted a much more expansive repertoire of cellular functions that are monitored by innate immunity. Here, we summarize several examples in which an innate immune response is triggered when homeostasis of macromolecule in the cell is disrupted in non-infectious or sterile settings. These ever-growing sensing mechanisms expand the repertoire of innate immune recognition, positioning it not only as a key player in host defense but also as a gatekeeper of cellular homeostasis. Therapeutics inspired by these advances to restore cellular homeostasis and correct the immune system could have far-reaching implications.

RNA modifications in cellular metabolism: implications for metabolism-targeted therapy and immunotherapy

Cellular metabolism is an intricate network satisfying bioenergetic and biosynthesis requirements of cells. Relevant studies have been constantly making inroads in our understanding of pathophysiology, and inspiring development of therapeutics. As a crucial component of epigenetics at post-transcription level, RNA modification significantly determines RNA fates, further affecting various biological processes and cellular phenotypes. To be noted, immunometabolism defines the metabolic alterations occur on immune cells in different stages and immunological contexts. In this review, we characterize the distribution features, modifying mechanisms and biological functions of 8 RNA modifications, including N6-methyladenosine (m6A), N6,2'-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), 5-methylcytosine (m5C), N4-acetylcytosine (ac4C), N7-methylguanosine (m7G), Pseudouridine (Ψ), adenosine-to-inosine (A-to-I) editing, which are relatively the most studied types. Then regulatory roles of these RNA modification on metabolism in diverse health and disease contexts are comprehensively described, categorized as glucose, lipid, amino acid, and mitochondrial metabolism. And we highlight the regulation of RNA modifications on immunometabolism, further influencing immune responses. Above all, we provide a thorough discussion about clinical implications of RNA modification in metabolism-targeted therapy and immunotherapy, progression of RNA modification-targeted agents, and its potential in RNA-targeted therapeutics. Eventually, we give legitimate perspectives for future researches in this field from methodological requirements, mechanistic insights, to therapeutic applications.

Malignant A-to-I RNA editing by ADAR1 drives T cell acute lymphoblastic leukemia relapse via attenuating dsRNA sensing

Leukemia-initiating cells (LICs) are regarded as the origin of leukemia relapse and therapeutic resistance. Identifying direct stemness determinants that fuel LIC self-renewal is critical for developing targeted approaches. Here, we show that the RNA-editing enzyme ADAR1 is a crucial stemness factor that promotes LIC self-renewal by attenuating aberrant double-stranded RNA (dsRNA) sensing. Elevated adenosine-to-inosine editing is a common attribute of relapsed T cell acute lymphoblastic leukemia (T-ALL) regardless of molecular subtype. Consequently, knockdown of ADAR1 severely inhibits LIC self-renewal capacity and prolongs survival in T-ALL patient-derived xenograft models. Mechanistically, ADAR1 directs hyper-editing of immunogenic dsRNA to avoid detection by the innate immune sensor melanoma differentiation-associated protein 5 (MDA5). Moreover, we uncover that the cell-intrinsic level of MDA5 dictates the dependency on the ADAR1-MDA5 axis in T-ALL. Collectively, our results show that ADAR1 functions as a self-renewal factor that limits the sensing of endogenous dsRNA. Thus, targeting ADAR1 presents an effective therapeutic strategy for eliminating T-ALL LICs.

The role of ADAR1 through and beyond its editing activity in cancer

Adenosine-to-inosine (A-to-I) editing of RNA, catalyzed by adenosine deaminase acting on RNA (ADAR) enzymes, is a prevalent RNA modification in mammals. It has been shown that A-to-I editing plays a critical role in multiple diseases, such as cardiovascular disease, neurological disorder, and particularly cancer. ADARs are the family of enzymes, including ADAR1, ADAR2, and ADAR3, that catalyze the occurrence of A-to-I editing. Notably, A-to-I editing is mainly catalyzed by ADAR1. Given the significance of A-to-I editing in disease development, it is important to unravel the complex roles of ADAR1 in cancer for the development of novel therapeutic interventions.In this review, we briefly describe the progress of research on A-to-I editing and ADARs in cancer, mainly focusing on the role of ADAR1 in cancer from both editing-dependent and independent perspectives. In addition, we also summarized the factors affecting the expression and editing activity of ADAR1 in cancer.

RMBase v3.0: decode the landscape, mechanisms and functions of RNA modifications

Although over 170 chemical modifications have been identified, their prevalence, mechanism and function remain largely unknown. To enable integrated analysis of diverse RNA modification profiles, we have developed RMBase v3.0 (http://bioinformaticsscience.cn/rmbase/), a comprehensive platform consisting of eight modules. These modules facilitate the exploration of transcriptome-wide landscape, biogenesis, interactome and functions of RNA modifications. By mining thousands of epitranscriptome datasets with novel pipelines, the 'RNA Modifications' module reveals the map of 73 RNA modifications of 62 species. the 'Genes' module allows to retrieve RNA modification profiles and clusters by gene and transcript. The 'Mechanisms' module explores 23 382 enzyme-catalyzed or snoRNA-guided modified sites to elucidate their biogenesis mechanisms. The 'Co-localization' module systematically formulates potential correlations between 14 histone modifications and 6 RNA modifications in various cell-lines. The 'RMP' module investigates the differential expression profiles of 146 RNA-modifying proteins (RMPs) in 18 types of cancers. The 'Interactome' integrates the interactional relationships between 73 RNA modifications with RBP binding events, miRNA targets and SNPs. The 'Motif' illuminates the enriched motifs for 11 types of RNA modifications identified from epitranscriptome datasets. The 'Tools' introduces a novel web-based 'modGeneTool' for annotating modifications. Overall, RMBase v3.0 provides various resources and tools for studying RNA modifications.

Liver cancer-specific mutations in functional domains of ADAR2 lead to the elevation of coding and non-coding RNA editing in multiple tumor-related genes

Mutation is the major cause of phenotypic innovations. Apart from DNA mutations, the alteration on RNA such as the ADAR-mediated A-to-I RNA editing could also shape the phenotype. These two layers of variations have not been systematically combined to study their collective roles in cancers. We collected the high-quality transcriptomes of ten hepatocellular carcinoma (HCC) and the matched control samples. We systematically identified HCC-specific mutations in the exonic regions and profiled the A-to-I RNA editome in each sample. All ten HCC samples had mutations in the CDS of ADAR2 gene (dsRNA-binding domain or catalytic domain). The consequence of these mutations converged to the elevation of ADAR2 efficiency as reflected by the global increase of RNA editing levels in HCC. The up-regulated editing sites (UES) were enriched in the CDS and UTR of oncogenes and tumor suppressor genes (TSG), indicating the possible roles of these target genes in HCC oncogenesis. We present the mutation-ADAR2-UES-oncogene/TSG-HCC axis that explains how mutations at different layers would finally lead to abnormal phenotype. In the light of central dogma, our work provides novel insights into how to fully take advantage of the transcriptome data to decipher the consequence of mutations.

Harnessing ADAR-Mediated Site-Specific RNA Editing in Immune-Related Disease: Prediction and Therapeutic Implications

ADAR (Adenosine Deaminases Acting on RNA) proteins are a group of enzymes that play a vital role in RNA editing by converting adenosine to inosine in RNAs. This process is a frequent post-transcriptional event observed in metazoan transcripts. Recent studies indicate widespread dysregulation of ADAR-mediated RNA editing across many immune-related diseases, such as human cancer. We comprehensively review ADARs' function as pattern recognizers and their capability to contribute to mediating immune-related pathways. We also highlight the potential role of site-specific RNA editing in maintaining homeostasis and its relationship to various diseases, such as human cancers. More importantly, we summarize the latest cutting-edge computational approaches and data resources for predicting and analyzing RNA editing sites. Lastly, we cover the recent advancement in site-directed ADAR editing tool development. This review presents an up-to-date overview of ADAR-mediated RNA editing, how site-specific RNA editing could potentially impact disease pathology, and how they could be harnessed for therapeutic applications.

Dissecting the basis for differential substrate specificity of ADAR1 and ADAR2

Millions of adenosines are deaminated throughout the transcriptome by ADAR1 and/or ADAR2 at varying levels, raising the question of what are the determinants guiding substrate specificity and how these differ between the two enzymes. We monitor how secondary structure modulates ADAR2 vs ADAR1 substrate selectivity, on the basis of systematic probing of thousands of synthetic sequences transfected into cell lines expressing exclusively ADAR1 or ADAR2. Both enzymes induce symmetric, strand-specific editing, yet with distinct offsets with respect to structural disruptions: -26 nt for ADAR2 and -35 nt for ADAR1. We unravel the basis for these differences in offsets through mutants, domain-swaps, and ADAR homologs, and find it to be encoded by the differential RNA binding domain (RBD) architecture. Finally, we demonstrate that this offset-enhanced editing can allow an improved design of ADAR2-recruiting therapeutics, with proof-of-concept experiments demonstrating increased on-target and potentially decreased off-target editing.

Identification and Interpretation of A-to-I RNA Editing Events in Insect Transcriptomes

Adenosine-to-inosine (A-to-I) RNA editing is the most prevalent RNA modification in the nervous systems of metazoans. To study the biological significance of RNA editing, we first have to accurately identify these editing events from the transcriptome. The genome-wide identification of RNA editing sites remains a challenging task. In this review, we will first introduce the occurrence, regulation, and importance of A-to-I RNA editing and then describe the established bioinformatic procedures and difficulties in the accurate identification of these sit esespecially in small sized non-model insects. In brief, (1) to obtain an accurate profile of RNA editing sites, a transcriptome coupled with the DNA resequencing of a matched sample is favorable; (2) the single-cell sequencing technique is ready to be applied to RNA editing studies, but there are a few limitations to overcome; (3) during mapping and variant calling steps, various issues, like mapping and base quality, soft-clipping, and the positions of mismatches on reads, should be carefully considered; (4) Sanger sequencing of both RNA and the matched DNA is the best verification of RNA editing sites, but other auxiliary evidence, like the nonsynonymous-to-synonymous ratio or the linkage information, is also helpful for judging the reliability of editing sites. We have systematically reviewed the understanding of the biological significance of RNA editing and summarized the methodology for identifying such editing events. We also raised several promising aspects and challenges in this field. With insightful perspectives on both scientific and technical issues, our review will benefit the researchers in the broader RNA editing community.

ADAR1p150 prevents MDA5 and PKR activation via distinct mechanisms to avert fatal autoinflammation

Effective immunity requires the innate immune system to distinguish foreign nucleic acids from cellular ones. Cellular double-stranded RNAs (dsRNAs) are edited by the RNA-editing enzyme ADAR1 to evade being recognized as viral dsRNA by cytoplasmic dsRNA sensors, including MDA5 and PKR. The loss of ADAR1-mediated RNA editing of cellular dsRNA activates MDA5. Additional RNA-editing-independent functions of ADAR1 have been proposed, but a specific mechanism has not been delineated. We now demonstrate that the loss of ADAR1-mediated RNA editing specifically activates MDA5, whereas loss of the cytoplasmic ADAR1p150 isoform or its dsRNA-binding activity enabled PKR activation. Deleting both MDA5 and PKR resulted in complete rescue of the embryonic lethality of Adar1p150-/- mice to adulthood, contrasting with the limited or no rescue by removing MDA5 or PKR alone. Our findings demonstrate that MDA5 and PKR are the primary in vivo effectors of fatal autoinflammation following the loss of ADAR1p150.

dsRID: in silico identification of dsRNA regions using long-read RNA-seq data

MOTIVATION

Double-stranded RNAs (dsRNAs) are potent triggers of innate immune responses upon recognition by cytosolic dsRNA sensor proteins. Identification of endogenous dsRNAs helps to better understand the dsRNAome and its relevance to innate immunity related to human diseases.

RESULTS

Here, we report dsRID (double-stranded RNA identifier), a machine-learning-based method to predict dsRNA regions in silico, leveraging the power of long-read RNA-sequencing (RNA-seq) and molecular traits of dsRNAs. Using models trained with PacBio long-read RNA-seq data derived from Alzheimer's disease (AD) brain, we show that our approach is highly accurate in predicting dsRNA regions in multiple datasets. Applied to an AD cohort sequenced by the ENCODE consortium, we characterize the global dsRNA profile with potentially distinct expression patterns between AD and controls. Together, we show that dsRID provides an effective approach to capture global dsRNA profiles using long-read RNA-seq data.

AVAILABILITY AND IMPLEMENTATION

Software implementation of dsRID, and genomic coordinates of regions predicted by dsRID in all samples are available at the GitHub repository: https://github.com/gxiaolab/dsRID.

Library Screening Reveals Sequence Motifs That Enable ADAR2 Editing at Recalcitrant Sites

Adenosine deaminases acting on RNA (ADARs) catalyze the hydrolytic deamination of adenosine to inosine in duplex RNA. The inosine product preferentially base pairs with cytidine resulting in an effective A-to-G edit in RNA. ADAR editing can result in a recoding event alongside other alterations to RNA function. A consequence of ADARs' selective activity on duplex RNA is that guide RNAs (gRNAs) can be designed to target an adenosine of interest and promote a desired recoding event. One of ADAR's main limitations is its preference to edit adenosines with specific 5' and 3' nearest neighbor nucleotides (e.g., 5' U, 3' G). Current rational design approaches are well-suited for this ideal sequence context, but limited when applied to difficult-to-edit sites. Here we describe a strategy for the in vitro evaluation of very large libraries of ADAR substrates (En Masse Evaluation of RNA Guides, EMERGe). EMERGe allows for a comprehensive screening of ADAR substrate RNAs that complements current design approaches. We used this approach to identify sequence motifs for gRNAs that enable editing in otherwise difficult-to-edit target sites. A guide RNA bearing one of these sequence motifs enabled the cellular repair of a premature termination codon arising from mutation of the MECP2 gene associated with Rett Syndrome. EMERGe provides an advancement in screening that not only allows for novel gRNA design, but also furthers our understanding of ADARs' specific RNA-protein interactions.