Systematic identification of tRNAome and its dynamics in Lactococcus lactis

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.

Alternate routes to mnm5s2U synthesis in Gram-positive bacteria

The wobble bases of tRNAs that decode split codons are often heavily modified. In bacteria, tRNAGlu, Gln, Asp contains a variety of xnm5s2U derivatives. The synthesis pathway for these modifications is complex and fully elucidated only in a handful of organisms, including the Gram-negative Escherichia coli K12 model. Despite the ubiquitous presence of mnm5s2U modification, genomic analysis shows the absence of mnmC orthologous genes, suggesting the occurrence of alternate biosynthetic schemes for the conversion of cmnm5s2U to mnm5s2U. Using a combination of comparative genomics and genetic studies, a member of the YtqA subgroup of the radical Sam superfamily was found to be involved in the synthesis of mnm5s2U in both Bacillus subtilis and Streptococcus mutans. This protein, renamed MnmL, is encoded in an operon with the recently discovered MnmM methylase involved in the methylation of the pathway intermediate nm5s2U into mnm5s2U in B. subtilis. Analysis of tRNA modifications of both S. mutans and Streptococcus pneumoniae shows that growth conditions and genetic backgrounds influence the ratios of pathway intermediates owing to regulatory loops that are not yet understood. The MnmLM pathway is widespread along the bacterial tree, with some phyla, such as Bacilli, relying exclusively on these two enzymes. Although mechanistic details of these newly discovered components are not fully resolved, the occurrence of fusion proteins, alternate arrangements of biosynthetic components, and loss of biosynthetic branches provide examples of biosynthetic diversity to retain a conserved tRNA modification in Nature.IMPORTANCEThe xnm5s2U modifications found in several tRNAs at the wobble base position are widespread in bacteria where they have an important role in decoding efficiency and accuracy. This work identifies a novel enzyme (MnmL) that is a member of a subgroup of the very versatile radical SAM superfamily and is involved in the synthesis of mnm5s2U in several Gram-positive bacteria, including human pathogens. This is another novel example of a non-orthologous displacement in the field of tRNA modification synthesis, showing how different solutions evolve to retain U34 tRNA modifications.

Biosynthesis and function of 7-deazaguanine derivatives in bacteria and phages

SUMMARYDeazaguanine modifications play multifaceted roles in the molecular biology of DNA and tRNA, shaping diverse yet essential biological processes, including the nuanced fine-tuning of translation efficiency and the intricate modulation of codon-anticodon interactions. Beyond their roles in translation, deazaguanine modifications contribute to cellular stress resistance, self-nonself discrimination mechanisms, and host evasion defenses, directly modulating the adaptability of living organisms. Deazaguanine moieties extend beyond nucleic acid modifications, manifesting in the structural diversity of biologically active natural products. Their roles in fundamental cellular processes and their presence in biologically active natural products underscore their versatility and pivotal contributions to the intricate web of molecular interactions within living organisms. Here, we discuss the current understanding of the biosynthesis and multifaceted functions of deazaguanines, shedding light on their diverse and dynamic roles in the molecular landscape of life.

Analyzing RNA posttranscriptional modifications to decipher the epitranscriptomic code

The discovery of RNA silencing has revealed that non-protein-coding sequences (ncRNAs) can cover essential roles in regulatory networks and their malfunction may result in severe consequences on human health. These findings have prompted a general reassessment of the significance of RNA as a key player in cellular processes. This reassessment, however, will not be complete without a greater understanding of the distribution and function of the over 170 variants of the canonical ribonucleotides, which contribute to the breathtaking structural diversity of natural RNA. This review surveys the analytical approaches employed for the identification, characterization, and detection of RNA posttranscriptional modifications (rPTMs). The merits of analyzing individual units after exhaustive hydrolysis of the initial biopolymer are outlined together with those of identifying their position in the sequence of parent strands. Approaches based on next generation sequencing and mass spectrometry technologies are covered in depth to provide a comprehensive view of their respective merits. Deciphering the epitranscriptomic code will require not only mapping the location of rPTMs in the various classes of RNAs, but also assessing the variations of expression levels under different experimental conditions. The fact that no individual platform is currently capable of meeting all such demands implies that it will be essential to capitalize on complementary approaches to obtain the desired information. For this reason, the review strived to cover the broadest possible range of techniques to provide readers with the fundamental elements necessary to make informed choices and design the most effective possible strategy to accomplish the task at hand.

Alternate routes to mnm 5 s 2 U synthesis in Gram-positive bacteria

The wobble bases of tRNAs that decode split codons are often heavily modified. In Bacteria tRNA Glu, Gln, Asp contain a variety of xnm 5 s 2 U derivatives. The synthesis pathway for these modifications is complex and fully elucidated only in a handful of organisms, including the Gram-negative Escherichia coli K12 model. Despite the ubiquitous presence of mnm 5 s 2 U modification, genomic analysis shows the absence of mnmC orthologous genes, suggesting the occurrence of alternate biosynthetic schemes for the installation of this modification. Using a combination of comparative genomics and genetic studies, a member of the YtqA subgroup of the Radical Sam superfamily was found to be involved in the synthesis of mnm 5 s 2 U in both Bacillus subtilis and Streptococcus mutans . This protein, renamed MnmL, is encoded in an operon with the recently discovered MnmM methylase involved in the methylation of the pathway intermediate nm 5 s 2 U into mnm 5 s 2 U in B. subtilis . Analysis of tRNA modifications of both S. mutans and Streptococcus pneumoniae shows that growth conditions and genetic backgrounds influence the ratios of pathways intermediates in regulatory loops that are not yet understood. The MnmLM pathway is widespread along the bacterial tree, with some phyla, such as Bacilli, relying exclusively on these two enzymes. The occurrence of fusion proteins, alternate arrangements of biosynthetic components, and loss of biosynthetic branches provide examples of biosynthetic diversity to retain a conserved tRNA modification in nature.

Importance

The xnm 5 s 2 U modifications found in several tRNAs at the wobble base position are widespread in Bacteria where they have an important role in decoding efficiency and accuracy. This work identifies a novel enzyme (MnmL) that is a member of a subgroup of the very versatile Radical SAM superfamily and is involved in the synthesis of mnm 5 s 2 U in several Gram-positive bacteria, including human pathogens. This is another novel example of a non-orthologous displacement in the field of tRNA modification synthesis, showing how different solutions evolve to retain U34 tRNA modifications.

Direct sequencing of total Saccharomyces cerevisiae tRNAs by LC-MS/MS

Among RNAs, transfer RNAs (tRNAs) contain the widest variety of abundant posttranscriptional chemical modifications. These modifications are crucial for tRNAs to participate in protein synthesis, promoting proper tRNA structure and aminoacylation, facilitating anticodon:codon recognition, and ensuring the reading frame maintenance of the ribosome. While tRNA modifications were long thought to be stoichiometric, it is becoming increasingly apparent that these modifications can change dynamically in response to the cellular environment. The ability to broadly characterize the fluctuating tRNA modification landscape will be essential for establishing the molecular level contributions of individual sites of tRNA modification. The locations of modifications within individual tRNA sequences can be mapped using liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS). In this approach, a single tRNA species is purified, treated with ribonucleases, and the resulting single-stranded RNA products are subject to LC-MS/MS analysis. The application of LC-MS/MS to study tRNAs is limited by the necessity of analyzing one tRNA at a time, because the digestion of total tRNA mixtures by commercially available ribonucleases produces many short digestion products unable to be uniquely mapped back to a single site within a tRNA. We overcame these limitations by taking advantage of the highly structured nature of tRNAs to prevent the full digestion by single-stranded RNA-specific ribonucleases. Folding total tRNA prior to digestion allowed us to sequence Saccharomyces cerevisiae tRNAs with up to 97% sequence coverage for individual tRNA species by LC-MS/MS. This method presents a robust avenue for directly detecting the distribution of modifications in total tRNAs.

High-throughput profiling of RNA modifications by ultra-performance liquid chromatography coupled to complementary mass spectrometry: Methods, quality control, and applications

Although next-generation sequencing technology has been used to delineate RNA modifications in recent years, the paucity of appropriate converting reactions or specific antibodies impedes the accurate characterization and quantification of numerous RNA modifications, especially when these modifications demonstrate wide variations across developmental stages and cell types. In this study, we developed a high-throughput analytical platform coupling ultra-performance liquid chromatograph (UPLC) with complementary mass spectrometry (MS) to identify and quantify RNA modifications in both synthetic and biological samples. Sixty-four types of RNA modifications, including positional isomers and hypermodified ribonucleosides, were successfully monitored within a 16-min single run of UPLC-MS. Two independent methods to cross-validate the purity of RNA extracted from Caenorhabditis elegans (C. elegans) were developed using the coexisting C. elegans and Escherichia coli (E. coli) as a surveillance system. To test the validity of the method, we investigated the RNA modification landscape of three model organisms, C. elegans, E. coli, and Arabidopsis thaliana (A. thaliana). Both the identity and molarity of modified ribonucleosides markedly varied among the species. Moreover, our platform is not only useful for exploring the dynamics of RNA modifications in response to environmental cues (e.g., cold shock) but can also help with the identification of RNA-modifying enzymes in genetic studies. Cumulatively, our method presents a novel platform for the comprehensive analysis of RNA modifications, which will be of benefit to both analytical chemists involved in biomarker discovery and biologists conducting functional studies of RNA modifications.

Detection technologies for RNA modifications

To date, more than 170 chemical modifications have been characterized in RNA, providing a new layer of gene expression regulation termed the 'epitranscriptome'. RNA modification detection methods and tools advance the functional studies of the epitranscriptome. According to the detection throughput and principles, existing RNA modification detection technologies can be categorized into four classes, including quantification methods, locus-specific detection methods, next-generation sequencing-based detection technologies and nanopore direct RNA sequencing-based technologies. In this review, we summarize the current knowledge about these RNA modification detection technologies and discuss the challenges for the existing detection tools, providing information for a comprehensive understanding of the epitranscriptome.

The Ramp Atlas: facilitating tissue and cell-specific ramp sequence analyses through an intuitive web interface

Ramp sequences occur when the average translational efficiency of codons near the 5′ end of highly expressed genes is significantly lower than the rest of the gene sequence, which counterintuitively increases translational efficiency by decreasing downstream ribosomal collisions. Here, we show that the relative codon adaptiveness within different tissues changes the existence of a ramp sequence without altering the underlying genetic code. We present the first comprehensive analysis of tissue and cell type-specific ramp sequences and report 3108 genes with ramp sequences that change between tissues and cell types, which corresponds with increased gene expression within those tissues and cells. The Ramp Atlas (https://ramps.byu.edu/) allows researchers to query precomputed ramp sequences in 18 388 genes across 62 tissues and 66 cell types and calculate tissue-specific ramp sequences from user-uploaded FASTA files through an intuitive web interface. We used The Ramp Atlas to identify seven SARS-CoV-2 genes and seven human SARS-CoV-2 entry factor genes with tissue-specific ramp sequences that may help explain viral proliferation within those tissues. We anticipate that The Ramp Atlas will facilitate personalized and creative tissue-specific ramp sequence analyses for both human and viral genes that will increase our ability to utilize this often-overlooked regulatory region.

Omics Approaches to Assess Flavor Development in Cheese

Cheese is characterized by a rich and complex microbiota that plays a vital role during both production and ripening, contributing significantly to the safety, quality, and sensory characteristics of the final product. In this context, it is vital to explore the microbiota composition and understand its dynamics and evolution during cheese manufacturing and ripening. Application of high-throughput DNA sequencing technologies have facilitated the more accurate identification of the cheese microbiome, detailed study of its potential functionality, and its contribution to the development of specific organoleptic properties. These technologies include amplicon sequencing, whole-metagenome shotgun sequencing, metatranscriptomics, and, most recently, metabolomics. In recent years, however, the application of multiple meta-omics approaches along with data integration analysis, which was enabled by advanced computational and bioinformatics tools, paved the way to better comprehension of the cheese ripening process, revealing significant associations between the cheese microbiota and metabolites, as well as their impact on cheese flavor and quality.

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

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

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.

Role of a cryptic tRNA gene operon in survival under translational stress

As compared to eukaryotes, bacteria have a reduced tRNA gene set encoding between 30 and 220 tRNAs. Although in most bacterial phyla tRNA genes are dispersed in the genome, many species from distinct phyla also show genes forming arrays. Here, we show that two types of arrays with distinct evolutionary origins exist. This work focuses on long tRNA gene arrays (L-arrays) that encompass up to 43 genes, which disseminate by horizontal gene transfer and contribute supernumerary tRNA genes to the host. Although in the few cases previously studied these arrays were reported to be poorly transcribed, here we show that the L-array of the model cyanobacterium Anabaena sp. PCC 7120, encoding 23 functional tRNAs, is largely induced upon impairment of the translation machinery. The cellular response to this challenge involves a global reprogramming of the transcriptome in two phases. tRNAs encoded in the array are induced in the second phase of the response, directly contributing to cell survival. Results presented here show that in some bacteria the tRNA gene set may be partitioned between a housekeeping subset, which constantly sustains translation, and an inducible subset that is generally silent but can provide functionality under particular conditions.

Transcriptomic time-series analysis of cold- and heat-shock response in psychrotrophic lactic acid bacteria

BACKGROUND

Psychrotrophic lactic acid bacteria (LAB) species are the dominant species in the microbiota of cold-stored modified-atmosphere-packaged food products and are the main cause of food spoilage. Despite the importance of psychrotrophic LAB, their response to cold or heat has not been studied. Here, we studied the transcriptome-level cold- and heat-shock response of spoilage lactic acid bacteria with time-series RNA-seq for Le. gelidum, Lc. piscium, and P. oligofermentans at 0 °C, 4 °C, 14 °C, 25 °C, and 28 °C.

RESULTS

We observed that the cold-shock protein A (cspA) gene was the main cold-shock protein gene in all three species. Our results indicated that DEAD-box RNA helicase genes (cshA, cshB) also play a critical role in cold-shock response in psychrotrophic LAB. In addition, several RNase genes were involved in cold-shock response in Lc. piscium and P. oligofermentans. Moreover, gene network inference analysis provided candidate genes involved in cold-shock response. Ribosomal proteins, tRNA modification, rRNA modification, and ABC and efflux MFS transporter genes clustered with cold-shock response genes in all three species, indicating that these genes could be part of the cold-shock response machinery. Heat-shock treatment caused upregulation of Clp protease and chaperone genes in all three species. We identified transcription binding site motifs for heat-shock response genes in Le. gelidum and Lc. piscium. Finally, we showed that food spoilage-related genes were upregulated at cold temperatures.

CONCLUSIONS

The results of this study provide new insights on the cold- and heat-shock response of psychrotrophic LAB. In addition, candidate genes involved in cold- and heat-shock response predicted using gene network inference analysis could be used as targets for future studies.

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

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

Functions of Bacterial tRNA Modifications: From Ubiquity to Diversity

Modified nucleotides in tRNA are critical components of the translation apparatus, but their importance in the process of translational regulation had until recently been greatly overlooked. Two breakthroughs have recently allowed a fuller understanding of the importance of tRNA modifications in bacterial physiology. One is the identification of the full set of tRNA modification genes in model organisms such as Escherichia coli K12. The second is the improvement of available analytical tools to monitor tRNA modification patterns. The role of tRNA modifications varies greatly with the specific modification within a given tRNA and with the organism studied. The absence of these modifications or reductions can lead to cell death or pleiotropic phenotypes or may have no apparent visible effect. By linking translation through their decoding functions to metabolism through their biosynthetic pathways, tRNA modifications are emerging as important components of the bacterial regulatory toolbox.

Survey and Validation of tRNA Modifications and Their Corresponding Genes in Bacillus subtilis sp Subtilis Strain 168

Extensive knowledge of both the nature and position of tRNA modifications in all cellular tRNAs has been limited to two bacteria, Escherichia coli and Mycoplasma capricolum. Bacillus subtilis sp subtilis strain 168 is the model Gram-positive bacteria and the list of the genes involved in tRNA modifications in this organism is far from complete. Mass spectrometry analysis of bulk tRNA extracted from B. subtilis, combined with next generation sequencing technologies and comparative genomic analyses, led to the identification of 41 tRNA modification genes with associated confidence scores. Many differences were found in this model Gram-positive bacteria when compared to E. coli. In general, B. subtilis tRNAs are less modified than those in E. coli, even if some modifications, such as m1A22 or ms2t6A, are only found in the model Gram-positive bacteria. Many examples of non-orthologous displacements and of variations in the most complex pathways are described. Paralog issues make uncertain direct annotation transfer from E. coli to B. subtilis based on homology only without further experimental validation. This difficulty was shown with the identification of the B. subtilis enzyme that introduces ψ at positions 31/32 of the tRNAs. This work presents the most up to date list of tRNA modification genes in B. subtilis, identifies the gaps in knowledge, and lays the foundation for further work to decipher the physiological role of tRNA modifications in this important model organism and other bacteria.

Identification of the 3-amino-3-carboxypropyl (acp) transferase enzyme responsible for acp3U formation at position 47 in Escherichia coli tRNAs

tRNAs from all domains of life contain modified nucleotides. However, even for the experimentally most thoroughly characterized model organism Escherichia coli not all tRNA modification enzymes are known. In particular, no enzyme has been found yet for introducing the acp³U modification at position 47 in the variable loop of eight E. coli tRNAs. Here we identify the so far functionally uncharacterized YfiP protein as the SAM-dependent 3-amino-3-carboxypropyl transferase catalyzing this modification and thereby extend the list of known tRNA modification enzymes in E. coli. Similar to the Tsr3 enzymes that introduce acp modifications at U or m¹Ψ nucleotides in rRNAs this protein contains a DTW domain suggesting that acp transfer reactions to RNA nucleotides are a general function of DTW domain containing proteins. The introduction of the acp³U-47 modification in E. coli tRNAs is promoted by the presence of the m⁷G-46 modification as well as by growth in rich medium. However, a deletion of the enzymes responsible for the modifications at position 46 and 47 in the variable loop of E. coli tRNAs did not lead to a clearly discernible phenotype suggesting that these two modifications play only a minor role in ensuring the proper function of tRNAs in E. coli.

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.

MODOMICS: a database of RNA modification pathways. 2017 update

MODOMICS is a database of RNA modifications that provides comprehensive information concerning the chemical structures of modified ribonucleosides, their biosynthetic pathways, the location of modified residues in RNA sequences, and RNA-modifying enzymes. In the current database version, we included the following new features and data: extended mass spectrometry and liquid chromatography data for modified nucleosides; links between human tRNA sequences and MINTbase - a framework for the interactive exploration of mitochondrial and nuclear tRNA fragments; new, machine-friendly system of unified abbreviations for modified nucleoside names; sets of modified tRNA sequences for two bacterial species, updated collection of mammalian tRNA modifications, 19 newly identified modified ribonucleosides and 66 functionally characterized proteins involved in RNA modification. Data from MODOMICS have been linked to the RNAcentral database of RNA sequences. MODOMICS is available at http://modomics.genesilico.pl.

To be or not to be modified: Miscellaneous aspects influencing nucleotide modifications in tRNAs

Transfer RNAs (tRNAs) are essential components of the cellular protein synthesis machineries, but are also implicated in many roles outside translation. To become functional, tRNAs, initially transcribed as longer precursor tRNAs, undergo a tightly controlled biogenesis process comprising the maturation of their extremities, removal of intronic sequences if present, addition of the 3'-CCA amino-acid accepting sequence, and aminoacylation. In addition, the most impressive feature of tRNA biogenesis consists in the incorporation of a large number of posttranscriptional chemical modifications along its sequence. The chemical nature of these modifications is highly diverse, with more than hundred different modifications identified in tRNAs to date. All functions of tRNAs in cells are controlled and modulated by modifications, making the understanding of the mechanisms that determine and influence nucleotide modifications in tRNAs an essential point in tRNA biology. This review describes the different aspects that determine whether a certain position in a tRNA molecule is modified or not. We describe how sequence and structural determinants, as well as the presence of prior modifications control modification processes. We also describe how environmental factors and cellular stresses influence the level and/or the nature of certain modifications introduced in tRNAs, and report situations where these dynamic modulations of tRNA modification levels are regulated by active demodification processes. © 2019 IUBMB Life, 71(8):1126-1140, 2019.

Mapping post-transcriptional modifications in Staphylococcus aureus tRNAs by nanoLC/MSMS

RNA modifications are involved in numerous biological processes. These modifications are constitutive or modulated in response to adaptive processes and can impact RNA base-pairing formation, protein recognition, RNA structure and stability. tRNAs are the most abundantly modified RNA molecules. Analysis of the roles of their modifications in response to stress, environmental changes, and infections caused by pathogens, has fueled new research areas. Nevertheless, the detection of modified nucleotides in RNAs is still a challenging task. We present here a reliable method to identify and localize tRNA modifications, which was applied to the human pathogenic bacteria, Staphyloccocus aureus. The method is based on a separation of tRNA species on a two-dimensional polyacrylamide gel electrophoresis followed by nano liquid chromatography-mass spectrometry. We provided a list of modifications mapped on 25 out of the 40 tRNA species (one isoacceptor for each amino acid). This method can be easily used to monitor the dynamics of tRNA modifications in S. aureus in response to stress adaptation and during infection of the host, a relatively unexplored field.

tRNA Modification Profiles and Codon-Decoding Strategies in Methanocaldococcus jannaschii

tRNAs play a critical role in mRNA decoding, and posttranscriptional modifications within tRNAs drive decoding efficiency and accuracy. The types and positions of tRNA modifications in model bacteria have been extensively studied, and tRNA modifications in a few eukaryotic organisms have also been characterized and localized to particular tRNA sequences. However, far less is known regarding tRNA modifications in archaea. While the identities of modifications have been determined for multiple archaeal organisms, Haloferax volcanii is the only organism for which modifications have been extensively localized to specific tRNA sequences. To improve our understanding of archaeal tRNA modification patterns and codon-decoding strategies, we have used liquid chromatography and tandem mass spectrometry to characterize and then map posttranscriptional modifications on 34 of the 35 unique tRNA sequences of Methanocaldococcus jannaschii A new posttranscriptionally modified nucleoside, 5-cyanomethyl-2-thiouridine (cnm⁵s²U), was discovered and localized to position 34. Moreover, data consistent with wyosine pathway modifications were obtained beyond the canonical tRNA^Phe as is typical for eukaryotes. The high-quality mapping of tRNA anticodon loops enriches our understanding of archaeal tRNA modification profiles and decoding strategies.IMPORTANCE While many posttranscriptional modifications in M. jannaschii tRNAs are also found in bacteria and eukaryotes, several that are unique to archaea were identified. By RNA modification mapping, the modification profiles of M. jannaschii tRNA anticodon loops were characterized, allowing a comparative analysis with H. volcanii modification profiles as well as a general comparison with bacterial and eukaryotic decoding strategies. This general comparison reveals that M. jannaschii, like H. volcanii, follows codon-decoding strategies similar to those used by bacteria, although position 37 appears to be modified to a greater extent than seen in H. volcanii.

RNA ribose methylation (2'-O-methylation): Occurrence, biosynthesis and biological functions

Methylation of riboses at 2'-OH group is one of the most common RNA modifications found in number of cellular RNAs from almost any species which belong to all three life domains. This modification was extensively studied for decades in rRNAs and tRNAs, but recent data revealed the presence of 2'-O-methyl groups also in low abundant RNAs, like mRNAs. Ribose methylation is formed in RNA by two alternative enzymatic mechanisms: either by stand-alone protein enzymes or by complex assembly of proteins associated with snoRNA guides (sno(s)RNPs). In that case one catalytic subunit acts at various RNA sites, the specificity is provided by base pairing of the sno(s)RNA guide with the target RNA. In this review we compile available information on 2'-OH ribose methylation in different RNAs, enzymatic machineries involved in their biosynthesis and dynamics, as well as on the physiological functions of these modified residues.

Microbiome characterization by high-throughput transfer RNA sequencing and modification analysis

Advances in high-throughput sequencing have facilitated remarkable insights into the diversity and functioning of naturally occurring microbes; however, current sequencing strategies are insufficient to reveal physiological states of microbial communities associated with protein translation dynamics. Transfer RNAs (tRNAs) are core components of protein synthesis machinery, present in all living cells, and are phylogenetically tractable, which make them ideal targets to gain physiological insights into environmental microbes. Here we report a direct sequencing approach, tRNA-seq, and a software suite, tRNA-seq-tools, to recover sequences, abundance profiles, and post-transcriptional modifications of microbial tRNA transcripts. Our analysis of cecal samples using tRNA-seq distinguishes high-fat- and low-fat-fed mice in a comparable fashion to 16S ribosomal RNA gene amplicons, and reveals taxon- and diet-dependent variations in tRNA modifications. Our results provide taxon-specific in situ insights into the dynamics of tRNA gene expression and post-transcriptional modifications within complex environmental microbiomes.

Improving RNA modification mapping sequence coverage by LC-MS through a nonspecific RNase U2-E49A mutant

We report the identification and use of a mutant of the purine selective ribonuclease RNase U2 that randomly cleaves RNA in a manner that is directly compatible with RNA modification mapping by mass spectrometry. A number of RNase U2 mutants were generated using site-saturation mutagenesis. The enzyme activity and specificity were tested using oligonucleotide substrates, which revealed an RNase U2 E49A mutant with limited specificity and a tendency to undercut RNA. Using this mutant, RNA digestion conditions were optimized to yield long, overlapping digestion products, which improve sequence coverage in RNA modification mapping experiments. The analytical utility of this mutant was demonstrated by liquid chromatography tandem mass spectrometry (LC-MS/MS) mapping of several modified RNAs where 100% sequence coverage could be obtained using only a single enzymatic digestion. This new mutant facilitates more accurate and efficient RNA modification mapping than traditional highly base-specific RNases that are currently used.

Detection of ribonucleoside modifications by liquid chromatography coupled with mass spectrometry

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

Cells alter their tRNA abundance to selectively regulate protein synthesis during stress conditions

Decoding the information in mRNA during protein synthesis relies on tRNA adaptors, the abundance of which can affect the decoding rate and translation efficiency. To determine whether cells alter tRNA abundance to selectively regulate protein expression, we quantified changes in the abundance of individual tRNAs at different time points in response to diverse stress conditions in Saccharomyces cerevisiae. We found that the tRNA pool was dynamic and rearranged in a manner that facilitated selective translation of stress-related transcripts. Through genomic analysis of multiple data sets, stochastic simulations, and experiments with designed sequences of proteins with identical amino acids but altered codon usage, we showed that changes in tRNA abundance affected protein expression independently of factors such as mRNA abundance. We suggest that cells alter their tRNA abundance to selectively affect the translation rates of specific transcripts to increase the amounts of required proteins under diverse stress conditions.

Acetate-dependent tRNA acetylation required for decoding fidelity in protein synthesis

Modification of tRNA anticodons plays a critical role in ensuring accurate translation. N⁴-acetylcytidine (ac⁴C) is present at the anticodon first position (position 34) of bacterial elongator tRNA^Met. Herein, we identified Bacillus subtilis ylbM (renamed tmcAL) as a novel gene responsible for ac⁴C34 formation. Unlike general acetyltransferases that use acetyl-CoA, TmcAL activates an acetate ion to form acetyladenylate and then catalyzes ac⁴C34 formation through a mechanism similar to tRNA aminoacylation. The crystal structure of TmcAL with an ATP analog reveals the molecular basis of ac⁴C34 formation. The ΔtmcAL strain displayed a cold-sensitive phenotype and a strong genetic interaction with tilS that encodes the enzyme responsible for synthesizing lysidine (L) at position 34 of tRNA^Ile to facilitate AUA decoding. Mistranslation of the AUA codon as Met in the ΔtmcAL strain upon tilS repression suggests that ac⁴C34 modification of tRNA^Met and L34 modification of tRNA^Ile act cooperatively to prevent misdecoding of the AUA codon.

Early Transcriptome Response of Lactococcus lactis to Environmental Stresses Reveals Differentially Expressed Small Regulatory RNAs and tRNAs

Bacteria can deploy various mechanisms to combat environmental stresses. Many genes have previously been identified in Lactococcus lactis that are involved in sensing the stressors and those that are involved in regulating and mounting a defense against the stressful conditions. However, the expression of small regulatory RNAs (sRNAs) during industrially relevant stress conditions has not been assessed yet in L. lactis, while sRNAs have been shown to be involved in many stress responses in other bacteria. We have previously reported the presence of hundreds of putative regulatory RNAs in L. lactis, and have used high-throughput RNA sequencing (RNA-seq) in this study to assess their expression under six different stress conditions. The uniformly designed experimental set-up enabled a highly reliable comparison between the different stress responses and revealed that many sRNAs are differentially expressed under the conditions applied. The primary stress responses of L. lactis NCDO712 was benchmarked to earlier work and, for the first time, the differential expression was assessed of transfer RNAs (tRNAs) and the genes from the six recently sequenced plasmids of NCDO712. Although, we only applied stresses for 5 min, the majority of the well-known specific stress-induced genes are already differentially expressed. We find that most tRNAs decrease after all stresses applied, except for a small number, which are increased upon cold stress. Starvation was shown to induce the highest differential response, both in terms of number and expression level of genes. Our data pinpoints many novel stress-related uncharacterized genes and sRNAs, which calls for further assessment of their molecular and cellular function. These insights furthermore could impact the way parameters are designed for bacterial culture production and milk fermentation, as we find that very short stress conditions already greatly alter gene expression.

Stable Isotope Labeling for Improved Comparative Analysis of RNA Digests by Mass Spectrometry

Even with the advent of high throughput methods to detect modified ribonucleic acids (RNAs), mass spectrometry remains a reliable method to detect, characterize, and place post-transcriptional modifications within an RNA sequence. Here we have developed a stable isotope labeling comparative analysis of RNA digests (SIL-CARD) approach, which improves upon the original ¹⁸O/¹⁶O labeling CARD method. Like the original, SIL-CARD allows sequence or modification information from a previously uncharacterized in vivo RNA sample to be obtained by direct comparison with a reference RNA, the sequence of which is known. This reference is in vitro transcribed using a ¹³C/¹⁵N isotopically enriched nucleoside triphosphate (NTP). The two RNAs are digested with an endonuclease, the specificity of which matches the labeled NTP used for transcription. As proof of concept, several transfer RNAs (tRNAs) were characterized by SIL-CARD, where labeled guanosine triphosphate was used for the reference in vitro transcription. RNase T1 digestion products from the in vitro transcript will be 15 Da higher in mass than the same digestion products from the in vivo tRNA that are unmodified, leading to a doublet in the mass spectrum. Singlets, rather than doublets, arise if a sequence variation or a post-transcriptional modification is present that results in a relative mass shift different from 15 Da. Moreover, the use of the in vitro synthesized tRNA transcript allows for quantitative measurement of RNA abundance. Overall, SIL-CARD simplifies data analysis and enhances quantitative RNA modification mapping by mass spectrometry. Graphical Abstract ᅟ.

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.

Transfer RNA is highly unstable during early amino acid starvation in Escherichia coli

Due to its long half-life compared to messenger RNA, bacterial transfer RNA is known as stable RNA. Here, we show that tRNAs become highly unstable as part of Escherichia coli's response to amino acid starvation. Degradation of the majority of cellular tRNA occurs within twenty minutes of the onset of starvation for each of several amino acids. Both the non-cognate and cognate tRNA for the amino acid that the cell is starving for are degraded, and both charged and uncharged tRNA species are affected. The alarmone ppGpp orchestrates the stringent response to amino acid starvation. However, tRNA degradation occurs in a ppGpp-independent manner, as it occurs with similar kinetics in a relaxed mutant. Further, we also observe rapid tRNA degradation in response to rifampicin treatment, which does not induce the stringent response. We propose a unifying model for these observations, in which the surplus tRNA is degraded whenever the demand for protein synthesis is reduced. Thus, the tRNA pool is a highly regulated, dynamic entity. We propose that degradation of surplus tRNA could function to reduce mistranslation in the stressed cell, because it would reduce competition between cognate and near-cognate charged tRNAs at the ribosomal A-site.

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.

Comparative Analysis of Ribonucleic Acid Digests (CARD) by Mass Spectrometry

We describe the comparative analysis of ribonucleic acid digests (CARD) approach for RNA modification analysis. This approach employs isotope labeling during RNase digestion, which allows the direct comparison of a tRNA of unknown modification status against a reference tRNA, whose sequence or modification status is known. The reference sample is labeled with ¹⁸O during RNase digestion while the candidate (unknown) sample is labeled with ¹⁶O. These RNase digestion products are combined and analyzed by mass spectrometry. Identical RNase digestion products will appear in the mass spectrum as characteristic doublets, separated by 2 Da due to the ¹⁶O/¹⁸O mass difference. Singlets arise in the mass spectrum when the sequence or modification status of a particular RNase digestion product from the reference is not matched in the candidate (unknown) sample. This CARD approach for RNA modification analysis simplifies the determination of differences between reference and candidate samples, providing a route for higher throughput screening of samples for modification profiles, including determination of tRNA methylation patterns.

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.

Identification of the mRNA targets of tRNA-specific regulation using genome-wide simulation of translation

tRNA gene copy number is a primary determinant of tRNA abundance and therefore the rate at which each tRNA delivers amino acids to the ribosome during translation. Low-abundance tRNAs decode rare codons slowly, but it is unclear which genes might be subject to tRNA-mediated regulation of expression. Here, those mRNA targets were identified via global simulation of translation. In-silico mRNA translation rates were compared for each mRNA in both wild-type and a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\rm{tRNA}}_{{\rm{CUG}}}^{{\rm{Gln}}}$\end{document}sup70-65 mutant, which exhibits a pseudohyphal growth phenotype and a 75% slower CAG codon translation rate. Of 4900 CAG-containing mRNAs, 300 showed significantly reduced in silico translation rates in a simulated tRNA mutant. Quantitative immunoassay confirmed that the reduced translation rates of sensitive mRNAs were \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\rm{tRNA}}_{{\rm{CUG}}}^{{\rm{Gln}}}$\end{document} concentration-dependent. Translation simulations showed that reduced \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${\rm{tRNA}}_{{\rm{CUG}}}^{{\rm{Gln}}}$\end{document} concentrations triggered ribosome queues, which dissipated at reduced translation initiation rates. To validate this prediction experimentally, constitutive gcn2 kinase mutants were used to reduce in vivo translation initiation rates. This repaired the relative translational rate defect of target mRNAs in the sup70-65 background, and ameliorated sup70-65 pseudohyphal growth phenotypes. We thus validate global simulation of translation as a new tool to identify mRNA targets of tRNA-specific gene regulation.

tRNAmodpred: A computational method for predicting posttranscriptional modifications in tRNAs

tRNA molecules contain numerous chemically altered nucleosides, which are formed by enzymatic modification of the primary transcripts during the complex tRNA maturation process. Some of the modifications are introduced by single reactions, while other require complex series of reactions carried out by several different enzymes. The location and distribution of various types of modifications vary greatly between different tRNA molecules, organisms and organelles. We have developed a computational method tRNAmodpred, for predicting modifications in tRNA sequences. Briefly, our method takes as an input one or more unmodified tRNA sequences and a set of protein sequences corresponding to a proteome of a cell. Subsequently it identifies homologs of known tRNA modification enzymes in the proteome, predicts tRNA modification activities and maps them onto known pathways of RNA modification from the MODOMICS database. Thereby, theoretically possible modification pathways are identified, and products of these modification reactions are proposed for query tRNAs. This method allows for predicting modification patterns for newly sequenced genomes as well as for checking tentative modification status of tRNAs from one species treated with enzymes from another source, e.g. to predict the possible modifications of eukaryotic tRNAs expressed in bacteria. tRNAmodpred is freely available as a web server at http://genesilico.pl/trnamodpred/.

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.

Distribution and frequencies of post-transcriptional modifications in tRNAs

Functional tRNA molecules always contain a wide variety of post-transcriptionally modified nucleosides. These modifications stabilize tRNA structure, allow for proper interaction with other macromolecules and fine-tune the decoding of mRNAs during translation. Their presence in functionally important regions of tRNA is conserved in all domains of life. However, the identities of many of these modified residues depend much on the phylogeny of organisms the tRNAs are found in, attesting for domain-specific strategies of tRNA maturation. In this work we present a new tool, tRNAmodviz web server (http://genesilico.pl/trnamodviz) for easy comparative analysis and visualization of modification patterns in individual tRNAs, as well as in groups of selected tRNA sequences. We also present results of comparative analysis of tRNA sequences derived from 7 phylogenetically distinct groups of organisms: Gram-negative bacteria, Gram-positive bacteria, cytosol of eukaryotic single cell organisms, Fungi and Metazoa, cytosol of Viridiplantae, mitochondria, plastids and Euryarchaeota. These data update the study conducted 20 y ago with the tRNA sequences available at that time.

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.

Detection of RNA nucleoside modifications with the uridine-specific ribonuclease MC1 from Momordica charantia

A codon-optimized recombinant ribonuclease, MC1 is characterized for its uridine-specific cleavage ability to map nucleoside modifications in RNA. The published MC1 amino acid sequence, as noted in a previous study, was used as a template to construct a synthetic gene with a natural codon bias favoring expression in Escherichia coli. Following optimization of various expression conditions, the active recombinant ribonuclease was successfully purified as a C-terminal His-tag fusion protein from E. coli [Rosetta 2(DE3)] cells. The isolated protein was tested for its ribonuclease activity against oligoribonucleotides and commercially available E. coli tRNA(Tyr I). Analysis of MC1 digestion products by ion-pairing reverse phase liquid-chromatography coupled with mass spectrometry (IP-RP-LC-MS) revealed enzymatic cleavage of RNA at the 5'-termini of uridine and pseudouridine, but cleavage was absent if the uridine was chemically modified or preceded by a nucleoside with a bulky modification. Furthermore, the utility of this enzyme to generate complementary digestion products to other common endonucleases, such as RNase T1, which enables the unambiguous mapping of modified residues in RNA is demonstrated.

Enhanced detection of post-transcriptional modifications using a mass-exclusion list strategy for RNA modification mapping by LC-MS/MS

There has been a renewed appreciation for the dynamic nature of ribonucleic acid (RNA) modifications and for the impact of modified RNAs on organism health resulting in an increased emphasis on developing analytical methods capable of detecting modifications within specific RNA sequence contexts. Here we demonstrate that a DNA-based exclusion list enhances data dependent liquid chromatography tandem mass spectrometry (LC-MS/MS) detection of post-transcriptionally modified nucleosides within specific RNA sequences. This approach is possible because all post-transcriptional modifications of RNA, except pseudouridine, result in a mass increase in the canonical nucleoside undergoing chemical modification. Thus, DNA-based sequences reflect the state of the RNA prior to or in the absence of modification. The utility of this exclusion list strategy is demonstrated through the RNA modification mapping of total tRNAs from the bacteria Escherichia coli, Lactococcus lactis, and Streptomyces griseus. Creation of a DNA-based exclusion list is shown to consistently enhance the number of detected modified ribonuclease (RNase) digestion products by ∼20%. All modified RNase digestion products that were detected during standard data dependent acquisition (DDA) LC-MS/MS were also detected when the DNA-based exclusion list was used. Consequently, the increase in detected modified RNase digestion products is attributed to new experimental information only obtained when using the exclusion list. This exclusion list strategy should be broadly applicable to any class of RNA and improves the utility of mass spectrometry approaches for discovery-based analyses of RNA modifications, such as are required for studies of the epitranscriptome.

Transfer RNAs Mediate the Rapid Adaptation of Escherichia coli to Oxidative Stress

Translational systems can respond promptly to sudden environmental changes to provide rapid adaptations to environmental stress. Unlike the well-studied translational responses to oxidative stress in eukaryotic systems, little is known regarding how prokaryotes respond rapidly to oxidative stress in terms of translation. In this study, we measured protein synthesis from the entire Escherichia coli proteome and found that protein synthesis was severely slowed down under oxidative stress. With unchanged translation initiation, this slowdown was caused by decreased translation elongation speed. We further confirmed by tRNA sequencing and qRT-PCR that this deceleration was caused by a global, enzymatic downregulation of almost all tRNA species shortly after exposure to oxidative agents. Elevation in tRNA levels accelerated translation and protected E. coli against oxidative stress caused by hydrogen peroxide and the antibiotic ciprofloxacin. Our results showed that the global regulation of tRNAs mediates the rapid adjustment of the E. coli translation system for prompt adaptation to oxidative stress.

All organisms need to respond quickly to sudden environmental changes. Translational regulation can occur in response to environmental stresses within minutes, which is much faster than transcriptional regulation, and thus normally provides immediate adaptation. Eukaryotic cells can manipulate their tRNA molecules, mainly in a reversible manner, to suppress translation. Here, we showed for the first time that bacteria respond to oxidative stress by adjusting the translational system in a manner that differs from that of eukaryotes. The bacteria nonspecifically, irreversibly, and enzymatically degrade tRNAs to block protein synthesis. Interestingly, we showed that elevated tRNA concentrations lead to opposing effects by causing increased protein aggregation, which impairs fitness under normal conditions but facilitates adaptation under oxidative stress, including that caused by antibiotics. Our results provide a new understanding of the role of global adjustments to the entire translation system during stress adaptation in bacteria. This mechanism may also be involved in the development of antibiotic resistance in bacteria.

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.

Optimization of translation profiles enhances protein expression and solubility

mRNA is translated with a non-uniform speed that actively coordinates co-translational folding of protein domains. Using structure-based homology we identified the structural domains in epoxide hydrolases (EHs) and introduced slow-translating codons to delineate the translation of single domains. These changes in translation speed dramatically improved the solubility of two EHs of metagenomic origin in Escherichia coli. Conversely, the importance of transient attenuation for the folding, and consequently solubility, of EH was evidenced with a member of the EH family from Agrobacterium radiobacter, which partitions in the soluble fraction when expressed in E. coli. Synonymous substitutions of codons shaping the slow-transiting regions to fast-translating codons render this protein insoluble. Furthermore, we show that low protein yield can be enhanced by decreasing the free folding energy of the initial 5’-coding region, which can disrupt mRNA secondary structure and enhance ribosomal loading. This study provides direct experimental evidence that mRNA is not a mere messenger for translation of codons into amino acids but bears an additional layer of information for folding, solubility and expression level of the encoded protein. Furthermore, it provides a general frame on how to modulate and fine-tune gene expression of a target protein.