REDIportal: a comprehensive database of A-to-I RNA editing events in humans

Narrowing down the candidates of beneficial A-to-I RNA editing by comparing the recoding sites with uneditable counterparts

Adar-mediated adenosine-to-inosine (A-to-I) RNA editing mainly occurs in nucleus and diversifies the transcriptome in a flexible manner. It has been a challenging task to identify beneficial editing sites from the sea of total editing events. The functional Ser>Gly auto-recoding site in insect Adar gene has uneditable Ser codons in ancestral nodes, indicating the selective advantage to having an editable status. Here, we extended this case study to more metazoan species, and also looked for all Drosophila recoding events with potential uneditable synonymous codons. Interestingly, in D. melanogaster, the abundant nonsynonymous editing is enriched in the codons that have uneditable counterparts, but the Adar Ser>Gly case suggests that the editable orthologous codons in other species are not necessarily edited. The use of editable versus ancestral uneditable codon is a smart way to infer the selective advantage of RNA editing, and priority might be given to these editing sites for functional studies due to the feasibility to construct an uneditable allele. Our study proposes an idea to narrow down the candidates of beneficial recoding sites. Meanwhile, we stress that the matched transcriptomes are needed to verify the conservation of editing events during evolution.

Spatially exploring RNA biology in archival formalin-fixed paraffin-embedded tissues

The capability to spatially explore RNA biology in formalin-fixed paraffin-embedded (FFPE) tissues holds transformative potential for histopathology research. Here, we present pathology-compatible deterministic barcoding in tissue (Patho-DBiT) by combining in situ polyadenylation and computational innovation for spatial whole transcriptome sequencing, tailored to probe the diverse RNA species in clinically archived FFPE samples. It permits spatial co-profiling of gene expression and RNA processing, unveiling region-specific splicing isoforms, and high-sensitivity transcriptomic mapping of clinical tumor FFPE tissues stored for 5 years. Furthermore, genome-wide single-nucleotide RNA variants can be captured to distinguish malignant subclones from non-malignant cells in human lymphomas. Patho-DBiT also maps microRNA regulatory networks and RNA splicing dynamics, decoding their roles in spatial tumorigenesis. Single-cell level Patho-DBiT dissects the spatiotemporal cellular dynamics driving tumor clonal architecture and progression. Patho-DBiT stands poised as a valuable platform to unravel rich RNA biology in FFPE tissues to aid in clinical pathology evaluation.

Charting and probing the activity of ADARs in human development and cell-fate specification

Adenosine deaminases acting on RNA (ADARs) impact diverse cellular processes and pathological conditions, but their functions in early cell-fate specification remain less understood. To gain insights here, we began by charting time-course RNA editing profiles in human organs from fetal to adult stages. Next, we utilized hPSC differentiation to experimentally probe ADARs, harnessing brain organoids as neural specific, and teratomas as pan-tissue developmental models. We show that time-series teratomas faithfully recapitulate fetal developmental trends, and motivated by this, conducted pan-tissue, single-cell CRISPR-KO screens of ADARs in teratomas. Knocking out ADAR leads to a global decrease in RNA editing across all germ-layers. Intriguingly, knocking out ADAR leads to an enrichment of adipogenic cells, revealing a role for ADAR in human adipogenesis. Collectively, we present a multi-pronged framework charting time-resolved RNA editing profiles and coupled ADAR perturbations in developmental models, thereby shedding light on the role of ADARs in cell-fate specification.

Validation of the APOBEC3A-mediated RNA Single Base Substitution Signature and Proposal of Novel APOBEC1, APOBEC3B, and APOBEC3G RNA Signatures

Mutational signature analysis gained significant attention for providing critical insights into the underlying mutational processes for various DNA single base substitution (SBS) signatures and their associations with different cancer types. Recently, RNA single base substitution (RNA-SBS) signatures were defined and described by decomposing RNA variants found in non-small cell lung cancer. Through statistical association, they attributed Apolipoprotein B mRNA Editing Enzyme, Catalytic Polypeptide 3A (APOBEC3A) mutagenesis to the RNA-SBS2 signature. Here, we provide the first validation of an RNA-SBS mutational signature by decomposing novel exogenous and endogenous APOBEC3A RNA editing signatures into COSMICv3.4 RNA-SBS reference signatures. Additionally, we have identified novel RNA-SBS signatures for APOBEC1, APOBEC3B, and APOBEC3G.

LoDEI: a robust and sensitive tool to detect transcriptome-wide differential A-to-I editing in RNA-seq data

RNA editing is a highly conserved process. Adenosine deaminase acting on RNA (ADAR) mediated deamination of adenosine (A-to-I editing) is associated with human disease and immune checkpoint control. Functional implications of A-to-I editing are currently of broad interest to academic and industrial research as underscored by the fast-growing number of clinical studies applying base editors as therapeutic tools. Analyzing the dynamics of A-to-I editing, in a biological or therapeutic context, requires the sensitive detection of differential A-to-I editing, a currently unmet need. We introduce the local differential editing index (LoDEI) to detect differential A-to-I editing in RNA-seq datasets using a sliding-window approach coupled with an empirical q value calculation that detects more A-to-I editing sites at the same false-discovery rate compared to existing methods. LoDEI is validated on known and novel datasets revealing that the oncogene MYCN increases and that a specific small non-coding RNA reduces A-to-I editing.

Epitranscriptome in action: RNA modifications in the DNA damage response

Complex pathways involving the DNA damage response (DDR) contend with cell-intrinsic and -extrinsic sources of DNA damage. DDR mis-regulation results in genome instability that can contribute to aging and diseases including cancer and neurodegeneration. Recent studies have highlighted key roles for several RNA species in the DDR, including short RNAs and RNA/DNA hybrids (R-loops) at DNA break sites, all contributing to efficient DNA repair. RNAs can undergo more than 170 distinct chemical modifications. These RNA modifications have emerged as key orchestrators of the DDR. Here, we highlight the function of enzyme- and non-enzyme-induced RNA modifications in the DDR, with particular emphasis on m6A, m5C, and RNA editing. We also discuss stress-induced RNA damage, including RNA alkylation/oxidation, RNA-protein crosslinks, and UV-induced RNA damage. Uncovering molecular mechanisms that underpin the contribution of RNA modifications to DDR and genome stability will have direct application to disease and approaches for therapeutic intervention.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk postmortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA editing and alternative polyadenylation (APA) within a cell-type-specific population of human neural progenitors and neurons. More RNA editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting that genetically mediated post-transcriptional regulation during brain development leads to differences in brain function.

Inflammation primes the murine kidney for recovery by activating AZIN1 adenosine-to-inosine editing

The progression of kidney disease varies among individuals, but a general methodology to quantify disease timelines is lacking. Particularly challenging is the task of determining the potential for recovery from acute kidney injury following various insults. Here, we report that quantitation of post-transcriptional adenosine-to-inosine (A-to-I) RNA editing offers a distinct genome-wide signature, enabling the delineation of disease trajectories in the kidney. A well-defined murine model of endotoxemia permitted the identification of the origin and extent of A-to-I editing, along with temporally discrete signatures of double-stranded RNA stress and adenosine deaminase isoform switching. We found that A-to-I editing of antizyme inhibitor 1 (AZIN1), a positive regulator of polyamine biosynthesis, serves as a particularly useful temporal landmark during endotoxemia. Our data indicate that AZIN1 A-to-I editing, triggered by preceding inflammation, primes the kidney and activates endogenous recovery mechanisms. By comparing genetically modified human cell lines and mice locked in either A-to-I-edited or uneditable states, we uncovered that AZIN1 A-to-I editing not only enhances polyamine biosynthesis but also engages glycolysis and nicotinamide biosynthesis to drive the recovery phenotype. Our findings implicate that quantifying AZIN1 A-to-I editing could potentially identify individuals who have transitioned to an endogenous recovery phase. This phase would reflect their past inflammation and indicate their potential for future recovery.

Loss of ADAR1 protein induces changes in small RNA landscape in hepatocytes

In recent years, numerous evidence has been accumulated about the extent of A-to-I editing in human RNAs and the key role ADAR1 plays in the cellular editing machinery. It has been shown that A-to-I editing occurrence and frequency are tissue-specific and essential for some tissue development, such as the liver. To study the effect of ADAR1 function in hepatocytes, we have created Huh7.5 ADAR1 KO cell lines. Upon IFN treatment, the Huh7.5 ADAR1 KO cells show rapid arrest of growth and translation, from which they do not recover. We analyzed translatome changes by using a method based on sequencing of separate polysome profile RNA fractions. We found significant changes in the transcriptome and translatome of the Huh7.5 ADAR1 KO cells. The most prominent changes include negatively affected transcription by RNA polymerase III and the deregulation of snoRNA and Y RNA levels. Furthermore, we observed that ADAR1 KO polysomes are enriched in mRNAs coding for proteins pivotal in a wide range of biological processes such as RNA localization and RNA processing, whereas the unbound fraction is enriched mainly in mRNAs coding for ribosomal proteins and translational factors. This indicates that ADAR1 plays a more relevant role in small RNA metabolism and ribosome biogenesis.

RBPs: an RNA editor's choice

RNA-binding proteins (RBPs) play a key role in gene expression and post-transcriptional RNA regulation. As integral components of ribonucleoprotein complexes, RBPs are susceptible to genomic and RNA Editing derived amino acid substitutions, impacting functional interactions. This article explores the prevalent RNA Editing of RBPs, unravelling the complex interplay between RBPs and RNA Editing events. Emphasis is placed on their influence on single amino acid variants (SAAVs) and implications for disease development. The role of Proteogenomics in identifying SAAVs is briefly discussed, offering insights into the RBP landscape. RNA Editing within RBPs emerges as a promising target for precision medicine, reshaping our understanding of genetic and epigenetic variations in health and disease.

RNA Editing by ADAR Adenosine Deaminases in the Cell Models of CAG Repeat Expansion Diseases: Significant Effect of Differentiation from Stem Cells into Brain Organoids in the Absence of Substantial Influence of CAG Repeats on the Level of Editing

Expansion of CAG repeats in certain genes is a known cause of several neurodegenerative diseases, but exact mechanism behind this is not yet fully understood. It is believed that the double-stranded RNA regions formed by CAG repeats could be harmful to the cell. This study aimed to test the hypothesis that these RNA regions might potentially interfere with ADAR RNA editing enzymes, leading to the reduced A-to-I editing of RNA and activation of the interferon response. We studied induced pluripotent stem cells (iPSCs) derived from the patients with Huntington's disease or ataxia type 17, as well as midbrain organoids developed from these cells. A targeted panel for next-generation sequencing was used to assess editing in the specific RNA regions. Differentiation of iPSCs into brain organoids led to increase in the ADAR2 gene expression and decrease in the expression of protein inhibitors of RNA editing. As a result, there was increase in the editing of specific ADAR2 substrates, which allowed identification of differential substrates of ADAR isoforms. However, comparison of the pathology and control groups did not show differences in the editing levels among the iPSCs. Additionally, brain organoids with 42-46 CAG repeats did not exhibit global changes. On the other hand, brain organoids with the highest number of CAG repeats in the huntingtin gene (76) showed significant decrease in the level of RNA editing of specific transcripts, potentially involving ADAR1. Notably, editing of the long non-coding RNA PWAR5 was nearly absent in this sample. It could be stated in conclusion that in most cultures with repeat expansion, the hypothesized effect on RNA editing was not confirmed.

Reovirus infection induces transcriptome-wide unique A-to-I editing changes in the murine fibroblasts

The conversion of Adenosine (A) to Inosine (I), by Adenosine Deaminases Acting on RNA or ADARs, is an essential post-transcriptional modification that contributes to proteome diversity and regulation in metazoans including humans. In addition to its transcriptome-regulating role, ADARs also play a major part in immune response to viral infection, where an interferon response activates interferon-stimulated genes, such as ADARp150, in turn dynamically regulating host-virus interactions. A previous report has shown that infection from reoviruses, despite strong activation of ADARp150, does not influence the editing of some of the major known editing targets, while likely editing others, suggesting a potentially nuanced editing pattern that may depend on different factors. However, the results were based on a handful of selected editing sites and did not cover the entire transcriptome. Thus, to determine whether and how reovirus infection specifically affects host ADAR editing patterns, we analyzed a publicly available deep-sequenced RNA-seq dataset, from murine fibroblasts infected with wild-type and mutant reovirus strains that allowed us to examine changes in editing patterns on a transcriptome-wide scale. To the best of our knowledge, this is the first transcriptome-wide report on host editing changes after reovirus infection. Our results demonstrate that reovirus infection induces unique nuanced editing changes in the host, including introducing sites uniquely edited in infected samples. Genes with edited sites are overrepresented in pathways related to immune regulation, cellular signaling, metabolism, and growth. Moreover, a shift in editing targets has also been observed, where the same genes are edited in infection and control conditions but at different sites, or where the editing rate is increased for some and decreased for other differential targets, supporting the hypothesis of dynamic and condition-specific editing by ADARs.

Detection of transcriptome-wide microRNA-target interactions in single cells with agoTRIBE

MicroRNAs (miRNAs) exert their gene regulatory effects on numerous biological processes based on their selection of target transcripts. Current experimental methods available to identify miRNA targets are laborious and require millions of cells. Here we have overcome these limitations by fusing the miRNA effector protein Argonaute2 to the RNA editing domain of ADAR2, allowing the detection of miRNA targets transcriptome-wide in single cells. miRNAs guide the fusion protein to their natural target transcripts, causing them to undergo A>I editing, which can be detected by sensitive single-cell RNA sequencing. We show that agoTRIBE identifies functional miRNA targets, which are supported by evolutionary sequence conservation. In one application of the method we study microRNA interactions in single cells and identify substantial differential targeting across the cell cycle. AgoTRIBE also provides transcriptome-wide measurements of RNA abundance and allows the deconvolution of miRNA targeting in complex tissues at the single-cell level.

Bioinformatics for Inosine: Tools and Approaches to Trace This Elusive RNA Modification

Inosine is a nucleotide resulting from the deamination of adenosine in RNA. This chemical modification process, known as RNA editing, is typically mediated by a family of double-stranded RNA binding proteins named Adenosine Deaminase Acting on dsRNA (ADAR). While the presence of ADAR orthologs has been traced throughout the evolution of metazoans, the existence and extension of RNA editing have been characterized in a more limited number of animals so far. Undoubtedly, ADAR-mediated RNA editing plays a vital role in physiology, organismal development and disease, making the understanding of the evolutionary conservation of this phenomenon pivotal to a deep characterization of relevant biological processes. However, the lack of direct high-throughput methods to reveal RNA modifications at single nucleotide resolution limited an extended investigation of RNA editing. Nowadays, these methods have been developed, and appropriate bioinformatic pipelines are required to fully exploit this data, which can complement existing approaches to detect ADAR editing. Here, we review the current literature on the "bioinformatics for inosine" subject and we discuss future research avenues in the field.

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

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

Precise in vivo RNA base editing with a wobble-enhanced circular CLUSTER guide RNA

Recruiting the endogenous editing enzyme adenosine deaminase acting on RNA (ADAR) with tailored guide RNAs for adenosine-to-inosine (A-to-I) RNA base editing is promising for safely manipulating genetic information at the RNA level. However, the precision and efficiency of editing are often compromised by bystander off-target editing. Here, we find that in 5'-UAN triplets, which dominate bystander editing, G•U wobble base pairs effectively mitigate off-target events while maintaining high on-target efficiency. This strategy is universally applicable to existing A-to-I RNA base-editing systems and complements other suppression methods such as G•A mismatches and uridine (U) depletion. Combining wobble base pairing with a circularized format of the CLUSTER approach achieves highly precise and efficient editing (up to 87%) of a disease-relevant mutation in the Mecp2 transcript in cell culture. Virus-mediated delivery of the guide RNA alone realizes functional MeCP2 protein restoration in the central nervous system of a murine Rett syndrome model with editing yields of up to 19% and excellent bystander control in vivo.

CanCellVar: A database for single-cell variants map in human cancer

Numerous variants, including both single-nucleotide variants (SNVs) in DNA and A>G RNA edits in mRNA as essential drivers of cellular proliferation and tumorigenesis, are commonly associated with cancer progression and growth. Thus, mining and summarizing single-cell variants will provide a refined and higher-resolution view of cancer and further contribute to precision medicine. Here, we established a database, CanCellVar, which aims to provide and visualize the comprehensive atlas of single-cell variants in tumor microenvironment. The current CanCellVar identified ∼3 million variants (∼1.4 million SNVs and ∼1.4 million A>G RNA edits) involved in 2,754,531 cells of 5 major cell types across 37 cancer types. CanCellVar provides the basic annotation information as well as cellular and molecular function properties of variants. In addition, the clinical relevance of variants can be obtained including tumor grade, treatment, metastasis, and others. Several flexible tools were also developed to aid retrieval and to analyze cell-cell interactions, gene expression, cell-development trajectories, regulation, and molecular structure affected by variants. Collectively, CanCellVar will serve as a valuable resource for investigating the functions and characteristics of single-cell variations and their roles in human tumor evolution and treatment.

A-to-I RNA co-editing predicts clinical outcomes and is associated with immune cells infiltration in hepatocellular carcinoma

Aberrant RNA editing has emerged as a pivotal factor in the pathogenesis of hepatocellular carcinoma (HCC), but the impact of RNA co-editing within HCC remains underexplored. We used a multi-step algorithm to construct an RNA co-editing network in HCC, and found that HCC-related RNA editings are predominantly centralized within the network. Furthermore, five pairs of risk RNA co-editing events were significantly correlated with the overall survival in HCC. Based on presence of risk RNA co-editings resulted in the categorization of HCC patients into high-risk and low-risk groups. Disparities in immune cell infiltrations were observed between the two groups, with the high-risk group exhibiting a greater abundance of exhausted T cells. Additionally, seven genes associated with risk RNA co-editing pairs were identified, whose expression effectively differentiates HCC tumor samples from normal ones. Our research offers an innovative perspective on the etiology and potential therapeutics for HCC.

Conquering new grounds: plant organellar C-to-U RNA editing factors can be functional in the plant cytosol

Plant mitochondrial and chloroplast transcripts are subject to numerous events of specific cytidine-to-uridine (C-to-U) RNA editing to correct genetic information. Key protein factors for this process are specific RNA-binding pentatricopeptide repeat (PPR) proteins, which are encoded in the nucleus and post-translationally imported into the two endosymbiotic organelles. Despite hundreds of C-to-U editing sites in the plant organelles, no comparable editing has been found for nucleo-cytosolic mRNAs raising the question why plant RNA editing is restricted to chloroplasts and mitochondria. Here, we addressed this issue in the model moss Physcomitrium patens, where all PPR-type RNA editing factors comprise specific RNA-binding and cytidine deamination functionalities in single proteins. To explore whether organelle-type RNA editing can principally also take place in the plant cytosol, we expressed PPR56, PPR65 and PPR78, three editing factors recently shown to also function in a bacterial setup, together with cytosolic co-transcribed native targets in Physcomitrium. While we obtained unsatisfying results upon their constitutive expression, we found strong cytosolic RNA editing under hormone-inducible expression. Moreover, RNA-Seq analyses revealed varying numbers of up to more than 900 off-targets in other cytosolic transcripts. We conclude that PPR-mediated C-to-U RNA editing is not per se incompatible with the plant cytosol but that its limited target specificity has restricted its occurrence to the much less complex transcriptomes of mitochondria and chloroplast in the course of evolution.

A-to-I Editing Is Subtype-Specific in Non-Hodgkin Lymphomas

Cancer is a complex and heterogeneous disease, in which a number of genetic and epigenetic changes occur in tumor onset and progression. Recent studies indicate that changes at the RNA level are also involved in tumorigenesis, such as adenosine-to-inosine (A-to-I) RNA editing. Here, we systematically investigate transcriptome-wide A-to-I editing events in a large number of samples from Non-Hodgkin lymphomas (NHLs). Using a computational pipeline that determines significant differences in editing level between NHL and normal samples at known A-to-I editing sites, we identify a number of differentially edited editing sites between NHL subtypes and normal samples. Most of the differentially edited sites are located in non-coding regions, and many such sites show a strong correlation between gene expression level and editing efficiency, indicating that RNA editing might have direct consequences for the cancer cell's aberrant gene regulation status in these cases. Moreover, we establish a strong link between RNA editing and NHL by demonstrating that NHL and normal samples and even NHL subtypes can be distinguished based on genome-wide RNA editing profiles alone. Our study establishes a strong link between RNA editing, cancer and aberrant gene regulation in NHL.

Divergent landscapes of A-to-I editing in postmortem and living human brain

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries offer more nuanced and accurate insights into the regulatory mechanisms of RNA editing in the human brain.

Exploring functional conservation in silico: a new machine learning approach to RNA-editing

Around 50 years ago, molecular biology opened the path to understand changes in forms, adaptations, complexity, or the basis of human diseases through myriads of reports on gene birth, gene duplication, gene expression regulation, and splicing regulation, among other relevant mechanisms behind gene function. Here, with the advent of big data and artificial intelligence (AI), we focus on an elusive and intriguing mechanism of gene function regulation, RNA editing, in which a single nucleotide from an RNA molecule is changed, with a remarkable impact in the increase of the complexity of the transcriptome and proteome. We present a new generation approach to assess the functional conservation of the RNA-editing targeting mechanism using two AI learning algorithms, random forest (RF) and bidirectional long short-term memory (biLSTM) neural networks with an attention layer. These algorithms, combined with RNA-editing data coming from databases and variant calling from same-individual RNA and DNA-seq experiments from different species, allowed us to predict RNA-editing events using both primary sequence and secondary structure. Then, we devised a method for assessing conservation or divergence in the molecular mechanisms of editing completely in silico: the cross-testing analysis. This novel method not only helps to understand the conservation of the editing mechanism through evolution but could set the basis for achieving a better understanding of the adenosine-targeting mechanism in other fields.

Pan-cancer discovery of somatic mutations from RNA sequencing data

Identification of somatic mutations (SMs) is essential for characterizing cancer genomes. While DNA-seq is the prevalent method for identifying SMs, RNA-seq provides an alternative strategy to discover tumor mutations in the transcribed genome. Here, we have developed a machine learning based pipeline to discover SMs based on RNA-seq data (designated as RNA-SMs). Subsequently, we have conducted a pan-cancer analysis to systematically identify RNA-SMs from over 8,000 tumors in The Cancer Genome Atlas (TCGA). In this way, we have identified over 105,000 novel SMs that had not been reported in previous TCGA studies. These novel SMs have significant clinical implications in designing targeted therapy for improved patient outcomes. Further, we have combined the SMs identified by both RNA-seq and DNA-seq analyses to depict an updated mutational landscape across 32 cancer types. This new online SM atlas, OncoDB ( https://oncodb.org ), offers a more complete view of gene mutations that underline the development and progression of various cancers.

Deep Learning for Elucidating Modifications to RNA-Status and Challenges Ahead

RNA-binding proteins and chemical modifications to RNA play vital roles in the co- and post-transcriptional regulation of genes. In order to fully decipher their biological roles, it is an essential task to catalogue their precise target locations along with their preferred contexts and sequence-based determinants. Recently, deep learning approaches have significantly advanced in this field. These methods can predict the presence or absence of modification at specific genomic regions based on diverse features, particularly sequence and secondary structure, allowing us to decipher the highly non-linear sequence patterns and structures that underlie site preferences. This article provides an overview of how deep learning is being applied to this area, with a particular focus on the problem of mRNA-RBP binding, while also considering other types of chemical modification to RNA. It discusses how different types of model can handle sequence-based and/or secondary-structure-based inputs, the process of model training, including choice of negative regions and separating sets for testing and training, and offers recommendations for developing biologically relevant models. Finally, it highlights four key areas that are crucial for advancing the field.

Divergent landscapes of A-to-I editing in postmortem and living human brain

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR1 and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries illuminate the nuanced functions and intricate regulatory mechanisms of RNA editing within the human brain.

RNA modifications in pulmonary diseases

Threatening public health, pulmonary disease (PD) encompasses diverse lung injuries like chronic obstructive PD, pulmonary fibrosis, asthma, pulmonary infections due to pathogen invasion, and fatal lung cancer. The crucial involvement of RNA epigenetic modifications in PD pathogenesis is underscored by robust evidence. These modifications not only shape cell fates but also finely modulate the expression of genes linked to disease progression, suggesting their utility as biomarkers and targets for therapeutic strategies. The critical RNA modifications implicated in PDs are summarized in this review, including N6-methylation of adenosine, N1-methylation of adenosine, 5-methylcytosine, pseudouridine (5-ribosyl uracil), 7-methylguanosine, and adenosine to inosine editing, along with relevant regulatory mechanisms. By shedding light on the pathology of PDs, these summaries could spur the identification of new biomarkers and therapeutic strategies, ultimately paving the way for early PD diagnosis and treatment innovation.

Integrated multiplexed assays of variant effect reveal determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase or decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

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.

Identification of skewed X chromosome inactivation using exome and transcriptome sequencing in patients with suspected rare genetic disease

BACKGROUND

X-chromosome inactivation (XCI) is an epigenetic process that occurs during early development in mammalian females by randomly silencing one of two copies of the X chromosome in each cell. The preferential inactivation of either the maternal or paternal copy of the X chromosome in a majority of cells results in a skewed or non-random pattern of X inactivation and is observed in over 25% of adult females. Identifying skewed X inactivation is of clinical significance in patients with suspected rare genetic diseases due to the possibility of biased expression of disease-causing genes present on the active X chromosome. The current clinical test for the detection of skewed XCI relies on the methylation status of the methylation-sensitive restriction enzyme (Hpall) binding site present in proximity of short tandem polymorphic repeats on the androgen receptor (AR) gene. This approach using one locus results in uninformative or inconclusive data for 10-20% of tests. Further, recent studies have shown inconsistency between methylation of the AR locus and the state of inactivation of the X chromosome. Herein, we develop a method for estimating X inactivation status, using exome and transcriptome sequencing data derived from blood in 227 female samples. We built a reference model for evaluation of XCI in 135 females from the GTEx consortium. We tested and validated the model on 11 female individuals with different types of undiagnosed rare genetic disorders who were clinically tested for X-skew using the AR gene assay and compared results to our outlier-based analysis technique.

RESULTS

In comparison to the AR clinical test for identification of X inactivation, our method was concordant with the AR method in 9 samples, discordant in 1, and provided a measure of X inactivation in 1 sample with uninformative clinical results. We applied this method on an additional 81 females presenting to the clinic with phenotypes consistent with different hereditary disorders without a known genetic diagnosis.

CONCLUSIONS

This study presents the use of transcriptome and exome sequencing data to provide an accurate and complete estimation of X-inactivation and skew status in a cohort of female patients with different types of suspected rare genetic disease.

RNA editing enzymes: structure, biological functions and applications

With the advancement of sequencing technologies and bioinformatics, over than 170 different RNA modifications have been identified. However, only a few of these modifications can lead to base pair changes, which are called RNA editing. RNA editing is a ubiquitous modification in mammalian transcriptomes and is an important co/posttranscriptional modification that plays a crucial role in various cellular processes. There are two main types of RNA editing events: adenosine to inosine (A-to-I) editing, catalyzed by ADARs on double-stranded RNA or ADATs on tRNA, and cytosine to uridine (C-to-U) editing catalyzed by APOBECs. This article provides an overview of the structure, function, and applications of RNA editing enzymes. We discuss the structural characteristics of three RNA editing enzyme families and their catalytic mechanisms in RNA editing. We also explain the biological role of RNA editing, particularly in innate immunity, cancer biogenesis, and antiviral activity. Additionally, this article describes RNA editing tools for manipulating RNA to correct disease-causing mutations, as well as the potential applications of RNA editing enzymes in the field of biotechnology and therapy.

Recommendations for detection, validation, and evaluation of RNA editing events in cardiovascular and neurological/neurodegenerative diseases

RNA editing, a common and potentially highly functional form of RNA modification, encompasses two different RNA modifications, namely adenosine to inosine (A-to-I) and cytidine to uridine (C-to-U) editing. As inosines are interpreted as guanosines by the cellular machinery, both A-to-I and C-to-U editing change the nucleotide sequence of the RNA. Editing events in coding sequences have the potential to change the amino acid sequence of proteins, whereas editing events in noncoding RNAs can, for example, affect microRNA target binding. With advancing RNA sequencing technology, more RNA editing events are being discovered, studied, and reported. However, RNA editing events are still often overlooked or discarded as sequence read quality defects. With this position paper, we aim to provide guidelines and recommendations for the detection, validation, and follow-up experiments to study RNA editing, taking examples from the fields of cardiovascular and brain disease. We discuss all steps, from sample collection, storage, and preparation, to different strategies for RNA sequencing and editing-sensitive data analysis strategies, to validation and follow-up experiments, as well as potential pitfalls and gaps in the available technologies. This paper may be used as an experimental guideline for RNA editing studies in any disease context.

In search of critical dsRNA targets of ADAR1

Recent studies have underscored the pivotal role of adenosine-to-inosine RNA editing, catalyzed by ADAR1, in suppressing innate immune interferon responses triggered by cellular double-stranded RNA (dsRNA). However, the specific ADAR1 editing targets crucial for this regulatory function remain elusive. We review analyses of transcriptome-wide ADAR1 editing patterns and their evolutionary dynamics, which offer valuable insights into this unresolved query. The growing appreciation of the significance of immunogenic dsRNAs and their editing in inflammatory and autoimmune diseases and cancer calls for a more comprehensive understanding of dsRNA immunogenicity, which may promote our understanding of these diseases and open doors to therapeutic avenues.

Spatially Exploring RNA Biology in Archival Formalin-Fixed Paraffin-Embedded Tissues

Spatial transcriptomics has emerged as a powerful tool for dissecting spatial cellular heterogeneity but as of today is largely limited to gene expression analysis. Yet, the life of RNA molecules is multifaceted and dynamic, requiring spatial profiling of different RNA species throughout the life cycle to delve into the intricate RNA biology in complex tissues. Human disease-relevant tissues are commonly preserved as formalin-fixed and paraffin-embedded (FFPE) blocks, representing an important resource for human tissue specimens. The capability to spatially explore RNA biology in FFPE tissues holds transformative potential for human biology research and clinical histopathology. Here, we present Patho-DBiT combining in situ polyadenylation and deterministic barcoding for spatial full coverage transcriptome sequencing, tailored for probing the diverse landscape of RNA species even in clinically archived FFPE samples. It permits spatial co-profiling of gene expression and RNA processing, unveiling region-specific splicing isoforms, and high-sensitivity transcriptomic mapping of clinical tumor FFPE tissues stored for five years. Furthermore, genome-wide single nucleotide RNA variants can be captured to distinguish different malignant clones from non-malignant cells in human lymphomas. Patho-DBiT also maps microRNA-mRNA regulatory networks and RNA splicing dynamics, decoding their roles in spatial tumorigenesis trajectory. High resolution Patho-DBiT at the cellular level reveals a spatial neighborhood and traces the spatiotemporal kinetics driving tumor progression. Patho-DBiT stands poised as a valuable platform to unravel rich RNA biology in FFPE tissues to study human tissue biology and aid in clinical pathology evaluation.

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.

Transcriptomic analysis reveals associations of blood-based A-to-I editing with Parkinson's disease

BACKGROUND

Adenosine-to-inosine (A-to-I) editing is the most common type of RNA editing in humans and the role of A-to-I RNA editing remains unclear in Parkinson's disease (PD).

OBJECTIVE

We aimed to explore the potential causal association between A-to-I editing and PD, and to assess whether changes in A-to-I editing were associated with cognitive progression in PD.

METHODS

The RNA-seq data from 380 PD patients and 178 healthy controls in the Parkinson's Progression Marker Initiative cohort was used to quantify A-to-I editing sites. We performed cis-RNA editing quantitative trait loci analysis and a two-sample Mendelian Randomization (MR) study by integrating genome-wide association studies to infer the potential causality between A-to-I editing and PD pathogenesis. The potential causal A-to-I editing sites were further confirmed by Summary-data-based MR analysis. Spearman's correlation analysis was performed to characterize the association between longitudinal A-to-I editing and cognitive progression in patients with PD.

RESULTS

We identified 17 potential causal A-to-I editing sites for PD and indicated that genetic risk variants may contribute to the risk of PD through A-to-I editing. These A-to-I editing sites were located in genes NCOR1, KANSL1 and BST1. Moreover, we observed 57 sites whose longitudinal A-to-I editing levels correlated with cognitive progression in PD.

CONCLUSIONS

We found potential causal A-to-I editing sites for PD onset and longitudinal changes of A-to-I editing were associated with cognitive progression in PD. We anticipate this study will provide new biological insights and drive the discovery of the epitranscriptomic role underlying Parkinson's disease.

The role of ADAR1 through and beyond its editing activity in cancer

Adenosine-to-inosine (A-to-I) editing of RNA, catalyzed by adenosine deaminase acting on RNA (ADAR) enzymes, is a prevalent RNA modification in mammals. It has been shown that A-to-I editing plays a critical role in multiple diseases, such as cardiovascular disease, neurological disorder, and particularly cancer. ADARs are the family of enzymes, including ADAR1, ADAR2, and ADAR3, that catalyze the occurrence of A-to-I editing. Notably, A-to-I editing is mainly catalyzed by ADAR1. Given the significance of A-to-I editing in disease development, it is important to unravel the complex roles of ADAR1 in cancer for the development of novel therapeutic interventions.In this review, we briefly describe the progress of research on A-to-I editing and ADARs in cancer, mainly focusing on the role of ADAR1 in cancer from both editing-dependent and independent perspectives. In addition, we also summarized the factors affecting the expression and editing activity of ADAR1 in cancer.

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

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

Unraveling C-to-U RNA editing events from direct RNA sequencing

In mammals, RNA editing events involve the conversion of adenosine (A) in inosine (I) by ADAR enzymes or the hydrolytic deamination of cytosine (C) in uracil (U) by the APOBEC family of enzymes, mostly APOBEC1. RNA editing has a plethora of biological functions, and its deregulation has been associated with various human disorders. While the large-scale detection of A-to-I is quite straightforward using the Illumina RNAseq technology, the identification of C-to-U events is a non-trivial task. This difficulty arises from the rarity of such events in eukaryotic genomes and the challenge of distinguishing them from background noise. Direct RNA sequencing by Oxford Nanopore Technology (ONT) permits the direct detection of Us on sequenced RNA reads. Surprisingly, using ONT reads from wild-type (WT) and APOBEC1-knock-out (KO) murine cell lines as well as in vitro synthesized RNA without any modification, we identified a systematic error affecting the accuracy of the Cs call, thereby leading to incorrect identifications of C-to-U events. To overcome this issue in direct RNA reads, here we introduce a novel machine learning strategy based on the isolation Forest (iForest) algorithm in which C-to-U editing events are considered as sequencing anomalies. Using in vitro synthesized and human ONT reads, our model optimizes the signal-to-noise ratio improving the detection of C-to-U editing sites with high accuracy, over 90% in all samples tested. Our results suggest that iForest, known for its rapid implementation and minimal memory requirements, is a promising tool to denoise ONT reads and reliably identify RNA modifications.

Harnessing ADAR-Mediated Site-Specific RNA Editing in Immune-Related Disease: Prediction and Therapeutic Implications

ADAR (Adenosine Deaminases Acting on RNA) proteins are a group of enzymes that play a vital role in RNA editing by converting adenosine to inosine in RNAs. This process is a frequent post-transcriptional event observed in metazoan transcripts. Recent studies indicate widespread dysregulation of ADAR-mediated RNA editing across many immune-related diseases, such as human cancer. We comprehensively review ADARs' function as pattern recognizers and their capability to contribute to mediating immune-related pathways. We also highlight the potential role of site-specific RNA editing in maintaining homeostasis and its relationship to various diseases, such as human cancers. More importantly, we summarize the latest cutting-edge computational approaches and data resources for predicting and analyzing RNA editing sites. Lastly, we cover the recent advancement in site-directed ADAR editing tool development. This review presents an up-to-date overview of ADAR-mediated RNA editing, how site-specific RNA editing could potentially impact disease pathology, and how they could be harnessed for therapeutic applications.

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

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

Elevated A-to-I RNA editing in COVID-19 infected individuals

Given the current status of coronavirus disease 2019 (COVID-19) as a global pandemic, it is of high priority to gain a deeper understanding of the disease's development and how the virus impacts its host. Adenosine (A)-to-Inosine (I) RNA editing is a post-transcriptional modification, catalyzed by the ADAR family of enzymes, that can be considered part of the inherent cellular defense mechanism as it affects the innate immune response in a complex manner. It was previously reported that various viruses could interact with the host's ADAR enzymes, resulting in epigenetic changes both to the virus and the host. Here, we analyze RNA-seq of nasopharyngeal swab specimens as well as whole-blood samples of COVID-19 infected individuals and show a significant elevation in the global RNA editing activity in COVID-19 compared to healthy controls. We also detect specific coding sites that exhibit higher editing activity. We further show that the increment in editing activity during the disease is temporary and returns to baseline shortly after the symptomatic period. These significant epigenetic changes may contribute to the immune system response and affect adverse outcomes seen in post-viral cases.

Cell-free multi-omics analysis reveals potential biomarkers in gastrointestinal cancer patients' blood

During cancer progression, tumorigenic and immune signals are spread through circulating molecules, such as cell-free DNA (cfDNA) and cell-free RNA (cfRNA) in the blood. So far, they have not been comprehensively investigated in gastrointestinal cancers. Here, we profile 4 categories of cell-free omics data from patients with colorectal cancer and patients with stomach adenocarcinoma and then assay 15 types of genomic, epigenomic, and transcriptomic variations. We find that multi-omics data are more appropriate for detection of cancer genes compared with single-omics data. In particular, cfRNAs are more sensitive and informative than cfDNAs in terms of detection rate, enriched functional pathways, etc. Moreover, we identify several peripheral immune signatures that are suppressed in patients with cancer. Specifically, we establish a γδ-T cell score and a cancer-associated-fibroblast (CAF) score, providing insights into clinical statuses like cancer stage and survival. Overall, we reveal a cell-free multi-molecular landscape that is useful for blood monitoring in personalized cancer treatment.

Integrated multiplexed assays of variant effect reveal cis-regulatory determinants of catechol-O-methyltransferase gene expression

Multiplexed assays of variant effect are powerful methods to profile the consequences of rare variants on gene expression and organismal fitness. Yet, few studies have integrated several multiplexed assays to map variant effects on gene expression in coding sequences. Here, we pioneered a multiplexed assay based on polysome profiling to measure variant effects on translation at scale, uncovering single-nucleotide variants that increase and decrease ribosome load. By combining high-throughput ribosome load data with multiplexed mRNA and protein abundance readouts, we mapped the cis-regulatory landscape of thousands of catechol-O-methyltransferase (COMT) variants from RNA to protein and found numerous coding variants that alter COMT expression. Finally, we trained machine learning models to map signatures of variant effects on COMT gene expression and uncovered both directional and divergent impacts across expression layers. Our analyses reveal expression phenotypes for thousands of variants in COMT and highlight variant effects on both single and multiple layers of expression. Our findings prompt future studies that integrate several multiplexed assays for the readout of gene expression.

Inflammation primes the kidney for recovery by activating AZIN1 A-to-I editing

The progression of kidney disease varies among individuals, but a general methodology to quantify disease timelines is lacking. Particularly challenging is the task of determining the potential for recovery from acute kidney injury following various insults. Here, we report that quantitation of post-transcriptional adenosine-to-inosine (A-to-I) RNA editing offers a distinct genome-wide signature, enabling the delineation of disease trajectories in the kidney. A well-defined murine model of endotoxemia permitted the identification of the origin and extent of A-to-I editing, along with temporally discrete signatures of double-stranded RNA stress and Adenosine Deaminase isoform switching. We found that A-to-I editing of Antizyme Inhibitor 1 (AZIN1), a positive regulator of polyamine biosynthesis, serves as a particularly useful temporal landmark during endotoxemia. Our data indicate that AZIN1 A-to-I editing, triggered by preceding inflammation, primes the kidney and activates endogenous recovery mechanisms. By comparing genetically modified human cell lines and mice locked in either A-to-I edited or uneditable states, we uncovered that AZIN1 A-to-I editing not only enhances polyamine biosynthesis but also engages glycolysis and nicotinamide biosynthesis to drive the recovery phenotype. Our findings implicate that quantifying AZIN1 A-to-I editing could potentially identify individuals who have transitioned to an endogenous recovery phase. This phase would reflect their past inflammation and indicate their potential for future recovery.

A Comprehensive Evaluation of the Effects of RNA-Editing Proteins ADAR and ADARB1 on the Expression of the Drug-Metabolizing Enzymes in HepaRG Cells

Two RNA-editing proteins, the adenosine deaminase acting on RNA, ADAR, and ADARB1, broadly regulate gene expression in editing-dependent and editing-independent manners. Previous studies showed that the expression of the drug-metabolizing cytochrome P450s (P450s) and UDP glucuronosyltransferases (UGTs) changes upon knockdown (KD) of ADAR or ADARB1 in different hepatic cell lines. To systematically survey the effects of these two ADARs on the expression of P450s and UGTs, we used small interfering RNA in HepaRG cells and tested the association between the expression of the P450s and ADARs in a liver sample cohort (n = 246). KD of ADAR increased the expression of the CYP3As and CYP2C9 and reduced the expression of the others, whereas KD of ADARB1 reduced the expression of nearly all genes tested. ADAR KD also suppressed the induction of most P450s, whereas ADARB1 KD had mixed effects depending on the inducer/gene combination. P450 expression was positively associated with both ADARs in liver samples, consistent with the KD results. However, after adjusting for the expression of transcription factors (TFs) known to regulate P450 expression, the associations disappeared, indicating that the effects of ADAR or ADARB1 primarily occur through TFs. Moreover, we found that the expression of normally spliced CYP3A5 transcripts is increased in both KDs, indicating a direct effect of the ADARs on promoting the usage of the cryptic splice site generated by CYP3A5*3. Taken together, our results revealed the nonoverlapping regulatory effects of ADAR and ADARB1 and supported their broad roles in controlling the expression of drug-metabolizing enzymes in the liver. SIGNIFICANCE STATEMENT: Here, this study systematically surveyed the roles of ADAR and ADARB1 in both basal and induced expression of drug-metabolizing enzymes and assessed their coexpression in liver samples. This study's results support that ADAR and ADARB1 regulate the expression of the drug-metabolizing enzymes in the liver, suggesting that factors affecting ADAR expression also have the potential to impact drug metabolism.

Autorecoding A-to-I RNA editing sites in the Adar gene underwent compensatory gains and losses in major insect clades

As one of the most prevalent RNA modifications in animals, adenosine-to-inosine (A-to-I) RNA editing facilitates the environmental adaptation of organisms by diversifying the proteome in a temporal-spatial manner. In flies and bees, the editing enzyme Adar has independently gained two different autorecoding sites that form an autofeedback loop, stabilizing the overall editing efficiency. This ensures cellular homeostasis by keeping the normal function of target genes. However, in a broader range of insects, the evolutionary dynamics and significance of this Adar autoregulatory mechanism are unclear. We retrieved the genomes of 377 arthropod species covering the five major insect orders (Hemiptera, Hymenoptera, Coleoptera, Diptera, and Lepidoptera) and aligned the Adar autorecoding sites across all genomes. We found that the two autorecoding sites underwent compensatory gains and losses during the evolution of two orders with the most sequenced species (Diptera and Hymenoptera), and that the two editing sites were mutually exclusive among them: One editable site is significantly linked to another uneditable site. This autorecoding mechanism of Adar could flexibly diversify the proteome and stabilize global editing activity. Many insects independently selected different autorecoding sites to achieve a feedback loop and regulate the global RNA editome, revealing an interesting phenomenon during evolution. Our study reveals the evolutionary force acting on accurate regulation of RNA editing activity in insects and thus deepens our understanding of the functional importance of RNA editing in environmental adaptation and evolution.

Transcriptome-wide profiling of A-to-I RNA editing by Slic-seq

Adenosine-to-inosine (A-to-I) RNA editing is a post-transcriptional processing event involved in diversifying the transcriptome and is responsible for various biological processes. In this context, we developed a new method based on the highly selective cleavage activity of Endonuclease V against Inosine and the universal activity of sodium periodate against all RNAs to enrich the inosine-containing RNA and accurately identify the editing sites. We validated the reliability of our method in human brain in both Alu and non-Alu elements. The conserved sites of A-to-I editing in human cells (HEK293T, HeLa, HepG2, K562 and MCF-7) primarily occurs in the 3'UTR of the RNA, which are highly correlated with RNA binding and protein binding. Analysis of the editing sites between the human brain and mouse brain revealed that the editing of exons is more conserved than that in other regions. This method was applied to three neurological diseases (Alzheimer's, epilepsy and ageing) of mouse brain, reflecting that A-to-I editing sites significantly decreased in neuronal activity genes.

RNA in situ conformation sequencing reveals novel long-range RNA structures with impact on splicing

Over recent years, long-range RNA structure has emerged as a factor that is fundamental to alternative splicing regulation. An increasing number of human disorders are now being associated with splicing defects; hence it is essential to develop methods that assess long-range RNA structure experimentally. RNA in situ conformation sequencing (RIC-seq) is a method that recapitulates RNA structure within physiological RNA-protein complexes. In this work, we juxtapose pairs of conserved complementary regions (PCCRs) that were predicted in silico with the results of RIC-seq experiments conducted in seven human cell lines. We show statistically that RIC-seq support of PCCRs correlates with their properties, such as equilibrium free energy, presence of compensatory substitutions, and occurrence of A-to-I RNA editing sites and forked eCLIP peaks. Exons enclosed in PCCRs that are supported by RIC-seq tend to have weaker splice sites and lower inclusion rates, which is indicative of post-transcriptional splicing regulation mediated by RNA structure. Based on these findings, we prioritize PCCRs according to their RIC-seq support and show, using antisense nucleotides and minigene mutagenesis, that PCCRs in two disease-associated human genes, PHF20L1 and CASK, and also PCCRs in their murine orthologs, impact alternative splicing. In sum, we demonstrate how RIC-seq experiments can be used to discover functional long-range RNA structures, and particularly those that regulate alternative splicing.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk post-mortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA-editing and alternative polyadenylation (APA), within a cell-type-specific population of human neural progenitors and neurons. More RNA-editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting genetically mediated post-transcriptional regulation during brain development lead to differences in brain function.

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

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