Molecular diversity through RNA editing: a balancing act

Mechanism of Nucleic Acid Phosphodiester Bond Cleavage by Human Endonuclease V: MD and QM/MM Calculations Reveal a Versatile Metal Dependence

Human endonuclease V (EndoV) catalytically removes deaminated nucleobases by cleaving the phosphodiester bond as part of RNA metabolism. Despite being implicated in several diseases (cancers, cardiovascular diseases, and neurological disorders) and potentially being a useful tool in biotechnology, details of the human EndoV catalytic pathway remain unclear due to limited experimental information beyond a crystal structure of the apoenzyme and select mutational data. Since a mechanistic understanding is critical for further deciphering the central roles and expanding applications of human EndoV in medicine and biotechnology, molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations were used to unveil the atomistic details of the catalytic pathway. Due to controversies surrounding the number of metals required for nuclease activity, enzyme-substrate models with different numbers of active site metals and various metal-substrate binding configurations were built based on structural data for other nucleases. Subsequent MD simulations revealed the structure and stability of the human EndoV-substrate complex for a range of active site metal binding architectures. Four unique pathways were then characterized using QM/MM that vary in metal number (one versus two) and modes of substrate coordination [direct versus indirect (water-mediated)], with several mechanisms being fully consistent with experimental structural, kinetic, and mutational data for related nucleases, including members of the EndoV family. Beyond uncovering key roles for several active site amino acids (D240 and K155), our calculations highlight that while one metal is essential for human EndoV activity, the enzyme can benefit from using two metals due to the presence of two suitable metal binding sites. By directly comparing one- versus two-metal-mediated P-O bond cleavage reactions within the confines of the same active site, our work brings a fresh perspective to the "number of metals" controversy.

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.

Development and assessment of an RNA editing-based risk model for the prognosis of cervical cancer patients

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

RNA modifications in cellular metabolism: implications for metabolism-targeted therapy and immunotherapy

Cellular metabolism is an intricate network satisfying bioenergetic and biosynthesis requirements of cells. Relevant studies have been constantly making inroads in our understanding of pathophysiology, and inspiring development of therapeutics. As a crucial component of epigenetics at post-transcription level, RNA modification significantly determines RNA fates, further affecting various biological processes and cellular phenotypes. To be noted, immunometabolism defines the metabolic alterations occur on immune cells in different stages and immunological contexts. In this review, we characterize the distribution features, modifying mechanisms and biological functions of 8 RNA modifications, including N6-methyladenosine (m6A), N6,2'-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), 5-methylcytosine (m5C), N4-acetylcytosine (ac4C), N7-methylguanosine (m7G), Pseudouridine (Ψ), adenosine-to-inosine (A-to-I) editing, which are relatively the most studied types. Then regulatory roles of these RNA modification on metabolism in diverse health and disease contexts are comprehensively described, categorized as glucose, lipid, amino acid, and mitochondrial metabolism. And we highlight the regulation of RNA modifications on immunometabolism, further influencing immune responses. Above all, we provide a thorough discussion about clinical implications of RNA modification in metabolism-targeted therapy and immunotherapy, progression of RNA modification-targeted agents, and its potential in RNA-targeted therapeutics. Eventually, we give legitimate perspectives for future researches in this field from methodological requirements, mechanistic insights, to therapeutic applications.

Is Metal Stabilization of the Leaving Group Required or Can Lysine Facilitate Phosphodiester Bond Cleavage in Nucleic Acids? A Computational Study of EndoV

Endonuclease V (EndoV) is a single-metal-dependent enzyme that repairs deaminated DNA nucleobases in cells by cleaving the phosphodiester bond, and this enzyme has proven to be a powerful tool in biotechnology and medicine. The catalytic mechanism used by EndoV must be understood to design new disease detection and therapeutic solutions and further exploit the enzyme in interdisciplinary applications. This study has used a mixed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) approach to compare eight distinct catalytic pathways and provides the first proposed mechanism for bacterial EndoV. The calculations demonstrate that mechanisms involving either direct or indirect metal coordination to the leaving group of the substrate previously proposed for other nucleases are unlikely for EndoV, regardless of the general base (histidine, aspartate, and substrate phosphate moiety). Instead, distinct catalytic pathways are characterized for EndoV that involve K139 stabilizing the leaving group, a metal-coordinated water stabilizing the transition structure, and either H214 or a substrate phosphate group activating the water nucleophile. In silico K139A and H214A mutational results support the newly proposed roles of these residues. Although this is a previously unseen combination of general base, general acid, and metal-binding architecture for a one-metal-dependent endonuclease, our proposed catalytic mechanisms are fully consistent with experimental kinetic, structural, and mutational data. In addition to substantiating a growing body of literature, suggesting that one metal is enough to catalyze P-O bond cleavage in nucleic acids, this new fundamental understanding of the catalytic function will promote the exploration of new and improved applications of EndoV.

REDH: A database of RNA editome in hematopoietic differentiation and malignancy

BACKGROUND

The conversion of adenosine (A) to inosine (I) through deamination is the prevailing form of RNA editing, impacting numerous nuclear and cytoplasmic transcripts across various eukaryotic species. Millions of high-confidence RNA editing sites have been identified and integrated into various RNA databases, providing a convenient platform for the rapid identification of key drivers of cancer and potential therapeutic targets. However, the available database for integration of RNA editing in hematopoietic cells and hematopoietic malignancies is still lacking.

METHODS

We downloaded RNA sequencing (RNA-seq) data of 29 leukemia patients and 19 healthy donors from National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database, and RNA-seq data of 12 mouse hematopoietic cell populations obtained from our previous research were also used. We performed sequence alignment, identified RNA editing sites, and obtained characteristic editing sites related to normal hematopoietic development and abnormal editing sites associated with hematologic diseases.

RESULTS

We established a new database, "REDH", represents RNA editome in hematopoietic differentiation and malignancy. REDH is a curated database of associations between RNA editome and hematopoiesis. REDH integrates 30,796 editing sites from 12 murine adult hematopoietic cell populations and systematically characterizes more than 400,000 edited events in malignant hematopoietic samples from 48 cohorts (human). Through the Differentiation, Disease, Enrichment, and knowledge modules, each A-to-I editing site is systematically integrated, including its distribution throughout the genome, its clinical information (human sample), and functional editing sites under physiological and pathological conditions. Furthermore, REDH compares the similarities and differences of editing sites between different hematologic malignancies and healthy control.

CONCLUSIONS

REDH is accessible at http://www.redhdatabase.com/ . This user-friendly database would aid in understanding the mechanisms of RNA editing in hematopoietic differentiation and malignancies. It provides a set of data related to the maintenance of hematopoietic homeostasis and identifying potential therapeutic targets in malignancies.

Characterization and analysis of multi-organ full-length transcriptomes in Sphaeropteris brunoniana and Alsophila latebrosa highlight secondary metabolism and chloroplast RNA editing pattern of tree ferns

BACKGROUND

Sphaeropteris brunoniana and Alsophila latebrosa are both old relict and rare tree ferns, which have experienced the constant changes of climate and environment. However, little is known about their high-quality genetic information and related research on environmental adaptation mechanisms of them. In this study, combined with PacBio and Illumina platforms, transcriptomic analysis was conducted on the roots, rachis, and pinna of S. brunoniana and A. latebrosa to identify genes and pathways involved in environmental adaptation. Additionally, based on the transcriptomic data of tree ferns, chloroplast genes were mined to analyze their gene expression levels and RNA editing events.

RESULTS

In the study, we obtained 11,625, 14,391 and 10,099 unigenes of S. brunoniana root, rachis, and pinna, respectively. Similarly, a total of 13,028, 11,431 and 12,144 unigenes were obtained of A. latebrosa root, rachis, and pinna, respectively. According to the enrichment results of differentially expressed genes, a large number of differentially expressed genes were enriched in photosynthesis and secondary metabolic pathways of S. brunoniana and A. latebrosa. Based on gene annotation results and phenylpropanoid synthesis pathways, two lignin synthesis pathways (H-lignin and G-lignin) were characterized of S. brunoniana. Among secondary metabolic pathways of A. latebrosa, three types of WRKY transcription factors were identified. Additionally, based on transcriptome data obtained in this study, reported transcriptome data, and laboratory available transcriptome data, positive selection sites were identified from 18 chloroplast protein-coding genes of four tree ferns. Among them, RNA editing was found in positive selection sites of four tree ferns. RNA editing affected the protein secondary structure of the rbcL gene. Furthermore, the expression level of chloroplast genes indicated high expression of genes related to the chloroplast photosynthetic system in all four species.

CONCLUSIONS

Overall, this work provides a comprehensive transcriptome resource of S. brunoniana and A. latebrosa, laying the foundation for future tree fern research.

Withania somnifera mitochondrial atp4 gene editing alters the ATP synthase b subunit, independent of salt stress

Numerous studies have shown that stress in plant cells and organelles with transport electron chains is related to RNA editing. The ATP synthase complex present in mitochondria plays a crucial role in cellular respiration and consists of several subunits. Among them is the b subunit, which is encoded by the mitochondrial atp4 gene. Computing-based analysis of the effects of RNA editing of the Withania somnifera atp4 gene in mitochondria leading to alterations in the b subunit of ATP synthase. Using the CLC Genomic Workbench 3, RNA editing analysis between the control and salt stress conditions was not significantly different. Depending on RNA editing, the tertiary structure model revealed a change in the states of the b subunit, reflecting differences in the central stalk and F1-catalytic domain. The study found that polar edits in the N-terminus of the b subunit allow for efficient H + ion selectivity and introduce a new coiled-coil alpha-helical structure that may help stabilize the complex. The most noteworthy finding of this study was the strong impact of these editing events on the tertiary structure of the b subunit, which has the potential to affect the ATPase activity and indicate that the editing in this subunit aimed to restore the original active protein and not as a response to salt stress.

The Genome-Wide Characterization of Alternative Splicing and RNA Editing in the Development of Coprinopsis cinerea

Coprinopsis cinerea is one of the model species used in fungal developmental studies. This mushroom-forming Basidiomycetes fungus has several developmental destinies in response to changing environments, with dynamic developmental regulations of the organism. Although the gene expression in C. cinerea development has already been profiled broadly, previous studies have only focused on a specific stage or process of fungal development. A comprehensive perspective across different developmental paths is lacking, and a global view on the dynamic transcriptional regulations in the life cycle and the developmental paths is far from complete. In addition, knowledge on co- and post-transcriptional modifications in this fungus remains rare. In this study, we investigated the transcriptional changes and modifications in C. cinerea during the processes of spore germination, vegetative growth, oidiation, sclerotia formation, and fruiting body formation by inducing different developmental paths of the organism and profiling the transcriptomes using the high-throughput sequencing method. Transition in the identity and abundance of expressed genes drive the physiological and morphological alterations of the organism, including metabolism and multicellularity construction. Moreover, stage- and tissue-specific alternative splicing and RNA editing took place and functioned in C. cinerea. These modifications were negatively correlated to the conservation features of genes and could provide extra plasticity to the transcriptome during fungal development. We suggest that C. cinerea applies different molecular strategies in its developmental regulation, including shifts in expressed gene sets, diversifications of genetic information, and reversible diversifications of RNA molecules. Such features would increase the fungal adaptability in the rapidly changing environment, especially in the transition of developmental programs and the maintenance and balance of genetic and transcriptomic divergence. The multi-layer regulatory network of gene expression serves as the molecular basis of the functioning of developmental regulation.

The cellular and KSHV A-to-I RNA editome in primary effusion lymphoma and its role in the viral lifecycle

Adenosine-to-inosine RNA editing is a major contributor to transcriptome diversity in animals with far-reaching biological consequences. Kaposi's sarcoma-associated herpesvirus (KSHV) is the etiological agent of several human malignancies including primary effusion lymphoma (PEL). The extent of RNA editing within the KSHV transcriptome is unclear as is its contribution to the viral lifecycle. Here, we leverage a combination of biochemical and genomic approaches to determine the RNA editing landscape in host- and KSHV transcriptomes during both latent and lytic replication in PEL. Analysis of RNA editomes reveals it is dynamic, with increased editing upon reactivation and the potential to deregulate pathways critical for latency and tumorigenesis. In addition, we identify conserved RNA editing events within a viral microRNA and discover their role in miRNA biogenesis as well as viral infection. Together, these results describe the editome of PEL cells as well as a critical role for A-to-I editing in the KSHV lifecycle.

Characteristics of Adenosine-to-Inosine RNA editing-based subtypes and novel risk score for the prognosis and drug sensitivity in stomach adenocarcinoma

Stomach adenocarcinoma (STAD) is always characterized by high mortality and poor prognosis with drug resistance and recrudescence due to individual genetic heterogeneity. Adenosine-to-Inosine RNA editing (ATIRE) has been reported associated with multiple tumors but the potential connection between ATIRE-related signatures and STAD remains unclear. In this study, we comprehensively elevated the genetic characteristics of ATIRE in STAD patients and first screened five vital survival-related ATIRE sites to identify a novel ATIRE-Risk score. Based on the risk scores, we further divided the patients into two different subtypes with diverse clinical characteristics and immune landscapes including immune cell infiltration (ICI), tumor microenvironment (TME), and immune checkpoint expression analysis. The low-risk subgroups, associated with better survival prognosis, were characterized by activated immune-cells, higher immune scores in TME, and down-expression of immunotherapy checkpoints. Moreover, different expressional genes (DEGs) between the above subtypes were further identified and the activation of immune-related pathways were found in low-risk patients. The stratified survival analysis further indicated patients with low-risk and high-tumor mutation burden (TMB) exhibited the best prognosis outcomes, implying the role of TMB and ATIRE-Risk scores was synergistic for the prognosis of STAD. Interestingly, anti-tumor chemotherapeutic drugs all exhibited lower IC50 values in low-risk subgroups, suggesting these patients might obtain a better curative response from the combined chemotherapy of STAD. Finally, combined with classical clinical features and ATIRE-Risk scores, we successfully established a promising nomogram system to accurately predict the 1/3/5-years survival ratio of STAD and this model was also estimated with high diagnostic efficiency and stable C-index with calibration curves. These significant ATIRE sites are promising to be further explored and might serve as a novel therapeutic target for STAD treatment.

RNA Editing Alterations Define Disease Manifestations in the Progression of Experimental Autoimmune Encephalomyelitis (EAE)

RNA editing is an epitranscriptomic modification, leading to targeted changes in RNA transcripts. It is mediated by the action of ADAR (adenosine deaminases acting on double-stranded (ds) RNA and APOBEC (apolipoprotein B mRNA editing enzyme catalytic polypeptide-like) deaminases and appears to play a major role in the pathogenesis of many diseases. Here, we assessed its role in experimental autoimmune encephalomyelitis (EAE), a widely used non-clinical model of autoimmune inflammatory diseases of the central nervous system (CNS), which resembles many aspects of human multiple sclerosis (MS). We have analyzed in silico data from microglia isolated at different timepoints through disease progression to identify the global editing events and validated the selected targets in murine tissue samples. To further evaluate the functional role of RNA editing, we induced EAE in transgenic animals lacking expression of APOBEC-1. We found that RNA-editing events, mediated by the APOBEC and ADAR deaminases, are significantly reduced throughout the course of disease, possibly affecting the protein expression necessary for normal neurological function. Moreover, the severity of the EAE model was significantly higher in APOBEC-1 knock-out mice, compared to wild-type controls. Our results implicate regulatory epitranscriptomic mechanisms in EAE pathogenesis that could be extrapolated to MS and other neurodegenerative disorders (NDs) with common clinical and molecular features.

Building integrative functional maps of gene regulation

Every cell in the human body inherits a copy of the same genetic information. The three billion base pairs of DNA in the human genome, and the roughly 50 000 coding and non-coding genes they contain, must thus encode all the complexity of human development and cell and tissue type diversity. Differences in gene regulation, or the modulation of gene expression, enable individual cells to interpret the genome differently to carry out their specific functions. Here we discuss recent and ongoing efforts to build gene regulatory maps, which aim to characterize the regulatory roles of all sequences in a genome. Many researchers and consortia have identified such regulatory elements using functional assays and evolutionary analyses; we discuss the results, strengths and shortcomings of their approaches. We also discuss new techniques the field can leverage and emerging challenges it will face while striving to build gene regulatory maps of ever-increasing resolution and comprehensiveness.

RNA modifications: importance in immune cell biology and related diseases

RNA modifications have become hot topics recently. By influencing RNA processes, including generation, transportation, function, and metabolization, they act as critical regulators of cell biology. The immune cell abnormality in human diseases is also a research focus and progressing rapidly these years. Studies have demonstrated that RNA modifications participate in the multiple biological processes of immune cells, including development, differentiation, activation, migration, and polarization, thereby modulating the immune responses and are involved in some immune related diseases. In this review, we present existing knowledge of the biological functions and underlying mechanisms of RNA modifications, including N⁶-methyladenosine (m⁶A), 5-methylcytosine (m⁵C), N¹-methyladenosine (m¹A), N⁷-methylguanosine (m⁷G), N⁴-acetylcytosine (ac⁴C), pseudouridine (Ψ), uridylation, and adenosine-to-inosine (A-to-I) RNA editing, and summarize their critical roles in immune cell biology. Via regulating the biological processes of immune cells, RNA modifications can participate in the pathogenesis of immune related diseases, such as cancers, infection, inflammatory and autoimmune diseases. We further highlight the challenges and future directions based on the existing knowledge. All in all, this review will provide helpful knowledge as well as novel ideas for the researchers in this area.

Reduction of A-to-I RNA editing in the failing human heart regulates formation of circular RNAs

Alterations of RNA editing that affect the secondary structure of RNAs can cause human diseases. We therefore studied RNA editing in failing human hearts. Transcriptome sequencing showed that adenosine-to-inosine (A-to-I) RNA editing was responsible for 80% of the editing events in the myocardium. Failing human hearts were characterized by reduced RNA editing. This was primarily attributable to Alu elements in introns of protein-coding genes. In the failing left ventricle, 166 circRNAs were upregulated and 7 circRNAs were downregulated compared to non-failing controls. Most of the upregulated circRNAs were associated with reduced RNA editing in the host gene. ADAR2, which binds to RNA regions that are edited from A-to-I, was decreased in failing human hearts. In vitro, reduction of ADAR2 increased circRNA levels suggesting a causal effect of reduced ADAR2 levels on increased circRNAs in the failing human heart. To gain mechanistic insight, one of the identified upregulated circRNAs with a high reduction of editing in heart failure, AKAP13, was further characterized. ADAR2 reduced the formation of double-stranded structures in AKAP13 pre-mRNA, thereby reducing the stability of Alu elements and the circularization of the resulting circRNA. Overexpression of circAKAP13 impaired the sarcomere regularity of human induced pluripotent stem cell-derived cardiomyocytes. These data show that ADAR2 mediates A-to-I RNA editing in the human heart. A-to-I RNA editing represses the formation of dsRNA structures of Alu elements favoring canonical linear mRNA splicing and inhibiting the formation of circRNAs. The findings are relevant to diseases with reduced RNA editing and increased circRNA levels and provide insights into the human-specific regulation of circRNA formation.

Molecular Diversity and Phylogeny Reconstruction of Genus Colobanthus (Caryophyllaceae) Based on Mitochondrial Gene Sequences

Mitochondrial genomes have become an interesting object of evolutionary and systematic study both for animals and plants, including angiosperms. Although the framework of the angiosperm phylogeny was built on the information derived from chloroplast and nuclear genes, mitochondrial sequences also revealed their usefulness in solving the phylogenetic issues at different levels of plant systematics. Here, we report for the first time the complete sequences of 26 protein-coding genes of eight Colobanthus species (Caryophyllaceae). Of these, 23 of them represented core mitochondrial genes, which are directly associated with the primary function of that organelle, and the remaining three genes represented a facultative set of mitochondrial genes. Comparative analysis of the identified genes revealed a generally high degree of sequence conservation. The Ka/Ks ratio was <1 for most of the genes, which indicated purifying selection. Only for rps12 was Ka/Ks > 1 in all studied species, suggesting positive selection. We identified 146−165 potential RNA editing sites in genes of the studied species, which is lower than in most angiosperms. The reconstructed phylogeny based on mitochondrial genes was consistent with the taxonomic position of the studied species, showing the separate character of the family Caryophyllaceae and close relationships between all studied Colobanthus species, with C. lycopodioides sharing less similarity.

Genome-wide investigation and functional analysis of RNA editing sites in wheat

Wheat is an important cereal and half of the world population consumed it. Wheat faces environmental stresses and different techniques (CRISPR, gene silencing, GWAS, etc.) were used to enhance its production but RNA editing (RESs) is not fully explored in wheat. RNA editing has a special role in controlling environmental stresses. The genome-wide identification and functional characterization of RESs in different types of wheat genotypes was done. We employed six wheat genotypes by RNA-seq analyses to achieve RESs. The findings revealed that RNA editing events occurred on all chromosomes equally. RNA editing sites were distributed randomly and 10–12 types of RESs were detected in wheat genotypes. Higher number of RESs were detected in drought-tolerant genotypes. A-to-I RNA editing (2952, 2977, 1916, 2576, 3422, and 3459) sites were also identified in six wheat genotypes. Most of the genes were found to be engaged in molecular processes after a Gene Ontology analysis. PPR (pentatricopeptide repeat), OZ1 (organelle zinc-finger), and MORF/RIP gene expression levels in wheat were also examined. Normal growth conditions diverge gene expression of these three different gene families, implying that normal growth conditions for various genotypes can modify RNA editing events and have an impact on gene expression levels. While the expression of PPR genes was not change. We used Variant Effect Predictor (VEP) to annotate RNA editing sites, and Local White had the highest RESs in the CDS region of the protein. These findings will be useful for prediction of RESs in other crops and will be helpful in drought tolerance development in wheat.

A Systematic Review of Common and Brain-Disease-Specific RNA Editing Alterations Providing Novel Insights into Neurological and Neurodegenerative Disease Manifestations

RNA editing contributes to transcriptome diversification through RNA modifications in relation to genome-encoded information (RNA–DNA differences, RDDs). The deamination of Adenosine (A) to Inosine (I) or Cytidine (C) to Uridine (U) is the most common type of mammalian RNA editing. It occurs as a nuclear co- and/or post-transcriptional event catalyzed by ADARs (Adenosine deaminases acting on RNA) and APOBECs (apolipoprotein B mRNA editing enzyme catalytic polypeptide-like genes). RNA editing may modify the structure, stability, and processing of a transcript. This review focuses on RNA editing in psychiatric, neurological, neurodegenerative (NDs), and autoimmune brain disorders in humans and rodent models. We discuss targeted studies that focus on RNA editing in specific neuron-enriched transcripts with well-established functions in neuronal activity, and transcriptome-wide studies, enabled by recent technological advances. We provide comparative editome analyses between human disease and corresponding animal models. Data suggest RNA editing to be an emerging mechanism in disease development, displaying common and disease-specific patterns. Commonly edited RNAs represent potential disease-associated targets for therapeutic and diagnostic values. Currently available data are primarily descriptive, calling for additional research to expand global editing profiles and to provide disease mechanistic insights. The potential use of RNA editing events as disease biomarkers and available tools for RNA editing identification, classification, ranking, and functional characterization that are being developed will enable comprehensive analyses for a better understanding of disease(s) pathogenesis and potential cures.

Characterization of RNA Editome in the Mammary Gland of Yaks during the Lactation and Dry Periods

Simple Summary

In order to study the influence of RNA editing sites on lactation and mammary gland development process in yaks, we comprehensively characterized the RNA editome of the yak mammary gland during the lactation period and dry period by using the transcriptome and genome sequencing data. The results revealed 82,872 nonredundant RNA editing sites, 14,159 of which were differentially edited between the lactation period and dry period. Enrichment analysis showed that the genes harboring differential editing sites were mainly associated with mammary gland development-related pathways, such as MAPK pathway, PI3K-Akt pathway, FoxO signaling pathway, GnRH signaling pathway, and focal adhesion pathway. Our findings offer some novel insights into the RNA editing function in the mammary gland of yaks.

Abstract

The mammary gland is a complicated organ comprising several types of cells, and it undergoes extensive morphogenetic and metabolic changes during the female reproductive cycle. RNA editing is a posttranscriptional modification event occurring at the RNA nucleotide level, and it drives transcriptomic and proteomic diversities, with potential functional consequences. RNA editing in the mammary gland of yaks, however, remains poorly understood. Here, we used REDItools to identify RNA editing sites in mammary gland tissues in yaks during the lactation period (LP, n = 2) and dry period (DP, n = 3). Totally, 82,872 unique RNA editing sites were identified, most of which were detected in the noncoding regions with a low editing degree. In the coding regions (CDS), we detected 5235 editing sites, among which 1884 caused nonsynonymous amino acid changes. Of these RNA editing sites, 486 were found to generate novel possible miRNA target sites or interfere with the initial miRNA binding sites, indicating that RNA editing was related to gene regulation mediated by miRNA. A total of 14,159 RNA editing sites (involving 3238 common genes) showed a significant differential editing level in the LP when compared with that in the DP through Tukey’s Honest Significant Difference method (p < 0.05). According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, genes that showed different RNA editing levels mainly participated in pathways highly related to mammary gland development, including MAPK, PI3K-Akt, FoxO, and GnRH signaling pathways. Collectively, this work demonstrated for the first time the dynamic RNA editome profiles in the mammary gland of yaks and shed more light on the mechanism that regulates lactation together with mammary gland development.

Genetic architecture of microRNA expression and its link to complex diseases in the Japanese population

Understanding the genetic effects on non-coding RNA (ncRNA) expression facilitates functional characterization of disease-associated genetic loci. Among several classes of ncRNAs, microRNAs (miRNAs) are key post-transcriptional gene regulators. Despite its biological importance, previous studies on the genetic architecture of miRNA expression focused mostly on the European individuals, underrepresented in other populations. Here, we mapped miRNA expression quantitative trait loci (miRNA-eQTL) for 343 miRNAs in 141 Japanese using small RNA sequencing (sRNA-seq) and whole-genome sequencing (WGS), identifying 1275 cis-miRNA-eQTL variants for 40 miRNAs (false discovery rate < 0.2). Of these, 25 miRNAs having eQTL were unreported in the European studies, including 5 miRNAs with their lead variant monomorphic in the European populations, which demonstrates the value of miRNA-eQTL analysis in diverse ancestral populations. MiRNAs with eQTL effect showed allele-specific expression (ASE) (e.g. miR-146a-3p), and ASE analysis further detected cis-regulatory variants not captured by the conventional miRNA-eQTL mapping (e.g. miR-933). We identified a copy number variation (CNV) associated with miRNA expression (e.g. miR-570-3p, P = 7.2 × 10-6), which contributes to a more comprehensive landscape of miRNA-eQTLs. To elucidate a post-transcriptional modification in miRNAs, we created a catalog of miRNA-editing sites, including ten canonical and six non-canonical sites. Finally, by integrating the miRNA-eQTLs and Japanese genome-wide association studies of 25 complex traits (mean n = 192 833), we conducted a transcriptome-wide association study (TWAS), identifying miR-1908-5p as a potential mediator for adult height, colorectal cancer, and type 2 diabetes (P < 9.1 × 10-5). Our study broadens the population diversity in ncRNA-eQTL studies and contributes to functional annotation of disease-associated loci found in non-European populations.

Landscape of tissue-specific RNA Editome provides insight into co-regulated and altered gene expression in pigs (Sus-scrofa)

RNA editing generates genetic diversity in mammals by altering amino acid sequences, miRNA targeting site sequences, influencing the stability of targeted RNAs, and causing changes in gene expression. However, the extent to which RNA editing affect gene expression via modifying miRNA binding site remains unexplored. Here, we first profiled the dynamic A-to-I RNA editome across tissues of Duroc and Luchuan pigs. The RNA editing events at the miRNA binding sites were generated. The biological function of the differentially edited gene in skeletal muscle was further characterized in pig muscle-derived satellite cells. RNA editome analysis revealed a total of 171,909 A-to-I RNA editing sites (RESs), and examination of its features showed that these A-to-I editing sites were mainly located in SINE retrotransposons PRE-1/Pre0_SS element. Analysis of differentially edited sites (DESs) revealed a total of 4,552 DESs across tissues between Duroc and Luchuan pigs, and functional category enrichment analysis of differentially edited gene (DEG) sets highlighted a significant association and enrichment of tissue-developmental pathways including TGF-beta, PI3K-Akt, AMPK, and Wnt signaling pathways. Moreover, we found that RNA editing events at the miRNA binding sites in the 3'-UTR of HSPA12B mRNA could prevent the miRNA-mediated mRNA downregulation of HSPA12B in the muscle-derived satellite (MDS) cell, consistent with the results obtained from the Luchuan skeletal muscle. This study represents the most systematic attempt to characterize the significance of RNA editing in regulating gene expression, particularly in skeletal muscle, constituting a new layer of regulation to understand the genetic mechanisms behind phenotype variance in animals.Abbreviations: A-to-I: Adenosine-to-inosine; ADAR: Adenosine deaminase acting on RNA; RES: RNA editing site; DEG: Differentially edited gene; DES: Differentially edited site; FDR: False discovery rate; GO: Gene Ontology; KEGG: Kyoto Encyclopaedia of Genes and Genomes; MDS cell: musclederived satellite cell; RPKM: Reads per kilobase of exon model in a gene per million mapped reads; UTR: Untranslated coding regions.

The first report of RNA U to C or G editing in the mitochondrial NADH dehydrogenase subunit 5 (Nad5) transcript of wild barley

BACKGROUND

Nad dehydrogenase complex in mtDNA has a significant role in cellular respiration. One of the largest subunits in the complex is subunit 5 (Nad5).

METHODS AND RESULTS

Four cDNAs of the Hordeum vulgare subsp. spontaneum nad5 gene have been characterized and subjected to four phases of 0.5 M salinity, at 0 h (control, accession no. MT235236), after 2 h (acc. no. MT235237), after 12 h (acc. no. MT235238) and after 24 h (acc. no. MT235239). Utilizing raw data from RNA-seq, ten RNA editing sites were reported. Seven sites have common editing from C to U in positions (C1490, C1859, C1895, C1900, C1901, C1916, C1918). A rare editing event U to C was detected in two positions (U1650 and U1652) and a novel editing event U to G was for the first time in positions nad5-U231. The highest editing level was shown in 2 and 12 h after salinity exposure. After 24 h, these edits were disrupted, possibly due to the launch of the programed cell death mechanism. However, the RNA editing in positions U1650, U1652 and U231 was fixed at all exposure times.

CONCLUSIONS

Although study clarified the role of salinity stress in nad5 RNA editing sites, the main achievements are first report of U to G RNA editing in plants at position U231 and first report of U to C editing in the nad5 gene at U1650 and U1652.

Identification and Analysis of RNA Editing Events in Ovarian Serous Cystadenoma Using RNA-seq Data

BACKGROUND

Recent studies have revealed thousands of A-to-I RNA editing events in primates. These events are closely related to the occurrence and development of multiple cancers, but the origination and general functions of these events in ovarian cancer remain incompletely understood.

OBJECTIVE

To further the determination of molecular mechanisms of ovarian cancer from the perspective of RNA editing.

METHODS

Here, we used the SNP-free RNA editing Identification Toolkit (SPRINT) to detect RNA editing sites. These editing sites were then annotated, and related functional analysis was performed.

RESULTS

In this study, about 1.7 million RES were detected in each sample, and 98% of these sites were due to A-to-G editing and were mainly distributed in non-coding regions. More than 1,000 A-- to-G RES were detected in CDS regions, and nearly 700 could lead to amino acid changes. Our results also showed that editing in the 3'UTR regions could influence miRNA-target binding. We predicted the network of changed miRNA-mRNA interaction caused by the A-to-I RNA editing sites. We also screened the differential RNA editing sites between ovarian cancer and adjacent normal tissues. We then performed GO and KEGG pathway enrichment analysis on the genes that contained these differential RNA editing sites. Finally, we identified the potential dysregulated RNA editing events in ovarian cancer samples.

CONCLUSION

This study systematically identified and analyzed RNA editing events in ovarian cancer and laid a foundation to explore the regulatory mechanism of RNA editing and its function in ovarian cancer.

Long Non-Coding RNA Epigenetics

Long noncoding RNAs exceeding a length of 200 nucleotides play an important role in ensuring cell functions and proper organism development by interacting with cellular compounds such as miRNA, mRNA, DNA and proteins. However, there is an additional level of lncRNA regulation, called lncRNA epigenetics, in gene expression control. In this review, we describe the most common modified nucleosides found in lncRNA, 6-methyladenosine, 5-methylcytidine, pseudouridine and inosine. The biosynthetic pathways of these nucleosides modified by the writer, eraser and reader enzymes are important to understanding these processes. The characteristics of the individual methylases, pseudouridine synthases and adenine-inosine editing enzymes and the methods of lncRNA epigenetics for the detection of modified nucleosides, as well as the advantages and disadvantages of these methods, are discussed in detail. The final sections are devoted to the role of modifications in the most abundant lncRNAs and their functions in pathogenic processes.

RNA editing enzyme APOBEC3A promotes pro-inflammatory M1 macrophage polarization

Pro-inflammatory M1 macrophage polarization is associated with microbicidal and antitumor responses. We recently described APOBEC3A-mediated cytosine-to-uracil (C > U) RNA editing during M1 polarization. However, the functional significance of this editing is unknown. Here we find that APOBEC3A-mediated cellular RNA editing can also be induced by influenza or Maraba virus infections in normal human macrophages, and by interferons in tumor-associated macrophages. Gene knockdown and RNA_Seq analyses show that APOBEC3A mediates C>U RNA editing of 209 exonic/UTR sites in 203 genes during M1 polarization. The highest level of nonsynonymous RNA editing alters a highly-conserved amino acid in THOC5, which encodes a nuclear mRNA export protein implicated in M-CSF-driven macrophage differentiation. Knockdown of APOBEC3A reduces IL6, IL23A and IL12B gene expression, CD86 surface protein expression, and TNF-α, IL-1β and IL-6 cytokine secretion, and increases glycolysis. These results show a key role of APOBEC3A cytidine deaminase in transcriptomic and functional polarization of M1 macrophages.

Identification and Spatiotemporal Expression of Adenosine Deaminases Acting on RNA (ADAR) during Earthworm Regeneration: Its Possible Implication in Muscle Redifferentiation

Simple Summary

Among the animal species capable of regenerating missing body parts, a species of earthworm, Perionyx excavatus, has the most powerful regeneration capacity, which can completely and regenerate an amputated head and tail. Earthworm regeneration is a form of epimorphosis, a simple mode of development in adults that occurs around the sites of damage rather than throughout the body. In order to achieve this process, the earthworm must have molecular tools via which a variety of cell and tissue types can be precisely recovered from the pluripotent (or possibly totipotent) blastemal cells. Adenosine to inosine (A-to-I) RNA editing catalyzed by adenosine deaminases acting on RNA (ADAR) can generate substantial transcriptome and proteome variability and provide an ideal tool for cell and tissue re-specification. To understand the role of ADAR during earthworm regeneration, the molecular characteristics of an ADAR gene identified from P. excavatus (Pex-ADAR) were analyzed, and its spatial and temporal expression patterns were observed during regeneration. Domain analysis showed that Pex-ADAR is a member of the ADAR1 class. Its expression level primarily increases when and where muscle redifferentiation is actively taking place, suggesting that the RNA-editing enzyme Pex-ADAR is involved in muscle redifferentiation.

Abstract

Adenosine deaminases acting on RNA (ADAR) catalyze the hydrolytic deamination of adenosine (A) to produce inosine (I) in double-stranded RNA substrates. A-to-I RNA editing has increasingly broad physiological significance in development, carcinogenesis, and environmental adaptation. Perionyx excavatus is an earthworm with potent regenerative potential; it can regenerate the head and tail and is an advantageous model system to investigate the molecular mechanisms of regeneration. During RNA sequencing analysis of P. excavatus regenerates, we identified an ADAR homolog (Pex-ADAR), which led us to examine its spatial and temporal expression to comprehend how Pex-ADAR is linked to regeneration. At first, in domain analysis, we discovered that Pex-ADAR only has one double-stranded RNA-binding domain (dsRBD) and a deaminase domain without a Z-DNA-binding domain (ZBD). In addition, a comparison of the core deaminase domains of Pex-ADAR with those of other ADAR family members indicated that Pex-ADAR comprises the conserved three active-site motifs and a glutamate residue for catalytic activity. Pex-ADAR also shares 11 conserved residues, a characteristic of ADAR1, supporting that Pex-ADAR is a member of ADAR1 class. Its temporal expression was remarkably low in the early stages of regeneration before suddenly increasing at 10 days post amputation (dpa) when diverse cell types and tissues were being regenerated. In situ hybridization of Pex-ADAR messenger RNA (mRNA) indicated that the main expression was observed in regenerating muscle layers and related connective tissues. Taken together, the present results demonstrate that an RNA-editing enzyme, Pex-ADAR, is implicated in muscle redifferentiation during earthworm regeneration.

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.

Adenosine deaminases acting on RNA modulate the expression of the human pregnane X receptor

The pregnane X receptor (PXR) is one of the major transcription factors that regulate the expression of different drug-metabolizing enzymes and transporters. Adenosine-to-inosine RNA editing, the most frequent nucleotide conversion on RNA, which is catalyzed by adenosine deaminase acting on RNA (ADAR) enzymes, may modulate gene expression and function. Here, we investigated the potential regulation of human PXR expression by adenosine-to-inosine RNA editing. Knockdown of ADAR1 increased PXR mRNA level, and the knockdown of ADAR1 or ADAR2 significantly increased PXR protein level in HepaRG cells. In HepG2 cells, the knockdown of ADAR1 or ADAR2 significantly increased PXR mRNA and protein levels. The increase in the PXR protein by ADAR1 knockdown resulted in increased cytochrome P450 3A4 (CYP3A4) transactivity and CYP3A4 and UDP-glucuronosyltransferase 1A1 (UGT1A1) expression. A reporter assay revealed that the 3'-untranslated region (UTR) of PXR mRNA, especially from +3371 to +3440, is responsible for the ADAR-mediated post-transcriptional control of PXR expression, despite the lack of RNA edited sites in this region. Collectively, we found that PXR is negatively regulated by ADAR1 via an indirect mechanism, which facilitates the degradation of PXR mRNA. We could demonstrate that ADAR1 can cause interindividual variability in hepatic drug metabolism potencies.

Decrease in ADAR1 expression by exposure to cigarette smoke enhances susceptibility to oxidative stress

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by adenosine deaminase acting on RNA (ADAR) enzymes, is the most frequent type of post-transcriptional nucleotide conversion in humans. It is known that innate abnormalities of A-to-I RNA editing are associated with the risk of certain diseases, such as amyotrophic lateral sclerosis. Extrinsic factors that modulate ADAR-mediated RNA editing remain to be clarified. In this study, we investigated the possibility that cigarette smoking may influence the expression of ADAR and that the changes may be biologically significant. Treatment of human lung adenocarcinoma A549 cells with cigarette smoke extract (CSE) induced a significant 50% decrease in ADAR1 protein levels. Since the decrease was counteracted by cotreatment with chloroquine, the CSE-dependent decrease in the ADAR1 protein levels may be due to the activation of autophagy. In addition to the in vitro study, we performed an in vivo study in mice and found a decrease in pulmonary Adar1 protein expression induced by cigarette smoking. Then, we investigated the biological significance of decreased ADAR1 expression. We found that knockdown of ADAR1 in A549 cells by siRNA resulted in an increase in the levels of protein carbonyl, a marker of oxidative stress. Moreover, knockdown of ADAR1 triggered a decrease in super oxide dismutase activity and heme oxygenase-1 expression, suggesting that ADAR1 plays a role to suppress oxidative stress. In conclusion, we show that ADAR1 expression is decreased by cigarette smoking and is a factor that contributes to the enhanced intracellular oxidative stress caused by cigarette smoking.

Light/heat effects on RNA editing in chloroplast NADH-plastoquinone oxidoreductase subunit 2 (ndhB) gene of Calotropis (Calotropis procera)

Background

RNA editing is common in terrestrial plants, especially in mitochondria and chloroplast. In the photosynthesis process, NAD dehydrogenase plays a very important role. Subunit 2 of NADH-dehydrogenase is one of the major subunits in NAD dehydrogenase complex. Using desert plant Calotropis (Calotropis procera), this study focuses on the RNA editing activity of ndhB based on light time.

Results

NdhB (NADH-dehydrogenase subunit 2) gene accession no. MK144329 was isolated from Calotropis procera genomic data (PRJNA292713). Additionally, using RNA-seq data, the cDNA of the ndhB gene of C. procera was isolated at three daylight periods, i.e., dawn (accession no. MK165161), at midday (accession no. MK165160), and pre-dusk (accession no. MK165159). Seven RNA editing sites have been found in several different positions (nucleotide no. C467, C586, C611, C737, C746, C830, and C1481) within the ndhB coding region. The rate of these alterations was deferentially edited across the three daylight periods. RNA editing rate of ndhB gene was highest at dawn, (87.5, 79.6, 78.5, 76, 68.6, 39.3, and 96.9%, respectively), less in midday (74.8, 54.1, 62.6, 47.4, 45.5, 47.4, and 93.4%, respectively), and less at pre-dusk (67, 52.6, 56.9, 40.1, 40.7, 33.2, and 90%, respectively), also all these sites were validated by qRT-PCR.

Conclusion

The differential editing of chloroplast ndhB gene across light periods may be led to a somehow relations between the RNA editing and control of photosynthesis.

Sex-Dependent RNA Editing and N6-adenosine RNA Methylation Profiling in the Gonads of a Fish, the Olive Flounder (Paralichthys olivaceus)

Adenosine-to-inosine (A-to-I) editing and N6-methyladenosine (m6A) are two of the most abundant RNA modifications. Here, we examined the characteristics of the RNA editing and transcriptome-wide m6A modification profile in the gonads of the olive flounder, Paralichthys olivaceus, an important maricultured fish in Asia. The gonadal differentiation and development of the flounder are controlled by genetic as well as environmental factors, and the epigenetic mechanism may play an important role. In total, 742 RNA editing events were identified, 459 of which caused A to I conversion. Most A-to-I sites were located in 3′UTRs, while 61 were detected in coding regions (CDs). The number of editing sites in the testis was higher than that in the ovary. Transcriptome-wide analyses showed that more than one-half of the transcribed genes presented an m6A modification in the flounder gonads, and approximately 60% of the differentially expressed genes (DEGs) between the testis and ovary appeared to be negatively correlated with m6A methylation enrichment. Further analyses revealed that the mRNA expression of some sex-related genes (e.g., dmrt1 and amh) in the gonads may be regulated by changes in mRNA m6A enrichment. Functional enrichment analysis indicated that the RNA editing and m6A modifications were enriched in several canonical pathways (e.g., Wnt and MAPK signaling pathways) in fish gonads and in some pathways whose roles have not been investigated in relation to fish sex differentiation and gonadal development (e.g., PPAR and RNA degradation pathways). There were 125 genes that were modified by both A-to-I editing and m6A, but the two types of modifications mostly occurred at different sites. Our results suggested that the presence of sex-specific RNA modifications may be involved in the regulation of gonadal development and gametogenesis.

Novel Insights into the Role of UBE3A in Regulating Apoptosis and Proliferation

The UBE3A gene codes for a protein with two known functions, a ubiquitin E3-ligase which catalyzes ubiquitin binding to substrate proteins and a steroid hormone receptor coactivator. UBE3A is most famous for its critical role in neuronal functioning. Lack of UBE3A protein expression leads to Angelman syndrome (AS), while its overexpression is associated with autism. In spite of extensive research, our understanding of UBE3A roles is still limited. We investigated the cellular and molecular effects of Ube3a deletion in mouse embryonic fibroblasts (MEFs) and Angelman syndrome (AS) mouse model hippocampi. Cell cultures of MEFs exhibited enhanced proliferation together with reduced apoptosis when Ube3a was deleted. These findings were supported by transcriptome and proteome analyses. Furthermore, transcriptome analyses revealed alterations in mitochondria-related genes. Moreover, an analysis of adult AS model mice hippocampi also found alterations in the expression of apoptosis- and proliferation-associated genes. Our findings emphasize the role UBE3A plays in regulating proliferation and apoptosis and sheds light into the possible effects UBE3A has on mitochondrial involvement in governing this balance.

Editome Disease Knowledgebase (EDK): a curated knowledgebase of editome-disease associations in human

RNA editing, as an essential co-/post-transcriptional RNA modification type, plays critical roles in many biological processes and involves with a variety of human diseases. Although several databases have been developed to collect RNA editing data in both model and non-model animals, there still lacks a resource integrating associations between editome and human disease. In this study, we present Editome-Disease Knowledgebase (EDK; http://bigd.big.ac.cn/edk), an integrated knowledgebase of RNA editome-disease associations manually curated from published literatures. In the current version, EDK incorporates 61 diseases associated with 248 experimentally validated abnormal editing events located in 32 mRNAs, 16 miRNAs, 1 lncRNA and 11 viruses, and 44 aberrant activities involved with 6 editing enzymes, which together are curated from more than 200 publications. In addition, to facilitate standardization of editome-disease knowledge integration, we propose a data curation model in EDK, factoring an abundance of relevant information to fully capture the context of editome-disease associations. Taken together, EDK is a comprehensive collection of editome-disease associations and bears the great utility in aid of better understanding the RNA editing machinery and complex molecular mechanisms associated with human diseases.

Salinity effects on nad3 gene RNA editing of wild barley mitochondria

Nad complex plays a very important role during cellular respiration. nad3 (nad dehydrogenase subunit 3) is one of the biggest subunits in this complex. Four cDNAs of nad3 gene were characterized in Hordeum vulgare subsp. spontaneum at exposed to four periods of 500 mM salinity, 0 h or control (accession no. MN066165), after 2 h (accession no. MN066166), after 12 h (accession no. MN066167) and after 24 h (accession no. MN066168) using RNA-seq raw data. Seventeen RNA editing sites were found in positions (or nucleotide nos. C5, C39, C44, C61, C62, C79, C80, C147, C185, C190, C191, C208, C209, C275, C317, C344, C349) within the nad3 coding region. These alterations represent differential editing at four exposure times. The maximum editing rate was revealed 2 and 12 h after salinity exposure. However, these edits were disrupted after 24 h probably due to the initiation of program cell death machinery. We found that RNA editing not only improved protein function but also may improve codon bias by altering the nucleotide without any change in amino acid. Characterization of pentatricopeptide repeat-containing protein At4g13650 (PPRSp1) in wild barley helped us to understand the behavior of editing sites C190 and C191 under salinity. Position - 6 in cis-element upstream editing sites of C155, C190 and C191 may be vital to the editing process in these sites by PPRSp1 protein. The differential editing of this gene under salinity led to a relationship between RNA editing and cellular respiration regulation.

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.

FairBase: a comprehensive database of fungal A-to-I RNA editing

Frequent A-to-I RNA editing has recently been identified in fungi despite the absence of recognizable homologues of metazoan ADARs ("Adenosine Deaminases Acting on RNA"). In particular, there is emerging evidence showing that A-to-I editing is involved in sexual reproduction of filamentous fungi. Here, we report on the creation of FairBase - a fungal A-to-I RNA editing database that provides a platform for deep exploration of fungal RNA editing to relevant academic communities. This database includes a comprehensive collection of A-to-I editing sites in six filamentous fungal species, together with extensive annotations for each editing site. In FairBase, users can conveniently search editing sites and obtain editing levels for each editing site in various RNA-seq samples. In addition, the pathways involving RNA editing are built in FairBase to help users understand the functions of RNA editing. Furthermore, each fungal species has a genome browser (JBrowse) that allows users to explore A-to-I editing in a genomic context. FairBase is the first fungal RNA editing database.

Molecular model linking Th2 polarized M2 tumour-associated macrophages with deaminase-mediated cancer progression mutation signatures

A new and diverse range of somatic mutation signatures are observed in late-stage cancers, but the underlying reasons are not fully understood. We advance a "combinatorial association model" for deaminase binding domain (DBD) diversification to explain the generation of previously observed cancer-progression associated mutation signatures. We also propose that changes in the polarization of tumour-associated macrophages (TAMs) are accompanied by the expression of deaminases with a new and diverse range of DBDs, and thus accounting for the generation of new somatic mutation signatures. The mechanism proposed is molecularly reminiscent of combinatorial association of heavy (H) and light (L) protein chains following V(D)J recombination of immunoglobulin molecules (and similarly for protein chains in heterodimers α/β and γ/δ of V(D)Js of T Cell Receptors) required for pathogen antigen recognition by B cells and T cells, respectively. We also discuss whether extracellular vesicles (EVs) emanating from tumour enhancing M2-polarized macrophages represent a likely source of the de novo deaminase DBDs. We conclude that M2-polarized macrophages extruding EVs loaded with deaminase proteins or deaminase-specific transcription/translation regulatory factors and like information may directly trigger deaminase diversification within cancer cells, and thus account for the many new somatic mutation signatures that are indicative of cancer progression. This hypothesis now has a plausible evidentiary base, and it is worth direct testing in future investigations. A long-term objective would be to identify molecular biomarkers predicting cancer progression (or metastatic disease) and to support the development of new drug targets before metastatic pathways are activated.

Genome-wide identification of RNA editing in seven porcine tissues by matched DNA and RNA high-throughput sequencing

BACKGROUND

RNA editing is a co/posttranscriptional modification mechanism that increases the diversity of transcripts, with potential functional consequences. The advent of next-generation sequencing technologies has enabled the identification of RNA edits at unprecedented throughput and resolution. However, our knowledge of RNA editing in swine is still limited.

RESULTS

Here, we utilized RES-Scanner to identify RNA editing sites in the brain, subcutaneous fat, heart, liver, muscle, lung and ovary in three 180-day-old Large White gilts based on matched strand-specific RNA sequencing and whole-genome resequencing datasets. In total, we identified 74863 editing sites, and 92.1% of these sites caused adenosine-to-guanosine (A-to-G) conversion. Most A-to-G sites were located in noncoding regions and generally had low editing levels. In total, 151 A-to-G sites were detected in coding regions (CDS), including 94 sites that could lead to nonsynonymous amino acid changes. We provide further evidence supporting a previous observation that pig transcriptomes are highly editable at PRE-1 elements. The number of A-to-G editing sites ranged from 4155 (muscle) to 25001 (brain) across the seven tissues. The expression levels of the ADAR enzymes could explain some but not all of this variation across tissues. The functional analysis of the genes with tissue-specific editing sites in each tissue revealed that RNA editing might play important roles in tissue function. Specifically, more pathways showed significant enrichment in the fat and liver than in other tissues, while no pathway was enriched in the muscle.

CONCLUSIONS

This study identified a total of 74863 nonredundant RNA editing sites in seven tissues and revealed the potential importance of RNA editing in tissue function. Our findings largely extend the porcine editome and enhance our understanding of RNA editing in swine.

Human C-to-U Coding RNA Editing Is Largely Nonadaptive

C-to-U RNA editing enzymatically converts the base C to U in RNA molecules and could lead to nonsynonymous changes when occurring in coding regions. Hundreds to thousands of coding sites were recently found to be C-to-U edited or editable in humans, but the biological significance of this phenomenon is elusive. Here, we test the prevailing hypothesis that nonsynonymous editing is beneficial because it provides a means for tissue- or time-specific regulation of protein function that may be hard to accomplish by mutations due to pleiotropy. The adaptive hypothesis predicts that the fraction of sites edited and the median proportion of RNA molecules edited (i.e., editing level) are both higher for nonsynonymous than synonymous editing. However, our empirical observations are opposite to these predictions. Furthermore, the frequency of nonsynonymous editing, relative to that of synonymous editing, declines as genes become functionally more important or evolutionarily more constrained, and the nonsynonymous editing level at a site is negatively correlated with the evolutionary conservation of the site. Together, these findings refute the adaptive hypothesis; they instead indicate that the reported C-to-U coding RNA editing is mostly slightly deleterious or neutral, probably resulting from off-target activities of editing enzymes. Along with similar conclusions on the more prevalent A-to-I editing and m6A modification of coding RNAs, our study suggests that, at least in humans, most events of each type of posttranscriptional coding RNA modification likely manifest cellular errors rather than adaptations, demanding a paradigm shift in the research of posttranscriptional modification.

RNA Biology Provides New Therapeutic Targets for Human Disease

RNA is the messenger molecule that conveys information from the genome and allows the production of biomolecules required for life in a responsive and regulated way. Most genes are able to produce multiple mRNA products in response to different internal or external environmental signals, in different tissues and organs, and at specific times in development or later life. This fine tuning of gene expression is dependent on the coordinated effects of a large and intricate set of regulatory machinery, which together orchestrate the genomic output at each locus and ensure that each gene is expressed at the right amount, at the right time and in the correct location. This complexity of control, and the requirement for both sequence elements and the entities that bind them, results in multiple points at which errors may occur. Errors of RNA biology are common and found in association with both rare, single gene disorders, but also more common, chronic diseases. Fortunately, complexity also brings opportunity. The existence of many regulatory steps also offers multiple levels of potential therapeutic intervention which can be exploited. In this review, I will outline the specific points at which coding RNAs may be regulated, indicate potential means of intervention at each stage, and outline with examples some of the progress that has been made in this area. Finally, I will outline some of the remaining challenges with the delivery of RNA-based therapeutics but indicate why there are reasons for optimism.

Curcumin modulates the apolipoprotein B mRNA editing by coordinating the expression of cytidine deamination to uridine editosome components in primary mouse hepatocytes

Curcumin, an active ingredient of Curcuma longa L., can reduce the concentration of low-density lipoproteins in plasma, in different ways. We had first reported that curcumin exhibits hypocholesterolemic properties by improving the apolipoprotein B (apoB) mRNA editing in primary rat hepatocytes. However, the role of curcumin in the regulation of apoB mRNA editing is not clear. Thus, we investigated the effect of curcumin on the expression of multiple editing components of apoB mRNA cytidine deamination to uridine (C-to-U) editosome. Our results demonstrated that treatment with 50 µM curcumin markedly increased the amount of edited apoB mRNA in primary mouse hepatocytes from 5.13%–8.05% to 27.63%–35.61%, and significantly elevated the levels of the core components apoB editing catalytic polypeptide-1 (APOBEC-1), apobec-1 complementation factor (ACF), and RNA-binding-motif-protein-47 (RBM47), as well as suppressed the level of the inhibitory component glycine-arginine-tyrosine-rich RNA binding protein. Moreover, the increased apoB RNA editing by 50 µM curcumin was significantly reduced by siRNA-mediated APOBEC-1, ACF, and RBM47 knockdown. These findings suggest that curcumin modulates apoB mRNA editing by coordinating the multiple editing components of the editosome in primary hepatocytes. Our data provided evidence for curcumin to be used therapeutically to prevent atherosclerosis.

Enhanced AZIN1 RNA editing and overexpression of its regulatory enzyme ADAR1 are important prognostic biomarkers in gastric cancer

BACKGROUND

Adenosine-to-inosine (A-to-I) RNA editing is catalyzed by adenosine deaminases acting on RNA (ADAR) enzymes. Recent evidence suggests that RNA editing of antizyme inhibitor 1 (AZIN1) RNA is emerging as a key epigenetic alteration underlying cancer pathogenesis.

METHODS

We evaluated AZIN1 RNA editing levels, and the expression of its regulator, ADAR1, in 280 gastric tissues from 140 patients, using a RNA editing site-specific quantitative polymerase chain reaction assays. We also analyzed the clinical significance of these results as disease biomarkers in gastric cancer (GC) patients.

RESULTS

Both AZIN1 RNA editing levels and ADAR1 expression were significantly elevated in GC tissues compared with matched normal mucosa (P < 0.0001, 0.0008, respectively); and AZIN1 RNA editing was positively correlated with ADAR1 expression. Elevated expression of ADAR1 significantly correlated with poor overall survival (P = 0.034), while hyper-edited AZIN1 emerged as an independent prognostic factor for OS and disease-free survival in GC patients [odds ratio (OR):1.98, 95% CI 1.17-3.35, P = 0.011, OR: 4.55, 95% CI 2.12-9.78, P = 0.0001, respectively]. Increased AZIN1 RNA editing and ADAR1 over-expression were significantly correlated with key clinicopathological factors, such as advanced T stage, presence of lymph node metastasis, distant metastasis, and higher TNM stages in GC patients. Logistic regression analysis revealed that hyper-editing status of AZIN1 RNA was an independent risk factor for lymph node metastasis in GC patients [hazard ratio (HR):3.03, 95% CI 1.19-7.71, P = 0.02].

CONCLUSIONS

AZIN1 RNA editing levels may be an important prognostic biomarker in GC patients, and may serve as a key clinical decision-making tool for determining preoperative treatment strategies in GC patients.

SPRINT: an SNP-free toolkit for identifying RNA editing sites

Motivation

RNA editing generates post-transcriptional sequence alterations. Detection of RNA editing sites (RESs) typically requires the filtering of SNVs called from RNA-seq data using an SNP database, an obstacle that is difficult to overcome for most organisms.

Results

Here, we present a novel method named SPRINT that identifies RESs without the need to filter out SNPs. SPRINT also integrates the detection of hyper RESs from remapped reads, and has been fully automated to any RNA-seq data with reference genome sequence available. We have rigorously validated SPRINT’s effectiveness in detecting RESs using RNA-seq data of samples in which genes encoding RNA editing enzymes are knock down or over-expressed, and have also demonstrated its superiority over current methods. We have applied SPRINT to investigate RNA editing across tissues and species, and also in the development of mouse embryonic central nervous system. A web resource (http://sprint.tianlab.cn) of RESs identified by SPRINT has been constructed.

Availability and implementation

The software and related data are available at http://sprint.tianlab.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.

Combinatory RNA-Sequencing Analyses Reveal a Dual Mode of Gene Regulation by ADAR1 in Gastric Cancer

BACKGROUND

Adenosine deaminase acting on RNA 1 (ADAR1) is known to mediate deamination of adenosine-to-inosine through binding to double-stranded RNA, the phenomenon known as RNA editing. Currently, the function of ADAR1 in gastric cancer is unclear.

AIMS

This study was aimed at investigating RNA editing-dependent and editing-independent functions of ADAR1 in gastric cancer, especially focusing on its influence on editing of 3' untranslated regions (UTRs) and subsequent changes in expression of messenger RNAs (mRNAs) as well as microRNAs (miRNAs).

METHODS

RNA-sequencing and small RNA-sequencing were performed on AGS and MKN-45 cells with a stable ADAR1 knockdown. Changed frequencies of editing and mRNA and miRNA expression were then identified by bioinformatic analyses. Targets of RNA editing were further validated in patients' samples.

RESULTS

In the Alu region of both gastric cell lines, editing was most commonly of the A-to-I type in 3'-UTR or intron. mRNA and protein levels of PHACTR4 increased in ADAR1 knockdown cells, because of the loss of seed sequences in 3'-UTR of PHACTR4 mRNA that are required for miRNA-196a-3p binding. Immunohistochemical analyses of tumor and paired normal samples from 16 gastric cancer patients showed that ADAR1 expression was higher in tumors than in normal tissues and inversely correlated with PHACTR4 staining. On the other hand, decreased miRNA-148a-3p expression in ADAR1 knockdown cells led to increased mRNA and protein expression of NFYA, demonstrating ADAR1's editing-independent function.

CONCLUSIONS

ADAR1 regulates post-transcriptional gene expression in gastric cancer through both RNA editing-dependent and editing-independent mechanisms.

A-to-I RNA editing of BLCAP lost the inhibition to STAT3 activation in cervical cancer

Bladder cancer-associated protein (BLCAP) gene is a highly conserved gene with tumor-suppressor function in different carcinomas. It is also a novel ADAR-mediated editing substrate undergoes multiple A-to-I RNA editing events. Although the anti-tumorigenic role of BLCAP has been examined in preliminarily studies, the relationship between BLCAP function and A-to-I RNA editing in cervical carcinogenesis still require further exploration. Herein, we analyzed the coding sequence of BLCAP transcripts in 35 paired cervical cancer samples using high-throughput sequencing. Of note, editing levels of three novel editing sites were statistically different between cancerous and adjacent cervical tissues, and editing of these three sites was closely correlated. Moreover, two editing sites of BLCAP coding region were mapped-in the key YXXQ motif which can bind to SH2 domain of STAT3. Further studies revealed that BLCAP interacted with signal transducer and activator of transcription 3 (STAT3) and inhibited its phosphorylation, while A-to-I RNA editing of BLCAP lost the inhibition to STAT3 activation in cervical cancer cell lines. Our findings reveal that A-to-I RNA editing events alter the genetically coded amino acid in BLCAP YXXQ motif, which drive the progression of cervical carcinogenesis through regulating STAT3 signaling pathway.

Global analysis of A-to-I RNA editing reveals association with common disease variants

RNA editing modifies transcripts and may alter their regulation or function. In humans, the most common modification is adenosine to inosine (A-to-I). We examined the global characteristics of RNA editing in 4,301 human tissue samples. More than 1.6 million A-to-I edits were identified in 62% of all protein-coding transcripts. mRNA recoding was extremely rare; only 11 novel recoding sites were uncovered. Thirty single nucleotide polymorphisms from genome-wide association studies were associated with RNA editing; one that influences type 2 diabetes (rs2028299) was associated with editing in ARPIN. Twenty-five genes, including LRP11 and PLIN5, had editing sites that were associated with plasma lipid levels. Our findings provide new insights into the genetic regulation of RNA editing and establish a rich catalogue for further exploration of this process.

Genome-wide analysis of differential RNA editing in epilepsy

The recoding of genetic information through RNA editing contributes to proteomic diversity, but the extent and significance of RNA editing in disease is poorly understood. In particular, few studies have investigated the relationship between RNA editing and disease at a genome-wide level. Here, we developed a framework for the genome-wide detection of RNA sites that are differentially edited in disease. Using RNA-sequencing data from 100 hippocampi from mice with epilepsy (pilocarpine–temporal lobe epilepsy model) and 100 healthy control hippocampi, we identified 256 RNA sites (overlapping with 87 genes) that were significantly differentially edited between epileptic cases and controls. The degree of differential RNA editing in epileptic mice correlated with frequency of seizures, and the set of genes differentially RNA-edited between case and control mice were enriched for functional terms highly relevant to epilepsy, including “neuron projection” and “seizures.” Genes with differential RNA editing were preferentially enriched for genes with a genetic association to epilepsy. Indeed, we found that they are significantly enriched for genes that harbor nonsynonymous de novo mutations in patients with epileptic encephalopathy and for common susceptibility variants associated with generalized epilepsy. These analyses reveal a functional convergence between genes that are differentially RNA-edited in acquired symptomatic epilepsy and those that contribute risk for genetic epilepsy. Taken together, our results suggest a potential role for RNA editing in the epileptic hippocampus in the occurrence and severity of epileptic seizures.