A Platform for Discovery and Quantification of Modified Ribonucleosides in RNA: Application to Stress-Induced Reprogramming of tRNA Modifications

tRNA modification enzyme-dependent redox homeostasis regulates synapse formation and memory

Post-transcriptional modification of RNA regulates gene expression at multiple levels. ALKBH8 is a tRNA-modifying enzyme that methylates wobble uridines in a subset of tRNAs to modulate translation. Through methylation of tRNA-selenocysteine, ALKBH8 promotes selenoprotein synthesis and regulates redox homeostasis. Pathogenic variants in ALKBH8 have been linked to intellectual disability disorders in the human population, but the role of ALKBH8 in the nervous system is unknown. Through in vivo studies in Drosophila, we show that ALKBH8 controls oxidative stress in the brain to restrain synaptic growth and support learning and memory. ALKBH8 null animals lack wobble uridine methylation and exhibit reduced protein synthesis in the nervous system, including a specific decrease in selenoprotein levels. Either loss of ALKBH8 or independent disruption of selenoprotein synthesis results in ectopic synapse formation. Genetic expression of antioxidant enzymes fully suppresses synaptic overgrowth in ALKBH8 null animals, confirming oxidative stress as the underlying cause of dysregulation. ALKBH8 null animals also exhibit associative memory impairments that are reversed by pharmacological antioxidant treatment. Together, these findings demonstrate the critical role of tRNA wobble uridine modification in redox homeostasis in the developing nervous system and reveal antioxidants as a potential therapy for ALKBH8-associated intellectual disability.

Expanding the Available RNA Labeling Toolbox With CutA Nucleotidyltransferase for Efficient Transcript Labeling with Purine and Pyrimidine Nucleotide Analogs

RNA labeling is an invaluable tool for investigation of the function and localization of nucleic acids. Labels are commonly incorporated into 3' end of RNA and the primary enzyme used for this purpose is RNA poly(A) polymerase (PAP), which belongs to the class of terminal nucleotidyltransferases (NTases). However, PAP preferentially adds ATP analogs, thus limiting the number of available substrates. Here, we report the use of another NTase, CutA from the fungus Thielavia terrestris. Using this enzyme, we were able to incorporate into the 3' end of RNA not only purine analogs, but also pyrimidine analogs. We engaged strain-promoted azide-alkyl cycloaddition (SPAAC) to obtain fluorescently labeled or biotinylated transcripts from RNAs extended with azide analogs by CutA. Importantly, modified transcripts retained their biological properties. Furthermore, fluorescently labeled mRNAs were suitable for visualization in cultured mammalian cells. Finally, we demonstrate that either affinity studies or molecular dynamic (MD) simulations allow for rapid screening of NTase substrates, what opens up new avenues in the search for the optimal substrates for this class of enzymes.

Esterification of Cyclic N6-Threonylcarbamoyladenosine During RNA Sample Preparation

The continuous deciphering of crucial biological roles of RNA modifications and their involvement in various pathological conditions, together with their key roles in the use of RNA-based therapeutics, has reignited interest in studying the occurrence and identity of non-canonical ribonucleoside structures during the past years. Discovery and structural elucidation of new modified structures is usually achieved by combination of liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) at the nucleoside level and stable isotope labeling experiments. This approach, however, has its pitfalls as demonstrated in the course of the present study: we structurally elucidated a new nucleoside structure that showed significant similarities to the family of (c)t6A modifications and was initially considered a genuine modification, but subsequently turned out to be an in vitro formed glycerol ester of t6A. This artifact is generated from ct6A during RNA hydrolysis upon addition of enzymes stored in glycerol containing buffers in a mildly alkaline milieu, and was moreover shown to undergo an intramolecular transesterification reaction. Our results demand for extra caution, not only in the discovery of new RNA modifications, but also with regard to the quantification of known modified structures, in particular chemically labile modifications, such as ct6A, that might suffer from exposure to putatively harmless reagents during the diverse steps of sample preparation.

PT-seq: A method for metagenomic analysis of phosphorothioate epigenetics in complex microbial communities

Among dozens of known epigenetic marks, naturally occurring phosphorothioate (PT) DNA modifications are unique in replacing a non-bridging phosphate oxygen with redox-active sulfur and function in prokaryotic restriction-modification and transcriptional regulation. Interest in PTs has grown due to the widespread distribution of the dnd, ssp, and brx genes among bacteria and archaea, as well as the discovery of PTs in 5-10% of gut microbes. Efforts to map PTs in complex microbiomes using existing next-generation and direct sequencing technologies have failed due to poor sensitivity. Here we developed PT-seq as a high-sensitivity method to quantitatively map PTs across genomes and metagenomically identify PT-containing microbes in complex genomic mixtures. Like other methods for mapping PTs in individual genomes, PT-seq exploits targeted DNA strand cleavage at PTs by iodine, followed by sequencing library construction using ligation or template switching approaches. However, PT-specific sequencing reads are dramatically increased by adding steps to heat denature the DNA, block pre-existing 3'-ends, fragment DNA after T-tailing, and enrich iodine-induced breaks using biotin-labeling and streptavidin beads capture. Iterative optimization of the sensitivity and specificity of PT-seq is demonstrated with individual bacteria and human fecal DNA.

Proteolytic cleavage and inactivation of the TRMT1 tRNA modification enzyme by SARS-CoV-2 main protease

Nonstructural protein 5 (Nsp5) is the main protease of SARS-CoV-2 that cleaves viral polyproteins into individual polypeptides necessary for viral replication. Here, we show that Nsp5 binds and cleaves human tRNA methyltransferase 1 (TRMT1), a host enzyme required for a prevalent post-transcriptional modification in tRNAs. Human cells infected with SARS-CoV-2 exhibit a decrease in TRMT1 protein levels and TRMT1-catalyzed tRNA modifications, consistent with TRMT1 cleavage and inactivation by Nsp5. Nsp5 cleaves TRMT1 at a specific position that matches the consensus sequence of SARS-CoV-2 polyprotein cleavage sites, and a single mutation within the sequence inhibits Nsp5-dependent proteolysis of TRMT1. The TRMT1 cleavage fragments exhibit altered RNA binding activity and are unable to rescue tRNA modification in TRMT1-deficient human cells. Compared to wild-type human cells, TRMT1-deficient human cells infected with SARS-CoV-2 exhibit reduced levels of intracellular viral RNA. These findings provide evidence that Nsp5-dependent cleavage of TRMT1 and perturbation of tRNA modification patterns contribute to the cellular pathogenesis of SARS-CoV-2 infection.

Abundance of Modifications in Mature miRNAs Revealed by LC-MS/MS Method Coupled with a Two-Step Hybridization Purification Strategy

Accurate detection of endogenous miRNA modifications, such as N6-methyladenosine (m6A), 7-methylguanosine (m7G), and 5-methylcytidine (m5C), poses significant challenges, resulting in considerable uncertainty regarding their presence in mature miRNAs. In this study, we demonstrate for the first time that liquid chromatography coupled with a tandem mass spectrometry (LC-MS/MS) nucleoside analysis method is a practical tool for quantitatively analyzing human miRNA modifications. The newly designed liquid-solid two-step hybridization (LSTH) strategy enhances specificity for miRNA purification, while LC-MS/MS offers robust capability in recognizing modifications and sufficient sensitivity with detection limits ranging from attomoles to low femtomoles. Therefore, it provides a more reliable approach compared to existing techniques for revealing modifications in endogenous miRNAs. With this approach, we characterized m6A, m7G, and m5C modifications in miR-21-5p, Let-7a/e-5p, and miR-10a-5p isolated from cultured cells and observed unexpectedly low abundance (<1% at each site) of these modifications.

Proteolytic cleavage and inactivation of the TRMT1 tRNA modification enzyme by SARS-CoV-2 main protease

Nonstructural protein 5 (Nsp5) is the main protease of SARS-CoV-2 that cleaves viral polyproteins into individual polypeptides necessary for viral replication. Here, we show that Nsp5 binds and cleaves human tRNA methyltransferase 1 (TRMT1), a host enzyme required for a prevalent post-transcriptional modification in tRNAs. Human cells infected with SARS-CoV-2 exhibit a decrease in TRMT1 protein levels and TRMT1-catalyzed tRNA modifications, consistent with TRMT1 cleavage and inactivation by Nsp5. Nsp5 cleaves TRMT1 at a specific position that matches the consensus sequence of SARS-CoV-2 polyprotein cleavage sites, and a single mutation within the sequence inhibits Nsp5-dependent proteolysis of TRMT1. The TRMT1 cleavage fragments exhibit altered RNA binding activity and are unable to rescue tRNA modification in TRMT1-deficient human cells. Compared to wildtype human cells, TRMT1-deficient human cells infected with SARS-CoV-2 exhibit reduced levels of intracellular viral RNA. These findings provide evidence that Nsp5-dependent cleavage of TRMT1 and perturbation of tRNA modification patterns contribute to the cellular pathogenesis of SARS-CoV-2 infection.

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.

Pitfalls in RNA Modification Quantification Using Nucleoside Mass Spectrometry

In recent years, there has been a high interest in researching RNA modifications, as they are involved in many cellular processes and in human diseases. A substantial set of enzymes within the cell, called RNA writers, place RNA modifications selectively and site-specifically. Another set of enzymes, called readers, recognize these modifications which guide the fate of the modified RNA. Although RNA is a transient molecule and RNA modification could be removed by RNA degradation, a subclass of enzymes, called RNA erasers, remove RNA modifications selectively and site-specifically to alter the characteristics of the RNA. The detection of RNA modifications can be done by various methods including second and next generation sequencing but also mass spectrometry. An approach capable of both qualitative and quantitative RNA modification analysis is liquid chromatography coupled to mass spectrometry of enzymatic hydrolysates of RNA into nucleosides. However, for successful detection and quantification, various factors must be considered to avoid biased identification and inaccurate quantification. In this Account, we identify three classes of errors that may distort the analysis. These classes comprise (I) errors related to chemical instabilities, (II) errors revolving around enzymatic hydrolysis to nucleosides, and (III) errors arising from issues with chromatographic separation and/or subsequent mass spectrometric analysis.A prominent example for class 1 is Dimroth rearrangement of m1A to m6A, but class 1 also comprises hydrolytic reactions and reactions with buffer components. Here, we also present the conversion of m3C to m3U under mild alkaline conditions and propose a practical solution to overcome these instabilities. Class 2 errors-such as contaminations in hydrolysis reagents or nuclease specificities-have led to erroneous discoveries of nucleosides in the past and possess the potential for misquantification of nucleosides. Impurities in the samples may also lead to class 3 errors: For instance, issues with chromatographic separation may arise from residual organic solvents, and salt adducts may hamper mass spectrometric quantification. This Account aims to highlight various errors connected to mass spectrometry analysis of nucleosides and presents solutions for how to overcome or circumnavigate those issues. Therefore, the authors anticipate that many scientists, but especially those who plan on doing nucleoside mass spectrometry, will benefit from the collection of data presented in this Account as a raised awareness, toward the variety of potential pitfalls, may further enhance the quality of data.

tRNA modification enzyme-dependent redox homeostasis regulates synapse formation and memory

Post-transcriptional modification of RNA regulates gene expression at multiple levels. ALKBH8 is a tRNA modifying enzyme that methylates wobble uridines in specific tRNAs to modulate translation. Through methylation of tRNA-selenocysteine, ALKBH8 promotes selenoprotein synthesis and regulates redox homeostasis. Pathogenic variants in ALKBH8 have been linked to intellectual disability disorders in the human population, but the role of ALKBH8 in the nervous system is unknown. Through in vivo studies in Drosophila, we show that ALKBH8 controls oxidative stress in the brain to restrain synaptic growth and support learning and memory. ALKBH8 null animals lack wobble uridine methylation and exhibit a global reduction in protein synthesis, including a specific decrease in selenoprotein levels. Loss of ALKBH8 or independent disruption of selenoprotein synthesis results in ectopic synapse formation. Genetic expression of antioxidant enzymes fully suppresses synaptic overgrowth in ALKBH8 null animals, confirming oxidative stress as the underlying cause of dysregulation. ALKBH8 animals also exhibit associative learning and memory impairments that are reversed by pharmacological antioxidant treatment. Together, these findings demonstrate the critical role of tRNA modification in redox homeostasis in the nervous system and reveal antioxidants as a potential therapy for ALKBH8-associated intellectual disability.

Expanded tRNA methyltransferase family member TRMT9B regulates synaptic growth and function

Nervous system function rests on the formation of functional synapses between neurons. We have identified TRMT9B as a new regulator of synapse formation and function in Drosophila. TRMT9B has been studied for its role as a tumor suppressor and is one of two metazoan homologs of yeast tRNA methyltransferase 9 (Trm9), which methylates tRNA wobble uridines. Whereas Trm9 homolog ALKBH8 is ubiquitously expressed, TRMT9B is enriched in the nervous system. However, in the absence of animal models, TRMT9B's role in the nervous system has remained unstudied. Here, we generate null alleles of TRMT9B and find it acts postsynaptically to regulate synaptogenesis and promote neurotransmission. Through liquid chromatography-mass spectrometry, we find that ALKBH8 catalyzes canonical tRNA wobble uridine methylation, raising the question of whether TRMT9B is a methyltransferase. Structural modeling studies suggest TRMT9B retains methyltransferase function and, in vivo, disruption of key methyltransferase residues blocks TRMT9B's ability to rescue synaptic overgrowth, but not neurotransmitter release. These findings reveal distinct roles for TRMT9B in the nervous system and highlight the significance of tRNA methyltransferase family diversification in metazoans.

Small RNA modifications: regulatory molecules and potential applications

Small RNAs (also referred to as small noncoding RNAs, sncRNA) are defined as polymeric ribonucleic acid molecules that are less than 200 nucleotides in length and serve a variety of essential functions within cells. Small RNA species include microRNA (miRNA), PIWI-interacting RNA (piRNA), small interfering RNA (siRNA), tRNA-derived small RNA (tsRNA), etc. Current evidence suggest that small RNAs can also have diverse modifications to their nucleotide composition that affect their stability as well as their capacity for nuclear export, and these modifications are relevant to their capacity to drive molecular signaling processes relevant to biogenesis, cell proliferation and differentiation. In this review, we highlight the molecular characteristics and cellular functions of small RNA and their modifications, as well as current techniques for their reliable detection. We also discuss how small RNA modifications may be relevant to the clinical applications for the diagnosis and treatment of human health conditions such as cancer.

Identification of 2-Methylthio-methylenethio-N6 -(cis-4-hydroxyisopentenyl)adenosine (msms2 io6 A37 ) as a Novel Modification at Adenosine 37 of tRNAs from Salmonella typhimurium

Post-transcriptional modifications of tRNA nucleotide are important determinants in folding, structure and function. We have successfully identified and characterized a new modified base named 2-methylthio-methylenethio-N⁶ -(cis-4-hydroxyisopentenyl)adenosine, which is present at position 37 in some tRNAs. We also showed that this new modified adenosine is derived from the known 2-methylthio-methylenethio-N⁶ -(isopentenyl)adenosine nucleoside by a catalytic cycle of the tRNA-diiron monooxygenase, MiaE, present in Salmonella typhimurium.

RNA modifications: an overview of select web-based tools

The field of epitranscriptomics has expanded dramatically in recent years, both in the number of identified RNA modifications and the number of researchers studying them. As knowledge of post-transcriptional modifications continues to expand, numerous new methods have been developed to detect these modifications. Additionally, modifications are being extended to therapeutic settings, such as with recent mRNA vaccines. With this increase in knowledge and use, the community is recognizing the necessity for user-friendly databases to (i) store information from both high- and low-throughput studies and (ii) provide prediction software on how RNA modifications contribute to RNA function and disease. This mini-review highlights select RNA modification databases and their key attributes with the aim of providing a resource to researchers in the field of epitranscriptomics.

Dysfunctional tRNA reprogramming and codon-biased translation in cancer

Many cancers hijack translation to increase the synthesis of tumor-driving proteins, the messenger mRNAs of which have specific codon usage patterns. Termed 'codon-biased translation' and originally identified in stress response regulation, this mechanism is supported by diverse studies demonstrating how the 50 RNA modifications of the epitranscriptome, specific tRNAs, and codon-biased mRNAs are used by oncogenic programs to promote proliferation and chemoresistance. The epitranscriptome writers METTL1-WDR4, Elongator complex protein (ELP)1-6, CTU1-2, and ALKBH8-TRM112 illustrate the principal mechanism of codon-biased translation, with gene amplifications, increased RNA modifications, and enhanced tRNA stability promoting cancer proliferation. Furthermore, systems-level analyses of 34 tRNA writers and 493 tRNA genes highlight the theme of tRNA epitranscriptome dysregulation in many cancers and identify candidate tRNA writers, tRNA modifications, and tRNA molecules as drivers of pathological codon-biased translation.

Higher-Energy Collisional Dissociation Mass Spectral Networks for the Rapid, Semi-automated Characterization of Known and Unknown Ribonucleoside Modifications

Higher-energy collisional dissociation (HCD) of modified ribonucleosides generates characteristic and highly reproducible nucleoside-specific tandem mass spectra (MS/MS). Here, we demonstrate the capability of HCD spectra in combination with spectral matching for the semi-automated characterization of ribonucleosides. This process involved the generation of an HCD spectral library and the establishment of a mass spectral network for rapid detection with high sensitivity and specificity in a retention time-independent fashion. Systematic spectral matching analysis of the MS/MS spectra of tRNA hydrolysates from different organisms has helped us to uncover evidence for the existence of novel ribonucleoside modifications such as s²Cm and OHyW-14. Such an untargeted label-free approach has the potential to be integrated with other methods, including those that use isotope labeling, to simplify the characterization of unknown modified ribonucleosides. These findings suggest the compilation of a universal spectral network, for the characterization of known and unknown ribonucleosides, could accelerate discoveries in the epitranscriptome.

The absence of the queuosine tRNA modification leads to pleiotropic phenotypes revealing perturbations of metal and oxidative stress homeostasis in Escherichia coli K12

Queuosine (Q) is a conserved hypermodification of the wobble base of tRNA containing GUN anticodons but the physiological consequences of Q deficiency are poorly understood in bacteria. This work combines transcriptomic, proteomic and physiological studies to characterize a Q-deficient Escherichia coli K12 MG1655 mutant. The absence of Q led to an increased resistance to nickel and cobalt, and to an increased sensitivity to cadmium, compared to the wild-type (WT) strain. Transcriptomic analysis of the WT and Q-deficient strains, grown in the presence and absence of nickel, revealed that the nickel transporter genes (nikABCDE) are downregulated in the Q- mutant, even when nickel is not added. This mutant is therefore primed to resist to high nickel levels. Downstream analysis of the transcriptomic data suggested that the absence of Q triggers an atypical oxidative stress response, confirmed by the detection of slightly elevated reactive oxygen species (ROS) levels in the mutant, increased sensitivity to hydrogen peroxide and paraquat, and a subtle growth phenotype in a strain prone to accumulation of ROS.

Data-Independent Acquisition for the Detection of Mononucleoside RNA Modifications by Mass Spectrometry

RNA is dynamically modified in cells by a plethora of chemical moieties to modulate molecular functions and processes. Over 140 modifications have been identified across species and RNA types, with the highest density and diversity of modifications found in tRNA (tRNA). The methods used to identify and quantify these modifications have developed over recent years and continue to advance, primarily in the fields of next-generation sequencing (NGS) and mass spectrometry (MS). Most current NGS methods are limited to antibody-recognized or chemically derivatized modifications and have limitations in identifying multiple modifications simultaneously. Mass spectrometry can overcome both of these issues, accurately identifying a large number of modifications in a single run. Here, we present advances in MS data acquisition for the purpose of RNA modification identification and quantitation. Using this approach, we identified multiple tRNA wobble position modifications in Arabidopsis thaliana that are upregulated in salt-stressed growth conditions and may stabilize translation of salt stress induced proteins. This work presents improvements in methods for studying RNA modifications and introduces a possible regulatory role of wobble position modifications in A. thaliana translation.

Challenges with Simulating Modified RNA: Insights into Role and Reciprocity of Experimental and Computational Approaches

RNA is critical to a broad spectrum of biological and viral processes. This functional diversity is a result of their dynamic nature; the variety of three-dimensional structures that they can fold into; and a host of post-transcriptional chemical modifications. While there are many experimental techniques to study the structural dynamics of biomolecules, molecular dynamics simulations (MDS) play a significant role in complementing experimental data and providing mechanistic insights. The accuracy of the results obtained from MDS is determined by the underlying physical models i.e., the force-fields, that steer the simulations. Though RNA force-fields have received a lot of attention in the last decade, they still lag compared to their protein counterparts. The chemical diversity imparted by the RNA modifications adds another layer of complexity to an already challenging problem. Insight into the effect of RNA modifications upon RNA folding and dynamics is lacking due to the insufficiency or absence of relevant experimental data. This review provides an overview of the state of MDS of modified RNA, focusing on the challenges in parameterization of RNA modifications as well as insights into relevant reference experiments necessary for their calibration.

Epitranscriptomic Reprogramming Is Required to Prevent Stress and Damage from Acetaminophen

Epitranscriptomic marks, in the form of enzyme catalyzed RNA modifications, play important gene regulatory roles in response to environmental and physiological conditions. However, little is known with respect to how acute toxic doses of pharmaceuticals influence the epitranscriptome. Here we define how acetaminophen (APAP) induces epitranscriptomic reprogramming and how the writer Alkylation Repair Homolog 8 (Alkbh8) plays a key gene regulatory role in the response. Alkbh8 modifies tRNA selenocysteine (tRNA^Sec) to translationally regulate the production of glutathione peroxidases (Gpx’s) and other selenoproteins, with Gpx enzymes known to play protective roles during APAP toxicity. We demonstrate that APAP increases toxicity and markers of damage, and decreases selenoprotein levels in Alkbh8 deficient mouse livers, when compared to wildtype. APAP also promotes large scale reprogramming of many RNA marks comprising the liver tRNA epitranscriptome including: 5-methoxycarbonylmethyluridine (mcm⁵U), isopentenyladenosine (i⁶A), pseudouridine (Ψ), and 1-methyladenosine (m¹A) modifications linked to tRNA^Sec and many other tRNA’s. Alkbh8 deficiency also leads to wide-spread epitranscriptomic dysregulation in response to APAP, demonstrating that a single writer defect can promote downstream changes to a large spectrum of RNA modifications. Our study highlights the importance of RNA modifications and translational responses to APAP, identifies writers as key modulators of stress responses in vivo and supports the idea that the epitranscriptome may play important roles in responses to pharmaceuticals.

High-performance nano-flow liquid chromatography column combined with high- and low-collision energy data-independent acquisition enables targeted and discovery identification of modified ribonucleotides by mass spectrometry

Over 170 post-transcriptional RNA modifications have been described and are common in all kingdoms of life. These modifications range from methylation to complex chemical structures, with methylation being the most abundant. RNA modifications play a key role in RNA folding and function and their dysregulation in humans has been linked to several diseases such as cancer, metabolic diseases or neurological disorder. Nowadays, liquid chromatography-tandem mass spectrometry is considered the gold standard method for the identification and quantification of these modifications due to its sensitivity and accuracy. However, the analysis of modified ribonucleosides by mass spectrometry is complex due to the presence of positional isomers. In this scenario, optimal separation of these compounds by highly sensitive liquid chromatography combined with the generation of high-information spectra is critical to unequivocally identify them, especially in high-complex mixtures. Here we present an analytical method that comprises a new type of mixed-mode nano-flow liquid chromatography column combined with high- and low-collision energy data-independent mass spectrometric acquisition for the identification and quantitation of modified ribonucleosides. The method produces content-rich spectra and combines targeted and screening capabilities thus enabling the identification of a variety of modified nucleosides in biological matrices by single-shot liquid chromatographic analysis coupled to mass spectrometry.

Strategies to Avoid Artifacts in Mass Spectrometry-Based Epitranscriptome Analyses

In this report, we perform structure validation of recently reported RNA phosphorothioate (PT) modifications, a new set of epitranscriptome marks found in bacteria and eukaryotes including humans. By comparing synthetic PT-containing diribonucleotides with native species in RNA hydrolysates by high-resolution mass spectrometry (MS), metabolic stable isotope labeling, and PT-specific iodine-desulfurization, we disprove the existence of PTs in RNA from E. coli, S. cerevisiae, human cell lines, and mouse brain. Furthermore, we discuss how an MS artifact led to the initial misidentification of 2'-O-methylated diribonucleotides as RNA phosphorothioates. To aid structure validation of new nucleic acid modifications, we present a detailed guideline for MS analysis of RNA hydrolysates, emphasizing how the chosen RNA hydrolysis protocol can be a decisive factor in discovering and quantifying RNA modifications in biological samples.

Locating chemical modifications in RNA sequences through ribonucleases and LC-MS based analysis

Knowledge of the structural information is essential for understanding the functional details of modified RNA. Cellular non-coding RNA such as rRNA, tRNA and even viral RNAs contain a number of post-transcriptional modifications with varied degree of diversity and density. In this chapter, we discuss the use of a combination of biochemical and analytical tools such as ribonucleases and liquid chromatography coupled with mass spectrometry approaches for characterization of modified RNA. We present the protocols and alternate strategies for obtaining confident modified sequence information to facilitate the understanding of function.

Instrumental analysis of RNA modifications

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

Broad-range RNA modification analysis of complex biological samples using rapid C18-UPLC-MS

Post-transcriptional RNA modifications play an important role in cellular metabolism with homoeostatic disturbances manifesting as a wide repertoire of phenotypes, reduced stress tolerance and translational perturbation, developmental defects, and diseases, such as type II diabetes, leukaemia, and carcinomas. Hence, there has been an intense effort to develop various methods for investigating RNA modifications and their roles in various organisms, including sequencing-based approaches and, more frequently, liquid chromatography–mass spectrometry (LC-MS)-based methods. Although LC-MS offers numerous advantages, such as being highly sensitive and quantitative over a broad detection range, some stationary phase chemistries struggle to resolve positional isomers. Furthermore, the demand for detailed analyses of complex biological samples often necessitates long separation times, hampering sample-to-sample turnover and making multisample analyses time consuming. To overcome this limitation, we have developed an ultra-performance LC-MS (UPLC-MS) method that uses an octadecyl carbon chain (C18)-bonded silica matrix for the efficient separation of 50 modified ribonucleosides, including positional isomers, in a single 9-min sample-to-sample run. To validate the performance and versatility of our method, we analysed tRNA modification patterns of representative microorganisms from each domain of life, namely Archaea (Methanosarcina acetivorans), Bacteria (Pseudomonas syringae), and Eukarya (Saccharomyces cerevisiae). Additionally, our method is flexible and readily applicable for detection and relative quantification using stable isotope labelling and targeted approaches like multiple reaction monitoring (MRM). In conclusion, this method represents a fast and robust tool for broad-range exploration and quantification of ribonucleosides, facilitating future homoeostasis studies of RNA modification in complex biological samples.

Bacterial tRNA 2'-O-methylation is dynamically regulated under stress conditions and modulates innate immune response

RNA modifications are a well-recognized way of gene expression regulation at the post-transcriptional level. Despite the importance of this level of regulation, current knowledge on modulation of tRNA modification status in response to stress conditions is far from being complete. While it is widely accepted that tRNA modifications are rather dynamic, such variations are mostly assessed in terms of total tRNA, with only a few instances where changes could be traced to single isoacceptor species. Using Escherichia coli as a model system, we explored stress-induced modulation of 2'-O-methylations in tRNAs by RiboMethSeq. This analysis and orthogonal analytical measurements by LC-MS show substantial, but not uniform, increase of the Gm18 level in selected tRNAs under mild bacteriostatic antibiotic stress, while other Nm modifications remain relatively constant. The absence of Gm18 modification in tRNAs leads to moderate alterations in E. coli mRNA transcriptome, but does not affect polysomal association of mRNAs. Interestingly, the subset of motility/chemiotaxis genes is significantly overexpressed in ΔTrmH mutant, this corroborates with increased swarming motility of the mutant strain. The stress-induced increase of tRNA Gm18 level, in turn, reduced immunostimulation properties of bacterial tRNAs, which is concordant with the previous observation that Gm18 is a suppressor of Toll-like receptor 7 (TLR7)-mediated interferon release. This documents an effect of stress induced modulation of tRNA modification that acts outside protein translation.

Eukaryotic life without tQCUG: the role of Elongator-dependent tRNA modifications in Dictyostelium discoideum

In the Elongator-dependent modification pathway, chemical modifications are introduced at the wobble uridines at position 34 in transfer RNAs (tRNAs), which serve to optimize codon translation rates. Here, we show that this three-step modification pathway exists in Dictyostelium discoideum, model of the evolutionary superfamily Amoebozoa. Not only are previously established modifications observable by mass spectrometry in strains with the most conserved genes of each step deleted, but also additional modifications are detected, indicating a certain plasticity of the pathway in the amoeba. Unlike described for yeast, D. discoideum allows for an unconditional deletion of the single tQ^CUG gene, as long as the Elongator-dependent modification pathway is intact. In gene deletion strains of the modification pathway, protein amounts are significantly reduced as shown by flow cytometry and Western blotting, using strains expressing different glutamine leader constructs fused to GFP. Most dramatic are these effects, when the tQ^CUG gene is deleted, or Elp3, the catalytic component of the Elongator complex is missing. In addition, Elp3 is the most strongly conserved protein of the modification pathway, as our phylogenetic analysis reveals. The implications of this observation are discussed with respect to the evolutionary age of the components acting in the Elongator-dependent modification pathway.

METTL6 is a tRNA m3C methyltransferase that regulates pluripotency and tumor cell growth

Recently, covalent modifications of RNA, such as methylation, have emerged as key regulators of all aspects of RNA biology and have been implicated in numerous diseases, for instance, cancer. Here, we undertook a combination of in vitro and in vivo screens to test 78 potential methyltransferases for their roles in hepatocellular carcinoma (HCC) cell proliferation. We identified methyltransferase-like protein 6 (METTL6) as a crucial regulator of tumor cell growth. We show that METTL6 is a bona fide transfer RNA (tRNA) methyltransferase, catalyzing the formation of 3-methylcytidine at C32 of specific serine tRNA isoacceptors. Deletion of Mettl6 in mouse stem cells results in changes in ribosome occupancy and RNA levels, as well as impaired pluripotency. In mice, Mettl6 knockout results in reduced energy expenditure. We reveal a previously unknown pathway in the maintenance of translation efficiency with a role in maintaining stem cell self-renewal, as well as impacting tumor cell growth profoundly.

RNA Phosphorothioate Modification in Prokaryotes and Eukaryotes

RNA modifications play important roles in RNA structures and regulation of gene expression and translation. We report the first RNA modification on the phosphate, the RNA phosphorothioate (PS) modification, discovered in both prokaryotes and eukaryotes. The PS modification is also first reported on nucleic acids of eukaryotes. The GpsG modification exists in the Rp configuration and was quantified with liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). By knocking out the DndA gene in E. coli, we show the Dnd clusters that regulate DNA PS modification may also play roles in RNA PS modification. We also show that the GpsG modification locates on rRNA in E. coli, L. lactis, and HeLa cells, and it is not detected in rRNA-depleted total RNAs from these cells.

Characterization of Post-Transcriptional RNA Modifications by Sheathless Capillary Electrophoresis-High Resolution Mass Spectrometry

Over the past decade there has been a growing interest in RNA modification analysis. High performance liquid chromatography-tandem mass spectrometry coupling (HPLC-MS/MS) is classically used to characterize post-transcriptional modifications of ribonucleic acids (RNAs). Here we propose a novel and simple workflow based on capillary zone electrophoresis-tandem mass spectrometry (CE-MS/MS), in positive mode, to characterize RNA modifications at nucleoside and oligonucleotide levels. By first totally digesting the purified RNA, prior to CE-MS/MS analysis, we were able to identify the nucleoside modifications. Then, using a bottom-up approach, sequencing of the RNAs and mapping of the modifications were performed. Sequence coverages from 68% to 97% were obtained for four tRNAs. Furthermore, unambiguous identification and mapping of several modifications were achieved.

The rRNA m6A methyltransferase METTL5 is involved in pluripotency and developmental programs

Covalent chemical modifications of cellular RNAs directly impact all biological processes. However, our mechanistic understanding of the enzymes catalyzing these modifications, their substrates and biological functions, remains vague. Amongst RNA modifications N⁶-methyladenosine (m⁶A) is widespread and found in messenger (mRNA), ribosomal (rRNA), and noncoding RNAs. Here, we undertook a systematic screen to uncover new RNA methyltransferases. We demonstrate that the methyltransferase-like 5 (METTL5) protein catalyzes m⁶A in 18S rRNA at position A₁₈₃₂ We report that absence of Mettl5 in mouse embryonic stem cells (mESCs) results in a decrease in global translation rate, spontaneous loss of pluripotency, and compromised differentiation potential. METTL5-deficient mice are born at non-Mendelian rates and develop morphological and behavioral abnormalities. Importantly, mice lacking METTL5 recapitulate symptoms of patients with DNA variants in METTL5, thereby providing a new mouse disease model. Overall, our biochemical, molecular, and in vivo characterization highlights the importance of m⁶A in rRNA in stemness, differentiation, development, and diseases.

A computational platform for high-throughput analysis of RNA sequences and modifications by mass spectrometry

The field of epitranscriptomics continues to reveal how post-transcriptional modification of RNA affects a wide variety of biological phenomena. A pivotal challenge in this area is the identification of modified RNA residues within their sequence contexts. Mass spectrometry (MS) offers a comprehensive solution by using analogous approaches to shotgun proteomics. However, software support for the analysis of RNA MS data is inadequate at present and does not allow high-throughput processing. Existing software solutions lack the raw performance and statistical grounding to efficiently handle the numerous modifications found on RNA. We present a free and open-source database search engine for RNA MS data, called NucleicAcidSearchEngine (NASE), that addresses these shortcomings. We demonstrate the capability of NASE to reliably identify a wide range of modified RNA sequences in four original datasets of varying complexity. In human tRNA, we characterize over 20 different modification types simultaneously and find many cases of incomplete modification.

An intellectual disability-associated missense variant in TRMT1 impairs tRNA modification and reconstitution of enzymatic activity

The human TRMT1 gene encodes an RNA methyltransferase enzyme responsible for catalyzing dimethylguanosine (m2,2G) formation in transfer RNAs (tRNAs). Frameshift mutations in TRMT1 have been shown to cause autosomal-recessive intellectual disability (ID) in the human population but additional TRMT1 variants remain to be characterized. Here, we describe a homozygous TRMT1 missense variant in a patient displaying developmental delay, ID, and epilepsy. The missense variant changes an arginine residue to a cysteine (R323C) within the methyltransferase domain and is expected to perturb protein folding. Patient cells expressing TRMT1-R323C exhibit a deficiency in m2,2G modifications within tRNAs, indicating that the mutation causes loss of function. Notably, the TRMT1 R323C mutant retains tRNA binding but is unable to rescue m2,2G formation in TRMT1-deficient human cells. Our results identify a pathogenic point mutation in TRMT1 that perturbs tRNA modification activity and demonstrate that m2,2G modifications are disrupted in the cells of patients with TRMT1-associated ID disorders.

NAIL-MS reveals the repair of 2-methylthiocytidine by AlkB in E. coli

RNAs contain post-transcriptional modifications, which fulfill a variety of functions in translation, secondary structure stabilization and cellular stress survival. Here, 2-methylthiocytidine (ms²C) is identified in tRNA of E. coli and P. aeruginosa using NAIL-MS (nucleic acid isotope labeling coupled mass spectrometry) in combination with genetic screening experiments. ms²C is only found in 2-thiocytidine (s²C) containing tRNAs, namely tRNA^Arg_CCG, tRNA^Arg_ICG, tRNA^Arg_UCU and tRNA^Ser_GCU at low abundances. ms²C is not formed by commonly known tRNA methyltransferases. Instead, we observe its formation in vitro and in vivo during exposure to methylating agents. More than half of the s²C containing tRNA can be methylated to carry ms²C. With a pulse-chase NAIL-MS experiment, the repair mechanism by AlkB dependent sulfur demethylation is demonstrated in vivo. Overall, we describe ms²C as a bacterial tRNA modification and damage product. Its repair by AlkB and other pathways is demonstrated in vivo by our powerful NAIL-MS approach.

Bacterial tRNA is modified by thiolation of nucleosides. Here the authors identify 2-methylthiocytidine in bacterial tRNA using nucleic acid isotope labeling coupled mass spectrometry. Exposure to methylating agents converts 2-thiocytidine to 2-methylthiocytidine, which is repaired by demethylase AlkB in vivo.

Bacterial wobble modifications of NNA-decoding tRNAs

Nucleotides of transfer RNAs (tRNAs) are highly modified, particularly at the anticodon. Bacterial tRNAs that read A-ending codons are especially notable. The U34 nucleotide canonically present in these tRNAs is modified by a wide range of complex chemical constituents. An additional two A-ending codons are not read by U34-containing tRNAs but are accommodated by either inosine or lysidine at the wobble position (I34 or L34). The structural basis for many N34 modifications in both tRNA aminoacylation and ribosome decoding has been elucidated, and evolutionary conservation of modifying enzymes is also becoming clearer. Here we present a brief review of the structure, function, and conservation of wobble modifications in tRNAs that translate A-ending codons. © 2019 IUBMB Life, 2019 © 2019 IUBMB Life, 71(8):1158-1166, 2019.

Cytosine-5 RNA methylation links protein synthesis to cell metabolism

Posttranscriptional modifications in transfer RNA (tRNA) are often critical for normal development because they adapt protein synthesis rates to a dynamically changing microenvironment. However, the precise cellular mechanisms linking the extrinsic stimulus to the intrinsic RNA modification pathways remain largely unclear. Here, we identified the cytosine-5 RNA methyltransferase NSUN2 as a sensor for external stress stimuli. Exposure to oxidative stress efficiently repressed NSUN2, causing a reduction of methylation at specific tRNA sites. Using metabolic profiling, we showed that loss of tRNA methylation captured cells in a distinct catabolic state. Mechanistically, loss of NSUN2 altered the biogenesis of tRNA-derived noncoding fragments (tRFs) in response to stress, leading to impaired regulation of protein synthesis. The intracellular accumulation of a specific subset of tRFs correlated with the dynamic repression of global protein synthesis. Finally, NSUN2-driven RNA methylation was functionally required to adapt cell cycle progression to the early stress response. In summary, we revealed that changes in tRNA methylation profiles were sufficient to specify cellular metabolic states and efficiently adapt protein synthesis rates to cell stress.

Matching tRNA modifications in humans to their known and predicted enzymes

tRNA are post-transcriptionally modified by chemical modifications that affect all aspects of tRNA biology. An increasing number of mutations underlying human genetic diseases map to genes encoding for tRNA modification enzymes. However, our knowledge on human tRNA-modification genes remains fragmentary and the most comprehensive RNA modification database currently contains information on approximately 20% of human cytosolic tRNAs, primarily based on biochemical studies. Recent high-throughput methods such as DM-tRNA-seq now allow annotation of a majority of tRNAs for six specific base modifications. Furthermore, we identified large gaps in knowledge when we predicted all cytosolic and mitochondrial human tRNA modification genes. Only 48% of the candidate cytosolic tRNA modification enzymes have been experimentally validated in mammals (either directly or in a heterologous system). Approximately 23% of the modification genes (cytosolic and mitochondrial combined) remain unknown. We discuss these 'unidentified enzymes' cases in detail and propose candidates whenever possible. Finally, tissue-specific expression analysis shows that modification genes are highly expressed in proliferative tissues like testis and transformed cells, but scarcely in differentiated tissues, with the exception of the cerebellum. Our work provides a comprehensive up to date compilation of human tRNA modifications and their enzymes that can be used as a resource for further studies.

Production and Application of Stable Isotope-Labeled Internal Standards for RNA Modification Analysis

Post-transcriptional RNA modifications have been found to be present in a wide variety of organisms and in different types of RNA. Nucleoside modifications are interesting due to their already known roles in translation fidelity, enzyme recognition, disease progression, and RNA stability. In addition, the abundance of modified nucleosides fluctuates based on growth phase, external stress, or possibly other factors not yet explored. With modifications ever changing, a method to determine absolute quantities for multiple nucleoside modifications is required. Here, we report metabolic isotope labeling to produce isotopically labeled internal standards in bacteria and yeast. These can be used for the quantification of 26 different modified nucleosides. We explain in detail how these internal standards are produced and show their mass spectrometric characterization. We apply our internal standards and quantify the modification content of transfer RNA (tRNA) from bacteria and various eukaryotes. We can show that the origin of the internal standard has no impact on the quantification result. Furthermore, we use our internal standard for the quantification of modified nucleosides in mouse tissue messenger RNA (mRNA), where we find different modification profiles in liver and brain tissue.

Lifestyle modifications: coordinating the tRNA epitranscriptome with codon bias to adapt translation during stress responses

Cells adapt to stress by altering gene expression at multiple levels. Here, we propose a new mechanism regulating stress-dependent gene expression at the level of translation, with coordinated interplay between the tRNA epitranscriptome and biased codon usage in families of stress-response genes. In this model, auxiliary genetic information contained in synonymous codon usage enables regulation of codon-biased and functionally related transcripts by dynamic changes in the tRNA epitranscriptome. This model partly explains the association between synchronous stress-dependent epitranscriptomic marks and significant multi-codon usage skewing in families of translationally regulated transcripts. The model also predicts translational adaptation during viral infections.

NAIL-MS in E. coli Determines the Source and Fate of Methylation in tRNA

In all domains of life, the nucleobases of tRNA can be methylated. These methylations are introduced either by enzymes or by the reaction of methylating agents with the nucleophilic centers of the nucleobases. Herein, we present a systematic approach to identify the methylation sites within RNA in vitro and in vivo. For discrimination between enzymatic tRNA methylation and tRNA methylation damage in bacteria, we used nucleic acid isotope labeling coupled mass spectrometry (NAIL‐MS). With NAIL‐MS, we clearly observed the formation of 7‐methylguanosine, 3‐methyluridine, and 6‐methyladenosine during exposure of bacteria to the alkylating agent methyl methanesulfonate (MMS) in vivo. These damage products were not reported to form in tRNA in vivo, as they were masked by the enzymatically formed modified nucleosides in previous studies. In addition, we found formation of the known damage products 1‐methyladenosine and 3‐methylcytidine in vivo. With a dynamic NAIL‐MS setup, we observed tRNA repair by demethylation of these two RNA modifications in vivo. Furthermore, we saw the potential repair of 6‐methyladenosine but not 7‐methylguanosine in bacterial tRNA.

Detection of ribonucleoside modifications by liquid chromatography coupled with mass spectrometry

A small set of ribonucleoside modifications have been found in different regions of mRNA including the open reading frame. Accurate detection of these specific modifications is critical to understanding their modulatory roles in facilitating mRNA maturation, translation and degradation. While transcriptome-wide next-generation sequencing (NGS) techniques could provide exhaustive information about the sites of one specific or class of modifications at a time, recent investigations strongly indicate cautionary interpretation due to the appearance of false positives. Therefore, it is suggested that NGS-based modification data can only be treated as predicted sites and their existence need to be validated by orthogonal methods. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) is an analytical technique that can yield accurate and reproducible information about the qualitative and quantitative characteristics of ribonucleoside modifications. Here, we review the recent advancements in LC-MS/MS technology that could help in securing accurate, gold-standard quality information about the resident post-transcriptional modifications of mRNA.

Improved application of RNAModMapper - An RNA modification mapping software tool - For analysis of liquid chromatography tandem mass spectrometry (LC-MS/MS) data

Research into post-transcriptional processing and modification of RNA continues to speed forward, as their ever-emerging role in the regulation of gene expression in biological systems continues to unravel. Liquid chromatography tandem mass spectrometry (LC-MS/MS) has proven for over two decades to be a powerful ally in the elucidation of RNA modification identity and location, but the technique has not proceeded without its own unique technical challenges. The throughput of LC-MS/MS modification mapping experiments continues to be impeded by tedious and time-consuming spectral interpretation, particularly during for the analysis of complex RNA samples. RNAModMapper was recently developed as a tool to improve the interpretation and annotation of LC-MS/MS data sets from samples containing post-transcriptionally modified RNAs. Here, we delve deeper into the methodology and practice of RNAModMapper to provide greater insight into its utility, and remaining hurdles, in current RNA modification mapping experiments.

tRNA epitranscriptomics and biased codon are linked to proteome expression in Plasmodium falciparum

Among components of the translational machinery, ribonucleoside modifications on tRNAs are emerging as critical regulators of cell physiology and stress response. Here, we demonstrate highly coordinated behavior of the repertoire of tRNA modifications of Plasmodium falciparum throughout the intra-erythrocytic developmental cycle (IDC). We observed both a synchronized increase in 22 of 28 modifications from ring to trophozoite stage, consistent with tRNA maturation during translational up-regulation, and asynchronous changes in six modifications. Quantitative analysis of ~2,100 proteins across the IDC revealed that up- and down-regulated proteins in late but not early stages have a marked codon bias that directly correlates with parallel changes in tRNA modifications and enhanced translational efficiency. We thus propose a model in which tRNA modifications modulate the abundance of stage-specific proteins by enhancing translation efficiency of codon-biased transcripts for critical genes. These findings reveal novel epitranscriptomic and translational control mechanisms in the development and pathogenesis of Plasmodium parasites.

The role of RNA modifications in the regulation of tRNA cleavage

Transfer RNA (tRNA) have been harbingers of many paradigms in RNA biology. They are among the first recognized noncoding RNA (ncRNA) playing fundamental roles in RNA metabolism. Although mainly recognized for their role in decoding mRNA and delivering amino acids to the growing polypeptide chain, tRNA also serve as an abundant source of small ncRNA named tRNA fragments. The functional significance of these fragments is only beginning to be uncovered. Early on, tRNA were recognized as heavily post-transcriptionally modified, which aids in proper folding and modulates the tRNA:mRNA anticodon-codon interactions. Emerging data suggest that these modifications play critical roles in the generation and activity of tRNA fragments. Modifications can both protect tRNA from cleavage or promote their cleavage. Modifications to individual fragments may be required for their activity. Recent work has shown that some modifications are critical for stem cell development and that failure to deposit certain modifications has profound effects on disease. This review will discuss how tRNA modifications regulate the generation and activity of tRNA fragments.

Differentiating Positional Isomers of Nucleoside Modifications by Higher-Energy Collisional Dissociation Mass Spectrometry (HCD MS)

The analytical identification of positional isomers (e.g., 3-, N⁴-, 5-methylcytidine) within the > 160 different post-transcriptional modifications found in RNA can be challenging. Conventional liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS) approaches rely on chromatographic separation for accurate identification because the collision-induced dissociation (CID) mass spectra of these isomers nearly exclusively yield identical nucleobase ions (BH₂⁺) from the same molecular ion (MH⁺). Here, we have explored higher-energy collisional dissociation (HCD) as an alternative fragmentation technique to generate more informative product ions that can be used to differentiate positional isomers. LC-MS/MS of modified nucleosides characterized using HCD led to the creation of structure- and HCD energy-specific fragmentation patterns that generated unique fingerprints, which can be used to identify individual positional isomers even when they cannot be separated chromatographically. While particularly useful for identifying positional isomers, the fingerprinting capabilities enabled by HCD also offer the potential to generate HPLC-independent spectral libraries for the rapid analysis of modified ribonucleosides. Graphical Abstract ᅟ.