Identifying modifications in RNA by MALDI mass spectrometry

MazF-rolling circle amplification combined MALDI-TOF MS for site-specific detection of N6-methyladenosine RNA

N6-methyladenosine (m6A) is one of the most abundant chemical modifications in RNA and has vital significance in cellular processes and tumor development. However, the accurate analysis of site-specific m6A modification remains a challenge. In this work, a MazF endoribonuclease activated rolling circle amplification (MazF-RCA) combined MALDI-TOF MS assay is developed for the detection of site-specific m6A-RNA. MazF endoribonuclease can specifically cleave the ACA motif, leaving methylated (m6A)CA motif intact. The intact methylated RNA can then be amplified through rolling circle amplification, and the generated reporter oligonucleotides are detected by MALDI-TOF MS. The assay exhibits good quantification ability, presenting a wide linear range (100 fM to 10 nM) with the limit-of-detection lower than 100 fM. Additionally, the assay can accurately detect methylated RNA in the presence of large amount of non-methylated RNA with a relative abundance of methylated RNA down to 0.5%. The developed assay was further applied to detect m6A-RNA spiked in MCF-7 cell RNA extracts, with the recovery rates in the range of 90.64-106.93%. The present assay provides a novel platform for the analysis of site-specific m6A-RNA at high specificity and sensitivity, which can promote the study of RNA methylation in clinical and biomedical research.

An improved residual network using deep fusion for identifying RNA 5-methylcytosine sites

MOTIVATION

5-Methylcytosine (m5C) is a crucial post-transcriptional modification. With the development of technology, it is widely found in various RNAs. Numerous studies have indicated that m5C plays an essential role in various activities of organisms, such as tRNA recognition, stabilization of RNA structure, RNA metabolism and so on. Traditional identification is costly and time-consuming by wet biological experiments. Therefore, computational models are commonly used to identify the m5C sites. Due to the vast computing advantages of deep learning, it is feasible to construct the predictive model through deep learning algorithms.

RESULTS

In this study, we construct a model to identify m5C based on a deep fusion approach with an improved residual network. First, sequence features are extracted from the RNA sequences using Kmer, K-tuple nucleotide frequency component (KNFC), Pseudo dinucleotide composition (PseDNC) and Physical and chemical property (PCP). Kmer and KNFC extract information from a statistical point of view. PseDNC and PCP extract information from the physicochemical properties of RNA sequences. Then, two parts of information are fused with new features using bidirectional long- and short-term memory and attention mechanisms, respectively. Immediately after, the fused features are fed into the improved residual network for classification. Finally, 10-fold cross-validation and independent set testing are used to verify the credibility of the model. The results show that the accuracy reaches 91.87%, 95.55%, 92.27% and 95.60% on the training sets and independent test sets of Arabidopsis thaliana and M.musculus, respectively. This is a considerable improvement compared to previous studies and demonstrates the robust performance of our model.

AVAILABILITY AND IMPLEMENTATION

The data and code related to the study are available at https://github.com/alivelxj/m5c-DFRESG.

Dihydrouridine in the Transcriptome: New Life for This Ancient RNA Chemical Modification

Until recently, post-transcriptional modifications of RNA were largely restricted to noncoding RNA species. However, this belief seems to have quickly dissipated with the growing number of new modifications found in mRNA that were originally thought to be primarily tRNA-specific, such as dihydrouridine. Recently, transcriptomic profiling, metabolic labeling, and proteomics have identified unexpected dihydrouridylation of mRNAs, greatly expanding the catalog of novel mRNA modifications. These data also implicated dihydrouridylation in meiotic chromosome segregation, protein translation rates, and cell proliferation. Dihydrouridylation of tRNAs and mRNAs are introduced by flavin-dependent dihydrouridine synthases. In this review, we will briefly outline the current knowledge on the distribution of dihydrouridines in the transcriptome, their chemical labeling, and highlight structural and mechanistic aspects regarding the dihydrouridine synthases enzyme family. A special emphasis on important research directions to be addressed will also be discussed. This new entry of dihydrouridine into mRNA modifications has definitely added a new layer of information that controls protein synthesis.

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

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

Characterization and Sequence Mapping of Large RNA and mRNA Therapeutics Using Mass Spectrometry

Large RNA including mRNA (mRNA) has emerged as an important new class of therapeutics. Recently, this has been demonstrated by two highly efficacious vaccines based on mRNA sequences encoding for a modified version of the SARS-CoV-2 spike protein. There is currently significant demand for the development of new and improved analytical methods for the characterization of large RNA including mRNA therapeutics. In this study, we have developed an automated, high-throughput workflow for the rapid characterization and direct sequence mapping of large RNA and mRNA therapeutics. Partial RNase digestions using RNase T1 immobilized on magnetic particles were performed in conjunction with high-resolution liquid chromatography–mass spectrometry analysis. Sequence mapping was performed using automated oligoribonucleotide annotation and identifications based on MS/MS spectra. Using this approach, a >80% sequence of coverage of a range of large RNAs and mRNA therapeutics including the SARS-CoV-2 spike protein was obtained in a single analysis. The analytical workflow, including automated sample preparation, can be completed within 90 min. The ability to rapidly identify, characterize, and sequence map large mRNA therapeutics with high sequence coverage provides important information for identity testing, sequence validation, and impurity analysis.

Intrinsic Strand-Incision Activity of Human UNG: Implications for Nick Generation in Immunoglobulin Gene Diversification

Uracil arises in cellular DNA by cytosine (C) deamination and erroneous replicative incorporation of deoxyuridine monophosphate opposite adenine. The former generates C → thymine transition mutations if uracil is not removed by uracil-DNA glycosylase (UDG) and replaced by C by the base excision repair (BER) pathway. The primary human UDG is hUNG. During immunoglobulin gene diversification in activated B cells, targeted cytosine deamination by activation-induced cytidine deaminase followed by uracil excision by hUNG is important for class switch recombination (CSR) and somatic hypermutation by providing the substrate for DNA double-strand breaks and mutagenesis, respectively. However, considerable uncertainty remains regarding the mechanisms leading to DNA incision following uracil excision: based on the general BER scheme, apurinic/apyrimidinic (AP) endonuclease (APE1 and/or APE2) is believed to generate the strand break by incising the AP site generated by hUNG. We report here that hUNG may incise the DNA backbone subsequent to uracil excision resulting in a 3´-α,β-unsaturated aldehyde designated uracil-DNA incision product (UIP), and a 5´-phosphate. The formation of UIP accords with an elimination (E2) reaction where deprotonation of C2´ occurs via the formation of a C1´ enolate intermediate. UIP is removed from the 3´-end by hAPE1. This shows that the first two steps in uracil BER can be performed by hUNG, which might explain the significant residual CSR activity in cells deficient in APE1 and APE2.

Regulation of translation by site-specific ribosomal RNA methylation

Ribosomes are complex ribozymes that interpret genetic information by translating messenger RNA (mRNA) into proteins. Natural variation in ribosome composition has been documented in several organisms and can arise from several different sources. A key question is whether specific control over ribosome heterogeneity represents a mechanism by which translation can be regulated. We used RiboMeth-seq to demonstrate that differential 2'-O-methylation of ribosomal RNA (rRNA) represents a considerable source of ribosome heterogeneity in human cells, and that modification levels at distinct sites can change dynamically in response to upstream signaling pathways, such as MYC oncogene expression. Ablation of one prominent methylation resulted in altered translation of select mRNAs and corresponding changes in cellular phenotypes. Thus, differential rRNA 2'-O-methylation can give rise to ribosomes with specialized function. This suggests a broader mechanism where the specific regulation of rRNA modification patterns fine tunes translation.

The hyperthermophilic partners Nanoarchaeum and Ignicoccus stabilize their tRNA T-loops via different but structurally equivalent modifications

The universal L-shaped tertiary structure of tRNAs is maintained with the help of nucleotide modifications within the D- and T-loops, and these modifications are most extensive within hyperthermophilic species. The obligate-commensal Nanoarchaeum equitans and its phylogenetically-distinct host Ignicoccus hospitalis grow physically coupled under identical hyperthermic conditions. We report here two fundamentally different routes by which these archaea modify the key conserved nucleotide U54 within their tRNA T-loops. In N. equitans, this nucleotide is methylated by the S-adenosylmethionine-dependent enzyme NEQ053 to form m⁵U54, and a recombinant version of this enzyme maintains specificity for U54 in Escherichia coli. In N. equitans, m⁵U54 is subsequently thiolated to form m⁵s²U54. In contrast, I. hospitalis isomerizes U54 to pseudouridine prior to methylating its N1-position and thiolating the O4-position of the nucleobase to form the previously uncharacterized nucleotide m¹s⁴Ψ. The methyl and thiol groups in m¹s⁴Ψ and m⁵s²U are presented within the T-loop in a spatially identical manner that stabilizes the 3′-endo-anti conformation of nucleotide-54, facilitating stacking onto adjacent nucleotides and reverse-Hoogsteen pairing with nucleotide m¹A58. Thus, two distinct structurally-equivalent solutions have evolved independently and convergently to maintain the tertiary fold of tRNAs under extreme hyperthermic conditions.

Unraveling the RNA modification code with mass spectrometry

The discovery and analysis of modifications on proteins and nucleic acids has provided functional information that has rapidly accelerated the field of epigenetics. While protein post-translational modifications (PTMs), especially on histones, have been highlighted as critical components of epigenetics, the post-transcriptional modification of RNA has been a subject of more recently emergent interest. Multiple RNA modifications have been known to be present in tRNA and rRNA since the 1960s, but the exploration of mRNA, small RNA, and inducible tRNA modifications remains nascent. Sequencing-based methods have been essential to the field by creating the first epitranscriptome maps of m6A, m5C, hm5C, pseudouridine, and inosine; however, these methods possess significant limitations. Here, we discuss the past, present, and future of the application of mass spectrometry (MS) to the study of RNA modifications.

Probing the diversity and regulation of tRNA modifications

Transfer RNAs (tRNAs) are non-coding RNAs essential for protein synthesis. tRNAs are heavily decorated with a variety of post-transcriptional modifications (tRNA modifications). Recent methodological advances provide new tools for rapid profiling of tRNA modifications and have led to discoveries of novel modifications and their regulation. Here, we provide an overview of the techniques for investigating tRNA modifications and of the expanding knowledge of their chemistry and regulation.

Detection of internal N7-methylguanosine (m7G) RNA modifications by mutational profiling sequencing

Methylation of guanosine on position N7 (m⁷G) on internal RNA positions has been found in all domains of life and have been implicated in human disease. Here, we present m⁷G Mutational Profiling sequencing (m⁷G-MaP-seq), which allows high throughput detection of m⁷G modifications at nucleotide resolution. In our method, m⁷G modified positions are converted to abasic sites by reduction with sodium borohydride, directly recorded as cDNA mutations through reverse transcription and sequenced. We detect positions with increased mutation rates in the reduced and control samples taking the possibility of sequencing/alignment error into account and use replicates to calculate statistical significance based on log likelihood ratio tests. We show that m⁷G-MaP-seq efficiently detects known m⁷G modifications in rRNA with mutational rates up to 25% and we map a previously uncharacterised evolutionarily conserved rRNA modification at position 1581 in Arabidopsis thaliana SSU rRNA. Furthermore, we identify m⁷G modifications in budding yeast, human and arabidopsis tRNAs and demonstrate that m⁷G modification occurs before tRNA splicing. We do not find any evidence for internal m⁷G modifications being present in other small RNA, such as miRNA, snoRNA and sRNA, including human Let-7e. Likewise, high sequencing depth m⁷G-MaP-seq analysis of mRNA from E. coli or yeast cells did not identify any internal m⁷G modifications.

Reductive Evolution and Diversification of C5-Uracil Methylation in the Nucleic Acids of Mollicutes

The C5-methylation of uracil to form 5-methyluracil (m5U) is a ubiquitous base modification of nucleic acids. Four enzyme families have converged to catalyze this methylation using different chemical solutions. Here, we investigate the evolution of 5-methyluracil synthase families in Mollicutes, a class of bacteria that has undergone extensive genome erosion. Many mollicutes have lost some of the m5U methyltransferases present in their common ancestor. Cases of duplication and subsequent shift of function are also described. For example, most members of the Spiroplasma subgroup use the ancestral tetrahydrofolate-dependent TrmFO enzyme to catalyze the formation of m5U54 in tRNA, while a TrmFO paralog (termed RlmFO) is responsible for m5U1939 formation in 23S rRNA. RlmFO has replaced the S-adenosyl-L-methionine (SAM)-enzyme RlmD that adds the same modification in the ancestor and which is still present in mollicutes from the Hominis subgroup. Another paralog of this family, the TrmFO-like protein, has a yet unidentified function that differs from the TrmFO and RlmFO homologs. Despite having evolved towards minimal genomes, the mollicutes possess a repertoire of m5U-modifying enzymes that is highly dynamic and has undergone horizontal transfer.

The human methyltransferase ZCCHC4 catalyses N6-methyladenosine modification of 28S ribosomal RNA

RNA methylations are essential both for RNA structure and function, and are introduced by a number of distinct methyltransferases (MTases). In recent years, N⁶-methyladenosine (m⁶A) modification of eukaryotic mRNA has been subject to intense studies, and it has been demonstrated that m⁶A is a reversible modification that regulates several aspects of mRNA function. However, m⁶A is also found in other RNAs, such as mammalian 18S and 28S ribosomal RNAs (rRNAs), but the responsible MTases have remained elusive. 28S rRNA carries a single m⁶A modification, found at position A4220 (alternatively referred to as A4190) within a stem–loop structure, and here we show that the MTase ZCCHC4 is the enzyme responsible for introducing this modification. Accordingly, we found that ZCCHC4 localises to nucleoli, the site of ribosome assembly, and that proteins involved in RNA metabolism are overrepresented in the ZCCHC4 interactome. Interestingly, the absence of m⁶A4220 perturbs codon-specific translation dynamics and shifts gene expression at the translational level. In summary, we establish ZCCHC4 as the enzyme responsible for m⁶A modification of human 28S rRNA, and demonstrate its functional significance in mRNA translation.

Excision of uracil from DNA by hSMUG1 includes strand incision and processing

Uracil arises in DNA by hydrolytic deamination of cytosine (C) and by erroneous incorporation of deoxyuridine monophosphate opposite adenine, where the former event is devastating by generation of C → thymine transitions. The base excision repair (BER) pathway replaces uracil by the correct base. In human cells two uracil-DNA glycosylases (UDGs) initiate BER by excising uracil from DNA; one is hSMUG1 (human single-strand-selective mono-functional UDG). We report that repair initiation by hSMUG1 involves strand incision at the uracil site resulting in a 3′-α,β-unsaturated aldehyde designated uracil-DNA incision product (UIP), and a 5′-phosphate. UIP is removed from the 3′-end by human apurinic/apyrimidinic (AP) endonuclease 1 preparing for single-nucleotide insertion. hSMUG1 also incises DNA or processes UIP to a 3′-phosphate designated uracil-DNA processing product (UPP). UIP and UPP were indirectly identified and quantified by polyacrylamide gel electrophoresis and chemically characterised by matrix-assisted laser desorption/ionisation time-of-flight mass-spectrometric analysis of DNA from enzyme reactions using ¹⁸O- or ¹⁶O-water. The formation of UIP accords with an elimination (E2) reaction where deprotonation of C2′ occurs via the formation of a C1′ enolate intermediate. A three-phase kinetic model explains rapid uracil excision in phase 1, slow unspecific enzyme adsorption/desorption to DNA in phase 2 and enzyme-dependent AP site incision in phase 3.

RNA 2'-O-Methylation (Nm) Modification in Human Diseases

Nm (2'-O-methylation) is one of the most common modifications in the RNA world. It has the potential to influence the RNA molecules in multiple ways, such as structure, stability, and interactions, and to play a role in various cellular processes from epigenetic gene regulation, through translation to self versus non-self recognition. Yet, building scientific knowledge on the Nm matter has been hampered for a long time by the challenges in detecting and mapping this modification. Today, with the latest advancements in the area, more and more Nm sites are discovered on RNAs (tRNA, rRNA, mRNA, and small non-coding RNA) and linked to normal or pathological conditions. This review aims to synthesize the Nm-associated human diseases known to date and to tackle potential indirect links to some other biological defects.

High resolution fingerprinting of single and double-stranded RNA using ion-pair reverse-phase chromatography

The emergence of new sustainable approaches for insect management using RNA interference (RNAi) based insecticides has created the demand for high throughput analytical techniques to fully characterise and accurately quantify double stranded RNA (dsRNA) prior to downstream RNAi applications. In this study we have developed a method for the rapid characterisation of single stranded and double stranded RNA using high resolution RNase mapping in conjunction with ion-pair reverse-phase chromatography utilising a column with superficially porous particles. The high resolution oligoribonucleotide map provides an important 'fingerprint' for identity testing and bioprocess monitoring. Reproducible RNA mapping chromatograms were generated from replicate analyses. Moreover, this approach was used to provide a method to rapidly distinguish different RNA sequences of the same size, based on differences in the resulting chromatograms. Principal components analysis of the high resolution RNA mapping data enabled us to rapidly compare multiple HPLC chromatograms and distinguish two dsRNA sequences of different size which share 72% sequence homology. We used the high resolution RNase mapping method to rapidly fingerprint biomanufactured dsRNA across a number of different batches. The resulting chromatograms in conjunction with principal components analysis demonstrated high similarity in the dsRNA produced across the different batches highlighting the potential ability of this method to provide information for batch release in a high throughput manner.

Mapping of ribosomal 23S ribosomal RNA modifications in Clostridium sporogenes

All organisms contain RNA modifications in their ribosomal RNA (rRNA), but the importance, positions and exact function of these are still not fully elucidated. Various functions such as stabilizing structures, controlling ribosome assembly and facilitating interactions have been suggested and in some cases substantiated. Bacterial rRNA contains much fewer modifications than eukaryotic rRNA. The rRNA modification patterns in bacteria differ from each other, but too few organisms have been mapped to draw general conclusions. This study maps 23S ribosomal RNA modifications in Clostridium sporogenes that can be characterized as a non-toxin producing Clostridium botulinum. Clostridia are able to sporulate and thereby survive harsh conditions, and are in general considered to be resilient to antibiotics. Selected regions of the 23S rRNA were investigated by mass spectrometry and by primer extension analysis to pinpoint modified sites and the nature of the modifications. Apparently, C. sporogenes 23S rRNA contains few modifications compared to other investigated bacteria. No modifications were identified in domain II and III of 23S rRNA. Three modifications were identified in domain IV, all of which have also been found in other organisms. Two unusual modifications were identified in domain V, methylated dihydrouridine at position U2449 and dihydrouridine at position U2500 (Escherichia coli numbering), in addition to four previously known modified positions. The enzymes responsible for the modifications were searched for in the C. sporogenes genome using BLAST with characterized enzymes as query. The search identified genes potentially coding for RNA modifying enzymes responsible for most of the found modifications.

Biofilm Formation and Motility Are Promoted by Cj0588-Directed Methylation of rRNA in Campylobacter jejuni

Numerous bacterial pathogens express an ortholog of the enzyme TlyA, which is an rRNA 2′-O-methyltransferase associated with resistance to cyclic peptide antibiotics such as capreomycin. Several other virulence traits have also been attributed to TlyA, and these appear to be unrelated to its methyltransferase activity. The bacterial pathogen Campylobacter jejuni possesses the TlyA homolog Cj0588, which has been shown to contribute to virulence. Here, we investigate the mechanism of Cj0588 action and demonstrate that it is a type I homolog of TlyA that 2′-O-methylates 23S rRNA nucleotide C1920. This same specific function is retained by Cj0588 both in vitro and also when expressed in Escherichia coli. Deletion of the cj0588 gene in C. jejuni or substitution with alanine of K⁸⁰, D¹⁶², or K¹⁸⁸ in the catalytic center of the enzyme cause complete loss of 2′-O-methylation activity. Cofactor interactions remain unchanged and binding affinity to the ribosomal substrate is only slightly reduced, indicating that the inactivated proteins are folded correctly. The substitution mutations thus dissociate the 2′-O-methylation function of Cj0588/TlyA from any other putative roles that the protein might play. C. jejuni strains expressing catalytically inactive versions of Cj0588 have the same phenotype as cj0588-null mutants, and show altered tolerance to capreomycin due to perturbed ribosomal subunit association, reduced motility and impaired ability to form biofilms. These functions are reestablished when methyltransferase activity is restored and we conclude that the contribution of Cj0588 to virulence in C. jejuni is a consequence of the enzyme's ability to methylate its rRNA.

Bypassing rRNA methylation by RsmA/Dim1during ribosome maturation in the hyperthermophilic archaeon Nanoarchaeum equitans

In all free-living organisms a late-stage checkpoint in the biogenesis of the small ribosomal subunit involves rRNA modification by an RsmA/Dim1 methyltransferase. The hyperthermophilic archaeon Nanoarchaeum equitans, whose existence is confined to the surface of a second archaeon, Ignicoccus hospitalis, lacks an RsmA/Dim1 homolog. We demonstrate here that the I. hospitalis host possesses the homolog Igni_1059, which dimethylates the N6-positions of two invariant adenosines within helix 45 of 16S rRNA in a manner identical to other RsmA/Dim1 enzymes. However, Igni_1059 is not transferred from I. hospitalis to N. equitans across their fused cell membrane structures and the corresponding nucleotides in N. equitans 16S rRNA remain unmethylated. An alternative mechanism for ribosomal subunit maturation in N. equitans is suggested by sRNA interactions that span the redundant RsmA/Dim1 site to introduce 2΄-O-ribose methylations within helices 44 and 45 of the rRNA.

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.

Profiling of 2'-O-Me in human rRNA reveals a subset of fractionally modified positions and provides evidence for ribosome heterogeneity

Ribose methylation is one of the two most abundant modifications in human ribosomal RNA and is believed to be important for ribosome biogenesis, mRNA selectivity and translational fidelity. We have applied RiboMeth-seq to rRNA from HeLa cells for ribosome-wide, quantitative mapping of 2'-O-Me sites and obtained a comprehensive set of 106 sites, including two novel sites, and with plausible box C/D guide RNAs assigned to all but three sites. We find approximately two-thirds of the sites to be fully methylated and the remainder to be fractionally modified in support of ribosome heterogeneity at the level of RNA modifications. A comparison to HCT116 cells reveals similar 2'-O-Me profiles with distinct differences at several sites. This study constitutes the first comprehensive mapping of 2'-O-Me sites in human rRNA using a high throughput sequencing approach. It establishes the existence of a core of constitutively methylated positions and a subset of variable, potentially regulatory positions, and paves the way for experimental analyses of the role of variations in rRNA methylation under different physiological or pathological settings.

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.

Identification of the methyltransferase targeting C2499 in Deinococcus radiodurans 23S ribosomal RNA

The bacterium Deinococcus radiodurans-like all other organisms-introduces nucleotide modifications into its ribosomal RNA. We have previously found that the bacterium contains a Carbon-5 methylation on cytidine 2499 of its 23S ribosomal RNA, which is so far the only modified version of cytidine 2499 reported. Using homology search, we identified the open reading frame DR_0049 as the primary candidate gene for the methyltransferase that modifies cytidine 2499. Mass spectrometric analysis demonstrated that recombinantly expressed DR0049 protein methylates E. coli cytidine 2499 both in vitro and in vivo. We also inactivated the DR_0049 gene in D. radiodurans through insertion of a chloramphenicol resistance cassette. This resulted in complete absence of the cytidine 2499 methylation, which all together demonstrates that DR_0049 encodes the methyltransferase producing m(5)C2499 in D. radiodurans 23S rRNA. Growth experiments disclosed that inactivation of DR_0049 is associated with a severe growth defect, but available ribosome structures show that cytidine 2499 is positioned very similar in D. radiodurans harbouring the modification and E. coli without the modification. Hence there is no obvious structure-based explanation for the requirement for the C2499 posttranscriptional modification in D. radiodurans.

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.

Resistance to ketolide antibiotics by coordinated expression of rRNA methyltransferases in a bacterial producer of natural ketolides

Ketolides are promising new antimicrobials effective against a broad range of Gram-positive pathogens, in part because of the low propensity of these drugs to trigger the expression of resistance genes. A natural ketolide pikromycin and a related compound methymycin are produced by Streptomyces venezuelae strain ATCC 15439. The producer avoids the inhibitory effects of its own antibiotics by expressing two paralogous rRNA methylase genes pikR1 and pikR2 with seemingly redundant functions. We show here that the PikR1 and PikR2 enzymes mono- and dimethylate, respectively, the N6 amino group in 23S rRNA nucleotide A2058. PikR1 monomethylase is constitutively expressed; it confers low resistance at low fitness cost and is required for ketolide-induced activation of pikR2 to attain high-level resistance. The regulatory mechanism controlling pikR2 expression has been evolutionary optimized for preferential activation by ketolide antibiotics. The resistance genes and the induction mechanism remain fully functional when transferred to heterologous bacterial hosts. The anticipated wide use of ketolide antibiotics could promote horizontal transfer of these highly efficient resistance genes to pathogens. Taken together, these findings emphasized the need for surveillance of pikR1/pikR2-based bacterial resistance and the preemptive development of drugs that can remain effective against the ketolide-specific resistance mechanism.

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.

On the mechanism of RNA phosphodiester backbone cleavage in the absence of solvent

Ribonucleic acid (RNA) modifications play an important role in the regulation of gene expression and the development of RNA-based therapeutics, but their identification, localization and relative quantitation by conventional biochemical methods can be quite challenging. As a promising alternative, mass spectrometry (MS) based approaches that involve RNA dissociation in ‘top-down’ strategies are currently being developed. For this purpose, it is essential to understand the dissociation mechanisms of unmodified and posttranscriptionally or synthetically modified RNA. Here, we have studied the effect of select nucleobase, ribose and backbone modifications on phosphodiester bond cleavage in collisionally activated dissociation (CAD) of positively and negatively charged RNA. We found that CAD of RNA is a stepwise reaction that is facilitated by, but does not require, the presence of positive charge. Preferred backbone cleavage next to adenosine and guanosine in CAD of (M+nH)ⁿ⁺ and (M−nH)ⁿ⁻ ions, respectively, is based on hydrogen bonding between nucleobase and phosphodiester moieties. Moreover, CAD of RNA involves an intermediate that is sufficiently stable to survive extension of the RNA structure and intramolecular proton redistribution according to simple Coulombic repulsion prior to backbone cleavage into c and y ions from phosphodiester bond cleavage.

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.

Dynamics of RNA modification by a multi-site-specific tRNA methyltransferase

In most organisms, the widely conserved 1-methyl-adenosine58 (m1A58) tRNA modification is catalyzed by an S-adenosyl-L-methionine (SAM)-dependent, site-specific enzyme TrmI. In archaea, TrmI also methylates the adjacent adenine 57, m1A57 being an obligatory intermediate of 1-methyl-inosine57 formation. To study this multi-site specificity, we used three oligoribonucleotide substrates of Pyrococcus abyssi TrmI (PabTrmI) containing a fluorescent 2-aminopurine (2-AP) at the two target positions and followed the RNA binding kinetics and methylation reactions by stopped-flow and mass spectrometry. PabTrmI did not modify 2-AP but methylated the adjacent target adenine. 2-AP seriously impaired the methylation of A57 but not A58, confirming that PabTrmI methylates efficiently the first adenine of the A57A58A59 sequence. PabTrmI binding provoked a rapid increase of fluorescence, attributed to base unstacking in the environment of 2-AP. Then, a slow decrease was observed only with 2-AP at position 57 and SAM, suggesting that m1A58 formation triggers RNA release. A model of the protein-tRNA complex shows both target adenines in proximity of SAM and emphasizes no major tRNA conformational change except base flipping during the reaction. The solvent accessibility of the SAM pocket is not affected by the tRNA, thereby enabling S-adenosyl-L-homocysteine to be replaced by SAM without prior release of monomethylated tRNA.

The flavoprotein Mcap0476 (RlmFO) catalyzes m5U1939 modification in Mycoplasma capricolum 23S rRNA

Efficient protein synthesis in all organisms requires the post-transcriptional methylation of specific ribosomal ribonucleic acid (rRNA) and transfer RNA (tRNA) nucleotides. The methylation reactions are almost invariably catalyzed by enzymes that use S-adenosylmethionine (AdoMet) as the methyl group donor. One noteworthy exception is seen in some bacteria, where the conserved tRNA methylation at m5U54 is added by the enzyme TrmFO using flavin adenine dinucleotide together with N5,N10-methylenetetrahydrofolate as the one-carbon donor. The minimalist bacterium Mycoplasma capricolum possesses two homologs of trmFO, but surprisingly lacks the m5U54 tRNA modification. We created single and dual deletions of the trmFO homologs using a novel synthetic biology approach. Subsequent analysis of the M. capricolum RNAs by mass spectrometry shows that the TrmFO homolog encoded by Mcap0476 specifically modifies m5U1939 in 23S rRNA, a conserved methylation catalyzed by AdoMet-dependent enzymes in all other characterized bacteria. The Mcap0476 methyltransferase (renamed RlmFO) represents the first folate-dependent flavoprotein seen to modify ribosomal RNA.

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.

Indigenous and acquired modifications in the aminoglycoside binding sites of Pseudomonas aeruginosa rRNAs

Aminoglycoside antibiotics remain the drugs of choice for treatment of Pseudomonas aeruginosa infections, particularly for respiratory complications in cystic-fibrosis patients. Previous studies on other bacteria have shown that aminoglycosides have their primary target within the decoding region of 16S rRNA helix 44 with a secondary target in 23S rRNA helix 69. Here, we have mapped P. aeruginosa rRNAs using MALDI mass spectrometry and reverse transcriptase primer extension to identify nucleotide modifications that could influence aminoglycoside interactions. Helices 44 and 45 contain indigenous (housekeeping) modifications at m (4)Cm1402, m (3)U1498, m (2)G1516, m (6) 2A1518, and m (6) 2A1519; helix 69 is modified at m (3)Ψ1915, with m (5)U1939 and m (5)C1962 modification in adjacent sequences. All modifications were close to stoichiometric, with the exception of m (3)Ψ1915, where about 80% of rRNA molecules were methylated. The modification status of a virulent clinical strain expressing the acquired methyltransferase RmtD was altered in two important respects: RmtD stoichiometrically modified m (7)G1405 conferring high resistance to the aminoglycoside tobramycin and, in doing so, impeded one of the methylation reactions at C1402. Mapping the nucleotide methylations in P. aeruginosa rRNAs is an essential step toward understanding the architecture of the aminoglycoside binding sites and the rational design of improved drugs against this bacterial pathogen.

Transcriptome-wide mapping of 5-methylcytidine RNA modifications in bacteria, archaea, and yeast reveals m5C within archaeal mRNAs

The presence of 5-methylcytidine (m⁵C) in tRNA and rRNA molecules of a wide variety of organisms was first observed more than 40 years ago. However, detection of this modification was limited to specific, abundant, RNA species, due to the usage of low-throughput methods. To obtain a high resolution, systematic, and comprehensive transcriptome-wide overview of m⁵C across the three domains of life, we used bisulfite treatment on total RNA from both gram positive (B. subtilis) and gram negative (E. coli) bacteria, an archaeon (S. solfataricus) and a eukaryote (S. cerevisiae), followed by massively parallel sequencing. We were able to recover most previously documented m⁵C sites on rRNA in the four organisms, and identified several novel sites in yeast and archaeal rRNAs. Our analyses also allowed quantification of methylated m⁵C positions in 64 tRNAs in yeast and archaea, revealing stoichiometric differences between the methylation patterns of these organisms. Molecules of tRNAs in which m⁵C was absent were also discovered. Intriguingly, we detected m⁵C sites within archaeal mRNAs, and identified a consensus motif of AUCGANGU that directs methylation in S. solfataricus. Our results, which were validated using m⁵C-specific RNA immunoprecipitation, provide the first evidence for mRNA modifications in archaea, suggesting that this mode of post-transcriptional regulation extends beyond the eukaryotic domain.

Ribonucleic acids are universally used to express genetic information in the form of gene transcripts. Although we envision RNA as a mere copy of the DNA four-base code, modification of specific RNA bases can expand the information code. Such modifications are abundant in transfer RNA (tRNA) and ribosomal RNA (rRNA), where they contribute to translation fidelity and ribosome assembly. Recent studies in eukaryotes have shown that mRNA modifications such as RNA-editing (conversion of an adenosine base to inosine), N6-adenine methylation (m⁶A), and 5-methylcytidine (m⁵C) can change the coding sequence, alter splicing patterns, or change RNA stability. However, no mRNA modifications in bacteria or archaea have been documented to date. We have used an approach that enables mapping of the m⁵C modifications across all expressed genes in a given organism. Applying this approach on model bacterial, archaeal, and fungal microorganisms enabled us to reveal the modified RNA bases in these organisms, and to provide an accurate and sensitive map of these modifications. In archaea, we documented multiple genes whose mRNAs are subject to RNA modification, suggesting that similar to eukaryotes, these organisms may utilize mRNA modifications as a mechanism for gene regulation.

Distinction between the Cfr methyltransferase conferring antibiotic resistance and the housekeeping RlmN methyltransferase

The cfr gene encodes the Cfr methyltransferase that primarily methylates C-8 in A2503 of 23S rRNA in the peptidyl transferase region of bacterial ribosomes. The methylation provides resistance to six classes of antibiotics of clinical and veterinary importance. The rlmN gene encodes the RlmN methyltransferase that methylates C-2 in A2503 in 23S rRNA and A37 in tRNA, but RlmN does not significantly influence antibiotic resistance. The enzymes are homologous and use the same mechanism involving radical S-adenosyl methionine to methylate RNA via an intermediate involving a methylated cysteine in the enzyme and a transient cross-linking to the RNA, but they differ in which carbon atom in the adenine they methylate. Comparative sequence analysis identifies differentially conserved residues that indicate functional sequence divergence between the two classes of Cfr- and RlmN-like sequences. The differentiation between the two classes is supported by previous and new experimental evidence from antibiotic resistance, primer extensions, and mass spectrometry. Finally, evolutionary aspects of the distribution of Cfr- and RlmN-like enzymes are discussed.

Association of the novel aminoglycoside resistance determinant RmtF with NDM carbapenemase in Enterobacteriaceae isolated in India and the UK

OBJECTIVES

16S rRNA methyltransferases are an emerging mechanism conferring high-level resistance to clinically relevant aminoglycosides and have been associated with important mechanisms such as NDM-1. We sought genes encoding these enzymes in isolates highly resistant (MIC >200 mg/L) to gentamicin and amikacin from an Indian hospital and we additionally screened for the novel RmtF enzyme in 132 UK isolates containing NDM.

METHODS

All highly aminoglycoside-resistant isolates were screened for armA and rmtA-E by PCR, with cloning experiments performed for isolates negative for these genes. Matrix-assisted laser desorption/ionization (MALDI) mass spectrometry was used to determine the methylation target of the novel RmtF methyltransferase. RmtF-bearing strains were characterized further, including susceptibility testing, PFGE, electroporation, PCR-based replicon typing and multilocus sequence typing of rmtF-bearing plasmids.

RESULTS

High-level aminoglycoside resistance was detected in 140/1000 (14%) consecutive isolates of Enterobacteriaceae from India. ArmA, RmtB and RmtC were identified among 46%, 20% and 27% of these isolates, respectively. The novel rmtF gene was detected in 34 aminoglycoside-resistant isolates (overall prevalence 3.4%), most (59%) of which also possessed a bla(NDM) gene; rmtF was detected in 6 NDM producers from the UK. It was found on different plasmid backbones. Four and two isolates showed resistance to tigecycline and colistin, respectively.

CONCLUSIONS

RmtF was often found in association with NDM in members of the Enterobacteriaceae and on diverse plasmids. It is of clinical concern that the RmtF- and NDM-positive strains identified here show additional resistance to tigecycline and colistin, current drugs of last resort for the treatment of serious bacterial infections.

Method for comparative analysis of ribonucleic acids using isotope labeling and mass spectrometry

Here, we describe a method for the comparative analysis of ribonucleic acids (RNAs). This method allows sequence or modification information from a previously uncharacterized RNA to be obtained by direct comparison with a reference RNA, whose sequence or modification information is known. This simple and rapid method is enabled by the differential labeling of two RNA samples. One sample, the reference RNA, is labeled with (16)O during enzymatic digestion. The second sample, the candidate or unknown RNA, is labeled with (18)O. By combining the two digests, digestion products that share the same sequence or post-transcriptional modification(s) between the reference and candidate will appear as doublets separated by 2 Da. Sequence or modification differences between the two will generate singlets that can be further characterized to identify how the candidate sequence differs from the reference. We illustrate the application of this approach for sequencing individual RNAs and demonstrate how this method can be used to identify sequence-specific differences in RNA modification. This comparative analysis of RNA digests (CARD) approach is scalable to multiple candidate RNAs using one or multiple reference RNAs and is compatible with existing methods for quantitative analysis of RNAs.

The human tRNA m (5) C methyltransferase Misu is multisite-specific

The human tRNA m ( 5) C methyltransferase Misu is a novel downstream target of the proto-oncogene Myc that participates in controlling cell division and proliferation. Misu catalyzes the transfer of a methyl group from S-adenosyl-L-methionine to carbon 5 of cytosines in tRNAs. It was previously shown to catalyze in vitro the intron-dependent formation of m ( 5) C at the first position of the anticodon (position 34) within the human pre-tRNA (Leu) (CAA). In addition, it was recently reported that C48 and C49 are methylated in vivo by Misu. We report here the expression of hMisu in Escherichia coli and its purification to homogeneity. We show that this enzyme methylates position 48 in tRNA (Leu) (CAA) with or without intron and positions 48, 49 and 50 in tRNA (Gly2) (GCC) in vitro. Therefore, hMisu is the enzyme responsible for the methylation of at least four cytosines in human tRNAs. By comparison, the orthologous yeast enzyme Trm4 catalyzes the methylation of carbon 5 of cytosine at positions 34, 40, 48 or 49 depending on the tRNAs.

Recognition of guanosine by dissimilar tRNA methyltransferases

Guanosines are important for biological activities through their specific functional groups that are recognized for RNA or protein interactions. One example is recognition of N(1) of G37 in tRNA by S-adenosyl-methionine (AdoMet)-dependent tRNA methyltransferases to synthesize m(1)G37-tRNA, which is essential for translational fidelity in all biological domains. Synthesis of m(1)G37-tRNA is catalyzed by TrmD in bacteria and by Trm5 in eukarya and archaea, using unrelated and dissimilar structural folds. This raises the question of how dissimilar proteins recognize the same guanosine. Here we probe the mechanism of discrimination among functional groups of guanosine by TrmD and Trm5. Guanosine analogs were systematically introduced into tRNA through a combination of chemical and enzymatic synthesis. Single turnover kinetic assays and thermodynamic analysis of the effect of each analog on m(1)G37-tRNA synthesis reveal that TrmD and Trm5 discriminate functional groups differently. While both recognize N(1) and O(6) of G37, TrmD places a much stronger emphasis on these functional groups than Trm5. While the exocyclic 2-amino group of G37 is important for TrmD, it is dispensable for Trm5. In addition, while an adjacent G36 is obligatory for TrmD, it is nonessential for Trm5. These results depict a more rigid requirement of guanosine functional groups for TrmD than for Trm5. However, the sensitivity of both enzymes to analog substitutions, together with an experimental revelation of their low cellular concentrations relative to tRNA substrates, suggests a model in which these enzymes rapidly screen tRNA by direct recognition of G37 in order to monitor the global state of m(1)G37-tRNA.

RTL-P: a sensitive approach for detecting sites of 2'-O-methylation in RNA molecules

2′-O-methylation is present within various cellular RNAs and is essential to RNA biogenesis and functionality. Several methods have been developed for the identification and localization of 2′-O-methylated sites in RNAs; however, the detection of RNA modifications, especially in low-abundance RNAs and small non-coding RNAs with a 2′-O-methylation at the 3′-end, remains a difficult task. Here, we introduce a new method to detect 2′-O-methylated sites in diverse RNA species, referred to as RTL-P [Reverse Transcription at Low deoxy-ribonucleoside triphosphate (dNTP) concentrations followed by polymerase chain reaction (PCR)] that demonstrates precise mapping and superior sensitivity compared with previous techniques. The main procedures of RTL-P include a site-specific primer extension by reverse transcriptase at a low dNTP concentration and a semi-quantitative PCR amplification step. No radiolabeled or fluorescent primers are required. By designing specific RT primers, we used RTL-P to detect both previously identified and novel 2′-O-methylated sites in human and yeast ribosomal RNAs (rRNAs), as well as mouse piwi-interacting RNAs (piRNAs). These results demonstrate the powerful application of RTL-P for the systematic analysis of fully or partially methylated residues in diverse RNA species, including low-abundance RNAs or small non-coding RNAs such as piRNAs and microRNAs (miRNAs).

Identification and characterization of the Thermus thermophilus 5-methylcytidine (m5C) methyltransferase modifying 23 S ribosomal RNA (rRNA) base C1942

Methylation of cytidines at carbon-5 is a common posttranscriptional RNA modification encountered across all domains of life. Here, we characterize the modifications of C1942 and C1962 in Thermus thermophilus 23 S rRNA as 5-methylcytidines (m(5)C) and identify the two associated methyltransferases. The methyltransferase modifying C1942, named RlmO, has not been characterized previously. RlmO modifies naked 23 S rRNA, but not the assembled 50 S subunit or 70 S ribosomes. The x-ray crystal structure of this enzyme in complex with the S-adenosyl-l-methionine cofactor at 1.7 Å resolution confirms that RlmO is structurally related to other m(5)C rRNA methyltransferases. Key residues in the active site are located similar to the further distant 5-methyluridine methyltransferase RlmD, suggestive of a similar enzymatic mechanism. RlmO homologues are primarily found in mesophilic bacteria related to T. thermophilus. In accordance, we find that growth of the T. thermophilus strain with an inactivated C1942 methyltransferase gene is not compromised at non-optimal temperatures.

Removal of 3'-phosphate group by bacterial alkaline phosphatase improves oligonucleotide sequence coverage of RNase digestion products analyzed by collision-induced dissociation mass spectrometry

RNase mapping by nucleobase-specific endonucleases combined with liquid chromatography/tandem mass spectrometry (LC/MS/MS) is a powerful analytical method for characterizing ribonucleic acids (RNAs). Endonuclease digestion of RNA yields products that contain a 3'-terminal phosphate group. MS/MS via collision-induced dissociation (CID) of these digestion products on a linear ion trap generates fragmentation pathways that include the loss of phosphoric acid (-H(3)PO(4); -98 u), which does not provide information about the sequence of the digestion products and can reduce ion abundance from other pathways that provide sequence information. Here we investigate the use of bacterial alkaline phosphatase (BAP) after RNase digestion to remove the 3'-terminal phosphate from all RNase digestion products prior to LC/MS/MS analysis. RNase digestion products lacking the 3'-phosphate were found to produce CID spectra with more consistent, high-abundance c- and y-type fragment ions as well as significantly more a-Base and w-type ions than digestion products retaining the 3'-phosphate. In this manner, RNase mapping with LC/MS/MS can provide more complete RNA sequence information from fragment ions of higher abundance that are easier to interpret and identify.

3-(3-amino-3-carboxypropyl)-5,6-dihydrouridine is one of two novel post-transcriptional modifications in tRNALys(UUU) from Trypanosoma brucei

tRNA is the most heavily modified of all RNA types, with typically 10-20% of the residues being post-transcriptionally altered. Unravelling the modification pattern of a tRNA is a challenging task; there are 92 currently known tRNA modifications, many of which are chemically similar. Furthermore, the tRNA has to be investigated with single-nucleotide resolution in order to ensure complete mapping of all modifications. In the present work, we characterized tRNA(Lys)(UUU) from Trypanosoma brucei, and provide a complete overview of its post-transcriptional modifications. The first step was MALDI-TOF MS of two independent digests of the tRNA, with RNase A and RNase T1, respectively. This revealed digestion products harbouring mass-changing modifications. Next, the modifications were mapped at the nucleotide level in the RNase products by tandem MS. Comparison with the sequence of the unmodified tRNA revealed the modified residues. The modifications were further characterized at the nucleoside level by chromatographic retention time and fragmentation pattern upon higher-order tandem MS. Phylogenetic comparison with modifications in tRNA(Lys) from other organisms was used through the entire analysis. We identified modifications on 12 nucleosides in tRNA(Lys)(UUU), where U47 exhibited a novel modification, 3-(3-amino-3-carboxypropyl)-5,6-dihydrouridine, based on identical chromatographic retention and MS fragmentation as the synthetic nucleoside. A37 was observed in two versions: a minor fraction with the previously described 2-methylthio-N(6)-threonylcarbamoyl-modification, and a major fraction with A37 being modified by a 294.0-Da moiety. The latter product is the largest adenosine modification reported so far, and we discuss its nature and origin.