Human coding RNA editing is generally nonadaptive

Genetic and selective constraints on the optimization of gene product diversity

RNA and protein expressed from the same gene can have diverse isoforms due to various post-transcriptional and post-translational modifications. For the vast majority of alternative isoforms, It is unknown whether they are adaptive or simply biological noise. As we cannot experimentally probe the function of each isoform, we can ask whether the distribution of isoforms across genes and across species is consistent with expectations from different evolutionary processes. However, there is currently no theoretical framework that can generate such predictions. To address this, we developed a mathematical model where isoform abundances are determined collectively by cis -acting loci, trans -acting factors, gene expression levels, and isoform decay rates to predict isoform abundance distributions across species and genes in the face of mutation, genetic drift, and selection. We found that factors beyond selection, such as effective population size and the number of cis -acting loci, significantly influence evolutionary outcomes. Notably, suboptimal phenotypes are more likely to evolve when the population is small and/or when the number of cis -loci is large. We also explored scenarios where modification processes have both beneficial and detrimental effects, revealing a non-monotonic relationship between effective population size and optimization, demonstrating how opposing selection pressures on cis - and trans -acting loci can constrain the optimization of gene product diversity. As a demonstration of the power of our theory, we compared the expected distribution of A-to-I RNA editing levels in coleoids and found this to be largely consistent with non-adaptive explanations.

Learning from the Codon Table: Convergent Recoding Provides Novel Understanding on the Evolution of A-to-I RNA Editing

Adenosine-to-inosine (A-to-I) RNA editing recodes the genetic information. Apart from diversifying the proteome, another tempting advantage of RNA recoding is to correct deleterious DNA mutation and restore ancestral allele. Solid evidences for beneficial restorative editing are very rare in animals. By searching for "convergent recoding" under a phylogenetic context, we proposed this term for judging the potential restorative functions of particular editing site. For the well-known mammalian Gln>Arg (Q>R) recoding site, its ancestral state in vertebrate genomes was the pre-editing Gln, and all 470 available mammalian genomes strictly avoid other three equivalent ways to achieve Arg in protein. The absence of convergent recoding from His>Arg, or synonymous mutations on Gln codons, could be attributed to the strong maintenance on editing motif and structure, but the absence of direct A-to-G mutation is extremely unexpected. With similar ideas, we found cases of convergent recoding in Drosophila genus, reducing the possibility of their restorative function. In summary, we defined an interesting scenario of convergent recoding, the occurrence of which could be used as preliminary judgements for whether a recoding site has a sole restorative role. Our work provides novel insights to the natural selection and evolution of RNA editing.

GTDrift: a resource for exploring the interplay between genetic drift, genomic and transcriptomic characteristics in eukaryotes

We present GTDrift, a comprehensive data resource that enables explorations of genomic and transcriptomic characteristics alongside proxies of the intensity of genetic drift in individual species. This resource encompasses data for 1506 eukaryotic species, including 1413 animals and 93 green plants, and is organized in three components. The first two components contain approximations of the effective population size, which serve as indicators of the extent of random genetic drift within each species. In the first component, we meticulously investigated public databases to assemble data on life history traits such as longevity, adult body length and body mass for a set of 979 species. The second component includes estimations of the ratio between the rate of non-synonymous substitutions and the rate of synonymous substitutions (dN/dS) in protein-coding sequences for 1324 species. This ratio provides an estimate of the efficiency of natural selection in purging deleterious substitutions. Additionally, we present polymorphism-derived N e estimates for 66 species. The third component encompasses various genomic and transcriptomic characteristics. With this component, we aim to facilitate comparative transcriptomics analyses across species, by providing easy-to-use processed data for more than 16 000 RNA-seq samples across 491 species. These data include intron-centered alternative splicing frequencies, gene expression levels and sequencing depth statistics for each species, obtained with a homogeneous analysis protocol. To enable cross-species comparisons, we provide orthology predictions for conserved single-copy genes based on BUSCO gene sets. To illustrate the possible uses of this database, we identify the most frequently used introns for each gene and we assess how the sequencing depth available for each species affects our power to identify major and minor splice variants.

New comparative genomic evidence supporting the proteomic diversification role of A-to-I RNA editing in insects

Adenosine-to-inosine (A-to-I) RNA editing, resembling A-to-G mutation, confers adaptiveness by increasing proteomic diversity in a temporal-spatial manner. This evolutionary theory named "proteomic diversifying hypothesis" has only partially been tested in very few organisms like Drosophila melanogaster, mainly by observing the positive selection on nonsynonymous editing events. To find additional genome-wide evidences supporting this interesting assumption, we retrieved the genomes of four Drosophila species and collected 20 deep-sequenced transcriptomes of different developmental stages and neuron populations of D. melanogaster. We systematically profiled the RNA editomes in these samples and performed meticulous comparative genomic analyses. Further evidences were found to support the diversifying hypothesis. (1) None of the nonsynonymous editing sites in D. melanogaster had ancestral G-alleles, while the silent editing sites had an unignorable fraction of ancestral G-alleles; (2) Only very few nonsynonymous editing sites in D. melanogaster had corresponding G-alleles derived in the genomes of sibling species, and the fraction of such situation was significantly lower than that of silent editing sites; (3) The few nonsynonymous editing with corresponding G-alleles had significantly more variable editing levels (across samples) than other nonsynonymous editing sites in D. melanogaster. The proteomic diversifying nature of RNA editing in Drosophila excludes the restorative role which favors an ancestral G-allele. The few fixed G-alleles in sibling species might facilitate the adaptation to particular environment and the corresponding nonsynonymous editing in D. melanogaster would introduce stronger advantage of flexible proteomic diversification. With multi-Omics data, our study consolidates the nature of evolutionary significance of A-to-I RNA editing sites in model insects.

Random genetic drift sets an upper limit on mRNA splicing accuracy in metazoans

Most eukaryotic genes undergo alternative splicing (AS), but the overall functional significance of this process remains a controversial issue. It has been noticed that the complexity of organisms (assayed by the number of distinct cell types) correlates positively with their genome-wide AS rate. This has been interpreted as evidence that AS plays an important role in adaptive evolution by increasing the functional repertoires of genomes. However, this observation also fits with a totally opposite interpretation: given that 'complex' organisms tend to have small effective population sizes (Ne), they are expected to be more affected by genetic drift, and hence more prone to accumulate deleterious mutations that decrease splicing accuracy. Thus, according to this 'drift barrier' theory, the elevated AS rate in complex organisms might simply result from a higher splicing error rate. To test this hypothesis, we analyzed 3496 transcriptome sequencing samples to quantify AS in 53 metazoan species spanning a wide range of Ne values. Our results show a negative correlation between Ne proxies and the genome-wide AS rates among species, consistent with the drift barrier hypothesis. This pattern is dominated by low abundance isoforms, which represent the vast majority of the splice variant repertoire. We show that these low abundance isoforms are depleted in functional AS events, and most likely correspond to errors. Conversely, the AS rate of abundant isoforms, which are relatively enriched in functional AS events, tends to be lower in more complex species. All these observations are consistent with the hypothesis that variation in AS rates across metazoans reflects the limits set by drift on the capacity of selection to prevent gene expression errors.

Liver cancer-specific mutations in functional domains of ADAR2 lead to the elevation of coding and non-coding RNA editing in multiple tumor-related genes

Mutation is the major cause of phenotypic innovations. Apart from DNA mutations, the alteration on RNA such as the ADAR-mediated A-to-I RNA editing could also shape the phenotype. These two layers of variations have not been systematically combined to study their collective roles in cancers. We collected the high-quality transcriptomes of ten hepatocellular carcinoma (HCC) and the matched control samples. We systematically identified HCC-specific mutations in the exonic regions and profiled the A-to-I RNA editome in each sample. All ten HCC samples had mutations in the CDS of ADAR2 gene (dsRNA-binding domain or catalytic domain). The consequence of these mutations converged to the elevation of ADAR2 efficiency as reflected by the global increase of RNA editing levels in HCC. The up-regulated editing sites (UES) were enriched in the CDS and UTR of oncogenes and tumor suppressor genes (TSG), indicating the possible roles of these target genes in HCC oncogenesis. We present the mutation-ADAR2-UES-oncogene/TSG-HCC axis that explains how mutations at different layers would finally lead to abnormal phenotype. In the light of central dogma, our work provides novel insights into how to fully take advantage of the transcriptome data to decipher the consequence of mutations.

Identification and Interpretation of A-to-I RNA Editing Events in Insect Transcriptomes

Adenosine-to-inosine (A-to-I) RNA editing is the most prevalent RNA modification in the nervous systems of metazoans. To study the biological significance of RNA editing, we first have to accurately identify these editing events from the transcriptome. The genome-wide identification of RNA editing sites remains a challenging task. In this review, we will first introduce the occurrence, regulation, and importance of A-to-I RNA editing and then describe the established bioinformatic procedures and difficulties in the accurate identification of these sit esespecially in small sized non-model insects. In brief, (1) to obtain an accurate profile of RNA editing sites, a transcriptome coupled with the DNA resequencing of a matched sample is favorable; (2) the single-cell sequencing technique is ready to be applied to RNA editing studies, but there are a few limitations to overcome; (3) during mapping and variant calling steps, various issues, like mapping and base quality, soft-clipping, and the positions of mismatches on reads, should be carefully considered; (4) Sanger sequencing of both RNA and the matched DNA is the best verification of RNA editing sites, but other auxiliary evidence, like the nonsynonymous-to-synonymous ratio or the linkage information, is also helpful for judging the reliability of editing sites. We have systematically reviewed the understanding of the biological significance of RNA editing and summarized the methodology for identifying such editing events. We also raised several promising aspects and challenges in this field. With insightful perspectives on both scientific and technical issues, our review will benefit the researchers in the broader RNA editing community.

High-level RNA editing diversifies the coleoid cephalopod brain proteome

Coleoid cephalopods (octopus, squid and cuttlefish) have unusually complex nervous systems. The coleoid nervous system is also the only one currently known to recode the majority of expressed proteins through A-to-I RNA editing. The deamination of adenosine by adenosine deaminase acting on RNA (ADAR) enzymes produces inosine, which is interpreted as guanosine during translation. If this occurs in an open reading frame, which is the case for tens of thousands of editing sites in coleoids, it can recode the encoded protein. Here, we describe recent findings aimed at deciphering the mechanisms underlying high-level recoding and its adaptive potential. We describe the complement of ADAR enzymes in cephalopods, including a recently discovered novel domain in sqADAR1. We further summarize current evidence supporting an adaptive role of high-level RNA recoding in coleoids, and review recent studies showing that a large proportion of recoding sites is temperature-sensitive. Despite these new findings, the mechanisms governing the high level of RNA recoding in coleoid cephalopods remain poorly understood. Recent advances using genome editing in squid may provide useful tools to further study A-to-I RNA editing in these animals.

Genome-Wide Analysis on Driver and Passenger RNA Editing Sites Suggests an Underestimation of Adaptive Signals in Insects

Adenosine-to-inosine (A-to-I) RNA editing leads to a similar effect to A-to-G mutations. RNA editing provides a temporo-spatial flexibility for organisms. Nonsynonymous (Nonsyn) RNA editing in insects is over-represented compared with synonymous (Syn) editing, suggesting adaptive signals of positive selection on Nonsyn editing during evolution. We utilized the brain RNA editome of Drosophila melanogaster to systematically study the LD (r2) between editing sites and infer its impact on the adaptive signals of RNA editing. Pairs of editing sites (PESs) were identified from the transcriptome. For CDS PESs of two consecutive editing sites, their occurrence was significantly biased to type-3 PES (Syn-Nonsyn). The haplotype frequency of type-3 PES exhibited a significantly higher abundance of AG than GA, indicating that the rear Nonsyn site is the driver that promotes the editing of the front Syn site (passenger). The exclusion of passenger Syn sites dramatically amplifies the adaptive signal of Nonsyn RNA editing. Our study for the first time quantitatively demonstrates that the linkage between RNA editing events comes from hitchhiking effects and leads to the underestimation of adaptive signals for Nonsyn editing. Our work provides novel insights for studying the evolutionary significance of RNA editing events.

Genes for highly abundant proteins in Escherichia coli avoid 5' codons that promote ribosomal initiation

In many species highly expressed genes (HEGs) over-employ the synonymous codons that match the more abundant iso-acceptor tRNAs. Bacterial transgene codon randomization experiments report, however, that enrichment with such "translationally optimal" codons has little to no effect on the resultant protein level. By contrast, consistent with the view that ribosomal initiation is rate limiting, synonymous codon usage following the 5' ATG greatly influences protein levels, at least in part by modifying RNA stability. For the design of bacterial transgenes, for simple codon based in silico inference of protein levels and for understanding selection on synonymous mutations, it would be valuable to computationally determine initiation optimality (IO) scores for codons for any given species. One attractive approach is to characterize the 5' codon enrichment of HEGs compared with the most lowly expressed genes, just as translational optimality scores of codons have been similarly defined employing the full gene body. Here we determine the viability of this approach employing a unique opportunity: for Escherichia coli there is both the most extensive protein abundance data for native genes and a unique large-scale transgene codon randomization experiment enabling objective definition of the 5' codons that cause, rather than just correlate with, high protein abundance (that we equate with initiation optimality, broadly defined). Surprisingly, the 5' ends of native genes that specify highly abundant proteins avoid such initiation optimal codons. We find that this is probably owing to conflicting selection pressures particular to native HEGs, including selection favouring low initiation rates, this potentially enabling high efficiency of ribosomal usage and low noise. While the classical HEG enrichment approach does not work, rendering simple prediction of native protein abundance from 5' codon content futile, we report evidence that initiation optimality scores derived from the transgene experiment may hold relevance for in silico transgene design for a broad spectrum of bacteria.

A-to-I RNA editing shows dramatic up-regulation in osteosarcoma and broadly regulates tumor-related genes by altering microRNA target regions

A-to-I RNA editing is a prevalent type of RNA modification in animals. The dysregulation of RNA editing has led to multiple human cancers. However, the role of RNA editing has never been studied in osteosarcoma, a complex bone cancer with unknown molecular basis. We retrieved the RNA-sequencing data from 24 primary osteosarcoma patients and 3 healthy controls. We systematically profiled the RNA editomes in these samples and quantitatively identified reliable differential editing sites (DES) between osteosarcoma and normal samples. RNA editing efficiency is dramatically increased in osteosarcoma, presumably due to the significant up-regulation of editing enzymes ADAR1 and ADAR2. Up-regulated DES in osteosarcoma are enriched in 3'UTRs. Strikingly, such 3'UTR sites are further enriched in microRNA binding regions of gene EMP2 and other oncogenes, abolishing the microRNA suppression on target genes. Accordingly, the expression of these tumor-promoting genes is elevated in osteosarcoma. There might be an RNA editing-dependent pathway leading to osteosarcoma. We expanded our knowledge on the potential roles of RNA editing in oncogenesis. Based on these molecular features, our work is valuable for future prognosis and diagnosis of osteosarcoma.

Temperature-dependent RNA editing in octopus extensively recodes the neural proteome

In poikilotherms, temperature changes challenge the integration of physiological function. Within the complex nervous systems of the behaviorally sophisticated coleoid cephalopods, these problems are substantial. RNA editing by adenosine deamination is a well-positioned mechanism for environmental acclimation. We report that the neural proteome of Octopus bimaculoides undergoes massive reconfigurations via RNA editing following a temperature challenge. Over 13,000 codons are affected, and many alter proteins that are vital for neural processes. For two highly temperature-sensitive examples, recoding tunes protein function. For synaptotagmin, a key component of Ca2+-dependent neurotransmitter release, crystal structures and supporting experiments show that editing alters Ca2+ binding. For kinesin-1, a motor protein driving axonal transport, editing regulates transport velocity down microtubules. Seasonal sampling of wild-caught specimens indicates that temperature-dependent editing occurs in the field as well. These data show that A-to-I editing tunes neurophysiological function in response to temperature in octopus and most likely other coleoids.

Experimental evidence for the functional importance and adaptive advantage of A-to-I RNA editing in fungi

Adenosine-to-inosine (A-to-I) editing is the most prevalent type of RNA editing in animals, and it occurs in fungi specifically during sexual reproduction. However, it is debatable whether A-to-I RNA editing is adaptive. Deciphering the functional importance of individual editing sites is essential for the mechanistic understanding of the adaptive advantages of RNA editing. Here, by performing gene deletion for 17 genes with conserved missense editing (CME) sites and engineering underedited (ue) and overedited (oe) mutants for 10 CME sites using site-specific mutagenesis at the native locus in Fusarium graminearum, we demonstrated that two CME sites in CME5 and CME11 genes are functionally important for sexual reproduction. Although the overedited mutant was normal in sexual reproduction, the underedited mutant of CME5 had severe defects in ascus and ascospore formation like the deletion mutant, suggesting that the CME site of CME5 is co-opted for sexual development. The preediting residue of Cme5 is evolutionarily conserved across diverse classes of Ascomycota, while the postediting one is rarely hardwired into the genome, implying that editing at this site leads to higher fitness than a genomic A-to-G mutation. More importantly, mutants expressing only the underedited or the overedited allele of CME11 are defective in ascosporogenesis, while those expressing both alleles displayed normal phenotypes, indicating that concurrently expressing edited and unedited versions of Cme11 is more advantageous than either. Our study provides convincing experimental evidence for the long-suspected adaptive advantages of RNA editing in fungi and likely in animals.

Evolutionary honing in and mutational replacement: how long-term directed mutational responses to specific environmental pressures are possible

Recent results have shown that the human malaria-resistant hemoglobin S mutation originates de novo more frequently in the gene and in the population where it is of adaptive significance, namely, in the hemoglobin subunit beta gene compared to the nonresistant but otherwise identical 20A[Formula: see text]T mutation in the hemoglobin subunit delta gene, and in sub-Saharan Africans, who have been subject to intense malarial pressure for many generations, compared to northern Europeans, who have not. This finding raises a fundamental challenge to the traditional notion of accidental mutation. Here, we address this finding with the replacement hypothesis, according to which preexisting genetic interactions can lead directly and mechanistically to mutations that simplify and replace them. Thus, an evolutionary process under selection can gradually hone in on interactions of importance for the currently evolving adaptations, from which large-effect mutations follow that are relevant to these adaptations. We exemplify this hypothesis using multiple types of mutation, including gene fusion mutations, gene duplication mutations, A[Formula: see text]G mutations in RNA-edited sites and transcription-associated mutations, and place it in the broader context of a system-level view of mutation origination called interaction-based evolution. Potential consequences include that similarity of mutation pressures may contribute to parallel evolution in genetically related species, that the evolution of genome organization may be driven by mutational mechanisms, that transposable element movements may also be explained by replacement, and that long-term directed mutational responses to specific environmental pressures are possible. Such mutational phenomena need to be further tested by future studies in natural and artificial settings.

On the origin and evolution of RNA editing in metazoans

Extensive adenosine-to-inosine (A-to-I) editing of nuclear-transcribed mRNAs is the hallmark of metazoan transcriptional regulation. Here, by profiling the RNA editomes of 22 species that cover major groups of Holozoa, we provide substantial evidence supporting A-to-I mRNA editing as a regulatory innovation originating in the last common ancestor of extant metazoans. This ancient biochemistry process is preserved in most extant metazoan phyla and primarily targets endogenous double-stranded RNA (dsRNA) formed by evolutionarily young repeats. We also find intermolecular pairing of sense-antisense transcripts as an important mechanism for forming dsRNA substrates for A-to-I editing in some but not all lineages. Likewise, recoding editing is rarely shared across lineages but preferentially targets genes involved in neural and cytoskeleton systems in bilaterians. We conclude that metazoan A-to-I editing might first emerge as a safeguard mechanism against repeat-derived dsRNA and was later co-opted into diverse biological processes due to its mutagenic nature.

Chloroplast C-to-U RNA editing in vascular plants is adaptive due to its restorative effect: testing the restorative hypothesis

The adaptiveness of nonsynonymous RNA editing (recoding) could be conferred by the flexibility of the temporal-spatially controllable proteomic diversity, or by its restorative effect which fixes unfavorable genomic mutations at the RNA level. These two complementary hypotheses, namely, the diversifying hypothesis and the restorative hypothesis, have distinct predictions on the landscape of RNA editing sites. We collected the chloroplast C-to-U RNA editomes of 21 vascular plants (11 angiosperms, four gymnosperms, and six ferns) from a previous study, aiming to testify whether the plant editomes typically conform to the restorative hypothesis. All predictions made by the restorative hypothesis are verified: (i) nonsynonymous editing sites are more frequent and have higher editing levels than synonymous sites; (ii) nonsynonymous editing levels are extremely high and show weak tissue-specificity in plants; (iii) on the inferred genomic sites with recent T-to-C mutations, nonsynonymous sites but not synonymous sites are compensated by C-to-U RNA editing. In conclusion, nonsynonymous C-to-U RNA editing in plants is adaptive due to its restorative effects. The recoding levels are high and are constantly required across the whole plant so that the recoding events could perfectly mimic DNA mutations. The evolutionary significance of plant RNA editing is systematically demonstrated at the genome-wide level.

Gene product diversity: adaptive or not?

One gene does not equal one RNA or protein. The genomic revolution has revealed numerous different RNA and protein molecules that can be produced from one gene, such as circular RNAs generated by back-splicing, proteins with residues mismatching the genomic encoding because of RNA editing, and proteins extended in the C terminus via stop codon readthrough in translation. Are these diverse products results of exquisite gene regulations or imprecise biological processes? While there are cases where the gene product diversity appears beneficial, genome-scale patterns suggest that much of this diversity arises from nonadaptive, molecular errors. This finding has important implications for studying the functions of diverse gene products and for understanding the fundamental properties and evolution of cellular life.

Loss of CaV1.3 RNA editing enhances mouse hippocampal plasticity, learning, and memory

L-type Ca_V1.3 calcium channels are expressed on the dendrites and soma of neurons, and there is a paucity of information about its role in hippocampal plasticity. Here, by genetic targeting to ablate Ca_V1.3 RNA editing, we demonstrate that unedited Ca_V1.3^ΔECS mice exhibited improved learning and enhanced long-term memory, supporting a functional role of RNA editing in behavior. Significantly, the editing paradox that functional recoding of Ca_V1.3 RNA editing sites slows Ca²⁺-dependent inactivation to increase Ca²⁺ influx but reduces channel open probability to decrease Ca²⁺ influx was resolved. Mechanistically, using hippocampal slice recordings, we provide evidence that unedited Ca_V1.3 channels permitted larger Ca²⁺ influx into the hippocampal pyramidal neurons to bolster neuronal excitability, synaptic transmission, late long-term potentiation, and increased dendritic arborization. Of note, RNA editing of the Ca_V1.3 IQ-domain was found to be evolutionarily conserved in mammals, which lends support to the importance of the functional recoding of the Ca_V1.3 channel in brain function.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Interplay between A-to-I Editing and Splicing of RNA: A Potential Point of Application for Cancer Therapy

Adenosine-to-inosine RNA editing is a system of post-transcriptional modification widely distributed in metazoans which is catalyzed by ADAR enzymes and occurs mostly in double-stranded RNA (dsRNA) before splicing. This type of RNA editing changes the genetic code, as inosine generally pairs with cytosine in contrast to adenosine, and this expectably modulates RNA splicing. We review the interconnections between RNA editing and splicing in the context of human cancer. The editing of transcripts may have various effects on splicing, and resultant alternatively spliced isoforms may be either tumor-suppressive or oncogenic. Dysregulated RNA splicing in cancer often causes the release of excess amounts of dsRNA into cytosol, where specific dsRNA sensors provoke antiviral-like responses, including type I interferon signaling. These responses may arrest cell division, causing apoptosis and, externally, stimulate antitumor immunity. Thus, small-molecule spliceosome inhibitors have been shown to facilitate the antiviral-like signaling and are considered to be potential cancer therapies. In turn, a cytoplasmic isoform of ADAR can deaminate dsRNA in cytosol, thereby decreasing its levels and diminishing antitumor innate immunity. We propose that complete or partial inhibition of ADAR may enhance the proapoptotic and cytotoxic effects of splicing inhibitors and that it may be considered a promising addition to cancer therapies targeting RNA splicing.

Landscape of adenosine-to-inosine RNA recoding across human tissues

RNA editing by adenosine deaminases changes the information encoded in the mRNA from its genomic blueprint. Editing of protein-coding sequences can introduce novel, functionally distinct, protein isoforms and diversify the proteome. The functional importance of a few recoding sites has been appreciated for decades. However, systematic methods to uncover these sites perform poorly, and the full repertoire of recoding in human and other mammals is unknown. Here we present a new detection approach, and analyze 9125 GTEx RNA-seq samples, to produce a highly-accurate atlas of 1517 editing sites within the coding region and their editing levels across human tissues. Single-cell RNA-seq data shows protein recoding contributes to the variability across cell subpopulations. Most highly edited sites are evolutionary conserved in non-primate mammals, attesting for adaptation. This comprehensive set can facilitate understanding of the role of recoding in human physiology and diseases.

Gabay et al. provide a highly-accurate atlas of recoding by A-to-I RNA editing in human, profiled across tissues and cell subpopulations. Most highly edited sites are evolutionary conserved in non-primate mammals, attesting for adaptation.

A hierarchy in clusters of cephalopod mRNA editing sites

RNA editing in the form of substituting adenine with inosine (A-to-I editing) is the most frequent type of RNA editing in many metazoan species. In most species, A-to-I editing sites tend to form clusters and editing at clustered sites depends on editing of the adjacent sites. Although functionally important in some specific cases, A-to-I editing usually is rare. The exception occurs in soft-bodied coleoid cephalopods, where tens of thousands of potentially important A-to-I editing sites have been identified, making coleoids an ideal model for studying of properties and evolution of A-to-I editing sites. Here, we apply several diverse techniques to demonstrate a strong tendency of coleoid RNA editing sites to cluster along the transcript. We show that clustering of editing sites and correlated editing substantially contribute to the transcriptome diversity that arises due to extensive RNA editing. Moreover, we identify three distinct types of editing site clusters, varying in size, and describe RNA structural features and mechanisms likely underlying formation of these clusters. In particular, these observations may explain sequence conservation at large distances around editing sites and the observed dependency of editing on mutations in the vicinity of editing sites.

Repeated Winning and Losing Experiences in Chronic Social Conflicts Are Linked to RNA Editing Pattern Difference

Winner-loser effects influence subsequent agonistic interactions between conspecifics. Previous winning experiences could strengthen future aggression and increase the chance of winning the next agonistic interaction, while previous losing experiences could have the opposite effect. Although the role of A-to-I RNA editing has been recently implicated in chronic social defeat stress and aggressive behavior, it remains to be further elucidated in chronic social conflicts in agonistic interactions, especially in the repeated aggression (winners) and repeated defeat (losers) resulted from these conflicts. In the current study, transcriptome-wide A-to-I RNA editing in the dorsal striatum was investigated in a mouse model of chronic social conflicts, and compared between mice repeatedly winning and losing daily agonistic interactions. Our analysis identified 622 A-to-I RNA editing sites in the mouse dorsal striatum, with 23 to be differentially edited in 22 genes, most of which had been previously associated with neurological, psychiatric, or immune disorders. Among these differential RNA editing (DRE) sites four missense variants were observed in neuroligin 2 (Nlgn2), Cdc42 guanine nucleotide exchange factor 9 (Arhgef9) BLCAP apoptosis inducing factor (Blcap), and cytoplasmic FMR1 interacting protein 2 (Cyfip2), as well as two noncoding RNA sites in small nucleolar RNA host gene 11 (Snhg11) and the maternally expressed 3 (Meg3) gene. Moreover, significant changes were observed in gene functions and pathways enriched by genes with A-to-I RNA editing in losers and especially winners compared to controls. Our results demonstrate that repeated winning and losing experiences in chronic social conflicts are linked to A-to-I RNA editing pattern difference, underlining its role in the molecular mechanism of agonistic interactions between conspecifics.

Fast evolution of SARS-CoV-2 driven by deamination systems in hosts

As an RNA virus, the fast evolution of SARS-CoV-2 is driven by the extensive RNA deamination by the host cells.

Adaptive Proteome Diversification by Nonsynonymous A-to-I RNA Editing in Coleoid Cephalopods

RNA editing by the ADAR enzymes converts selected adenosines into inosines, biological mimics for guanosines. By doing so, it alters protein-coding sequences, resulting in novel protein products that diversify the proteome beyond its genomic blueprint. Recoding is exceptionally abundant in the neural tissues of coleoid cephalopods (octopuses, squids, and cuttlefishes), with an over-representation of nonsynonymous edits suggesting positive selection. However, the extent to which proteome diversification by recoding provides an adaptive advantage is not known. It was recently suggested that the role of evolutionarily conserved edits is to compensate for harmful genomic substitutions, and that there is no added value in having an editable codon as compared with a restoration of the preferred genomic allele. Here, we show that this hypothesis fails to explain the evolutionary dynamics of recoding sites in coleoids. Instead, our results indicate that a large fraction of the shared, strongly recoded, sites in coleoids have been selected for proteome diversification, meaning that the fitness of an editable A is higher than an uneditable A or a genomically encoded G.

Mammalian circular RNAs result largely from splicing errors

Ubiquitous in eukaryotes, circular RNAs (circRNAs) comprise a large class of mostly non-coding RNAs produced by back-splicing. Although some circRNAs have demonstrated biochemical activities, whether most circRNAs are functional is unknown. Here, we test the hypothesis that circRNA production primarily results from splicing error and so is deleterious instead of beneficial. In support of the error hypothesis, our analysis of RNA sequencing data from 11 shared tissues of humans, macaques, and mice finds that (1) back-splicing is much rarer than linear-splicing, (2) the rate of back-splicing diminishes with the splicing amount, (3) the overall prevalence of back-splicing in a species declines with its effective population size, and (4) circRNAs are overall evolutionarily unconserved. We estimate that more than 97% of the observed circRNA production is deleterious. We identify a small number of functional circRNA candidates, and the genome-wide trend strongly suggests that circRNAs are largely non-functional products of splicing errors.

A porcine brain-wide RNA editing landscape

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by ADAR enzymes, is an essential post-transcriptional modification. Although hundreds of thousands of RNA editing sites have been reported in mammals, brain-wide analysis of the RNA editing in the mammalian brain remains rare. Here, a genome-wide RNA-editing investigation is performed in 119 samples, representing 30 anatomically defined subregions in the pig brain. We identify a total of 682,037 A-to-I RNA editing sites of which 97% are not identified before. Within the pig brain, cerebellum and olfactory bulb are regions with most edited transcripts. The editing level of sites residing in protein-coding regions are similar across brain regions, whereas region-distinct editing is observed in repetitive sequences. Highly edited conserved recoding events in pig and human brain are found in neurotransmitter receptors, demonstrating the evolutionary importance of RNA editing in neurotransmission functions. Although potential data biases caused by age, sex or health status are not considered, this study provides a rich resource to better understand the evolutionary importance of post-transcriptional RNA editing.

Genome-wide expression changes mediated by A-to-I RNA editing correlate with hepatic oncogenesis

Background

Adenosine-to-inosine (A-to-I) RNA editing is one of the most prevalent RNA modifications in the animal kingdom. Since inosine is recognized as guanosines, the A-to-I process mimics A-to-G DNA mutations but can be controlled in a more flexible manner compared to DNA alterations.

Methods

We parsed the transcriptomes and translatomes of liver cancer and normal tissues from ten patients. We profiled the landscape of the A-to-I RNA editome in these samples and interrogated whether the A-to-I processes participated in the gene expression regulation in oncogenesis.

Results

Globally, editing activity was enhanced in all tumor samples compared to that in normal samples. Accordingly, expression of the gene encoding the RNA editing enzyme ADAR (adenosine deaminase acting on RNA) was elevated. Two intronic self-editing sites in ADAR mRNAs controlled its splicing pattern and may regulate its translation efficiency (TE). Moreover, the expression of oncogenes was generally upregulated in tumors, whereas tumor suppressor genes (TSG) were downregulated, possibly due to alterations to microRNA binding sites or RNA splicing defects caused by A-to-I editing.

Conclusions

A-to-I RNA editing plays a crucial role in the oncogenesis of liver cancer. ADAR regulates its own expression via self-editing, and it also affects global transcriptomes and translatomes involving cancer-related genes by RNA editing and changing their expression patterns.

Evolutionary driving forces of A-to-I editing in metazoans

Adenosine-to-inosine (A-to-I) RNA editing is an evolutionarily conserved mechanism that converts adenosines to inosines in metazoans' transcriptomes. However, the landscapes of editomes have considerably changed during evolution. Here, we review some of our current knowledge of A-to-I editing in the metazoan transcriptomes, focusing on the possible evolutionary driving forces underlying the editing events. First, we review the evolution of ADAR gene family in animals. Then, we summarize the recent advances in characterizing the editomes of various metazoan species. Next, we highlight several factors contributing to the interspecies differences in editomes, including variations in copy number and expression patterns of ADAR genes, the differences in genomic architectures and contents, and the differences in the efficacy of natural selection. After that, we review the possible diversifying and restorative effects of the editing (recoding) events that change the protein sequences. Finally, we discuss the possible convergent evolution of RNA editing in distantly related clades. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Processing > RNA Editing and Modification.

Mammalian Alternative Translation Initiation Is Mostly Nonadaptive

Alternative translation initiation (ATLI) refers to the existence of multiple translation initiation sites per gene and is a widespread phenomenon in eukaryotes. ATLI is commonly assumed to be advantageous through creating proteome diversity or regulating protein synthesis. We here propose an alternative hypothesis that ATLI arises primarily from nonadaptive initiation errors presumably due to the limited ability of ribosomes to distinguish sequence motifs truly signaling translation initiation from similar sequences. Our hypothesis, but not the adaptive hypothesis, predicts a series of global patterns of ATLI, all of which are confirmed at the genomic scale by quantitative translation initiation sequencing in multiple human and mouse cell lines and tissues. Similarly, although many codons differing from AUG by one nucleotide can serve as start codons, our analysis suggests that using non-AUG start codons is mostly disadvantageous. These and other findings strongly suggest that ATLI predominantly results from molecular error, requiring a major revision of our understanding of the precision and regulation of translation initiation.

Proteome Diversification by RNA Editing

RNA editing is an RNA modification that alters the RNA sequence relative to its genomic blueprint. The most common type of RNA editing is A-to-I editing by double-stranded RNA-specific adenosine deaminase (ADAR) enzymes. Editing of a protein-coding region within the RNA molecule may result in non-synonymous substitutions, leading to a modified protein product. These editing sites, also known as "recoding" sites, contribute to the complexity and diversification of the proteome. Recent computational transcriptomic studies have identified thousands of recoding sites in multiple species, many of which are conserved within (but not usually across) lineages and have functional and evolutionary importance. In this chapter we describe the recoding phenomenon across species, consider its potential utility for diversity and adaptation, and discuss its evolution.

Specificity of RNA Folding and Its Association with Evolutionarily Adaptive mRNA Secondary Structures

The secondary structure is a fundamental feature of both noncoding and messenger RNAs. However, our understanding of the secondary structure of mRNA, especially that of the coding regions, remains elusive, likely due to translation and the lack of RNA-binding proteins that sustain the consensus structure, such as those that bind to noncoding RNA. Indeed, mRNA has recently been found to adopt diverse alternative structures, the overall functional significance of which remains untested. We hereby approached this problem by estimating the folding specificity, i.e., the probability that a fragment of RNA folds back to the same partner once refolded. We showed that the folding specificity of mRNA is lower than that of noncoding RNA and exhibits moderate evolutionary conservation. Notably, we found that specific rather than alternative folding is likely evolutionarily adaptive since specific folding is frequently associated with functionally important genes or sites within a gene. Additional analysis in combination with ribosome density suggests the ability to modulate ribosome movement as one potential functional advantage provided by specific folding. Our findings revealed a novel facet of the RNA structurome with important functional and evolutionary implications and indicated a potential method for distinguishing the mRNA secondary structures maintained by natural selection from molecular noise.

A-to-I RNA editing in honeybees shows signals of adaptation and convergent evolution

Social insects exhibit extensive phenotypic diversities among the genetically similar individuals, suggesting a role for the epigenetic regulations beyond the genome level. The ADAR-mediated adenosine-to-inosine (A-to-I) RNA editing, an evolutionarily conserved mechanism, facilitates adaptive evolution by expanding proteomic diversities. Here, we characterize the A-to-I RNA editome of honeybees (Apis mellifera), identifying 407 high-confidence A-to-I editing sites. Editing is most abundant in the heads and shows signatures for positive selection. Editing behavior differs between foragers and nurses, suggesting a role for editing in caste differentiation. Although only five sites are conserved between bees and flies, an unexpectedly large number of genes exhibit editing in both species, albeit at different locations, including the nonsynonymous auto-editing of Adar. This convergent evolution, where the same target genes independently acquire recoding events in distant diverged clades, together with the signals of adaptation observed in honeybees alone, further supports the notion of recoding being adaptive.

•

Nonsynonymous editing sites in honeybees were under positive selection

•

Differential editing may contribute to the phenotypic diversity between sub-castes

•

Target genes acquire editing in different clades, suggesting convergent evolution

Human A-to-I RNA editing SNP loci are enriched in GWAS signals for autoimmune diseases and under balancing selection

BACKGROUND

Adenosine-to-inosine (A-to-I) RNA editing plays important roles in diversifying the transcriptome and preventing MDA5 sensing of endogenous dsRNA as nonself. To date, few studies have investigated the population genomic signatures of A-to-I editing due to the lack of editing sites overlapping with SNPs.

RESULTS

In this study, we applied a pipeline to robustly identify SNP editing sites from population transcriptomic data and combined functional genomics, GWAS, and population genomics approaches to study the function and evolution of A-to-I editing. We find that the G allele, which is equivalent to edited I, is overrepresented in editing SNPs. Functionally, A/G editing SNPs are highly enriched in GWAS signals of autoimmune and immune-related diseases. Evolutionarily, derived allele frequency distributions of A/G editing SNPs for both A and G alleles as the ancestral alleles are skewed toward intermediate frequency alleles relative to neutral SNPs, a hallmark of balancing selection, suggesting that both A and G alleles are functionally important. The signal of balancing selection is confirmed by a number of additional population genomic analyses.

CONCLUSIONS

We uncovered a hidden layer of A-to-I RNA editing SNP loci as a common target of balancing selection, and we propose that the maintenance of such editing SNP variations may be at least partially due to constraints on the resolution of the balance between immune activity and self-tolerance.

Adaptive evolution at mRNA editing sites in soft-bodied cephalopods

BACKGROUND

The bulk of variability in mRNA sequence arises due to mutation-change in DNA sequence which is heritable if it occurs in the germline. However, variation in mRNA can also be achieved by post-transcriptional modification including mRNA editing, changes in mRNA nucleotide sequence that mimic the effect of mutations. Such modifications are not inherited directly; however, as the processes affecting them are encoded in the genome, they have a heritable component, and therefore can be shaped by selection. In soft-bodied cephalopods, adenine-to-inosine RNA editing is very frequent, and much of it occurs at nonsynonymous sites, affecting the sequence of the encoded protein.

METHODS

We study selection regimes at coleoid A-to-I editing sites, estimate the prevalence of positive selection, and analyze interdependencies between the editing level and contextual characteristics of editing site.

RESULTS

Here, we show that mRNA editing of individual nonsynonymous sites in cephalopods originates in evolution through substitutions at regions adjacent to these sites. As such substitutions mimic the effect of the substitution at the edited site itself, we hypothesize that they are favored by selection if the inosine is selectively advantageous to adenine at the edited position. Consistent with this hypothesis, we show that edited adenines are more frequently substituted with guanine, an informational analog of inosine, in the course of evolution than their unedited counterparts, and for heavily edited adenines, these transitions are favored by positive selection. Our study shows that coleoid editing sites may enhance adaptation, which, together with recent observations on Drosophila and human editing sites, points at a general role of RNA editing in the molecular evolution of metazoans.

Portuguese wild grapevine genome re-sequencing (Vitis vinifera sylvestris)

The first genome of Vitis vinifera vinifera (PN40024), published in 2007, boosted grapevine related studies. While this reference genome is a suitable tool for the overall studies in the field, it lacks the ability to unveil changes accumulated during V. v. vinifera domestication. The subspecies V. v. sylvestris preserves wild characteristics, making it a good material to provide insights into V. v. vinifera domestication. The difference in the reproductive strategy between both subspecies is one of the characteristics that set them apart. While V. v. vinifera flowers are hermaphrodite, V. v. sylvestris is mostly dioecious. In this paper, we compare the re-sequencing of the genomes from a male and a female individual of the wild sylvestris, against the reference vinifera genome (PN40024). Variant analysis reveals a low number but with high impact modifications in coding regions, essentially non-synonymous single nucleotide polymorphisms and frame shifts caused by insertions and deletions. The sex-locus was manually inspected, and the results obtained are in line with the most recent works related with wild grapevine sex. In this paper we also describe for the first time RNA editing in transcripts of 14 genes in the sex-determining region, including VviYABBY and VviPLATZ.

ADARs, RNA editing and more in hematological malignancies

Adenosine-to-inosine (A-to-I) editing is the most prevalent type of RNA editing in humans, mediated by the adenosine deaminases acting on RNA (ADARs). Physiologically, these enzymes are present in the nucleus and/or the cytoplasm, where they catalyze the conversion of adenosines (A) to inosines (I) on double-stranded mRNA molecules. Aberrant ADAR-mediated-editing is a prominent feature in a variety of cancers. Importantly, the biological functions of ADARs and its functional implications in hematological malignancies have recently been unraveled. In this review, we will highlight the functions of ADARs and their involvements in cancer, specifically in hematological malignancies. RNA editing-independent function of cellular processes by ADARs and the potential of developing novel therapeutic approaches revolving RNA editing will also be discussed.

Inosine RNA modifications are enriched at the codon wobble position in mouse oocytes and eggs†

Mammalian oocytes and eggs are transcriptionally quiescent and depend on post-transcriptional mechanisms for proper maturation. Post-transcriptional mRNA modifications comprise an important regulatory mechanism that can alter protein and miRNA recognition sites, splicing, stability, secondary structure, and protein coding. We discovered that fully grown mouse germinal vesicle oocytes and metaphase II eggs display abundant inosine mRNA modifications compared to growing oocytes from postnatal day 12 oocytes. These inosines were enriched in mRNA protein coding regions (CDS) and specifically located at the third codon base, or wobble position. Inosines, observed at lower frequencies in CDS of somatic tissues, were similarly enriched at the codon wobble position. In oocytes and eggs, inosine modifications lead primarily to synonymous changes in mRNA transcripts. Inosines may ultimately affect maternal mRNA stability by changing codon usage, thereby altering translational efficiency and translationally coupled mRNA degradation. These important observations advance our understanding of post-transcriptional mechanisms contributing to mammalian oocyte maturation.

Dynamic landscape and evolution of m6A methylation in human

m⁶A is a prevalent internal modification in mRNAs and has been linked to the diverse effects on mRNA fate. To explore the landscape and evolution of human m⁶A, we generated 27 m⁶A methylomes across major adult tissues. These data reveal dynamic m⁶A methylation across tissue types, uncover both broadly or tissue-specifically methylated sites, and identify an unexpected enrichment of m⁶A methylation at non-canonical cleavage sites. A comparison of fetal and adult m⁶A methylomes reveals that m⁶A preferentially occupies CDS regions in fetal tissues. Moreover, the m⁶A sub-motifs vary between fetal and adult tissues or across tissue types. From the evolutionary perspective, we uncover that the selection pressure on m⁶A sites varies and depends on their genic locations. Unexpectedly, we found that ∼40% of the 3′UTR m⁶A sites are under negative selection, which is higher than the evolutionary constraint on miRNA binding sites, and much higher than that on A-to-I RNA modification. Moreover, the recently gained m⁶A sites in human populations are clearly under positive selection and associated with traits or diseases. Our work provides a resource of human m⁶A profile for future studies of m⁶A functions, and suggests a role of m⁶A modification in human evolutionary adaptation and disease susceptibility.

Pan-RNA editing analysis of the bovine genome

RNA editing is an essential process for modifying nucleotides at specific RNA sites during post-transcription in many species. However, its genomic landscape and characters have not been systematically explored in the bovine genome. In the present study, we characterized global RNA editing profiles from 50 samples of cattle and revealed a range of RNA editing profiles in different tissues. Most editing sites were significantly enriched in specific BovB-derived SINEs, especially the dispersed Bov-tAs, which likely forms dsRNA structures similar to the primate-specific Alu elements. Interestingly, ADARB1 (ADAR2) was observed to be predominant in determining global editing in the bovine genome. Common RNA editing sites among similar tissues were associated with tissue-specific biological functions. Taken together, the wide distribution of RNA editing sites and their tissue-specific characters implied the bovine RNA editome should be further explored.

Natural Selection on Exonic SNPs Shapes Allelic Expression Imbalance (AEI) Adaptability in Lung Cancer Progression

Tumors are driven by a sequence of genetic and epigenetic alterations. Previous studies have mostly focused on the roles of somatic mutations in tumorigenesis, but how germline variants act is largely unknown. In this study, we hypothesized that allelic expression imbalance (AEI) participated in the process of germline variants on tumorigenesis. We screened single-nucleotide polymorphisms (SNPs) as representative germline variants. By using 127 patients’ RNA sequencing data from paired lung cancer and adjacent normal tissues from public databases, we analyzed the effects of the functional consequence of SNPs, function and conservativeness on genes with AEI. We found that natural selection can affect AEI. Functional adaptability of genes with a high frequency of AEI and a correlation of the incidence of AEI with conservativeness were observed in both adjacent tissues and tumor tissues. Moreover, we observed a higher incidence of AEI in genes with non-synonymous SNPs than in those with synonymous SNPs. However, we also found that AEI was affected by allele expression noise, especially in tumor tissues, which led to an increased proportion of AEI, weakened the effect of natural selection and eliminated the influence of the functional consequence of SNPs on AEI. We unveiled a previously unknown adaptive regulatory mechanism in which the effect of natural selection on SNPs can be reflected in allelic expression, which provides insight into a better understanding of cancer evolution.

Pig-specific RNA editing during early embryo development revealed by genome-wide comparisons

Posttranscriptional modification of mRNA sequences through RNA editing can increase transcriptome and proteome diversity in eukaryotes. Studies of fetal and adult tissues showed that adenosine‐to‐inosine RNA editing plays a crucial role in early human development, but there is a lack of global understanding of dynamic RNA editing during mammalian early embryonic development. Therefore, here we used RNA sequencing data from human, pig and mouse during early embryonic development to detect edited genes that may regulate stem cell pluripotency. We observed that although most of the RNA editing sites are located in intergenic, intron and UTR, a few editing sites are in coding regions and may result in nonsynonymous amino acid changes. Some editing sites are predicted to change the structure of a protein. We also report that HNF1A, TBX3, ACLY, ECI1 and ERDR1 are related to embryonic development and cell division.

Mitochondrial genome of the nonphotosynthetic mycoheterotrophic plant Hypopitys monotropa, its structure, gene expression and RNA editing

Heterotrophic plants—plants that have lost the ability to photosynthesize—are characterized by a number of changes at all levels of organization. Heterotrophic plants are divided into two large categories—parasitic and mycoheterotrophic (MHT). The question of to what extent such changes are similar in these two categories is still open. The plastid genomes of nonphotosynthetic plants are well characterized, and they exhibit similar patterns of reduction in the two groups. In contrast, little is known about the mitochondrial genomes of MHT plants. We report the structure of the mitochondrial genome of Hypopitys monotropa, a MHT member of Ericaceae, and the expression of its genes. In contrast to its highly reduced plastid genome, the mitochondrial genome of H. monotropa is larger than that of its photosynthetic relative Vaccinium macrocarpon, and its complete size is ~810 Kb. We observed an unusually long repeat-rich structure of the genome that suggests the existence of linear fragments. Despite this unique feature, the gene content of the H. monotropa mitogenome is typical of flowering plants. No acceleration of substitution rates is observed in mitochondrial genes, in contrast to previous observations in parasitic non-photosynthetic plants. Transcriptome sequencing revealed the trans-splicing of several genes and RNA editing in 33 of 38 genes. Notably, we did not find any traces of horizontal gene transfer from fungi, in contrast to plant parasites, which extensively integrate genetic material from their hosts.

A-to-I RNA Editing Uncovers Hidden Signals of Adaptive Genome Evolution in Animals

In animals, the most common type of RNA editing is the deamination of adenosines (A) into inosines (I). Because inosines basepair with cytosines (C), they are interpreted as guanosines (G) by the cellular machinery and genomically encoded G alleles at edited sites mimic the function of edited RNAs. The contribution of this hardwiring effect on genome evolution remains obscure. We looked for population genomics signatures of adaptive evolution associated with A-to-I RNA edited sites in humans and Drosophila melanogaster. We found that single nucleotide polymorphisms at edited sites occur 3 (humans) to 15 times (Drosophila) more often than at unedited sites, the nucleotide G is virtually the unique alternative allele at edited sites and G alleles segregate at higher frequency at edited sites than at unedited sites. Our study reveals that a significant fraction of coding synonymous and nonsynonymous as well as silent and intergenic A-to-I RNA editing sites are likely adaptive in the distantly related human and Drosophila lineages.

The preponderance of nonsynonymous A-to-I RNA editing in coleoids is nonadaptive

A-to-I editing enzymatically converts the base adenosine (A) in RNA molecules to inosine (I), which is recognized as guanine (G) in translation. Exceptionally abundant A-to-I editing was recently discovered in the neural tissues of coleoids (octopuses, squids, and cuttlefishes), with a greater fraction of nonsynonymous sites than synonymous sites subject to high levels of editing. Although this phenomenon is thought to indicate widespread adaptive editing, its potential advantage is unknown. Here we propose an alternative, nonadaptive explanation. Specifically, increasing the cellular editing activity permits some otherwise harmful G-to-A nonsynonymous substitutions, because the As are edited to Is at sufficiently high levels. These high editing levels are constrained upon substitutions, resulting in the predominance of nonsynonymous editing at highly edited sites. Our evidence for this explanation suggests that the prevalent nonsynonymous editing in coleoids is generally nonadaptive, as in species with much lower editing activities.

The neural tissues of coleoids have a greater fraction of nonsynonymous sites than synonymous sites subject to high levels of A-to-I RNA editing, a pattern thought to indicate widespread adaptive editing. Here the authors propose and provide evidence for an alternative, nonadaptive explanation.

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.