Data from Single-Cell Analysis in Lung Adenocarcinoma Implicates RNA Editing in Cancer Innate Immunity and Patient Prognosis

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

OP0297 ABERRANT ADENOSINE TO INOSINE RNA EDITING IN ACTIVE RHEUMATOID ARTHRITIS

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. See related commentary by Luo and Liang, p. 351.

Background: RNA modification plays an important role in many diseases. A comprehensive study of tumor microenvironment (TME) characteristics mediated by RNA modification regulators will improve the understanding of TME immune regulation.Methods: We selected 26 RNA modification “writers” of lung adenocarcinoma (LUAD) samples and performed unsupervised clustering analysis to explore RNA modification patterns in LUAD. Differentially expressed genes (DEGs) with RNA modification patterns were screened to develop a “writers” of RNA modification score (WM score) system. The infiltration ratio of TME cell subsets was analyzed by CIBERSORT.Results: We identified two RNA modification modes showing different characteristics of overall survival (OS) and TME cell infiltration. According to WM score, LUAD patients were divided into a high-WM score group and a low-WM score group. High-scored patients had a poor prognosis and higher tumor mutation burden (TMB), they were more sensitive to four LUAD therapies (erlotinib, XA V939, gefitinib, and KU-55933) and more clinically responsive to PD-L1 treatment. Those with a low WM score showed higher stromal scores, ESTIMATE scores, and survival chance.Conclusion: Our work revealed the potential role of RNA modification patterns in TME, genetic variation, targeted inhibitor therapy, and immunotherapy. Identifying RNA modification pattern of LUAD patients help understand the characteristics of TME and may promote the development of immunotherapy strategies.

Many Zika virus (ZIKV) pathogenesis-related studies have focused primarily on virus-driven pathology and neurotoxicity, instead of considering the possibility of pathogenesis as an (unintended) consequence of host innate immunity: specifically, as the side-effect of an otherwise well-functioning machine. The hypothesis presented here suggests a new way of thinking about the role of host immune mechanisms in disease pathogenesis, focusing on dysregulation of post-transcriptional RNA editing as a candidate driver of a broad range of observed neurodevelopmental defects and neurodegenerative clinical symptoms in both infants and adults linked with ZIKV infections. We collect and synthesize existing evidence of ZIKV-mediated changes in expression of adenosine deaminases that act on RNA (ADARs), known links between abnormal RNA editing and pathogenesis, as well as ideas for potential translational applications, including genomic profile-based molecular diagnostic tools and/or treatment strategies.

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

Background: Ulcerative colitis (UC) is a chronic inflammatory disease continuous from the rectum. We previously reported that chronic inflammation may induce the RNA-editing enzyme ADAR1 and promote carcinogenesis (Journal of Crohn's and Colitis, 2023). However, the relationship between RNA editing and carcinogenesis in UC remains unresolved because we could not identify specific RNA editing sites. Therefore, in this study, we performed a new comprehensive RNA analysis using RNA sequencing to clarify the role of RNA editing in the carcinogenesis process in UC and to construct a prediction model of carcinogenic risk. Methods: We performed immunostaining of ADAR1, an RNA-editing enzyme, in 139 cases of ulcerative colitis. In addition, normal colorectal tissue, UC inflammatory phase, normal rectal mucosa of UC carcinoma cases, and UC carcinoma sites were submitted to RNA-Seq testing in 6 cases each (24 cases total). Bioinformatics analysis was performed to compare RNA editing sites in each group. Results: RNA-seq analysis showed a significant increase in global RNA editing sites in the inflammatory phase compared to normal colon tissue (p = 0.02). This phenomenon was due to increased ADAR1 expression associated with type 1 interferon. We hypothesized that global RNA editing was decreased in normal colonic mucosa during the remission phase but that some sites of prolonged RNA editing might promote carcinogenesis. We attempted to identify carcinogenesis-specific RNA editing sites. We identified 72 cancer-specific RNA editing sites, mainly in long noncoding RNAs such as NEAT1 and in the 3'UTR region of mRNAs such as CTSS. Assuming that the accumulation of these RNA editing leads to carcinogenesis, pathway analysis revealed that pathways promoting cancer invasion and stemness are activated. When ROC curves were generated using these editing sites as markers, the AUC was 0.935, suggesting that analysis of RNA editing in the colonic mucosa of UC patients may be able to predict the risk of carcinogenesis. Conclusions: We identified potential editing sites for carcinogenesis in UC, allowing the identification of carcinogenic mechanisms in UC and the establishment of predictive markers for carcinogenic risk. Citation Format: Toshiaki Takahashi, Kunitoshi Shigeyasu, Yoshitaka Kondo, Hibiki Umeda, Kazuya Moriwake, Masashi Kayano, Kaori Nitta, Fuminori Teraishi, Yuzo Umeda, Hiroshi Tazawa, Toshiyoshi Fujiwara. Establishment of predictive biomarkers for carcinogenesis in ulcerative colitis using RNA editing analyses [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 2 (Late-Breaking, Clinical Trial, and Invited Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(7_Suppl):Abstract nr LB360.

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.

Impairment of RNA editing at a handful of coding sites causes severe disorders, prompting the view that coding RNA editing is highly advantageous. Recent genomic studies have expanded the list of human coding RNA editing sites by more than 100 times, raising the question of how common advantageous RNA editing is. Analyzing 1,783 human coding A-to-G editing sites, we show that both the frequency and level of RNA editing decrease as the importance of a site or gene increases; that during evolution, edited As are more likely than unedited As to be replaced with Gs but not with Ts or Cs; and that among nonsynonymously edited As, those that are evolutionarily least conserved exhibit the highest editing levels. These and other observations reveal the overall nonadaptive nature of coding RNA editing, despite the presence of a few sites in which editing is clearly beneficial. We propose that most observed coding RNA editing results from tolerable promiscuous targeting by RNA editing enzymes, the original physiological functions of which remain elusive.

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.

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by adenosine deaminases acting on RNA (ADARs), is a widespread modification in metazoans. Cumulative evidence has revealed the altered A-to-I editing profiles in cancers, but the underlying mechanism remains unclear. Here, we discover the well-known histone lysine methyltransferase enhancer of zeste homologue 2 (EZH2) as an unexplored ADAR interactor and editing regulator in prostate cancer (PCa). Through competing with interleukin enhancer binding factor 2 (ILF2) for ADAR1 binding, EZH2 reshapes the substrate selectivity of ADAR1 and thus exhibits a bidirectional role in editing regulation. Moreover, EZH2 depletion induces the translational repression of transportin-1 (TRN1), which further results in the accumulation of cytoplasmic ADAR1p110 isoform to protect many oncogenic transcripts from degradation. Consistently, depletion of ADAR1 dramatically enhances the sensitivity of cancer cells and tumors to EZH2 selective degraders. Collectively, our study sheds new light on a link between two layers of epigenetic regulations at histone modification and RNA editing levels, demonstrates a previously uncharacterized role of EZH2 in RNA editing and mRNA stability independently of its lysine methyltransferase activity, and reveals the significance of EZH2-ADAR1 cascade in governing RNA editing and mRNA stability, which may provide additional perspectives for the advancement of EZH2-targeting cancer therapies.

Activation of AZIN1 RNA editing is a novel mechanism that promotes invasive potential of cancer-associated fibroblasts in colorectal cancer

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.

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.

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