Next Article in Journal
High KRT8 Expression Independently Predicts Poor Prognosis for Lung Adenocarcinoma Patients
Next Article in Special Issue
m6A mRNA Destiny: Chained to the rhYTHm by the YTH-Containing Proteins
Previous Article in Journal
Using a Chemical Genetic Screen to Enhance Our Understanding of the Antimicrobial Properties of Gallium against Escherichia coli
Previous Article in Special Issue
m6A-Dependent RNA Dynamics in T Cell Differentiation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 10
Issue 1
10.3390/genes10010035
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

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

by
Yuri Motorin
1,* and
Mark Helm
2,*
1
Laboratoire IMoPA, UMR7365 National Centre for Scientific Research (CNRS)-Lorraine University, Biopôle, 9 Avenue de la Forêt de Haye, 54505 Vandoeuvre-les-Nancy, France
2
Institute of Pharmacy and Biochemistry, Johannes Gutenberg-University Mainz, Staudingerweg 5, 55128 Mainz, Germany
*
Authors to whom correspondence should be addressed.
Genes 2019, 10(1), 35; https://doi.org/10.3390/genes10010035
Submission received: 16 November 2018 / Revised: 3 January 2019 / Accepted: 4 January 2019 / Published: 9 January 2019
(This article belongs to the Special Issue RNA Modifications)
Browse Figures
Versions Notes

Abstract

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

1. Introduction

Precise and reliable detection of the numerous modified nucleotides in RNA is still a highly challenging task in the epitranscriptomics field. In addition to traditional detection and quantification strategies, including various types of chromatography such as TLC, HPLC, and/or LC-MS(MS), approaches relying on deep sequencing have emerged in force, and numerous applications to large transcriptomes were reported, starting in 2012 [1,2,3,4,5], with a particularly dynamic development during the past two–three years. These methods provided an unprecedented amount of information and depth, in certain cases even allowing the analysis of low abundant RNA species. In this brief overview, we consider techniques that utilize deep sequencing for modification calling, including both well-established and innovative methods. The current portfolio of concepts for the selective detection of RNA modifications ranges from antibody-(Ab)-based anti-methylated/modified nucleoside immunoprecipitation (MeRIP) and variants of individual-nucleotide-resolution crosslinking and immunoprecipitation (iCLIP) method (also called miCLIP) [6,7], over the use of particular enzymes conferring specificity, to applications of known and newly developed organic chemistry of nucleic acids.

2. Antibody-Based Technologies (MeRIP-Seq, i/miCLIP)

The first application of m6A-Ab-based enrichment dates back to the 1980s [8], but only its coupling to deep sequencing brought on the revolution in the field that we are seeing today. The elucidation of human and mouse m6A-methylomes by MeRIP-Seq immunocapturing and sequencing [2,9], was, and still is, being followed by a vast number of studies using this technology, or a derivative thereof. Indeed, despite many openly discussed shortcomings of the various available anti-m6A antibodies, mapping of this modification, especially in mRNA, is arguably the most frequent and visible application in the field. The initial embodiments of MeRIP-Seq had only low resolution and provided broad peaks of ~100–200 nt in which the actual m6A modification site was supposed to be contained. Further improvements of this technique came from the combination of MeRIP with Ab cross-linking, allowing single nucleotide resolution. The latter technique is a specialized derivative of CLIP, termed miCLIP (methylation individual-nucleotide resolution UV cross-linking and immunoprecipitation, [6]).
Thus far, MeRIP-Seq or Ab-related techniques were also applied to the detection of hm5C [10], ac4C [11,12], m5C [13], and of m1A [14,15,16,17], even though m1A data remain highly controversial [18,19]. Despite their massive use in different studies, all Ab-related approaches suffer from poorly characterized, and thus unpredictable specificity of the Ab used for enrichment [20]. These uncertain biases may include sequence context, nucleobase specificity, and a presumed cooperativity in the binding to multiple methylation sites in the same fragment (discussed in Reference [21]), plausibly caused by the presence of two Fab copies in a typical IgG. In addition, the exact stoichiometry of the modification remains largely unknown, since peak enrichment produces signals whose intensities are not likely to be linearly correlated to the original level of modification. More precise quantification was introduced by LAIC-Seq [22], where both m6A-positive and negative mRNA fractions from the same sample were converted to libraries and sequenced in parallel. Careful calibration with spiked-in mRNA standard allows more exact approximation of the m6A stoichiometry in mRNAs. Data obtained in two human cell lines showed that molar modification levels varied considerably from one mRNA site to another, and covered the entire range from almost zero to 1 mol m6A/mol mRNA.
As mentioned above, miCLIP applications to MeRIP-Seq brought considerable improvement, allowing detection via analysis of reverse transcriptase (RT)-stops and/or misincorporations, which occur as a consequence of the cross-link events either at the m6A position or at the +1 site [7,23]. Similarly, m6Am, which is frequently present at m7G-mRNA cap, was detected in the form of a mutational signature [7]. Overall, miCLIP is thought to provide better specificity in mapping, since the simultaneous presences of an enrichment peak and the RT-signature improve confidence in modification calling. Another variant of miCLIP (PA-m6A-Seq) used s4U incorporation for PhotoActivatable-Ribonucleoside-Enhanced Cross-Linking and Immunoprecipitation (PAR-CLIP, [24]).

3. Detection Using Natural or Induced RT-Arrests and Signatures

3.1. Naturally Existing RT-Signatures

Modifications in an RNA template are known to potentially cause reverse transcriptases to either arrest, to incorporate a dNTP cognate to the unmodified template, or to misincorporate non-Watson-Crick compatible dNTPs into the nascent cDNA. The combination of these events constitutes the so-called RT-signature, which varies with the type of modification, and likely with reaction conditions including the type of enzyme used [25]. One of the first attempts to use RT-misincorporation events caused by modifications in the RNA template was termed HAMR for High-throughput Annotation of Modified Ribonucleotides [26,27]. This approach was used to map some potential m3C, m1A, m1I, m22G, and m1G sites in human tRNAs, of which a selection of predicted m3C sites were experimentally validated.
The idea of a specific RT-signature was developed in application to m1A detection in rRNA and tRNAs [25,28]. It was demonstrated that an RT-signature is sequence context dependent, but may be used for mapping of unknown modification sites, and, after thorough calibration, even for their relative quantification.
More recently, the analysis of m1A-generated RT-signatures was extended to analysis of the human transcriptome [17]. Non-enriched mRNA samples show only a very few high-confidence sites with m1A-induced RT-signature, but this level can be increased by prior MeRIP type enrichment with anti-m1A-Ab, indicating that modification levels at the relevant sites were most probably sub-stoichiometric.
The detection of modifications via their RT-signature can only be applied to a limited subset of known RNA modifications, since many do not affect Watson-Crick base pairing during cDNA synthesis and thus are RT-silent. Interestingly, for some presumed RT-silent modifications, particular conditions were identified in which they become visible [29,30]. For example, it was observed early on that ribose 2′-O-methylations in RNA could be detected as RT-stops at reduced dNTP concentrations (1–5 µM final concentration) during primer extension. This method was extensively used in the past for mapping of 2′-O-Me groups in various rRNAs [31]. This methodology was also coupled to NGS for mapping of 2′-O-Me [32], termed 2OMe-seq. Applications of this techniques were only described for validation of methylation sites in rRNA, since RT-primer extension at low [dNTP] is rather prone to generate false-positive hits. In addition, the intensity of RT-stop varies for different sequence contexts, leading to uncertainties in quantification.

3.2. Enhancing RT-Signatures with Engineered Enzymes and Substrates

In addition to manipulating reaction conditions during RT, RT-silent RNA modifications can also be detected using particular, notably engineered RT-enzymes. A first example was an engineered version of a thermostable KlenTaq DNA polymerase, which displayed significant intrinsic RT activity. Mutant KlenTaq DNA polymerases were generated that displayed sensitivity to 2′-O-Me RNA residues even at normal dNTP concentrations (Figure 1a), and these were used for RTL-P-like [33] detection of modifications [34]. Applications of this enzyme to high throughput RNA-Seq have not yet been described, but appear promising in 2′-O-methylation detection. Additionally, starting from a KlenTaq DNA polymerase, protein engineering yielded an enzyme sensitive to m6A (Figure 1b), which exhibits increased misincorporation at respective sites of m6A in the template [35]. This enzyme, too, still awaits application to transcriptome wide analysis.
On the other hand, chemical manipulation at the level of the dNTPs used for primer extension can also exacerbate the RT-signature for specific RNA modifications (Figure 1b). One specific dNTP featuring a selenium substitution of the O4 of thymidine (4Se-dTTP) was found to induce RT-stops and used to generate a specific signature for m6A [36]. Combined with FTO-assisted RNA de-methylation, this approach may potentially become an alternative to Ab-dependent m6A detection strategies with single nucleotide resolution.
Finally, in addition to the manipulation of enzyme and dNTPs, modulation of an RT-signature can also be achieved by manipulation of the RNA template. Enzymatic removal of methyl groups by demethylase enzymes of the AlkB-type was demonstrated for m1A, m3C, and m1G RNA modifications, which are commonly found in tRNAs (Figure 1c). A coupling of this approach to deep sequencing was developed into a method termed ARM-Seq [37], which used pretreatment of RNA with Escherichia coli AlkB to remove RT-blocking modifications. Another method using the same basic idea, in which WT AlkB was combined with an AlkB D135S mutant of improved activity, was called DM-tRNA-Seq [38].
More recently, another AlkB variant (D135S/L118V) was shown to efficiently and selectively de-methylate m22G in tRNAs, thus improving tRNA deep sequencing and allowing more specific detection [39].

4. Exploiting Specific Chemical Reactivity of Modified Nucleobases

4.1. Chemically Induced Alteration of RT-Profiles

The property of naturally occurring RNA modifications to affect the performance of RT enzymes during cDNA synthesis in terms of misincorporation or RT-arrest can be modulated by chemical treatment prior to the RT reaction. Various reagents are known to either selectively react with modified nucleosides, or to selectively affect non-modified nucleotides. An example for the latter is bisulfite sequencing, which was adapted from 5mC detection in DNA to application in RNA by Schäfer and Lyko [40]. The bisulfite reagent is a nucleophile that performs a Michael addition to cytidines, ultimately leading to their deamination. The methyl group in m5C is thought to increase the electron density of the nucleobase, thus leading to repulsion of the nucleophile and thus inertness of m5C to deamination under these specific conditions. Consequently, the readout consists in detection of cytidines in deep seq data sets, and is therefore, in principle, well suited for quantification and single nucleotide resolution mapping. The method is generally considered robust for abundant RNAs such as tRNA and rRNA [40,41,42,43], but results on lower abundant RNAs, including mRNAs and lncRNAs [4,13,44,45], are contested [46].
Further, building on bisulfite chemistry, procedures for mapping oxidation products of m5C were developed, which were, again, first applied on the respective DNA modifications. While hm5C yields signals that are indistinguishable from m5C in bisulfite sequencing, f5C remains invisible. However, the two species can be converted into each other by mild redox agents: f5C can be reduced to hm5C by NaBH4 (RedBS-Seq [47] in DNA, application to RNA: [48]), or be protected from conversion by treatment with O-ethylhydroxylamine [49], application to RNA [48]. Permutations of these treatments coupled to deep seq analysis have been used for mapping in both DNA and RNA. In DNA, hm5C was oxidized to f5C with K2RuO4, followed by an aldol-type addition-elimination-cyclization sequence [50], and the ensuing alteration of base-pairing properties was exploited for sequencing. However, this technique still awaits potential adaptation to RNA-seq.
From PAR-CLIP experiments [51], it was known that the reverse transcription profile of s4U displays a certain percentage of misincorporation resulting in an U- > C transition in the cDNA. SLAM-seq [52] exploits an altered hydrogen binding pattern of s4U after alkylation with iodoacetamide (ICH2-CONH2), which causes s4U to be retrotrancribed as a C, increasing the apparent transition rate to near-quantitative (Figure 2a). While this method was mostly applied to trace the incorporation of s4U nucleoside into nascent RNA in metabolic labeling experiments, it is equally suited to trace posttranscriptional s4U modifications that occur, e.g., at position 8 of bacterial tRNAs.
Inversely, inosine, which is read like a guanosine by known RT-enzymes, does not normally display any tendency to cause RT-arrest. Treatment with acrylonitrile was shown to chemically modify inosines by alkylation, leading to RT-arrest at inosine sites (Figure 2b), and ultimately allowing identification of inosines in a transcriptome-wide search with increased confidence (ICE-Seq [53]).
In contrast, detection of pseudouridine must, so far, entirely rely on chemical treatment for the generation of RT-stops, as the nucleoside itself leaves no exploitable RT-signature. To this end, several labs have performed mapping experiments using soluble carbodiimide (N-Cyclohexyl-N’-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate, usually abbreviated as CMCT) [3,54,55,56,57,58], an agent developed by Bakin and Ofengand [59] for mapping of pseudouridine in ribosomal RNA. A promising variant used a “clickable” derivative, which was conjugated with an affinity tag for the enrichment of derivatized, pseudouridine containing RNA [60,61]. Open questions include the moderate overlap of pseudouridine sites from the different studies [62], potentially a consequence of variable sequencing depth [63].
The reverse transcription signature of an RNA modification is, in principle, understood to also contain abortive RT events resulting in truncated cDNA, for which library preparation protocols must be specifically tailored. The composite nature of the RT-signature of m1A, which was already mentioned [25], can be specifically altered by alkaline treatment, causing a Dimroth rearrangement of m1A, whose product m6A does not display any significant signature. This reaction was used as a means of validation, whereby the RT-signature of m1A was altered by incubation at high temperature under alkaline conditions, leading to the m1A- > m6A conversion [14,17]. The resulting disappearance of m1A signals would be expected to coincide with the appearance of an m6A signal in, e.g., MeRIP-Seq data. Since the alkaline treatment causes significant RNA degradation, and because m6A detection is still inaccurate, applications of the Dimroth rearrangement so far have been limited to validation of predicted m1A sites, rather than for their global profiling.
One of the oldest combination of chemical treatments known to specifically react with a modified nucleotide [30,64] has only recently been combined with an RNA-Seq approach. Reduction of m7G with NaBH4 renders the corresponding site in the RNA chain susceptible to aniline induced cleavage, which was then analyzed as RT-stop in RNA-Seq data in a combination termed TRAC-seq [65].
A chemically related treatment is the reduction of ac4C with NaBH4, leading to saturation of the 5,6 double-bond in acetylated cytidines. This, in turn, produced ~20–30% of misincorporation signals in cDNA, as determined by Sanger sequencing [66]. This method, which still awaits its application to RNA-Seq, might profit from validation of acetylated sites by chemical deacetylation with hydroxylamine, which would be expected to erase signals [11].

4.2. Protection of RNA from Cleavage

Variations of the so-called RiboMethSeq approach, first published by the Nielsen lab and further independently developed by two other groups [67,68,69,70], are based on relative protection of phosphodiester bond in RNA when 5′-neighboring ribose is 2′-O-methylated. This protection can be converted into a signal after random alkaline hydrolysis, conversion of the resulting fragments to the sequencing library, and counting of 5′- and 3′-extremities after alignment to the reference sequence. Modern variants of RiboMethSeq require very low quantities of input RNA (10 ng is generally sufficient), thus opening the possibility for large-scale biomedical applications on clinical samples. Compared to other methods for 2′-O-Me detection (2′OMe-seq [32] and Nm(RibOxi)-Seq [71,72]), RiboMethSeq requires higher sequencing depth (minimum 100 reads/nt position in reference sequence), but provides precise site-by-site quantification of RNA modification, which is difficult to achieve by other approaches. This requirement for significant sequencing depth makes RiboMethSeq less appropriate for low abundant RNAs, but this limitation can be overcome by preliminary enrichment of target RNA or simply by increased sequencing depth for analysis [73].

4.3. Specific Enrichment of RNA Fragments by Selective Ligation

A novel trend in the field is the development of library preparation protocols which selectively produce RNA fragments with modified nucleotides at the 5′- or 3′-end extremity, properties which can then be exploited for enrichment of these fragments in RNA Seq libraries. Two recently developed methods, Nm-Seq [71] and RibOxi-Seq [72], exploit the stability of 2′-O-Me riboses against NaIO4 oxidation, under conditions where the unmodified cis-diol of a normal ribose is rapidly oxidized and cleaved, thus making it inaccessible to 3′-adapter ligation. In this protocol, RNA was randomly cleaved by nuclease to release 5′-phosphates and 3′-riboses, and fragments were subjected to multiple rounds of a sequence consisting of NaIO4 oxidation treatments, subsequent removal of the 3′-end nucleotide, and de-phosphorylation. After several repetitions, an accumulation of Nm at the 3′-terminal position was to be expected. Positive enrichment of 3′-Nm-containing fragments, and thus improvement of Signal/noise ratio, is a very valuable feature for transcriptome-wide RNA analysis. However, the Nm-Seq protocol with repetitive oxidation/cleavage/de-phosphorylation cycles requires considerable amounts of input RNA, due to losses during such treatment.
Another recently established technique, named AlkAnilineSeq [74], relies on an alternative approach for selective fragment enrichment (Figure 2c). The method is based on the properties of RNA abasic sites, which are readily cleaved by aniline via β- and δ-elimination. These two consecutive reactions leave a 5′-phosphate at the N + 1 nucleotide in the sequence, which renders the corresponding fragments uniquely competent for subsequent adapter ligation. Thus, these abasic sites, which were created by alkaline hydrolysis at positions containing m7G and m3C, ultimately led to a positive selection of fragments resulting from modifications, and thus to exquisite sensitivity, because large populations of irrelevant RNA sequences were excluded from the library. This procedure, which was also shown to be partially sensitive to D and ho5C residues, thus required very limited amounts of input RNA and showed exceptional selectivity and specificity because of extremely low background.
The intrinsic reactivity of RNA (or DNA) abasic sites can be explored for detection of RNA modifications. For DNA, formation of abasic site is generally catalyzed by specific DNA-glycosylase, which removes damaged DNA bases, and subsequently the aldehyde group of abasic site can be chemically ligated to, e.g., a reactive amine, which in turn is coupled to biotin [75,76]. A similar approach was also proposed for the detection of oxidized sites in RNA [77,78]. Up to now, aldehyde reactivity at the RNA abasic site or decomposed RNA modified nucleotide was not yet explored for coupling to deep sequencing.

5. Conclusions

At first glance, this overview might convey the impression that RNA-Seq based detection methods can be successfully applied for mapping of almost any desired RNA modification throughout the transcriptome. However, while ever-new concepts for modification mapping get published, the community looks to the biology emerging, based on the methods published since 2012. Now, about six years after the furious resurrection of research on m6A in mRNA, that field is in full swing, with high impact papers appearing on a weekly basis, despite the fact that comprehensive mapping of m6A at single nucleotide resolution in a quantifiable manner is still a largely unresolved problem.
Detection of 2′-O-methylations in abundant RNAs like rRNA and tRNA is quite robust and is now applied even to quantitative analysis of RNA modification modulation under different conditions [79,80,81], but high-confidence transcriptome-wide 2′-O-Me screening has not yet been reported, even if some unpublished observations were made available to epitranscriptomic community (https://doi.org/10.1101/271916).
In contrast, although the performance, at least on paper, of published methods for the mapping of m5C and pseudouridine in mRNA are much better, papers on this topic issue from only a handful of research groups, and especially the numbers of reported modification sites have been met with skepticism [46,62,63]. Similarly, the occurrence of m1A in mRNA, as such, is agreed upon, but the opinions on numbers differ by two orders of magnitude [18,19].
Clearly, such controversies present a formidable obstacle for new groups to engage in targeted research on the effect of mRNA modifications. Research on the impact of defined modification sites is likely to focus on those sites identified by several labs in agreement. This strongly suggests that conservative “candidate-calling” might benefit the field. This means avoidance of false positives at the risk of overlooking some residues, i.e., producing false negatives out of caution.

Funding

Y.M. was supported by ANR-DFG grant HTRNA-Mod ANR-13-ISV8-0001 and Research funds from Lorraine Région (France). M.H. was supported by the DFG (HE3397/8-1, HE3397/13-2, SPP1784).

Acknowledgments

This work was performed in the COST EPITRAN (CA16120) framework.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Squires, J.E.; Patel, H.R.; Nousch, M.; Sibbritt, T.; Humphreys, D.T.; Parker, B.J.; Suter, C.M.; Preiss, T. Widespread occurrence of 5-methylcytosine in human coding and non-coding RNA. Nucleic Acids Res. 2012, 40, 5023–5033. [Google Scholar] [CrossRef] [Green Version]
  2. Dominissini, D.; Moshitch-Moshkovitz, S.; Schwartz, S.; Salmon-Divon, M.; Ungar, L.; Osenberg, S.; Cesarkas, K.; Jacob-Hirsch, J.; Amariglio, N.; Kupiec, M.; et al. Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq. Nature 2012, 485, 201–206. [Google Scholar] [CrossRef] [PubMed]
  3. Carlile, T.M.; Rojas-Duran, M.F.; Zinshteyn, B.; Shin, H.; Bartoli, K.M.; Gilbert, W.V. Pseudouridine profiling reveals regulated mRNA pseudouridylation in yeast and human cells. Nature 2014, 515, 143–146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Edelheit, S.; Schwartz, S.; Mumbach, M.R.; Wurtzel, O.; Sorek, R. Transcriptome-wide mapping of 5-methylcytidine RNA modifications in bacteria, archaea, and yeast reveals m5C within archaeal mRNAs. PLoS Genet. 2013, 9, e1003602. [Google Scholar] [CrossRef] [PubMed]
  5. Meyer, K.D.; Saletore, Y.; Zumbo, P.; Elemento, O.; Mason, C.E.; Jaffrey, S.R. Comprehensive analysis of mRNA methylation reveals enrichment in 3′ UTRs and near stop codons. Cell 2012, 149, 1635–1646. [Google Scholar] [CrossRef] [PubMed]
  6. Hussain, S.; Sajini, A.A.; Blanco, S.; Dietmann, S.; Lombard, P.; Sugimoto, Y.; Paramor, M.; Gleeson, J.G.; Odom, D.T.; Ule, J.; et al. NSun2-mediated cytosine-5 methylation of vault noncoding RNA determines its processing into regulatory small RNAs. Cell Rep. 2013, 4, 255–261. [Google Scholar] [CrossRef] [PubMed]
  7. Linder, B.; Grozhik, A.V.; Olarerin-George, A.O.; Meydan, C.; Mason, C.E.; Jaffrey, S.R. Single-nucleotide-resolution mapping of m6A and m6Am throughout the transcriptome. Nat. Methods 2015, 12, 767–772. [Google Scholar] [CrossRef] [PubMed]
  8. Horowitz, S.; Horowitz, A.; Nilsen, T.W.; Munns, T.W.; Rottman, F.M. Mapping of N6-methyladenosine residues in bovine prolactin mRNA. Proc. Natl. Acad. Sci. USA 1984, 81, 5667–5671. [Google Scholar] [CrossRef]
  9. Dominissini, D.; Moshitch-Moshkovitz, S.; Salmon-Divon, M.; Amariglio, N.; Rechavi, G. Transcriptome-wide mapping of N6-methyladenosine by m6A-seq based on immunocapturing and massively parallel sequencing. Nat. Protoc. 2013, 8, 176–189. [Google Scholar] [CrossRef]
  10. Delatte, B.; Wang, F.; Ngoc, L.V.; Collignon, E.; Bonvin, E.; Deplus, R.; Calonne, E.; Hassabi, B.; Putmans, P.; Awe, S.; et al. RNA biochemistry. Transcriptome-wide distribution and function of RNA hydroxymethylcytosine. Science 2016, 351, 282–285. [Google Scholar] [CrossRef]
  11. Sinclair, W.R.; Arango, D.; Shrimp, J.H.; Zengeya, T.T.; Thomas, J.M.; Montgomery, D.C.; Fox, S.D.; Andresson, T.; Oberdoerffer, S.; Meier, J.L. Profiling cytidine acetylation with specific affinity and reactivity. ACS Chem. Biol. 2017, 12, 2922–2926. [Google Scholar] [CrossRef] [PubMed]
  12. Arango, D.; Sturgill, D.; Alhusaini, N.; Dillman, A.A.; Sweet, T.J.; Hanson, G.; Hosogane, M.; Sinclair, W.R.; Nanan, K.K.; Mandler, M.D.; et al. Acetylation of cytidine in mRNA promotes translation efficiency. Cell 2018, 175, 1872–1886.e24. [Google Scholar] [CrossRef] [PubMed]
  13. Amort, T.; Rieder, D.; Wille, A.; Khokhlova-Cubberley, D.; Riml, C.; Trixl, L.; Jia, X.-Y.; Micura, R.; Lusser, A. Distinct 5-methylcytosine profiles in poly(A) RNA from mouse embryonic stem cells and brain. Genome Biol. 2017, 18, 1. [Google Scholar] [CrossRef]
  14. Dominissini, D.; Nachtergaele, S.; Moshitch-Moshkovitz, S.; Peer, E.; Kol, N.; Ben-Haim, M.S.; Dai, Q.; Di Segni, A.; Salmon-Divon, M.; Clark, W.C.; et al. The dynamic N(1)-methyladenosine methylome in eukaryotic messenger RNA. Nature 2016, 530, 441–446. [Google Scholar] [CrossRef] [PubMed]
  15. Li, X.; Xiong, X.; Zhang, M.; Wang, K.; Chen, Y.; Zhou, J.; Mao, Y.; Lv, J.; Yi, D.; Chen, X.-W.; et al. Base-resolution mapping reveals distinct m1A methylome in nuclear- and mitochondrial-encoded transcripts. Mol. Cell 2017, 68, 993–1005.e9. [Google Scholar] [CrossRef]
  16. Li, X.; Xiong, X.; Wang, K.; Wang, L.; Shu, X.; Ma, S.; Yi, C. Transcriptome-wide mapping reveals reversible and dynamic N(1)-methyladenosine methylome. Nat. Chem. Biol. 2016, 12, 311–316. [Google Scholar] [CrossRef] [PubMed]
  17. Safra, M.; Sas-Chen, A.; Nir, R.; Winkler, R.; Nachshon, A.; Bar-Yaacov, D.; Erlacher, M.; Rossmanith, W.; Stern-Ginossar, N.; Schwartz, S. The m1A landscape on cytosolic and mitochondrial mRNA at single-base resolution. Nature 2017, 551, 251–255. [Google Scholar] [CrossRef]
  18. Xiong, X.; Li, X.; Wang, K.; Yi, C. Perspectives on topology of the human m1A methylome at single nucleotide resolution. RNA 2018, 24, 1437–1442. [Google Scholar] [CrossRef]
  19. Schwartz, S. m1A within cytoplasmic mRNAs at single nucleotide resolution: A reconciled transcriptome-wide map. RNA 2018, 24, 1427–1436. [Google Scholar] [CrossRef]
  20. Mishima, E.; Jinno, D.; Akiyama, Y.; Itoh, K.; Nankumo, S.; Shima, H.; Kikuchi, K.; Takeuchi, Y.; Elkordy, A.; Suzuki, T.; et al. Immuno-Northern Blotting: Detection of RNA modifications by using antibodies against modified nucleosides. PLoS ONE 2015, 10, e0143756. [Google Scholar] [CrossRef]
  21. Grozhik, A.V.; Jaffrey, S.R. Distinguishing RNA modifications from noise in epitranscriptome maps. Nat. Chem. Biol. 2018, 14, 215–225. [Google Scholar] [CrossRef] [PubMed]
  22. Molinie, B.; Wang, J.; Lim, K.S.; Hillebrand, R.; Lu, Z.-X.; Van Wittenberghe, N.; Howard, B.D.; Daneshvar, K.; Mullen, A.C.; Dedon, P.; et al. m6A-LAIC-seq reveals the census and complexity of the m6A epitranscriptome. Nat. Methods 2016, 13, 692–698. [Google Scholar] [CrossRef] [PubMed]
  23. Grozhik, A.V.; Linder, B.; Olarerin-George, A.O.; Jaffrey, S.R. Mapping m6A at individual-nucleotide resolution using crosslinking and immunoprecipitation (miCLIP). Methods Mol. Biol. 2017, 1562, 55–78. [Google Scholar] [PubMed]
  24. Chen, K.; Lu, Z.; Wang, X.; Fu, Y.; Luo, G.-Z.; Liu, N.; Han, D.; Dominissini, D.; Dai, Q.; Pan, T.; et al. High-resolution N(6)-methyladenosine (m6A) map using photo-crosslinking-assisted m6A sequencing. Angew. Chem. Int. Ed. Engl. 2015, 54, 1587–1590. [Google Scholar] [CrossRef] [PubMed]
  25. Hauenschild, R.; Tserovski, L.; Schmid, K.; Thüring, K.; Winz, M.-L.; Sharma, S.; Entian, K.-D.; Wacheul, L.; Lafontaine, D.L.J.; Anderson, J.; et al. The reverse transcription signature of N-1-methyladenosine in RNA-Seq is sequence dependent. Nucleic Acids Res. 2015, 43, 9950–9964. [Google Scholar] [CrossRef] [PubMed]
  26. Ryvkin, P.; Leung, Y.Y.; Silverman, I.M.; Childress, M.; Valladares, O.; Dragomir, I.; Gregory, B.D.; Wang, L.-S. HAMR: High-throughput annotation of modified ribonucleotides. RNA 2013, 19, 1684–1692. [Google Scholar] [CrossRef] [PubMed]
  27. Kuksa, P.P.; Leung, Y.Y.; Vandivier, L.E.; Anderson, Z.; Gregory, B.D.; Wang, L.-S. In silico identification of RNA modifications from high-throughput sequencing data using HAMR. Methods Mol. Biol. 2017, 1562, 211–229. [Google Scholar] [PubMed]
  28. Tserovski, L.; Marchand, V.; Hauenschild, R.; Blanloeil-Oillo, F.; Helm, M.; Motorin, Y. High-throughput sequencing for 1-methyladenosine (m1A) mapping in RNA. Methods 2016, 107, 110–121. [Google Scholar] [CrossRef]
  29. Motorin, Y.; Muller, S.; Behm-Ansmant, I.; Branlant, C. Identification of modified residues in RNAs by reverse transcription-based methods. Meth. Enzymol. 2007, 425, 21–53. [Google Scholar]
  30. Behm-Ansmant, I.; Helm, M.; Motorin, Y. Use of specific chemical reagents for detection of modified nucleotides in RNA. J. Nucleic Acids 2011, 2011, 408053. [Google Scholar] [CrossRef]
  31. Maden, B.E. Mapping 2′-O-methyl groups in ribosomal RNA. Methods 2001, 25, 374–382. [Google Scholar] [CrossRef] [PubMed]
  32. Incarnato, D.; Anselmi, F.; Morandi, E.; Neri, F.; Maldotti, M.; Rapelli, S.; Parlato, C.; Basile, G.; Oliviero, S. High-throughput single-base resolution mapping of RNA 2′-O-methylated residues. Nucleic Acids Res. 2017, 45, 1433–1441. [Google Scholar] [CrossRef] [PubMed]
  33. Dong, Z.-W.; Shao, P.; Diao, L.-T.; Zhou, H.; Yu, C.-H.; Qu, L.-H. RTL-P: A sensitive approach for detecting sites of 2′-O-methylation in RNA molecules. Nucleic Acids Res. 2012, 40, e157. [Google Scholar] [CrossRef] [PubMed]
  34. Aschenbrenner, J.; Marx, A. Direct and site-specific quantification of RNA 2′-O-methylation by PCR with an engineered DNA polymerase. Nucleic Acids Res. 2016, 44, 3495–3502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Aschenbrenner, J.; Werner, S.; Marchand, V.; Adam, M.; Motorin, Y.; Helm, M.; Marx, A. Engineering of a DNA polymerase for direct m6A sequencing. Angew. Chem. Int. Ed. Engl. 2018, 57, 417–421. [Google Scholar] [CrossRef] [PubMed]
  36. Hong, T.; Yuan, Y.; Chen, Z.; Xi, K.; Wang, T.; Xie, Y.; He, Z.; Su, H.; Zhou, Y.; Tan, Z.-J.; et al. Precise antibody-independent m6A identification via 4SedTTP-involved and FTO-assisted strategy at single-nucleotide resolution. J. Am. Chem. Soc. 2018, 140, 5886–5889. [Google Scholar] [CrossRef]
  37. Cozen, A.E.; Quartley, E.; Holmes, A.D.; Hrabeta-Robinson, E.; Phizicky, E.M.; Lowe, T.M. ARM-seq: AlkB-facilitated RNA methylation sequencing reveals a complex landscape of modified tRNA fragments. Nat. Methods 2015, 12, 879–884. [Google Scholar] [CrossRef] [Green Version]
  38. Zheng, G.; Qin, Y.; Clark, W.C.; Dai, Q.; Yi, C.; He, C.; Lambowitz, A.M.; Pan, T. Efficient and quantitative high-throughput tRNA sequencing. Nat. Methods 2015, 12, 835–837. [Google Scholar] [CrossRef] [Green Version]
  39. Dai, Q.; Zheng, G.; Schwartz, M.H.; Clark, W.C.; Pan, T. Selective enzymatic demethylation of N2,N2-dimethylguanosine in RNA and its application in high-throughput tRNA sequencing. Angew. Chem. Int. Ed. Engl. 2017, 56, 5017–5020. [Google Scholar] [CrossRef]
  40. Schaefer, M.; Pollex, T.; Hanna, K.; Lyko, F. RNA cytosine methylation analysis by bisulfite sequencing. Nucleic Acids Res. 2009, 37, e12. [Google Scholar] [CrossRef]
  41. Müller, M.; Hartmann, M.; Schuster, I.; Bender, S.; Thüring, K.L.; Helm, M.; Katze, J.R.; Nellen, W.; Lyko, F.; Ehrenhofer-Murray, A.E. Dynamic modulation of Dnmt2-dependent tRNA methylation by the micronutrient queuine. Nucleic Acids Res. 2015, 43, 10952–10962. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Tuorto, F.; Liebers, R.; Musch, T.; Schaefer, M.; Hofmann, S.; Kellner, S.; Frye, M.; Helm, M.; Stoecklin, G.; Lyko, F. RNA cytosine methylation by Dnmt2 and NSun2 promotes tRNA stability and protein synthesis. Nat. Struct. Mol. Biol. 2012, 19, 900–905. [Google Scholar] [CrossRef] [PubMed]
  43. Bourgeois, G.; Ney, M.; Gaspar, I.; Aigueperse, C.; Schaefer, M.; Kellner, S.; Helm, M.; Motorin, Y. Eukaryotic rRNA modification by yeast 5-methylcytosine-methyltransferases and human proliferation-associated antigen p120. PLoS ONE 2015, 10, e0133321. [Google Scholar] [CrossRef] [PubMed]
  44. David, R.; Burgess, A.; Parker, B.; Li, J.; Pulsford, K.; Sibbritt, T.; Preiss, T.; Searle, I.R. Transcriptome-wide mapping of RNA 5-Methylcytosine in Arabidopsis mRNAs and noncoding RNAs. Plant Cell 2017, 29, 445–460. [Google Scholar] [CrossRef] [PubMed]
  45. Wei, Z.; Panneerdoss, S.; Timilsina, S.; Zhu, J.; Mohammad, T.A.; Lu, Z.-L.; de Magalhães, J.P.; Chen, Y.; Rong, R.; Huang, Y.; et al. Topological characterization of human and mouse m5C epitranscriptome revealed by bisulfite sequencing. Int. J. Genom. 2018, 2018, 1351964. [Google Scholar] [CrossRef] [PubMed]
  46. Legrand, C.; Tuorto, F.; Hartmann, M.; Liebers, R.; Jacob, D.; Helm, M.; Lyko, F. Statistically robust methylation calling for whole-transcriptome bisulfite sequencing reveals distinct methylation patterns for mouse RNAs. Genome Res. 2017, 27, 1589–1596. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Booth, M.J.; Marsico, G.; Bachman, M.; Beraldi, D.; Balasubramanian, S. Quantitative sequencing of 5-formylcytosine in DNA at single-base resolution. Nat. Chem. 2014, 6, 435–440. [Google Scholar] [CrossRef] [Green Version]
  48. Van Haute, L.; Dietmann, S.; Kremer, L.; Hussain, S.; Pearce, S.F.; Powell, C.A.; Rorbach, J.; Lantaff, R.; Blanco, S.; Sauer, S.; et al. Deficient methylation and formylation of mt-tRNAMet wobble cytosine in a patient carrying mutations in NSUN3. Nat. Commun. 2016, 7, 12039. [Google Scholar] [CrossRef]
  49. Song, C.-X.; Szulwach, K.E.; Dai, Q.; Fu, Y.; Mao, S.-Q.; Lin, L.; Street, C.; Li, Y.; Poidevin, M.; Wu, H.; et al. Genome-wide profiling of 5-formylcytosine reveals its roles in epigenetic priming. Cell 2013, 153, 678–691. [Google Scholar] [CrossRef]
  50. Zeng, H.; He, B.; Xia, B.; Bai, D.; Lu, X.; Cai, J.; Chen, L.; Zhou, A.; Zhu, C.; Meng, H.; et al. Bisulfite-Free, Nanoscale analysis of 5-hydroxymethylcytosine at single base resolution. J. Am. Chem. Soc. 2018, 140, 13190–13194. [Google Scholar] [CrossRef]
  51. Hafner, M.; Landthaler, M.; Burger, L.; Khorshid, M.; Hausser, J.; Berninger, P.; Rothballer, A.; Ascano, M.; Jungkamp, A.-C.; Munschauer, M.; et al. Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP. Cell 2010, 141, 129–141. [Google Scholar] [CrossRef] [PubMed]
  52. Herzog, V.A.; Reichholf, B.; Neumann, T.; Rescheneder, P.; Bhat, P.; Burkard, T.R.; Wlotzka, W.; von Haeseler, A.; Zuber, J.; Ameres, S.L. Thiol-linked alkylation of RNA to assess expression dynamics. Nat. Methods 2017, 14, 1198–1204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Suzuki, T.; Ueda, H.; Okada, S.; Sakurai, M. Transcriptome-wide identification of adenosine-to-inosine editing using the ICE-seq method. Nat. Protoc. 2015, 10, 715–732. [Google Scholar] [CrossRef] [PubMed]
  54. Carlile, T.M.; Rojas-Duran, M.F.; Gilbert, W.V. Pseudo-Seq: Genome-wide detection of pseudouridine modifications in RNA. Meth. Enzymol. 2015, 560, 219–245. [Google Scholar] [PubMed]
  55. Carlile, T.M.; Rojas-Duran, M.F.; Gilbert, W.V. Transcriptome-wide identification of pseudouridine modifications using pseudo-seq. Curr. Protoc. Mol. Biol. 2015, 112, 1–24. [Google Scholar]
  56. Schwartz, S.; Bernstein, D.A.; Mumbach, M.R.; Jovanovic, M.; Herbst, R.H.; León-Ricardo, B.X.; Engreitz, J.M.; Guttman, M.; Satija, R.; Lander, E.S.; et al. Transcriptome-wide mapping reveals widespread dynamic-regulated pseudouridylation of ncRNA and mRNA. Cell 2014, 159, 148–162. [Google Scholar] [CrossRef] [PubMed]
  57. Lovejoy, A.F.; Riordan, D.P.; Brown, P.O. Transcriptome-wide mapping of pseudouridines: Pseudouridine synthases modify specific mRNAs in S. cerevisiae. PLoS ONE 2014, 9, e110799. [Google Scholar] [CrossRef] [PubMed]
  58. Nakamoto, M.A.; Lovejoy, A.F.; Cygan, A.M.; Boothroyd, J.C. mRNA pseudouridylation affects RNA metabolism in the parasite Toxoplasma gondii. RNA 2017, 23, 1834–1849. [Google Scholar] [CrossRef]
  59. Bakin, A.; Ofengand, J. Four newly located pseudouridylate residues in Escherichia coli 23S ribosomal RNA are all at the peptidyltransferase center: Analysis by the application of a new sequencing technique. Biochemistry 1993, 32, 9754–9762. [Google Scholar] [CrossRef]
  60. Li, X.; Zhu, P.; Ma, S.; Song, J.; Bai, J.; Sun, F.; Yi, C. Chemical pulldown reveals dynamic pseudouridylation of the mammalian transcriptome. Nat. Chem. Biol. 2015, 11, 592–597. [Google Scholar] [CrossRef]
  61. Li, X.; Ma, S.; Yi, C. Pseudouridine chemical labeling and profiling. Meth. Enzymol. 2015, 560, 247–272. [Google Scholar] [PubMed]
  62. Zaringhalam, M.; Papavasiliou, F.N. Pseudouridylation meets next-generation sequencing. Methods 2016, 107, 63–72. [Google Scholar] [CrossRef] [PubMed]
  63. Safra, M.; Nir, R.; Farouq, D.; Vainberg Slutskin, I.; Schwartz, S. TRUB1 is the predominant pseudouridine synthase acting on mammalian mRNA via a predictable and conserved code. Genome Res. 2017, 27, 393–406. [Google Scholar] [CrossRef] [Green Version]
  64. Wintermeyer, W.; Zachau, H.G. Tertiary structure interactions of 7-methylguanosine in yeast tRNAPhe as studied by borohydride reduction. FEBS Lett. 1975, 58, 306–309. [Google Scholar] [CrossRef]
  65. Lin, S.; Liu, Q.; Lelyveld, V.S.; Choe, J.; Szostak, J.W.; Gregory, R.I. Mettl1/Wdr4-mediated m7G tRNA methylome is required for normal mRNA translation and embryonic stem cell self-renewal and differentiation. Mol. Cell 2018, 71, 244–255.e5. [Google Scholar] [CrossRef] [PubMed]
  66. Thomas, J.M.; Briney, C.A.; Nance, K.D.; Lopez, J.E.; Thorpe, A.L.; Fox, S.D.; Bortolin-Cavaille, M.-L.; Sas-Chen, A.; Arango, D.; Oberdoerffer, S.; et al. A chemical signature for cytidine acetylation in RNA. J. Am. Chem. Soc. 2018, 140, 12667–12670. [Google Scholar] [CrossRef] [PubMed]
  67. Birkedal, U.; Christensen-Dalsgaard, M.; Krogh, N.; Sabarinathan, R.; Gorodkin, J.; Nielsen, H. Profiling of ribose methylations in RNA by high-throughput sequencing. Angew. Chem. Int. Ed. Engl. 2015, 54, 451–455. [Google Scholar] [CrossRef]
  68. Marchand, V.; Blanloeil-Oillo, F.; Helm, M.; Motorin, Y. Illumina-based RiboMethSeq approach for mapping of 2′-O-Me residues in RNA. Nucleic Acids Res. 2016, 44, e135. [Google Scholar] [CrossRef]
  69. Marchand, V.; Pichot, F.; Thüring, K.; Ayadi, L.; Freund, I.; Dalpke, A.; Helm, M.; Motorin, Y. Next-generation sequencing-based RiboMethSeq protocol for analysis of tRNA 2′-O-methylation. Biomolecules 2017, 7, 13. [Google Scholar] [CrossRef]
  70. Gumienny, R.; Jedlinski, D.J.; Schmidt, A.; Gypas, F.; Martin, G.; Vina-Vilaseca, A.; Zavolan, M. High-throughput identification of C/D box snoRNA targets with CLIP and RiboMeth-seq. Nucleic Acids Res. 2017, 45, 2341–2353. [Google Scholar] [CrossRef]
  71. Dai, Q.; Moshitch-Moshkovitz, S.; Han, D.; Kol, N.; Amariglio, N.; Rechavi, G.; Dominissini, D.; He, C. Nm-seq maps 2′-O-methylation sites in human mRNA with base precision. Nat. Methods 2017, 14, 695–698. [Google Scholar] [CrossRef] [PubMed]
  72. Zhu, Y.; Pirnie, S.P.; Carmichael, G.G. High-throughput and site-specific identification of 2′-O-methylation sites using ribose oxidation sequencing (RibOxi-seq). RNA 2017, 23, 1303–1314. [Google Scholar] [CrossRef] [PubMed]
  73. Krogh, N.; Kongsbak-Wismann, M.; Geisler, C.; Nielsen, H. Substoichiometric ribose methylations in spliceosomal snRNAs. Org. Biomol. Chem. 2017, 15, 8872–8876. [Google Scholar] [CrossRef] [PubMed]
  74. Marchand, V.; Ayadi, L.; Ernst, F.G.M.; Hertler, J.; Bourguignon-Igel, V.; Galvanin, A.; Kotter, A.; Helm, M.; Lafontaine, D.L.J.; Motorin, Y. AlkAniline-Seq: Profiling of m7G and m3C RNA modifications at single nucleotide resolution. Angew. Chem. Int. Ed. Engl. 2018, 57, 16785–16790. [Google Scholar] [CrossRef]
  75. Kubo, K.; Ide, H.; Wallace, S.S.; Kow, Y.W. A novel, sensitive, and specific assay for abasic sites, the most commonly produced DNA lesion. Biochemistry 1992, 31, 3703–3708. [Google Scholar] [CrossRef] [PubMed]
  76. Ide, H.; Akamatsu, K.; Kimura, Y.; Michiue, K.; Makino, K.; Asaeda, A.; Takamori, Y.; Kubo, K. Synthesis and damage specificity of a novel probe for the detection of abasic sites in DNA. Biochemistry 1993, 32, 8276–8283. [Google Scholar] [CrossRef] [PubMed]
  77. Tanaka, M.; Han, S.; Song, H.; Küpfer, P.A.; Leumann, C.J.; Sonntag, W.E. An assay for RNA oxidation induced abasic sites using the Aldehyde Reactive Probe. Free Radic. Res. 2011, 45, 237–247. [Google Scholar] [CrossRef]
  78. Tanaka, M.; Han, S.; Küpfer, P.A.; Leumann, C.J.; Sonntag, W.E. Quantification of oxidized levels of specific RNA species using an aldehyde reactive probe. Anal. Biochem. 2011, 417, 142–148. [Google Scholar] [CrossRef] [Green Version]
  79. Krogh, N.; Jansson, M.D.; Häfner, S.J.; Tehler, D.; Birkedal, U.; Christensen-Dalsgaard, M.; Lund, A.H.; Nielsen, H. Profiling of 2′-O-Me in human rRNA reveals a subset of fractionally modified positions and provides evidence for ribosome heterogeneity. Nucleic Acids Res. 2016, 44, 7884–7895. [Google Scholar] [CrossRef] [Green Version]
  80. Sharma, S.; Marchand, V.; Motorin, Y.; Lafontaine, D.L.J. Identification of sites of 2′-O-methylation vulnerability in human ribosomal RNAs by systematic mapping. Sci. Rep. 2017, 7, 11490. [Google Scholar] [CrossRef]
  81. Erales, J.; Marchand, V.; Panthu, B.; Gillot, S.; Belin, S.; Ghayad, S.E.; Garcia, M.; Laforêts, F.; Marcel, V.; Baudin-Baillieu, A.; et al. Evidence for rRNA 2′-O-methylation plasticity: Control of intrinsic translational capabilities of human ribosomes. Proc. Natl. Acad. Sci. USA 2017, 114, 12934–12939. [Google Scholar] [CrossRef] [PubMed]
Genes 10 00035 g001 550
Figure 1. Detection of 2′-O-methylation, m6A and AlkB-sensitive modifications in RNA (a) Engineered polymerase arrests reverse transcription at 2′-O-methylated residue at normal [dNTP]; (b) Ab-independent detection of m6A residue using an engineered polymerase or a 4-Se derivative of dTTP (4Se-dTTP). The presence of an m6A residue in the RNA template is reflected by a misincorporation signature in the cDNA; (c) Improved confidence in the detection of AlkB-sensitive RNA modifications (m3C, m1A, m1G, and m22G) comes from comparison of RT-signatures (including misincorporation and RT-arrest) before and after AlkB treatment.
Figure 1. Detection of 2′-O-methylation, m6A and AlkB-sensitive modifications in RNA (a) Engineered polymerase arrests reverse transcription at 2′-O-methylated residue at normal [dNTP]; (b) Ab-independent detection of m6A residue using an engineered polymerase or a 4-Se derivative of dTTP (4Se-dTTP). The presence of an m6A residue in the RNA template is reflected by a misincorporation signature in the cDNA; (c) Improved confidence in the detection of AlkB-sensitive RNA modifications (m3C, m1A, m1G, and m22G) comes from comparison of RT-signatures (including misincorporation and RT-arrest) before and after AlkB treatment.
Genes 10 00035 g001
Genes 10 00035 g002 550
Figure 2. Use of specific chemical reagents for detection of s4U, inosine (I) and m7G/m3C: (a) 4-thiouridine (s4U) readily reacts with iodoacetamide, with subsequent change of base-paring properties. After iodoacetamide treatment, derivatized s4U is read as C during RT step and this U- > C transition is detected in the sequencing data; (b) Detection of inosine (I) in RNA by acrylonitrile derivatization (ICE-Seq). Unmodified inosine gives an A- > G misincorporation signature, which can be distinguished from SNPs by additional treatment with acrylonitrile. The resulting derivatized inosine creates abortive products during RT-primer extension; (c) AlkAniline-Seq protocol for detection of m7G and m3C residues in RNA. Mild alkaline hydrolysis creates an RNA abasic site at the sensitive residue and subsequent aniline treatment releases 5′-phosphate at the N + 1 nucleotide. These 5′-phosphates are used for selective adapter ligation during library preparation step. AlkAniline-Seq signal is detected as a substantial proportion of sequencing reads starting at the same position.
Figure 2. Use of specific chemical reagents for detection of s4U, inosine (I) and m7G/m3C: (a) 4-thiouridine (s4U) readily reacts with iodoacetamide, with subsequent change of base-paring properties. After iodoacetamide treatment, derivatized s4U is read as C during RT step and this U- > C transition is detected in the sequencing data; (b) Detection of inosine (I) in RNA by acrylonitrile derivatization (ICE-Seq). Unmodified inosine gives an A- > G misincorporation signature, which can be distinguished from SNPs by additional treatment with acrylonitrile. The resulting derivatized inosine creates abortive products during RT-primer extension; (c) AlkAniline-Seq protocol for detection of m7G and m3C residues in RNA. Mild alkaline hydrolysis creates an RNA abasic site at the sensitive residue and subsequent aniline treatment releases 5′-phosphate at the N + 1 nucleotide. These 5′-phosphates are used for selective adapter ligation during library preparation step. AlkAniline-Seq signal is detected as a substantial proportion of sequencing reads starting at the same position.
Genes 10 00035 g002

Share and Cite

MDPI and ACS Style

Motorin, Y.; Helm, M. Methods for RNA Modification Mapping Using Deep Sequencing: Established and New Emerging Technologies. Genes 2019, 10, 35. https://doi.org/10.3390/genes10010035

AMA Style

Motorin Y, Helm M. Methods for RNA Modification Mapping Using Deep Sequencing: Established and New Emerging Technologies. Genes. 2019; 10(1):35. https://doi.org/10.3390/genes10010035

Chicago/Turabian Style

Motorin, Yuri, and Mark Helm. 2019. "Methods for RNA Modification Mapping Using Deep Sequencing: Established and New Emerging Technologies" Genes 10, no. 1: 35. https://doi.org/10.3390/genes10010035

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop