Elucidating the inosinome: global approaches to adenosine-to-inosine RNA editing

RNA modifications identification based on chemical reactions

RNA modification identification is an emerging field in epigenetics due to its indispensable regulatory role in the cell life cycle. With advancements in identification methods, an increasing number of RNA modifications has been discovered, thereby driving the development of more efficient and accurate techniques for localizing modified RNAs and elucidating their functions. High-throughput sequencing approaches for modified RNA detection can be categorized into antibody-based, enzymatic-based, and chemical-labeling-based methods. Given the intrinsic chemical reactions involved in all biochemical processes, we provide a comprehensive review of recent advancements in artificial chemical labeling and transformations of ten distinct RNA modifications and their applications in sequencing. Our aim is to contribute to a deeper understanding of the mechanisms underlying these modifications. We focus on the chemical reactions associated with RNA modifications and briefly compare the advantages and disadvantages of detection methods based on these reactions. Additionally, we introduce several approaches that identify multiple modifications through chemical labeling. As the field of RNA modification research continues to expand, we anticipate that the techniques and insights presented in this review will serve as a valuable resource for future studies aimed at further elucidating the functional roles of RNA modifications in biological processes.

Temperature-Dependent tRNA Modifications in Bacillales

Transfer RNA (tRNA) modifications are essential for the temperature adaptation of thermophilic and psychrophilic organisms as they control the rigidity and flexibility of transcripts. To further understand how specific tRNA modifications are adjusted to maintain functionality in response to temperature fluctuations, we investigated whether tRNA modifications represent an adaptation of bacteria to different growth temperatures (minimal, optimal, and maximal), focusing on closely related psychrophilic (P. halocryophilus and E. sibiricum), mesophilic (B. subtilis), and thermophilic (G. stearothermophilus) Bacillales. Utilizing an RNA sequencing approach combined with chemical pre-treatment of tRNA samples, we systematically profiled dihydrouridine (D), 4-thiouridine (s4U), 7-methyl-guanosine (m7G), and pseudouridine (Ψ) modifications at single-nucleotide resolution. Despite their close relationship, each bacterium exhibited a unique tRNA modification profile. Our findings revealed increased tRNA modifications in the thermophilic bacterium at its optimal growth temperature, particularly showing elevated levels of s4U8 and Ψ55 modifications compared to non-thermophilic bacteria, indicating a temperature-dependent regulation that may contribute to thermotolerance. Furthermore, we observed higher levels of D modifications in psychrophilic and mesophilic bacteria, indicating an adaptive strategy for cold environments by enhancing local flexibility in tRNAs. Our method demonstrated high effectiveness in identifying tRNA modifications compared to an established tool, highlighting its potential for precise tRNA profiling studies.

[Advances in mapping analysis of ribonucleic acid modifications through sequencing]

Over 170 chemical modifications have been discovered in various types of ribonucleic acids (RNAs), including messenger RNA (mRNA), ribosomal RNA (rRNA), transfer RNA (tRNA), and small nuclear RNA (snRNA). These RNA modifications play crucial roles in a wide range of biological processes such as gene expression regulation, RNA stability maintenance, and protein translation. RNA modifications represent a new dimension of gene expression regulation known as the "epitranscriptome". The discovery of RNA modifications and the relevant writers, erasers, and readers provides an important basis for studies on the dynamic regulation and physiological functions of RNA modifications. Owing to the development of detection technologies for RNA modifications, studies on RNA epitranscriptomes have progressed to the single-base resolution, multilayer, and full-coverage stage. Transcriptome-wide methods help discover new RNA modification sites and are of great importance for elucidating the molecular regulatory mechanisms of epitranscriptomics, exploring the disease associations of RNA modifications, and understanding their clinical applications. The existing RNA modification sequencing technologies can be categorized according to the pretreatment approach and sequencing principle as direct high-throughput sequencing, antibody-enrichment sequencing, enzyme-assisted sequencing, chemical labeling-assisted sequencing, metabolic labeling sequencing, and nanopore sequencing technologies. These methods, as well as studies on the functions of RNA modifications, have greatly expanded our understanding of epitranscriptomics. In this review, we summarize the recent progress in RNA modification detection technologies, focusing on the basic principles, advantages, and limitations of different methods. Direct high-throughput sequencing methods do not require complex RNA pretreatment and allow for the mapping of RNA modifications using conventional RNA sequencing methods. However, only a few RNA modifications can be analyzed by high-throughput sequencing. Antibody enrichment followed by high-throughput sequencing has emerged as a crucial approach for mapping RNA modifications, significantly advancing the understanding of RNA modifications and their regulatory functions in different species. However, the resolution of antibody-enrichment sequencing is limited to approximately 100-200 bp. Although chemical crosslinking techniques can achieve single-base resolution, these methods are often complex, and the specificity of the antibodies used in these methods has raised concerns. In particular, the issue of off-target binding by the antibodies requires urgent attention. Enzyme-assisted sequencing has improved the accuracy of the localization analysis of RNA modifications and enables stoichiometric detection with single-base resolution. However, the enzymes used in this technique show poor reactivity, specificity, and sequence preference. Chemical labeling sequencing has become a widely used approach for profiling RNA modifications, particularly by altering reverse transcription (RT) signatures such as RT stops, misincorporations, and deletions. Chemical-assisted sequencing provides a sequence-independent RNA modification detection strategy that enables the localization of multiple RNA modifications. Additionally, when combined with the biotin-streptavidin affinity method, low-abundance RNA modifications can be enriched and detected. Nevertheless, the specificity of many chemical reactions remains problematic, and the development of specific reaction probes for particular modifications should continue in the future to achieve the precise localization of RNA modifications. As an indirect localization method, metabolic labeling sequencing specifically localizes the sites at which modifying enzymes act, which is of great significance in the study of RNA modification functions. However, this method is limited by the intracellular labeling of RNA and cannot be applied to biological samples such as clinical tissues and blood samples. Nanopore sequencing is a direct RNA-sequencing method that does not require RT or the polymerase chain reaction (PCR). However, challenges in analyzing the data obtained from nanopore sequencing, such as the high rate of false positives, must be resolved. Discussing sequencing analysis methods for various types of RNA modifications is instructive for the future development of novel RNA modification mapping technologies, and will aid studies on the functions of RNA modifications across the entire transcriptome.

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.

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

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

Detection technologies for RNA modifications

To date, more than 170 chemical modifications have been characterized in RNA, providing a new layer of gene expression regulation termed the 'epitranscriptome'. RNA modification detection methods and tools advance the functional studies of the epitranscriptome. According to the detection throughput and principles, existing RNA modification detection technologies can be categorized into four classes, including quantification methods, locus-specific detection methods, next-generation sequencing-based detection technologies and nanopore direct RNA sequencing-based technologies. In this review, we summarize the current knowledge about these RNA modification detection technologies and discuss the challenges for the existing detection tools, providing information for a comprehensive understanding of the epitranscriptome.

RNA modifications: importance in immune cell biology and related diseases

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

Crucial Roles of miR-625 in Human Cancer

Genetic and epigenetic characteristics are core factors of cancer. MicroRNAs (miRNAs) are small non-coding RNAs which regulate gene expression at the post-transcriptional level via binding to corresponding mRNAs. Recently, increasing evidence has proven that miRNAs regulate the occurrence and development of human cancer. Here, we mainly review the abnormal expression of miR-625 in a variety of cancers. In summarizing the role and potential molecular mechanisms of miR-625 in various tumors in detail, we reveal that miR-625 is involved in a variety of biological processes, such as cell proliferation, invasion, migration, apoptosis, cell cycle regulation, and drug resistance. In addition, we discuss the lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA networks and briefly explain the specific mechanisms of competing endogenous RNAs. In conclusion, we reveal the potential value of miR-625 in cancer diagnosis, treatment, and prognosis and hope to provide new ideas for the clinical application of miR-625.

RNA modifications as emerging therapeutic targets

The field of epitranscriptome, posttranscriptional modifications to RNAs, is still growing up and has presented substantial evidences for the role of RNA modifications in diseases. In terms of new drug development, RNA modifications have a great promise for therapy. For example, more than 170 type of modifications exist in various types of RNAs. Regulatory genes and their roles in critical biological process have been identified and they are associated with several diseases. Current data, for example, identification of inhibitors targeting RNA modifications regulatory genes, strongly support the idea that RNA modifications have potential as emerging therapeutic targets. Therefore, in this review, RNA modifications and regulatory genes were comprehensively documented in terms of drug development by summarizing the findings from previous studies. It was discussed how RNA modifications or regulatory genes can be targeted by altering molecular mechanisms. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > RNA Editing and Modification.

The transcribed ultraconserved region uc.160+ enhances processing and A-to-I editing of the miR-376 cluster: hypermethylation improves glioma prognosis

Transcribed ultraconserved regions (T‐UCRs) are noncoding RNAs derived from DNA sequences that are entirely conserved across species. Their expression is altered in many tumor types, and, although a role for T‐UCRs as regulators of gene expression has been proposed, their functions remain largely unknown. Herein, we describe the epigenetic silencing of the uc.160+ T‐UCR in gliomas and mechanistically define a novel RNA–RNA regulatory network in which uc.160+ modulates the biogenesis of several members of the miR‐376 cluster. This includes the positive regulation of primary microRNA (pri‐miRNA) cleavage and an enhanced A‐to‐I editing on its mature sequence. As a consequence, the expression of uc.160+ affects the downstream, miR‐376‐regulated genes, including the transcriptional coregulators RING1 and YY1‐binding protein (RYBP) and forkhead box P2 (FOXP2). Finally, we elucidate the clinical impact of our findings, showing that hypermethylation of the uc.160+ CpG island is an independent prognostic factor associated with better overall survival in lower‐grade gliomas, highlighting the importance of T‐UCRs in cancer pathophysiology.

RNA A-to-I editing, environmental exposure, and human diseases

Epigenetic modifications have gained attention since they can be potentially changed with environmental stimuli and can be associated with adverse health outcomes. Epitranscriptome field has begun to attract attention with several aspects since RNA modifications have been linked with critical biological processes and implicated in diseases. Several RNA modifications have been identified as reversible indicating the dynamic features of modification which can be altered by environmental cues. Currently, we know more than 150 RNA modifications in different organisms and on different bases which are modified by various chemical groups. RNA editing, which is one of the RNA modifications, occurs after transcription, which results in RNA sequence different from its corresponding DNA sequence. Emerging evidence reveals the functions of RNA editing as well as the association between RNA editing and diseases. However, the RNA editing field is beginning to grow up and needs more empirical evidence in regard to disease and toxicology. Thus, this review aims to provide the current evidence-based studies on RNA editing modifying genes for genotoxicity and cancer. The review presented the association between environmental xenobiotics exposure and RNA editing modifying genes and focused on the association between the expression of RNA editing modifying genes and cancer. Furthermore, we discussed the future directions of scientific studies in the area of RNA modifications, especially in the RNA editing field, and provided a knowledge-based framework for further studies.

RNA modifications in cardiovascular diseases, the potential therapeutic targets

More than one hundred RNA modifications decorate the chemical and topological properties of these ribose nucleotides, thereby executing their biological functions through post-transcriptional regulation. In cardiovascular diseases, a wide range of RNA modifications including m⁶A (N6-adenosine methylation), m⁵C (5-methylcytidin), Nm (2'-O-ribose-methylation), Ψ (pseudouridine), m⁷G (N7-methylguanosine), and m¹A (N1-adenosine methylation) have been found in tRNA, rRNA, mRNA and other noncoding RNA, which can function as a novel mechanism in metabolic syndrome, heart failure, coronary heart disease, and hypertension. In this review, we will summarize the current understanding of the regulatory roles and significance of several types of RNA modifications in CVDs (cardiovascular diseases) and the interplay between RNA modifications and noncoding RNA, epigenetics. Finally, we will focus on the potential therapeutic strategies by using RNA modifications.

A Bayesian Framework for Inferring the Influence of Sequence Context on Point Mutations

The probability of point mutations is expected to be highly influenced by the flanking nucleotides that surround them, known as the sequence context. This phenomenon may be mainly attributed to the enzyme that modifies or mutates the genetic material, because most enzymes tend to have specific sequence contexts that dictate their activity. Here, we develop a statistical model that allows for the detection and evaluation of the effects of different sequence contexts on mutation rates from deep population sequencing data. This task is computationally challenging, as the complexity of the model increases exponentially as the context size increases. We established our novel Bayesian method based on sparse model selection methods, with the leading assumption that the number of actual sequence contexts that directly influence mutation rates is minuscule compared with the number of possible sequence contexts. We show that our method is highly accurate on simulated data using pentanucleotide contexts, even when accounting for noisy data. We next analyze empirical population sequencing data from polioviruses and HIV-1 and detect a significant enrichment in sequence contexts associated with deamination by the cellular deaminases ADAR 1/2 and APOBEC3G, respectively. In the current era, where next-generation sequencing data are highly abundant, our approach can be used on any population sequencing data to reveal context-dependent base alterations and may assist in the discovery of novel mutable sites or editing sites.

Global RNA editome landscape discovers reduced RNA editing in glioma: loss of editing of gamma-amino butyric acid receptor alpha subunit 3 (GABRA3) favors glioma migration and invasion

Background

Gliomas are the most common and lethal type of intracranial tumors. With the current treatment regime, the median survival of patients with grade IV glioma (glioblastoma/GBM) remains at 14-16 months. RNA editing modifies the function and regulation of transcripts. The development of glial tumors may be caused by altered RNA editing events.

Methods

In this study, we uncover the global RNA editome landscape of glioma patients from RNA-seq data of control, lower grade glioma (LGG) and GBM samples (n = 1,083).

Results

A-to-I editing events were found to comprise 80% of the total editing events of which 96% were located in the Alu regions. The total RNA editing events were found to be reduced in glioma compared to control samples. More specifically, we found Gamma-aminobutyric acid type A receptor alpha3 (GABRA3) to be edited (c.1026 A-to-G; pI343M) in 73% (editing ratio 0.8) of control samples compared to LGG (28.96%; 0.47) and GBM (5.2%; 0.53) samples. GABRA3 transcript level was found to be downregulated in glioma compared to control in a grade-specific manner with GBMs having the lowest level of the transcript. Further, GABRA3 transcripts were observed to be higher in edited compared to unedited glioma samples. The transcript and protein levels of exogenously expressed gene were found to be higher for edited compared to unedited GABRA3 in glioma cells. Further, exogenously expressed edited GABRA3 inhibited migration and invasion of glioma cells efficiently but not the unedited GABRA3.

Conclusion

Collectively, our study discovered a reduction in RNA editing during glioma development. We further demonstrate that elevated RNA editing maintains a high level of GABRA3 RNA and protein in normal glial cells which provides a less migratory environment for the normal functioning of the brain. In contrast, the reduction in GABRA3 protein levels, due to lower stability of unedited RNA, results in the loss of function which confers an aggressive phenotype to GBM tumor.

Epitranscriptomics: RNA Modifications in Bacteria and Archaea

The increasingly complex functionality of RNA is contrasted by its simple chemical composition. RNA is generally built from only four different nucleotides (adenine, guanine, cytosine, and uracil). To date, >160 chemical modifications are known to decorate RNA molecules and thereby alter their function or stability. Many RNA modifications are conserved throughout bacteria, archaea, and eukaryotes, while some are unique to each branch of life. Most known modifications occur at internal positions, while there is limited diversity at the termini. The dynamic nature of RNA modifications and newly discovered regulatory functions of some of these RNA modifications gave birth to a new field, now often referred to as "epitranscriptomics." This review highlights the major developments in this field and summarizes detection principles for internal as well as 5'-terminal mRNA modifications in prokaryotes and archaea to investigate their biological significance.

The nucleotides they are a-changin': function of RNA binding proteins in post-transcriptional messenger RNA editing and modification in Arabidopsis

During and after transcription, the fate of an RNA molecule is almost entirely directed by the cohorts of interacting RNA-binding proteins (RBPs). RBPs regulate all stages of the life cycle of a messenger RNA (mRNA) molecule, including splicing, polyadenylation, transport out of the nucleus, RNA stability, and translation. In addition to these functions, RBPs can function to modify or edit the sequences encoded by the RNA. While the sequence for each transcript is determined in the genome, by the time an RNA reaches its final fate, the sequence may have been edited, where one nucleotide is converted to another, or modified, where a chemical group, or sometimes others moieties, are covalently linked to a nucleotide base. These changes to the RNA sequence have major consequences on the function of the RNA. Additionally, variation in the levels of the RBPs that perform the editing or modification can drastically affect the fitness of an organism. Here, we review RBPs that are known to edit or modify RNA ribonucleotides, focusing on the RNA editing ability of the pentatricopeptide repeat (PPR) proteins and the RBPs that modify adenosine to N⁶- methyladenosine.

Adenosine-to-Inosine Editing of MicroRNA-487b Alters Target Gene Selection After Ischemia and Promotes Neovascularization

Rationale:

Adenosine-to-inosine editing of microRNAs has the potential to cause a shift in target site selection. 2′-O-ribose-methylation of adenosine residues, however, has been shown to inhibit adenosine-to-inosine editing.

Objective:

To investigate whether angiomiR miR487b is subject to adenosine-to-inosine editing or 2′-O-ribose-methylation during neovascularization.

Methods and Results:

Complementary DNA was prepared from C57BL/6-mice subjected to hindlimb ischemia. Using Sanger sequencing and endonuclease digestion, we identified and validated adenosine-to-inosine editing of the miR487b seed sequence. In the gastrocnemius muscle, pri-miR487b editing increased from 6.7±0.4% before to 11.7±1.6% ( P =0.02) 1 day after ischemia. Edited pri-miR487b is processed into a novel microRNA, edited miR487b, which is also upregulated after ischemia. We confirmed editing of miR487b in multiple human primary vascular cell types. Short interfering RNA–mediated knockdown demonstrated that editing is adenosine deaminase acting on RNA 1 and 2 dependent. Using reverse-transcription at low dNTP concentrations followed by quantitative-PCR, we found that the same adenosine residue is methylated in mice and human primary cells. In the murine gastrocnemius, the estimated methylation fraction increased from 32.8±14% before to 53.6±12% 1 day after ischemia. Short interfering RNA knockdown confirmed that methylation is fibrillarin dependent. Although we could not confirm that methylation directly inhibits editing, we do show that adenosine deaminase acting on RNA 1 and 2 and fibrillarin negatively influence each other’s expression. Using multiple luciferase reporter gene assays, we could demonstrate that editing results in a complete switch of target site selection. In human primary cells, we confirmed the shift in miR487b targeting after editing, resulting in a edited miR487b targetome that is enriched for multiple proangiogenic pathways. Furthermore, overexpression of edited miR487b, but not wild-type miR487b, stimulates angiogenesis in both in vitro and ex vivo assays.

Conclusions:

MiR487b is edited in the seed sequence in mice and humans, resulting in a novel, proangiogenic microRNA with a unique targetome. The rate of miR487b editing, as well as 2′-O-ribose-methylation, is increased in murine muscle tissue during postischemic neovascularization. Our findings suggest miR487b editing plays an intricate role in postischemic neovascularization.

Integrated sequencing of exome and mRNA of large-sized single cells

Current approaches of single cell DNA-RNA integrated sequencing are difficult to call SNPs, because a large amount of DNA and RNA is lost during DNA-RNA separation. Here, we performed simultaneous single-cell exome and transcriptome sequencing on individual mouse oocytes. Using microinjection, we kept the nuclei intact to avoid DNA loss, while retaining the cytoplasm inside the cell membrane, to maximize the amount of DNA and RNA captured from the single cell. We then conducted exome-sequencing on the isolated nuclei and mRNA-sequencing on the enucleated cytoplasm. For single oocytes, exome-seq can cover up to 92% of exome region with an average sequencing depth of 10+, while mRNA-sequencing reveals more than 10,000 expressed genes in enucleated cytoplasm, with similar performance for intact oocytes. This approach provides unprecedented opportunities to study DNA-RNA regulation, such as RNA editing at single nucleotide level in oocytes. In future, this method can also be applied to other large cells, including neurons, large dendritic cells and large tumour cells for integrated exome and transcriptome sequencing.

CAPS1 RNA Editing Promotes Dense Core Vesicle Exocytosis

Calcium-dependent activator protein for secretion 1 (CAPS1) plays a distinct role in the priming step of dense core vesicle (DCV) exocytosis. CAPS1 pre-mRNA is known to undergo adenosine-to-inosine RNA editing in its coding region, which results in a glutamate-to-glycine conversion at a site in its C-terminal region. However, the physiological significance of CAPS1 RNA editing remains elusive. Here, we created mutant mice in which edited CAPS1 was solely expressed. These mice were lean due to increased energy expenditure caused by physical hyperactivity. Electrophysiological and biochemical analyses demonstrated that the exocytosis of DCVs was upregulated in the chromaffin cells and neurons of these mice. Furthermore, we showed that edited CAPS1 bound preferentially to the activated form of syntaxin-1A, a component of the exocytotic fusion complex. These findings suggest that RNA editing regulates DCV exocytosis in vivo, affecting physical activity.

Epitranscriptome sequencing technologies: decoding RNA modifications

In recent years, major breakthroughs in RNA-modification-mediated regulation of gene expression have been made, leading to the emerging field of epitranscriptomics.Our understanding of the distribution, regulation and function of these dynamic RNA modifications is based on sequencing technologies. In this Review, we focus on the major mRNA modifications in the transcriptome of eukaryotic cells: N6-methyladenosine, N6, 2'-O-dimethyladenosine, 5-methylcytidine, 5-hydroxylmethylcytidine, inosine, pseudouridine and N¹-methyladenosine. We discuss the sequencing technologies used to profile these epitranscriptomic marks, including scale, resolution, quantitative feature, pre-enrichment capability and the corresponding bioinformatics tools. We also discuss the challenges of epitranscriptome profiling and highlight the prospect of future detection tools. We aim to guide the choice of different detection methods and inspire new ideas in RNA biology.

RNA editing signature during myeloid leukemia cell differentiation

Adenosine deaminases acting on RNA (ADARs) are key proteins for hematopoietic stem cell self-renewal and for survival of differentiating progenitor cells. However, their specific role in myeloid cell maturation has been poorly investigated. Here we show that ADAR1 is present at basal level in the primary myeloid leukemia cells obtained from patients at diagnosis as well as in myeloid U-937 and THP1 cell lines and its expression correlates with the editing levels. Upon phorbol-myristate acetate or Vitamin D3/granulocyte macrophage colony-stimulating factor (GM-CSF)-driven differentiation, both ADAR1 and ADAR2 enzymes are upregulated, with a concomitant global increase of A-to-I RNA editing. ADAR1 silencing caused an editing decrease at specific ADAR1 target genes, without, however, interfering with cell differentiation or with ADAR2 activity. Remarkably, ADAR2 is absent in the undifferentiated cell stage, due to its elimination through the ubiquitin-proteasome pathway, being strongly upregulated at the end of the differentiation process. Of note, peripheral blood monocytes display editing events at the selected targets similar to those found in differentiated cell lines. Taken together, the data indicate that ADAR enzymes play important and distinct roles in myeloid cells.

Chemical Modifications to RNA: A New Layer of Gene Expression Regulation

The first chemical modification to RNA was discovered nearly 60 years ago; to date, more than 100 chemically distinct modifications have been identified in cellular RNA. With the recent development of novel chemical and/or biochemical methods, dynamic modifications to RNA have been identified in the transcriptome, including N⁶-methyladenosine (m⁶A), inosine (I), 5-methylcytosine (m⁵C), pseudouridine (Ψ), 5-hydroxymethylcytosine (hm⁵C), and N¹-methyladenosine (m¹A). Collectively, the multitude of RNA modifications are termed epitranscriptome, leading to the emerging field of epitranscriptomics. In this review, we primarily focus on recently reported chemical modifications to mRNA; we discuss their chemical properties, biological functions, and mechanisms with an emphasis on their high-throughput detection methods. We also envision that future tools, particularly novel chemical biology methods, could further facilitate and enable studies in the field of epitranscriptomics.

Elucidating the cancer-specific genetic alteration spectrum of glioblastoma derived cell lines from whole exome and RNA sequencing

Cell lines derived from tumor tissues have been used as a valuable system to study gene regulation and cancer development. Comprehensive characterization of the genetic background of cell lines could provide clues on novel genes responsible for carcinogenesis and help in choosing cell lines for particular studies. Here, we have carried out whole exome and RNA sequencing of commonly used glioblastoma (GBM) cell lines (U87, T98G, LN229, U343, U373 and LN18) to unearth single nucleotide variations (SNVs), indels, differential gene expression, gene fusions and RNA editing events. We obtained an average of 41,071 SNVs out of which 1,594 (3.88%) were potentially cancer-specific. The cell lines showed frequent SNVs and indels in some of the genes that are known to be altered in GBM- EGFR, TP53, PTEN, SPTA1 and NF1. Chromatin modifying genes- ATRX, MLL3, MLL4, SETD2 and SRCAP also showed alterations. While no cell line carried IDH1 mutations, five cell lines showed hTERT promoter activating mutations with a concomitant increase in hTERT transcript levels. Five significant gene fusions were found of which NUP93-CYB5B was validated. An average of 18,949 RNA editing events was also obtained. Thus we have generated a comprehensive catalogue of genetic alterations for six GBM cell lines.

Adenosine-to-inosine RNA editing controls cathepsin S expression in atherosclerosis by enabling HuR-mediated post-transcriptional regulation

Adenosine-to-inosine (A-to-I) RNA editing, which is catalyzed by a family of adenosine deaminase acting on RNA (ADAR) enzymes, is important in the epitranscriptomic regulation of RNA metabolism. However, the role of A-to-I RNA editing in vascular disease is unknown. Here we show that cathepsin S mRNA (CTSS), which encodes a cysteine protease associated with angiogenesis and atherosclerosis, is highly edited in human endothelial cells. The 3' untranslated region (3' UTR) of the CTSS transcript contains two inverted repeats, the AluJo and AluSx⁺ regions, which form a long stem-loop structure that is recognized by ADAR1 as a substrate for editing. RNA editing enables the recruitment of the stabilizing RNA-binding protein human antigen R (HuR; encoded by ELAVL1) to the 3' UTR of the CTSS transcript, thereby controlling CTSS mRNA stability and expression. In endothelial cells, ADAR1 overexpression or treatment of cells with hypoxia or with the inflammatory cytokines interferon-γ and tumor-necrosis-factor-α induces CTSS RNA editing and consequently increases cathepsin S expression. ADAR1 levels and the extent of CTSS RNA editing are associated with changes in cathepsin S levels in patients with atherosclerotic vascular diseases, including subclinical atherosclerosis, coronary artery disease, aortic aneurysms and advanced carotid atherosclerotic disease. These results reveal a previously unrecognized role of RNA editing in gene expression in human atherosclerotic vascular diseases.

Reconciling proteomics with next generation sequencing

Both genomics and proteomics technologies have matured in the last decade to a level where they are able to deliver system-wide data on the qualitative and quantitative abundance of their respective molecular entities, that is DNA/RNA and proteins. A next logical step is the collective use of these technologies, ideally gathering data on matching samples. The first large scale so-called proteogenomics studies are emerging, and display the benefits each of these layers of analysis has on the other layers to together generate more meaningful insight into the connection between the phenotype/physiology and genotype of the system under study. Here we review a selected number of these studies, highlighting what they can uniquely deliver. We also discuss the future potential and remaining challenges, from a somewhat proteome biased perspective.

Chemical Modifications Mark Alternatively Spliced and Uncapped Messenger RNAs in Arabidopsis

Posttranscriptional chemical modification of RNA bases is a widespread and physiologically relevant regulator of RNA maturation, stability, and function. While modifications are best characterized in short, noncoding RNAs such as tRNAs, growing evidence indicates that mRNAs and long noncoding RNAs (lncRNAs) are likewise modified. Here, we apply our high-throughput annotation of modified ribonucleotides (HAMR) pipeline to identify and classify modifications that affect Watson-Crick base pairing at three different levels of the Arabidopsis thaliana transcriptome (polyadenylated, small, and degrading RNAs). We find this type of modifications primarily within uncapped, degrading mRNAs and lncRNAs, suggesting they are the cause or consequence of RNA turnover. Additionally, modifications within stable mRNAs tend to occur in alternatively spliced introns, suggesting they regulate splicing. Furthermore, these modifications target mRNAs with coherent functions, including stress responses. Thus, our comprehensive analysis across multiple RNA classes yields insights into the functions of covalent RNA modifications in plant transcriptomes.

Transcriptome-wide measurement of plant RNA secondary structure

RNAs fold into intricate and precise secondary structures. These structural patterns regulate multiple steps of the RNA lifecycle, while also conferring catalytic and scaffolding functions to certain transcripts. Therefore, a full understanding of RNA posttranscriptional regulation requires a comprehensive picture of secondary structure. Here, we review several high throughput sequencing-based methods to globally survey plant RNA secondary structure. These methods are more accurate than computational prediction, and more scalable than physical techniques such as crystallography. We note hurdles to reliably measuring secondary structure, including RNA-binding proteins, RNA base modifications, and intramolecular duplexes. Finally, we survey the functional knowledge that has been gleaned from each of these methods, and identify some unanswered questions that remain.

Measles Virus Defective Interfering RNAs Are Generated Frequently and Early in the Absence of C Protein and Can Be Destabilized by Adenosine Deaminase Acting on RNA-1-Like Hypermutations

Defective interfering RNAs (DI-RNAs) of the viral genome can form during infections of negative-strand RNA viruses and outgrow full-length viral genomes, thereby modulating the severity and duration of infection. Here we document the frequent de novo generation of copy-back DI-RNAs from independent rescue events both for a vaccine measles virus (vac2) and for a wild-type measles virus (IC323) as early as passage 1 after virus rescue. Moreover, vaccine and wild-type C-protein-deficient (C-protein-knockout [CKO]) measles viruses generated about 10 times more DI-RNAs than parental virus, suggesting that C enhances the processivity of the viral polymerase. We obtained the nucleotide sequences of 65 individual DI-RNAs, identified breakpoints and reinitiation sites, and predicted their structural features. Several DI-RNAs possessed clusters of A-to-G or U-to-C transitions. Sequences flanking these mutation sites were characteristic of those favored by adenosine deaminase acting on RNA-1 (ADAR1), which catalyzes in double-stranded RNA the C-6 deamination of adenosine to produce inosine, which is recognized as guanosine, a process known as A-to-I RNA editing. In individual DI-RNAs the transitions were of the same type and occurred on both sides of the breakpoint. These patterns of mutations suggest that ADAR1 edits unencapsidated DI-RNAs that form double-strand RNA structures. Encapsidated DI-RNAs were incorporated into virus particles, which reduced the infectivity of virus stocks. The CKO phenotype was dominant: DI-RNAs derived from vac2 with a CKO suppressed the replication of vac2, as shown by coinfections of interferon-incompetent lymphatic cells with viruses expressing different fluorescent reporter proteins. In contrast, coinfection with a C-protein-expressing virus did not counteract the suppressive phenotype of DI-RNAs.

IMPORTANCE

Recombinant measles viruses (MVs) are in clinical trials as cancer therapeutics and as vectored vaccines for HIV-AIDS and other infectious diseases. The efficacy of MV-based vectors depends on their replication proficiency and immune activation capacity. Here we document that copy-back defective interfering RNAs (DI-RNAs) are generated by recombinant vaccine and wild-type MVs immediately after rescue. The MV C protein interferes with DI-RNA generation and may enhance the processivity of the viral polymerase. We frequently detected clusters of A-to-G or U-to-C transitions and noted that sequences flanking individual mutations contain motifs favoring recognition by the adenosine deaminase acting on RNA-1 (ADAR1). The consistent type of transitions on the DI-RNAs indicates that these are direct substrates for editing by ADAR1. The ADAR1-mediated biased hypermutation events are consistent with the protein kinase R (PKR)-ADAR1 balancing model of innate immunity activation. We show by coinfection that the C-defective phenotype is dominant.

To translate, or not to translate: viral and host mRNA regulation by interferon-stimulated genes

•

Translational regulation of viral and host mRNA is an important function of a subset of interferon-stimulated genes (ISGs).

•

These ISGs repress translation through binding to viral RNA, direct interaction with, or perturbation of the translation machinery components.

•

Several ISGs localize to cytoplasmic granules such as stress granules (SGs) and processing bodies (PBs), and interfere with the processing or function of microRNAs (miRNAs).

Type I interferon (IFN) is one of the first lines of cellular defense against viral pathogens. As a result of IFN signaling, a wide array of IFN-stimulated gene (ISG) products is upregulated to target different stages of the viral life cycle. We review recent findings implicating a subset of ISGs in translational regulation of viral and host mRNAs. Translation inhibition is mediated either by binding to viral RNA or by disrupting physiological interactions or levels of the translation complex components. In addition, many of these ISGs localize to translationally silent cytoplasmic granules, such as stress granules and processing bodies, and intersect with the microRNA (miRNA)-mediated silencing pathway to regulate translation of cellular mRNAs.

Evidence for multiple, distinct ADAR-containing complexes in Xenopus laevis

ADAR (adenosine deaminase acting on RNA) is an RNA-editing enzyme present in most metazoans that converts adenosines in double-stranded RNA targets into inosines. Although the RNA targets of ADAR-mediated editing have been extensively cataloged, our understanding of the cellular function of such editing remains incomplete. We report that long, double-stranded RNA added to Xenopus laevis egg extract is incorporated into an ADAR-containing complex whose protein components resemble those of stress granules. This complex localizes to microtubules, as assayed by accumulation on meiotic spindles. We observe that the length of a double-stranded RNA influences its incorporation into the microtubule-localized complex. ADAR forms a similar complex with endogenous RNA, but the endogenous complex fails to localize to microtubules. In addition, we characterize the endogenous, ADAR-associated RNAs and discover that they are enriched for transcripts encoding transcriptional regulators, zinc-finger proteins, and components of the secretory pathway. Interestingly, association with ADAR correlates with previously reported translational repression in early embryonic development. This work demonstrates that ADAR is a component of two, distinct ribonucleoprotein complexes that contain different types of RNAs and exhibit diverse cellular localization patterns. Our findings offer new insight into the potential cellular functions of ADAR.

De novo identification and quantification of single amino-acid variants in human brain

The detection of single amino-acid variants (SAVs) usually depends on single-nucleotide polymorphisms (SNPs) database. Here, we describe a novel method that discovers SAVs at proteome level independent of SNPs data. Using mass spectrometry-based de novo sequencing algorithm, peptide-candidates are identified and compared with theoretical protein database to generate SAVs under pairing strategy, which is followed by database re-searching to control false discovery rate. In human brain tissues, we can confidently identify known and novel protein variants with diverse origins. Combined with DNA/RNA sequencing, we verify SAVs derived from DNA mutations, RNA alternative splicing, and unknown post-transcriptional mechanisms. Furthermore, quantitative analysis in human brain tissues reveals several tissue-specific differential expressions of SAVs. This approach provides a novel access to high-throughput detection of protein variants, which may offer the potential for clinical biomarker discovery and mechanistic research.

The emerging role of RNA and DNA editing in cancer

Nucleotide sequence modification through single base editing in animals is emerging as an important player in tumorigenesis. RNA editing especially has increased greatly during mammalian evolution and modulates diverse cellular functions presumably in a context-dependent manner. Sequence editing impacts development, including pluripotency and hematopoiesis, and multiple recent studies have shown that dysregulation of editing is associated with tumor biology. Much is yet to be learned about the role of sequence editing in human biology but this process is a critical modulator of cell regulation and may present an attractive option for therapeutic intervention in cancer in the future.

SIGNIFICANCE

Sequence editing provides an additional regulatory layer of cancer initiation and progression that may be amenable to therapeutic design. Although editing of both RNA and DNA substrates has been known to occur for some time, the extent and implications of these modifications have been grossly underappreciated until recent genome-wide and disease-association studies were reported. This review highlights the cellular processes controlled by sequence editing, their implications in normal and cancerous states and considers potential targeted therapeutic strategies.

A biochemical landscape of A-to-I RNA editing in the human brain transcriptome

Inosine is an abundant RNA modification in the human transcriptome and is essential for many biological processes in modulating gene expression at the post-transcriptional level. Adenosine deaminases acting on RNA (ADARs) catalyze the hydrolytic deamination of adenosines to inosines (A-to-I editing) in double-stranded regions. We previously established a biochemical method called "inosine chemical erasing" (ICE) to directly identify inosines on RNA strands with high reliability. Here, we have applied the ICE method combined with deep sequencing (ICE-seq) to conduct an unbiased genome-wide screening of A-to-I editing sites in the transcriptome of human adult brain. Taken together with the sites identified by the conventional ICE method, we mapped 19,791 novel sites and newly found 1258 edited mRNAs, including 66 novel sites in coding regions, 41 of which cause altered amino acid assignment. ICE-seq detected novel editing sites in various repeat elements as well as in short hairpins. Gene ontology analysis revealed that these edited mRNAs are associated with transcription, energy metabolism, and neurological disorders, providing new insights into various aspects of human brain functions.

RNA-protein interactions: an overview

RNA binding proteins (RBPs) are key players in the regulation of gene expression. In this chapter we discuss the main protein-RNA recognition modes used by RBPs in order to regulate multiple steps of RNA processing. We discuss traditional and state-of-the-art technologies that can be used to study RNAs bound by individual RBPs, or vice versa, for both in vitro and in vivo methodologies. To help highlight the biological significance of RBP mediated regulation, online resources on experimentally verified protein-RNA interactions are briefly presented. Finally, we present the major tools to computationally infer RNA binding sites according to the modeling features and to the unsupervised or supervised frameworks that are adopted. Since some RNA binding site search algorithms are derived from DNA binding site search algorithms, we discuss the commonalities and novelties introduced to handle both sequence and structural features uniquely characterizing protein-RNA interactions.

Genome-wide analysis of A-to-I RNA editing by single-molecule sequencing in Drosophila

The accurate and thorough genome-wide detection of adenosine-to-inosine editing, a biologically indispensable process, has proven challenging. Here, we present a discovery pipeline in adult Drosophila, with 3,581 high-confidence editing sites identified with an estimated accuracy of 87%. The target genes and specific sites highlight global biological properties and functions of RNA editing, including hitherto-unknown editing in well-characterized classes of noncoding RNAs and 645 sites that cause amino acid substitutions, usually at conserved positions. The spectrum of functions that these gene targets encompass suggests that editing participates in a diverse set of cellular processes. Editing sites in Drosophila exhibit sequence-motif preferences and tend to be concentrated within a small subset of total RNAs. Finally, editing regulates expression levels of target mRNAs and strongly correlates with alternative splicing.

Tertiary structural elements determine the extent and specificity of messenger RNA editing

The specificity and extent of RNA editing by ADAR enzymes is determined largely by local primary sequence and secondary structural imperfections in duplex RNA. Here we surgically alter conserved cis elements associated with a cluster of ADAR modification sites within the endogenous Drosophila paralytic transcript. In addition to the local requirement for a central imperfect RNA duplex containing the modified adenosines, we demonstrate that a secondary RNA duplex containing splicing signals strongly modulates RNA editing. A subtle non-coding mutation, extending base pairing of this accessory helix, confers significant phenotypic consequences via effects on splicing. Through mutation/counter-mutation, we also uncover and functionally replace a highly conserved intronic long-range tertiary pseudoknot that is absolutely required for deamination of one particular adenosine in the central duplex. Our results demonstrate that complex RNA tertiary structures, which may be difficult to predict computationally, form in vivo and can regulate RNA-editing events.

Antagonistic and stimulative roles of ADAR1 in RNA silencing

Adenosine deaminases acting on RNA (ADARs) are involved in RNA editing that converts adenosine residues to inosine specifically in double-stranded RNAs (dsRNA). This A-to-I RNA editing pathway and the RNA interference (RNAi) pathway seem to interact antagonistically by competing for their common dsRNA substrates. For instance, A-to-I editing of certain microRNA (miRNA) precursors by ADAR1 and ADAR2 inhibits their processing to mature miRNAs. Recent studies unexpectedly revealed the presence of a completely different type of interaction between the RNA editing mechanism and the RNAi machinery. ADAR1 forms a complex via direct protein-protein interaction with Dicer, an RNase III gene family member involved in the RNAi mechanism. ADAR1 in the Dicer complex promotes pre-miRNA cleavage by Dicer and facilitates loading of miRNA onto RNA-induced silencing complexes, giving rise to an unsuspected stimulative function of ADAR1 on miRNA processing and RNAi mechanisms. ADAR1 differentiates its functions in RNA editing and RNAi by formation of either ADAR1-ADAR1 homodimer or Dicer-ADAR1 heterodimer complexes. Expression of miRNAs is globally inhibited in ADAR1-null mouse embryos, which, in turn, alters expression of their target genes and may contribute to their embryonic lethal phenotype.

Identification of the nature of reading frame transitions observed in prokaryotic genomes

Our goal was to identify evolutionary conserved frame transitions in protein coding regions and to uncover an underlying functional role of these structural aberrations. We used the ab initio frameshift prediction program, GeneTack, to detect reading frame transitions in 206 991 genes (fs-genes) from 1106 complete prokaryotic genomes. We grouped 102 731 fs-genes into 19 430 clusters based on sequence similarity between protein products (fs-proteins) as well as conservation of predicted position of the frameshift and its direction. We identified 4010 pseudogene clusters and 146 clusters of fs-genes apparently using recoding (local deviation from using standard genetic code) due to possessing specific sequence motifs near frameshift positions. Particularly interesting was finding of a novel type of organization of the dnaX gene, where recoding is required for synthesis of the longer subunit, τ. We selected 20 clusters of predicted recoding candidates and designed a series of genetic constructs with a reporter gene or affinity tag whose expression would require a frameshift event. Expression of the constructs in Escherichia coli demonstrated enrichment of the set of candidates with sequences that trigger genuine programmed ribosomal frameshifting; we have experimentally confirmed four new families of programmed frameshifts.

New insights from existing sequence data: generating breakthroughs without a pipette

With the rapidly declining cost of data generation and the accumulation of massive data sets, molecular biology is entering an era in which incisive analysis of existing data will play an increasingly prominent role in the discovery of new biological phenomena and the elucidation of molecular mechanisms. Here, we discuss resources of publicly available sequencing data most useful for interrogating the mechanisms of gene expression. Existing next-generation sequence data sets, however, come with significant challenges in the form of technical and bioinformatic artifacts, which we discuss in detail. We also recount several breakthroughs made largely through the analysis of existing data, primarily in the RNA field.

Reprogramming of microRNAs by adenosine-to-inosine editing and the selective elimination of edited microRNA precursors in mouse oocytes and preimplantation embryos

Adenosine deaminases-acting-on-RNA (ADAR) proteins induce adenosine-to-inosine editing in double-stranded RNA molecules. This editing generates RNA diversity at the post-transcriptional level, and it has been implicated in the control of cell differentiation and development. The editing of microRNA (miRNA) precursors, along with Tudor-SN (Snd1) activity, could lead to the elimination of selected miRNAs and reprogram miRNA activity. Here, we report the dynamics of adenosine-to-inosine editing in miRNA precursors and their selected elimination during mouse preimplantation development. Adar1p110 and Snd1 were found to be strongly but differentially expressed in oocytes and zygotes with respect to later pre-implantation stages. When the biogenesis of miR-151 was assessed, the majority of miR-151 precursors was edited and subsequently eliminated during early development. Deep sequencing of this and other miRNAs confirmed that, in general, edited precursors were selectively eliminated at early post-zygotic stages. Moreover, in oocytes and throughout the zygote-to-blastocyst stages, Tudor-SN accumulated in newly discovered aggregates termed ‘T bodies’. These results provide new insight into how editing and Tudor-SN-mediated elimination of miRNA precursors is regulated during early development.

Development of strategies for SNP detection in RNA-seq data: application to lymphoblastoid cell lines and evaluation using 1000 Genomes data

Next-generation RNA sequencing (RNA-seq) maps and analyzes transcriptomes and generates data on sequence variation in expressed genes. There are few reported studies on analysis strategies to maximize the yield of quality RNA-seq SNP data. We evaluated the performance of different SNP-calling methods following alignment to both genome and transcriptome by applying them to RNA-seq data from a HapMap lymphoblastoid cell line sample and comparing results with sequence variation data from 1000 Genomes. We determined that the best method to achieve high specificity and sensitivity, and greatest number of SNP calls, is to remove duplicate sequence reads after alignment to the genome and to call SNPs using SAMtools. The accuracy of SNP calls is dependent on sequence coverage available. In terms of specificity, 89% of RNA-seq SNPs calls were true variants where coverage is >10X. In terms of sensitivity, at >10X coverage 92% of all expected SNPs in expressed exons could be detected. Overall, the results indicate that RNA-seq SNP data are a very useful by-product of sequence-based transcriptome analysis. If RNA-seq is applied to disease tissue samples and assuming that genes carrying mutations relevant to disease biology are being expressed, a very high proportion of these mutations can be detected.

The ADAR protein family

Adenosine to inosine (A-to-I) RNA editing is a post-transcriptional process by which adenosines are selectively converted to inosines in double-stranded RNA (dsRNA) substrates. A highly conserved group of enzymes, the adenosine deaminase acting on RNA (ADAR) family, mediates this reaction. All ADARs share a common domain architecture consisting of a variable number of amino-terminal dsRNA binding domains (dsRBDs) and a carboxy-terminal catalytic deaminase domain. ADAR family members are highly expressed in the metazoan nervous system, where these enzymes predominantly localize to the neuronal nucleus. Once in the nucleus, ADARs participate in the modification of specific adenosines in pre-mRNAs of proteins involved in electrical and chemical neurotransmission, including pre-synaptic release machineries, and voltage- and ligand-gated ion channels. Most RNA editing sites in these nervous system targets result in non-synonymous codon changes in functionally important, usually conserved, residues and RNA editing deficiencies in various model organisms bear out a crucial role for ADARs in nervous system function. Mutation or deletion of ADAR genes results in striking phenotypes, including seizure episodes, extreme uncoordination, and neurodegeneration. Not only does the process of RNA editing alter important nervous system peptides, but ADARs also regulate gene expression through modification of dsRNA substrates that enter the RNA interference (RNAi) pathway and may then act at the chromatin level. Here, we present a review on the current knowledge regarding the ADAR protein family, including evolutionary history, key structural features, localization, function and mechanism.