Structure-mediated modulation of mRNA abundance by A-to-I editing

Functions and mechanisms of A-to-I RNA editing in filamentous ascomycetes

Although lack of ADAR (adenosine deaminase acting on RNA) orthologs, genome-wide A-to-I editing occurs specifically during sexual reproduction in a number of filamentous ascomycetes, including Fusarium graminearum and Neurospora crassa. Unlike ADAR-mediated editing in animals, fungal A-to-I editing has a strong preference for hairpin loops and U at -1 position, which leads to frequent editing of UAG and UAA stop codons. Majority of RNA editing events in fungi are in the coding region and cause amino acid changes. Some of these editing events have been experimentally characterized for providing heterozygote and adaptive advantages in F. graminearum. Recent studies showed that FgTad2 and FgTad3, 2 ADAT (adenosine deaminase acting on tRNA) enzymes that normally catalyze the editing of A34 in the anticodon of tRNA during vegetative growth mediate A-to-I mRNA editing during sexual reproduction. Stage specificity of RNA editing is conferred by stage-specific expression of short transcript isoforms of FgTAD2 and FgTAD3 as well as cofactors such as AME1 and FIP5 that facilitate the editing of mRNA in perithecia. Taken together, fungal A-to-I RNA editing during sexual reproduction is catalyzed by ADATs and it has the same sequence and structural preferences with editing of A34 in tRNA.

Navigating the landscape of epitranscriptomics and host immunity

RNA modifications, also termed epitranscriptomic marks, encompass chemical alterations to individual nucleotides, including processes such as methylation and editing. These marks contribute to a wide range of biological processes, many of which are related to host immune system defense. The functions of immune-related RNA modifications can be categorized into three main groups: regulation of immunogenic RNAs, control of genes involved in innate immune response, and facilitation of adaptive immunity. Here, we provide an overview of recent research findings that elucidate the contributions of RNA modifications to each of these processes. We also discuss relevant methods for genome-wide identification of RNA modifications and their immunogenic substrates. Finally, we highlight recent advances in cancer immunotherapies that aim to reduce cancer cell viability by targeting the enzymes responsible for RNA modifications. Our presentation of these dynamic research avenues sets the stage for future investigations in this field.

Structural and functional effects of inosine modification in mRNA

Inosine (I), resulting from the deamination of adenosine (A), is a prominent modification in the human transcriptome. The enzymes responsible for the conversion of adenosine to inosine in human mRNAs are the ADARs (adenosine deaminases acting on RNA). Inosine modification introduces a layer of complexity to mRNA processing and function, as it can impact various aspects of RNA biology, including mRNA stability, splicing, translation, and protein binding. The relevance of this process is emphasized in the growing number of human disorders associated with dysregulated A-to-I editing pathways. Here, we describe the impact of the A-to-I conversion on the structure and stability of duplex RNA and on the consequences of this modification at different locations in mRNAs. Furthermore, we highlight specific open questions regarding the interplay between inosine formation in duplex RNA and the innate immune response.

The miRNA-target interactions: An underestimated intricacy

MicroRNAs (miRNAs) play indispensable roles in posttranscriptional gene regulation. Their cellular regulatory impact is determined not solely by their sheer number, which likely amounts to >2000 individual miRNAs in human, than by the regulatory effectiveness of single miRNAs. Although, one begins to develop an understanding of the complex mechanisms underlying miRNA-target interactions (MTIs), the overall knowledge of MTI functionality is still rather patchy. In this critical review, we summarize key features of mammalian MTIs. We especially highlight latest insights on (i) the dynamic make-up of miRNA binding sites including non-canonical binding sites, (ii) the cooperativity between miRNA binding sites, (iii) the adaptivity of MTIs through sequence modifications, (iv) the bearing of intra-cellular miRNA localization changes and (v) the role of cell type and cell status specific miRNA interaction partners. The MTI biology is discussed against the background of state-of-the-art approaches with particular emphasis on experimental strategies for evaluating miRNA functionality.

Malignant A-to-I RNA editing by ADAR1 drives T cell acute lymphoblastic leukemia relapse via attenuating dsRNA sensing

Leukemia-initiating cells (LICs) are regarded as the origin of leukemia relapse and therapeutic resistance. Identifying direct stemness determinants that fuel LIC self-renewal is critical for developing targeted approaches. Here, we show that the RNA-editing enzyme ADAR1 is a crucial stemness factor that promotes LIC self-renewal by attenuating aberrant double-stranded RNA (dsRNA) sensing. Elevated adenosine-to-inosine editing is a common attribute of relapsed T cell acute lymphoblastic leukemia (T-ALL) regardless of molecular subtype. Consequently, knockdown of ADAR1 severely inhibits LIC self-renewal capacity and prolongs survival in T-ALL patient-derived xenograft models. Mechanistically, ADAR1 directs hyper-editing of immunogenic dsRNA to avoid detection by the innate immune sensor melanoma differentiation-associated protein 5 (MDA5). Moreover, we uncover that the cell-intrinsic level of MDA5 dictates the dependency on the ADAR1-MDA5 axis in T-ALL. Collectively, our results show that ADAR1 functions as a self-renewal factor that limits the sensing of endogenous dsRNA. Thus, targeting ADAR1 presents an effective therapeutic strategy for eliminating T-ALL LICs.

Epitranscriptomics as a New Layer of Regulation of Gene Expression in Skeletal Muscle: Known Functions and Future Perspectives

Epitranscriptomics refers to post-transcriptional regulation of gene expression via RNA modifications and editing that affect RNA functions. Many kinds of modifications of mRNA have been described, among which are N6-methyladenosine (m6A), N1-methyladenosine (m1A), 7-methylguanosine (m7G), pseudouridine (Ψ), and 5-methylcytidine (m5C). They alter mRNA structure and consequently stability, localization and translation efficiency. Perturbation of the epitranscriptome is associated with human diseases, thus opening the opportunity for potential manipulations as a therapeutic approach. In this review, we aim to provide an overview of the functional roles of epitranscriptomic marks in the skeletal muscle system, in particular in embryonic myogenesis, muscle cell differentiation and muscle homeostasis processes. Further, we explored high-throughput epitranscriptome sequencing data to identify RNA chemical modifications in muscle-specific genes and we discuss the possible functional role and the potential therapeutic applications.

Genetic pathways regulating the longitudinal acquisition of cocaine self-administration in a panel of inbred and recombinant inbred mice

To identify addiction genes, we evaluate intravenous self-administration of cocaine or saline in 84 inbred and recombinant inbred mouse strains over 10 days. We integrate the behavior data with brain RNA-seq data from 41 strains. The self-administration of cocaine and that of saline are genetically distinct. We maximize power to map loci for cocaine intake by using a linear mixed model to account for this longitudinal phenotype while correcting for population structure. A total of 15 unique significant loci are identified in the genome-wide association study. A transcriptome-wide association study highlights the Trpv2 ion channel as a key locus for cocaine self-administration as well as identifying 17 additional genes, including Arhgef26, Slc18b1, and Slco5a1. We find numerous instances where alternate splice site selection or RNA editing altered transcript abundance. Our work emphasizes the importance of Trpv2, an ionotropic cannabinoid receptor, for the response to cocaine.

RNA modifications in cancer

Currently, more than 170 modifications have been identified on RNA. Among these RNA modifications, various methylations account for two-thirds of total cases and exist on almost all RNAs. Roles of RNA modifications in cancer are garnering increasing interest. The research on m6A RNA methylation in cancer is in full swing at present. However, there are still many other popular RNA modifications involved in the regulation of gene expression post-transcriptionally besides m6A RNA methylation. In this review, we focus on several important RNA modifications including m1A, m5C, m7G, 2'-O-Me, Ψ and A-to-I editing in cancer, which will provide a new perspective on tumourigenesis by peeking into the complex regulatory network of epigenetic RNA modifications, transcript processing, and protein translation.

L-GIREMI uncovers RNA editing sites in long-read RNA-seq

Although long-read RNA-seq is increasingly applied to characterize full-length transcripts it can also enable detection of nucleotide variants, such as genetic mutations or RNA editing sites, which is significantly under-explored. Here, we present an in-depth study to detect and analyze RNA editing sites in long-read RNA-seq. Our new method, L-GIREMI, effectively handles sequencing errors and read biases. Applied to PacBio RNA-seq data, L-GIREMI affords a high accuracy in RNA editing identification. Additionally, our analysis uncovered novel insights about RNA editing occurrences in single molecules and double-stranded RNA structures. L-GIREMI provides a valuable means to study nucleotide variants in long-read RNA-seq.

The Potential Regulation of A-to-I RNA Editing on Genes in Parkinson's Disease

Parkinson's disease (PD) is characterized by dopaminergic neurodegeneration and an abnormal accumulation of α-synuclein aggregates. A number of genetic factors have been shown to increase the risk of PD. Exploring the underlying molecular mechanisms that mediate PD's transcriptomic diversity can help us understand neurodegenerative pathogenesis. In this study, we identified 9897 A-to-I RNA editing events associated with 6286 genes across 372 PD patients. Of them, 72 RNA editing events altered miRNA binding sites and this may directly affect miRNA regulations of their host genes. However, RNA editing effects on the miRNA regulation of genes are more complex. They can (1) abolish existing miRNA binding sites, which allows miRNAs to regulate other genes; (2) create new miRNA binding sites that may sequester miRNAs from regulating other genes; or (3) occur in the miRNA seed regions and change their targets. The first two processes are also referred to as miRNA competitive binding. In our study, we found 8 RNA editing events that may alter the expression of 1146 other genes via miRNA competition. We also found one RNA editing event that modified a miRNA seed region, which was predicted to disturb the regulation of four genes. Considering the PD-related functions of the affected genes, 25 A-to-I RNA editing biomarkers for PD are proposed, including the 3 editing events in the EIF2AK2, APOL6, and miR-4477b seed regions. These biomarkers may alter the miRNA regulation of 133 PD-related genes. All these analyses reveal the potential mechanisms and regulations of RNA editing in PD pathogenesis.

Widespread RNA hypoediting in schizophrenia and its relevance to mitochondrial function

RNA editing, the endogenous modification of nucleic acids, is known to be altered in genes with important neurological function in schizophrenia (SCZ). However, the global profile and molecular functions of disease-associated RNA editing remain unclear. Here, we analyzed RNA editing in postmortem brains of four SCZ cohorts and uncovered a significant and reproducible trend of hypoediting in patients of European descent. We report a set of SCZ-associated editing sites via WGCNA analysis, shared across cohorts. Using massively parallel reporter assays and bioinformatic analyses, we observed that differential 3' untranslated region (3'UTR) editing sites affecting host gene expression were enriched for mitochondrial processes. Furthermore, we characterized the impact of two recoding sites in the mitofusin 1 (MFN1) gene and showed their functional relevance to mitochondrial fusion and cellular apoptosis. Our study reveals a global reduction of editing in SCZ and a compelling link between editing and mitochondrial function in the disease.

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

RNA editing modifies single nucleotides of RNAs, regulating primary protein structure and protein abundance. In recent years, the diversity of proteins and complexity of gene regulation associated with RNA editing dysregulation has been increasingly appreciated in oncology. Large-scale shifts in editing have been observed in bulk tumors across various cancer types. However, RNA editing in single cells and individual cell types within tumors has not been explored. By profiling editing in single cells from lung adenocarcinoma biopsies, we found that the increased editing trend of bulk lung tumors was unique to cancer cells. Elevated editing levels were observed in cancer cells resistant to targeted therapy, and editing sites associated with drug response were enriched. Consistent with the regulation of antiviral pathways by RNA editing, higher editing levels in cancer cells were associated with reduced antitumor innate immune response, especially levels of natural killer cell infiltration. In addition, the level of RNA editing in cancer cells was positively associated with somatic point mutation burden. This observation motivated the definition of a new metric, RNA editing load, reflecting the amount of RNA mutations created by RNA editing. Importantly, in lung cancer, RNA editing load was a stronger predictor of patient survival than DNA mutations. This study provides the first single cell dissection of editing in cancer and highlights the significance of RNA editing load in cancer prognosis.

SIGNIFICANCE

RNA editing analysis in single lung adenocarcinoma cells uncovers RNA mutations that correlate with tumor mutation burden and cancer innate immunity and reveals the amount of RNA mutations that strongly predicts patient survival. See related commentary by Luo and Liang, p. 351.

Dynamic regulation and key roles of ribonucleic acid methylation

Ribonucleic acid (RNA) methylation is the most abundant modification in biological systems, accounting for 60% of all RNA modifications, and affects multiple aspects of RNA (including mRNAs, tRNAs, rRNAs, microRNAs, and long non-coding RNAs). Dysregulation of RNA methylation causes many developmental diseases through various mechanisms mediated by N ⁶-methyladenosine (m⁶A), 5-methylcytosine (m⁵C), N ¹-methyladenosine (m¹A), 5-hydroxymethylcytosine (hm⁵C), and pseudouridine (Ψ). The emerging tools of RNA methylation can be used as diagnostic, preventive, and therapeutic markers. Here, we review the accumulated discoveries to date regarding the biological function and dynamic regulation of RNA methylation/modification, as well as the most popularly used techniques applied for profiling RNA epitranscriptome, to provide new ideas for growth and development.

Spatiotemporal and genetic regulation of A-to-I editing throughout human brain development

Posttranscriptional RNA modifications by adenosine-to-inosine (A-to-I) editing are abundant in the brain, yet elucidating functional sites remains challenging. To bridge this gap, we investigate spatiotemporal and genetically regulated A-to-I editing sites across prenatal and postnatal stages of human brain development. More than 10,000 spatiotemporally regulated A-to-I sites were identified that occur predominately in 3' UTRs and introns, as well as 37 sites that recode amino acids in protein coding regions with precise changes in editing levels across development. Hyper-edited transcripts are also enriched in the aging brain and stabilize RNA secondary structures. These features are conserved in murine and non-human primate models of neurodevelopment. Finally, thousands of cis-editing quantitative trait loci (edQTLs) were identified with unique regulatory effects during prenatal and postnatal development. Collectively, this work offers a resolved atlas linking spatiotemporal variation in editing levels to genetic regulatory effects throughout distinct stages of brain maturation.

Chemical Probe-Based Nanopore Sequencing to Selectively Assess the RNA Modifications

Nanopore direct RNA sequencing (dRNA-Seq) reads reveal RNA modifications through consistent error profiles specific to a modified nucleobase. However, a null data set is required to identify actual RNA modification-associated errors for distinguishing it from confounding highly intrinsic sequencing errors. Here, we reveal that inosine creates a signature mismatch error in dRNA-Seq reads and obviates the need for a null data set by harnessing the selective reactivity of acrylonitrile for validating the presence of actual inosine modifications. Selective reactivity of acrylonitrile toward inosine altered multiple dRNA-Seq parameters like signal intensity and trace value. We also deduced the stoichiometry of inosine modification through deviation in signal intensity and trace value using this chemical biology approach. Furthermore, we devised Nano ICE-Seq, a protocol to overcome the low coverage issue associated with direct RNA sequencing. Taken together, our chemical probe-based approach may facilitate the knockout-free detection of disease-associated RNA modifications in clinical scenarios.

Multifaceted role of RNA editing in promoting loss-of-function of PODXL in cancer

PODXL, a protein that is dysregulated in multiple cancers, plays an important role in promoting cancer metastasis. In this study, we report that RNA editing promotes the inclusion of a PODXL alternative exon. The resulting edited PODXL long isoform is more prone to protease digestion and has the strongest effects on reducing cell migration and cisplatin chemoresistance among the three PODXL isoforms (short, unedited long, and edited long isoforms). Importantly, the editing level of the PODXL recoding site and the inclusion level of the PODXL alternative exon are strongly associated with overall patient survival in Kidney Renal Clear Cell Carcinoma (KIRC). Supported by significant enrichment of exonic RNA editing sites in alternatively spliced exons, we hypothesize that exonic RNA editing sites may enhance proteomic diversity through alternative splicing, in addition to amino acid changes, a previously under-appreciated aspect of RNA editing function.

ADAR1 Isoforms Regulate Let-7d Processing in Idiopathic Pulmonary Fibrosis

Double-stranded RNA adenosine deaminase 1 (ADAR1) is significantly down-regulated in fibroblasts derived from Idiopathic Pulmonary Fibrosis (IPF) patients, and its overexpression restored levels of miRNA-21, PELI1, and SPRY2. There are two ADAR1 isoforms in humans, ADAR1-p110 and ADAR1-p150, generated by an alternative promoter. Let-7d is considered an essential microRNA in Pulmonary Fibrosis (PF). In silico analysis revealed COL3A1 and SMAD2, proteins involved in the development of IPF, as Let-7d targets. We analyzed the role of ADAR1-p110 and ADAR1-p150 isoforms in the regulation of Let-7d maturation and the effect of this regulation on the expression of COL3A1 and SMAD2 in IPF fibroblast. We demonstrated that differential expression and subcellular distribution of ADAR1 isoforms in fibroblasts contribute to the up-regulation of pri-miR-Let-7d and down-regulation of mature Let-7d. Induction of overexpression of ADAR1 reestablishes the expression of pri-miR-Let-7d and Let-7d in lung fibroblasts. The reduction of mature Let-7d upregulates the expression of COL3A1 and SMAD2. Thus, ADAR1 isoforms and Let-7d could have a synergistic role in IPF, which is a promising explanation in the mechanisms of fibrosis development, and the regulation of both molecules could be used as a therapeutic approach in IPF.

Survey of the binding preferences of RNA-binding proteins to RNA editing events

Background

Adenosine-to-inosine (A-to-I) editing is an important RNA posttranscriptional process related to a multitude of cellular and molecular activities. However, systematic characterizations of whether and how the events of RNA editing are associated with the binding preferences of RNA sequences to RNA-binding proteins (RBPs) are still lacking.

Results

With the RNA-seq and RBP eCLIP-seq datasets from the ENCODE project, we quantitatively survey the binding preferences of 150 RBPs to RNA editing events, followed by experimental validations. Such analyses of the RBP-associated RNA editing at nucleotide resolution and genome-wide scale shed light on the involvement of RBPs specifically in RNA editing-related processes, such as RNA splicing, RNA secondary structures, RNA decay, and other posttranscriptional processes.

Conclusions

These results highlight the relevance of RNA editing in the functions of many RBPs and therefore serve as a resource for further characterization of the functional associations between various RNA editing events and RBPs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02741-8.

RNA and neuronal function: the importance of post-transcriptional regulation

The brain represents an organ with a particularly high diversity of genes that undergo post-transcriptional gene regulation through multiple mechanisms that affect RNA metabolism and, consequently, brain function. This vast regulatory process in the brain allows for a tight spatiotemporal control over protein expression, a necessary factor due to the unique morphologies of neurons. The numerous mechanisms of post-transcriptional regulation or translational control of gene expression in the brain include alternative splicing, RNA editing, mRNA stability and transport. A large number of trans-elements such as RNA-binding proteins and micro RNAs bind to specific cis-elements on transcripts to dictate the fate of mRNAs including its stability, localization, activation and degradation. Several trans-elements are exemplary regulators of translation, employing multiple cofactors and regulatory machinery so as to influence mRNA fate. Networks of regulatory trans-elements exert control over key neuronal processes such as neurogenesis, synaptic transmission and plasticity. Perturbations in these networks may directly or indirectly cause neuropsychiatric and neurodegenerative disorders. We will be reviewing multiple mechanisms of gene regulation by trans-elements occurring specifically in neurons.

Brain Epitranscriptomic Analysis Revealed Altered A-to-I RNA Editing in Septic Patients

Recent studies suggest that RNA editing is associated with impaired brain function and neurological and psychiatric disorders. However, the role of A-to-I RNA editing during sepsis-associated encephalopathy (SAE) remains unclear. In this study, we analyzed adenosine-to-inosine (A-to-I) RNA editing in postmortem brain tissues from septic patients and controls. A total of 3024 high-confidence A-to-I RNA editing sites were identified. In sepsis, there were fewer A-to-I RNA editing genes and editing sites than in controls. Among all A-to-I RNA editing sites, 42 genes showed significantly differential RNA editing, with 23 downregulated and 19 upregulated in sepsis compared to controls. Notably, more than 50% of these genes were highly expressed in the brain and potentially related to neurological diseases. Notably, cis-regulatory analysis showed that the level of RNA editing in six differentially edited genes was significantly correlated with the gene expression, including HAUS augmin-like complex subunit 2 (HAUS2), protein phosphatase 3 catalytic subunit beta (PPP3CB), hook microtubule tethering protein 3 (HOOK3), CUB and Sushi multiple domains 1 (CSMD1), methyltransferase-like 7A (METTL7A), and kinesin light chain 2 (KLC2). Furthermore, enrichment analysis showed that fewer gene functions and KEGG pathways were enriched by edited genes in sepsis compared to controls. These results revealed alteration of A-to-I RNA editing in the human brain associated with sepsis, thus providing an important basis for understanding its role in neuropathology in SAE.

Inosine and its methyl derivatives: Occurrence, biogenesis, and function in RNA

Inosine is one of the most common post-transcriptional modifications. Since its discovery, it has been noted for its ability to contribute to non-Watson-Crick interactions within RNA. Rapidly accumulating evidence points to the widespread generation of inosine through hydrolytic deamination of adenosine to inosine by different classes of adenosine deaminases. Three naturally occurring methyl derivatives of inosine, i.e., 1-methylinosine, 2'-O-methylinosine and 1,2'-O-dimethylinosine are currently reported in RNA modification databases. These modifications are expected to lead to changes in the structure, folding, dynamics, stability and functions of RNA. The importance of the modifications is indicated by the strong conservation of the modifying enzymes across organisms. The structure, binding and catalytic mechanism of the adenosine deaminases have been well-studied, but the underlying mechanism of the catalytic reaction is not very clear yet. Here we extensively review the existing data on the occurrence, biogenesis and functions of inosine and its methyl derivatives in RNA. We also included the structural and thermodynamic aspects of these modifications in our review to provide a detailed and integrated discussion on the consequences of A-to-I editing in RNA and the contribution of different structural and thermodynamic studies in understanding its role in RNA. We also highlight the importance of further studies for a better understanding of the mechanisms of the different classes of deamination reactions. Further investigation of the structural and thermodynamic consequences and functions of these modifications in RNA should provide more useful information about their role in different diseases.

RNA editing regulates lncRNA splicing in human early embryo development

RNA editing is a co- or post-transcriptional modification through which some cells can make discrete changes to specific nucleotide sequences within an RNA molecule after transcription. Previous studies found that RNA editing may be critically involved in cancer and aging. However, the function of RNA editing in human early embryo development is still unclear. In this study, through analyzing single cell RNA sequencing data, 36.7% RNA editing sites were found to have a have differential editing ratio among early embryo developmental stages, and there was a great reprogramming of RNA editing rates at the 8-cell stage, at which most of the differentially edited RNA editing sites (99.2%) had a decreased RNA editing rate. In addition, RNA editing was more likely to occur on RNA splicing sites during human early embryo development. Furthermore, long non-coding RNA (lncRNA) editing sites were found more likely to be on RNA splicing sites (odds ratio = 2.19, P = 1.37×10-8), while mRNA editing sites were less likely (odds ratio = 0.22, P = 8.38×10-46). Besides, we found that the RNA editing rate on lncRNA had a significantly higher correlation coefficient with the percentage spliced index (PSI) of lncRNA exons (R = 0.75, P = 4.90×10-16), which indicated that RNA editing may regulate lncRNA splicing during human early embryo development. Finally, functional analysis revealed that those RNA editing-regulated lncRNAs were enriched in signal transduction, the regulation of transcript expression, and the transmembrane transport of mitochondrial calcium ion. Overall, our study might provide a new insight into the mechanism of RNA editing on lncRNAs in human developmental biology and common birth defects.

Primate-specific retrotransposons and the evolution of circadian networks in the human brain

The circadian rhythm of the human brain is attuned to sleep-wake cycles that entail global alterations in neuronal excitability. This periodicity involves a highly coordinated regulation of gene expression. A growing number of studies are documenting a fascinating connection between primate-specific retrotransposons (Alu elements) and key epigenetic regulatory processes in the primate brain. Collectively, these studies indicate that Alu elements embedded in the human neuronal genome mediate post-transcriptional processes that unite human-specific neuroepigenetic landscapes and circadian rhythm. Here, we review evidence linking Alu retrotransposon-mediated posttranscriptional pathways to circadian gene expression. We hypothesize that Alu retrotransposons participate in the organization of circadian brain function through multidimensional neuroepigenetic pathways. We anticipate that these pathways are closely tied to the evolution of human cognition and their perturbation contributes to the manifestation of human-specific neurological diseases. Finally, we address current challenges and accompanying opportunities in studying primate- and human-specific transposable elements.

Inosine Substitutions in RNA Activate Latent G-Quadruplexes

It is well-accepted that gene expression is heavily influenced by RNA structure. For instance, stem-loops and G-quadruplexes (rG4s) are dynamic motifs in mRNAs that influence gene expression. Adenosine-to-inosine (A-to-I) editing is a common chemical modification of RNA which introduces a nucleobase that is iso-structural with guanine, thereby changing RNA base-pairing properties. Here, we provide biophysical, chemical, and biological evidence that A-to-I exchange can activate latent rG4s by filling incomplete G-quartets with inosine. We demonstrate the formation of inosine-containing rG4s (GI-quadruplexes) in vitro and verify their activity in cells. GI-quadruplexes adopt parallel topologies, stabilized by potassium ions. They exhibit moderately reduced thermal stability compared to conventional G-quadruplexes. To study inosine-induced structural changes in a naturally occurring RNA, we use a synthetic approach that enables site-specific inosine incorporation in long RNAs. In summary, RNA GI-quadruplexes are a previously unrecognized structural motif that may contribute to the regulation of gene expression in vivo.

RDDSVM: accurate prediction of A-to-I RNA editing sites from sequence using support vector machines

Adenosine to inosine (A-to-I) editing in RNA is involved in various biological processes like gene expression, alternative splicing, and mRNA degradation associated with carcinogenesis and various human diseases. Therefore, accurate identification of RNA editing sites in transcriptome is valuable for research and medicine. RNA-seq is very useful for the detection of RNA editing events in condition-specific cells. However, computational analysis methods of RNA-seq data have considerable false-positive risks due to mapping errors. In this study, we developed a simple machine learning method using support vector machines to train sequence and structure information derived from flanking sequences of experimentally verified A-to-I editing sites to predict new A-to-I editing sites in RNA. The highest performance results were obtained by the model that utilizes the composition of the triplet sequence elements in the flanking regions of the in A-to-I editing sites. Using this model, the SVM classifier also showed high performance on experimentally verified data providing a sensitivity of 92.8%, specificity of 77.1%, and accuracy of 90.2%. To compare the predictive capacity of our method with other classifiers that use sequence information, we have used validated human A-to-I RNA editing sites by Sanger sequencing. Out of 58 validated editing sites, our method recognized 53 of them correctly with an accuracy of 91.4% outperforming other classifiers. As to our knowledge, this is the first case of utilization of the composition of the triplet sequence elements neighboring A-to-I editing sites for the prediction of new A-to-I editing sites in RNA. The methodology is very easy to perform and computationally low demanding making it a convenient and valuable choice for facilities with low sources. To facilitate the usage of the method publicly, we developed an open-source program called RDDSVM to perform prediction on candidate A-to-I RNA editing sites using support vector machines.

ADeditome provides the genomic landscape of A-to-I RNA editing in Alzheimer's disease

A-to-I RNA editing, contributing to nearly 90% of all editing events in human, has been reported to involve in the pathogenesis of Alzheimer's disease (AD) due to its roles in brain development and immune regulation, such as the deficient editing of GluA2 Q/R related to cell death and memory loss. Currently, there are urgent needs for the systematic annotations of A-to-I RNA editing events in AD. Here, we built ADeditome, the annotation database of A-to-I RNA editing in AD available at https://ccsm.uth.edu/ADeditome, aiming to provide a resource and reference for functional annotation of A-to-I RNA editing in AD to identify therapeutically targetable genes in an individual. We detected 1676 363 editing sites in 1524 samples across nine brain regions from ROSMAP, MayoRNAseq and MSBB. For these editing events, we performed multiple functional annotations including identification of specific and disease stage associated editing events and the influence of editing events on gene expression, protein recoding, alternative splicing and miRNA regulation for all the genes, especially for AD-related genes in order to explore the pathology of AD. Combing all the analysis results, we found 108 010 and 26 168 editing events which may promote or inhibit AD progression, respectively. We also found 5582 brain region-specific editing events with potentially dual roles in AD across different brain regions. ADeditome will be a unique resource for AD and drug research communities to identify therapeutically targetable editing events. Significance: ADeditome is the first comprehensive resource of the functional genomics of individual A-to-I RNA editing events in AD, which will be useful for many researchers in the fields of AD pathology, precision medicine, and therapeutic researches.

MiREDiBase, a manually curated database of validated and putative editing events in microRNAs

MicroRNAs (miRNAs) are regulatory small non-coding RNAs that function as translational repressors. MiRNAs are involved in most cellular processes, and their expression and function are presided by several factors. Amongst, miRNA editing is an epitranscriptional modification that alters the original nucleotide sequence of selected miRNAs, possibly influencing their biogenesis and target-binding ability. A-to-I and C-to-U RNA editing are recognized as the canonical types, with the A-to-I type being the predominant one. Albeit some bioinformatics resources have been implemented to collect RNA editing data, it still lacks a comprehensive resource explicitly dedicated to miRNA editing. Here, we present MiREDiBase, a manually curated catalog of editing events in miRNAs. The current version includes 3,059 unique validated and putative editing sites from 626 pre-miRNAs in humans and three primates. Editing events in mature human miRNAs are supplied with miRNA-target predictions and enrichment analysis, while minimum free energy structures are inferred for edited pre-miRNAs. MiREDiBase represents a valuable tool for cell biology and biomedical research and will be continuously updated and expanded at https://ncrnaome.osumc.edu/miredibase .

Epitranscriptomics: A New Layer of microRNA Regulation in Cancer

MicroRNAs are pervasive regulators of gene expression at the post-transcriptional level in metazoan, playing key roles in several physiological and pathological processes. Accordingly, these small non-coding RNAs are also involved in cancer development and progression. Furthermore, miRNAs represent valuable diagnostic and prognostic biomarkers in malignancies. In the last twenty years, the role of RNA modifications in fine-tuning gene expressions at several levels has been unraveled. All RNA species may undergo post-transcriptional modifications, collectively referred to as epitranscriptomic modifications, which, in many instances, affect RNA molecule properties. miRNAs are not an exception, in this respect, and they have been shown to undergo several post-transcriptional modifications. In this review, we will summarize the recent findings concerning miRNA epitranscriptomic modifications, focusing on their potential role in cancer development and progression.

Regulatory Mechanisms of the RNA Modification m6A and Significance in Brain Function in Health and Disease

RNA modifications have emerged as an additional layer of regulatory complexity governing the function of almost all species of RNA. N⁶-methyladenosine (m⁶A), the addition of methyl groups to adenine residues, is the most abundant and well understood RNA modification. The current review discusses the regulatory mechanisms governing m⁶A, how this influences neuronal development and function and how aberrant m⁶A signaling may contribute to neurological disease. M⁶A is known to regulate the stability of mRNA, the processing of microRNAs and function/processing of tRNAs among other roles. The development of antibodies against m⁶A has facilitated the application of next generation sequencing to profile methylated RNAs in both health and disease contexts, revealing the extent of this transcriptomic modification. The mechanisms by which m⁶A is deposited, processed, and potentially removed are increasingly understood. Writer enzymes include METTL3 and METTL14 while YTHDC1 and YTHDF1 are key reader proteins, which recognize and bind the m⁶A mark. Finally, FTO and ALKBH5 have been identified as potential erasers of m⁶A, although there in vivo activity and the dynamic nature of this modification requires further study. M⁶A is enriched in the brain and has emerged as a key regulator of neuronal activity and function in processes including neurodevelopment, learning and memory, synaptic plasticity, and the stress response. Changes to m⁶A have recently been linked with Schizophrenia and Alzheimer disease. Elucidating the functional consequences of m⁶A changes in these and other brain diseases may lead to novel insight into disease pathomechanisms, molecular biomarkers and novel therapeutic targets.

Deciphering the principles of the RNA editing code via large-scale systematic probing

Adenosine-to-inosine editing is catalyzed by ADAR1 at thousands of sites transcriptome-wide. Despite intense interest in ADAR1 from physiological, bioengineering, and therapeutic perspectives, the rules of ADAR1 substrate selection are poorly understood. Here, we used large-scale systematic probing of ∼2,000 synthetic constructs to explore the structure and sequence context determining editability. We uncover two structural layers determining the formation and propagation of A-to-I editing, independent of sequence. First, editing is robustly induced at fixed intervals of 35 bp upstream and 30 bp downstream of structural disruptions. Second, editing is symmetrically introduced on opposite sites on a double-stranded structure. Our findings suggest a recursive model for RNA editing, whereby the structural alteration induced by the editing at one site iteratively gives rise to the formation of an additional editing site at a fixed periodicity, serving as a basis for the propagation of editing along and across both strands of double-stranded RNA structures.

Detecting and Characterizing A-To-I microRNA Editing in Cancer

Adenosine to inosine (A-to-I) editing consists of an RNA modification where single adenosines along the RNA sequence are converted into inosines. Such a biochemical transformation is catalyzed by enzymes belonging to the family of adenosine deaminases acting on RNA (ADARs) and occurs either co- or post-transcriptionally. The employment of powerful, high-throughput detection methods has recently revealed that A-to-I editing widely occurs in non-coding RNAs, including microRNAs (miRNAs). MiRNAs are a class of small regulatory non-coding RNAs (ncRNAs) acting as translation inhibitors, known to exert relevant roles in controlling cell cycle, proliferation, and cancer development. Indeed, a growing number of recent researches have evidenced the importance of miRNA editing in cancer biology by exploiting various detection and validation methods. Herein, we briefly overview early and currently available A-to-I miRNA editing detection and validation methods and discuss the significance of A-to-I miRNA editing in human cancer.

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.

RNA editing in cancer impacts mRNA abundance in immune response pathways

BACKGROUND

RNA editing generates modifications to the RNA sequences, thereby increasing protein diversity and shaping various layers of gene regulation. Recent studies have revealed global shifts in editing levels across many cancer types, as well as a few specific mechanisms implicating individual sites in tumorigenesis or metastasis. However, most tumor-associated sites, predominantly in noncoding regions, have unknown functional relevance.

RESULTS

Here, we carry out integrative analysis of RNA editing profiles between epithelial and mesenchymal tumors, since epithelial-mesenchymal transition is a key paradigm for metastasis. We identify distinct editing patterns between epithelial and mesenchymal tumors in seven cancer types using TCGA data, an observation further supported by single-cell RNA sequencing data and ADAR perturbation experiments in cell culture. Through computational analyses and experimental validations, we show that differential editing sites between epithelial and mesenchymal phenotypes function by regulating mRNA abundance of their respective genes. Our analysis of RNA-binding proteins reveals ILF3 as a potential regulator of this process, supported by experimental validations. Consistent with the known roles of ILF3 in immune response, epithelial-mesenchymal differential editing sites are enriched in genes involved in immune and viral processes. The strongest target of editing-dependent ILF3 regulation is the transcript encoding PKR, a crucial player in immune and viral response.

CONCLUSIONS

Our study reports widespread differences in RNA editing between epithelial and mesenchymal tumors and a novel mechanism of editing-dependent regulation of mRNA abundance. It reveals the broad impact of RNA editing in cancer and its relevance to cancer-related immune pathways.

Statistical inference of differential RNA-editing sites from RNA-sequencing data by hierarchical modeling

MOTIVATION

RNA-sequencing (RNA-seq) enables global identification of RNA-editing sites in biological systems and disease. A salient step in many studies is to identify editing sites that statistically associate with treatment (e.g. case versus control) or covary with biological factors, such as age. However, RNA-seq has technical features that incumbent tests (e.g. t-test and linear regression) do not consider, which can lead to false positives and false negatives.

RESULTS

In this study, we demonstrate the limitations of currently used tests and introduce the method, RNA-editing tests (REDITs), a suite of tests that employ beta-binomial models to identify differential RNA editing. The tests in REDITs have higher sensitivity than other tests, while also maintaining the type I error (false positive) rate at the nominal level. Applied to the GTEx dataset, we unveil RNA-editing changes associated with age and gender, and differential recoding profiles between brain regions.

AVAILABILITY AND IMPLEMENTATION

REDITs are implemented as functions in R and freely available for download at https://github.com/gxiaolab/REDITs. The repository also provides a code example for leveraging parallelization using multiple cores.

The development and controversy of competitive endogenous RNA hypothesis in non-coding genes

As a momentous post-transcriptional regulator, microRNAs (miRNAs) are attracting more and more attention. The classical miRNAs regulated mechanism shows it binds to the targets' 3'UTR thus play the role in post-transcription. Meanwhile, single miRNA can target multiple genes, so those should compete to bind that miRNA. Vice versa, single gene can sponge mass of miRNAs as well. Thus the competitive endogenous RNAs (ceRNAs) hypothesis was put forward in 2011. The ceRNA hypothesis has made huge achievements, in particular in non-coding genes, which including long non-coding RNAs (lncRNAs), circle RNAs (circRNAs) and pseudogenes, even viral transcripts. It also contributed greatly to epigenetics development. However, an increasing number of controversies have occurred with applause. Based on this situation, this review introduces something in detail about the ceRNAs hypothesis achieved in lncRNAs, circRNAs, pseudogenes and viral transcripts, respectively. Meanwhile, it also covers controversy of the ceRNAs hypothesis.

A-to-I RNA editing as a tuner of noncoding RNAs in cancer

Recent advancement in RNA technology and computation biology shows the abundance and impact of RNA editing at the genome-wide level. Of RNA editing events, Adenosine-to-inosine (A-to-I) RNA editing is one of the most frequent types of RNA editing catalyzed by ADAR proteins. Indeed, A-to-I RNA editing occurs at the various coding and noncoding regions, triggering abnormal signaling pathways involved in cancer pathogenesis. Noncoding RNAs such as microRNA and long noncoding RNA have emerged as key regulators of pathways in cancer. The RNA editing including A-to-I editing is enriched in noncoding regions because of the abundance of noncoding RNAs accounting for 99% of total transcripts in the human genome. The effects of A-to-I editing in coding genes have been investigated and reported. However, those in noncoding RNAs have been less known in spite of the high frequency of editing events in noncoding regions. In this review, we will briefly discuss current findings and potential directions of A-to-I RNA editing research of noncoding RNAs and cancer. We will also introduce the concept of A-to-I editing, ADAR proteins, RNA editing technologies and databases.

Identification of genetic variants controlling RNA editing and their effect on RNA structure stabilization

Post-transcriptional modification of RNA (RNA editing, RNAe) results in differences between the RNA transcript and the genomic DNA sequence (RDD). Enzymatic modification of adenosine to inosine (A2I) by ADAR is the most studied type of RNAe. However, few genetic association studies with A2I RNAe events have been conducted. Some studies have analyzed the inter-population RNAe-QTL diversity in humans, but the sample size of these studies was limited. Other types of RNA and DNA differences have been reported but are largely understudied. Here, we report a comprehensive analysis of all types of RDD, based on two independent datasets. We found that A2I was by far the most observed type of RDD. Moreover, manual curation suggests that A2I is likely the only enzymatically driven RNAe type observed in blood derived DNA, all other non-A2I RDD could either be attributed to sequencing and processing artifacts, or are a result of somatic DNA rearrangements. We then conducted an in-cis genetic association study and identified 472 genetic associations (RNAe-QTL), that were replicated in both datasets. We confirm the potential effect of the RNAe-QTL on RNA structure by showing that allele specific RNAe occurs in heterozygotes. Although the generally assumed function of RNAe is to destabilize double stranded RNA structure, we found clear evidence for the potential additional involvement of RNAe in maintaining RNA hairpin that has been altered by the RNAe-QTL. Our study confirms, in two independent datasets, the potential role of RNAe in maintaining RNA structure in the presence of genetic variation.

ADAR1-Dependent RNA Editing Promotes MET and iPSC Reprogramming by Alleviating ER Stress

RNA editing of adenosine to inosine (A to I) is catalyzed by ADAR1 and dramatically alters the cellular transcriptome, although its functional roles in somatic cell reprogramming are largely unexplored. Here, we show that loss of ADAR1-mediated A-to-I editing disrupts mesenchymal-to-epithelial transition (MET) during induced pluripotent stem cell (iPSC) reprogramming and impedes acquisition of induced pluripotency. Using chemical and genetic approaches, we show that absence of ADAR1-dependent RNA editing induces aberrant innate immune responses through the double-stranded RNA (dsRNA) sensor MDA5, unleashing endoplasmic reticulum (ER) stress and hindering epithelial fate acquisition. We found that A-to-I editing impedes MDA5 sensing and sequestration of dsRNAs encoding membrane proteins, which promote ER homeostasis by activating the PERK-dependent unfolded protein response pathway to consequently facilitate MET. This study therefore establishes a critical role for ADAR1 and its A-to-I editing activity during cell fate transitions and delineates a key regulatory layer underlying MET to control efficient reprogramming.

ADAR1 Facilitates KSHV Lytic Reactivation by Modulating the RLR-Dependent Signaling Pathway

Kaposi's sarcoma-associated herpesvirus (KSHV) is a double-stranded DNA virus that exhibits two alternative life cycles: latency and lytic reactivation. During lytic reactivation, host innate immune responses are activated to restrict viral replication. Here, we report that adenosine deaminase acting on RNA 1 (ADAR1) is required for optimal KSHV lytic reactivation from latency. Knockdown of ADAR1 in KSHV latently infected cells inhibits viral gene transcription and viral replication during KSHV lytic reactivation. ADAR1 deficiency also significantly increases type I interferon production during KSHV reactivation. This increased interferon response is dependent on activation of the RIG-I-like receptor (RLR) pathway. Depletion of ADAR1 together with either RIG-I, MDA5, or MAVS reverses the increased IFNβ production and rescues KSHV lytic replication. These data suggest that ADAR1 serves as a proviral factor for KSHV lytic reactivation and facilitates DNA virus reactivation by dampening the RLR pathway-mediated innate immune response.

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.

Genome-wide quantification of ADAR adenosine-to-inosine RNA editing activity

Adenosine-to-inosine (A-to-I) RNA editing by the adenosine deaminase that acts on RNA (ADAR) enzymes is a common RNA modification, preventing false activation of the innate immune system by endogenous double-stranded RNAs. Methods for quantification of ADAR activity are sought after, due to an increasing interest in the role of ADARs in cancer and autoimmune disorders, as well as attempts to harness the ADAR enzymes for RNA engineering. Here, we present the Alu editing index (AEI), a robust and simple-to-use computational tool devised for this purpose. We describe its properties and demonstrate its superiority to current quantification methods of ADAR activity. The AEI is used to map global editing across a large dataset of healthy human samples and identify putative regulators of ADAR, as well as previously unknown factors affecting the observed Alu editing levels. These should be taken into account in future comparative studies of ADAR activity. The AEI tool is available at https://github.com/a2iEditing/RNAEditingIndexer.

Epigenetic Components of Myalgic Encephalomyelitis/Chronic Fatigue Syndrome Uncover Potential Transposable Element Activation

PURPOSE

Studies to determine epigenetic changes associated with myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) remain scarce; however, current evidence clearly shows that methylation patterns of genomic DNA and noncoding RNA profiles of immune cells differ between patients and healthy subjects, suggesting an active role of these epigenetic mechanisms in the disease. The present study compares and contrasts the available ME/CFS epigenetic data in an effort to evidence overlapping pathways capable of explaining at least some of the dysfunctional immune parameters linked to this disease.

METHODS

A systematic search of the literature evaluating the ME/CFS epigenome landscape was performed following the Preferred Reporting Items for Systematic Reviews and Meta-Analyses criteria. Differential DNA methylation and noncoding RNA differential expression patterns associated with ME/CFS were used to screen for the presence of transposable elements using the Dfam browser, a search program nurtured with the Repbase repetitive sequence database and the RepeatMasker annotation tool.

FINDINGS

Unexpectedly, particular associations of transposable elements and ME/CFS epigenetic hallmarks were uncovered. A model for the disease emerged involving transcriptional induction of endogenous dormant transposons and structured cellular RNA interactions, triggering the activation of the innate immune system without a concomitant active infection.

IMPLICATIONS

Repetitive sequence filters (ie, RepeatMasker) should be avoided when analyzing transcriptomic data to assess the potential participation of repetitive sequences ("junk repetitive DNA"), representing >45% of the human genome, in the onset and evolution of ME/CFS. In addition, transposable element screenings aimed at designing cost-effective, focused empirical assays that can confirm or disprove the suspected involvement of transposon transcriptional activation in this disease, following the pilot strategy presented here, will require databases gathering large ME/CFS epigenetic datasets.

Genomic Positional Dissection of RNA Editomes in Tumor and Normal Samples

RNA editing is phenomenon that occurs in both protein coding and non-coding RNAs. Increasing evidence have shown that adenosine-to-inosine RNA editing can potentially rendering substantial functional effects throughout the genome. Using RNA editing datasets from two large consortiums: The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) project, we quantitatively analyzed human genome-wide RNA editing events derived from tumor or normal tissues. Generally, a common RNA editing site tends to have a higher editing level in tumors as compared to normal samples. Of the 14 tumor-normal-paired cancer types examined, Eleven of the 14 cancers tested had overall increased RNA editing levels in the tumors. The editomes in cancer or normal tissues were dissected by genomic locations, and significant RNA editing locational difference was found between cancerous and healthy subjects. Additionally, our results indicated a significant correlation between the RNA editing rate and the gene density across chromosomes, highlighted hyper RNA editing clusters through visualization of running RNA editing rates along chromosomes, and identified hyper RNA edited genes (protein-coding genes, lincRNAs, and pseudogenes) that embody a large portion of cancer prognostic predictors. This study reinforces the potential functional effects of RNA editing in protein-coding genes, and also makes a strong foundation for further exploration of RNA editing’s roles in non-coding regions.

Widespread cis-regulation of RNA editing in a large mammal

Post-transcriptional RNA editing may regulate transcript expression and diversity in cells, with potential impacts on various aspects of physiology and environmental adaptation. A small number of recent genome-wide studies in Drosophila, mouse, and human have shown that RNA editing can be genetically modulated, highlighting loci that quantitatively impact editing of transcripts. The potential gene expression and physiological consequences of these RNA-editing quantitative trait loci (edQTL), however, are almost entirely unknown. Here, we present analyses of RNA editing in a large domestic mammal (Bos taurus), where we use whole-genome and high-depth RNA sequencing to discover, characterize, and conduct genetic mapping studies of novel transcript edits. Using a discovery population of nine deeply sequenced cows, we identify 2413 edit sites in the mammary transcriptome, the majority of which are adenosine to inosine edits (98.6%). Most sites are predicted to reside in double-stranded secondary structures (85.1%), and quantification of the rates of editing in an additional 355 cows reveals editing is negatively correlated with gene expression in the majority of cases. Genetic analyses of RNA editing and gene expression highlight 152 cis-regulated edQTL, of which 15 appear to cosegregate with expression QTL effects. Trait association analyses in a separate population of 9989 lactating cows also shows 12 of the cis-edQTL coincide with at least one cosegregating lactation QTL. Together, these results enhance our understanding of RNA-editing dynamics in mammals, and suggest mechanistic links by which loci may impact phenotype through RNA editing mediated processes.

Evolutionarily significant A-to-I RNA editing events originated through G-to-A mutations in primates

BACKGROUND

Recent studies have revealed thousands of A-to-I RNA editing events in primates, but the origination and general functions of these events are not well addressed.

RESULTS

Here, we perform a comparative editome study in human and rhesus macaque and uncover a substantial proportion of macaque A-to-I editing sites that are genomically polymorphic in some animals or encoded as non-editable nucleotides in human. The occurrence of these recent gain and loss of RNA editing through DNA point mutation is significantly more prevalent than that expected for the nearby regions. Ancestral state analyses further demonstrate that an increase in recent gain of editing events contribute to the over-representation, with G-to-A mutation site as a favorable location for the origination of robust A-to-I editing events. Population genetics analyses of the focal editing sites further reveal that a portion of these young editing events are evolutionarily significant, indicating general functional relevance for at least a fraction of these sites.

CONCLUSIONS

Overall, we report a list of A-to-I editing events that recently originated through G-to-A mutations in primates, representing a valuable resource to investigate the features and evolutionary significance of A-to-I editing events at the population and species levels. The unique subset of primate editome also illuminates the general functions of RNA editing by connecting it to particular gene regulatory processes, based on the characterized outcome of a gene regulatory level in different individuals or primate species with or without these editing events.

Multiple QTL underlie milk phenotypes at the CSF2RB locus

Background

Over many years, artificial selection has substantially improved milk production by cows. However, the genes that underlie milk production quantitative trait loci (QTL) remain relatively poorly characterised. Here, we investigate a previously reported QTL located at the CSF2RB locus on chromosome 5, for several milk production phenotypes, to better understand its underlying genetic and molecular causes.

Results

Using a population of 29,350 taurine dairy cows, we conducted association analyses for milk yield and composition traits, and identified highly significant QTL for milk yield, milk fat concentration, and milk protein concentration. Strikingly, protein concentration and milk yield appear to show co-located yet genetically distinct QTL. To attempt to understand the molecular mechanisms that might be mediating these effects, gene expression data were used to investigate eQTL for 11 genes in the broader interval. This analysis highlighted genetic impacts on CSF2RB and NCF4 expression that share similar association signatures to those observed for lactation QTL, strongly implicating one or both of these genes as responsible for these effects. Using the same gene expression dataset representing 357 lactating cows, we also identified 38 novel RNA editing sites in the 3′ UTR of CSF2RB transcripts. The extent to which two of these sites were edited also appears to be genetically co-regulated with lactation QTL, highlighting a further layer of regulatory complexity that involves the CSF2RB gene.

Conclusions

This locus presents a diversity of molecular and lactation QTL, likely representing multiple overlapping effects that, at a minimum, highlight the CSF2RB gene as having a causal role in these processes.

Electronic supplementary material

The online version of this article (10.1186/s12711-019-0446-x) contains supplementary material, which is available to authorized users.

Regulation of RNA editing by RNA-binding proteins in human cells

Adenosine-to-inosine (A-to-I) editing, mediated by the ADAR enzymes, diversifies the transcriptome by altering RNA sequences. Recent studies reported global changes in RNA editing in disease and development. Such widespread editing variations necessitate an improved understanding of the regulatory mechanisms of RNA editing. Here, we study the roles of >200 RNA-binding proteins (RBPs) in mediating RNA editing in two human cell lines. Using RNA-sequencing and global protein-RNA binding data, we identify a number of RBPs as key regulators of A-to-I editing. These RBPs, such as TDP-43, DROSHA, NF45/90 and Ro60, mediate editing through various mechanisms including regulation of ADAR1 expression, interaction with ADAR1, and binding to Alu elements. We highlight that editing regulation by Ro60 is consistent with the global up-regulation of RNA editing in systemic lupus erythematosus. Additionally, most key editing regulators act in a cell type-specific manner. Together, our work provides insights for the regulatory mechanisms of RNA editing.

ncRNA Editing: Functional Characterization and Computational Resources

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