Detection of canonical A-to-G editing events at 3' UTRs and microRNA target sites in human lungs using next-generation sequencing

A-Kinase Anchor Protein 1 deficiency causes mitochondrial dysfunction in mouse model of hyperoxia induced acute lung injury

Background: Critically ill patients on supplemental oxygen therapy eventually develop acute lung injury (ALI). Reactive oxygen species (ROS) produced during ALI perturbs the mitochondrial dynamics resulting in cellular damage. Genetic deletion of the mitochondrial A-kinase anchoring protein 1 (Akap1) in mice resulted in mitochondrial damage, Endoplasmic reticulum (ER) stress, increased expression of mitophagy proteins and pro-inflammatory cytokines, exacerbating hyperoxia-induced Acute Lung Injury (HALI).

Objective: Despite a strong causal link between mitochondrial dysfunction and HALI, the mechanisms governing the disease progression at the transcriptome level is unknown.

Methods: In this study, RNA sequencing (RNA-seq) analysis was carried out using the lungs of Akap1 knockout (Akap1−/−) mice exposed to normoxia or 48 h of hyperoxia followed by quantitative real time PCR and Ingenuity pathway analysis (IPA). Western blot analysis assessed mitochondrial dysfunction, OXPHOS complex (I-V), apoptosis and antioxidant proteins. Mitochondrial enzymatic assays was used to measure the aconitase, fumarase, citrate synthase activities in isolated mitochondria from Akap1−/− vs. Wt mice exposed to hyperoxia.

Results: Transcriptome analysis of Akap1−/− exposed to hyperoxia reveals increases in transcripts encoding electron transport chain (ETC) and tricarboxylic acid cycle (TCA) proteins. Ingenuity pathway analysis (IPA) shows enrichment of mitochondrial dysfunction and oxidative phosphorylation in Akap1−/− mice. Loss of AKAP1, coupled with oxidant injury, significantly decreases the activities of TCA enzymes. Mechanistically, a significant loss of dynamin-related protein 1 (Drp1) phosphorylation at the protein kinase A (PKA) site Serine 637 (Ser637), decreases in Akt phosphorylation at Serine 437 (Ser47) and increase in the expression of pro-apoptotic protein Bax indicate mitochondrial dysfunction. Heme oxygenase-1 (HO-1) levels significantly increased in CD68 positive alveolar macrophages in Akap1−/− lungs, suggesting a strong antioxidant response to hyperoxia.

Conclusion: Overall these results suggest that AKAP1 overexpression and modulation of Drp1 phosphorylation at Ser637 is an important therapeutic strategy for acute lung injury.

Dynamic Variations of 3'UTR Length Reprogram the mRNA Regulatory Landscape

This paper concerns 3′-untranslated regions (3′UTRs) of mRNAs, which are non-coding regulatory platforms that control stability, fate and the correct spatiotemporal translation of mRNAs. Many mRNAs have polymorphic 3′UTR regions. Controlling 3′UTR length and sequence facilitates the regulation of the accessibility of functional effectors (RNA binding proteins, miRNAs or other ncRNAs) to 3′UTR functional boxes and motifs and the establishment of different regulatory landscapes for mRNA function. In this context, shortening of 3′UTRs would loosen miRNA or protein-based mechanisms of mRNA degradation, while 3′UTR lengthening would strengthen accessibility to these effectors. Alterations in the mechanisms regulating 3′UTR length would result in widespread deregulation of gene expression that could eventually lead to diseases likely linked to the loss (or acquisition) of specific miRNA binding sites. Here, we will review the mechanisms that control 3′UTR length dynamics and their alterations in human disorders. We will discuss, from a mechanistic point of view centered on the molecular machineries involved, the generation of 3′UTR variability by the use of alternative polyadenylation and cleavage sites, of mutually exclusive terminal alternative exons (exon skipping) as well as by the process of exonization of Alu cassettes to generate new 3′UTRs with differential functional features.

Genome-Wide Characterization of RNA Editing Sites in Primary Gastric Adenocarcinoma through RNA-seq Data Analysis

RNA editing is a posttranscriptional nucleotide modification in humans. Of the various types of RNA editing, the adenosine to inosine substitution is the most widespread in higher eukaryotes, which is mediated by the ADAR family enzymes. Inosine is recognized by the biological machinery as guanosine; therefore, editing could have substantial functional effects throughout the genome. RNA editing could contribute to cancer either by exclusive editing of tumor suppressor/promoting genes or by introducing transcriptomic diversity to promote cancer progression. Here, we provided a comprehensive overview of the RNA editing sites in gastric adenocarcinoma and highlighted some of their possible contributions to gastric cancer. RNA-seq data corresponding to 8 gastric adenocarcinoma and their paired nontumor counterparts were retrieved from the GEO database. After preprocessing and variant calling steps, a stringent filtering pipeline was employed to distinguish potential RNA editing sites from SNPs. The identified potential editing sites were annotated and compared with those in the DARNED database. Totally, 12362 high-confidence adenosine to inosine RNA editing sites were detected across all samples. Of these, 12105 and 257 were known and novel editing events, respectively. These editing sites were unevenly distributed across genomic regions, and nearly half of them were located in 3′UTR. Our results revealed that 4868 editing sites were common in both normal and cancer tissues. From the remaining sites, 3985 and 3509 were exclusive to normal and cancer tissues, respectively. Further analysis revealed a significant number of differentially edited events among these sites, which were located in protein coding genes and microRNAs. Given the distinct pattern of RNA editing in gastric adenocarcinoma and adjacent normal tissue, edited sites have the potential to serve as the diagnostic biomarkers and therapeutic targets in gastric cancer.

RNA editing alterations define manifestation of prion diseases

Prion diseases are fatal neurodegenerative disorders characterized by rapidly progressive dementia. Sporadic Creutzfeldt–Jakob disease (sCJD) is the most prevalent. We report that, specific gene-expression alterations utilizing a reliable in vivo mouse model (tg340-PRNP129MM) with sCJD MM1 subtype, correlate with human disease manifestations in the brain cortex related to disease progression. RNA-editing functions mediated by the APOBEC and ADAR deaminases possibly affecting protein expression necessary for normal brain function, are altered in disease stages. Our data provide powerful evidence, derived from a humanized sCJD mouse model and human autopsy material, discerning the critical role of gene expression and RNA-editing signatures, introducing disease-associated targets that can be extrapolated in other neurodegenerative disorders with common clinical and molecular features.

Prion diseases are fatal neurodegenerative disorders caused by misfolding of the normal prion protein into an infectious cellular pathogen. Clinically characterized by rapidly progressive dementia and accounting for 85% of human prion disease cases, sporadic Creutzfeldt–Jakob disease (sCJD) is the prevalent human prion disease. Although sCJD neuropathological hallmarks are well-known, associated molecular alterations are elusive due to rapid progression and absence of preclinical stages. To investigate transcriptome alterations during disease progression, we utilized tg340-PRNP129MM mice infected with postmortem material from sCJD patients of the most susceptible genotype (MM1 subtype), a sCJD model that faithfully recapitulates the molecular and pathological alterations of the human disease. Here we report that transcriptomic analyses from brain cortex in the context of disease progression, reveal epitranscriptomic alterations (specifically altered RNA edited pathway profiles, eg., ER stress, lysosome) that are characteristic and possibly protective mainly for preclinical and clinical disease stages. Our results implicate regulatory epitranscriptomic mechanisms in prion disease neuropathogenesis, whereby RNA-editing targets in a humanized sCJD mouse model were confirmed in pathological human autopsy material.

Human cancer tissues exhibit reduced A-to-I editing of miRNAs coupled with elevated editing of their targets

A-to-I RNA editing is an important post-transcriptional modification, known to be altered in tumors. It targets dozens of sites within miRNAs, some of which impact miRNA biogenesis and function, as well as many miRNA recognition sites. However, the full extent of the effect of editing on regulation by miRNAs and its behavior in human cancers is still unknown. Here we systematically characterized miRNA editing in 10 593 human samples across 32 cancer types and normal controls. We find that the majority of previously reported sites show little to no evidence for editing in this dataset, compile a list of 58 reliable miRNA editing sites, and study them across normal and cancer samples. Edited miRNA versions tend to suppress expression of known oncogenes, and, consistently, we observe a clear global tendency for hypo-editing in tumors, in strike contrast to the behavior for mRNA editing, allowing an accurate classification of normal/tumor samples based on their miRNA editing profile. In many cancers this profile correlates with patients' survival. Finally, thousands of miRNA binding sites are differentially edited in cancer. Our study thus establishes the important effect of RNA editing on miRNA-regulation in the tumor cell, with prospects for diagnostic and prognostic applications.

ncRNA Editing: Functional Characterization and Computational Resources

Noncoding RNAs (ncRNAs) have received much attention due to their central role in gene expression and translational regulation as well as due to their involvement in several biological processes and disease development. Small noncoding RNAs (sncRNAs), such as microRNAs and piwiRNAs, have been thoroughly investigated and functionally characterized. Long noncoding RNAs (lncRNAs), known to play an important role in chromatin-interacting transcription regulation, posttranscriptional regulation, cell-to-cell signaling, and protein regulation, are also being investigated to further elucidate their functional roles.Next-generation sequencing (NGS) technologies have greatly aided in characterizing the ncRNAome. Moreover, the coupling of NGS technology together with bioinformatics tools has been essential to the genome-wide detection of RNA modifications in ncRNAs. RNA editing, a common human co-transcriptional and posttranscriptional modification, is a dynamic biological phenomenon able to alter the sequence and the structure of primary transcripts (both coding and noncoding RNAs) during the maturation process, consequently influencing the biogenesis, as well as the function, of ncRNAs. In particular, the dysregulation of the RNA editing machineries have been associated with the onset of human diseases.In this chapter we discuss the potential functions of ncRNA editing and describe the knowledge base and bioinformatics resources available to investigate such phenomenon.

Genome-wide analysis of consistently RNA edited sites in human blood reveals interactions with mRNA processing genes and suggests correlations with cell types and biological variables

BACKGROUND

A-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.

RESULTS

We 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.

CONCLUSIONS

Our 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.

Comprehensive analysis of the tumor immune micro-environment in non-small cell lung cancer for efficacy of checkpoint inhibitor

Characterizing the molecular immune subtype and micro-environment of lung cancer is necessary to understand immunogenic interactions between infiltrating immune and stromal cells, and how tumor cells overcome immune checkpoint blockades. This study seeks to identify computational methodologies for subtyping gene expression-based tumor-immune micro-environment interactions, which differentiate non-small cell lung cancer (NSCLC) into immune-defective and immune-competent subtypes. Here, 101 lung squamous cell carcinomas (LUSCs) and 87 lung adenocarcinomas (LUADs) tumor samples have been analyzed. Several micro-environmental factors differentially induce LUAD or LUSC immune subtypes, as well as immune checkpoint expression. In particular, tumor-associated macrophages (TAMs) are key immune cells play a vital role in inflammation and cancer micro-environments of LUSCs; whereas, regulatory B cells are immunosuppressive and tumorigenic in LUADs. Additionally, cytolytic activity upon CD8⁺ T cell activation is decreased by the abundance of B cells and macrophages in immune-competent subtypes. Therefore, identifying immune subtypes in lung cancer and their impact on tumor micro-environment will lead to clinical tools for assessing LUADs and LUSCs in patients, as well as maximize the efficacy of immune checkpoint inhibitors.

Whole Exome and Transcriptome Analyses Integrated with Microenvironmental Immune Signatures of Lung Squamous Cell Carcinoma

The immune microenvironment in lung squamous cell carcinoma (LUSC) is not well understood, with interactions between the host immune system and the tumor, as well as the molecular pathogenesis of LUSC, awaiting better characterization. To date, no molecularly targeted agents have been developed for LUSC treatment. Identification of predictive and prognostic biomarkers for LUSC could help optimize therapy decisions. We sequenced whole exomes and RNA from 101 tumors and matched noncancer control Korean samples. We used the information to predict subtype-specific interactions within the LUSC microenvironment and to connect genomic alterations with immune signatures. Hierarchical clustering based on gene expression and mutational profiling revealed subtypes that were either immune defective or immune competent. We analyzed infiltrating stromal and immune cells to further characterize the tumor microenvironment. Elevated expression of macrophage 2 signature genes in the immune competent subtype confirmed that tumor-associated macrophages (TAM) linked inflammation and mutation-driven cancer. A negative correlation was evident between the immune score and the amount of somatic copy-number variation (SCNV) of immune genes (r = -0.58). The SCNVs showed a potential detrimental effect on immunity in the immune-deficient subtype. Knowledge of the genomic alterations in the tumor microenvironment could be used to guide design of immunotherapy options that are appropriate for patients with certain cancer subtypes. Cancer Immunol Res; 6(7); 848-59. ©2018 AACR.

ADAR Mediated RNA Editing Modulates MicroRNA Targeting in Human Breast Cancer

RNA editing by RNA specific adenosine deaminase acting on RNA (ADAR) is increasingly being found to alter microRNA (miRNA) regulation. Editing of miRNA transcripts can affect their processing, as well as which messenger RNAs (mRNAs) they target. Further, editing of target mRNAs can also affect their complementarity to miRNAs. Notably, ADAR editing is often increased in malignancy with the effect of these RNA changes being largely unclear. In addition, numerous reports have now identified an array of miRNAs that directly contribute to various malignancies although the majority of their targets remain largely undefined. Here we propose that modulating the targets of miRNAs via mRNA editing is a frequent occurrence in cancer and an underappreciated participant in pathology. In order to more accurately characterize the relationship between these two regulatory processes, this study examined RNA editing events within mRNA sequences of two breast cancer cell lines (MCF-7 and MDA-MB-231) and determined whether or not these edits could modulate miRNA associations. Computational analyses of RNA-Seq data from these two cell lines identified over 50,000 recurrent editing sites within human mRNAs, and many of these were located in 3’ untranslated regions (UTRs). When these locations were screened against the list of currently-annotated miRNAs we discovered that editing caused a subset (~9%) to have significant alterations to mRNA complementarity. One miRNA in particular, miR-140–3p, is known to be misexpressed in many breast cancers, and we found that mRNA editing allowed this miRNA to directly target the apoptosis inducing gene DFFA in MCF-7, but not in MDA-MB-231 cells. As these two cell lines are known to have distinct characteristics in terms of morphology, invasiveness and physiological responses, we hypothesized that the differential RNA editing of DFFA in these two cell lines could contribute to their phenotypic differences. Indeed, we confirmed through western blotting that inhibiting miR-140–3p increases expression of the DFFA protein product in MCF-7, but not MDA-MB-231, and further that inhibition of miR-140–3p also increases cellular growth in MCF-7, but not MDA-MB-231. Broadly, these results suggest that the creation of miRNA targets may be an underappreciated function of ADAR and may help further elucidate the role of RNA editing in tumor pathogenicity.

Aberrant overexpression of ADAR1 promotes gastric cancer progression by activating mTOR/p70S6K signaling

ADAR1, one of adenosine deaminases acting on RNA, modulates RNA transcripts through converting adenosine (A) to inosine (I) by deamination. Emerging evidence has implicated that ADAR1 plays an important role in a few of human cancers, however, its expression and physiological significance in gastric cancer remain undefined. In the present study, we demonstrated that ADAR1 was frequently overexpressed in gastric cancer samples by quantitative real-time PCR analysis. In a gastric cancer tissue microarray, ADAR1 staining was closely correlated with tumor stage (P < 0.001) and N classification (P < 0.001). Functional analysis indicated that ADAR1 overexpression promoted cell proliferation and migration in vitro, whereas ADAR1 knockdown resulted in an opposite phenotypes. Furthermore, ADAR1 knockdown also inhibited tumorigenicity and lung metastasis potential of gastric cancer cells in nude mice models. Mechanistically, ADAR1 expression had a significant effect on phosphorylation level of mTOR, p70S kinase, and S6 ribosomal protein, implying its involvement in the regulation of mTOR signaling pathway. We conclude that ADAR1 contributes to gastric cancer development and progression via activating mTOR/p70S6K/S6 ribosomal protein signaling axis. Our findings suggest that ADAR1 may be a valuable biomarker for GC diagnosis and prognosis and may represent a new novel therapeutic opportunities.

An RNA editing/dsRNA binding-independent gene regulatory mechanism of ADARs and its clinical implication in cancer

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by Adenosine DeAminases acting on double-stranded RNA(dsRNA) (ADAR), occurs predominantly in the 3′ untranslated regions (3′UTRs) of spliced mRNA. Here we uncover an unanticipated link between ADARs (ADAR1 and ADAR2) and the expression of target genes undergoing extensive 3′UTR editing. Using METTL7A (Methyltransferase Like 7A), a novel tumor suppressor gene with multiple editing sites at its 3′UTR, we demonstrate that its expression could be repressed by ADARs beyond their RNA editing and double-stranded RNA (dsRNA) binding functions. ADARs interact with Dicer to augment the processing of pre-miR-27a to mature miR-27a. Consequently, mature miR-27a targets the METTL7A 3′UTR to repress its expression level. In sum, our study unveils that the extensive 3′UTR editing of METTL7A is merely a footprint of ADAR binding, and there are a subset of target genes that are equivalently regulated by ADAR1 and ADAR2 through their non-canonical RNA editing and dsRNA binding-independent functions, albeit maybe less common. The functional significance of ADARs is much more diverse than previously appreciated and this gene regulatory function of ADARs is most likely to be of high biological importance beyond the best-studied editing function. This non-editing side of ADARs opens another door to target cancer.

ADAR2 regulates RNA stability by modifying access of decay-promoting RNA-binding proteins

Adenosine deaminases acting on RNA (ADARs) catalyze the editing of adenosine residues to inosine (A-to-I) within RNA sequences, mostly in the introns and UTRs (un-translated regions). The significance of editing within non-coding regions of RNA is poorly understood. Here, we demonstrate that association of ADAR2 with RNA stabilizes a subset of transcripts. ADAR2 interacts with and edits the 3΄UTR of nuclear-retained Cat2 transcribed nuclear RNA (Ctn RNA). In absence of ADAR2, the abundance and half-life of Ctn RNA are significantly reduced. Furthermore, ADAR2-mediated stabilization of Ctn RNA occurred in an editing-independent manner. Unedited Ctn RNA shows enhanced interaction with the RNA-binding proteins HuR and PARN [Poly(A) specific ribonuclease deadenylase]. HuR and PARN destabilize Ctn RNA in absence of ADAR2, indicating that ADAR2 stabilizes Ctn RNA by antagonizing its degradation by PARN and HuR. Transcriptomic analysis identified other RNAs that are regulated by a similar mechanism. In summary, we identify a regulatory mechanism whereby ADAR2 enhances target RNA stability by limiting the interaction of RNA-destabilizing proteins with their cognate substrates.