MODOMICS: a database of RNA modification pathways. 2021 update

Utilization of nanopore direct RNA sequencing to analyze viral RNA modifications

Modifications on viral RNAs (vRNAs), either genomic RNAs or RNA transcripts, have complex effects on the viral life cycle and cellular responses to viral infection. The advent of Oxford Nanopore Technologies Direct RNA Sequencing provides a new strategy for studying RNA modifications. To this end, multiple computational tools have been developed, but a systemic evaluation of their performance in mapping vRNA modifications is lacking. Here, 10 computational tools were tested using the Sindbis virus (SINV) RNAs isolated from infected mammalian (BHK-21) or mosquito (C6/36) cells, with in vitro-transcribed RNAs serving as modification-free control. Three single-mode approaches were shown to be inapplicable in the viral context, and three out of seven comparative methods required cutoff adjustments to reduce false-positive predictions. Utilizing optimized cutoffs, an integrated analysis of comparative tools suggested that the intersected predictions of Tombo_com and xPore were significantly enriched compared with the background. Consequently, a pipeline integrating Tombo_com and xPore was proposed for vRNA modification detection; the performance of which was supported by N6-methyladenosine prediction in severe acute respiratory syndrome coronavirus 2 RNAs using publicly available data. When applied to SINV RNAs, this pipeline revealed more intensive modifications in subgenomic RNAs than in genomic RNAs. Modified uridines were frequently identified, exhibiting substantive overlapping between vRNAs generated in different cell lines. On the other hand, the interpretation of other modifications remained unclear, underlining the limitations of the current computational tools despite their notable potential.IMPORTANCEComputational approaches utilizing Oxford Nanopore Technologies Direct RNA Sequencing data were almost exclusively designed to map eukaryotic epitranscriptomes. Therefore, extra caution must be exercised when using these tools to detect vRNA modifications, as in most cases, vRNA modification profiles should be regarded as unknown epitranscriptomes without prior knowledge. Here, we comprehensively evaluated the performance of 10 computational tools in detecting vRNA modification sites. All tested single-mode methods failed to differentiate native and in vitro-transcribed samples. Using optimized cutoff values, seven tested comparative tools generated very different predictions. An integrated analysis showed significant enrichment of Tombo_com and xPore predictions against the background. A pipeline for vRNA modification detection was proposed accordingly and applied to Sindbis virus RNAs. In conclusion, our study underscores the need for the careful application of computational tools to analyze viral epitranscriptomics. It also offers insights into alphaviral RNA modifications, although further validation is required.

Dissection of mRNA ac4C acetylation modifications in AC and Nr fruits: insights into the regulation of fruit ripening by ethylene

N4-acetylcytidine (ac4C) modification of mRNA has been shown to be present in plant RNAs, but its regulatory function in plant remains largely unexplored. In this study, we investigated the differentially expressed mRNAs, lncRNAs and acetylation modifications of mRNAs in tomato fruits from both genotypes. By comparing wild-type (AC) tomato and the ethylene receptor-mutant (Nr) tomato from mature green (MG) to six days after the breaker (Br6) stage, we identified differences in numerous key genes related to fruit ripening and observed the corresponding lncRNAs positively regulated the target genes expression. At the post-transcriptional level, the acetylation level decreased and increased in AC and Nr tomatoes from MG to Br6 stage, respectively. The integrated analysis of RNA-seq and ac4C-seq data revealed the potential positive role of acetylation modification in regulating gene expression. Furthermore, we found differential acetylation modifications of certain transcripts (ACO, ETR, ERF, PG, CesA, β-Gal, GAD, AMY, and SUS) in AC and Nr fruits which may explain the differences in ethylene production, fruit texture, and flavor during their ripening processes. The present study provides new insights into the molecular mechanisms by which acetylation modification differentially regulates the ripening process of wild-type and mutant tomato fruits deficient in ethylene signaling.

Inosine-Induced Base Pairing Diversity during Reverse Transcription

A-to-I editing catalyzed by adenosine deaminase acting on RNAs impacts numerous physiological and biochemical processes that are essential for cellular functions and is a big contributor to the infectivity of certain RNA viruses. The outcome of this deamination leads to changes in the eukaryotic transcriptome functionally resembling A-G transitions since inosine preferentially pairs with cytosine. Moreover, hyper-editing or multiple A to G transitions in clusters were detected in measles virus. Inosine modifications either directly on viral RNA or on cellular RNA can have antiviral or pro-viral repercussions. While many of the significant roles of inosine in cellular RNAs are well understood, the effects of hyper-editing of A to I on viral polymerase activity during RNA replication remain elusive. Moreover, biological strategies such as molecular cloning and RNA-seq for transcriptomic interrogation rely on RT-polymerase chain reaction with little to no emphasis placed on the first step, reverse transcription, which may reshape the sequencing results when hypermodification is present. In this study, we systematically explore the influence of inosine modification, varying the number and position of inosines, on decoding outcomes using three different reverse transcriptases (RTs) followed by standard Sanger sequencing. We find that inosine alone or in clusters can differentially affect the RT activity. To gain structural insights into the accommodation of inosine in the polymerase site of HIV-1 reverse transcriptase (HIV-1-RT) and how this structural context affects the base pairing rules for inosine, we performed molecular dynamics simulations of the HIV-1-RT. The simulations highlight the importance of the protein-nucleotide interaction as a critical factor in deciphering the base pairing behavior of inosine clusters. This effort sets the groundwork for decrypting the physiological significance of inosine and linking the fidelity of reverse transcriptase and the possible diverse transcription outcomes of cellular RNAs and/or viral RNAs where hyper-edited inosines are present in the transcripts.

Recent Advances in RNA-Targeted Cancer Therapy

Ribonucleic acid (RNA) plays a pivotal role in gene regulation and protein biosynthesis. Interfering the physiological function of key RNAs to induce cell apoptosis holds great promise for cancer treatment. Many RNA-targeted anti-cancer strategies have emerged continuously. Among them, RNA interference (RNAi) has been recognized as a promising therapeutic modality for various disease treatments. Nevertheless, the primary obstacle in siRNA delivery-escaping the endosome and crossing the plasma membrane severely impedes its therapeutic potential. Thus far, a variety of nanosystems as well as carrier-free bioconjugation for siRNA delivery have been developed and employed to enhance the drug delivery and anti-tumor efficiency. Besides, the use of small molecules to target specific RNA structures and disrupt their function, along with the covalent modification of RNA, has also drawn tremendous attention recently owing to high therapeutic efficacy. In this review, we will provide an overview of recent progress in RNA-targeted cancer therapy including various siRNA delivery strategies, RNA-targeting small molecules, and newly emerged covalent RNA modification. Finally, challenges and future perspectives faced in this research field will be discussed.

Cell death or survival: Insights into the role of mRNA translational control

Cellular stress is an intrinsic part of cell physiology that underlines cell survival or death. The ability of mammalian cells to regulate global protein synthesis (aka translational control) represents a critical, yet underappreciated, layer of regulation during the stress response. Various cellular stress response pathways monitor conditions of cell growth and subsequently reshape the cellular translatome to optimize translational outputs. On the molecular level, such translational reprogramming involves an intricate network of interactions between translation machinery, RNA-binding proteins, mRNAs, and non-protein coding RNAs. In this review, we will discuss molecular mechanisms, signaling pathways, and targets of translational control that contribute to cellular adaptation to stress and to cell survival or death.

Contribution of tRNA sequence and modifications to the decoding preferences of E. coli and M. mycoides tRNAGlyUCC for synonymous glycine codons

tRNA superwobbling, used by certain bacteria and organelles, is an intriguing decoding concept in which a single tRNA isoacceptor is used to decode all synonymous codons of a four-fold degenerate codon box. While Escherichia coli relies on three tRNAGly isoacceptors to decode the four glycine codons (GGN), Mycoplasma mycoides requires only a single tRNAGly. Both organisms express tRNAGly with the anticodon UCC, which are remarkably similar in sequence but different in their decoding ability. By systematically introducing mutations and altering the number and type of tRNA modifications using chemically synthesized tRNAs, we elucidated the contribution of individual nucleotides and chemical groups to decoding by the E. coli and M. mycoides tRNAGly. The tRNA sequence was identified as the key factor for superwobbling, revealing the T-arm sequence as a novel pivotal element. In addition, the presence of tRNA modifications, although not essential for providing superwobbling, was shown to delicately fine-tune and balance the decoding of synonymous codons. This emphasizes that the tRNA sequence and its modifications together form an intricate system of high complexity that is indispensable for accurate and efficient decoding.

Internal RNA 2'-O-methylation on the HIV-1 genome impairs reverse transcription

Viral RNA genomes are modified by epitranscriptomic marks, including 2'-O-methylation that is added by cellular or viral methyltransferases. 2'-O-Methylation modulates RNA structure, function and discrimination between self- and non-self-RNA by innate immune sensors such as RIG-I-like receptors. This is illustrated by human immunodeficiency virus type-1 (HIV-1) that decorates its RNA genome through hijacking the cellular FTSJ3 2'-O-methyltransferase, thereby limiting immune sensing and interferon production. However, the impact of such an RNA modification during viral genome replication is poorly understood. Here we show by performing endogenous reverse transcription on methylated or hypomethylated HIV-1 particles, that 2'-O-methylation negatively affects HIV-1 reverse transcriptase activity. Biochemical assays confirm that RNA 2'-O-methylation impedes reverse transcriptase activity, especially at low dNTP concentrations reflecting those in quiescent cells, by reducing nucleotide incorporation efficiency and impairing translocation. Mutagenesis highlights K70 as a critical amino acid for the reverse transcriptase to bypass 2'-O-methylation. Hence, the observed antiviral effect due to viral RNA 2'-O-methylation antagonizes the FTSJ3-mediated proviral effects, suggesting the fine-tuning of RNA methylation during viral replication.

Advances in the Structural and Functional Understanding of m1A RNA Modification

ConspectusRNA modification is a co- or post-transcriptional process by which specific nucleotides are chemically altered by enzymes after their initial incorporation into the RNA chain, expanding the chemical and functional diversity of RNAs. Our understanding of RNA modifications has changed dramatically in recent years. In the past decade, RNA methyltransferases (MTases) have been highlighted in numerous clinical studies and disease models, modifications have been found to be dynamically regulated by demodification enzymes, and significant technological advances have been made in the fields of RNA sequencing, mass spectrometry, and structural biology. Among RNAs, transfer RNAs (tRNAs) exhibit the greatest diversity and density of post-transcriptional modifications, which allow for potential cross-talks and regulation during their incorporation. N1-methyladenosine (m1A) modification is found in tRNAs at positions 9, 14, 16, 22, 57, and 58, depending on the tRNA and organism.Our laboratory has used and developed a large panel of tools to decipher the different mechanisms used by m1A tRNA MTases to recognize and methylate tRNA. We have solved the structures of TrmI from Thermus thermophilus (m1A58), TrmK from Bacillus subtilis (m1A22), and human TRMT10C (m1A9). These MTases do not share the same structure or organization to recognize tRNAs, but they all modify an adenosine, forming a non-Watson-Crick (WC) interaction. For TrmK, nuclear magnetic resonance (NMR) chemical shift mapping of the binding interface between TrmK and tRNASer was invaluable to build a TrmK/tRNA model, where both domains of TrmK participate in the binding of a full-length L-shaped tRNA and where the non-WC purine 13-A22 base pair positions the A22 N1-atom close to the methyl of the S-adenosyl-l-methionine (SAM) TrmK cofactor. For TRMT10C, cryoEM structures showed the MTase poised to N1-methylate A9 or G9 in tRNA and revealed different steps of tRNA maturation, where TRMT10C acts as a tRNA binding platform for sequential docking of each maturation enzyme. This work confers a role for TRMT10C in tRNA quality control and provides a framework to understand the link between mitochondrial tRNA maturation dysfunction and diseases.Methods to directly detect the incorporation of modifications during tRNA biosynthesis are rare and do not provide easy access to the temporality of their introduction. To this end, we have introduced time-resolved NMR to monitor tRNA maturation in the cellular environment. Combined with genetic and biochemical approaches involving the synthesis of specifically modified tRNAs, our methodology revealed that some modifications are incorporated in a defined sequential order, controlled by cross-talks between modification events. In particular, a strong modification circuit, namely Ψ55 → m5U54 → m1A58, controls the modification process in the T-arm of yeast elongator tRNAs. Conversely, we showed that m1A58 is efficiently introduced on unmodified initiator tRNAiMet without the need of any prior modification. Two distinct pathways are therefore followed for m1A58 incorporation in elongator and initiator tRNAs.We are undoubtedly entering an exciting period for the elucidation of the functions of RNA modifications and the intricate mechanisms by which modification enzymes identify and alter their RNA substrates. These are promising directions for the field of epitranscriptomics.

General Principles and Limitations for Detection of RNA Modifications by Sequencing

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

Mass Spectrometry-Based Pipeline for Identifying RNA Modifications Involved in a Functional Process: Application to Cancer Cell Adaptation

Cancer onset and progression are known to be regulated by genetic and epigenetic events, including RNA modifications (a.k.a. epitranscriptomics). So far, more than 150 chemical modifications have been described in all RNA subtypes, including messenger, ribosomal, and transfer RNAs. RNA modifications and their regulators are known to be implicated in all steps of post-transcriptional regulation. The dysregulation of this complex yet delicate balance can contribute to disease evolution, particularly in the context of carcinogenesis, where cells are subjected to various stresses. We sought to discover RNA modifications involved in cancer cell adaptation to inhospitable environments, a peculiar feature of cancer stem cells (CSCs). We were particularly interested in the RNA marks that help the adaptation of cancer cells to suspension culture, which is often used as a surrogate to evaluate the tumorigenic potential. For this purpose, we designed an experimental pipeline consisting of four steps: (1) cell culture in different growth conditions to favor CSC survival; (2) simultaneous RNA subtype (mRNA, rRNA, tRNA) enrichment and RNA hydrolysis; (3) the multiplex analysis of nucleosides by LC-MS/MS followed by statistical/bioinformatic analysis; and (4) the functional validation of identified RNA marks. This study demonstrates that the RNA modification landscape evolves along with the cancer cell phenotype under growth constraints. Remarkably, we discovered a short epitranscriptomic signature, conserved across colorectal cancer cell lines and associated with enrichment in CSCs. Functional tests confirmed the importance of selected marks in the process of adaptation to suspension culture, confirming the validity of our approach and opening up interesting prospects in the field.

Stack-DHUpred: Advancing the accuracy of dihydrouridine modification sites detection via stacking approach

Dihydrouridine (DHU, D) is one of the most abundant post-transcriptional uridine modifications found in tRNA, mRNA, and snoRNA, closely associated with disease pathogenesis and various biological processes in eukaryotes. Identifying D sites is important for understanding the modification mechanisms and/or epigenetic regulation. However, biological experiments for detecting D sites are time-consuming and expensive. Given these challenges, computational methods have been developed for accurately identifying the D sites in genome-wide datasets. However, existing methods have some limitations, and their prediction performance needs to be improved. In this work, we have developed a new computational predictor for accurately identifying D sites called Stack-DHUpred. Briefly, we trained 66 baseline models or single-feature models by connecting six machine learning classifiers with eleven different feature encoding methods and stacked different baseline models to build stacked ensemble learning models. Subsequently, the optimal combination of the baseline models was identified for the construction of the final stacked model. Remarkably, the Stack-DHUpred outperformed the existing predictors on our new independent dataset, indicating that the stacking approach significantly improved the prediction performance. We have made Stack-DHUpred available to the public through a web server (http://kurata35.bio.kyutech.ac.jp/Stack-DHUpred) and a standalone program (https://github.com/kuratahiroyuki/Stack-DHUpred). We believe that Stack-DHUpred will be a valuable tool for accelerating the discovery of D modifications and understanding their role in post-transcriptional regulation.

Heterogeneity of chemical modifications on RNA

The chemical modifications of RNAs broadly impact almost all cellular events and influence various diseases. The rapid advance of sequencing and other technologies opened the door to global methods for profiling all RNA modifications, namely the "epitranscriptome." The mapping of epitranscriptomes in different cells and tissues unveiled that RNA modifications exhibit extensive heterogeneity, in type, amount, and in location. In this mini review, we first introduce the current understanding of modifications on major types of RNAs and the methods that enabled their discovery. We next discuss the tissue and cell heterogeneity of RNA modifications and briefly address the limitations of current technologies. With much still remaining unknown, the development of the epitranscriptomic field lies in the further developments of novel technologies.

Tissue-specific RNA methylation prediction from gene expression data using sparse regression models

N6-methyladenosine (m6A) is a highly prevalent and conserved post-transcriptional modification observed in mRNA and long non-coding RNA (lncRNA). Identifying potential m6A sites within RNA sequences is crucial for unraveling the potential influence of the epitranscriptome on biological processes. In this study, we introduce Exp2RM, a novel approach that formulates single-site-based tissue-specific elastic net models for predicting tissue-specific methylation levels utilizing gene expression data. The resulting ensemble model demonstrates robust predictive performance for tissue-specific methylation levels, with an average R-squared value of 0.496 and a median R-squared value of 0.482 across all 22 human tissues. Since methylation distribution varies among tissues, we trained the model to incorporate similar patterns, significantly improves accuracy with the median R-squared value increasing to 0.728. Additonally, functional analysis reveals Exp2RM's ability to capture coefficient genes in relevant biological processes. This study emphasizes the importance of tissue-specific methylation distribution in enhancing prediction accuracy and provides insights into the functional implications of methylation sites.

RNA modifications in physiology and disease: towards clinical applications

The ability of chemical modifications of single nucleotides to alter the electrostatic charge, hydrophobic surface and base pairing of RNA molecules is exploited for the clinical use of stable artificial RNAs such as mRNA vaccines and synthetic small RNA molecules - to increase or decrease the expression of therapeutic proteins. Furthermore, naturally occurring biochemical modifications of nucleotides regulate RNA metabolism and function to modulate crucial cellular processes. Studies showing the mechanisms by which RNA modifications regulate basic cell functions in higher organisms have led to greater understanding of how aberrant RNA modification profiles can cause disease in humans. Together, these basic science discoveries have unravelled the molecular and cellular functions of RNA modifications, have provided new prospects for therapeutic manipulation and have led to a range of innovative clinical approaches.

Update on N6-methyladenosine methylation in obesity-related diseases

Obesity is a chronic metabolic disease that is closely related to type 2 diabetes mellitus, cardiovascular diseases, nonalcoholic fatty liver disease, obstructive sleep apnea, and osteoarthritis. The prevalence of obesity is increasing rapidly every year and is recognized as a global public health problem. In recent years, the role of epigenetics in the development of obesity and related diseases has been recognized and is currently a research hotspot. N6-methyladenosine (m6A) methylation is the most abundant epigenetic modification in the eukaryotic RNA, including mRNA and noncoding RNA. Several studies have shown that the m6A modifications in the target mRNA and the corresponding m6A regulators play a significant role in lipid metabolism and are strongly associated with the pathogenesis of obesity-related diseases. In this review, the latest research findings regarding the role of m6A methylation in obesity and related metabolic diseases are summarized. The authors' aim is to highlight evidence that suggests the clinical utility of m6A modifications and the m6A regulators as novel early prediction biomarkers and precision therapeutics for obesity and obesity-related diseases.

Protection-Free, Two-step Synthesis of C5-C Functionalized Pyrimidine Nucleosides

A simple, reliable, and efficient method for the gram-scale chemical synthesis of pyrimidine nucleosides functionalized with C5-carboxyl, nitrile, ester, amide, or amidine, starting from unprotected uridine and cytidine, is described. The protocol involves the synthesis of 5-trifluoromethyluridine and 5-trifluoromethylcytidine with Langlois reagent (CF3 SO2 Na) in the presence of tert-butyl hydroperoxide and subsequent transformation of the CF3 group to the C5-C 'carbon substituents' under alkaline conditions. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Synthesis and characterization of 5-trifluoromethyluridine (5-CF3 U) and 5-trifluoromethylcytidine (5-CF3 C) Basic Protocol 2: Conversion of 5-CF3 U and 5-CF3 C to several C5-substituted ribonucleosides.

The enigmatic epitranscriptome of bacteriophages: putative RNA modifications in viral infections

RNA modifications play essential roles in modulating RNA function, stability, and fate across all kingdoms of life. The entirety of the RNA modifications within a cell is defined as the epitranscriptome. While eukaryotic RNA modifications are intensively studied, understanding bacterial RNA modifications remains limited, and knowledge about bacteriophage RNA modifications is almost nonexistent. In this review, we shed light on known mechanisms of bacterial RNA modifications and propose how this knowledge might be extended to bacteriophages. We build hypotheses on enzymes potentially responsible for regulating the epitranscriptome of bacteriophages and their host. This review highlights the exciting prospects of uncovering the unexplored field of bacteriophage epitranscriptomics and its potential role to shape bacteriophage-host interactions.

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.

A low-cost, low-input method establishment for m6A MeRIP-seq

N6-methyladenosine (m6A) is a highly prevalent modification found in mammal mRNA molecules that plays a crucial role in the regulation of cellular function. m6A RNA immunoprecipitation sequencing (MeRIP-seq) has been frequently used in transcriptomics research to identify the location of m6A. MABE572 (Millipore) is the most widely utilized and efficient anti-m6A antibody for MeRIP-seq. However, due to the high dose and price of this antibody, which has also been taken off the market, we discovered that CST's anti-m6A antibody can be used instead of MABE572 to map the m6A transcriptome. In the present study, we performed different concentrations of the CST anti-m6A antibodies with the corresponding initiation RNA of HEK293T cells, 2.5 μg antibody with 1 μg total RNA, 1.25 μg antibody with 0.5 μg total RNA, and 1.25 μg antibody with 0.1 μg total RNA. By comparing the m6A peak calling, enriched motifs, alternative splicing events, and nuclear transcripts modified by m6A between the CST and Millipore libraries, it was found that the CST library presented similar data to Millipore, even at incredibly low doses. The volume and cost of antibodies are significantly reduced by this refined MeRIP-seq using CST antibody, making it convenient to map future large-scale sample m6A methylation.

Effects of N6-methyladenosine modification on metabolic reprogramming in digestive tract tumors

N6-methyladenosine (m6A), the most abundant RNA modification within cells, participates in various biological and pathological processes, including self-renewal, invasion and proliferation, drug resistance, and stem cell characteristics. The m6A methylation plays a crucial role in tumors by regulating multiple RNA processes such as transcription, processing, and translation. Three protein types are primarily involved in m6A methylation: methyltransferases (such as METTL3, METTL14, ZC3H13, and KIAA1429), demethylases (such as FTO, ALKBH5), and RNA-binding proteins (such as the family of YTHDF, YTHDC1, YTHDC2, and IGF2BPs). Various metabolic pathways are reprogrammed in digestive tumors to meet the heightened growth demands and sustain cellular functionality. Recent studies have highlighted the extensive impact of m6A on the regulation of digestive tract tumor metabolism, further modulating tumor initiation and progression. Our review aims to provide a comprehensive understanding of the expression patterns, functional roles, and regulatory mechanisms of m6A in digestive tract tumor metabolism-related molecules and pathways. The characterization of expression profiles of m6A regulatory factors and in-depth studies on m6A methylation in digestive system tumors may provide new directions for clinical prediction and innovative therapeutic interventions.

RNA N6-Methyladenosine Affects Copper-Induced Oxidative Stress Response in Arabidopsis thaliana

Recently, post-transcriptional regulation of mRNA mediated by N6-methyladenosine (m6A) has been found to have profound effects on transcriptome regulation during plant responses to various abiotic stresses. However, whether this RNA modification can affect an oxidative stress response in plants has not been studied. To assess the role of m6A modifications during copper-induced oxidative stress responses, m6A-IP-seq was performed in Arabidopsis seedlings exposed to high levels of copper sulfate. This analysis revealed large-scale shifts in this modification on the transcripts most relevant for oxidative stress. This altered epitranscriptomic mark is known to influence transcript abundance and translation; therefore we scrutinized these possibilities. We found an increased abundance of copper-enriched m6A-containing transcripts. Similarly, we also found increased ribosome occupancy of copper-enriched m6A-containing transcripts, specifically those encoding proteins involved with stress responses relevant to oxidative stressors. Furthermore, the significance of the m6A epitranscriptome on plant oxidative stress tolerance was uncovered by assessing germination and seedling development of the mta (N6-methyladenosine RNA methyltransferase A mutant complemented with ABI3:MTA) mutant exposed to high copper treatment. These analyses suggested hypersensitivity of the mta mutant compared to the wild-type plants in response to copper-induced oxidative stress. Overall, our findings suggest an important role for m6A in the oxidative stress response of Arabidopsis.

The Bacillus subtilis ywbD gene encodes RlmQ, the 23S rRNA methyltransferase forming m7G2574 in the A-site of the peptidyl transferase center

Ribosomal RNA contains many posttranscriptionally modified nucleosides, particularly in the functional parts of the ribosome. The distribution of these modifications varies from one organism to another. In Bacillus subtilis, the model organism for Gram-positive bacteria, mass spectrometry experiments revealed the presence of 7-methylguanosine (m7G) at position 2574 of the 23S rRNA, which lies in the A-site of the peptidyl transferase center of the large ribosomal subunit. Testing several m7G methyltransferase candidates allowed us to identify the RlmQ enzyme, encoded by the ywbD open reading frame, as the MTase responsible for this modification. The enzyme methylates free RNA and not ribosomal 50S or 70S particles, suggesting that modification occurs in the early steps of ribosome biogenesis.

AlkB-Facilitated Demethylation Enables Quantitative and Site-Specific Detection of Dual Methylation of Adenosine in RNA

RNA molecules undergo various chemical modifications that play critical roles in a wide range of biological processes. N6,N6-Dimethyladenosine (m6,6A) is a conserved RNA modification and is essential for the processing of rRNA. To gain a deeper understanding of the functions of m6,6A, site-specific and accurate quantification of this modification in RNA is indispensable. In this study, we developed an AlkB-facilitated demethylation (AD-m6,6A) method for the site-specific detection and quantification of m6,6A in RNA. The N6,N6-dimethyl groups in m6,6A can cause reverse transcription to stall at the m6,6A site, resulting in truncated cDNA. However, we found that Escherichia coli AlkB demethylase can effectively demethylate m6,6A in RNA, generating full-length cDNA from AlkB-treated RNA. By quantifying the amount of full-length cDNA produced using quantitative real-time PCR, we were able to achieve site-specific detection and quantification of m6,6A in RNA. Using the AD-m6,6A method, we successfully detected and quantified m6,6A at position 1851 of 18S rRNA and position 937 of mitochondrial 12S rRNA in human cells. Additionally, we found that the level of m6,6A at position 1007 of mitochondrial 12S rRNA was significantly reduced in lung tissues from sleep-deprived mice compared with control mice. Overall, the AD-m6,6A method provides a valuable tool for easy, accurate, quantitative, and site-specific detection of m6,6A in RNA, which can aid in uncovering the functions of m6,6A in human diseases.

Role of RNA Modifications, Especially m6A, in Aflatoxin Biosynthesis of Aspergillus flavus

RNA modifications play key roles in eukaryotes, but the functions in Aspergillus flavus are still unknown. Temperature has been reported previously to be a critical environmental factor that regulates the aflatoxin production of A. flavus, but much remains to be learned about the molecular networks. Here, we demonstrated that 12 kinds of RNA modifications in A. flavus were significantly changed under 29 °C compared to 37 °C incubation; among them, m6A was further verified by a colorimetric method. Then, the transcriptome-wide m6A methylome and m6A-altered genes were comprehensively illuminated through methylated RNA immunoprecipitation sequencing and RNA sequencing, from which 22 differentially methylated and expressed transcripts under 29 °C were screened out. It is especially notable that AFCA_009549, an aflatoxin biosynthetic pathway gene (aflQ), and the m6A methylation of its 332nd adenine in the mRNA significantly affect aflatoxin biosynthesis in A. flavus both on media and crop kernels. The content of sterigmatocystin in both ΔaflQ and aflQA332C strains was significantly higher than that in the WT strain. Together, these findings reveal that RNA modifications are associated with secondary metabolite biosynthesis of A. flavus.

Mapping m6A Sites on HIV-1 RNA Using Oligonucleotide LC-MS/MS

The biological significance of chemical modifications to the ribonucleic acid (RNA) of human immunodeficiency virus type-1 (HIV-1) has been recognized. However, our understanding of the site-specific and context-dependent roles of these chemical modifications remains limited, primarily due to the absence of nucleotide-resolution mapping of modification sites. In this study, we present a method for achieving nucleotide-resolution mapping of chemical modification sites on HIV-1 RNA using liquid chromatography and tandem mass spectrometry (LC-MS/MS). LC-MS/MS, a powerful tool capable of directly analyzing native RNAs, has proven effective for mapping RNA modifications in small RNA molecules, including ribosomal RNA and transfer RNA. However, longer RNAs have posed challenges, such as the 9 Kb HIV-1 virion RNA, due to the complexity of and ambiguity in mass differences among RNase T1-cleaved RNA fragments in LC-MS/MS data. Here, we introduce a new target RNA enrichment method to isolate small local RNA fragments of HIV-1 RNA that potentially harbor site-specific N6-methyladenosine (m6A) modifications. In our initial trial, we used target-specific DNA probes only and encountered insufficient RNA fragmentation due to inefficient S1 digestion near the target site. Recognizing that inefficient S1 digestion by HIV-1 RNA is likely due to the formation of secondary structures in proximity to the target site, we designed multiple DNA probes annealing to various sites of HIV-1 RNA to better control the structures of RNA substrates for S1 digestion. The use of these non-target DNA probes significantly improved the isolation of more homogeneous target RNA fragments of approximately 50 bases in length. Oligonucleotide LC-MS/MS analysis of these isolated target RNA fragments successfully separated and detected both m6A-methylated and non-methylated oligomers at the two m6A-predicted sites. The principle of this new target enrichment strategy holds promise and should be broadly applicable to the analysis of any lengthy RNA that was previously deemed infeasible for investigation using oligonucleotide LC-MS/MS.

RMBase v3.0: decode the landscape, mechanisms and functions of RNA modifications

Although over 170 chemical modifications have been identified, their prevalence, mechanism and function remain largely unknown. To enable integrated analysis of diverse RNA modification profiles, we have developed RMBase v3.0 (http://bioinformaticsscience.cn/rmbase/), a comprehensive platform consisting of eight modules. These modules facilitate the exploration of transcriptome-wide landscape, biogenesis, interactome and functions of RNA modifications. By mining thousands of epitranscriptome datasets with novel pipelines, the 'RNA Modifications' module reveals the map of 73 RNA modifications of 62 species. the 'Genes' module allows to retrieve RNA modification profiles and clusters by gene and transcript. The 'Mechanisms' module explores 23 382 enzyme-catalyzed or snoRNA-guided modified sites to elucidate their biogenesis mechanisms. The 'Co-localization' module systematically formulates potential correlations between 14 histone modifications and 6 RNA modifications in various cell-lines. The 'RMP' module investigates the differential expression profiles of 146 RNA-modifying proteins (RMPs) in 18 types of cancers. The 'Interactome' integrates the interactional relationships between 73 RNA modifications with RBP binding events, miRNA targets and SNPs. The 'Motif' illuminates the enriched motifs for 11 types of RNA modifications identified from epitranscriptome datasets. The 'Tools' introduces a novel web-based 'modGeneTool' for annotating modifications. Overall, RMBase v3.0 provides various resources and tools for studying RNA modifications.

PRMD: an integrated database for plant RNA modifications

The scope and function of RNA modifications in model plant systems have been extensively studied, resulting in the identification of an increasing number of novel RNA modifications in recent years. Researchers have gradually revealed that RNA modifications, especially N6-methyladenosine (m6A), which is one of the most abundant and commonly studied RNA modifications in plants, have important roles in physiological and pathological processes. These modifications alter the structure of RNA, which affects its molecular complementarity and binding to specific proteins, thereby resulting in various of physiological effects. The increasing interest in plant RNA modifications has necessitated research into RNA modifications and associated datasets. However, there is a lack of a convenient and integrated database with comprehensive annotations and intuitive visualization of plant RNA modifications. Here, we developed the Plant RNA Modification Database (PRMD; http://bioinformatics.sc.cn/PRMD and http://rnainformatics.org.cn/PRMD) to facilitate RNA modification research. This database contains information regarding 20 plant species and provides an intuitive interface for displaying information. Moreover, PRMD offers multiple tools, including RMlevelDiff, RMplantVar, RNAmodNet and Blast (for functional analyses), and mRNAbrowse, RNAlollipop, JBrowse and Integrative Genomics Viewer (for displaying data). Furthermore, PRMD is freely available, making it useful for the rapid development and promotion of research on plant RNA modifications.

MODOMICS: a database of RNA modifications and related information. 2023 update

The MODOMICS database was updated with recent data and now includes new data types related to RNA modifications. Changes to the database include an expanded modification catalog, encompassing both natural and synthetic residues identified in RNA structures. This addition aids in representing RNA sequences from the RCSB PDB database more effectively. To manage the increased number of modifications, adjustments to the nomenclature system were made. Updates in the RNA sequences section include the addition of new sequences and the reintroduction of sequence alignments for tRNAs and rRNAs. The protein section was updated and connected to structures from the RCSB PDB database and predictions by AlphaFold. MODOMICS now includes a data annotation system, with 'Evidence' and 'Estimated Reliability' features, offering clarity on data support and accuracy. This system is open to all MODOMICS entries, enhancing the accuracy of RNA modification data representation. MODOMICS is available at https://iimcb.genesilico.pl/modomics/.

m6A-Atlas v2.0: updated resources for unraveling the N6-methyladenosine (m6A) epitranscriptome among multiple species

N 6-Methyladenosine (m6A) is one of the most abundant internal chemical modifications on eukaryote mRNA and is involved in numerous essential molecular functions and biological processes. To facilitate the study of this important post-transcriptional modification, we present here m6A-Atlas v2.0, an updated version of m6A-Atlas. It was expanded to include a total of 797 091 reliable m6A sites from 13 high-resolution technologies and two single-cell m6A profiles. Additionally, three methods (exomePeaks2, MACS2 and TRESS) were used to identify >16 million m6A enrichment peaks from 2712 MeRIP-seq experiments covering 651 conditions in 42 species. Quality control results of MeRIP-seq samples were also provided to help users to select reliable peaks. We also estimated the condition-specific quantitative m6A profiles (i.e. differential methylation) under 172 experimental conditions for 19 species. Further, to provide insights into potential functional circuitry, the m6A epitranscriptomics were annotated with various genomic features, interactions with RNA-binding proteins and microRNA, potentially linked splicing events and single nucleotide polymorphisms. The collected m6A sites and their functional annotations can be freely queried and downloaded via a user-friendly graphical interface at: http://rnamd.org/m6a.

The role of the methyltransferase METTL3 in prostate cancer: a potential therapeutic target

The incidence of prostate cancer (PCa), the most prevalent malignancy, is currently at the forefront. RNA modification is a subfield of the booming field of epigenetics. To date, more than 170 types of RNA modifications have been described, and N6-methyladenosine (m6A) is the most abundant and well-characterized internal modification of mRNAs involved in various aspects of cancer progression. METTL3, the first identified key methyltransferase, regulates human mRNA and non-coding RNA expression in an m6A-dependent manner. This review elucidates the biological function and role of METTL3 in PCa and discusses the implications of METTL3 as a potential therapeutic target for future research directions and clinical applications.

M6A plays a potential role in carotid atherosclerosis by modulating immune cell modification and regulating aging-related genes

RNA N6-methyladenosine (m6A) regulators play essential roles in diverse biological processes, including immune responses. Mounting evidence suggests that their dysregulation is intricately linked to numerous diseases. However, the role of m6A-associated genes in carotid atherosclerosis and their relationship with aging and immune cells remain unclear. Analyze the expression profiles of m6A-related genes in carotid atherosclerosis-related datasets. Based on the expression patterns of m6A-related genes, perform consistent clustering analysis of carotid atherosclerosis samples and investigate associated immune cell infiltration patterns and aging characteristics. Develop an m6A prediction model specific to carotid atherosclerosis and analyze the relationships between immune cells infiltration and aging features. The m6A methylation modification level exhibited a substantial decrease in early-stage carotid atherosclerosis samples compared to late-stage carotid atherosclerosis samples. Subsequently, two distinct m6A subtypes were defined through consensus clustering analysis, with the lower m6A modification level group showing associations with heightened immune cell infiltration and increased expression of aging-related genes. A model composed of five m6A-related genes was formulated, and the results indicated that this model possesses effective predictive and therapeutic capabilities for carotid atherosclerosis. Furthermore, the downregulation of YTHDC1 expression resulted in elevated expression of inflammatory factors and a decrease in the expression of the aging-related gene RGN. Single-cell data analysis suggests that the reduced expression of YTHDC1 may decrease the degradation of inflammation-related factors in macrophages, leading to a highly inflammatory state in the carotid artery wall. Furthermore, the sustained release of inflammatory factors may increase the expression of the aging-related gene RGN in vascular smooth muscle cells, further exacerbating the progression of atherosclerosis. A reduced level of m6A methylation modification could enhance inflammation and expedite cellular aging, thereby contributing to the development of carotid atherosclerosis.

The rise of epitranscriptomics: recent developments and future directions

The epitranscriptomics field has undergone tremendous growth since the discovery that the RNA N6-methyladenosine (m6A) modification is reversible and is distributed throughout the transcriptome. Efforts to map RNA modifications transcriptome-wide and reshape the epitranscriptome in disease settings have facilitated mechanistic understanding and drug discovery in the field. In this review we discuss recent advancements in RNA modification detection methods and consider how these developments can be applied to gain novel insights into the epitranscriptome. We also highlight drug discovery efforts aimed at developing epitranscriptomic therapeutics for cancer and other diseases. Finally, we consider engineering of the epitranscriptome as an emerging direction to investigate RNA modifications and their causal effects on RNA processing at high specificity.

The role of m6A RNA methylation in autoimmune diseases: Novel therapeutic opportunities

N6-methyladenosine (m6A) modifications, as one of the most common forms of internal RNA chemical modifications in eukaryotic cells, have gained increasing attention in recent years. The m6A RNA modifications exert various crucial roles in various biological processes, such as embryonic development, neurogenesis, circadian rhythms, and tumorigenesis. Recent advances have highlighted that m6A RNA modification plays an important role in immune response, especially in the initiation and progression of autoimmune diseases. In this review, we summarized the regulatory mechanisms of m6A methylation and its biological functions in the immune system and mainly focused on recent progress in research on the potential role of m6A RNA methylation in the pathogenesis of autoimmune diseases, thus providing possible biomarkers and potential targets for the prevention and treatment of autoimmune diseases.

Global and single-nucleotide resolution detection of 7-methylguanosine in RNA

RNA modifications, including N-7-methylguanosine (m7G), are pivotal in governing RNA stability and gene expression regulation. The accurate detection of internal m7G modifications is of paramount significance, given recent associations between altered m7G deposition and elevated expression of the methyltransferase METTL1 in various human cancers. The development of robust m7G detection techniques has posed a significant challenge in the field of epitranscriptomics. In this study, we introduce two methodologies for the global and accurate identification of m7G modifications in human RNA. We introduce borohydride reduction sequencing (Bo-Seq), which provides base resolution mapping of m7G modifications. Bo-Seq achieves exceptional performance through the optimization of RNA depurination and scission, involving the strategic use of high concentrations of NaBH4, neutral pH and the addition of 7-methylguanosine monophosphate (m7GMP) during the reducing reaction. Notably, compared to NaBH4-based methods, Bo-Seq enhances the m7G detection performance, and simplifies the detection process, eliminating the necessity for intricate chemical steps and reducing the protocol duration. In addition, we present an antibody-based approach, which enables the assessment of m7G relative levels across RNA molecules and biological samples, however it should be used with caution due to limitations associated with variations in antibody quality between batches. In summary, our novel approaches address the pressing need for reliable and accessible methods to detect RNA m7G methylation in human cells. These advancements hold the potential to catalyse future investigations in the critical field of epitranscriptomics, shedding light on the complex regulatory roles of m7G in gene expression and its implications in cancer biology.

In-Gel Cyanoethylation for Pseudouridines Mass Spectrometry Detection of Bacterial Regulatory RNA

Regulatory RNAs, as well as many RNA families, contain chemically modified nucleotides, including pseudouridines (ψ). To map nucleotide modifications, approaches based on enzymatic digestion of RNA followed by nano liquid chromatography-tandem mass spectrometry (nanoLC-MS/MS) analysis were implemented several years ago. However, detection of ψ by mass spectrometry (MS) is challenging as ψ exhibits the same mass as uridine. Thus, a chemical labeling strategy using acrylonitrile was developed to detect this mass-silent modification. Acrylonitrile reacts specifically to ψ to form 1-cyanoethylpseudouridine (Ceψ), resulting in a mass shift of ψ detectable by MS. Here, a protocol detailing the steps from the purification of RNA by polyacrylamide gel electrophoresis, including in-gel labeling of ψ, to MS data interpretation to map ψ and other modifications is proposed. To demonstrate its efficiency, the protocol was applied to bacterial regulatory RNAs from E. coli: 6S RNA and transfer-messenger RNA (tmRNA, also known as 10Sa RNA). Moreover, ribonuclease P (RNase P) was also mapped using this approach. This method enabled the detection of several ψ at single nucleotide resolution.

A review on the role of RNA methylation in aging-related diseases

Senescence is the underlying mechanism of organism aging and is robustly regulated at the post-transcriptional level. This regulation involves the chemical modifications, of which the RNA methylation is the most common. Recently, a rapidly growing number of studies have demonstrated that methylation is relevant to aging and aging-associated diseases. Owing to the rapid development of detection methods, the understanding on RNA methylation has gone deeper. In this review, we summarize the current understanding on the influence of RNA modification on cellular senescence, with a focus on mRNA methylation in aging-related diseases, and discuss the emerging potential of RNA modification in diagnosis and therapy.

Interconnections between m6A RNA modification, RNA structure, and protein-RNA complex assembly

Protein-RNA complexes exist in many forms within the cell, from stable machines such as the ribosome to transient assemblies like the spliceosome. All protein-RNA assemblies rely on spatially and temporally coordinated interactions between specific proteins and RNAs to achieve a functional form. RNA folding and structure are often critical for successful protein binding and protein-RNA complex formation. RNA modifications change the chemical nature of a given RNA and often alter its folding kinetics. Both these alterations can affect how and if proteins or other RNAs can interact with the modified RNA and assemble into complexes. N6-methyladenosine (m6A) is the most common base modification on mRNAs and regulatory noncoding RNAs and has been shown to impact RNA structure and directly modulate protein-RNA interactions. In this review, focusing on the mechanisms and available quantitative information, we discuss first how the METTL3/14 m6A writer complex is specifically targeted to RNA assisted by protein-RNA and other interactions to enable site-specific and co-transcriptional RNA modification and, once introduced, how the m6A modification affects RNA folding and protein-RNA interactions.

Exploring urinary modified nucleosides as biomarkers for diabetic retinopathy: Development and validation of a ultra performance liquid chromatography-tandem mass spectrometry method

The dynamic modification of RNA plays a crucial role in biological regulation and is strongly linked to human disease development and progression. Notably, modified nucleosides in urine have shown promising potential as early diagnostic biomarkers for various conditions. In this study, we developed and validated a rapid, sensitive, and accurate UPLC-MS/MS method for quantifying eight types of modified nucleosides (N1-methyladenosine (m1A), N6-methyladenosine (m6A), 5-methyluridine (m5U), 5-taurinomethyl-2-thiouridine (τm5s2U), 5-methylcytidine (m5C), 2'-O-methylcytidine (Cm), N1-methylguanosine (m1G), and N7-methylguanosine (m7G) in human urine. Using the method, we measured the urinary concentrations of m1A, m6A, m5U, τm5s2U, m5C, Cm, m1G, and m7G in a total of 21 control individuals and 23 patients diagnosed with diabetic retinopathy (DR). Cm levels showed promise as a diagnostic marker for diabetic retinopathy (DR), with a significant value (P < 0.01) and an AUC of 0.735. Other modified nucleosides also exhibited significant differences within specific subpopulations. As non-proliferative diabetic retinopathy (NPDR) signifies the latent early stage of diabetic retinopathy, we developed a multivariate linear model that integrates patients' sex, age, height, and urinary concentration of modified nucleosides which aims to predict and differentiate between healthy individuals, NPDR patients, and proliferative diabetic retinopathy (PDR) patients. Encouragingly, the model achieved satisfactory accuracy rates: healthy (81%), NPDR (75%), and PDR (80%). Our findings provide valuable insights into the development of an early, cost-effective, and noninvasive diagnostic approach for diabetic retinopathy.

Epitranscriptomic Mass Spectrometry

Every chemical group that is added to any one of the canonical ribonucleotides in a transcript would create a specific RNA modification. Currently, 170+ RNA modifications have been identified. A specific epitranscriptome refers to all the RNA modifications in a given biological system and is considered to play an important role in the regulations of cellular activities. Mass spectrometry-based methods have proven to be the most accurate way to identify RNA modifications and determine the amount of each detectable modification. Relating to the recent development of mapping specific RNA modifications within a transcriptome, the profiling of all RNA modifications can serve as a prescreening tool for mapping and provides support for analyzing the data obtained from mapping. In this chapter, the details for setting up a commonly used mass spectrometry-based method to profile all the RNA modifications in specific epitranscriptomes are described, and the possible options if available are discussed.

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.

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

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

RNA m6A modifications regulate crosstalk between tumor metabolism and immunity

In recent years, m6A modifications in RNA transcripts have arisen as a hot topic in cancer research. Indeed, a number of independent studies have elaborated that the m6A modification impacts the behavior of tumor cells and tumor-infiltrating immune cells, altering tumor cell metabolism along with the differentiation and functional activity of immune cells. This review elaborates on the links between RNA m6A modifications, tumor cell metabolism, and immune cell behavior, discussing this topic from the viewpoint of reciprocal regulation through "RNA m6A-tumor cell metabolism-immune cell behavior" and "RNA m6A-immune cell behavior-tumor cell metabolism" axes. In addition, we discuss the various factors affecting RNA m6A modifications in the tumor microenvironment, particularly the effects of hypoxia associated with cancer cell metabolism along with immune cell-secreted cytokines. Our analysis proposes the conclusion that RNA m6A modifications support widespread interactions between tumor metabolism and tumor immunity. With the current viewpoint that long-term cancer control must tackle cancer cell malignant behavior while strengthening anti-tumor immunity, the recognition of RNA m6A modifications as a key factor provides a new direction for the targeted therapy of tumors. This article is categorized under: RNA Processing > RNA Editing and Modification RNA in Disease and Development > RNA in Disease RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

The catalytic mechanism of the RNA methyltransferase METTL3

The complex of methyltransferase-like proteins 3 and 14 (METTL3-14) is the major enzyme that deposits N6-methyladenosine (m6A) modifications on mRNA in humans. METTL3-14 plays key roles in various biological processes through its methyltransferase (MTase) activity. However, little is known about its substrate recognition and methyl transfer mechanism from its cofactor and methyl donor S-adenosylmethionine (SAM). Here, we study the MTase mechanism of METTL3-14 by a combined experimental and multiscale simulation approach using bisubstrate analogues (BAs), conjugates of a SAM-like moiety connected to the N6-atom of adenosine. Molecular dynamics simulations based on crystal structures of METTL3-14 with BAs suggest that the Y406 side chain of METTL3 is involved in the recruitment of adenosine and release of m6A. A crystal structure with a bisubstrate analogue representing the transition state of methyl transfer shows a direct involvement of the METTL3 side chains E481 and K513 in adenosine binding which is supported by mutational analysis. Quantum mechanics/molecular mechanics (QM/MM) free energy calculations indicate that methyl transfer occurs without prior deprotonation of adenosine-N6. Furthermore, the QM/MM calculations provide further support for the role of electrostatic contributions of E481 and K513 to catalysis. The multidisciplinary approach used here sheds light on the (co)substrate binding mechanism, catalytic step, and (co)product release catalysed by METTL3, and suggests that the latter step is rate-limiting. The atomistic information on the substrate binding and methyl transfer reaction of METTL3 can be useful for understanding the mechanisms of other RNA MTases and for the design of transition state analogues as their inhibitors.

The role of RNA modification in urological cancers: mechanisms and clinical potential

RNA modification is a post-transcriptional level of regulation that is widely distributed in all types of RNAs, including mRNA, tRNA, rRNA, miRNA, and lncRNA, where N6-methyladenine (m6A) is the most abundant mRNA methylation modification. Significant evidence has depicted that m6A modifications are closely related to human diseases, especially cancer, and play pivotal roles in RNA transcription, splicing, stabilization, and translation processes. The most common urological cancers include prostate, bladder, kidney, and testicular cancers, accounting for a certain proportion of human cancers, with an ever-increasing incidence and mortality. The recurrence, systemic metastasis, poor prognosis, and drug resistance of urologic tumors have prompted the identification of new therapeutic targets and mechanisms. Research on m6A modifications may provide new solutions to the current puzzles. In this review, we provide a comprehensive overview of the key roles played by RNA modifications, especially m6A modifications, in urologic cancers, as well as recent research advances in diagnostics and molecularly targeted therapies.

Diversity in Biological Function and Mechanism of the tRNA Methyltransferase Trm10

Transfer ribonucleic acid (tRNA) is the most highly modified RNA species in the cell, and loss of tRNA modifications can lead to growth defects in yeast as well as metabolic, neurological, and mitochondrial disorders in humans. Significant progress has been made toward identifying the enzymes that are responsible for installing diverse modifications in tRNA, revealing a landscape of fascinating biological and mechanistic diversity that remains to be fully explored. Most early discoveries of tRNA modification enzymes were in model systems, where many enzymes were not strictly required for viability, an observation somewhat at odds with the extreme conservation of many of the same enzymes throughout multiple domains of life. Moreover, many tRNA modification enzymes act on more than one type of tRNA substrate, which is not necessarily surprising given the similar overall secondary and tertiary structures of tRNA, yet biochemical characterization has revealed interesting patterns of substrate specificity that can be challenging to rationalize on a molecular level. Questions about how many enzymes efficiently select a precise set of target tRNAs from among a structurally similar pool of molecules persist.The tRNA methyltransferase Trm10 provides an exciting paradigm to study the biological and mechanistic questions surrounding tRNA modifications. Even though the enzyme was originally characterized in Saccharomyces cerevisiae where its deletion causes no detectable phenotype under standard lab conditions, several more recently identified phenotypes provide insight into the requirement for this modification in the overall quality control of the tRNA pool. Studies of Trm10 in yeast also revealed another characteristic feature that has turned out to be a conserved feature of enzymes throughout the Trm10 family tree. We were initially surprised to see that purified S. cerevisiae Trm10 was capable of modifying tRNA substrates that were not detectably modified by the enzyme in vivo in yeast. This pattern has continued to emerge as we and others have studied Trm10 orthologs from Archaea and Eukarya, with enzymes exhibiting in vitro substrate specificities that can differ significantly from in vivo patterns of modification. While this feature complicates efforts to predict substrate specificities of Trm10 enzymes in the absence of appropriate genetic systems, it also provides an exciting opportunity for studying how enzyme activities can be regulated to achieve dynamic patterns of biological tRNA modification, which have been shown to be increasingly important for stress responses and human disease. Finally, the intriguing diversity in target nucleotide modification that has been revealed among Trm10 orthologs is distinctive among known tRNA modifying enzymes and necessitates unusual and likely novel catalytic strategies for methylation that are being revealed by biochemical and structural studies directed toward various family members. These efforts will no doubt yield more surprising discoveries in terms of tRNA modification enzymology.

Roles of RNA m6A modifications in plant-virus interactions

Viral RNAs have been known to contain N6-methyladenosine (m6A) modifications since the 1970s. The function of these modifications remained unknown until the development of genome-wide methods to map m6A residues. Increasing evidence has recently revealed a strong association between m6A modifications and plant viral infection. This highlight introduces advances in the roles of RNA m6A modifications in plant-virus interactions.

Extending the toolbox for RNA biology with SegModTeX: a polymerase-driven method for site-specific and segmental labeling of RNA

RNA performs a wide range of functions regulated by its structure, dynamics, and often post-transcriptional modifications. While NMR is the leading method for understanding RNA structure and dynamics, it is currently limited by the inability to reduce spectral crowding by efficient segmental labeling. Furthermore, because of the challenging nature of RNA chemistry, the tools being developed to introduce site-specific modifications are increasingly complex and laborious. Here we use a previously designed Tgo DNA polymerase mutant to present SegModTeX - a versatile, one-pot, copy-and-paste approach to address these challenges. By precise, stepwise construction of a diverse set of RNA molecules, we demonstrate the technique to be superior to RNA polymerase driven and ligation methods owing to its substantially high yield, fidelity, and selectivity. We also show the technique to be useful for incorporating some fluorescent- and a wide range of other probes, which significantly extends the toolbox of RNA biology in general.

The role of m6A epigenetic modifications in tumor coding and non-coding RNA processing

BACKGROUND

Epigenetic modifications of RNA significantly contribute to the regulatory processes in tumors and have, thus, received considerable attention. The m6A modification, known as N6-methyladenosine, is the predominant epigenetic alteration found in both eukaryotic mRNAs and ncRNAs.

MAIN BODY

m6A methylation modifications are dynamically reversible and are catalyzed, removed, and recognized by the complex of m6A methyltransferase (MTases), m6A demethylase, and m6A methyl recognition proteins (MRPs). Published evidence suggests that dysregulated m6A modification results in abnormal biological behavior of mature mRNA, leading to a variety of abnormal physiological processes, with profound implications for tumor development in particular.

CONCLUSION

Abnormal RNA processing due to dysregulation of m6A modification plays an important role in tumor pathogenesis and potential mechanisms of action. In this review, we comprehensively explored the mechanisms by which m6A modification regulates mRNA and ncRNA processing, focusing on their roles in tumors, and aiming to understand the important regulatory function of m6A modification, a key RNA epigenetic modification, in tumor cells, with a view to providing theoretical support for tumor diagnosis and treatment. Video Abstract.

DPred_3S: identifying dihydrouridine (D) modification on three species epitranscriptome based on multiple sequence-derived features

Introduction: Dihydrouridine (D) is a conserved modification of tRNA among all three life domains. D modification enhances the flexibility of a single nucleotide base in the spatial structure and is disease- and evolution-associated. Recent studies have also suggested the presence of dihydrouridine on mRNA. Methods: To identify D in epitranscriptome, we provided a prediction framework named "DPred_3S" based on the machine learning approach for three species D epitranscriptome, which used epitranscriptome sequencing data as training data for the first time. Results: The optimal features were evaluated by the F-score and integration of different features; our model achieved area under the receiver operating characteristic curve (AUROC) scores 0.955, 0.946, and 0.905 for Saccharomyces cerevisiae, Escherichia coli, and Schizosaccharomyces pombe, respectively. The performances of different machine learning algorithms were also compared in this study. Discussion: The high performances of our model suggest the D sites can be distinguished based on their surrounding sequence, but the lower performance of cross-species prediction may be limited by technique preferences.

Insights into the Mechanism of Installation of 5-Carboxymethylaminomethyl Uridine Hypermodification by tRNA-Modifying Enzymes MnmE and MnmG

The evolutionarily conserved bacterial proteins MnmE and MnmG (and their homologues in Eukarya) install a 5-carboxymethylaminomethyl (cmnm5) or a 5-taurinomethyl (τm5) group onto wobble uridines of several tRNA species. The Escherichia coli MnmE binds guanosine-5'-triphosphate (GTP) and methylenetetrahydrofolate (CH2THF), while MnmG binds flavin adenine dinucleotide (FAD) and a reduced nicotinamide adenine dinucleotide (NADH). Together with glycine, MnmEG catalyzes the installation of cmnm5 in a reaction that also requires hydrolysis of GTP. In this letter, we investigated key steps of the MnmEG reaction using a combination of biochemical techniques. We show multiple lines of evidence supporting flavin-iminium FADH[N5═CH2]+ as a central intermediate in the MnmEG reaction. Using a synthetic FADH[N5═CD2]+ analogue, the intermediacy of the FAD in the transfer of the methylene group from CH2THF to the C5 position of U34 was unambiguously demonstrated. Further, MnmEG reactions containing the deuterated flavin-iminium intermediate and alternate nucleophiles such as taurine and ammonia also led to the formation of the anticipated U34-modified tRNAs, showing FAD[N5═CH2]+ as the universal intermediate for all MnmEG homologues. Additionally, an RNA-protein complex stable to urea-denaturing polyacrylamide gel electrophoresis was identified. Studies involving a series of nuclease (RNase T1) and protease (trypsin) digestions along with reverse transcription experiments suggest that the complex may be noncovalent. While the conserved MnmG cysteine C47 and C277 mutant variants were shown to reduce FAD, they were unable to promote the modified tRNA formation. Overall, this study provides critical insights into the biochemical mechanism underlying tRNA modification by the MnmEG.