Wobbling Forth and Drifting Back: The Evolutionary History and Impact of Bacterial tRNA Modifications

tRNA modifications: greasing the wheels of translation and beyond

Transfer RNA (tRNA) is one of the most abundant RNA types in cells, acting as an adaptor to bridge the genetic information in mRNAs with the amino acid sequence in proteins. Both tRNAs and small fragments processed from them play many nonconventional roles in addition to translation. tRNA molecules undergo various types of chemical modifications to ensure the accuracy and efficiency of translation and regulate their diverse functions beyond translation. In this review, we discuss the biogenesis and molecular mechanisms of tRNA modifications, including major tRNA modifications, writer enzymes, and their dynamic regulation. We also summarize the state-of-the-art technologies for measuring tRNA modification, with a particular focus on 2'-O-methylation (Nm), and discuss their limitations and remaining challenges. Finally, we highlight recent discoveries linking dysregulation of tRNA modifications with genetic diseases.

The Evolution and Implications of the Inosine tRNA Modification

Ever since the legendary publication by Francis Crick in JMB introducing the wobble hypothesis in 1966, inosine has been a permanent part of molecular biology. This review aims to integrate the rich array of novel insights emerging from subsequent research on the adenine-to-inosine modification of tRNA, with an emphasis on the results obtained during the last 5 years. Both the grand panorama of 4 billion years of evolution of life and the medical implications of defects in inosine modification will be reviewed. The most salient insights are that: (1) inosine at position 34 (the first position in the anticodon) is not universally present in the tree of life; (2) in many bacteria just a single homodimeric enzyme (TadA) is responsible for both tRNA inosine modification and mRNA inosine modification; (3) rapid progress is currently being made both in the molecular understanding of the heterodimeric ADAT2/ADAT3 enzyme responsible for inosine modifications in eukaryotes and in experimental capabilities for monitoring both the cytoplasmic tRNA pool and their modifications; (4) for selected tRNAs, inosine modification at position 37 has been demonstrated but this modification remains under-studied; (5) modification of tRNAs known to contain inosine can be incomplete; (6) the GC content of the T-stem is of great importance for wobble behavior, including wobbling behavior of inosine; and (7) the tRNA inosine modification is of direct relevance to human disease. In summary, research on inosine continues to yield important novel insights.

Distal Domains of the Bacterial-Exclusive Wobble-Modifying Enzyme TilS Contribute to Catalysis

tRNAIle2 lysidine synthetase (TilS) is a bacterial-specific wobble-modifying enzyme that acts on the isoleucine-accepting tRNAIle2 CAU. TilS installs a lysine at the C34 position of the anticodon, generating the lysidine modification. The resulting LAU anticodon enables exclusive decoding of infrequently used AUA isoleucine codons, rejecting AUG methionine codons. Compared to other wobble-modifying enzymes that contact only the anticodon arm of their cognate tRNAs, TilS is distinct in containing additional domains outside of the N-terminal active site. For type I TilS enzymes such as the B. cenocepacia TilS (BcTilS) investigated here, appended domains contact the tRNAIle2 substrate along the body and through the acceptor stem, up to 60 Å away from the target C34. Among bacterial tRNAs, only unmodified tRNAIle2 and tRNAMet share an anticodon, suggesting that the appended domains of TilS provide substrate recognition strategies that other wobble-modifying enzymes do not need. Here, we investigate both protein and tRNA elements to understand the strategy by which TilS accepts its cognate tRNAIle2 substrate and rejects the near-cognate tRNAMet.

Genome-wide profiling of tRNA modifications by Induro-tRNAseq reveals coordinated changes

While all native tRNAs undergo extensive post-transcriptional modifications as a mechanism to regulate gene expression, mapping these modifications remains challenging. The critical barrier is the difficulty of readthrough of modifications by reverse transcriptases (RTs). Here we use Induro-a new group-II intron-encoded RT-to map and quantify genome-wide tRNA modifications in Induro-tRNAseq. We show that Induro progressively increases readthrough over time by selectively overcoming RT stops without altering the misincorporation frequency. In a parallel analysis of Induro vs. a related RT, we provide comparative datasets to facilitate the prediction of each modification. We assess tRNA modifications across five human cell lines and three mouse tissues and show that, while the landscape of modifications is highly variable throughout the tRNA sequence framework, it is stabilized for modifications that are required for reading of the genetic code. The coordinated changes have fundamental importance for development of tRNA modifications in protein homeostasis.

Impact of the chemical modification of tRNAs anticodon loop on the variability and evolution of codon usage in proteobacteria

Despite the highly conserved nature of the genetic code, the frequency of usage of each codon can vary significantly. The evolution of codon usage is shaped by two main evolutionary forces: mutational bias and selection pressures. These pressures can be driven by environmental factors, but also by the need for efficient translation, which depends heavily on the concentration of transfer RNAs (tRNAs) within the cell. The data presented here supports the proposal that tRNA modifications play a key role in shaping the overall preference of codon usage in proteobacteria. Interestingly, some codons, such as CGA and AGG (encoding arginine), exhibit a surprisingly low level of variation in their frequency of usage, even across genomes with differing GC content. These findings suggest that the evolution of GC content in proteobacterial genomes might be primarily driven by changes in the usage of a specific subset of codons, whose usage is itself influenced by tRNA modifications.

A tRNA modification pattern that facilitates interpretation of the genetic code

Interpretation of the genetic code from triplets of nucleotides to amino acids is fundamental to life. This interpretation is achieved by cellular tRNAs, each reading a triplet codon through its complementary anticodon (positions 34-36) while delivering the amino acid charged to its 3'-end. This amino acid is then incorporated into the growing polypeptide chain during protein synthesis on the ribosome. The quality and versatility of the interpretation is ensured not only by the codon-anticodon pairing, but also by the post-transcriptional modifications at positions 34 and 37 of each tRNA, corresponding to the wobble nucleotide at the first position of the anticodon and the nucleotide on the 3'-side of the anticodon, respectively. How each codon is read by the matching anticodon, and which modifications are required, cannot be readily predicted from the codon-anticodon pairing alone. Here we provide an easily accessible modification pattern that is integrated into the genetic code table. We focus on the Gram-negative bacterium Escherichia coli as a model, which is one of the few organisms whose entire set of tRNA modifications and modification genes is identified and mapped. This work provides an important reference tool that will facilitate research in protein synthesis, which is at the core of the cellular life.

Examining the Effects of Temperature on the Evolution of Bacterial tRNA Pools

The genetic code consists of 61 codons coding for 20 amino acids. These codons are recognized by transfer RNAs (tRNAs) that bind to specific codons during protein synthesis. All organisms utilize less than all 61 possible anticodons due to base pair wobble: the ability to have a mismatch with a codon at its third nucleotide. Previous studies observed a correlation between the tRNA pool of bacteria and the temperature of their respective environments. However, it is unclear if these patterns represent biological adaptations to maintain the efficiency and accuracy of protein synthesis in different environments. A mechanistic mathematical model of mRNA translation is used to quantify the expected elongation rates and error rate for each codon based on an organism's tRNA pool. A comparative analysis across a range of bacteria that accounts for covariance due to shared ancestry is performed to quantify the impact of environmental temperature on the evolution of the tRNA pool. We find that thermophiles generally have more anticodons represented in their tRNA pool than mesophiles or psychrophiles. Based on our model, this increased diversity is expected to lead to increased missense errors. The implications of this for protein evolution in thermophiles are discussed.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Determining the effects of temperature on the evolution of bacterial tRNA pools

The genetic code consists of 61 codon coding for 20 amino acids. These codons are recognized by transfer RNAs (tRNA) that bind to specific codons during protein synthesis. Most organisms utilize less than all 61 possible anticodons due to base pair wobble: the ability to have a mismatch with a codon at its third nucleotide. Previous studies observed a correlation between the tRNA pool of bacteria and the temperature of their respective environments. However, it is unclear if these patterns represent biological adaptations to maintain the efficiency and accuracy of protein synthesis in different environments. A mechanistic mathematical model of mRNA translation is used to quantify the expected elongation rates and error rate for each codon based on an organism's tRNA pool. A comparative analysis across a range of bacteria that accounts for covariance due to shared ancestry is performed to quantify the impact of environmental temperature on the evolution of the tRNA pool. We find that thermophiles generally have more anticodons represented in their tRNA pool than mesophiles or psychrophiles. Based on our model, this increased diversity is expected to lead to increased missense errors. The implications of this for protein evolution in thermophiles are discussed.

Shifts in mutation spectra enhance access to beneficial mutations

Biased mutation spectra are pervasive, with wide variation in the magnitude of mutational biases that influence genome evolution and adaptation. How do such diverse biases evolve? Our experiments show that changing the mutation spectrum allows populations to sample previously undersampled mutational space, including beneficial mutations. The resulting shift in the distribution of fitness effects is advantageous: Beneficial mutation supply and beneficial pleiotropy both increase, while deleterious load reduces. More broadly, simulations indicate that reducing or reversing the direction of a long-term bias is always selectively favored. Such changes in mutation bias can occur easily via altered function of DNA repair genes. A phylogenetic analysis shows that these genes are repeatedly gained and lost in bacterial lineages, leading to frequent bias shifts in opposite directions. Thus, shifts in mutation spectra may evolve under selection and can directly alter the outcome of adaptive evolution by facilitating access to beneficial mutations.

The layered costs and benefits of translational redundancy

The rate and accuracy of translation hinges upon multiple components - including transfer RNA (tRNA) pools, tRNA modifying enzymes, and rRNA molecules - many of which are redundant in terms of gene copy number or function. It has been hypothesized that the redundancy evolves under selection, driven by its impacts on growth rate. However, we lack empirical measurements of the fitness costs and benefits of redundancy, and we have poor a understanding of how this redundancy is organized across components. We manipulated redundancy in multiple translation components of Escherichia coli by deleting 28 tRNA genes, 3 tRNA modifying systems, and 4 rRNA operons in various combinations. We find that redundancy in tRNA pools is beneficial when nutrients are plentiful and costly under nutrient limitation. This nutrient-dependent cost of redundant tRNA genes stems from upper limits to translation capacity and growth rate, and therefore varies as a function of the maximum growth rate attainable in a given nutrient niche. The loss of redundancy in rRNA genes and tRNA modifying enzymes had similar nutrient-dependent fitness consequences. Importantly, these effects are also contingent upon interactions across translation components, indicating a layered hierarchy from copy number of tRNA and rRNA genes to their expression and downstream processing. Overall, our results indicate both positive and negative selection on redundancy in translation components, depending on a species' evolutionary history with feasts and famines.

TadA orthologs enable both cytosine and adenine editing of base editors

Cytidine and adenosine deaminases are required for cytosine and adenine editing of base editors respectively, and no single deaminase could enable concurrent and comparable cytosine and adenine editing. Additionally, distinct properties of cytidine and adenosine deaminases lead to various types of off-target effects, including Cas9-indendepent DNA off-target effects for cytosine base editors (CBEs) and RNA off-target effects particularly severe for adenine base editors (ABEs). Here we demonstrate that 25 TadA orthologs could be engineered to generate functional ABEs, CBEs or ACBEs via single or double mutations, which display minimized Cas9-independent DNA off-target effects and genotoxicity, with orthologs B5ZCW4, Q57LE3, E8WVH3, Q13XZ4 and B3PCY2 as promising candidates for further engineering. Furthermore, RNA off-target effects of TadA ortholog-derived base editors could be further reduced or even eliminated by additional single mutation. Taken together, our work expands the base editing toolkits, and also provides important clues for the potential evolutionary process of deaminases.

On distinguishing between canonical tRNA genes and tRNA gene fragments in prokaryotes

Automated genome annotation is essential for extracting biological information from sequence data. The identification and annotation of tRNA genes is frequently performed by the software package tRNAscan-SE, the output of which is listed for selected genomes in the Genomic tRNA database (GtRNAdb). Here, we highlight a pervasive error in prokaryotic tRNA gene sets on GtRNAdb: the mis-categorization of partial, non-canonical tRNA genes as standard, canonical tRNA genes. Firstly, we demonstrate the issue using the tRNA gene sets of 20 organisms from the archaeal taxon Thermococcaceae. According to GtRNAdb, these organisms collectively deviate from the expected set of tRNA genes in 15 instances, including the listing of eleven putative canonical tRNA genes. However, after detailed manual annotation, only one of these eleven remains; the others are either partial, non-canonical tRNA genes resulting from the integration of genetic elements or CRISPR-Cas activity (seven instances), or attributable to ambiguities in input sequences (three instances). Secondly, we show that similar examples of the mis-categorization of predicted tRNA sequences occur throughout the prokaryotic sections of GtRNAdb. While both canonical and non-canonical prokaryotic tRNA gene sequences identified by tRNAscan-SE are biologically interesting, the challenge of reliably distinguishing between them remains. We recommend employing a combination of (i) screening input sequences for the genetic elements typically associated with non-canonical tRNA genes, and ambiguities, (ii) activating the tRNAscan-SE automated pseudogene detection function, and (iii) scrutinizing predicted tRNA genes with low isotype scores. These measures greatly reduce manual annotation efforts, and lead to improved prokaryotic tRNA gene set predictions.

A large-scale analysis of codon usage bias in 4868 bacterial genomes shows association of codon adaptation index with GC content, protein functional domains and bacterial phenotypes

Multiple synonymous codons code for the same amino acid, resulting in the degeneracy of the genetic code and in the preferred used of some codons called codon bias usage (CBU). We performed a large-scale analysis of codon usage bias analysing the distribution of the codon adaptation index (CAI) and the codon relative adaptiveness index (RA) in 4868 bacterial genomes. We found that CAI values differ significantly between protein functional domains and part of the protein outside domains and show how CAI, GC content and preferred usage of polymerase III alpha subunits are related. Additionally, we give evidence of the association between CAI and bacterial phenotypes.

Evolutionary jumps in bacterial GC content

Genomic GC (Guanine-Cytosine) content is a fundamental molecular trait linked with many key genomic features such as codon and amino acid use. Across bacteria, GC content is surprisingly diverse and has been studied for many decades; yet its evolution remains incompletely understood. Since it is difficult to observe GC content evolve on laboratory time scales, phylogenetic comparative approaches are instrumental; but this dimension is rarely studied systematically in the case of bacterial GC content. We applied phylogenetic comparative models to analyze GC content evolution in multiple bacterial groups across 2 major bacterial phyla. We find that GC content diversifies via a combination of gradual evolution and evolutionary "jumps." Surprisingly, unlike prior reports that solely focused on reductions in GC, we found a comparable number of jumps with both increased and decreased GC content. Overall, many of the identified jumps occur in lineages beyond the well-studied peculiar examples of endosymbiotic and AT-rich marine bacteria and do not support the predicted role of oxygen dependence. Our analysis of rapid and large shifts in GC content thus identifies new clades and novel contexts to further understand the ecological and evolutionary drivers of this important genomic trait.

Human tRNAs with inosine 34 are essential to efficiently translate eukarya-specific low-complexity proteins

The modification of adenosine to inosine at the wobble position (I34) of tRNA anticodons is an abundant and essential feature of eukaryotic tRNAs. The expansion of inosine-containing tRNAs in eukaryotes followed the transformation of the homodimeric bacterial enzyme TadA, which generates I34 in tRNAArg and tRNALeu, into the heterodimeric eukaryotic enzyme ADAT, which modifies up to eight different tRNAs. The emergence of ADAT and its larger set of substrates, strongly influenced the tRNA composition and codon usage of eukaryotic genomes. However, the selective advantages that drove the expansion of I34-tRNAs remain unknown. Here we investigate the functional relevance of I34-tRNAs in human cells and show that a full complement of these tRNAs is necessary for the translation of low-complexity protein domains enriched in amino acids cognate for I34-tRNAs. The coding sequences for these domains require codons translated by I34-tRNAs, in detriment of synonymous codons that use other tRNAs. I34-tRNA-dependent low-complexity proteins are enriched in functional categories related to cell adhesion, and depletion in I34-tRNAs leads to cellular phenotypes consistent with these roles. We show that the distribution of these low-complexity proteins mirrors the distribution of I34-tRNAs in the phylogenetic tree.

Inosine in Biology and Disease

The nucleoside inosine plays an important role in purine biosynthesis, gene translation, and modulation of the fate of RNAs. The editing of adenosine to inosine is a widespread post-transcriptional modification in transfer RNAs (tRNAs) and messenger RNAs (mRNAs). At the wobble position of tRNA anticodons, inosine profoundly modifies codon recognition, while in mRNA, inosines can modify the sequence of the translated polypeptide or modulate the stability, localization, and splicing of transcripts. Inosine is also found in non-coding and exogenous RNAs, where it plays key structural and functional roles. In addition, molecular inosine is an important secondary metabolite in purine metabolism that also acts as a molecular messenger in cell signaling pathways. Here, we review the functional roles of inosine in biology and their connections to human health.

On the Track of the Missing tRNA Genes: A Source of Non-Canonical Functions?

Cellular tRNAs appear today as a diverse population of informative macromolecules with conserved general elements ensuring essential common functions and different and distinctive features securing specific interactions and activities. Their differential expression and the variety of post-transcriptional modifications they are subject to, lead to the existence of complex repertoires of tRNA populations adjusted to defined cellular states. Despite the tRNA-coding genes redundancy in prokaryote and eukaryote genomes, it is surprising to note the absence of genes coding specific translational-active isoacceptors throughout the phylogeny. Through the analysis of different releases of tRNA databases, this review aims to provide a general summary about those “missing tRNA genes.” This absence refers to both tRNAs that are not encoded in the genome, as well as others that show critical sequence variations that would prevent their activity as canonical translation adaptor molecules. Notably, while a group of genes are universally missing, others are absent in particular kingdoms. Functional information available allows to hypothesize that the exclusion of isodecoding molecules would be linked to: 1) reduce ambiguities of signals that define the specificity of the interactions in which the tRNAs are involved; 2) ensure the adaptation of the translational apparatus to the cellular state; 3) divert particular tRNA variants from ribosomal protein synthesis to other cellular functions. This leads to consider the “missing tRNA genes” as a source of putative non-canonical tRNA functions and to broaden the concept of adapter molecules in ribosomal-dependent protein synthesis.

The birth of a bacterial tRNA gene by large-scale, tandem duplication events

Organisms differ in the types and numbers of tRNA genes that they carry. While the evolutionary mechanisms behind tRNA gene set evolution have been investigated theoretically and computationally, direct observations of tRNA gene set evolution remain rare. Here, we report the evolution of a tRNA gene set in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The growth defect caused by deleting the single-copy tRNA gene, serCGA, is rapidly compensated by large-scale (45–290 kb) duplications in the chromosome. Each duplication encompasses a second, compensatory tRNA gene (serTGA) and is associated with a rise in tRNA-Ser(UGA) in the mature tRNA pool. We postulate that tRNA-Ser(CGA) elimination increases the translational demand for tRNA-Ser(UGA), a pressure relieved by increasing serTGA copy number. This work demonstrates that tRNA gene sets can evolve through duplication of existing tRNA genes, a phenomenon that may contribute to the presence of multiple, identical tRNA gene copies within genomes.

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.

Intelligent Design of 14-3-3 Docking Proteins Utilizing Synthetic Evolution Artificial Intelligence (SYN-AI)

The ability to write DNA code from scratch will allow for the discovery of new and interesting chemistries as well as allowing the rewiring of cell signal pathways. Herein, we have utilized synthetic evolution artificial intelligence (SYN-AI) to intelligently design a set of 14-3-3 docking genes. SYN-AI engineers synthetic genes utilizing a parental gene as an evolution template. Wherein, evolution is fast-forwarded by transforming template gene sequences to DNA secondary and tertiary codes based upon gene hierarchical structural levels. The DNA secondary code allows identification of genomic building blocks across an orthologous sequence space comprising multiple genomes. Where, the DNA tertiary code allows engineering of supersecondary structures. SYN-AI constructed a library of 10 million genes that was reduced to three structurally functional 14-3-3 docking genes by applying natural selection protocols. Synthetic protein identity was verified utilizing Clustal Omega sequence alignments and Phylogeny.fr phylogenetic analysis. Wherein, we were able to confirm the three-dimensional structure utilizing I-TASSER and protein-ligand interactions utilizing COACH and Cofactor. The conservation of allosteric communications was confirmed utilizing elastic and anisotropic network models. Whereby, we utilized elNemo and ANM2.1 to confirm conservation of the 14-3-3 ζ amphipathic groove. Notably, to the best of our knowledge, we report the first 14-3-3 docking genes to be written from scratch.

Bacterial wobble modifications of NNA-decoding tRNAs

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

The Expansion of Inosine at the Wobble Position of tRNAs, and Its Role in the Evolution of Proteomes

The modification of adenosine to inosine at the first position of transfer RNA (tRNA) anticodons (I34) is widespread among bacteria and eukaryotes. In bacteria, the modification is found in tRNAArg and is catalyzed by tRNA adenosine deaminase A, a homodimeric enzyme. In eukaryotes, I34 is introduced in up to eight different tRNAs by the heterodimeric adenosine deaminase acting on tRNA. This substrate expansion significantly influenced the evolution of eukaryotic genomes in terms of codon usage and tRNA gene composition. However, the selective advantages driving this process remain unclear. Here, we have studied the evolution of I34, tRNA adenosine deaminase A, adenosine deaminase acting on tRNA, and their relevant codons in a large set of bacterial and eukaryotic species. We show that a functional expansion of I34 to tRNAs other than tRNAArg also occurred within bacteria, in a process likely initiated by the emergence of unmodified A34-containing tRNAs. In eukaryotes, we report on a large variability in the use of I34 in protists, in contrast to a more uniform presence in fungi, plans, and animals. Our data support that the eukaryotic expansion of I34-tRNAs was driven by the improvement brought by these tRNAs to the synthesis of proteins highly enriched in certain amino acids.

Can Protein Expression Be Regulated by Modulation of tRNA Modification Profiles?

tRNAs are the central adaptor molecules in translation. Their decoding properties are influenced by post-transcriptional modifications, particularly in the critical anticodon-stem-loop (ASL) region. Synonymous codon choice, also called codon usage bias, affects both translation efficiency and accuracy, and ASL modifications play key roles in both of these processes. In combination with a handful of historical examples, recent studies integrating ribosome profiling, proteomics, codon-usage analyses, and modification quantifications show that levels of tRNA modifications can change under stress, during development, or under specific metabolic conditions and can modulate the expression of specific genes. Deconvoluting the different responses (global or specific) to tRNA modification deficiencies can be difficult because of pleiotropic effects, but, as more cases emerge, it does seem that tRNA modification changes could add another layer of regulation in the transfer of information from DNA to protein.