NOseq: amplicon sequencing evaluation method for RNA m6A sites after chemical deamination

General Principles and Limitations for Detection of RNA Modifications by Sequencing

ConspectusAmong the many analytical methods applied to RNA modifications, a particularly pronounced surge has occurred in the past decade in the field of modification mapping. The occurrence of modifications such as m6A in mRNA, albeit known since the 1980s, became amenable to transcriptome-wide analyses through the advent of next-generation sequencing techniques in a rather sudden manner. The term "mapping" here refers to detection of RNA modifications in a sequence context, which has a dramatic impact on the interpretation of biological functions. As a consequence, an impressive number of mapping techniques were published, most in the perspective of what now has become known as "epitranscriptomics". While more and more different modifications were reported to occur in mRNA, conflicting reports and controversial results pointed to a number of technical and theoretical problems rooted in analytics, statistics, and reagents. Rather than finding the proverbial needle in a haystack, the tasks were to determine how many needles of what color in what size of a haystack one was looking at.As the authors of this Account, we think it important to outline the limitations of different mapping methods since many life scientists freshly entering the field confuse the accuracy and precision of modification mapping with that of normal sequencing, which already features numerous caveats by itself. Indeed, we propose here to qualify a specific mapping method by the size of the transcriptome that can be meaningfully analyzed with it.We here focus on high throughput sequencing by Illumina technology, referred to as RNA-Seq. We noted with interest the development of methods for modification detection by other high throughput sequencing platforms that act directly on RNA, e.g., PacBio SMRT and nanopore sequencing, but those are not considered here.In contrast to approaches relying on direct RNA sequencing, current Illumina RNA-Seq protocols require prior conversion of RNA into DNA. This conversion relies on reverse transcription (RT) to create cDNA; thereafter, the cDNA undergoes a sequencing-by-synthesis type of analysis. Thus, a particular behavior of RNA modified nucleotides during the RT-step is a prerequisite for their detection (and quantification) by deep sequencing, and RT properties have great influence on the detection efficiency and reliability. Moreover, the RT-step requires annealing of a synthetic primer, a prerequisite with a crucial impact on library preparation. Thus, all RNA-Seq protocols must feature steps for the introduction of primers, primer landing sites, or adapters on both the RNA 3'- and 5'-ends.

Exploring N6-methyladenosine (m6A) modification in tree species: opportunities and challenges

N 6-methyladenosine (m6A) in eukaryotes is the most common and widespread internal modification in mRNA. The modification regulates mRNA stability, translation efficiency, and splicing, thereby fine-tuning gene regulation. In plants, m6A is dynamic and critical for various growth stages, embryonic development, morphogenesis, flowering, stress response, crop yield, and biomass. Although recent high-throughput sequencing approaches have enabled the rapid identification of m6A modification sites, the site-specific mechanism of this modification remains unclear in trees. In this review, we discuss the functional significance of m6A in trees under different stress conditions and discuss recent advancements in the quantification of m6A. Quantitative and functional insights into the dynamic aspect of m6A modification could assist researchers in engineering tree crops for better productivity and resistance to various stress conditions.

Experience with German Research Consortia in the Field of Chemical Biology of Native Nucleic Acid Modifications

The chemical biology of native nucleic acid modifications has seen an intense upswing, first concerning DNA modifications in the field of epigenetics and then concerning RNA modifications in a field that was correspondingly rebaptized epitranscriptomics by analogy. The German Research Foundation (DFG) has funded several consortia with a scientific focus in these fields, strengthening the traditionally well-developed nucleic acid chemistry community and inciting it to team up with colleagues from the life sciences and data science to tackle interdisciplinary challenges. This Perspective focuses on the genesis, scientific outcome, and downstream impact of the DFG priority program SPP1784 and offers insight into how it fecundated further consortia in the field. Pertinent research was funded from mid-2015 to 2022, including an extension related to the coronavirus pandemic. Despite being a detriment to research activity in general, the pandemic has resulted in tremendously boosted interest in the field of RNA and RNA modifications as a consequence of their widespread and successful use in vaccination campaigns against SARS-CoV-2. Funded principal investigators published over 250 pertinent papers with a very substantial impact on the field. The program also helped to redirect numerous laboratories toward this dynamic field. Finally, SPP1784 spawned initiatives for several funded consortia that continue to drive the fields of nucleic acid modification.

Quantifying m6A and Ψ Modifications in the Transcriptome via Chemical-Assisted Approaches

Since the discovery of the first chemically modified RNA nucleotide in 1951, more than 170 types of chemical modifications have been characterized in RNA so far. Since the discovery of the reversible and dynamic nature of N6-methyladenosine (m6A) in mRNA modification, researchers have identified about ten modifications in eukaryotic mRNA; together with modifications on the noncoding RNAs, the term "epitranscriptome" has been coined to describe the ensemble of various chemical RNA modifications. The past decade has witnessed the discovery of many novel molecular functions of mRNA modifications, demonstrating their crucial roles in gene expression regulation. As the most abundant modifications in mRNA, the study of m6A and Ψ has been facilitated by innovative high-throughput sequencing technologies, which can be based on antibodies, enzymes, or novel chemistry. Among them, chemical-assisted methods utilize selective chemistry that can discriminate modified versus unmodified nucleotides, enabling the transcriptome-wide mapping of m6A and Ψ modifications and functional studies.Our group has developed several sequencing technologies to investigate these epitranscriptomic marks including m6A, Ψ, m1A, and m6Am. Among them, we have recently developed two methods for absolute quantification of m6A and Ψ in the transcriptome based on chemical reactivity to distinguish and measure the two modifications. In GLORI, we used glyoxal and nitrite to mediate efficient deamination of regular adenosine, while m6A remained unaffected, thereby enabling efficient and unbiased detection of single-base resolution and absolute quantification of m6A modification. In CeU-seq and PRAISE, we used different chemistry to achieve selective labeling and detection of transcriptome-wide Ψ. CeU-seq is based on an azido-derivatized carbodiimide compound, while PRAISE utilizes the unique activity of bisulfite to Ψ. PRAISE results in the formation of ring-opening Ψ-bisulfite adduct and quantitatively detects Ψ as 1-2 nt deletion signatures during sequencing. The resulting base-resolution and stoichiometric information expanded our understanding to the profiles of RNA modifications in the transcriptome. In particular, the quantitative information on RNA methylome is critical for characterizing the dynamic and reversible nature of RNA modifications, for instance, under environmental stress or during development. Additionally, base-resolution and stoichiometric information can greatly facilitate the analysis and characterization of functional modification sites that are important for gene expression regulation, especially when one modification type may have multiple or even opposing functions within a specific transcript. Together, the quantitative profiling methods provide the modification stoichiometry information, which is critical to study the regulatory roles of RNA modifications.In this Account, we will focus on the quantitative sequencing technologies of m6A and Ψ developed in our group, review recent advances in chemical-assisted reactions for m6A and Ψ detection, and discuss the challenges and future opportunities of transcriptome-wide mapping technologies for RNA modifications.

Aberrant RNA m6A modification in gastrointestinal malignancies: versatile regulators of cancer hallmarks and novel therapeutic opportunities

Gastrointestinal (GI) cancer is one of the most common malignancies, and a leading cause of cancer-related death worldwide. However, molecular targeted therapies are still lacking, leading to poor treatment efficacies. As an important layer of epigenetic regulation, RNA N6-Methyladenosine (m⁶A) modification is recently linked to various biological hallmarks of cancer by orchestrating RNA metabolism, including RNA splicing, export, translation, and decay, which is partially involved in a novel biological process termed phase separation. Through these regulatory mechanisms, m⁶A dictates gene expression in a dynamic and reversible manner and may play oncogenic, tumor suppressive or context-dependent roles in GI tumorigenesis. Therefore, regulators and effectors of m⁶A, as well as their modified substrates, represent a novel class of molecular targets for cancer treatments. In this review, we comprehensively summarize recent advances in this field and highlight research findings that documented key roles of RNA m⁶A modification in governing hallmarks of GI cancers. From a historical perspective, milestone findings in m⁶A machinery are integrated with a timeline of developing m⁶A targeting compounds. These available chemical compounds, as well as other approaches that target core components of the RNA m⁶A pathway hold promises for clinical translational to treat human GI cancers. Further investigation on several outstanding issues, e.g. how oncogenic insults may disrupt m⁶A homeostasis, and how m⁶A modification impacts on the tumor microenvironment, may dissect novel mechanisms underlying human tumorigenesis and identifies next-generation anti-cancer therapeutics. In this review, we discuss advances in our understanding of m⁶A RNA modification since its discovery in the 1970s to the latest progress in defining its potential clinic relevance. We summarize the molecular basis and roles of m⁶A regulators in the hallmarks of GI cancer and discuss their context-dependent functions. Furthermore, the identification and characterization of inhibitors or activators of m⁶A regulators and their potential anti-cancer effects are discussed. With the rapid growth in this field there is significant potential for developing m⁶A targeted therapy in GI cancers.

Absolute quantification of single-base m6A methylation in the mammalian transcriptome using GLORI

N6-methyladenosine (m6A) is the most abundant RNA modification in mammalian cells and the best-studied epitranscriptomic mark. Despite the development of various tools to map m6A, a transcriptome-wide method that enables absolute quantification of m6A at single-base resolution is lacking. Here we use glyoxal and nitrite-mediated deamination of unmethylated adenosines (GLORI) to develop an absolute m6A quantification method that is conceptually similar to bisulfite-sequencing-based quantification of DNA 5-methylcytosine. We apply GLORI to quantify the m6A methylomes of mouse and human cells and reveal clustered m6A modifications with differential distribution and stoichiometry. In addition, we characterize m6A dynamics under stress and examine the quantitative landscape of m6A modification in gene expression regulation. GLORI is an unbiased, convenient method for the absolute quantification of the m6A methylome.

Mutual regulation between N6-methyladenosine (m6A) modification and circular RNAs in cancer: impacts on therapeutic resistance

The resistance of tumor cells to therapy severely impairs the efficacy of treatment, leading to recurrence and metastasis of various cancers. Clarifying the underlying mechanisms of therapeutic resistance may provide new strategies for overcoming cancer resistance. N6-methyladenosine (m6A) is the most prevalent RNA modification in eukaryotes, and is involved in the regulation of RNA splicing, translation, transport, degradation, stability and processing, thus affecting several physiological processes and cancer progression. As a novel type of multifunctional non-coding RNAs (ncRNAs), circular RNAs (circRNAs) have been demonstrated to play vital roles in anticancer therapy. Currently, accumulating studies have revealed the mutual regulation of m6A modification and circRNAs, and their interaction can further influence the sensitivity of cancer treatment. In this review, we mainly summarized the recent advances of m6A modification and circRNAs in the modulation of cancer therapeutic resistance, as well as their interplay and potential mechanisms, providing promising insights and future directions in reversal of therapeutic resistance in cancer.

Towards a comprehensive understanding of RNA deamination: synthesis and properties of xanthosine-modified RNA

Nucleobase deamination, such as A-to-I editing, represents an important posttranscriptional modification of RNA. When deamination affects guanosines, a xanthosine (X) containing RNA is generated. However, the biological significance and chemical consequences on RNA are poorly understood. We present a comprehensive study on the preparation and biophysical properties of X-modified RNA. Thermodynamic analyses revealed that base pairing strength is reduced to a level similar to that observed for a G•U replacement. Applying NMR spectroscopy and X-ray crystallography, we demonstrate that X can form distinct wobble geometries with uridine depending on the sequence context. In contrast, X pairing with cytidine occurs either through wobble geometry involving protonated C or in Watson–Crick-like arrangement. This indicates that the different pairing modes are of comparable stability separated by low energetic barriers for switching. Furthermore, we demonstrate that the flexible pairing properties directly affect the recognition of X-modified RNA by reverse transcription enzymes. Primer extension assays and PCR-based sequencing analysis reveal that X is preferentially read as G or A and that the ratio depends on the type of reverse transcriptase. Taken together, our results elucidate important properties of X-modified RNA paving the way for future studies on its biological significance.

Rapid Determination of RNA Modifications in Consensus Motifs by Nuclease Protection with Ion-Tagged Oligonucleotide Probes and Matrix-Assisted Laser Desorption Ionization Mass Spectrometry

The reversible and substoichiometric modification of RNA has recently emerged as an additional layer of translational regulation in normal biological function and disease. Modifications are often enzymatically deposited in and removed from short (~5 nt) consensus motif sequences to carefully control the translational output of the cell. Although characterization of modification occupancy at consensus motifs can be accomplished using RNA sequencing methods, these approaches are generally time-consuming and do not directly detect post-transcriptional modifications. Here, we present a nuclease protection assay coupled with matrix-assisted laser desorption ionization mass spectrometry (MALDI-MS) to rapidly characterize modifications in consensus motifs, such as GGACU, which frequently harbor N6-methyladenosine (m⁶A). While conventional nuclease protection methods rely on long (~30 nt) oligonucleotide probes that preclude the global assessment of consensus motif modification stoichiometry, we investigated a series of ion-tagged oligonucleotide (ITO) probes and found that a benzylimidazolium-functionalized ITO (ABzIM-ITO) conferred significantly improved nuclease resistance for GGACU targets. After optimizing the conditions of the nuclease protection assay, we applied the ITO and MALDI-MS-based method for determining the stoichiometry of GG(m⁶A)CU and GGACU in RNA mixtures. Overall, the ITO-based nuclease protection and MALDI-MS method constitutes a rapid and promising approach for determining modification stoichiometries of consensus motifs.

NT-seq: a chemical-based sequencing method for genomic methylome profiling

DNA methylation plays vital roles in both prokaryotes and eukaryotes. There are three forms of DNA methylation in prokaryotes: N⁶-methyladenine (6mA), N⁴-methylcytosine (4mC), and 5-methylcytosine (5mC). Although many sequencing methods have been developed to sequence specific types of methylation, few technologies can be used for efficiently mapping multiple types of methylation. Here, we present NT-seq for mapping all three types of methylation simultaneously. NT-seq reliably detects all known methylation motifs in two bacterial genomes and can be used for identifying de novo methylation motifs. NT-seq provides a simple and efficient solution for detecting multiple types of DNA methylation.

RNA marker modifications reveal the necessity for rigorous preparation protocols to avoid artifacts in epitranscriptomic analysis

The accurate definition of an epitranscriptome is endangered by artefacts resulting from RNA degradation after cell death, a ubiquitous yet little investigated process. By tracing RNA marker modifications through tissue preparation protocols, we identified a major blind spot from daily lab routine, that has massive impact on modification analysis in small RNAs. In particular, m6,6A and Am as co-varying rRNA marker modifications, appeared in small RNA fractions following rRNA degradation in vitro and in cellulo. Analysing mouse tissue at different time points post mortem, we tracked the progress of intracellular RNA degradation after cell death, and found it reflected in RNA modification patterns. Differences were dramatic between liver, where RNA degradation commenced immediately after death, and brain, yielding essentially undamaged RNA. RNA integrity correlated with low amounts of co-varying rRNA markers. Thus validated RNA preparations featured differentially modified tRNA populations whose information content allowed a distinction even among the related brain tissues cortex, cerebellum and hippocampus. Inversely, advanced cell death correlated with high rRNA marker content, and correspondingly little with the naïve state of living tissue. Therefore, unless RNA and tissue preparations are executed with utmost care, interpretation of modification patterns in tRNA and small RNA are prone to artefacts.

scDART-seq: Mapping m6A at the single-cell level

In this issue of Molecular Cell, Tegowski et al. established a single-cell DART-seq (scDART-seq) method for m⁶A site profiling and further revealed the fundamental features of m⁶A biology in single cells. The findings improve our understating of the intrinsic regulatory and functional mechanisms of m⁶A in a single-cell basis.

Transcriptome-wide profiling of N6-methyladenosine via a selective chemical labeling method

Studies of chemical modifications on RNA have ushered in the field of epitranscriptomics. N⁶-Methyladenosine (m⁶A) is the most typical RNA modification and is indispensable for basic biological processes. This study presents a chemical pulldown method (m⁶A-ORL-Seq) for transcriptome-wide profiling of m⁶A. Moreover, chemical labeling results in a specific reverse transcription (RT) truncation signature. This study has identified four thousand high-confidence m⁶A sites at single-base resolution in the human transcriptome. Unlike previously reported methods based on m⁶A-antibody or m⁶A-sensitive enzymes, the antibody/enzyme-free chemical method provides a new perspective for single-base m⁶A detection at the transcriptome level.

m⁶A-ORL-Seq: A chemical labeling method for transcriptome-wide m⁶A profiling.

RNA N6-methyladenosine in nonocular and ocular disease

N⁶ -methyladenosine (m⁶ A), the sixth N methylation of adenylate (A) in RNA, is the most abundant transcriptome modification in eukaryotic messenger RNA (mRNAs). m⁶ A modification exists in both coding mRNA and noncoding RNAs, and its functions are controlled by methyltransferase, demethylase, and m⁶ A reading proteins. Methylation modification of m⁶ A can regulate RNA cleavage, transport, stability, and expression. This review summarizes the enzymes involved in RNA m⁶ A methylation and the commonly used detection methods. The role of m⁶ A modification in physiological processes is described, and its impact on tumorigenesis, viral infection, and diabetes is further highlighted. Moreover, up-to-date knowledge of the implications of RNA m⁶ A modification in ocular diseases such as uveal melanoma and diabetic retinopathy is introduced. Clarifying the mechanism of RNA m⁶ A methylation will help elucidate the pathogenesis of various diseases, providing options for subsequent treatment.

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.

Long Non-Coding RNA Epigenetics

Long noncoding RNAs exceeding a length of 200 nucleotides play an important role in ensuring cell functions and proper organism development by interacting with cellular compounds such as miRNA, mRNA, DNA and proteins. However, there is an additional level of lncRNA regulation, called lncRNA epigenetics, in gene expression control. In this review, we describe the most common modified nucleosides found in lncRNA, 6-methyladenosine, 5-methylcytidine, pseudouridine and inosine. The biosynthetic pathways of these nucleosides modified by the writer, eraser and reader enzymes are important to understanding these processes. The characteristics of the individual methylases, pseudouridine synthases and adenine-inosine editing enzymes and the methods of lncRNA epigenetics for the detection of modified nucleosides, as well as the advantages and disadvantages of these methods, are discussed in detail. The final sections are devoted to the role of modifications in the most abundant lncRNAs and their functions in pathogenic processes.

Elucidating the Functions of Non-Coding RNAs from the Perspective of RNA Modifications

It is now commonly accepted that most of the mammalian genome is transcribed as RNA, yet less than 2% of such RNA encode for proteins. A majority of transcribed RNA exists as non-protein-coding RNAs (ncRNAs) with various functions. Because of the lack of sequence homologies among most ncRNAs species, it is difficult to infer the potential functions of ncRNAs by examining sequence patterns, such as catalytic domains, as in the case of proteins. Added to the existing complexity of predicting the functions of the ever-growing number of ncRNAs, increasing evidence suggests that various enzymes modify ncRNAs (e.g., ADARs, METTL3, and METTL14), which has opened up a new field of study called epitranscriptomics. Here, we examine the current status of ncRNA research from the perspective of epitranscriptomics.

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.