Human coding RNA editing is generally nonadaptive

RNA editing is phenomenon that occurs in both protein coding and non-coding RNAs. Increasing evidence have shown that adenosine-to-inosine RNA editing can potentially rendering substantial functional effects throughout the genome. Using RNA editing datasets from two large consortiums: The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) project, we quantitatively analyzed human genome-wide RNA editing events derived from tumor or normal tissues. Generally, a common RNA editing site tends to have a higher editing level in tumors as compared to normal samples. Of the 14 tumor-normal-paired cancer types examined, Eleven of the 14 cancers tested had overall increased RNA editing levels in the tumors. The editomes in cancer or normal tissues were dissected by genomic locations, and significant RNA editing locational difference was found between cancerous and healthy subjects. Additionally, our results indicated a significant correlation between the RNA editing rate and the gene density across chromosomes, highlighted hyper RNA editing clusters through visualization of running RNA editing rates along chromosomes, and identified hyper RNA edited genes (protein-coding genes, lincRNAs, and pseudogenes) that embody a large portion of cancer prognostic predictors. This study reinforces the potential functional effects of RNA editing in protein-coding genes, and also makes a strong foundation for further exploration of RNA editing’s roles in non-coding regions.

Background. Despite recent advances in surgical techniques and treatment options, gastric cancer (GC) remains the third most common cause of cancer-related deaths worldwide. Besides DNA sequence mutations, epigenetic alterations have emerged as significantly major players in cancer development. A-to-I RNA editing is a post-transcriptional modification that converts adenosines to inosines in both coding and noncoding RNA transcripts, and has recently been recognized has a novel epigenetic mechanism in GC pathogenesis. More specifically, A-to-I editing of AZIN1 transcripts was shown to be regulated by an adenosine deaminase acting on RNA-1 (ADAR1), and edited AZIN1 resulted in an aggressive phenotype during disease progression in some human cancers. The aim of this study was to clarify the clinical consequences of RNA editing status of AZIN1 and the RNA editing enzymes (ADAR1 and 2) in GC patients. Methods. Two hundred eighty-eight gastric specimens from one hundred forty-four patients who underwent surgery for GC were evaluated. We analyzed the RNA editing status of AZIN1 by RNA editing site-specific quantitative PCR (RESSqPCR). RESSqPCR allows quantitation of RNA editing levels using wild-type or edited AZIN1-specific primers, and RNA editing levels are calculated by determining the ratios of Ct values between edited vs. wild-type transcript expression levels. Furthermore, expression levels of ADAR1 and ADAR2 were evaluated by qPCR in GC tissues. Results. We observed a higher frequency of the AZIN1 RNA editing in tumors compared with normal mucosa in GC. RNA editing status of AZIN1 was significantly increased in a stage-dependent manner, and high frequency of RNA editing in AZIN1 significantly correlated with advanced T stage and presence of lymph node metastasis in GC patients. Multivariate analysis revealed that high frequency of RNA editing in the AZIN1 gene was an independent prognostic factor for poor disease free survival and overall survival. Furthermore, significant upregulation of ADAR1 and downregulation of ADAR2 was also observed in GC tissues compared with matched normal mucosa, and these alterations also significantly correlated with disease progression factors and poor prognosis in GC patients. Interestingly, ADAR1 expression level was positively correlated with RNA editing status of AZIN1 in GC tissues. Conclusions. Our findings revealed that altered gene-specific A-to-I editing events mediated by ADAR1 promote disease progression in GC. We conclude that assessment of RNA editing status of the AZIN1 gene might offer a novel and superior biomarker for a more accurate diagnosis of disease staging and may help predict clinical outcomes in GC patients. Citation Format: Yoshinaga Okugawa, Yuji Toiyama, Kunitoshi Shigeyasu, Takashi Ichikawa, Satoshi Oki, Koichiro Mori, Yuka Nagano, Hiromi Yasuda, Shigeyuki Yoshiyama, Masaki Ohi, Koji Tanaka, Yasuhiro Inoue, Toshimitsu Araki, Yasuhiko Mohri, Motoyoshi Tanaka, Chikao Miki, Ajay Goel, Masato Kusunoki. Novel findings for the clinical significance of RNA editing status of AZIN1 and ADAR 1 and 2 expression levels in gastric cancer patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4720. doi:10.1158/1538-7445.AM2017-4720

<div>Abstract<p>RNA editing modifies single nucleotides of RNAs, regulating primary protein structure and protein abundance. In recent years, the diversity of proteins and complexity of gene regulation associated with RNA editing dysregulation has been increasingly appreciated in oncology. Large-scale shifts in editing have been observed in bulk tumors across various cancer types. However, RNA editing in single cells and individual cell types within tumors has not been explored. By profiling editing in single cells from lung adenocarcinoma biopsies, we found that the increased editing trend of bulk lung tumors was unique to cancer cells. Elevated editing levels were observed in cancer cells resistant to targeted therapy, and editing sites associated with drug response were enriched. Consistent with the regulation of antiviral pathways by RNA editing, higher editing levels in cancer cells were associated with reduced antitumor innate immune response, especially levels of natural killer cell infiltration. In addition, the level of RNA editing in cancer cells was positively associated with somatic point mutation burden. This observation motivated the definition of a new metric, RNA editing load, reflecting the amount of RNA mutations created by RNA editing. Importantly, in lung cancer, RNA editing load was a stronger predictor of patient survival than DNA mutations. This study provides the first single cell dissection of editing in cancer and highlights the significance of RNA editing load in cancer prognosis.</p>Significance:<p>RNA editing analysis in single lung adenocarcinoma cells uncovers RNA mutations that correlate with tumor mutation burden and cancer innate immunity and reveals the amount of RNA mutations that strongly predicts patient survival.</p><p><i><a href="https://aacrjournals.org/cancerres/article/doi/10.1158/0008-5472.CAN-22-3537" target="_blank">See related commentary by Luo and Liang, p. 351</a></i></p></div>

BackgroundThe hair coat is available for the yak to live in the harsh environment of the plateau. Besides, improving the hair production of yak is necessary for its textile industry development. Hair grows from hair follicles (HFs). The HFs undergo periodic growth after birth and are regulated by the complex gene regulatory network. However, the molecular mechanism of HFs regeneration in the Tianzhu white yak remains unclear. RNA editing is a post-transcriptional mechanism that regulates gene expression and produces new transcripts. Hence, we investigated the influence of the A-to-I RNA editing events on the HFs cycle of the Tianzhu white yak.ResultsWe finally identified 54,707 adenosine-to-inosine (A-to-I) RNA editing sites (RESs) from RNA sequencing data of the HFs cycle in the Tianzhu white yak. Annotation results showed RESs caused missense amino acid changes in 7 known genes. And 202 A-to-I editing sites altered 23 target genes of 140 microRNAs. A total of 1,722 differential RESs were identified during the HFs cycle of Tianzhu white yak. GO and KEGG enrichment analysis revealed several signaling pathways and GO terms involved skin development, hair growth, and HFs cycle. Such as genes with differential RNA editing levels were significantly enriched in the peroxisome, metabolic pathways, Notch signaling pathway, and PPAR signaling pathway. Besides, the editing sites in HFs development-related genes FAS, APCDD1, WWOX, MPZL3, RUNX1, KANK2, DCN, DSC2, LEPR, HEPHL1, and PTK2B were suggested as the potential RESs involving HFs development.ConclusionThis study investigated the global A-to-I RNA editing events during the HFs cycle of yak skin tissue and expanded the knowledge of A-to-I RNA editing on the HFs cycle. Furthermore, this study revealed that RNA editing-influenced genes may regulate the HFs cycle by participating in the HFs development-related pathways. The findings might provide new insight into the regulation of RNA editing in hair growth.

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

RNA editing is a post-transcriptional alteration of RNA sequences that, via insertions, deletions or base substitutions, can affect protein structure as well as RNA and protein expression. Recently, it has been suggested that RNA editing may be more frequent than previously thought. A great impediment, however, to a deeper understanding of this process is the paramount sequencing effort that needs to be undertaken to identify RNA editing events. Here, we describe an in silico approach, based on machine learning, that ameliorates this problem. Using 41 nucleotide long DNA sequences, we show that novel A-to-I RNA editing events can be predicted from known A-to-I RNA editing events intra- and interspecies. The validity of the proposed method was verified in an independent experimental dataset. Using our approach, 203 202 putative A-to-I RNA editing events were predicted in the whole human genome. Out of these, 9% were previously reported. The remaining sites require further validation, e.g., by targeted deep sequencing. In conclusion, the approach described here is a useful tool to identify potential A-to-I RNA editing events without the requirement of extensive RNA sequencing.

OP0297 ABERRANT ADENOSINE TO INOSINE RNA EDITING IN ACTIVE RHEUMATOID ARTHRITIS

RNA editing, as an epigenetic mechanism, exhibits a strong correlation with the occurrence and development of cancers. Nevertheless, few studies have been conducted to investigate the impact of RNA editing on cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC). In order to study the connection between RNA editing and CESC patients' prognoses, we obtained CESC-related information from The Cancer Genome Atlas (TCGA) database and randomly allocated the patients into the training group or testing group. An RNA editing-based risk model for CESC patients was established by Cox regression analysis and least absolute shrinkage and selection operator (LASSO). According to the median score generated by this RNA editing-based risk model, patients were categorized into subgroups with high and low risks. We further constructed the nomogram by risk scores and clinical characteristics and analyzed the impact of RNA editing levels on host gene expression levels and adenosine deaminase acting on RNA. Finally, we also compared the biological functions and pathways of differentially expressed genes (DEGs) between different subgroups by enrichment analysis. In this risk model, we screened out 6 RNA editing sites with significant prognostic value. The constructed nomogram performed well in forecasting patients' prognoses. Furthermore, the level of RNA editing at the prognostic site exhibited a strong correlation with host gene expression. In the high-risk subgroup, we observed multiple biological functions and pathways associated with immune response, cell proliferation, and tumor progression. This study establishes an RNA editing-based risk model that helps forecast patients' prognoses and offers a new understanding of the underlying mechanism of RNA editing in CESC.

RNA editing is a crucial post-transcriptional modification process in plant organellar RNA metabolism. rRNA removal-based total RNA-seq is one of the most common methods to study this event. However, the lack of commercial kits to remove rRNAs limits the usage of this method, especially for non-model plant species. DSN-seq is a transcriptome sequencing method utilizing duplex-specific nuclease (DSN) to degrade highly abundant cDNA species especially those from rRNAs while keeping the robustness of transcript levels of the majority of other mRNAs, and has not been applied to study RNA editing in plants before. In this study, we evaluated the capability of DSN-seq to reduce rRNA content and profile organellar RNA editing events in plants, as well we used commercial Ribo-off-seq and standard mRNA-seq as comparisons. Our results demonstrated that DSN-seq efficiently reduced rRNA content and enriched organellar transcriptomes in rice. With high sensitivity to RNA editing events, DSN-seq and Ribo-off-seq provided a more complete and accurate RNA editing profile of rice, which was further validated by Sanger sequencing. Furthermore, DSN-seq also demonstrated efficient organellar transcriptome enrichment and high sensitivity for profiling RNA editing events in Arabidopsis thaliana. Our study highlights the capability of rRNA removal-based total RNA-seq for profiling RNA editing events in plant organellar transcriptomes and also suggests DSN-seq as a widely accessible RNA editing profiling method for various plant species.

BackgroundA-to-I RNA editing is a co−/post-transcriptional modification catalyzed by ADAR enzymes, that deaminates Adenosines (A) into Inosines (I). Most of known editing events are located within inverted ALU repeats, but they also occur in coding sequences and may alter the function of encoded proteins. RNA editing contributes to generate transcriptomic diversity and it is found altered in cancer, autoimmune and neurological disorders. Emerging evidences indicate that editing process could be influenced by genetic variations, biological and environmental variables.ResultsWe analyzed RNA editing levels in human blood using RNA-seq data from 459 healthy individuals and identified 2079 sites consistently edited in this tissue. As expected, analysis of gene expression revealed that ADAR is the major contributor to editing on these sites, explaining ~ 13% of observed variability. After removing ADAR effect, we found significant associations for 1122 genes, mainly involved in RNA processing. These genes were significantly enriched in genes encoding proteins interacting with ADARs, including 276 potential ADARs interactors and 9 ADARs direct partners. In addition, our analysis revealed several factors potentially influencing RNA editing in blood, including cell composition, age, Body Mass Index, smoke and alcohol consumption. Finally, we identified genetic loci associated with editing levels, including known ADAR eQTLs and a small region on chromosome 7, containing LOC730338, a lincRNA gene that appears to modulate ADARs mRNA expression.ConclusionsOur data provides a detailed picture of the most relevant RNA editing events and their variability in human blood, giving interesting insights on potential mechanisms behind this post-transcriptional modification and its regulation in this tissue.

Rett syndrome is a neurological disease due to loss-of-function mutations in the transcription factor, Methyl CpG binding protein 2 (MECP2). Because overexpression of endogenous MECP2 also causes disease, we have exploited a targeted RNA-editing approach to repair patient mutations where levels of MECP2 protein will never exceed endogenous levels. Here, we have constructed adeno-associated viruses coexpressing a bioengineered wild-type ADAR2 catalytic domain (Editasewt) and either Mecp2-targeting or nontargeting gfp RNA guides. The viruses are introduced systemically into male mice containing a guanosine to adenosine mutation that eliminates MeCP2 protein and causes classic Rett syndrome in humans. We find that in the mutant mice injected with the Mecp2-targeting virus, the brainstem exhibits the highest RNA-editing frequency compared to other brain regions. The efficiency is sufficient to rescue MeCP2 expression and function in the brainstem of mice expressing the Mecp2-targeting virus. Correspondingly, we find that abnormal Rett-like respiratory patterns are alleviated, and survival is prolonged, compared to mice injected with the control gfp guide virus. The levels of RNA editing among most brain regions corresponds to the distribution of guide RNA rather than Editasewt. Our results provide evidence that a targeted RNA-editing approach can alleviate a hallmark symptom in a mouse model of human disease.

Genome-wide detection of conserved stage-specific RNA editing events for nuclear genes in the Plasmodium falciparum 3D7 malaria parasite.

Hundreds to millions of adenosine (A)-to-inosine (I) modifications are present in eukaryotic transcriptomes and play an essential role in the creation of proteomic and phenotypic diversity. As adenosine and inosine have different base-pairing properties, the functional consequences of these modifications or 'edits' include altering coding potential, splicing, and miRNA-mediated gene silencing of transcripts. However, rather than serving as a static control of gene expression, A-to-I editing provides a means to dynamically rewire the genetic code during development and in a cell-type specific manner. Interestingly, during normal development, in specific cells, and in both neuropathological diseases and cancers, the extent of RNA editing does not directly correlate with levels of the substrate mRNA or the adenosine deaminase that act on RNA (ADAR) editing enzymes, implying that cellular factors are required for spatiotemporal regulation of A-to-I editing. The factors that affect the specificity and extent of ADAR activity have been thoroughly dissected in vitro. Yet, we still lack a complete understanding of how specific ADAR family members can selectively deaminate certain adenosines while others cannot. Additionally, in the cellular environment, ADAR specificity and editing efficiency is likely to be influenced by cellular factors, which is currently an area of intense investigation. Data from many groups have suggested two main mechanisms for controlling A-to-I editing in the cell: (1) regulating ADAR accessibility to target RNAs and (2) protein-protein interactions that directly alter ADAR enzymatic activity. Recent studies suggest cis- and trans-acting RNA elements, heterodimerization and RNA-binding proteins play important roles in regulating RNA editing levels in vivo. WIREs RNA 2016, 7:113-127. doi: 10.1002/wrna.1319.

BackgroundAs the number of RNA-seq datasets that become available to explore transcriptome diversity increases, so does the need for easy-to-use comprehensive computational workflows. Many available tools facilitate analyses of one of the two major mechanisms of transcriptome diversity, namely, differential expression of isoforms due to alternative splicing, while the second major mechanism—RNA editing due to post-transcriptional changes of individual nucleotides—remains under-appreciated. Both these mechanisms play an essential role in physiological and diseases processes, including cancer and neurological disorders. However, elucidation of RNA editing events at transcriptome-wide level requires increasingly complex computational tools, in turn resulting in a steep entrance barrier for labs who are interested in high-throughput variant calling applications on a large scale but lack the manpower and/or computational expertise.ResultsHere we present an easy-to-use, fully automated, computational pipeline (Automated Isoform Diversity Detector, AIDD) that contains open source tools for various tasks needed to map transcriptome diversity, including RNA editing events. To facilitate reproducibility and avoid system dependencies, the pipeline is contained within a pre-configured VirtualBox environment. The analytical tasks and format conversions are accomplished via a set of automated scripts that enable the user to go from a set of raw data, such as fastq files, to publication-ready results and figures in one step. A publicly available dataset of Zika virus-infected neural progenitor cells is used to illustrate AIDD’s capabilities.ConclusionsAIDD pipeline offers a user-friendly interface for comprehensive and reproducible RNA-seq analyses. Among unique features of AIDD are its ability to infer RNA editing patterns, including ADAR editing, and inclusion of Guttman scale patterns for time series analysis of such editing landscapes. AIDD-based results show importance of diversity of ADAR isoforms, key RNA editing enzymes linked with the innate immune system and viral infections. These findings offer insights into the potential role of ADAR editing dysregulation in the disease mechanisms, including those of congenital Zika syndrome. Because of its automated all-inclusive features, AIDD pipeline enables even a novice user to easily explore common mechanisms of transcriptome diversity, including RNA editing landscapes.

Dihydrofolate reductase (DHFR) plays a key role in folate metabolism and is a target molecule of methotrexate. An increase in the cellular expression level of DHFR is one of the mechanisms of tumor resistance to methotrexate. The present study investigated the possibility that adenosine-to-inosine RNA editing, which causes nucleotide conversion by adenosine deaminase acting on RNA (ADAR) enzymes, might modulate DHFR expression. In human breast adenocarcinoma-derived MCF-7 cells, 26 RNA editing sites were identified in the 3'-UTR of DHFR. Knockdown of ADAR1 decreased the RNA editing levels of DHFR and resulted in a decrease in the DHFR mRNA and protein levels, indicating that ADAR1 up-regulates DHFR expression. Using a computational analysis, miR-25-3p and miR-125a-3p were predicted to bind to the non-edited 3'-UTR of DHFR but not to the edited sequence. The decrease in DHFR expression by the knockdown of ADAR1 was restored by transfection of antisense oligonucleotides for these miRNAs, suggesting that RNA editing mediated up-regulation of DHFR requires the function of these miRNAs. Interestingly, we observed that the knockdown of ADAR1 decreased cell viability and increased the sensitivity of MCF-7 cells to methotrexate. ADAR1 expression levels and the RNA editing levels in the 3'-UTR of DHFR in breast cancer tissues were higher than those in adjacent normal tissues. Collectively, the present study demonstrated that ADAR1 positively regulates the expression of DHFR by editing the miR-25-3p and miR-125a-3p binding sites in the 3'-UTR of DHFR, enhancing cellular proliferation and resistance to methotrexate.