Characterizing of functional human coding RNA editing from evolutionary, structural, and dynamic perspectives

A-to-I RNA editing contributes to the persistence of predicted damaging mutations in populations

Adenosine-to-inosine (A-to-I) RNA editing is a very common co-/posttranscriptional modification that can lead to A-to-G changes at the RNA level and compensate for G-to-A genomic changes to a certain extent. It has been shown that each healthy individual can carry dozens of missense variants predicted to be severely deleterious. Why strongly detrimental variants are preserved in a population and not eliminated by negative natural selection remains mostly unclear. Here, we ask if RNA editing correlates with the burden of deleterious A/G polymorphisms in a population. Integrating genome and transcriptome sequencing data from 447 human lymphoblastoid cell lines, we show that nonsynonymous editing activities (prevalence/level) are negatively correlated with the deleteriousness of A-to-G genomic changes and positively correlated with that of G-to-A genomic changes within the population. We find a significantly negative correlation between nonsynonymous editing activities and allele frequency of A within the population. This negative editing-allele frequency correlation is particularly strong when editing sites are located in highly important genes/loci. Examinations of deleterious missense variants from the 1000 Genomes Project further show a significantly higher proportion of rare missense mutations for G-to-A changes than for other types of changes. The proportion for G-to-A changes increases with increasing deleterious effects of the changes. Moreover, the deleteriousness of G-to-A changes is significantly positively correlated with the percentage of editing enzyme binding motifs at the variants. Overall, we show that nonsynonymous editing is associated with the increased burden of G-to-A missense mutations in healthy individuals, expanding RNA editing in pathogenomics studies.

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.

RNA editing derived epitopes function as cancer antigens to elicit immune responses

In addition to genomic mutations, RNA editing is another major mechanism creating sequence variations in proteins by introducing nucleotide changes in mRNA sequences. Deregulated RNA editing contributes to different types of human diseases, including cancers. Here we report that peptides generated as a consequence of RNA editing are indeed naturally presented by human leukocyte antigen (HLA) molecules. We provide evidence that effector CD8⁺ T cells specific for edited peptides derived from cyclin I are present in human tumours and attack tumour cells that are presenting these epitopes. We show that subpopulations of cancer patients have increased peptide levels and that levels of edited RNA correlate with peptide copy numbers. These findings demonstrate that RNA editing extends the classes of HLA presented self-antigens and that these antigens can be recognised by the immune system.

RNA editing is a biological process that creates sequence variation. Here the authors show that peptides generated as a consequence of RNA editing are naturally presented by human leukocyte antigen (HLA) and serve as antigens to elicit anti-tumour immune responses.

Decreased A-to-I RNA editing as a source of keratinocytes' dsRNA in psoriasis

Recognition of dsRNA molecules activates the MDA5-MAVS pathway and plays a critical role in stimulating type-I interferon responses in psoriasis. However, the source of the dsRNA accumulation in psoriatic keratinocytes remains largely unknown. A-to-I RNA editing is a common co- or post-transcriptional modification that diversifies adenosine in dsRNA, and leads to unwinding of dsRNA structures. Thus, impaired RNA editing activity can result in an increased load of endogenous dsRNAs. Here we provide a transcriptome-wide analysis of RNA editing across dozens of psoriasis patients, and we demonstrate a global editing reduction in psoriatic lesions. In addition to the global alteration, we also detect editing changes in functional recoding sites located in the IGFBP7, COPA, and FLNA genes. Accretion of dsRNA activates autoimmune responses, and therefore the results presented here, linking for the first time an autoimmune disease to reduction in global editing level, are relevant to a wide range of autoimmune diseases.

Genome-wide identification and analysis of A-to-I RNA editing events in bovine by transcriptome sequencing

RNA editing increases the diversity of the transcriptome and proteome. Adenosine-to-inosine (A-to-I) editing is the predominant type of RNA editing in mammals and it is catalyzed by the adenosine deaminases acting on RNA (ADARs) family. Here, we used a largescale computational analysis of transcriptomic data from brain, heart, colon, lung, spleen, kidney, testes, skeletal muscle and liver, from three adult animals in order to identify RNA editing sites in bovine. We developed a computational pipeline and used a rigorous strategy to identify novel editing sites from RNA-Seq data in the absence of corresponding DNA sequence information. Our methods take into account sequencing errors, mapping bias, as well as biological replication to reduce the probability of obtaining a false-positive result. We conducted a detailed characterization of sequence and structural features related to novel candidate sites and found 1,600 novel canonical A-to-I editing sites in the nine bovine tissues analyzed. Results show that these sites 1) occur frequently in clusters and short interspersed nuclear elements (SINE) repeats, 2) have a preference for guanines depletion/enrichment in the flanking 5′/3′ nucleotide, 3) occur less often in coding sequences than other regions of the genome, and 4) have low evolutionary conservation. Further, we found that a positive correlation exists between expression of ADAR family members and tissue-specific RNA editing. Most of the genes with predicted A-to-I editing in each tissue were significantly enriched in biological terms relevant to the function of the corresponding tissue. Lastly, the results highlight the importance of the RNA editome in nervous system regulation. The present study extends the list of RNA editing sites in bovine and provides pipelines that may be used to investigate the editome in other organisms.

Cloning and expression of Malabar grouper (Epinephelus malabaricus) ADAR1 gene in response to immune stimulants and nervous necrosis virus

ADARs are RNA editing catalysts that bind double-stranded RNA and convert adenosine to inosine, a process that can lead to destabilization of dsRNA structures and suppression of mRNA translation. In mammals, ADAR1 genes are involved in various cellular pathways, including interferon (IFN)-mediated response. However, the function of fish ADAR1 remains unclear. We report here the cloning of ADAR1 in Malabar grouper (Epinephelus malabaricus) (MgADAR1) and its response to various immune stimulants. The MgADAR1 cDNA is 5371-bp long, consisting of an open reading frame encoding a putative protein of 1381 amino acids, a 235-nt 5'-terminal untranslated region (UTR), and a 990-nt 3'-UTR. The deduced amino acid sequence exhibits signature features of a chitin synthesis regulation domain, two Z-DNA-binding domains (Z alpha), three dsRNA binding motifs (DSRM) and one tRNA-specific and dsRNA adenosine deaminase domain (ADEAMc). MgADAR1 mRNA expressed ubiquitously in tissues of healthy Malabar grouper, with elevated levels in the brain, gills and eyes. In response to poly (I: C), the MgADAR1 mRNA level was significantly up-regulated in the brain and spleen, but not head kidney. Upon nervous necrosis virus (NNV) infection the level of MgADAR1 increased in the brain, whereas Mx increased in the brain, spleen and head kidney. Induction of MgADAR1 by poly (I: C) and NNV was also observed in vitro. Additionally, the expression of MgADAR1 was upregulated by recombinant grouper IFN in grouper cells. These data indicate an intricate interplay between ADAR1 and NNV infection in grouper as MgADAR1 might be regulated in a tissue-specific manner.

The evolution and adaptation of A-to-I RNA editing

Adenosine-to-inosine (A-to-I) RNA editing is an important post-transcriptional modification that affects the information encoded from DNA to RNA to protein. RNA editing can generate a multitude of transcript isoforms and can potentially be used to optimize protein function in response to varying conditions. In light of this and the fact that millions of editing sites have been identified in many different species, it is interesting to examine the extent to which these sites have evolved to be functionally important. In this review, we discuss results pertaining to the evolution of RNA editing, specifically in humans, cephalopods, and Drosophila. We focus on how comparative genomics approaches have aided in the identification of sites that are likely to be advantageous. The use of RNA editing as a mechanism to adapt to varying environmental conditions will also be reviewed.

Adenosine Deaminases That Act on RNA (ADARs)

Inosine is one of the most common modifications found in human RNAs and the Adenosine Deaminases that act on RNA (ADARs) are the main enzymes responsible for its production. ADARs were first discovered in the 1980s and since then our understanding of ADARs has advanced tremendously. For instance, it is now known that defective ADAR function can cause human diseases. Furthermore, recently solved crystal structures of the human ADAR2 deaminase bound to RNA have provided insights regarding the catalytic and substrate recognition mechanisms. In this chapter, we describe the occurrence of inosine in human RNAs and the newest perspective on the ADAR family of enzymes, including their substrate recognition, catalytic mechanism, regulation as well as the consequences of A-to-I editing, and their relation to human diseases.

Evolutionary analysis reveals regulatory and functional landscape of coding and non-coding RNA editing

Adenosine-to-inosine RNA editing diversifies the transcriptome and promotes functional diversity, particularly in the brain. A plethora of editing sites has been recently identified; however, how they are selected and regulated and which are functionally important are largely unknown. Here we show the cis-regulation and stepwise selection of RNA editing during Drosophila evolution and pinpoint a large number of functional editing sites. We found that the establishment of editing and variation in editing levels across Drosophila species are largely explained and predicted by cis-regulatory elements. Furthermore, editing events that arose early in the species tree tend to be more highly edited in clusters and enriched in slowly-evolved neuronal genes, thus suggesting that the main role of RNA editing is for fine-tuning neurological functions. While nonsynonymous editing events have been long recognized as playing a functional role, in addition to nonsynonymous editing sites, a large fraction of 3’UTR editing sites is evolutionarily constrained, highly edited, and thus likely functional. We find that these 3’UTR editing events can alter mRNA stability and affect miRNA binding and thus highlight the functional roles of noncoding RNA editing. Our work, through evolutionary analyses of RNA editing in Drosophila, uncovers novel insights of RNA editing regulation as well as its functions in both coding and non-coding regions.

Many important modifications are made to RNA to fine-tune genomic information. One type, Adenosine-to-Inosine (A-to-I) RNA editing, changes certain adenosines to inosines and is essential for the neurological well-being of many animals. Although RNA editing occurs at thousands of sites across the genomes of various animals, the functions of nearly all editing events–particularly those in non-coding regions–have not been studied, and what determines whether particular adenosines across the genome are edited has not been fully explored. Here, using the Drosophila genus as model organisms, we analyze the evolution of A-to-I RNA editing to identify a large fraction of both coding and non-coding editing events that are under evolutionary constraint and therefore likely functionally important. We find that non-coding editing events in the 3’UTRs of genes could affect miRNA binding and are associated with a decrease in gene expression levels.

A Common Variant of PROK1 (V67I) Acts as a Genetic Modifier in Early Human Pregnancy through Down-Regulation of Gene Expression

PROK1-V67I has been shown to play a role as a modifier gene in the PROK1-PROKR system of human early pregnancy. To explore the related modifier mechanism of PROK1-V67I, we carried out a comparison study at the gene expression level and the cell function alternation of V67I, and its wild-type (WT), in transiently-transfected cells. We, respectively, performed quantitative RT-PCR and ELISA assays to evaluate the protein and/or transcript level of V67I and WT in HTR-8/SV neo, JAR, Ishikawa, and HEK293 cells. Transiently V67I- or WT-transfected HTR-8/SV neo and HEK293 cells were used to investigate cell function alternations. The transcript and protein expressions were down-regulated in all cell lines, ranging from 20% to 70%, compared with WT. There were no significant differences in the ligand activities of V67I and WT with regard to cell proliferation, cell invasion, calcium influx, and tubal formation. Both PROK1 alleles promoted cell invasion and intracellular calcium mobilization, whereas they had no significant effects on cell proliferation and tubal formation. In conclusion, the biological effects of PROK1-V67I on cell functions are similar to those of WT, and the common variant of V67I may act as a modifier in the PROK1-PROKR system through down-regulation of PROK1 expression. This study may provide a general mechanism that the common variant of V67I, modifying the disease severity of PROK1-related pathophysiologies.

Genome organization and transcriptional regulation of Adenosine Deaminase Acting on RNA gene 1 (ADAR1) in grass carp (Ctenopharyngodon idella)

ADAR1, involved in A-to-I RNA editing, belongs to adenosine deaminase acting on RNA (ADAR) family. A-to-I RNA editing is the most widespread editing phenomenon in higher eukaryotes. In the present study, we cloned and identified the full-length cDNA, complete genomic sequence and the promoter sequence of grass carp (Ctenopharyngodon idella) ADAR1 (CiADAR1) by homology cloning strategy and genome walking. CiADAR1 full-length cDNA is comprised of a 5'UTR (43  bp), a 3'UTR (229 bp) and a 4179 bp ORF encoding a polypeptide of 1392 amino acids. The deduced amino acid sequence of CiADAR1 contains two Z-DNA binding domains, three dsRNA binding motifs and a conserved catalytic domain. The complete genomic CiADAR1 has 16 exons and 15 introns. Phylogenetic tree analysis revealed that CiADAR1 shared high homology with Danio rerio ADAR1 (DrADAR1). RT-PCR showed that CiADAR1 were ubiquitously expressed and significantly up-regulated after stimulation with poly I:C. In spleen and liver, CiADAR1 mRNA reached the peak at 12 h and maintained the highest level during 12-24 h post-injection. CiADAR1 promoter was found to be 1102 bp in length and divided into two distinct regions, the proximal region containing three putative interferon regulatory factor binding elements (IRF-E) and the distal region containing only one IRF-E. To further study the transcriptional regulatory mechanism of CiADAR1, grass carp IRF1 (CiIRF1) and IRF3 (CiIRF3) were expressed in Escherichia coli BL21 and purified by affinity chromatography with the Ni-NTA His-Bind resin. Then, gel mobility shift assay was employed to analyze the affinity of CiADAR1 promoter sequence with CiIRF1 and CiIRF3 in vitro. The result revealed that CiIRF1 and CiIRF3 bound to CiADAR1 promoter with high affinity, indicating that IRF1 and IRF3 could be the potential transcriptional regulatory factor for CiADAR1. Co-transfection of pcDNA3.1-IRF1 (or pcDNA3.1-IRF3) with pGL3-CiADAR1 into C. idella kidney (CIK) cells showed that both IRF1 and IRF3 played a positive role in CiADAR1 transcription. In addition, the mutant assay revealed that the proximal region of CiADAR1 promoter is the main regulatory region in CiADAR1 transcription.

A-to-I RNA Editing: Current Knowledge Sources and Computational Approaches with Special Emphasis on Non-Coding RNA Molecules

RNA editing is a dynamic mechanism for gene regulation attained through the alteration of the sequence of primary RNA transcripts. A-to-I (adenosine-to-inosine) RNA editing, which is catalyzed by members of the adenosine deaminase acting on RNA (ADAR) family of enzymes, is the most common post-transcriptional modification in humans. The ADARs bind double-stranded regions and deaminate adenosine (A) into inosine (I), which in turn is interpreted by the translation and splicing machineries as guanosine (G). In recent years, this modification has been discovered to occur not only in coding RNAs but also in non-coding RNAs (ncRNA), such as microRNAs, small interfering RNAs, transfer RNAs, and long non-coding RNAs. This may have several consequences, such as the creation or disruption of microRNA/mRNA binding sites, and thus affect the biogenesis, stability, and target recognition properties of ncRNAs. The malfunction of the editing machinery is not surprisingly associated with various human diseases, such as neurodegenerative, cardiovascular, and carcinogenic diseases. Despite the enormous efforts made so far, the real biological function of this phenomenon, as well as the features of the ADAR substrate, in particular in non-coding RNAs, has still not been fully understood. In this work, we focus on the current knowledge of RNA editing on ncRNA molecules and provide a few examples of computational approaches to elucidate its biological function.

Knowledge in the Investigation of A-to-I RNA Editing Signals

RNA editing is a post-transcriptional alteration of RNA sequences that is able to affect protein structure as well as RNA and protein expression. Adenosine-to-inosine (A-to-I) RNA editing is the most frequent and common post-transcriptional modification in human, where adenosine (A) deamination produces its conversion into inosine (I), which in turn is interpreted by the translation and splicing machineries as guanosine (G). The disruption of the editing machinery has been associated to various human diseases such as cancer or neurodegenerative diseases. This biological phenomenon is catalyzed by members of the adenosine deaminase acting on RNA (ADAR) family of enzymes and occurs on dsRNA structures. Despite the enormous efforts made in the last decade, the real biological function underlying such a phenomenon, as well as ADAR's substrate features still remain unknown. In this work, we summarize the major computational aspects of predicting and understanding RNA editing events. We also investigate the detection of short motif sequences potentially characterizing RNA editing signals and the use of a logistic regression technique to model a predictor of RNA editing events. The latter, named AIRlINER, an algorithmic approach to assessment of A-to-I RNA editing sites in non-repetitive regions, is available as a web app at: http://alpha.dmi.unict.it/airliner/. Results and comparisons with the existing methods encourage our findings on both aspects.