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

Nucleo-SAFARI: Automated Identification of Fragment Ions in Top-Down MS/MS Spectra of Nucleic Acids

Recent progress in top-down mass spectrometry analysis of progressively larger nucleic acids has enabled in-depth characterization of intact, modified RNA molecules. Development of methods for desalting and MS/MS fragmentation allows rapid acquisition of high-quality top-down MS/MS spectra of nucleic acids up to 100 nt, which has spurred the need for development of software approaches to identify and validate nucleic acid fragment ions. We have implemented an R-based approach to aid in analysis of MS/MS spectra of nucleic acids based on fragment ions observed directly in the m/z domain. This program, entitled Shiny Application for Fragment Assignment by Relative Isotopes (Nucleo-SAFARI), utilizes the Shiny HTML framework for deployment of a user-friendly application for automated annotation of top-down MS/MS spectra of nucleic acids recorded on Orbitrap mass spectrometer platforms. This approach proceeds through in silico generation of fragment ions and their isotopic distributions, followed by algorithmic assessment of the experimental isotopic distributions. Nucleo-SAFARI is available for download at https://github.com/mblanzillotti/Nucleo-SAFARI.

Characterization of mRNA therapeutics

Therapeutic messenger RNAs (mRNAs) have emerged as powerful tools in the treatment of complex diseases, especially for conditions that lack efficacious treatment. The successful application of this modality can be attributed to its ability to encode entire proteins. While the large nature of these molecules has supported their success as therapeutics, its extended size creates several analytical challenges. To further support therapeutic mRNA development and its deployment in clinical trials, appropriate methods to support their characterization must be developed. In this review, we describe current analytical methods that have been used in the characterization of RNA quality, identity, and integrity. Advantages and limitations from several analytical techniques ranging from gel electrophoresis to liquid chromatography-mass spectrometry and from shotgun sequencing to intact mass measurements are discussed. We comprehensively describe the application of analytical methods in the measurements of capping efficiency, poly A tail analysis, as well as their applicability in stability studies.

Analysis of RNA and Its Modifications

Ribonucleic acids (RNAs) are key biomolecules responsible for the transmission of genetic information, the synthesis of proteins, and modulation of many biochemical processes. They are also often the key components of viruses. Synthetic RNAs or oligoribonucleotides are becoming more widely used as therapeutics. In many cases, RNAs will be chemically modified, either naturally via enzymatic systems within a cell or intentionally during their synthesis. Analytical methods to detect, sequence, identify, and quantify RNA and its modifications have demands that far exceed requirements found in the DNA realm. Two complementary platforms have demonstrated their value and utility for the characterization of RNA and its modifications: mass spectrometry and next-generation sequencing. This review highlights recent advances in both platforms, examines their relative strengths and weaknesses, and explores some alternative approaches that lie at the horizon.

The role of post-transcriptional modification on a new tRNAIle(GAU) identified from Ganoderma lucidum in its fragments' cytotoxicity on cancer cells

Ganoderma lucidum (Ganoderma) is a famous Chinese herbal medicine which has been used clinically for thousands of years in China. Despite numerous studies on triterpenes and polysaccharides, the bioactivity of RNAs abundant in Ganoderma remains unknown. Here, based on LC-MS techniques, dihydrouracil, 5-methyluridine (m⁵U) and pseudouridine were identified at position 19, 52 and 53 of a new tRNA^Ile(GAU) which was isolated as the most abundant tRNA species in Ganoderma, and is the first purified tRNA from fungus. Cytotoxic screening of tRNA-half (t-half) and tRNA fragment (tRF) derived from this tRNA, as well as their mimics (t-half or tRF as antisense strand), demonstrated that the double-stranded form, i.e., tRF and t-halve mimics, exhibited stronger cytotoxicity than their single-stranded form, and the cytotoxicity of t-half mimic is significantly stronger than that of tRF mimic. Notably, the cytotoxicity of 3'-t-half mimic is not only much more potent than that of taxol, but also is much more potent than that of ganoderic acids, the major bioactive components in Ganoderma. Furthermore, 3'-t-half mimic_M2 (m⁵U modified) exhibited significantly stronger cytotoxicity than unmodified 3'-t-half mimic, which is consistent with the computational simulation showing that m⁵U modification enhances the stability of the tertiary structure of 3'-t-half mimic. Overall, the present study not only indicates t-halves are bioactive components in Ganoderma which should not be neglected, but also reveals an important role of post-transcriptional modification on tRNA in its fragments' cytotoxicity against cancer cells, which benefits the design and development of RNAi drugs from natural resource.

Characterizing RNA modifications in the central nervous system and single cells by RNA sequencing and liquid chromatography-tandem mass spectrometry techniques

Post-transcriptional modifications to RNA constitute a newly appreciated layer of translation regulation in the central nervous system (CNS). The identity, stoichiometric quantity, and sequence position of these unusual epitranscriptomic marks are central to their function, making analytical methods that are capable of accurate and reproducible measurements paramount to the characterization of the neuro-epitranscriptome. RNA sequencing-based methods and liquid chromatography-tandem mass spectrometry (LC-MS/MS) techniques have been leveraged to provide an early glimpse of the landscape of RNA modifications in bulk CNS tissues. However, recent advances in sample preparation, separations, and detection methods have revealed that individual cells display remarkable heterogeneity in their RNA modification profiles, raising questions about the prevalence and function of cell-specific distributions of post-transcriptionally modified nucleosides in the brain. In this Trends article, we present an overview of RNA sequencing and LC-MS/MS methodologies for the analysis of RNA modifications in the CNS with special emphasis on recent advancements in techniques that facilitate single-cell and subcellular detection.

Challenges and emerging trends in liquid chromatography-based analyses of mRNA pharmaceuticals

Lipid encapsulated messenger RNA (LNP mRNA) has garnered a significant amount of interest from the pharmaceutical industry and general public alike. This attention has been catalyzed by the clinical success of LNP mRNA for SARS-CoV-2 vaccination as well as future promises that might be fulfilled by the biotechnology pipeline, such as the in vivo delivery of a CRISPR/Cas9 complex that can edit patient cells to reduce levels of low-density lipoprotein. LNP mRNAs are comprised of various chemically diverse molecules brought together in a sophisticated intermolecular complex. This can make it challenging to achieve thorough analytical characterization. Nevertheless, liquid chromatography is becoming an increasingly relied upon technique for LNP mRNA analyses. Although there have been significant advances in all types of LNP mRNA analyses, this review focuses on recent developments and the possibilities of applying anion exchange (AEX) and ion pairing reversed phase (IP-RP) liquid chromatography for intact mRNAs as well as techniques for oligo mapping analysis, 5' endcap testing and lipid compositional assays.

RNA Cleavage Properties of Nucleobase-Specific RNase MC1 and Cusativin Are Determined by the Dinucleotide-Binding Interactions in the Enzyme-Active Site

Knowledge of the cleavage specificity of ribonucleases is critical for their application in RNA modification mapping or RNA-protein binding studies. Here, we detail the cleavage specificity and efficiency of ribonuclease MC1 and cusativin using a customized RNA sequence that contained all dinucleotide combinations and homopolymer sequences. The sequencing of the oligonucleotide digestion products by a semi-quantitative liquid chromatography coupled with mass spectrometry (LC-MS) analysis documented as little as 0.5-1% cleavage levels for a given dinucleotide sequence combination. While RNase MC1 efficiently cleaved the [A/U/C]pU dinucleotide bond, no cleavage was observed for the GpU bond. Similarly, cusativin efficiently cleaved Cp[U/A/G] dinucleotide combinations along with UpA and [A/U]pU, suggesting a broader specificity of dinucleotide preferences. The molecular interactions between the substrate and active site as determined by the dinucleotide docking studies of protein models offered additional evidence and support for the observed substrate specificity. Targeted alteration of the key amino acid residues in the nucleotide-binding site confirms the utility of this in silico approach for the identification of key interactions. Taken together, the use of bioanalytical and computational approaches, involving LC-MS and ligand docking of tertiary structural models, can form a powerful combination to help explain the RNA cleavage behavior of RNases.

Technical challenges in defining RNA modifications

It is well established that DNA base modifications play a key role in gene regulation during development and in response to environmental stress. This type of epigenetic control of development and environmental responses has been intensively studied over the past few decades. Similar to DNA, various RNA species also undergo modifications that play important roles in, for example, RNA splicing, protein translation, and the avoidance of immune surveillance by host. More than 160 different types of RNA modifications have been identified. In addition to base modifications, RNA modification also involves splicing of pre-mRNAs, leading to as many as tens of transcript isoforms from a single pre-RNA, especially in higher organisms. However, the function, prevalence and distribution of RNA modifications are poorly understood. The lack of a suitable method for the reliable identification of RNA modifications constitutes a significant challenge to studying their functions. This review focuses on the technologies that enable de novo identification of RNA base modifications and the alternatively spliced mRNA transcripts.

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.

Analysis of the epitranscriptome with ion-pairing reagent free oligonucleotide mass spectrometry

RNA modifications gain growing attention as a new frontier in the life sciences but with the rise of RNA vaccines also in biomedical drug development. Impeccable characterization of RNA modifications within their sequence context remains an analytical challenge. Oligonucleotide mass spectrometry (ON-MS), an approach similar to bottom-up proteome analysis, is capable of defining a short 5-15 nucleotide sequence context of an RNA modification while delivering information on the chemical character of the modified nucleotide. Commonly, ON-MS requires the use of ion pairing reagents for ON separation which is not compatible with most other MS-based applications and only few laboratories run dedicated MS instruments for the task. Here, we present an ON-MS technique which is independent of ion pairing reagents and can be used on any available mass spectrometer without risking its sensitivity for other analytes. In this chapter, we describe the experiments necessary for ON-MS method development, ON-MS application to native and synthetic RNAs and finally a guideline for data analysis.

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.

Full-Range Profiling of tRNA Modifications Using LC-MS/MS at Single-Base Resolution through a Site-Specific Cleavage Strategy

Transfer RNAs (tRNAs) are the most heavily modified RNA species. Liquid chromatography coupled with mass spectrometry (LC-MS/MS) is a powerful tool for characterizing tRNA modifications, which involves pretreating tRNAs with base-specific ribonucleases to produce smaller oligonucleotides amenable to MS. However, the quality and quantity of products from base-specific digestions are severely impacted by the base composition of tRNAs. This often leads to a loss of sequence information. Here, we report a method for the full-range profiling of tRNA modifications at single-base resolution by combining site-specific RNase H digestion with the LC-MS/MS and RNA-seq techniques. The key steps were designed to generate high-quality products of optimal lengths and ionization properties. A linear correlation between collision energies and the m/z of oligonucleotides significantly improved the information content of collision-induced dissociation (CID) spectra. False positives were eliminated by up to 95% using novel inclusion criteria for collecting a census of modifications. This method is illustrated by the mapping of mouse mitochondrial tRNA^His(GUG) and tRNA^Val(UAC), which were hitherto not investigated. The identities and locations of the five species of modifications on these tRNAs were fully characterized. This approach is universally applicable to any tRNA species and provides an experimentally realizable pathway to the de novo sequencing of post-transcriptionally modified tRNAs with high sequence coverage.

Aom2 S: A new web-based application for DNA/RNA tandem mass spectrometry data interpretation

RATIONALE

The Analysis of Oligonucleotide Modifications from Mass Spectra (Aom² S) was created to support the analysis of oligonucleotide mass spectra. This application complements the existing software tools by providing a comprehensive analysis of oligonucleotide fragments from high-resolution tandem mass spectrometry (HR-MS/MS) data in a flexible and user-friendly manner, directly accessible through a web browser without any need for installation.

METHODS

MS measurements of aminoC6-DNA and inosine-RNA were performed using an LTQ Orbitrap FT-MS instrument. The obtained data were analyzed by our newly developed open-source package Aom² S accessible from the ms.epfl.ch web page or directly at https://mstools.epfl.ch/am2s/ to demonstrate the various functionalities of this tool, notably the possibility to identify different product ions from a nucleotide sequence with any fixed/variable modification by matching theoretical isotopic patterns to any experimental mass spectra with similarity scores ranking.

RESULTS

A detailed description of the Aom² S tool with its user-friendly interface is exemplified using HR-MS/MS data of modified DNA and RNA oligonucleotides. Explanations of analysis parameters and tool workflow, as well as multiple options for viewing and exporting the results, are provided. Product ion assignment and modification localization can be achieved in seconds, and results can be exported as tables, matched mass spectra, and fragmentation maps.

CONCLUSIONS

A new open source tool (Aom² S) for the analysis of HR-MS/MS data for modified DNA and RNA oligonucleotides is described. Aom² S is fast, highly flexible, and versatile, allowing automatic precursor and product ion assignment in a comprehensive manner, including internal fragments and variable modification localization, with clear graphical representation of the results.

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.

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.

Improved RNA modification mapping of cellular non-coding RNAs using C- and U-specific RNases

Locating ribonucleoside modifications within an RNA sequence requires digestion of the RNA into oligoribonucleotides of amenable size for subsequent analysis by LC-MS (liquid chromatography-mass spectrometry). This approach, widely referred to as RNA modification mapping, is facilitated through ribonucleases (RNases) such as T1 (guanosine-specific), U2 (purine-selective) and A (pyrimidine-specific) among others. Sequence coverage by these enzymes depends on positioning of the recognized nucleobase (such as guanine or purine or pyrimidine) in the sequence and its ribonucleotide composition. Using E. coli transfer RNA (tRNA) and ribosomal RNA (rRNA) as model samples, we demonstrate the ability of complementary nucleobase-specific ribonucleases cusativin (C-specific) and MC1 (U-specific) to generate digestion products that facilitate confident mapping of modifications in regions such as G-rich and pyrimidine-rich segments of RNA, and to distinguish C to U sequence differences. These enzymes also increase the number of oligonucleotide digestion products that are unique to a specific RNA sequence. Further, with these additional RNases, multiple modifications can be localized with high confidence in a single set of experiments with minimal dependence on the individual tRNA abundance in a mixture. The sequence overlaps observed with these complementary digestion products and that of RNase T1 improved sequence coverage to 75% or above. A similar level of sequence coverage was also observed for the 2904 nt long 23S rRNA indicating their utility has no dependence on RNA size. Wide-scale adoption of these additional modification mapping tools could help expedite the characterization of modified RNA sequences to understand their structural and functional role in various living systems.

Methods and Challenges for Computational Data Analysis for DNA Adductomics

Frequent exposure to chemicals in the environment, diet, and endogenous electrophiles leads to chemical modification of DNA and the formation of DNA adducts. Some DNA adducts can induce mutations during cell division and, when occurring in critical regions of the genome, can lead to the onset of disease, including cancer. The targeted analysis of DNA adducts over the past 30 years has revealed that the human genome contains many types of DNA damages. However, a long-standing limitation in conducting DNA adduct measurements has been the inability to screen for the total complement of DNA adducts derived from a wide range of chemicals in a single assay. With the advancement of high-resolution mass spectrometry (MS) instrumentation and new scanning technologies, nontargeted "omics" approaches employing data-dependent acquisition and data-independent acquisition methods have been established to simultaneously screen for multiple DNA adducts, a technique known as DNA adductomics. However, notable challenges in data processing must be overcome for DNA adductomics to become a mature technology. DNA adducts occur at low abundance in humans, and current softwares do not reliably detect them when using common MS data acquisition methods. In this perspective, we discuss contemporary computational tools developed for feature finding of MS data widely utilized in the disciplines of proteomics and metabolomics and highlight their limitations for conducting nontargeted DNA-adduct biomarker discovery. Improvements to existing MS data processing software and new algorithms for adduct detection are needed to develop DNA adductomics into a powerful tool for the nontargeted identification of potential cancer-causing agents.

Using spectral matching to interpret LC-MS/MS data during RNA modification mapping

While a number of approaches have been developed to analyze liquid chromatography tandem mass spectrometry (LC-MS/MS) data obtained from modified oligonucleotides, the majority of these methods require analyzing every MS/MS spectrum de novo to sequence the oligonucleotide and place the modification. Spectral matching is an alternative approach for analyzing MS/MS data that is based on creating a library of annotated MS/MS spectra against which individual MS/MS data can be searched. Here, we have adapted the existing NIST spectral matching software to enable its use in the interpretation of MS/MS data obtained from modified oligonucleotides. In particular, we demonstrate the utility of this approach to identify specific post-transcriptionally modified nucleosides in particular transfer RNAs (tRNAs) obtained through a conventional RNA modification mapping experimental protocol. Spectral matching was found to be an efficient approach for screening for known modified tRNAs by using the experimental data as the library and previously annotated RNase T1 digestion products of tRNAs as the reference spectra. The utility of spectral matching for rapid analysis of multiple LC-MS/MS analyses was demonstrated by screening mutant strains of Streptococcus mutans to identify the enzyme(s) responsible for synthesizing the tRNA position 37 modification threonylcarbamoyladenosine (t⁶ A).

Distinct substrate specificities of the human tRNA methyltransferases TRMT10A and TRMT10B

The tRNA m1R9 methyltransferase (Trm10) family is conserved throughout Eukarya and Archaea. Despite the presence of a single Trm10 gene in Archaea and most single-celled eukaryotes, metazoans encode up to three homologs of Trm10. Several disease states correlate with a deficiency in the human homolog TRMT10A, despite the presence of another cytoplasmic enzyme, TRMT10B. Here we investigate these phenomena and demonstrate that human TRMT10A (hTRMT10A) and human TRMT10B (hTRMT10B) are not biochemically redundant. In vitro activity assays with purified hTRMT10A and hTRMT10B reveal a robust activity for hTRMT10B as a tRNAAsp-specific m1A9 methyltransferase and suggest that it is the relevant enzyme responsible for this newly discovered m1A9 modification in humans. Moreover, a comparison of the two cytosolic enzymes with multiple tRNA substrates exposes the enzymes' distinct substrate specificities, and suggests that hTRMT10B exhibits a restricted selectivity hitherto unseen in the Trm10 enzyme family. Single-turnover kinetics and tRNA binding assays highlight further differences between the two enzymes and eliminate overall tRNA affinity as a primary determinant of substrate specificity for either enzyme. These results increase our understanding of the important biology of human tRNA modification systems, which can aid in understanding the molecular basis for diseases in which their aberrant function is increasingly implicated.

Emerging Technologies in Mass Spectrometry-Based DNA Adductomics

The measurement of DNA adducts, the covalent modifications of DNA upon the exposure to the environmental and dietary genotoxicants and endogenously produced electrophiles, provides molecular evidence for DNA damage. With the recent improvements in the sensitivity and scanning speed of mass spectrometry (MS) instrumentation, particularly high-resolution MS, it is now feasible to screen for the totality of DNA damage in the human genome through DNA adductomics approaches. Several MS platforms have been used in DNA adductomic analysis, each of which has its strengths and limitations. The loss of 2′-deoxyribose from the modified nucleoside upon collision-induced dissociation is the main transition feature utilized in the screening of DNA adducts. Several advanced data-dependent and data-independent scanning techniques originated from proteomics and metabolomics have been tailored for DNA adductomics. The field of DNA adductomics is an emerging technology in human exposure assessment. As the analytical technology matures and bioinformatics tools become available for analysis of the MS data, DNA adductomics can advance our understanding about the role of chemical exposures in DNA damage and disease risk.

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.