Computational approaches for detection and quantification of A-to-I RNA-editing

Site-specific quantification of Adenosine-to-Inosine RNA editing by Endonuclease-Mediated qPCR

RNA molecules contain diverse modified nucleobases that play pivotal roles in numerous biological processes. Adenosine-to-inosine (A-to-I) RNA editing, one of the most prevalent RNA modifications in mammalian cells, is linked to a multitude of human diseases. To unveil the functions of A-to-I RNA editing, accurate quantification of inosine at specific sites is essential. In this study, we developed an endonuclease-mediated cleavage and real-time fluorescence quantitative PCR method for A-to-I RNA editing (EM-qPCR) to quantitatively analyze A-to-I RNA editing at a single site. By employing this method, we successfully quantified the levels of A-to-I RNA editing on various transfer RNA (tRNA) molecules at position 34 (I34) in mammalian cells with precision. Subsequently, this method was applied to tissues from sleep-deprived mice, revealing a notable alteration in the levels of I34 between sleep-deprived and control mice. The proposed method sets a precedent for the quantitative analysis of A-to-I RNA editing at specific sites, facilitating a deeper understanding of the biological implications of A-to-I RNA editing.

Spatial Visualization of A-to-I Editing in Cells Using Endonuclease V Immunostaining Assay (EndoVIA)

Adenosine-to-inosine (A-to-I) editing is one of the most widespread post-transcriptional RNA modifications and is catalyzed by adenosine deaminases acting on RNA (ADARs). Varying across tissue types, A-to-I editing is essential for numerous biological functions, and dysregulation leads to autoimmune and neurological disorders, as well as cancer. Recent evidence has also revealed a link between RNA localization and A-to-I editing, yet understanding of the mechanisms underlying this relationship and its biological impact remains limited. Current methods rely primarily on in vitro characterization of extracted RNA that ultimately erases subcellular localization and cell-to-cell heterogeneity. To address these challenges, we have repurposed endonuclease V (EndoV), a magnesium-dependent ribonuclease that cleaves inosine bases in edited RNA, to selectively bind and detect A-to-I edited RNA in cells. The work herein introduces an endonuclease V immunostaining assay (EndoVIA), a workflow that provides spatial visualization of edited transcripts, enables rapid quantification of overall inosine abundance, and maps the landscape of A-to-I editing within the transcriptome at the nanoscopic level.

Transcriptome analysis of Citrus sinensis reveals potential responsive events triggered by Candidatus Liberibacter asiaticus

Citrus Huanglongbing (HLB), caused by Candidatus Liberibacter asiaticus (CLas), is a devastating immune-mediated disorder that has a detrimental effect on the citrus industry, with the distinguishing feature being an eruption of reactive oxygen species (ROS). This study explored the alterations in antioxidant enzyme activity, transcriptome, and RNA editing events of organelles in C. sinensis during CLas infection. Results indicated that there were fluctuations in the performance of antioxidant enzymes, such as ascorbate peroxidase (APX), catalase (CAT), glutathione reductase (GR), peroxidase (POD), and superoxide dismutase (SOD), in plants affected by HLB. Transcriptome analysis revealed 3604 genes with altered expression patterns between CLas-infected and healthy samples, including those associated with photosynthesis, biotic interactions, and phytohormones. Samples infected with CLas showed a decrease in the expression of most genes associated with photosynthesis and gibberellin metabolism. It was discovered that RNA editing frequency and the expression level of various genes in the chloroplast and mitochondrion genomes were affected by CLas infection. Our findings provide insights into the inhibition of photosynthesis, gibberellin metabolism, and antioxidant enzymes during CLas infection in C. sinensis.

Spatial visualization of A-to-I Editing in cells using Endonuclease V Immunostaining Assay (EndoVIA)

Adenosine-to-Inosine (A-to-I) editing is one of the most widespread post-transcriptional RNA modifications and is catalyzed by adenosine deaminases acting on RNA (ADARs). Varying across tissue types, A-to-I editing is essential for numerous biological functions and dysregulation leads to autoimmune and neurological disorders, as well as cancer. Recent evidence has also revealed a link between RNA localization and A-to-I editing, yet understanding of the mechanisms underlying this relationship and its biological impact remains limited. Current methods rely primarily on in vitro characterization of extracted RNA that ultimately erases subcellular localization and cell-to-cell heterogeneity. To address these challenges, we have repurposed Endonuclease V (EndoV), a magnesium dependent ribonuclease that cleaves inosine bases in edited RNA, to selectively bind and detect A-to-I edited RNA in cells. The work herein introduces Endonuclease V Immunostaining Assay (EndoVIA), a workflow that provides spatial visualization of edited transcripts, enables rapid quantification of overall inosine abundance, and maps the landscape of A-to-I editing within the transcriptome at the nanoscopic level.

Transcriptome-wide profiling of A-to-I RNA editing by Slic-seq

Adenosine-to-inosine (A-to-I) RNA editing is a post-transcriptional processing event involved in diversifying the transcriptome and is responsible for various biological processes. In this context, we developed a new method based on the highly selective cleavage activity of Endonuclease V against Inosine and the universal activity of sodium periodate against all RNAs to enrich the inosine-containing RNA and accurately identify the editing sites. We validated the reliability of our method in human brain in both Alu and non-Alu elements. The conserved sites of A-to-I editing in human cells (HEK293T, HeLa, HepG2, K562 and MCF-7) primarily occurs in the 3'UTR of the RNA, which are highly correlated with RNA binding and protein binding. Analysis of the editing sites between the human brain and mouse brain revealed that the editing of exons is more conserved than that in other regions. This method was applied to three neurological diseases (Alzheimer's, epilepsy and ageing) of mouse brain, reflecting that A-to-I editing sites significantly decreased in neuronal activity genes.

Profiling A-to-I RNA editing during mouse somatic reprogramming at the single-cell level

Mouse somatic cells can be reprogrammed into induced pluripotent stem cells through a highly heterogeneous process regulated by numerous biological factors, including adenosine-to-inosine (A-to-I) RNA editing. In this study, we analyzed A-to-I RNA editing sites using a single-cell RNA sequencing (scRNA-seq) dataset with high-depth and full-length coverage. Our method revealed that A-to-I RNA editing frequency varied widely at the single-cell level and underwent dynamic changes. We also found that A-to-I RNA editing level was correlated with the expression of the RNA editing enzyme ADAR1. The analysis combined with gene ontology (GO) enrichment revealed that ADAR1-dependent A-to-I editing may downregulate the expression levels of Igtp, Irgm2, Mndal, Ifi202b, and Tapbp in the early stage, to inhibit the pathways of cellular response to interferon-beta and regulation of protein complex stability to promote mesenchymal-epithelial transition (MET). Notably, we identified a negative correlation between A-to-I RNA editing frequency and the expression of certain genes, such as Nras, Ube2l6, Zfp987, and Adsl.

Quantitative and site-specific detection of inosine modification in RNA by acrylonitrile labeling-mediated elongation stalling

RNA molecules contain diverse modifications that play crucial roles in a wide variety of biological processes. Inosine is one of the most prevalent modifications in RNA and dysregulation of inosine is correlated with many human diseases. Herein, we established an acrylonitrile labeling-mediated elongation stalling (ALES) method for quantitative and site-specific detection of inosine in RNA from biological samples. In ALES method, inosine is selectively cyanoethylated with acrylonitrile to form N1-cyanoethylinosine (ce1I) through a Michael addition reaction. The N1-cyanoethyl group of ce1I compromises the hydrogen bond between ce1I and other nucleobases, leading to the stalling of reverse transcription at original inosine site. This specific property of stalling at inosine site could be evaluated by subsequent real-time quantitative PCR (qPCR). With the proposed ALES method, we found the significantly increased level of inosine at position Chr1:63117284 of Ino80dos RNA of multiple tissues from sleep-deprived mice compared to the control mice. This is the first report on the investigation of inosine modification in sleep-deprived mice, which may open up new direction for deciphering insomnia from RNA modifications. In addition, we found the decreased level of inosine at GluA2 Q/R site (Chr4:157336723) in glioma tissues, indicating the decreased level of inosine at GluA2 Q/R site may serve as potential indicator for the diagnosis of glioma. Taken together, the proposed ALES method is capable of quantitative and site-specific detection of inosine in RNA, which provides a valuable tool to uncover the functions of inosine in human diseases.

Single-Base Resolution Detection of Adenosine-to-Inosine RNA Editing by Endonuclease-Mediated Sequencing

RNA molecules contain diverse modifications that play crucial roles in a wide variety of biological processes. Adenosine-to-inosine (A-to-Ino) RNA editing is one of the most prevalent modifications among all types of RNA. Abnormal A-to-InoRNA editing has been demonstrated to be associated with many human diseases. Identification of A-to-Ino editing sites is indispensable to deciphering their biological roles. Herein, by employing the unique property of human endonuclease V (hEndoV), we proposed a hEndoV-mediated sequencing (hEndoV-seq) method for the single-base resolution detection of A-to-InoRNA editing sites. In this approach, the terminal 3'OH of RNA is first blocked by 3'-deoxyadenosine (3'-deoxy-A). Specific cleavage of Ino sites by hEndoV protein produces new terminal 3'OH, which can be identified by sequencing analysis, and therefore offers the site-specific detection of Ino in RNA. The principle of hEndoV-seq is straightforward and the analytical procedure is simple. No chemical reaction is involved in the sequencing library preparation. The whole procedure in hEndoV-seq is carried out under mild conditions and RNA is not prone to degradation. Taken together, the proposed hEndoV-seq method is capable of site-specific identification of A-to-Ino editing in RNA, which provides a valuable tool for elucidating the functions of A-to-Ino editing in RNA.

Direct identification of A-to-I editing sites with nanopore native RNA sequencing

Inosine is a prevalent RNA modification in animals and is formed when an adenosine is deaminated by the ADAR family of enzymes. Traditionally, inosines are identified indirectly as variants from Illumina RNA-sequencing data because they are interpreted as guanosines by cellular machineries. However, this indirect method performs poorly in protein-coding regions where exons are typically short, in non-model organisms with sparsely annotated single-nucleotide polymorphisms, or in disease contexts where unknown DNA mutations are pervasive. Here, we show that Oxford Nanopore direct RNA sequencing can be used to identify inosine-containing sites in native transcriptomes with high accuracy. We trained convolutional neural network models to distinguish inosine from adenosine and guanosine, and to estimate the modification rate at each editing site. Furthermore, we demonstrated their utility on the transcriptomes of human, mouse and Xenopus. Our approach expands the toolkit for studying adenosine-to-inosine editing and can be further extended to investigate other RNA modifications.

Virus-specific editing identification approach reveals the landscape of A-to-I editing and its impacts on SARS-CoV-2 characteristics and evolution

Upon SARS-CoV-2 infection, viral intermediates specifically activate the IFN response through MDA5-mediated sensing and accordingly induce ADAR1 p150 expression, which might lead to viral A-to-I RNA editing. Here, we developed an RNA virus-specific editing identification pipeline, surveyed 7622 RNA-seq data from diverse types of samples infected with SARS-CoV-2, and constructed an atlas of A-to-I RNA editing sites in SARS-CoV-2. We found that A-to-I editing was dynamically regulated, varied between tissue and cell types, and was correlated with the intensity of innate immune response. On average, 91 editing events were deposited at viral dsRNA intermediates per sample. Moreover, editing hotspots were observed, including recoding sites in the spike gene that affect viral infectivity and antigenicity. Finally, we provided evidence that RNA editing accelerated SARS-CoV-2 evolution in humans during the epidemic. Our study highlights the ability of SARS-CoV-2 to hijack components of the host antiviral machinery to edit its genome and fuel its evolution, and also provides a framework and resource for studying viral RNA editing.

Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA

Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.

General Principles for the Detection of Modified Nucleotides in RNA by Specific Reagents

Epitranscriptomics heavily rely on chemical reagents for the detection, quantification, and localization of modified nucleotides in transcriptomes. Recent years have seen a surge in mapping methods that use innovative and rediscovered organic chemistry in high throughput approaches. While this has brought about a leap of progress in this young field, it has also become clear that the different chemistries feature variegated specificity and selectivity. The associated error rates, e.g., in terms of false positives and false negatives, are in large part inherent to the chemistry employed. This means that even assuming technically perfect execution, the interpretation of mapping results issuing from the application of such chemistries are limited by intrinsic features of chemical reactivity. An important but often ignored fact is that the huge stochiometric excess of unmodified over-modified nucleotides is not inert to any of the reagents employed. Consequently, any reaction aimed at chemical discrimination of modified versus unmodified nucleotides has optimal conditions for selectivity that are ultimately anchored in relative reaction rates, whose ratio imposes intrinsic limits to selectivity. Here chemical reactivities of canonical and modified ribonucleosides are revisited as a basis for an understanding of the limits of selectivity achievable with chemical methods.

RESIC: A Tool for Comprehensive Adenosine to Inosine RNA Editing Site Identification and Classification

Adenosine to inosine (A-to-I) RNA editing, the most prevalent type of RNA editing in metazoans, is carried out by adenosine deaminases (ADARs) in double-stranded RNA regions. Several computational approaches have been recently developed to identify A-to-I RNA editing sites from sequencing data, each addressing a particular issue. Here, we present RNA Editing Sites Identification and Classification (RESIC), an efficient pipeline that combines several approaches for the detection and classification of RNA editing sites. The pipeline can be used for all organisms and can use any number of RNA-sequencing datasets as input. RESIC provides (1) the detection of editing sites in both repetitive and non-repetitive genomic regions; (2) the identification of hyper-edited regions; and (3) optional exclusion of polymorphism sites to increase reliability, based on DNA, and ADAR-mutant RNA sequencing datasets, or SNP databases. We demonstrate the utility of RESIC by applying it to human, successfully overlapping and extending the list of known putative editing sites. We further tested changes in the patterns of A-to-I RNA editing, and RNA abundance of ADAR enzymes, following SARS-CoV-2 infection in human cell lines. Our results suggest that upon SARS-CoV-2 infection, compared to mock, the number of hyper editing sites is increased, and in agreement, the activity of ADAR1, which catalyzes hyper-editing, is enhanced. These results imply the involvement of A-to-I RNA editing in conceiving the unpredicted phenotype of COVID-19 disease. RESIC code is open-source and is easily extendable.

Systematic identification of A-to-I RNA editing in zebrafish development and adult organs

A-to-I RNA editing is a common post transcriptional mechanism, mediated by the Adenosine deaminase that acts on RNA (ADAR) enzymes, that increases transcript and protein diversity. The study of RNA editing is limited by the absence of editing maps for most model organisms, hindering the understanding of its impact on various physiological conditions. Here, we mapped the vertebrate developmental landscape of A-to-I RNA editing, and generated the first comprehensive atlas of editing sites in zebrafish. Tens of thousands unique editing events and 149 coding sites were identified with high-accuracy. Some of these edited sites are conserved between zebrafish and humans. Sequence analysis of RNA over seven developmental stages revealed high levels of editing activity in early stages of embryogenesis, when embryos rely on maternal mRNAs and proteins. In contrast to the other organisms studied so far, the highest levels of editing were detected in the zebrafish ovary and testes. This resource can serve as the basis for understanding of the role of editing during zebrafish development and maturity.

Discovering RNA Editing Events in Fungi

RNA editing is an important posttranscriptional process that alters the genetic information of RNA encoded by genomic DNA. Adenosine-to-inosine (A-to-I) editing is the most prevalent type of RNA editing in animal kingdom, catalyzed by adenosine deaminases acting on RNA (ADARs). Recently, genome-wide A-to-I RNA editing is discovered in fungi, involving adenosine deamination mechanisms distinct from animals. Aiming to draw more attention to RNA editing in fungi, here we discuss the considerations for deep sequencing data preparation and the available various methods for detecting RNA editing, with a special emphasis on their usability for fungal RNA editing detection. We describe computational protocols for the identification of candidate RNA editing sites in fungi by using two software packages REDItools and RES-Scanner with RNA sequencing (RNA-Seq) and genomic DNA sequencing (DNA-Seq) data.

RNA modifications as a common denominator between tRNA and mRNA

Recent studies underscore RNA modifications as a novel mechanism to coordinate expression and function of different genes. While modifications on the sugar or base moieties of tRNA are well known, their roles in mRNA regulation are only starting to emerge. Interestingly, some modifications are present in both tRNA and mRNA, and here we discuss the functional significance of these common features. We describe key modifications that are present in both RNA types, elaborate on proteins that interact with them, and indicate recent works that identify roles in communicating tRNA processes and mRNA regulation. We propose that as tools are developed, the shortlist of features that are common between types of RNA will greatly expand and proteins that interact with them will be identified. In conclusion, the presence of the same modification in both RNA types provides an intersect between tRNA processes and mRNA regulation and implies a novel mechanism for connecting diverse cellular processes.

A-to-I RNA Editing in Cancer: From Evaluating the Editing Level to Exploring the Editing Effects

As an important regulatory mechanism at the posttranscriptional level in metazoans, adenosine deaminase acting on RNA (ADAR)-induced A-to-I RNA editing modification of double-stranded RNA has been widely detected and reported. Editing may lead to non-synonymous amino acid mutations, RNA secondary structure alterations, pre-mRNA processing changes, and microRNA-mRNA redirection, thereby affecting multiple cellular processes and functions. In recent years, researchers have successfully developed several bioinformatics software tools and pipelines to identify RNA editing sites. However, there are still no widely accepted editing site standards due to the variety of parallel optimization and RNA high-seq protocols and programs. It is also challenging to identify RNA editing by normal protocols in tumor samples due to the high DNA mutation rate. Numerous RNA editing sites have been reported to be located in non-coding regions and can affect the biosynthesis of ncRNAs, including miRNAs and circular RNAs. Predicting the function of RNA editing sites located in non-coding regions and ncRNAs is significantly difficult. In this review, we aim to provide a better understanding of bioinformatics strategies for human cancer A-to-I RNA editing identification and briefly discuss recent advances in related areas, such as the oncogenic and tumor suppressive effects of RNA editing.

Profiling neural editomes reveals a molecular mechanism to regulate RNA editing during development

Adenosine (A) to inosine (I) RNA editing contributes to transcript diversity and modulates gene expression in a dynamic, cell type–specific manner. During mammalian brain development, editing of specific adenosines increases, whereas the expression of A-to-I editing enzymes remains unchanged, suggesting molecular mechanisms that mediate spatiotemporal regulation of RNA editing exist. Herein, by using a combination of biochemical and genomic approaches, we uncover a molecular mechanism that regulates RNA editing in a neural- and development-specific manner. Comparing editomes during development led to the identification of neural transcripts that were edited only in one life stage. The stage-specific editing is largely regulated by differential gene expression during neural development. Proper expression of nearly one-third of the neurodevelopmentally regulated genes is dependent on adr-2, the sole A-to-I editing enzyme in C. elegans. However, we also identified a subset of neural transcripts that are edited and expressed throughout development. Despite a neural-specific down-regulation of adr-2 during development, the majority of these sites show increased editing in adult neural cells. Biochemical data suggest that ADR-1, a deaminase-deficient member of the adenosine deaminase acting on RNA (ADAR) family, is competing with ADR-2 for binding to specific transcripts early in development. Our data suggest a model in which during neural development, ADR-2 levels overcome ADR-1 repression, resulting in increased ADR-2 binding and editing of specific transcripts. Together, our findings reveal tissue- and development-specific regulation of RNA editing and identify a molecular mechanism that regulates ADAR substrate recognition and editing efficiency.

The ADAR Family in Amphioxus: RNA Editing and Conserved Orthologous Site Predictions

RNA editing is a relatively unexplored process in which transcribed RNA is modified at specific nucleotides before translation, adding another level of regulation of gene expression. Cephalopods use it extensively to increase the regulatory complexity of their nervous systems, and mammals use it too, but less prominently. Nevertheless, little is known about the specifics of RNA editing in most of the other clades and the relevance of RNA editing from an evolutionary perspective remains unknown. Here we analyze a key element of the editing machinery, the ADAR (adenosine deaminase acting on RNA) gene family, in an animal with a key phylogenetic position at the root of chordates: the cephalochordate amphioxus. We show, that as in cephalopods, ADAR genes in amphioxus are predominantly expressed in the nervous system; we identify a number of RNA editing events in amphioxus; and we provide a newly developed method to identify RNA editing events in highly polymorphic genomes using orthology as a guide. Overall, our work lays the foundations for future comparative analysis of RNA-editing events across the metazoan tree.

Chemical Profiling of A-to-I RNA Editing Using a Click-Compatible Phenylacrylamide

Straightforward methods for detecting adenosine-to-inosine (A-to-I) RNA editing are key to a better understanding of its regulation, function, and connection with disease. We address this need by developing a novel reagent, N-(4-ethynylphenyl)acrylamide (EPhAA), and illustrating its ability to selectively label inosine in RNA. EPhAA is synthesized in a single step, reacts rapidly with inosine, and is "click"-compatible, enabling flexible attachment of fluorescent probes at editing sites. We first validate EPhAA reactivity and selectivity for inosine in both ribonucleosides and RNA substrates, and then apply our approach to directly monitor in vitro A-to-I RNA editing activity using recombinant ADAR enzymes. This method improves upon existing inosine chemical-labeling techniques and provides a cost-effective, rapid, and non-radioactive approach for detecting inosine formation in RNA. We envision this method will improve the study of A-to-I editing and enable better characterization of RNA modification patterns in different settings.

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.

Reading Chemical Modifications in the Transcriptome

Diverse chemical modifications have been identified in the transcriptome, leading to the emerging field of epitranscriptomics. In eukaryotic mRNA, the 5' cap and 3' poly(A) tail play important roles in regulation, and multiple internal modifications have also been revealed to participate in RNA metabolism. In this review, we focus on internal modifications in eukaryotic mRNA, including modifications to A/U/C/G bases and to ribose as well. We provide an overview of their biogenesis, high-throughput detection methods, biological functions, and regulatory mechanisms, with an emphasis on their reported reader proteins (RNA-binding proteins that specifically bind to modified RNA). We also briefly discuss the current problems in the investigation of mRNA modifications that need to be solved.