Use of specific chemical reagents for detection of modified nucleotides in RNA

[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.

Mapping the tRNA modification landscape of Bartonella henselae Houston I and Bartonella quintana Toulouse

Transfer RNA (tRNA) modifications play a crucial role in maintaining translational fidelity and efficiency, and they may function as regulatory elements in stress response and virulence. Despite their pivotal roles, a comprehensive mapping of tRNA modifications and their associated synthesis genes is still limited, with a predominant focus on free-living bacteria. In this study, we employed a multidisciplinary approach, incorporating comparative genomics, mass spectrometry, and next-generation sequencing, to predict the set of tRNA modification genes responsible for tRNA maturation in two intracellular pathogens-Bartonella henselae Houston I and Bartonella quintana Toulouse, which are causative agents of cat-scratch disease and trench fever, respectively. This analysis presented challenges, particularly because of host RNA contamination, which served as a potential source of error. However, our approach predicted 26 genes responsible for synthesizing 23 distinct tRNA modifications in B. henselae and 22 genes associated with 23 modifications in B. quintana. Notably, akin to other intracellular and symbiotic bacteria, both Bartonella species have undergone substantial reductions in tRNA modification genes, mostly by simplifying the hypermodifications present at positions 34 and 37. Bartonella quintana exhibited the additional loss of four modifications and these were linked to examples of gene decay, providing snapshots of reductive evolution.

Mapping the tRNA Modification Landscape of Bartonella henselae Houston I and Bartonella quintana Toulouse

Transfer RNA (tRNA) modifications play a crucial role in maintaining translational fidelity and efficiency, and they may function as regulatory elements in stress response and virulence. Despite their pivotal roles, a comprehensive mapping of tRNA modifications and their associated synthesis genes is still limited, with a predominant focus on free-living bacteria. In this study, we employed a multidisciplinary approach, incorporating comparative genomics, mass spectrometry, and next-generation sequencing, to predict the set of tRNA modification genes responsible for tRNA maturation in two intracellular pathogens- Bartonella henselae Houston I and Bartonella quintana Toulouse, which are causative agents of cat-scratch disease and trench fever, respectively. This analysis presented challenges, particularly because of host RNA contamination, which served as a potential source of error. However, our approach predicted 26 genes responsible for synthesizing 23 distinct tRNA modifications in B. henselae and 22 genes associated with 23 modifications in B. quintana . Notably, akin to other intracellular and symbiotic bacteria, both Bartonella species have undergone substantial reductions in tRNA modification genes, mostly by simplifying the hypermodifications present at positions 34 and 37. B. quintana exhibited the additional loss of four modifications and these were linked to examples of gene decay, providing snapshots of reductive evolution.

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

Analysis approaches for the identification and prediction of N6-methyladenosine sites

The global dynamics in a variety of biological processes can be revealed by mapping transcriptional m⁶A sites, in particular full-transcriptome m⁶A. And individual m⁶A sites have contributed to biological function, which can be evaluated by stoichiometric information obtained from the single nucleotide resolution. Currently, the identification of m⁶A sites is mainly carried out by experiment and prediction methods, based on high-throughput sequencing and machine learning model respectively. This review summarizes the recent topics and progress made in bioinformatics methods of deciphering the m⁶A methylation, including the experimental detection of m⁶A methylation sites, techniques of data analysis, the way of predicting m⁶A methylation sites, m⁶A methylation databases, and detection of m⁶A modification in circRNA. At the end, the essay makes a brief discussion for the development perspective in this area.

Advantages and challenges associated with bisulfite-assisted nanopore direct RNA sequencing for modifications

Nanopore direct RNA sequencing is a technology that allows sequencing for epitranscriptomic modifications with the possibility of a quantitative assessment. In the present work, pseudouridine (Ψ) was sequenced with the nanopore before and after the pH 7 bisulfite reaction that yields stable ribose adducts at C1' of Ψ. The adducted sites produced greater base call errors in the form of deletion signatures compared to Ψ. Sequencing studies on E. coli rRNA and tmRNA before and after the pH 7 bisulfite reaction demonstrated that using chemically-assisted nanopore sequencing has distinct advantages for minimization of false positives and false negatives in the data. The rRNA from E. coli has 19 known U/C sequence variations that give similar base call signatures as Ψ, and therefore, are false positives when inspecting base call data; however, these sites are refractory to reacting with bisulfite as is easily observed in nanopore data. The E. coli tmRNA has a low occupancy Ψ in a pyrimidine-rich sequence context that is called a U representing a false negative; partial occupancy by Ψ is revealed after the bisulfite reaction. In a final study, 5-methylcytidine (m5C) in RNA can readily be observed after the pH 5 bisulfite reaction in which the parent C deaminates to U and the modified site does not react. This locates m5C when using bisulfite-assisted nanopore direct RNA sequencing, which is otherwise challenging to observe. The advantages and challenges of the overall approach are discussed.

Direct Nanopore Sequencing for the 17 RNA Modification Types in 36 Locations in the E. coli Ribosome Enables Monitoring of Stress-Dependent Changes

The bacterium Escherichia coli possesses 16S and 23S rRNA strands that have 36 chemical modification sites with 17 different structures. Nanopore direct RNA sequencing using a protein nanopore sensor and helicase brake, which is also a sensor, was applied to the rRNAs. Nanopore current levels, base calling profile, and helicase dwell times for the modifications relative to unmodified synthetic rRNA controls found signatures for nearly all modifications. Signatures for clustered modifications were determined by selective sequencing of writer knockout E. coli and sequencing of synthetic RNAs utilizing some custom-synthesized nucleotide triphosphates for their preparation. The knowledge of each modification's signature, apart from 5-methylcytidine, was used to determine how metabolic and cold-shock stress impact rRNA modifications. Metabolic stress resulted in either no change or a decrease, and one site increased in modification occupancy, while cold-shock stress led to either no change or a decrease. The double modification m4Cm1402 resides in 16S rRNA, and it decreased with both stressors. Using the helicase dwell time, it was determined that the N4 methyl group is lost during both stressors, and the 2'-OMe group remained. In the ribosome, this modification stabilizes binding to the mRNA codon at the P-site resulting in increased translational fidelity that is lost during stress. The E. coli genome has seven rRNA operons (rrn), and the earlier studies aligned the nanopore reads to a single operon (rrnA). Here, the reads were aligned to all seven operons to identify operon-specific changes in the 11 pseudouridines. This study demonstrates that direct sequencing for >16 different RNA modifications in a strand is achievable.

Oligodendrocyte differentiation alters tRNA modifications and codon optimality-mediated mRNA decay

Oligodendrocytes are specialized cells that confer neuronal myelination in the central nervous system. Leukodystrophies associated with oligodendrocyte deficits and hypomyelination are known to result when a number of tRNA metabolism genes are mutated. Thus, for unknown reasons, oligodendrocytes may be hypersensitive to perturbations in tRNA biology. In this study, we survey the tRNA transcriptome in the murine oligodendrocyte cell lineage and find that specific tRNAs are hypomodified in oligodendrocytes within or near the anticodon compared to oligodendrocyte progenitor cells (OPCs). This hypomodified state may be the result of differential expression of key modification enzymes during oligodendrocyte differentiation. Moreover, we observe a concomitant relationship between tRNA hypomodification and tRNA decoding potential; observing oligodendrocyte specific alterations in codon optimality-mediated mRNA decay and ribosome transit. Our results reveal that oligodendrocytes naturally maintain a delicate, hypersensitized tRNA/mRNA axis. We suggest this axis is a potential mediator of pathology in leukodystrophies and white matter disease when further insult to tRNA metabolism is introduced.

RNA methylation and cancer treatment

To this date, over 100 different types of RNA modification have been identified. Methylation of different RNA species has emerged as a critical regulator of transcript expression. RNA methylation and its related downstream signaling pathways are involved in plethora biological processes, including cell differentiation, sex determination and stress response, and others. It is catalyzed by the RNA methyltransferases, is demethylated by the demethylases (FTO and ALKBH5) and read by methylation binding protein (YTHDF1 and IGF2BP1). Increasing evidence indicates that this process closely connected to cancer cell proliferation, cellular stress, metastasis, immune response. And RNA methylation related protein has been becoming a promising targets of cancer therapy. This review outlines the relationship between different types of RNA methylation and cancer, and some FTO inhibitors in cancer treatment.

A chemical probe based on the PreQ1 metabolite enables transcriptome-wide mapping of binding sites

The role of metabolite-responsive riboswitches in regulating gene expression in bacteria is well known and makes them useful systems for the study of RNA-small molecule interactions. Here, we study the PreQ₁ riboswitch system, assessing sixteen diverse PreQ₁-derived probes for their ability to selectively modify the class-I PreQ₁ riboswitch aptamer covalently. For the most active probe (11), a diazirine-based photocrosslinking analog of PreQ₁, X-ray crystallography and gel-based competition assays demonstrated the mode of binding of the ligand to the aptamer, and functional assays demonstrated that the probe retains activity against the full riboswitch. Transcriptome-wide mapping using Chem-CLIP revealed a highly selective interaction between the bacterial aptamer and the probe. In addition, a small number of RNA targets in endogenous human transcripts were found to bind specifically to 11, providing evidence for candidate PreQ₁ aptamers in human RNA. This work demonstrates a stark influence of linker chemistry and structure on the ability of molecules to crosslink RNA, reveals that the PreQ₁ aptamer/ligand pair are broadly useful for chemical biology applications, and provides insights into how PreQ₁, which is similar in structure to guanine, interacts with human RNAs.

General Principles for the Detection of Modified Nucleotides in RNA by Specific Reagents

Epitranscriptomics heavily rely on chemical reagents for the detection, quantification, and localization of modified nucleotides in transcriptomes. Recent years have seen a surge in mapping methods that use innovative and rediscovered organic chemistry in high throughput approaches. While this has brought about a leap of progress in this young field, it has also become clear that the different chemistries feature variegated specificity and selectivity. The associated error rates, e.g., in terms of false positives and false negatives, are in large part inherent to the chemistry employed. This means that even assuming technically perfect execution, the interpretation of mapping results issuing from the application of such chemistries are limited by intrinsic features of chemical reactivity. An important but often ignored fact is that the huge stochiometric excess of unmodified over-modified nucleotides is not inert to any of the reagents employed. Consequently, any reaction aimed at chemical discrimination of modified versus unmodified nucleotides has optimal conditions for selectivity that are ultimately anchored in relative reaction rates, whose ratio imposes intrinsic limits to selectivity. Here chemical reactivities of canonical and modified ribonucleosides are revisited as a basis for an understanding of the limits of selectivity achievable with chemical methods.

Analysis of pseudouridines and other RNA modifications using HydraPsiSeq protocol

Detection of RNA modified nucleotides using deep sequencing can be performed by several approaches, including antibody-driven enrichment and natural or chemically induced RT signatures. However, only very few RNA modified nucleotides generate natural RT signatures and antibody-driven enrichment heavily depends on the quality of antibodies used and may be highly biased. Thus, the use of chemically-induced RT signatures is now considered as the most trusted experimental approach. In addition, the use of chemical reagents allows inclusion of simple "mock-treated" controls, to exclude spontaneous RT arrests, SNPs and other misincorporation-prone sites. Hydrazine is a well-known RNA-specific reagent, already extensively used in the past for RNA sequencing and structural probing. Hydrazine is highly reactive to U and shows low reaction rates with ψ residues, allowing their distinction by deep sequencing-based protocols. However, other modified RNA residues also show particular behavior upon hydrazine treatment. Here we present methodological developments allowing to use HydraPsiSeq for precise quantification of RNA pseudouridylation and also detection and quantification of some other RNA modifications, in addition to ψ residues.

Arrow pushing in RNA modification sequencing

Methods to accurately determine the location and abundance of RNA modifications are critical to understanding their functional role. In this review, we describe recent efforts in which chemical reactivity and next-generation sequencing have been integrated to detect modified nucleotides in RNA. For eleven exemplary modifications, we detail chemical, enzymatic, and metabolic labeling protocols that can be used to differentiate them from canonical nucleobases. By emphasizing the molecular rationale underlying these detection methods, our survey highlights new opportunities for chemistry to define the role of RNA modifications in disease.

Mapping of 7-methylguanosine (m7G), 3-methylcytidine (m3C), dihydrouridine (D) and 5-hydroxycytidine (ho5C) RNA modifications by AlkAniline-Seq

Precise and reliable mapping of modified nucleotides in RNA is a challenging task in epitranscriptomics analysis. Only deep sequencing-based methods are able to provide both, a single-nucleotide resolution and sufficient selectivity and sensitivity. A number of protocols employing specific chemical reagents to distinguish modified RNA nucleotides from canonical parental residues have already proven their performance. We developed a deep-sequencing analytical pipeline for simultaneous detection of several modified nucleotides of different nature (methylation, hydroxylation, reduction) in RNA. The AlkAniline-Seq protocol uses intrinsic fragility of the N-glycosidic bond present in certain modified residues (7-methylguanosine (m⁷G), 3-methylcytidine (m³C), dihydrouridine (D) and 5-hydroxycytidine (ho⁵C)) to induce cleavage under heat combined with alkaline conditions. The resulting RNA abasic site is decomposed by aniline-driven β-elimination and creates a 5'-phosphate (5'-P) at the adjacent N+1 residue. This 5'-P is the crucial entry point for a highly selective ligation of sequencing adapters during the subsequent Illumina library preparation protocol. AlkAniline-Seq protocol has a very low background, and is both highly sensitive and specific. Applications of AlkAniline-Seq include mapping of m⁷G, m³C, D, and ho⁵C in variety of cellular RNAs, including in particular rRNAs and tRNAs.

RNA Secondary Structure Study by Chemical Probing Methods Using DMS and CMCT

RNA folds into secondary structures that can serve in understanding various RNA functions (Weeks KM. Curr Opin Struct Biol 20(3):295-304, 2010). Chemical probing is a method that enables the characterization of RNA secondary structures using chemical reagents that specifically modify RNA nucleotides that are located in single-stranded areas. In our protocol, we used Dimethyl Sulfate (DMS) and Cyclohexyl-3-(2-Morpholinoethyl) Carbodiimide metho-p-Toluene sulfonate (CMCT) that are both base-specific modifying reagents (Behm-Ansmant I, et al. J Nucleic Acids 2011:408053, 2011). These modifications are mapped by primer extension arrests using 5' fluorescently labeled primers. In this protocol, we show a comprehensive method to identify RNA secondary structures in vitro using fluorescently labeled oligos. To demonstrate the efficiency of the method, we give an example of a structure we have designed which corresponds to a part of the 5'-UTR regulatory element called Translation Inhibitory Element (TIE) from Hox a3 mRNA (Xue S, et al. Nature 517(7532):33-38, 2015).

RNA Structure Analysis by Chemical Probing with DMS and CMCT

RNA structure is important for understanding RNA function and stability within a cell. Chemical probing is a well-established and convenient method to evaluate the structure of an RNA. Several structure-sensitive chemicals can differentiate paired and unpaired nucleotides. This chapter specifically addresses the use of DMS and CMCT. Although exhibiting different affinities, the combination of these two chemical reagents enables screening of all four nucleobases. DMS and CMCT are only reactive with exposed unpaired nucleotides. We have used this method to analyze the effect of the RNA chaperone Hfq on the conformation of the 16S rRNA. The strategy here described may be applied for the study of many other RNA-binding proteins and RNAs.

Changes of the tRNA Modification Pattern during the Development of Dictyostelium discoideum

Dictyostelium discoideum is a social amoeba, which on starvation develops from a single-cell state to a multicellular fruiting body. This developmental process is accompanied by massive changes in gene expression, which also affect non-coding RNAs. Here, we investigate how tRNAs as key regulators of the translation process are affected by this transition. To this end, we used LOTTE-seq to sequence the tRNA pool of D. discoideum at different developmental time points and analyzed both tRNA composition and tRNA modification patterns. We developed a workflow for the specific detection of modifications from reverse transcriptase signatures in chemically untreated RNA-seq data at single-nucleotide resolution. It avoids the comparison of treated and untreated RNA-seq data using reverse transcription arrest patterns at nucleotides in the neighborhood of a putative modification site as internal control. We find that nucleotide modification sites in D. discoideum tRNAs largely conform to the modification patterns observed throughout the eukaroytes. However, there are also previously undescribed modification sites. We observe substantial dynamic changes of both expression levels and modification patterns of certain tRNA types during fruiting body development. Beyond the specific application to D. discoideum our results demonstrate that the developmental variability of tRNA expression and modification can be traced efficiently with LOTTE-seq.

DeepOMe: A Web Server for the Prediction of 2'-O-Me Sites Based on the Hybrid CNN and BLSTM Architecture

2'-O-methylations (2'-O-Me or Nm) are one of the most important layers of regulatory control over gene expression. With increasing attentions focused on the characteristics, mechanisms and influences of 2'-O-Me, a revolutionary technique termed Nm-seq were established, allowing the identification of precise 2'-O-Me sites in RNA sequences with high sensitivity. However, as the costs and complexities involved with this new method, the large-scale detection and in-depth study of 2'-O-Me is still largely limited. Therefore, the development of a novel computational method to identify 2'-O-Me sites with adequate reliability is urgently needed at the current stage. To address the above issue, we proposed a hybrid deep-learning algorithm named DeepOMe that combined Convolutional Neural Networks (CNN) and Bidirectional Long Short-term Memory (BLSTM) to accurately predict 2'-O-Me sites in human transcriptome. Validating under 4-, 6-, 8-, and 10-fold cross-validation, we confirmed that our proposed model achieved a high performance (AUC close to 0.998 and AUPR close to 0.880). When testing in the independent data set, DeepOMe was substantially superior to NmSEER V2.0. To facilitate the usage of DeepOMe, a user-friendly web-server was constructed, which can be freely accessed at http://deepome.renlab.org.

Nucleotide resolution profiling of m3C RNA modification by HAC-seq

Cellular RNAs are subject to a myriad of different chemical modifications that play important roles in controlling RNA expression and function. Dysregulation of certain RNA modifications, the so-called 'epitranscriptome', contributes to human disease. One limitation in studying the functional, physiological, and pathological roles of the epitranscriptome is the availability of methods for the precise mapping of individual RNA modifications throughout the transcriptome. 3-Methylcytidine (m3C) modification of certain tRNAs is well established and was also recently detected in mRNA. However, methods for the specific mapping of m3C throughout the transcriptome are lacking. Here, we developed a m3C-specific technique, Hydrazine-Aniline Cleavage sequencing (HAC-seq), to profile the m3C methylome at single-nucleotide resolution. We applied HAC-seq to analyze ribosomal RNA (rRNA)-depleted total RNAs in human cells. We found that tRNAs are the predominant m3C-modified RNA species, with 17 m3C modification sites on 11 cytoplasmic and 2 mitochondrial tRNA isoacceptors in MCF7 cells. We found no evidence for m3C-modification of mRNA or other non-coding RNAs at comparable levels to tRNAs in these cells. HAC-seq provides a novel method for the unbiased, transcriptome-wide identification of m3C RNA modification at single-nucleotide resolution, and could be widely applied to reveal the m3C methylome in different cells and tissues.

Instrumental analysis of RNA modifications

Organisms from all domains of life invest a substantial amount of energy for the introduction of RNA modifications into nearly all transcripts studied to date. Instrumental analysis of RNA can focus on the modified residues and reveal the function of these epitranscriptomic marks. Here, we will review recent advances and breakthroughs achieved by NMR spectroscopy, sequencing, and mass spectrometry of the epitranscriptome.

Analysis of RNA Modifications by Second- and Third-Generation Deep Sequencing: 2020 Update

The precise mapping and quantification of the numerous RNA modifications that are present in tRNAs, rRNAs, ncRNAs/miRNAs, and mRNAs remain a major challenge and a top priority of the epitranscriptomics field. After the keystone discoveries of massive m⁶A methylation in mRNAs, dozens of deep sequencing-based methods and protocols were proposed for the analysis of various RNA modifications, allowing us to considerably extend the list of detectable modified residues. Many of the currently used methods rely on the particular reverse transcription signatures left by RNA modifications in cDNA; these signatures may be naturally present or induced by an appropriate enzymatic or chemical treatment. The newest approaches also include labeling at RNA abasic sites that result from the selective removal of RNA modification or the enhanced cleavage of the RNA ribose-phosphate chain (perhaps also protection from cleavage), followed by specific adapter ligation. Classical affinity/immunoprecipitation-based protocols use either antibodies against modified RNA bases or proteins/enzymes, recognizing RNA modifications. In this survey, we review the most recent achievements in this highly dynamic field, including promising attempts to map RNA modifications by the direct single-molecule sequencing of RNA by nanopores.

An epigenetic 'extreme makeover': the methylation of flaviviral RNA (and beyond)

Beyond their high clinical relevance worldwide, flaviviruses (comprising dengue and Zika viruses) are of particular interest to understand the spatiotemporal control of RNA metabolism. Indeed, their positive single-stranded viral RNA genome (vRNA) undergoes in the cytoplasm replication, translation and encapsidation, three steps of the flavivirus life cycle that are coordinated through a fine-tuned equilibrium. Over the last years, RNA methylation has emerged as a powerful mechanism to regulate messenger RNA metabolism at the posttranscriptional level. Not surprisingly, flaviviruses exploit RNA epigenetic strategies to control crucial steps of their replication cycle as well as to evade sensing by the innate immune system. This review summarizes the current knowledge about vRNA methylation events and their impacts on flavivirus replication and pathogenesis. We also address the important challenges that the field of epitranscriptomics faces in reliably and accurately identifying RNA methylation sites, which should be considered in future studies on viral RNA modifications.

Deciphering RNA modifications at base resolution: from chemistry to biology

Nearly 200 distinct chemical modifications of RNAs have been discovered to date. Their analysis via direct methods has been possible in abundant RNA species, such as ribosomal, transfer or viral RNA, since several decades. However, their analysis in less abundant RNAs species, especially cellular messenger RNAs, was rendered possible only recently with the advent of high throughput sequencing techniques. Given the growing biomedical interest of the proteins that write, erase and read RNA modifications, ingenious new methods to enrich and identify RNA modifications at base resolution have been implemented, and more efforts are underway to render them more quantitative. Here, we review several crucial modification-specific (bio)chemical approaches and discuss their advantages and shortcomings for exploring the epitranscriptome.

Detection of Modified Bases in Bacteriophage Genomic DNA

Collectively, the dsDNA tailed bacteriophages (Caudovirales) contain the largest chemical diversity of naturally occurring deoxynucleotides in DNA observed to date. The continuing discovery of new modifications in phages suggest many more are waiting to be found. Thus, methods for the observation and characterization of noncanonical nucleosides are timely. We present here protocols for extraction of genomic DNA from bacteriophage particles, enzymatic hydrolysis of DNA to free nucleosides, and examination of nucleoside composition by HPLC and mass spectrometry.

HydraPsiSeq: a method for systematic and quantitative mapping of pseudouridines in RNA

Developing methods for accurate detection of RNA modifications remains a major challenge in epitranscriptomics. Next-generation sequencing-based mapping approaches have recently emerged but, often, they are not quantitative and lack specificity. Pseudouridine (ψ), produced by uridine isomerization, is one of the most abundant RNA modification. ψ mapping classically involves derivatization with soluble carbodiimide (CMCT), which is prone to variation making this approach only semi-quantitative. Here, we developed 'HydraPsiSeq', a novel quantitative ψ mapping technique relying on specific protection from hydrazine/aniline cleavage. HydraPsiSeq is quantitative because the obtained signal directly reflects pseudouridine level. Furthermore, normalization to natural unmodified RNA and/or to synthetic in vitro transcripts allows absolute measurements of modification levels. HydraPsiSeq requires minute amounts of RNA (as low as 10-50 ng), making it compatible with high-throughput profiling of diverse biological and clinical samples. Exploring the potential of HydraPsiSeq, we profiled human rRNAs, revealing strong variations in pseudouridylation levels at ∼20-25 positions out of total 104 sites. We also observed the dynamics of rRNA pseudouridylation throughout chondrogenic differentiation of human bone marrow stem cells. In conclusion, HydraPsiSeq is a robust approach for the systematic mapping and accurate quantification of pseudouridines in RNAs with applications in disease, aging, development, differentiation and/or stress response.

Survey and Validation of tRNA Modifications and Their Corresponding Genes in Bacillus subtilis sp Subtilis Strain 168

Extensive knowledge of both the nature and position of tRNA modifications in all cellular tRNAs has been limited to two bacteria, Escherichia coli and Mycoplasma capricolum. Bacillus subtilis sp subtilis strain 168 is the model Gram-positive bacteria and the list of the genes involved in tRNA modifications in this organism is far from complete. Mass spectrometry analysis of bulk tRNA extracted from B. subtilis, combined with next generation sequencing technologies and comparative genomic analyses, led to the identification of 41 tRNA modification genes with associated confidence scores. Many differences were found in this model Gram-positive bacteria when compared to E. coli. In general, B. subtilis tRNAs are less modified than those in E. coli, even if some modifications, such as m1A22 or ms2t6A, are only found in the model Gram-positive bacteria. Many examples of non-orthologous displacements and of variations in the most complex pathways are described. Paralog issues make uncertain direct annotation transfer from E. coli to B. subtilis based on homology only without further experimental validation. This difficulty was shown with the identification of the B. subtilis enzyme that introduces ψ at positions 31/32 of the tRNAs. This work presents the most up to date list of tRNA modification genes in B. subtilis, identifies the gaps in knowledge, and lays the foundation for further work to decipher the physiological role of tRNA modifications in this important model organism and other bacteria.

Detection of internal N7-methylguanosine (m7G) RNA modifications by mutational profiling sequencing

Methylation of guanosine on position N7 (m⁷G) on internal RNA positions has been found in all domains of life and have been implicated in human disease. Here, we present m⁷G Mutational Profiling sequencing (m⁷G-MaP-seq), which allows high throughput detection of m⁷G modifications at nucleotide resolution. In our method, m⁷G modified positions are converted to abasic sites by reduction with sodium borohydride, directly recorded as cDNA mutations through reverse transcription and sequenced. We detect positions with increased mutation rates in the reduced and control samples taking the possibility of sequencing/alignment error into account and use replicates to calculate statistical significance based on log likelihood ratio tests. We show that m⁷G-MaP-seq efficiently detects known m⁷G modifications in rRNA with mutational rates up to 25% and we map a previously uncharacterised evolutionarily conserved rRNA modification at position 1581 in Arabidopsis thaliana SSU rRNA. Furthermore, we identify m⁷G modifications in budding yeast, human and arabidopsis tRNAs and demonstrate that m⁷G modification occurs before tRNA splicing. We do not find any evidence for internal m⁷G modifications being present in other small RNA, such as miRNA, snoRNA and sRNA, including human Let-7e. Likewise, high sequencing depth m⁷G-MaP-seq analysis of mRNA from E. coli or yeast cells did not identify any internal m⁷G modifications.

High-resolution ion mobility spectrometry-mass spectrometry of isomeric/isobaric ribonucleotide variants

In this report, we explored the benefits of cyclic ion mobility (cIM) mass spectrometry in the analysis of isomeric post-transcriptional modifications of RNA. Standard methyl-cytidine samples were initially utilized to test the ability to correctly distinguish different structures sharing the same elemental composition and thus molecular mass. Analyzed individually, the analytes displayed characteristic arrival times (tD ) determined by the different positions of the modifying methyl groups onto the common cytidine scaffold. Analyzed in mixture, the widths of the respective signals resulted in significant overlap that initially prevented their resolution on the tD scale. The separation of the four isomers was achieved by increasing the number of passes through the cIM device, which enabled to fully differentiate the characteristic ion mobility behaviors associated with very subtle structural variations. The placement of the cIM device between the mass-selective quadrupole and the time-of-flight analyzer allowed us to perform gas-phase activation of each of these ion populations, which had been first isolated according to a common mass-to-charge ratio and then separated on the basis of different ion mobility behaviors. The observed fragmentation patterns confirmed the structures of the various isomers thus substantiating the benefits of complementing unique tD information with specific fragmentation data to reach more stringent analyte identification. These capabilities were further tested by analyzing natural mono-nucleotide mixtures obtained by exonuclease digestion of total RNA extracts. In particular, the combination of cIM separation and post-mobility dissociation allowed us to establish the composition of methyl-cytidine and methyl-adenine components present in the entire transcriptome of HeLa cells. For this reason, we expect that this technique will benefit not only epitranscriptomic studies requiring the determination of identity and expression levels of RNA modifications, but also metabolomics investigations involving the analysis of natural extracts that may possibly contain subsets of isomeric/isobaric species.

Chemical tagging for sensitive determination of uridine modifications in RNA

The discovery of dynamic and reversible modifications in messenger RNA (mRNA) is opening new directions in RNA modification-mediated regulation of biological processes. Methylation is the most prevalent modification occurring in mRNA and the methyl group is mainly decorated in the adenine, cytosine, and guanine base or in the 2′-hydroxyl group of ribose. However, methylation of the uracil base (5-methyluridine, m⁵U) has not been discovered in mRNA of eukaryotes. In the current study, we established a method of N-cyclohexyl-N′-β-(4-methylmorpholinium) ethylcarbodiimide p-toluenesulfonate (CMCT) labelling coupled with liquid chromatography-electrospray ionization-mass spectrometry (LC-ESI-MS/MS) analysis for the sensitive determination of uridine modifications in RNA. Our results demonstrated that the detection sensitivities of uridine modifications in RNA increased up to 1408 fold upon CMCT labelling. Using the developed method, we identified the distinct existence of m⁵U in mRNA of various mammalian cells and tissues. In addition, the stable isotope tracing monitored by mass spectrometry revealed that the methyl group of m⁵U originated from S-adenosyl-l-methionine (SAM). Our study expanded the list of modifications occurring in mRNA of mammals. Future work on transcriptome-wide mapping of m⁵U will further uncover the functional roles of m⁵U in mRNA of mammals.

The discovery of dynamic and reversible modifications in messenger RNA is opening new directions in RNA modification-mediated regulation of biological processes.

Distinct methylation levels of mature microRNAs in gastrointestinal cancers

The biological significance of micro (mi)RNAs has traditionally been evaluated according to their RNA expression levels based on the assumption that miRNAs recognize and regulate their targets in an unvarying fashion. Here we show that a fraction of mature miRNAs including miR-17-5p, -21-5p, and -200c-3p and let-7a-5p harbor methyl marks that potentially alter their stability and target recognition. Importantly, methylation of these miRNAs was significantly increased in cancer tissues as compared to paired normal tissues. Furthermore, miR-17-5p methylation level in serum samples distinguished early pancreatic cancer patients from healthy controls with extremely high sensitivity and specificity. These findings provide a basis for diagnostic strategies for early-stage cancer and add a dimension to our understanding of miRNA biology.

In cancer it is assumed that microRNAs recognise and regulate their targets uniformly. Here, the authors show that in gastrointestinal cancers methylation of microRNAs may impact their stability, and that levels of microRNA methylation are distinct in pancreatic cancer patients compared to healthy controls with potential diagnostic implications.

RNA Framework: an all-in-one toolkit for the analysis of RNA structures and post-transcriptional modifications

RNA is emerging as a key regulator of a plethora of biological processes. While its study has remained elusive for decades, the recent advent of high-throughput sequencing technologies provided the unique opportunity to develop novel techniques for the study of RNA structure and post-transcriptional modifications. Nonetheless, most of the required downstream bioinformatics analyses steps are not easily reproducible, thus making the application of these techniques a prerogative of few laboratories. Here we introduce RNA Framework, an all-in-one toolkit for the analysis of most NGS-based RNA structure probing and post-transcriptional modification mapping experiments. To prove the extreme versatility of RNA Framework, we applied it to both an in-house generated DMS-MaPseq dataset, and to a series of literature available experiments. Notably, when starting from publicly available datasets, our software easily allows replicating authors' findings. Collectively, RNA Framework provides the most complete and versatile toolkit to date for a rapid and streamlined analysis of the RNA epistructurome. RNA Framework is available for download at: http://www.rnaframework.com.

Transcriptome-wide profiling of multiple RNA modifications simultaneously at single-base resolution

The field of RNA modification would be significantly advanced by the development of sensitive, accurate, single-base resolution methods for profiling multiple common RNA modifications in the same RNA molecule. Our work provides several advances toward that goal, including (i) quantitative methods for profiling Ψ sites at true base-pair resolution transcriptome-wide, (ii) a chemical understanding of our observed Ψ-dependent deletion signature, (iii) improved methods for profiling m⁵C and m¹A, and (iv) a coupling of these methods for the simultaneous detection of all three modifications in the same RNA. Together, the combinatorial ability and relative ease of execution provided by this procedure should greatly forward epitranscriptome studies involving these three very common RNA modifications.

The breadth and importance of RNA modifications are growing rapidly as modified ribonucleotides can impact the sequence, structure, function, stability, and fate of RNAs and their interactions with other molecules. Therefore, knowing cellular RNA modifications at single-base resolution could provide important information regarding cell status and fate. A current major limitation is the lack of methods that allow the reproducible profiling of multiple modifications simultaneously, transcriptome-wide and at single-base resolution. Here we developed RBS-Seq, a modification of RNA bisulfite sequencing that enables the sensitive and simultaneous detection of m⁵C, Ψ, and m¹A at single-base resolution transcriptome-wide. With RBS-Seq, m⁵C and m¹A are accurately detected based on known signature base mismatches and are detected here simultaneously along with Ψ sites that show a 1–2 base deletion. Structural analyses revealed the mechanism underlying the deletion signature, which involves Ψ-monobisulfite adduction, heat-induced ribose ring opening, and Mg²⁺-assisted reorientation, causing base-skipping during cDNA synthesis. Detection of each of these modifications through a unique chemistry allows high-precision mapping of all three modifications within the same RNA molecule, enabling covariation studies. Application of RBS-Seq on HeLa RNA revealed almost all known m⁵C, m¹A, and ψ sites in tRNAs and rRNAs and provided hundreds of new m⁵C and Ψ sites in noncoding RNAs and mRNAs. However, our results diverge greatly from earlier work, suggesting ∼10-fold fewer m⁵C sites in noncoding and coding RNAs and the absence of substantial m¹A in mRNAs. Taken together, the approaches and refined datasets in this work will greatly enable future epitranscriptome studies.

N1-Methyladenosine detection with CRISPR-Cas13a/C2c2

Recent studies suggested that the widespread presence of N1-methyladenosine (m¹A) plays a very important role in environmental stress, ribosome biogenesis and antibiotic resistance. The RNA-guided, RNA-targeting CRISPR Cas13a exhibits a "collateral effect" of promiscuous RNase activity upon target recognition with high sensitivity. Inspired by the advantage of CRISPR Cas13a, we designed a system to detect m¹A induced mismatch, providing a rapid, simple and fluorescence-based m¹A detection. For A-ssRNA, the Cas13a-based molecular detection platform showed a high fluorescence signal. For m¹A-ssRNA, there is an about 90% decline of fluorescence. Moreover, this approach can also be used to quantify m¹A in RNAs and applied for the analysis of dynamic m¹A demethylation of 28S rRNA with AlkB.

Methods for RNA Modification Mapping Using Deep Sequencing: Established and New Emerging Technologies

New analytics of post-transcriptional RNA modifications have paved the way for a tremendous upswing of the biological and biomedical research in this field. This especially applies to methods that included RNA-Seq techniques, and which typically result in what is termed global scale modification mapping. In this process, positions inside a cell’s transcriptome are receiving a status of potential modification sites (so called modification calling), typically based on a score of some kind that issues from the particular method applied. The resulting data are thought to represent information that goes beyond what is contained in typical transcriptome data, and hence the field has taken to use the term “epitranscriptome”. Due to the high rate of newly published mapping techniques, a significant number of chemically distinct RNA modifications have become amenable to mapping, albeit with variegated accuracy and precision, depending on the nature of the technique. This review gives a brief overview of known techniques, and how they were applied to modification calling.

Post-transcriptional regulation by cytosine-5 methylation of RNA

The recent advent of high-throughput sequencing technologies coupled with RNA modifications detection methods has allowed the detection of RNA modifications at single nucleotide resolution giving a more comprehensive landscape of post-transcriptional gene regulation pathways. In this review, we focus on the occurrence of 5-methylcytosine (m⁵C) in the transcriptome. We summarise the main findings of the molecular role in post-transcriptional regulation that governs m⁵C deposition in RNAs. Functionally, m⁵C deposition can regulate several cellular and physiological processes including development, differentiation and survival to stress stimuli. Despite many aspects concerning m⁵C deposition in RNA, such as position or sequence context and the fact that many readers and erasers still remain elusive, the overall recent findings indicate that RNA cytosine methylation is a powerful mechanism to post-transcriptionally regulate physiological processes. In addition, mutations in RNA cytosine-5 methyltransferases are associated to pathological processes ranging from neurological syndromes to cancer.

7-Methylguanosine Modifications in Transfer RNA (tRNA)

More than 90 different modified nucleosides have been identified in tRNA. Among the tRNA modifications, the 7-methylguanosine (m⁷G) modification is found widely in eubacteria, eukaryotes, and a few archaea. In most cases, the m⁷G modification occurs at position 46 in the variable region and is a product of tRNA (m⁷G46) methyltransferase. The m⁷G46 modification forms a tertiary base pair with C13-G22, and stabilizes the tRNA structure. A reaction mechanism for eubacterial tRNA m⁷G methyltransferase has been proposed based on the results of biochemical, bioinformatic, and structural studies. However, an experimentally determined mechanism of methyl-transfer remains to be ascertained. The physiological functions of m⁷G46 in tRNA have started to be determined over the past decade. For example, tRNA m⁷G46 or tRNA (m⁷G46) methyltransferase controls the amount of other tRNA modifications in thermophilic bacteria, contributes to the pathogenic infectivity, and is also associated with several diseases. In this review, information of tRNA m⁷G modifications and tRNA m⁷G methyltransferases is summarized and the differences in reaction mechanism between tRNA m⁷G methyltransferase and rRNA or mRNA m⁷G methylation enzyme are discussed.

An In Vitro Assay to Detect tRNA-Isopentenyl Transferase Activity

N⁶-isopentenyladenosine RNA modifications are functionally diverse and highly conserved among prokaryotes and eukaryotes. One of the most highly conserved N⁶-isopentenyladenosine modifications occurs at the A37 position in a subset of tRNAs. This modification improves translation efficiency and fidelity by increasing the affinity of the tRNA for the ribosome. Mutation of enzymes responsible for this modification in eukaryotes are associated with several disease states, including mitochondrial dysfunction and cancer. Therefore, understanding the substrate specificity and biochemical activities of these enzymes is important for understanding of normal and pathologic eukaryotic biology. A diverse array of methods has been employed to characterize i⁶A modifications. Herein is described a direct approach for the detection of isopentenylation by Mod5. This method utilizes incubation of RNAs with a recombinant isopentenyl transferase, followed by RNase T1 digestion, and 1-dimensional gel electrophoresis analysis to detect i⁶A modifications. In addition, the potential adaptability of this protocol to characterize other RNA-modifying enzymes is discussed.

Reversible RNA Modification N1-methyladenosine (m1A) in mRNA and tRNA

More than 100 modifications have been found in RNA. Analogous to epigenetic DNA methylation, epitranscriptomic modifications can be written, read, and erased by a complex network of proteins. Apart from N6-methyladenosine (m6A), N1-methyladenosine (m1A) has been found as a reversible modification in tRNA and mRNA. m1A occurs at positions 9, 14, and 58 of tRNA, with m1A58 being critical for tRNA stability. Other than the hundreds of m1A sites in mRNA and long non-coding RNA transcripts, transcriptome-wide mapping of m1A also identifies >20 m1A sites in mitochondrial genes. m1A in the coding region of mitochondrial transcripts can inhibit the translation of the corresponding proteins. In this review, we summarize the current understanding of m1A in mRNA and tRNA, covering high-throughput sequencing methods developed for m1A methylome, m1A-related enzymes (writers and erasers), as well as its functions in mRNA and tRNA.

Monitoring the 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2U) modification in eukaryotic tRNAs via the γ-toxin endonuclease

The post-transcriptional modification of tRNA at the wobble position is a universal process occurring in all domains of life. In eukaryotes, the wobble uridine of particular tRNAs is transformed to the 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2U) modification which is critical for proper mRNA decoding and protein translation. However, current methods to detect mcm5s2U are technically challenging and/or require specialized instrumental expertise. Here, we show that γ-toxin endonuclease from the yeast Kluyveromyces lactis can be used as a probe for assaying mcm5s2U status in the tRNA of diverse eukaryotic organisms ranging from protozoans to mammalian cells. The assay couples the mcm5s2U-dependent cleavage of tRNA by γ-toxin with standard molecular biology techniques such as northern blot analysis or quantitative PCR to monitor mcm5s2U levels in multiple tRNA isoacceptors. The results gained from the γ-toxin assay reveals the evolutionary conservation of the mcm5s2U modification across eukaryotic species. Moreover, we have used the γ-toxin assay to verify uncharacterized eukaryotic Trm9 and Trm112 homologs that catalyze the formation of mcm5s2U. These findings demonstrate the use of γ-toxin as a detection method to monitor mcm5s2U status in diverse eukaryotic cell types for cellular, genetic, and biochemical studies.

Understanding RNA modifications: the promises and technological bottlenecks of the 'epitranscriptome'

The discovery of mechanisms that alter genetic information via RNA editing or introducing covalent RNA modifications points towards a complexity in gene expression that challenges long-standing concepts. Understanding the biology of RNA modifications represents one of the next frontiers in molecular biology. To this date, over 130 different RNA modifications have been identified, and improved mass spectrometry approaches are still adding to this list. However, only recently has it been possible to map selected RNA modifications at single-nucleotide resolution, which has created a number of exciting hypotheses about the biological function of RNA modifications, culminating in the proposition of the ‘epitranscriptome’. Here, we review some of the technological advances in this rapidly developing field, identify the conceptual challenges and discuss approaches that are needed to rigorously test the biological function of specific RNA modifications.

In Silico Identification of RNA Modifications from High-Throughput Sequencing Data Using HAMR

RNA molecules are often altered post-transcriptionally by the covalent modification of their nucleotides. These modifications are known to modulate the structure, function, and activity of RNAs. When reverse transcribed into cDNA during RNA sequencing library preparation, atypical (modified) ribonucleotides that affect Watson-Crick base pairing will interfere with reverse transcriptase (RT), resulting in cDNA products with mis-incorporated bases or prematurely terminated RNA products. These interactions with RT can therefore be inferred from mismatch patterns in the sequencing reads, and are distinguishable from simple base-calling errors, single-nucleotide polymorphisms (SNPs), or RNA editing sites. Here, we describe a computational protocol for the in silico identification of modified ribonucleotides from RT-based RNA-seq read-out using the High-throughput Analysis of Modified Ribonucleotides (HAMR) software. HAMR can identify these modifications transcriptome-wide with single nucleotide resolution, and also differentiate between different types of modifications to predict modification identity. Researchers can use HAMR to identify and characterize RNA modifications using RNA-seq data from a variety of common RT-based sequencing protocols such as Poly(A), total RNA-seq, and small RNA-seq.

Accurate detection of chemical modifications in RNA by mutational profiling (MaP) with ShapeMapper 2

Mutational profiling (MaP) enables detection of sites of chemical modification in RNA as sequence changes during reverse transcription (RT), subsequently read out by massively parallel sequencing. We introduce ShapeMapper 2, which integrates careful handling of all classes of adduct-induced sequence changes, sequence variant correction, basecall quality filters, and quality-control warnings to now identify RNA adduct sites as accurately as achieved by careful manual analysis of electrophoresis data, the prior highest-accuracy standard. MaP and ShapeMapper 2 provide a robust, experimentally concise, and accurate approach for reading out nucleic acid chemical probing experiments.

Modern Approaches for Identification of Modified Nucleotides in RNA

This review considers approaches for detection of modified monomers in the RNA structure of living organisms. Recently, some data on dynamic alterations in the pool of modifications of the key RNA species that depend on external factors affecting the cells and physiological conditions of the whole organism have been accumulated. The recent studies have presented experimental data on relationship between the mechanisms of formation of modified/minor nucleotides of RNA in mammalian cells and the development of various pathologies. The development of novel methods for detection of chemical modifications of RNA nucleotides in the cells of living organisms and accumulation of knowledge on the contribution of modified monomers to metabolism and functioning of individual RNA species establish the basis for creation of novel diagnostic and therapeutic approaches. This review includes a short description of routine methods for determination of modified nucleotides in RNA and considers in detail modern approaches that enable not only detection but also quantitative assessment of the modification level of various nucleotides in individual RNA species.

Northern Blot Analysis of Circular RNAs

Northern blotting enables the specific detection and characterization of RNA molecules. Recently, circular RNAs (circRNAs) were described as a new class of cell type-specific noncoding RNAs. With the discovery of many novel circRNAs on the basis of high-throughput sequencing and bioinformatics, a solid biochemical approach is required to directly detect and validate specific circRNA species. Here we give a detailed overview of how different Northern blot methods can be employed to validate specific circRNAs. Different Northern gel and detection systems are introduced, in combination with additional tools for circRNA characterization, such as RNase R and RNase H treatments.

High-throughput single-base resolution mapping of RNA 2΄-O-methylated residues

Functional characterization of the transcriptome requires tools for the systematic investigation of RNA post-transcriptional modifications. 2΄-O-methylation (2΄-OMe) of the ribose moiety is one of the most abundant post-transcriptional modifications of RNA, although its systematic analysis is difficult due to the lack of reliable high-throughput mapping methods. We describe here a novel high-throughput approach, named 2OMe-seq, that enables fast and accurate mapping at single-base resolution, and relative quantitation, of 2΄-OMe modified residues. We compare our method to other state-of-art approaches, and show that it achieves higher sensitivity and specificity. By applying 2OMe-seq to HeLa cells, we show that it is able to recover the majority of the annotated 2΄-OMe sites on ribosomal RNA. By performing knockdown of the Fibrillarin methyltransferase in mouse embryonic stem cells (ESCs) we show the ability of 2OMe-seq to capture 2΄-O-Methylation level variations. Moreover, using 2OMe-seq data we here report the discovery of 12 previously unannotated 2΄-OMe sites across 18S and 28S rRNAs, 11 of which are conserved in both human and mouse cells, and assigned the respective snoRNAs for all sites. Our approach expands the repertoire of methods for transcriptome-wide mapping of RNA post-transcriptional modifications, and promises to provide novel insights into the role of this modification.

Detection of a Subset of Posttranscriptional Transfer RNA Modifications in Vivo with a Restriction Fragment Length Polymorphism-Based Method

Transfer RNAs (tRNAs) are among the most heavily modified RNA species. Posttranscriptional tRNA modifications (ptRMs) play fundamental roles in modulating tRNA structure and function and are being increasingly linked to human physiology and disease. Detection of ptRMs is often challenging, expensive, and laborious. Restriction fragment length polymorphism (RFLP) analyses study the patterns of DNA cleavage after restriction enzyme treatment and have been used for the qualitative detection of modified bases on mRNAs. It is known that some ptRMs induce specific and reproducible base "mutations" when tRNAs are reverse transcribed. For example, inosine, which derives from the deamination of adenosine, is detected as a guanosine when an inosine-containing tRNA is reverse transcribed, amplified via polymerase chain reaction (PCR), and sequenced. ptRM-dependent base changes on reverse transcription PCR amplicons generated as a consequence of the reverse transcription reaction might create or abolish endonuclease restriction sites. The suitability of RFLP for the detection and/or quantification of ptRMs has not been studied thus far. Here we show that different ptRMs can be detected at specific sites of different tRNA types by RFLP. For the examples studied, we show that this approach can reliably estimate the modification status of the sample, a feature that can be useful in the study of the regulatory role of tRNA modifications in gene expression.

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.

Properties and reactivity of nucleic acids relevant to epigenomics, transcriptomics, and therapeutics

Developments in epigenomics, toxicology, and therapeutic nucleic acids all rely on a precise understanding of nucleic acid properties and chemical reactivity. In this review we discuss the properties and chemical reactivity of each nucleobase and attempt to provide some general principles for nucleic acid targeting or engineering. For adenine-thymine and guanine-cytosine base pairs, we review recent quantum chemical estimates of their Watson-Crick interaction energy, π-π stacking energies, as well as the nuclear quantum effects on tautomerism. Reactions that target nucleobases have been crucial in the development of new sequencing technologies and we believe further developments in nucleic acid chemistry will be required to deconstruct the enormously complex transcriptome.

A to I editing in disease is not fake news

Adenosine deaminases acting on RNA (ADARs) are zinc-containing enzymes that deaminate adenosine bases to inosines within dsRNA regions in transcripts. In short, structured dsRNA hairpins individual adenosine bases may be targeted specifically and edited with up to one hundred percent efficiency, leading to the production of alternative protein variants. However, the majority of editing events occur within longer stretches of dsRNA formed by pairing of repetitive sequences. Here, many different adenosine bases are potential targets but editing efficiency is usually much lower.

Recent work shows that ADAR-mediated RNA editing is also required to prevent aberrant activation of antiviral innate immune sensors that detect viral dsRNA in the cytoplasm. Missense mutations in the ADAR1 RNA editing enzyme cause a fatal auto-inflammatory disease, Aicardi–Goutières syndrome (AGS) in affected children. In addition RNA editing by ADARs has been observed to increase in many cancers and also can contribute to vascular disease. Thus the role of RNA editing in the progression of various diseases can no longer be ignored.

The ability of ADARs to alter the sequence of RNAs has also been used to artificially target model RNAs in vitro and in cells for RNA editing. Potentially this approach may be used to repair genetic defects and to alter genetic information at the RNA level.

In this review we focus on the role of ADARs in disease development and progression and on their potential use to artificially modify RNAs in a targeted manner.