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.

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

The dynamic RNA modification 5-methylcytosine and its emerging role as an epitranscriptomic mark

It is a well‐known fact that RNA is the target of a plethora of modifications which currently amount to over a hundred. The vast majority of these modifications was observed in the two most abundant classes of RNA, rRNA and tRNA. With the recent advance in mapping technologies, modifications have been discovered also in mRNA and in less abundant non‐coding RNA species. These developments have sparked renewed interest in elucidating the nature and functions of those “epitransciptomic” modifications in RNA. N6‐methyladenosine (m⁶A) is the best understood and most frequent mark of mRNA with demonstrated functions ranging from pre‐mRNA processing, translation, miRNA biogenesis to mRNA decay. By contrast, much less research has been conducted on 5‐methylcytosine (m5C), which was detected in tRNAs and rRNAs and more recently in poly(A)RNAs. In this review, we discuss recent developments in the discovery of m5C RNA methylomes, the functions of m5C as well as the proteins installing, translating and manipulating this modification. Although our knowledge about m5C in RNA transcripts is just beginning to consolidate, it has become clear that cytosine methylation represents a powerful mechanistic strategy to regulate cellular processes on an epitranscriptomic level.

This article is categorized under:

RNA Processing > RNA Editing and Modification

RNA Interactions with Proteins and Other Molecules > Protein–RNA Interactions: Functional Implications

RNA Processing > tRNA Processing

RNA Turnover and Surveillance > Regulation of RNA Stability

Mapping Post-Transcriptional Modifications onto Transfer Ribonucleic Acid Sequences by Liquid Chromatography Tandem Mass Spectrometry

Liquid chromatography, coupled with tandem mass spectrometry, has become one of the most popular methods for the analysis of post-transcriptionally modified transfer ribonucleic acids (tRNAs). Given that the information collected using this platform is entirely determined by the mass of the analyte, it has proven to be the gold standard for accurately assigning nucleobases to the sequence. For the past few decades many labs have worked to improve the analysis, contiguous to instrumentation manufacturers developing faster and more sensitive instruments. With biological discoveries relating to ribonucleic acid happening more frequently, mass spectrometry has been invaluable in helping to understand what is happening at the molecular level. Here we present a brief overview of the methods that have been developed and refined for the analysis of modified tRNAs by liquid chromatography tandem mass spectrometry.

Going global: the new era of mapping modifications in RNA

The post-transcriptional modification of RNA by the addition of one or more chemical groups has been known for over 50 years. These chemical modifications, once thought to be static, are now being discovered to play key regulatory roles in gene expression. The advent of massive parallel sequencing of RNA (RNA-seq) now allows us to probe the complexity of cellular RNA and how chemically altering RNA structure expands the RNA vocabulary. Here we present an overview of the various strategies and technologies that are available to profile RNA chemical modifications at the cellular level. These strategies can be characterized as targeted and untargeted approaches: targeted strategies are developed for one single chemical modification while untargeted strategies are more broadly applicable to a range of such chemical changes. Key for all of these approaches is the ability to locate modifications within the RNA sequence. While most of these methods are built upon an RNA-Seq pipeline, alternative approaches based on mass spectrometry or conventional DNA sequencing retain value in the overall analysis process. We also look forward toward future opportunities and technologies that may expand the types of modifications that can be globally profiled. Given the ever increasing recognition that these RNA chemical modifications play important biological roles, a variety of methods, preferably orthogonal approaches, will be required to globally identify, validate and quantify RNA chemical modifications found in the transcriptome. WIREs RNA 2017, 8:e1367. doi: 10.1002/wrna.1367 For further resources related to this article, please visit the WIREs website.

tRNA-mediated codon-biased translation in mycobacterial hypoxic persistence

Microbial pathogens adapt to the stress of infection by regulating transcription, translation and protein modification. We report that changes in gene expression in hypoxia-induced non-replicating persistence in mycobacteria—which models tuberculous granulomas—are partly determined by a mechanism of tRNA reprogramming and codon-biased translation. Mycobacterium bovis BCG responded to each stage of hypoxia and aerobic resuscitation by uniquely reprogramming 40 modified ribonucleosides in tRNA, which correlate with selective translation of mRNAs from families of codon-biased persistence genes. For example, early hypoxia increases wobble cmo⁵U in tRNA^Thr(UGU), which parallels translation of transcripts enriched in its cognate codon, ACG, including the DosR master regulator of hypoxic bacteriostasis. Codon re-engineering of dosR exaggerates hypoxia-induced changes in codon-biased DosR translation, with altered dosR expression revealing unanticipated effects on bacterial survival during hypoxia. These results reveal a coordinated system of tRNA modifications and translation of codon-biased transcripts that enhance expression of stress response proteins in mycobacteria.

Mycobacteria can adapt to the stress of human infection by entering a dormant state. Here the authors show that hypoxia-induced dormancy in M. bovis BCG involves the reprogramming of tRNA wobble modifications and copy numbers, coupled with biased use of synonymous codons in survival genes.

Quantifying the 'escapers' among RNA species

tRNAs are fundamental components of translation and emerging evidence places them more centrally in various other cellular processes. However, rather than being uniformly conserved, tRNA abundance is instead highly variable and adaptable. The amount of tRNA genes greatly differs among species. Moreover, even within the same genome, tRNA abundance shapes the proteome in a tissue- and cell-specific manner and is dynamically regulated in response to stress. Here, we review approaches for identification and quantification of tRNAs and their functional integrity. We discuss the resolution of each method and highlight new approaches with cell-wide resolution based on deep-sequencing technologies.

Sequence mapping of transfer RNA chemical modifications by liquid chromatography tandem mass spectrometry

Mass spectrometry is a powerful analytical tool for identifying and characterizing structural modifications to the four canonical bases in RNA, information that is lost when using techniques such as PCR for RNA analysis. Here we described an updated method for sequence mapping of modified nucleosides in transfer RNA. This modification mapping approach utilizes knowledge of the modified nucleosides present in the sample along with the genome-derived tRNA sequence to readily locate modifications site-specifically in the tRNA sequence. The experimental approach involves isolation of the tRNA of interest followed by separate enzymatic digestion to nucleosides and oligonucleotides. Both samples are analyzed by liquid chromatography tandem mass spectrometry (LC-MS/MS) and the data sets are then combined to yield the modification profile of the tRNA. Data analysis is facilitated by the use of unmodified sequence exclusion lists and new developments in software that can automate MS/MS spectral annotation. The method is illustrated using tRNA-Asn isolated from Thermus thermophilus.

Mass spectrometry of modified RNAs: recent developments

A common feature of ribonucleic acids (RNAs) is that they can undergo a variety of chemical modifications. As nearly all of these chemical modifications result in an increase in the mass of the canonical nucleoside, mass spectrometry has long been a powerful approach for identifying and characterizing modified RNAs. Over the past several years, significant advances have been made in method development and software for interpreting tandem mass spectra resulting in approaches that can yield qualitative and quantitative information on RNA modifications, often at the level of sequence specificity. We discuss these advances along with instrumentation developments that have increased our ability to extract such information from relatively complex biological samples. With the increasing interest in how these modifications impact the epitranscriptome, mass spectrometry will continue to play an important role in bioanalytical investigations revolving around RNA.

The identification and characterization of non-coding and coding RNAs and their modified nucleosides by mass spectrometry

The analysis of ribonucleic acids (RNA) by mass spectrometry has been a valuable analytical approach for more than 25 years. In fact, mass spectrometry has become a method of choice for the analysis of modified nucleosides from RNA isolated out of biological samples. This review summarizes recent progress that has been made in both nucleoside and oligonucleotide mass spectral analysis. Applications of mass spectrometry in the identification, characterization and quantification of modified nucleosides are discussed. At the oligonucleotide level, advances in modern mass spectrometry approaches combined with the standard RNA modification mapping protocol enable the characterization of RNAs of varying lengths ranging from low molecular weight short interfering RNAs (siRNAs) to the extremely large 23 S rRNAs. New variations and improvements to this protocol are reviewed, including top-down strategies, as these developments now enable qualitative and quantitative measurements of RNA modification patterns in a variety of biological systems.

Essentiality of threonylcarbamoyladenosine (t(6)A), a universal tRNA modification, in bacteria

Threonylcarbamoyladenosine (t(6)A) is a modified nucleoside universally conserved in tRNAs in all three kingdoms of life. The recently discovered genes for t(6)A synthesis, including tsaC and tsaD, are essential in model prokaryotes but not essential in yeast. These genes had been identified as antibacterial targets even before their functions were known. However, the molecular basis for this prokaryotic-specific essentiality has remained a mystery. Here, we show that t(6)A is a strong positive determinant for aminoacylation of tRNA by bacterial-type but not by eukaryotic-type isoleucyl-tRNA synthetases and might also be a determinant for the essential enzyme tRNA(Ile)-lysidine synthetase. We confirm that t(6)A is essential in Escherichia coli and a survey of genome-wide essentiality studies shows that genes for t(6)A synthesis are essential in most prokaryotes. This essentiality phenotype is not universal in Bacteria as t(6)A is dispensable in Deinococcus radiodurans, Thermus thermophilus, Synechocystis PCC6803 and Streptococcus mutans. Proteomic analysis of t(6)A(-) D. radiodurans strains revealed an induction of the proteotoxic stress response and identified genes whose translation is most affected by the absence of t(6)A in tRNAs. Thus, although t(6)A is universally conserved in tRNAs, its role in translation might vary greatly between organisms.

RoboOligo: software for mass spectrometry data to support manual and de novo sequencing of post-transcriptionally modified ribonucleic acids

Ribosomal ribonucleic acid (RNA), transfer RNA and other biological or synthetic RNA polymers can contain nucleotides that have been modified by the addition of chemical groups. Traditional Sanger sequencing methods cannot establish the chemical nature and sequence of these modified-nucleotide containing oligomers. Mass spectrometry (MS) has become the conventional approach for determining the nucleotide composition, modification status and sequence of modified RNAs. Modified RNAs are analyzed by MS using collision-induced dissociation tandem mass spectrometry (CID MS/MS), which produces a complex dataset of oligomeric fragments that must be interpreted to identify and place modified nucleosides within the RNA sequence. Here we report the development of RoboOligo, an interactive software program for the robust analysis of data generated by CID MS/MS of RNA oligomers. There are three main functions of RoboOligo: (i) automated de novo sequencing via the local search paradigm. (ii) Manual sequencing with real-time spectrum labeling and cumulative intensity scoring. (iii) A hybrid approach, coined 'variable sequencing', which combines the user intuition of manual sequencing with the high-throughput sampling of automated de novo sequencing.

Identification of RNA sequence isomer by isotope labeling and LC-MS/MS

Recently, we developed a method for modified ribonucleic acid (RNA) analysis based on the comparative analysis of RNA digests (CARD). Within this CARD approach, sequence or modification differences between two samples are identified through differential isotopic labeling of two samples. Components present in both samples will each be labeled, yielding doublets in the CARD mass spectrum. Components unique to only one sample should be detected as singlets. A limitation of the prior singlet identification strategy occurs when the two samples contain components of unique sequence but identical base composition. At the first stage of mass spectrometry, these sequence isomers cannot be differentiated and would appear as doublets rather than singlets. However, underlying sequence differences should be detectable by collision-induced dissociation tandem mass spectrometry (CID MS/MS), as y-type product ions will retain the original enzymatically incorporated isotope label. Here, we determine appropriate instrumental conditions that enable CID MS/MS of isotopically labeled ribonuclease T1 (RNase T1) digestion products such that the original isotope label is maintained in the product ion mass spectrum. Next, we demonstrate how y-type product ions can be used to differentiate singlets and doublets from isomer sequences. We were then able to extend the utility of this approach by using CID MS/MS for the confirmation of an expected RNase T1 digestion product within the CARD analysis of an Escherichia coli mutant strain even in the presence of interfering and overlapping digestion products from other transfer RNAs.

Systematic identification of tRNAome and its dynamics in Lactococcus lactis

Transfer RNAs (tRNAs) through their abundance and modification pattern significantly influence protein translation. Here, we present a systematic analysis of the tRNAome of Lactococcus lactis. Using the next-generation sequencing approach, we identified 40 tRNAs which carry 16 different post-transcriptional modifications as revealed by mass spectrometry analysis. While small modifications are located in the tRNA body, hypermodified nucleotides are mainly present in the anticodon loop, which through wobbling expand the decoding potential of the tRNAs. Using tRNA-based microarrays, we also determined the dynamics in tRNA abundance upon changes in the growth rate and heterologous protein overexpression stress. With a fourfold increase in the growth rate, the relative abundance of tRNAs cognate to low abundance codons decrease, while the tRNAs cognate to major codons remain mostly unchanged. Significant changes in the tRNA abundances are observed upon protein overexpression stress, which does not correlate with the codon usage of the overexpressed gene but rather reflects the altered expression of housekeeping genes.

Metabolic de-isotoping for improved LC-MS characterization of modified RNAs

Mapping, sequencing, and quantifying individual noncoding ribonucleic acids (ncRNAs), including post-transcriptionally modified nucleosides, by mass spectrometry is a challenge that often requires rigorous sample preparation prior to analysis. Previously, we have described a simplified method for the comparative analysis of RNA digests (CARD) that is applicable to relatively complex mixtures of ncRNAs. In the CARD approach for transfer RNA (tRNA) analysis, two complete sets of digestion products from total tRNA are compared using the enzymatic incorporation of (16)O/(18)O isotopic labels. This approach allows one to rapidly screen total tRNAs from gene deletion mutants or comparatively sequence total tRNA from two related bacterial organisms. However, data analysis can be challenging because of convoluted mass spectra arising from the natural (13)C and (15) N isotopes present in the ribonuclease-digested tRNA samples. Here, we demonstrate that culturing in (12)C-enriched/(13)C-depleted media significantly reduces the isotope patterns that must be interpreted during the CARD experiment. Improvements in data quality yield a 35 % improvement in detection of tRNA digestion products that can be uniquely assigned to particular tRNAs. These mass spectral improvements lead to a significant reduction in data processing attributable to the ease of spectral identification of labeled digestion products and will enable improvements in the relative quantification of modified RNAs by the (16)O/(18)O differential labeling approach.