Permalink
Cannot retrieve contributors at this time
Fetching contributors…
Cannot retrieve contributors at this time
| <?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.1d1 20130915//EN" "JATS-archivearticle1.dtd"><article article-type="research-article" dtd-version="1.1d1" xmlns:xlink="http://www.w3.org/1999/xlink"><front><journal-meta><journal-id journal-id-type="nlm-ta">elife</journal-id><journal-id journal-id-type="hwp">eLife</journal-id><journal-id journal-id-type="publisher-id">eLife</journal-id><journal-title-group><journal-title>eLife</journal-title></journal-title-group><issn publication-format="electronic">2050-084X</issn><publisher><publisher-name>eLife Sciences Publications, Ltd</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">01339</article-id><article-id pub-id-type="doi">10.7554/eLife.01339</article-id><article-categories><subj-group subj-group-type="display-channel"><subject>Research article</subject></subj-group><subj-group subj-group-type="heading"><subject>Cell biology</subject></subj-group><subj-group subj-group-type="heading"><subject>Genomics and evolutionary biology</subject></subj-group></article-categories><title-group><article-title>tRNA genes rapidly change in evolution to meet novel translational demands</article-title></title-group><contrib-group><contrib contrib-type="author" equal-contrib="yes" id="author-7132"><name><surname>Yona</surname><given-names>Avihu H</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="equal-contrib">†</xref><xref ref-type="fn" rid="con1"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" equal-contrib="yes" id="author-7133"><name><surname>Bloom-Ackermann</surname><given-names>Zohar</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="equal-contrib">†</xref><xref ref-type="fn" rid="con3"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" equal-contrib="yes" id="author-7134"><name><surname>Frumkin</surname><given-names>Idan</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="equal-contrib">†</xref><xref ref-type="fn" rid="con2"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-7135"><name><surname>Hanson-Smith</surname><given-names>Victor</given-names></name><xref ref-type="aff" rid="aff2"/><xref ref-type="aff" rid="aff3"/><xref ref-type="fn" rid="con4"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-7136"><name><surname>Charpak-Amikam</surname><given-names>Yoav</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="con5"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-7137"><name><surname>Feng</surname><given-names>Qinghua</given-names></name><xref ref-type="aff" rid="aff4"/><xref ref-type="fn" rid="con6"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-6261"><name><surname>Boeke</surname><given-names>Jef D</given-names></name><xref ref-type="aff" rid="aff5"/><xref ref-type="fn" rid="pa1">‡</xref><xref ref-type="fn" rid="con7"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" id="author-7139"><name><surname>Dahan</surname><given-names>Orna</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="fn" rid="con8"/><xref ref-type="fn" rid="conf1"/></contrib><contrib contrib-type="author" corresp="yes" id="author-6809"><name><surname>Pilpel</surname><given-names>Yitzhak</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="corresp" rid="cor1">*</xref><xref ref-type="other" rid="par-1"/><xref ref-type="other" rid="par-2"/><xref ref-type="fn" rid="con9"/><xref ref-type="fn" rid="conf1"/></contrib><aff id="aff1"><institution content-type="dept">Department of Molecular Genetics</institution>, <institution>Weizmann Institute of Science</institution>, <addr-line><named-content content-type="city">Rehovot</named-content></addr-line>, <country>Israel</country></aff><aff id="aff2"><institution content-type="dept">Department of Microbiology</institution>, <institution>University of California, San Francisco</institution>, <addr-line><named-content content-type="city">San Francisco</named-content></addr-line>, <country>United States</country></aff><aff id="aff3"><institution content-type="dept">Department of Immunology</institution>, <institution>University of California, San Francisco</institution>, <addr-line><named-content content-type="city">San Francisco</named-content></addr-line>, <country>United States</country></aff><aff id="aff4"><institution content-type="dept">Department of Pathology</institution>, <institution>University of Washington</institution>, <addr-line><named-content content-type="city">Seattle</named-content></addr-line>, <country>United States</country></aff><aff id="aff5"><institution content-type="dept">Department of Molecular Biology and Genetics</institution>, <institution>The Johns Hopkins University School of Medicine</institution>, <addr-line><named-content content-type="city">Baltimore</named-content></addr-line>, <country>United States</country></aff></contrib-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Laub</surname><given-names>Michael</given-names></name><role>Reviewing editor</role><aff><institution>Massachusetts Institute of Technology</institution>, <country>United States</country></aff></contrib></contrib-group><author-notes><corresp id="cor1"><label>*</label>For correspondence: <email>Pilpel@weizmann.ac.il</email></corresp><fn fn-type="con" id="equal-contrib"><label>†</label><p>These authors contributed equally to this work</p></fn><fn fn-type="present-address" id="pa1"><label>‡</label><p>New York University Langone Medical Center, New York, United States</p></fn></author-notes><pub-date date-type="pub" publication-format="electronic"><day>20</day><month>12</month><year>2013</year></pub-date><pub-date pub-type="collection"><year>2013</year></pub-date><volume>2</volume><elocation-id>e01339</elocation-id><history><date date-type="received"><day>05</day><month>08</month><year>2013</year></date><date date-type="accepted"><day>30</day><month>10</month><year>2013</year></date></history><permissions><copyright-statement>© 2013, Yona et al</copyright-statement><copyright-year>2013</copyright-year><copyright-holder>Yona et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/3.0/"><license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/3.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="elife01339.pdf"/><abstract><object-id pub-id-type="doi">10.7554/eLife.01339.001</object-id><p>Changes in expression patterns may occur when organisms are presented with new environmental challenges, for example following migration or genetic changes. To elucidate the mechanisms by which the translational machinery adapts to such changes, we perturbed the tRNA pool of <italic>Saccharomyces cerevisiae</italic> by tRNA gene deletion. We then evolved the deletion strain and observed that the genetic adaptation was recurrently based on a strategic mutation that changed the anticodon of other tRNA genes to match that of the deleted one. Strikingly, a systematic search in hundreds of genomes revealed that anticodon mutations occur throughout the tree of life. We further show that the evolution of the tRNA pool also depends on the need to properly couple translation to protein folding. Together, our observations shed light on the evolution of the tRNA pool, demonstrating that mutation in the anticodons of tRNA genes is a common adaptive mechanism when meeting new translational demands.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.001">http://dx.doi.org/10.7554/eLife.01339.001</ext-link></p></abstract><abstract abstract-type="executive-summary"><object-id pub-id-type="doi">10.7554/eLife.01339.002</object-id><title>eLife digest</title><p>Genes contain the blueprints for the proteins that are essential for countless biological functions and processes, and the path that leads from a particular gene to the corresponding protein is long and complex. The genetic information stored in the DNA must first be transcribed to produce a messenger RNA molecule, which then has to be translated to produce a string of amino acids that fold to form a protein. The translation step is performed by a molecular machine called the ribosome, with transfer RNA molecules bringing the amino acids that are needed to make the protein.</p><p>The information in messenger RNA is stored as a series of letters, with groups of three letters called codons representing the different amino acids. Since there are four letters—A, C, G and U—it is possible to form 64 different codons. And since there are only 20 amino acids, two or more different codons can specify the same amino acid (for example, AGU and AGC both specify serine), and two or more different transfer RNA molecules can take this amino acid to the ribosome. Moreover, some codons are found more often than others in the messenger RNA molecules, so the genes that encode the related transfer RNA molecules are more common than the genes for other transfer RNA molecules.</p><p>Environmental pressures mean that organisms must adapt to survive, with some genes and proteins increasing in importance, and others becoming less important. Clearly the relative numbers of the different transfer RNA molecules will also need to change to reflect these evolutionary changes, but the details of how this happens were not understood.</p><p>Now Yona et al. have explored this issue by studying yeast cells that lack a gene for one of the less common transfer RNA molecules (corresponding to the codon AGG, which specifies the amino acid arginine). At first this mutation resulted in slower growth of the yeast cells, but after being allowed to evolve over 200 generations, the rate of growth matched that of a normal strain with all transfer RNA genes. Yona et al. found that the gene for a more common transfer RNA molecule, corresponding to the codon AGA, which also specifies arginine, had mutated to AGG. As a result, the mutated yeast was eventually able to produce proteins as quickly as wild type yeast. Moreover, further experiments showed that the levels of some transfer RNAs are kept deliberately low in order to slow down the production of proteins so as to ensure that the proteins assume their correct structure.</p><p>But does the way these cells evolved in the lab resemble what happened in nature? To address this question Yona et al. examined a database of transfer RNA sequences from more than 500 species, and found evidence for the same codon-based switching mechanism in many species across the tree of life.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.002">http://dx.doi.org/10.7554/eLife.01339.002</ext-link></p></abstract><kwd-group kwd-group-type="author-keywords"><title>Author keywords</title><kwd>tRNA genes</kwd><kwd>anticodon</kwd><kwd>translation</kwd><kwd>evolution</kwd><kwd>protein folding stress</kwd></kwd-group><kwd-group kwd-group-type="research-organism"><title>Research organism</title><kwd><italic>S. cerevisiae</italic></kwd></kwd-group><funding-group><award-group id="par-1"><funding-source><institution-wrap><institution>European Research Council</institution></institution-wrap></funding-source><award-id>ERC-2007-StG 205199-ERNBPTC</award-id><principal-award-recipient><name><surname>Pilpel</surname><given-names>Yitzhak</given-names></name></principal-award-recipient></award-group><award-group id="par-2"><funding-source><institution-wrap><institution>Ben-May Charitable Trust</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Pilpel</surname><given-names>Yitzhak</given-names></name></principal-award-recipient></award-group><funding-statement>The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.</funding-statement></funding-group><custom-meta-group><custom-meta><meta-name>elife-xml-version</meta-name><meta-value>2</meta-value></custom-meta><custom-meta specific-use="meta-only"><meta-name>Author impact statement</meta-name><meta-value>Experimental evolution and systematic sequence analysis of transfer RNA genes reveal that anticodon mutations provide adaptive plasticity to the translation machinery.</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="s1" sec-type="intro"><title>Introduction</title><p>The process of gene translation is fundamental to the function of living cells, and as such its apparatus is highly conserved across the tree of life (<xref ref-type="bibr" rid="bib40">Müller and Wittmann-Liebold, 1997</xref>; <xref ref-type="bibr" rid="bib26">Itoh et al., 1999</xref>; <xref ref-type="bibr" rid="bib72">Wolf et al., 2001</xref>). Yet, the capacity of the translation machinery to adaptively evolve is crucial in order to support life in changing environments. Therefore, a key open question is to identify the mechanisms by which the translation machinery adapts to changing conditions.</p><p>A thoroughly studied aspect of translation that demonstrates its adaptation capacities is the different proportions by which synonymous codons are used, a phenomenon known as ‘codon usage bias’. Although differential use of codons can be the result of neutral processes such as mutational biases and the genomic GC content (<xref ref-type="bibr" rid="bib66">Urrutia and Hurst, 2001</xref>; <xref ref-type="bibr" rid="bib48">Rao et al., 2011</xref>), natural selection also influences codon usage bias. Indeed, it has been demonstrated that codon choice affects expression level, protein folding, translational accuracy, and other translational features (<xref ref-type="bibr" rid="bib1">Akashi, 1994</xref>; <xref ref-type="bibr" rid="bib43">Parmley and Hurst, 2007</xref>; <xref ref-type="bibr" rid="bib74">Zhou et al., 2009</xref>; <xref ref-type="bibr" rid="bib25">Hudson et al., 2011</xref>). Since both neutral and selective processes govern codon usage bias, the balance between selection, mutational bias and drift is crucial in shaping the codon usage of each species (<xref ref-type="bibr" rid="bib9">Bulmer, 1991</xref>). Importantly, although the selective advantage offered by alternative synonymous codons is considered to be moderate, it was recently demonstrated that selection can still shape codon usage patterns in vertebrates even with their small effective population sizes (<xref ref-type="bibr" rid="bib17">Doherty and McInerney, 2013</xref>).</p><p>Notably, the differential usage of codons represents the evolution of the ‘demand’ aspect of translation, namely the codon usage of all expressed genes. Yet, the adaptation mechanisms of the ‘supply’, namely the expression level of each tRNA type that is loaded with an amino acid, are not fully understood. While ribosomal genes do not exhibit appreciable changes in response to environmental alterations (<xref ref-type="bibr" rid="bib40">Müller and Wittmann-Liebold, 1997</xref>; <xref ref-type="bibr" rid="bib26">Itoh et al., 1999</xref>; <xref ref-type="bibr" rid="bib72">Wolf et al., 2001</xref>), tRNA genes may provide an important source of evolutionary plasticity for fine tuning translation.</p><p>tRNAs constitute a fundamental component in the process of translation, linking codons to their corresponding amino acids (<xref ref-type="bibr" rid="bib69">Widmann et al., 2010</xref>). tRNA genes are classified into gene families according to their anticodon, with each gene family containing between one and several copies scattered throughout the genome. Importantly, it has been experimentally observed for <italic>Saccharomyces cerevisiae</italic> (<xref ref-type="bibr" rid="bib65">Tuller et al., 2010</xref>) and <italic>Escherichia coli</italic> (<xref ref-type="bibr" rid="bib18">Dong et al., 1996</xref>) that the cellular concentrations of each tRNA family in the cell (i.e., the tRNA pool) correlate with its genomic tRNA copy number (<xref ref-type="bibr" rid="bib47">Percudani et al., 1997</xref>; <xref ref-type="bibr" rid="bib28">Kanaya et al., 1999</xref>). Notably, the rate-limiting step of polypeptide elongation is the recruitment of a tRNA that matches the translated codon (<xref ref-type="bibr" rid="bib67">Varenne et al., 1984</xref>). Thus, the translation efficiency is defined as the extent to which the tRNA pool can accommodate the transcriptome (<xref ref-type="bibr" rid="bib58">Sharp and Li, 1987</xref>; <xref ref-type="bibr" rid="bib19">Dos Reis et al., 2004</xref>; <xref ref-type="bibr" rid="bib60">Stoletzki and Eyre-Walker, 2007</xref>), thereby affecting protein production and accuracy.</p><p>In general, highly expressed genes exhibit a marked codon usage bias toward ‘optimal’ codons, whose corresponding tRNA gene copy number is high (<xref ref-type="bibr" rid="bib56">Sharp and Li, 1986a</xref>, <xref ref-type="bibr" rid="bib57">1986b</xref>). The evolutionary force that acts to maintain optimal translation efficiency of such genes was coined ‘translational selection’ (<xref ref-type="bibr" rid="bib19">Dos Reis et al., 2004</xref>). It was previously suggested that translational selection acts to maintain a balance between codon usage and tRNA availability. On the one hand, there is a selective pressure to increase the frequency of preferred codons in highly expressed genes. On the other hand, changes in the tRNA pool may also occur, for example duplication of tRNA genes for which high codon demand exists. Thus, codon frequencies and tRNA copy numbers coevolve toward a supply versus demand balance that facilitates optimal protein production (<xref ref-type="bibr" rid="bib24">Higgs and Ran, 2008</xref>; <xref ref-type="bibr" rid="bib22">Gingold et al., 2012</xref>).</p><p>The fitness effects of an unmet translational demand and its potential role in shaping the tRNA pool are not fully characterized. Evolutionary changes to the tRNA pool were appreciated mainly via bioinformatics studies (<xref ref-type="bibr" rid="bib49">Rawlings et al., 2003</xref>; <xref ref-type="bibr" rid="bib71">Withers et al., 2006</xref>; <xref ref-type="bibr" rid="bib24">Higgs and Ran, 2008</xref>; <xref ref-type="bibr" rid="bib6">Bermudez-Santana et al., 2010</xref>; <xref ref-type="bibr" rid="bib50">Rogers et al., 2010</xref>) and only a handful of experimental findings have been reported, which rely on genetic manipulations (<xref ref-type="bibr" rid="bib10">Byström and Fink, 1989</xref>; <xref ref-type="bibr" rid="bib68">Von Pawel-Rammingen et al., 1992</xref>; <xref ref-type="bibr" rid="bib4">Aström et al., 1993</xref>) or direct mutagenesis (<xref ref-type="bibr" rid="bib53">Saks et al., 1998</xref>). Sequence analyses of divergent genomes have demonstrated that both the sequence and copy number of tRNA genes may change among various species or strains. However, it is unclear whether the observed variations in the tRNA pool are a consequence of an adaptive process due to unbalanced translational demand or the result of random genomic processes, as tRNA genes are a known source of genomic instability (<xref ref-type="bibr" rid="bib39">McFarlane and Whitehall, 2009</xref>).</p><p>Further, the forces that direct and maintain low copy tRNA families remain unclear. Specifically, it is not clear whether translational selection acts only to favor optimal codons or also acts in particular cases to keep other codons deliberately as ‘non-optimal’ by maintaining their tRNA supply at low level. Encoding genes with optimal codons might not always lead to higher protein expression levels (<xref ref-type="bibr" rid="bib34">Kudla et al., 2009</xref>). Similarly, the use of ‘slow codons’ may not always result in lower levels of protein expression as they could have functional roles in improving expression, for example when enriched at the beginning of a transcript in order to reduce the energy of the RNA structure (<xref ref-type="bibr" rid="bib23">Goodman et al., 2013</xref>) or to efficiently allocate ribosomes along the mRNA (<xref ref-type="bibr" rid="bib65">Tuller et al., 2010</xref>). Additionally, it has been proposed that non-optimal codons may play a role in governing the process of cotranslational folding by slowing down translation, which supports proper folding between domain boundaries (<xref ref-type="bibr" rid="bib63">Thanaraj, 1996</xref>; <xref ref-type="bibr" rid="bib33">Kramer et al., 2009</xref>; <xref ref-type="bibr" rid="bib11">Cabrita et al., 2010</xref>; <xref ref-type="bibr" rid="bib70">Wilke and Drummond, 2010</xref>; <xref ref-type="bibr" rid="bib46">Pechmann and Frydman, 2012</xref>). Yet, the contribution of non-optimal codons to proper protein folding was observed only for specific genes (<xref ref-type="bibr" rid="bib16">Crombie et al., 1992</xref>; <xref ref-type="bibr" rid="bib31">Komar et al., 1999</xref>; <xref ref-type="bibr" rid="bib15">Cortazzo et al., 2002</xref>; <xref ref-type="bibr" rid="bib64">Tsai et al., 2008</xref>; <xref ref-type="bibr" rid="bib73">Zhang et al., 2009</xref>; <xref ref-type="bibr" rid="bib75">Zhou et al., 2013</xref>). Thus, the extent and relevance of this phenomenon to the global folding state of the proteome needs to be substantiated.</p><p>To elucidate the importance of restoring translational equilibrium, we used an experimental evolution approach. To this end, we genetically perturbed the tRNA pool of the budding yeast <italic>S. cerevisiae</italic>. In this yeast, the genetic code is decoded by 42 different tRNA families that make up a total of 274 tRNA genes (<xref ref-type="bibr" rid="bib12">Chan and Lowe, 2009</xref>). Each tRNA gene family ranges from 1–16 copies, with 6 tRNA families consisting of only a single copy. In a recent study (<xref ref-type="bibr" rid="bib7">Bloom-Ackermann et al., In press</xref>), we have systematically manipulated the tRNA pool in <italic>S. cerevisiae</italic> by individually deleting most tRNA genes from its genome. Here, we focus on one particular deletion strain that showed the most extreme fitness reduction among the viable deletion mutants in this tRNA deletion library. This tRNA exists in only one copy in the genome, thus after its deletion the cell is left without a tRNA with the corresponding anticodon. Lab-evolution experiments performed on this strain demonstrated that the translational balance was rapidly restored by mutations in other tRNA genes that compensated for the tRNA deletion. An extensive bioinformatic analysis revealed that a similar evolutionary trend is widespread in nature too, suggesting that the anticodon mutations we observed in the lab recapitulate an existing mechanism that shapes the tRNA pool. To shed light on the constraints that shape the size of tRNA gene families, we artificially overexpressed singleton tRNAs, rather than deleting them. We found that when low copy tRNAs were overexpressed, the protein quality control machinery was challenged due to increased proteotoxic stress. This observation suggests that low tRNA availability for particular codon can serve an essential means to ensure proper cotranslation folding of proteins.</p></sec><sec id="s2" sec-type="results"><title>Results</title><sec id="s2-1"><title>Deletion of singleton tRNA gene breaks the translational balance</title><p>To demonstrate the importance of the balance between codon usage and the cellular tRNA pool we created a yeast strain in which the single copy of the arginine tRNA gene, <italic>tR(CCU)J</italic>, was deleted (designated ΔtRNA<sup>Arg</sup><sub>CCU</sub>). Consequently, in this deletion strain, the arginine codon AGG cannot be translated with its fully matched tRNA and it is presumably translated by another arginine tRNA, tRNA<sup>Arg</sup><sub>UCU</sub>, owing to a wobble interaction (<xref ref-type="bibr" rid="bib5">Begley et al., 2007</xref>). This shortage in tRNA supply is particularly evident given the demand: AGG is the second most highly used codon for arginine in the yeast genome (<xref ref-type="supplementary-material" rid="SD2-data">Supplementary file 1A</xref>). Indeed, the ΔtRNA<sup>Arg</sup><sub>CCU</sub> strain showed a severe growth defect compared to the wild-type strain (<xref ref-type="fig" rid="fig1">Figure 1A</xref>, blue and green curves, respectively). This growth difference demonstrates the effect of translational imbalance on cellular growth. Although the deletion mutant of this single copy tRNA is viable (<xref ref-type="bibr" rid="bib13">Clare et al., 1988</xref>; <xref ref-type="bibr" rid="bib29">Kawakami et al., 1993</xref>), its severe growth defect also reveals the inability of the wobble interactions to fully compensate for the tRNA gene deletion.<fig id="fig1" position="float"><object-id pub-id-type="doi">10.7554/eLife.01339.003</object-id><label>Figure 1.</label><caption><title>The growth defect associated with deletion of a singleton tRNA gene was rapidly rescued during the lab-evolution experiment.</title><p>(<bold>A</bold>) Growth curve measurements of wild-type (WT) (green), ΔtRNA<sup>Arg</sup><sub>CCU</sub> (blue) and the evolved deletion (red) are shown in optical density (OD) values over time during continuous growth on rich medium at 30°C. (<bold>B</bold>) The mutation that was found to recover the deletion phenotype in the evolved strains is shown on the secondary structure of tRNA<sup>Arg</sup><sub>UCU.</sub> The UCU anticodon nucleotides are marked with black circles, and the red circle indicates the mutation that occurred in the anticodon, that is T→C transition. (<bold>C</bold>) <italic>MutΔtRNA</italic><sup><italic>Arg</italic></sup><sub><italic>CCU</italic></sub> in which the same mutation that was found in the evolved strain was deliberately engineered, exhibits similar growth as the WT. Growth curve measurements of WT (green) and of <italic>MutΔtRNA</italic><sup><italic>Arg</italic></sup><sub><italic>CCU</italic></sub> (magenta) are shown in OD<sub>600</sub> values over time during continuous growth on rich medium at 30°C.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.003">http://dx.doi.org/10.7554/eLife.01339.003</ext-link></p></caption><graphic xlink:href="elife01339f001"/></fig></p></sec><sec id="s2-2"><title>The tRNA pool can rapidly evolve to meet translational demands</title><p>To learn how genomes adapt to translational imbalances, we performed lab-evolution experiments on the ΔtRNA<sup>Arg</sup><sub>CCU</sub> strain, employing the procedure of daily growth and dilution to a fresh medium (<xref ref-type="bibr" rid="bib35">Lenski et al., 1991</xref>). The deletion strain was grown under optimal laboratory conditions (rich medium at 30°C) and was diluted daily into a fresh medium by a factor of 120, corresponding to approximately 7 generations per cycle. Every 50 generations, the growth of the evolving population was compared to both the wild-type and the ancestor ΔtRNA<sup>Arg</sup><sub>CCU</sub> strains. Strikingly, after 200 generations we observed a full recovery of the growth defect of the ancestor strain ΔtRNA<sup>Arg</sup><sub>CCU</sub>, as the growth curve of the evolved population was indistinguishable from that of the wild-type strain (<xref ref-type="fig" rid="fig1">Figure 1A</xref>, red curve). Similar dynamics were observed in all four independent evolutionary lines.</p><p>In search of the potential genetic adaptations underlying this rapid recovery, we first looked for genetic alterations in other arginine tRNA genes. We found a single point mutation in another arginine tRNA gene that codes for tRNA<sup>Arg</sup><sub>UCU</sub>. This mutation changed the anticodon triplet of tRNA<sup>Arg</sup><sub>UCU</sub> from UCU to CCU (i.e., T→C transition). Consequently, the evolved tRNA<sup>Arg</sup><sub>UCU</sub> perfectly matched the AGG codon (<xref ref-type="fig" rid="fig1">Figure 1B</xref>). Unlike the singleton tRNA<sup>Arg</sup><sub>CCU</sub>, there are 11 copies of tRNA<sup>Arg</sup><sub>UCU</sub> in the yeast genome. Although each of the 4 independent lab-evolution experiments showed the exact same solution (that is, a mutation in the anticodon of a tRNA<sup>Arg</sup><sub>UCU</sub> gene), 3 different copies of this gene were changed among the 4 lines (i.e., 1 of the 11 copies of tRNA<sup>Arg</sup><sub>UCU</sub> was mutated in 2 repetitions; see ‘Materials and methods’). To confirm that a single point mutation in the anticodon of tRNA<sup>Arg</sup><sub>UCU</sub> is sufficient to fully compensate for the growth defect of ΔtRNA<sup>Arg</sup><sub>CCU</sub>, we artificially inserted the same T→C mutation into the deletion ΔtRNA<sup>Arg</sup><sub>CCU</sub> mutant. We inserted the mutation into 1 of the 11 copies of the tRNA<sup>Arg</sup><sub>UCU</sub> genes, a copy that resides on chromosome XI, and thus spontaneously mutated in 1 of the evolution lines. Indeed, the artificially mutated strain, termed as <italic>MutΔtRNA</italic><sup><italic>Arg</italic></sup><sub><italic>CCU</italic></sub>, showed a full recovery of the deletion adverse phenotype (<xref ref-type="fig" rid="fig1">Figure 1C</xref>). This indicates that the T→C mutation in the anticodon is sufficient for the full recovery of the tRNA<sup>Arg</sup><sub>CCU</sub> deletion phenotype.</p></sec><sec id="s2-3"><title>Mutated tRNA<sup>Arg</sup><sub>UCU</sub> is functional despite sequence dissimilarities with respect to the deleted tRNA<sup>Arg</sup><sub>CCU</sub></title><p>In general, all copies of each tRNA gene family tend to be highly similar in sequence in <italic>S. cerevisiae</italic> (<xref ref-type="bibr" rid="bib12">Chan and Lowe, 2009</xref>). In particular, the sequences of the 11 copies of tRNA<sup>Arg</sup><sub>UCU</sub> are 100% identical to each other. Yet, the 2 arginine tRNA, tRNA<sup>Arg</sup><sub>UCU</sub>, and tRNA<sup>Arg</sup><sub>CCU</sub>, differ in 21 of their 72 nucleotides (including the third anticodon position, <xref ref-type="fig" rid="fig2">Figure 2A</xref>). Thus, the evolutionary solution that occurred in our experiments created a ‘chimeric’ tRNA with a CCU anticodon, whereas the rest of the tRNA sequence (termed as the ‘tRNA scaffold’) remained as tRNA<sup>Arg</sup><sub>UCU</sub>. The sequence identity among all members of the tRNA<sup>Arg</sup><sub>UCU</sub> family suggests a functional role for the tRNA scaffold in addition to that of the anticodon (<xref ref-type="bibr" rid="bib55">Schultz and Yarus, 1994</xref>; <xref ref-type="bibr" rid="bib32">Konevega et al., 2004</xref>; <xref ref-type="bibr" rid="bib14">Cochella and Green, 2005</xref>; <xref ref-type="bibr" rid="bib42">Olejniczak et al., 2005</xref>; <xref ref-type="bibr" rid="bib52">Saks and Conery, 2007</xref>; <xref ref-type="bibr" rid="bib54">Schmeing et al., 2011</xref>). Therefore, it is surprising that the chimeric tRNA performed just as well as the deleted tRNA<sup>Arg</sup><sub>CCU</sub> in terms of cell growth, despite the major sequence differences between the two tRNA scaffolds. Thus, we raised the hypothesis that more challenging growth conditions may expose possible inadequacies in the chimeric tRNA. To test this notion, we compared the rescued strain, <italic>MutΔ</italic>tRNA<sup>Arg</sup><sub>CCU</sub>, which carries the chimeric tRNA, to the wild-type strain under an array of challenging conditions. Surprisingly, under all the examined conditions, we observed no significant growth difference between the two strains (<xref ref-type="fig" rid="fig2">Figure 2B</xref>). Hence, the chimeric tRNA provides a direct in vivo indication that the scaffolds of tRNAs, which encode for the same amino acid, may be interchangeable in terms of their effect on cellular growth under the conditions we tested.<fig-group><fig id="fig2" position="float"><object-id pub-id-type="doi">10.7554/eLife.01339.004</object-id><label>Figure 2.</label><caption><title>The growth of MutΔtRNA<sup>Arg</sup><sub>CCU</sub> carrying the chimeric tRNA compared to wild-type (WT) under different conditions.</title><p>(<bold>A</bold>) The sequence of the chimeric tRNA is drawn showing the scaffold of tRNA<sup>Arg</sup><sub>UCU</sub> with the mutated CCU anticodon. The anticodon triplet is marked with black circles. The evolved mutation is marked with a red circle. All 20 nucleotide differences between tRNA<sup>Arg</sup><sub>UCU</sub> and tRNA<sup>Arg</sup><sub>CCU</sub> are marked with blue background, next to which, in green letters, the original nucleotide of tRNA<sup>Arg</sup><sub>CCU</sub> are written. (<bold>B</bold>) Growth curve measurements of WT (green) and of <italic>MutΔtRNA</italic><sup><italic>Arg</italic></sup><sub><italic>CCU</italic></sub> (magenta) are shown in OD<sub>600</sub> values over time during continuous growth.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.004">http://dx.doi.org/10.7554/eLife.01339.004</ext-link></p></caption><graphic xlink:href="elife01339f002"/></fig><fig id="fig2s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01339.005</object-id><label>Figure 2—figure supplement 1.</label><caption><title>Quadruple deletion of tRNA<sup>ser</sup><sub>GCU</sub> is lethal.</title><p>We perturbed the tRNA pool in a wild-type strain by deletion of an entire serine tRNA family. Here, the supply of tRNA<sup>Ser</sup><sub>GCU</sub> was eliminated by deletion of all four identical gene copies of this family located on chromosomes IV, VI, XII and X. A complete deletion of this gene family is lethal, indicating that the tRNA<sup>Ser</sup><sub>GCU</sub> is essential in <italic>Saccharomyces cerevisiae</italic>. To validate that the lethality is indeed due to the deletion of the tRNA<sup>Ser</sup><sub>GCU</sub> genes and not due to an unintentional perturbation of other putative genetic features in the vicinity of the deleted tRNA<sup>Ser</sup><sub>GCU</sub> copies, we introduced a plasmid with the tRNA<sup>Ser</sup><sub>GCU</sub> gene. As expected, the quadruple deletion strain was viable when supplemented with a plasmid carrying the tRNA<sup>Ser</sup><sub>GCU</sub> gene. (<bold>A</bold> and <bold>C</bold>) BY384 and tRNA<sup>Ser</sup><sub>GCU</sub> quadruple deletion strains, respectively, with a plasmid harboring a tRNA<sup>Ser</sup><sub>GCU</sub> gene grown on yeast extract/peptone/dextrose (YPD) plates. (<bold>B</bold> and <bold>D</bold>) BY384 and tRNA<sup>Ser</sup><sub>GCU</sub> quadruple deletion strains, respectively, grown on liquid medium that contains uracil (YPD) to allow the growth of cells that lost the plasmid and then plated on 5-fluoro-orotic acid (5-FOA) plates (see ‘Materials and methods’). As expected, normal size colonies were only observed in <bold>A</bold>, <bold>B</bold> and <bold>C</bold>. This means that the quadruple deletion strain with the plasmid could not lose the tRNA gene because of its essentiality.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.005">http://dx.doi.org/10.7554/eLife.01339.005</ext-link></p></caption><graphic xlink:href="elife01339fs001"/></fig><fig id="fig2s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01339.006</object-id><label>Figure 2—figure supplement 2.</label><caption><title>A chimeric serine tRNA can rescue the lethality of the quadruple deletion.</title><p>Here, rather than inserting a plasmid with the original tRNA<sup>Ser</sup><sub>GCU</sub> gene, we complemented the tRNA<sup>Ser</sup><sub>GCU</sub> family deletion strain with a plasmid containing a chimeric serine tRNA with a GCU anticodon. The chimeric tRNA has an alternative scaffold derived from a different serine tRNA, tRNA<sup>Ser</sup><sub>CGA</sub>, that differs in 22 positions from the deleted serine tRNA family. This strain formed normal size colonies on a yeast extract/peptone/dextrose (YPD) plate.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.006">http://dx.doi.org/10.7554/eLife.01339.006</ext-link></p></caption><graphic xlink:href="elife01339fs002"/></fig></fig-group></p><p>To examine the generality of our observation, we once again perturbed the tRNA pool in a wild-type (WT) strain by deletion of an entire serine tRNA family, tRNA<sup>Ser</sup><sub>GCU</sub> that has four copies in the genome. A complete deletion of this gene family is lethal, indicating that the tRNA<sup>Ser</sup><sub>GCU</sub> is essential in <italic>S. cerevisiae</italic>. Although we could not evolve that strain, we did find that this quadruple deletion strain was viable when supplemented with a centromeric plasmid carrying the tRNA<sup>Ser</sup><sub>GCU</sub> gene. Thus, the lethality is conferred directly from the tRNA loss and is not due to other indirect effects (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref>). We hypothesized that, as with tRNA<sup>Arg</sup><sub>CCU</sub>, other chimeric serine tRNAs that carry a GCU anticodon, yet with the scaffold of another tRNA for serine, would also prevent the observed lethality. Indeed, a strain carrying a chimeric tRNA with a scaffold of tRNA<sup>Ser</sup><sub>CGA</sub> and the GCU anticodon is viable on the background of the tRNA<sup>Ser</sup><sub>GCU</sub>-family deletion (<xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2</xref>). Therefore, we concluded that the identity of the anticodon is essential for the function of the tRNA<sup>Ser</sup><sub>GCU</sub> gene family. Thus, it appears that for the examined tRNAs the anticodon is a dominant feature in terms of cellular fitness, overshadowing other sequence elements.</p></sec><sec id="s2-4"><title>Anticodon switching is a widespread phenomenon in nature</title><p>Although the anticodon of a tRNA gene was rapidly mutated under our laboratory conditions, thus regaining proper translational equilibrium, it is unclear to what extent this mechanism naturally occurs in species across the tree of life. To address this question, we performed a systematic bioinformatics screen for tRNA switching events in nature. We defined an anticodon-switching event as a case of a tRNA whose nucleotide sequence is closer to a tRNA gene with a different anticodon than to a tRNA gene with the same anticodon. To this end, we downloaded all the known tRNA sequences from the Genomic tRNA Database (<xref ref-type="bibr" rid="bib12">Chan and Lowe, 2009</xref>), a collection that stores the tRNA pools of 524 species. We masked the anticodon triplet as ‘NNN’ in all tRNA genes, aligned all tRNA sequences from each species individually and inferred a maximum likelihood phylogenetic tree for each alignment. For each tRNA sequence, we calculated the shortest phylogenetic distance to another tRNA with the same anticodon (designated d<sub>same</sub>) and the shortest distance to another tRNA with a different anticodon (designated d<sub>diff</sub>). For each species, we defined its set of tRNA switching events as those in which d<sub>diff</sub><d<sub>same</sub> (see ‘Materials and methods’, <xref ref-type="supplementary-material" rid="SD1-data">Figure 3—source data 1</xref>).</p><p>Our analysis included 416 eubacterial, 68 eukaryotic, and 40 archaeal species. We found that tRNA switching events are present in all domains of life, as we detected at least 1 tRNA switching event per species in 8 bacteria, 58 eukarya, and 1 archaeal species (<xref ref-type="fig" rid="fig3">Figure 3A</xref>). A retrospective counting revealed that most switching events occurred due to a mutation in the first position in the anticodon triplet that corresponds to the third codon position (see details in ‘Materials and methods’). For comparison, we masked as ‘NNN’ additional triplets of nucleotides within the tRNA molecule, and found a higher percentage of discrepancies compared to the anticodon triplet (<xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1</xref>; <xref ref-type="supplementary-material" rid="SD2-data">Supplementary file 1B</xref>).<fig-group><fig id="fig3" position="float"><object-id pub-id-type="doi">10.7554/eLife.01339.007</object-id><label>Figure 3.</label><caption><title>Anticodon switching is a widespread phenomenon in nature.</title><p>(<bold>A</bold>) Number of species with at least one tRNA switching event in each domain of life. (<bold>B</bold>) The anticodon UUC convergently evolved in <italic>Mus musculus</italic>. A maximum likelihood phylogeny of tRNA sequences in <italic>M. musculus</italic> that decode glutamic acid (Glu) codons. Branch lengths express average nucleotide substitutions per site. Decimals on internal branches express branch support. (<bold>C</bold>) A comparison of nucleotide sequences for glutamic acid tRNA genes in <italic>M. musculus</italic> with anticodon UUC (top, tRNA1547 and tRNA359), ‘switched’ UUC tRNAs (middle, tRNA286 and tRNA754), and CUC tRNAs (bottom, tRNA1002, tRNA745, tRNA303, tRNA999, tRNA996 tRNA709, tRNA1001, tRNA1912 and tRNA81). The anticodon triplet is boxed in gray. Red vertical bars indicate differences between sequences. (<bold>D</bold>) The anticodon UAC convergently evolved in <italic>Homo sapiens</italic>. A maximum likelihood phylogeny of tRNA sequences in <italic>H. sapiens</italic> encoding for valine (Val) is shown. (<bold>E</bold>) A comparison of nucleotide sequences for <italic>H. sapiens</italic> tRNAs with anticodons UAC (top, tRNA6), a ‘switched’ UAC tRNA (middle, tRNA40), and an AAC tRNA (bottom, tRNA136). The number of genes is according to the tRNA database.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.007">http://dx.doi.org/10.7554/eLife.01339.007</ext-link></p><p><supplementary-material id="SD1-data"><object-id pub-id-type="doi">10.7554/eLife.01339.008</object-id><label>Figure 3—source data 1.</label><caption><title>Table of anticodon switchings in different species across the tree of life.</title><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.008">http://dx.doi.org/10.7554/eLife.01339.008</ext-link></p></caption><media mime-subtype="xlsx" mimetype="application" xlink:href="elife01339s001.xlsx"/></supplementary-material></p></caption><graphic xlink:href="elife01339f003"/></fig><fig id="fig3s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01339.009</object-id><label>Figure 3—figure supplement 1.</label><caption><title>A comparison of discrepancy proportions at the anticodon triplet vs control triplets.</title><p>Six control triplets were chosen; their locations are marked on a tRNA structure. (<bold>A</bold>–<bold>F</bold>) A comparison of the proportion of switched tRNAs for each species when the anticodon was masked versus the proportion of alternative discrepancies when the control triplet was masked. Each point corresponds to one species. A tRNA was identified as ‘switched’ if its nucleotide sequence phylogenetically clustered with other tRNAs with dissimilar mask triplets (see ‘Materials and methods’).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.009">http://dx.doi.org/10.7554/eLife.01339.009</ext-link></p></caption><graphic xlink:href="elife01339fs003"/></fig></fig-group></p><p><xref ref-type="fig" rid="fig3">Figure 3</xref> demonstrates two examples of tRNA switching events, the first in <italic>Mus musculus</italic> and the second in <italic>Homo sapiens</italic>. In the first example, the phylogeny of tRNA sequences with glutamic acid anticodons is presented (<xref ref-type="fig" rid="fig3">Figure 3B</xref>). Notably, six out of the eight tRNAs with a UUC anticodon in <italic>M. musculus</italic> were clustered together in our analysis, while two other copies of the same anticodon identity were clustered closer to tRNA genes with a CUC anticodon (<xref ref-type="fig" rid="fig3">Figure 3C</xref>). The second example demonstrates a switching event for tRNA genes encoding for valine anticodons. In this study, a tRNA with a UAC anticodon was clustered with CAC and AAC tRNA genes and not with the other four UAC tRNAs (<xref ref-type="fig" rid="fig3">Figure 3D,E</xref>). Interestingly, the CAC and AAC tRNA genes are intermixed in the tree, suggesting that anticodon switching was prevalent in the evolution of CAC and AAC tRNA genes in <italic>H. sapiens</italic> (<xref ref-type="fig" rid="fig3">Figure 3D,E</xref>). Also of interest, the switching events shown in mouse were not found in human and vice versa. Thus, in each of these two mammals the switching examples shown here probably occurred after they split from their common ancestor. In general, inspecting the relationship across species between the size of the tRNA pool and the number of detected switching events revealed a modest correlation, and in particular species with same size of tRNA repertoire manifested tRNA switching to different extents (not shown). This analysis suggests that future examination of the tRNA switching phenomenon in individual species could be of interest.</p></sec><sec id="s2-5"><title>Multiple copies of rare tRNAs are deleterious to the cell</title><p>After demonstrating the prevalence of anticodon switching, we refocused on our lab-evolution results. The switching events that we observed (from tRNA<sup>Arg</sup><sub>UCU</sub> to tRNA<sup>Arg</sup><sub>CCU</sub>) suggest that the effective gene copy number of each tRNA anticodon set can change during evolution, presumably due to demand-to-supply changes. Given that a single point mutation can functionally convert one tRNA into another, an interesting question emerges: why does the genome maintain a single copy of tRNA<sup>Arg</sup><sub>CCU</sub>? T to C mutations must have occurred in evolution but they appear to have been selected against so as to preserve only a single copy of the CCU anticodon tRNA. Consistent with this hypothesis is the observation that other yeast species maintain tRNA<sup>Arg</sup><sub>CCU</sub> at a single copy (<xref ref-type="supplementary-material" rid="SD2-data">Supplementary file 1C</xref>). We thus reasoned that an artificial increase in the copy number of a rare tRNA, but not of an abundant one, might result in a deleterious effect on the cells.</p><p>Indeed, transformation of a multi-copy plasmid containing a tRNA<sup>Arg</sup><sub>CCU</sub> gene to a wild-type strain (termed as <italic>WTmultiCCU</italic>) resulted in a substantial growth reduction when compared to wild-type cells carrying an empty multi-copy plasmid (termed <italic>WTmultiControl</italic>). In contrast, when we created a strain with a similar multi-copy plasmid that contains the abundant tRNA<sup>Arg</sup><sub>UCU</sub> gene, designated as <italic>WTmultiUCU</italic>, a growth profile much closer to that of <italic>WTmultiControl</italic> was observed (<xref ref-type="fig" rid="fig4">Figure 4A</xref>). These findings are consistent with the evolutionary tendency for yeast to keep a low copy number of tRNA<sup>Arg</sup><sub>CCU</sub> and suggest that a high copy number of such rare tRNAs is deleterious to cells.<fig-group><fig id="fig4" position="float"><object-id pub-id-type="doi">10.7554/eLife.01339.010</object-id><label>Figure 4.</label><caption><title><italic>WTmultiCCU</italic> experiences a growth defect compared to <italic>WTmultiUCU</italic> and demonstrates higher levels of misfolded proteins.</title><p>(<bold>A</bold>) Growth curve measurements of <italic>WTmultiControl</italic> (blue), <italic>WTmultiUCU</italic> (brown) and <italic>WTmultiCCU</italic> (khaki) are shown in optical density (OD) values over time during continuous growth. The <italic>WTmultiCCU</italic> strain carrying a high copy number plasmid harboring tRNA<sup>Arg</sup><sub>CCU</sub> demonstrates slower growth compared to cells with an empty plasmid or with a tRNA<sup>Arg</sup><sub>UCU</sub> plasmid that is mainly characterized by a longer growth delay (lag phase). (<bold>B</bold>) A demonstration of a <italic>WTmultiCCU</italic> cell in which the mCherry-Von Hippel–Lindau (VHL) proteins appear with a punctum phenotype when the protein quality control machinery is saturated with misfolded proteins. (<bold>C</bold>) A demonstration <italic>WTmultiUCU</italic> cell in which the quality control machinery is not occupied with other proteins; mCherry-VHL is localized to the cytosol. (<bold>D</bold>) <italic>WTmultiCCU</italic>, <italic>WTmultiUCU</italic> and <italic>WTmultiControl</italic> were transformed with a VHL-mCherry containing plasmid and visualized under the microscope; 1000 cells per strain were counted for either cytosolic or punctum localization of the VHL protein. The fold change in the number of cells containing puncta was then deduced by normalization to the <italic>WTmultiControl</italic> population. The 95% confidence interval is indicated. (<bold>E</bold>) The mRNA fold change of six representative heat-shock genes measured by real-time quantitative PCR (RT-qPCR). Presented values are the mean of two biological repetitions ± SEM. The significance of the fold change differences was examined using a <italic>t</italic> test, with *p<0.001 or **p<0.0001.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.010">http://dx.doi.org/10.7554/eLife.01339.010</ext-link></p></caption><graphic xlink:href="elife01339f004"/></fig><fig id="fig4s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01339.011</object-id><label>Figure 4—figure supplement 1.</label><caption><title>Multiple copies of rare tRNA<sup>Ser</sup><sub>CGA</sub> gene are deleterious compared to abundant tRNA<sup>Ser</sup><sub>AGA</sub>.</title><p>Growth curve measurements of <italic>WTmultiControl</italic> (blue), <italic>WTmultiAGA</italic> (brown) and <italic>WTmultiCGA</italic> (khaki) are shown in optical density (OD) values over time during continuous growth on rich medium at 30°C. The <italic>WTmultiCGA</italic> strain carrying a high copy number plasmid harboring tRNA<sup>Ser</sup><sub>CGA</sub> demonstrates slower growth compared to cells with an empty plasmid or with a tRNA<sup>Ser</sup><sub>AGA</sub> plasmid.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.011">http://dx.doi.org/10.7554/eLife.01339.011</ext-link></p></caption><graphic xlink:href="elife01339fs004"/></fig><fig id="fig4s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01339.012</object-id><label>Figure 4—figure supplement 2.</label><caption><title>Multiple copies of the rare tRNA<sup>Gln</sup><sub>CUG</sub> gene are deleterious compared to abundant tRNA<sup>Gln</sup><sub>UUG</sub>.</title><p>Growth curve measurements of <italic>WTmultiControl</italic> (blue), <italic>WTmultiUUG</italic> (brown) and <italic>WTmultiCUG</italic> (khaki) are shown in optical density (OD) values over time during continuous growth on rich medium at 30°C. The <italic>WTmultiCUG</italic> strain carrying a high copy number plasmid harboring tRNA<sup>Gln</sup><sub>CUG</sub> demonstrates slower growth compared to cells with an empty plasmid or with a tRNA<sup>Gln</sup><sub>UUG</sub> plasmid.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.012">http://dx.doi.org/10.7554/eLife.01339.012</ext-link></p></caption><graphic xlink:href="elife01339fs005"/></fig><fig id="fig4s3" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.01339.013</object-id><label>Figure 4—figure supplement 3.</label><caption><title>Addition of low copy number tRNA<sup>Arg</sup><sub>CCU</sub> is deleterious compared to low copy number tRNA<sup>Arg</sup><sub>UCU</sub> when grown in heat.</title><p>We created two strains harboring low copy number (CEN) plasmids carrying either tRNA<sup>Arg</sup><sub>UCU</sub> or tRNA<sup>Arg</sup><sub>CCU</sub> in addition to the endogenous copy. We termed these strains <italic>WTextraUCU</italic> and <italic>WTextraCCU</italic>, respectively. Growth curve measurements of <italic>WTextraUCU</italic> (green) and <italic>WTextraCCU</italic> (red) are shown in optical density (OD) values over time during continuous growth on rich medium at 39°C. Cells with the centromeric tRNA<sup>Arg</sup><sub>CCU</sub> plasmid showed a modest growth defect compared to cells with the centromeric tRNA<sup>Arg</sup><sub>UCU</sub> plasmid.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.013">http://dx.doi.org/10.7554/eLife.01339.013</ext-link></p></caption><graphic xlink:href="elife01339fs006"/></fig></fig-group></p><p>To demonstrate the generality of our findings, we employed the same assays in two additional cases. First, we examined 2 serine tRNAs, the singleton tRNA<sup>Ser</sup><sub>CGA</sub> and tRNA<sup>Ser</sup><sub>AGA</sub> that is found in the genome in 11 copies. In the second case, we focused on two glutamine tRNAs, the singleton tRNA<sup>Gln</sup><sub>CUG</sub> and tRNA<sup>Gln</sup><sub>UUG</sub> that is found in the genome in nine copies. In both the cases, we observed that the wild-type strain supplemented with multiple copies of a singleton tRNA exhibit impaired growth compared to the same strain supplemented with the abundant tRNA for the same amino acid (<xref ref-type="fig" rid="fig4s1 fig4s2">Figure 4—figure supplements 1 and 2</xref>). Since the changes in tRNA family sizes during evolution likely occur gradually, perhaps one copy at a time, we also examined the effect of adding low copy number plasmids carrying either tRNA<sup>Arg</sup><sub>UCU</sub> or tRNA<sup>Arg</sup><sub>CCU</sub>. The cells with the tRNA<sup>Arg</sup><sub>CCU</sub> plasmid showed a modest growth defect compared to the cells with tRNA<sup>Arg</sup><sub>UCU</sub> plasmid, yet only when grown at 39°C (<xref ref-type="fig" rid="fig4s3">Figure 4—figure supplement 3</xref>).</p></sec><sec id="s2-6"><title>Multiple copies of the rare tRNA<sup>Arg</sup><sub>CCU</sub> induce proteotoxic stress</title><p>Why is it essential to keep certain tRNAs at a low level? One interesting possibility is that rare tRNAs are essential for the process of cotranslation folding, presumably because low abundance tRNAs provide a pause in translation that might be needed for proper folding (<xref ref-type="bibr" rid="bib63">Thanaraj, 1996</xref>; <xref ref-type="bibr" rid="bib20">Drummond and Wilke, 2008</xref>; <xref ref-type="bibr" rid="bib11">Cabrita et al., 2010</xref>; <xref ref-type="bibr" rid="bib46">Pechmann and Frydman, 2012</xref>). Other deleterious effects that may stem from a high copy number of tRNA<sup>Arg</sup><sub>CCU</sub> could be misincorporation of arginine into non-arginine codons, or the misloading of arginine tRNA molecules with other amino acids. These potential sources of errors are not mutually exclusive and can each contribute to the observed growth defect by exerting a protein folding stress.</p><p>To examine the possibility that the growth defect associated with multiple copies of tRNA<sup>Arg</sup><sub>CCU</sub> is indeed associated with such a proteotoxic stress, we used an established method that examines the load on the protein quality control machinery of the cell (see ‘Materials and methods’) (<xref ref-type="bibr" rid="bib27">Kaganovich et al., 2008</xref>). In this assay, we transformed cells with a plasmid harboring the human gene, von Hippel–Lindau (VHL), fused to a fluorescent tag (mCherry). Fluorescently tagged VHL that is present as aggregated puncta (<xref ref-type="fig" rid="fig4">Figure 4B</xref>), and not as a disperse cytosolic localization (<xref ref-type="fig" rid="fig4">Figure 4C</xref>), indicates that the protein quality control machinery is saturated due to high levels of misfolding in the cell’s endogenous proteins. We transformed the VHL-mCherry plasmid to each of the multi-copy tRNA strains, <italic>WTmultiCCU</italic>, <italic>WTmultiUCU</italic> and <italic>WTmultiControl,</italic> and monitored the level of proteotoxic stress by quantifying the number of cells with puncta phenotype in each population. The fold change in those cells was then deduced by normalization to the <italic>WTmultiControl</italic> population. We found that while <italic>WTmultiUCU</italic> exhibited similar amount of cells with puncta as the <italic>WTmultiControl</italic>, the <italic>WTmultiCCU</italic> exhibited a threefold increase (<xref ref-type="fig" rid="fig4">Figure 4D</xref>).</p><p>The proteotoxic stress experienced by the two strains overexpressing tRNAs was further assessed by measuring the induction level of an array of heat-shock proteins (HSPs) using real-time quantitative PCR (RT-qPCR). Since the HSPs have been shown to undergo induction under proteotoxic stress, they are an excellent indicator for this stress (<xref ref-type="bibr" rid="bib38">McClellan et al., 2005</xref>; <xref ref-type="bibr" rid="bib27">Kaganovich et al., 2008</xref>). Indeed, we found a significant upregulation in mRNA levels for all the examined HSP genes in the <italic>WTmultiCCU</italic> strain, but not in the <italic>WTmultiUCU</italic> strain (<xref ref-type="fig" rid="fig4">Figure 4E</xref>). These findings further demonstrate that increasing the copy number of a rare tRNA gene, but not of an already abundant one, results in proteotoxic stress in the cell.</p></sec></sec><sec id="s3" sec-type="discussion"><title>Discussion</title><p>Genomic duplications, deletions, and anticodon mutations shape tRNA gene families, yet the evolutionary scenarios that trigger changes in the tRNA pool have not been thoroughly explored. In our evolution experiments, a translational imbalance was imposed by a tRNA gene deletion that compromised growth and drove the tRNA pool to adapt to a novel translational demand. Importantly, organisms may experience equivalent imbalances when their gene expression changes due to altered environmental conditions or upon migrating to a new ecological niche (<xref ref-type="bibr" rid="bib22">Gingold et al., 2012</xref>). This scenario is particularly feasible given that the genes needed in various environments do show differences in codon usage, for example respiration as opposed to fermentation in yeast (<xref ref-type="bibr" rid="bib37">Man and Pilpel, 2007</xref>).</p><p>Indeed, when faced with different environmental challenges, transcriptional changes affect the codon usage of the transcriptome (<xref ref-type="bibr" rid="bib22">Gingold et al., 2012</xref>), and hence the demand for the various tRNAs, and may thus cause translational imbalances. To maintain optimal protein production, the tRNA pool is under pressure to restore the translational balance by accommodating the new translational demands. On a short timescale, the tRNA pool might respond non-genetically by changing expression profiles of the tRNAs (<xref ref-type="bibr" rid="bib65">Tuller et al., 2010</xref>; <xref ref-type="bibr" rid="bib51">Saikia et al., 2012</xref>; <xref ref-type="bibr" rid="bib44">Pavon-Eternod et al., 2013a</xref>, <xref ref-type="bibr" rid="bib45">2013b</xref>). Yet, if changes in demand-to-supply persist, a genetic change in the tRNA pool might become beneficial evolutionarily. In this work, we demonstrate how anticodon mutations provide a rapid mechanism to alter the tRNA pool. We propose that during evolution, novel translational requirements can be addressed by anticodon shifting of tRNA copies more readily than by duplications and deletions of tRNA genes. The tRNA pool can evolve to meet new translation demands by adjusting the ratios of tRNA families that code for the same amino acid. Within a single mutational event, anticodon switching holds the potential to rapidly change the ratios of tRNAs within the pool, by increasing the copy number of one tRNA family at the expense of a counterpart. A similar solution could be obtained by a sequence of genomic duplications and deletions of tRNA genes. These alternatives are likely to fixate less frequently than anticodon switching, as they may carry negative effects due to duplications or deletions of adjacent unrelated genetic features. Furthermore, our systematic search for tRNA switching events throughout the tree of life revealed the prevalence of tRNA anticodon mutations in nature. This observation is consistent with the results of our lab-evolution experiment and may be the evolutionary outcome to novel translational demands in the wild.</p><p>Studies on methionine tRNAs have previously shown that the scaffold sequences determine their function as either initiator or elongator Met tRNA (<xref ref-type="bibr" rid="bib68">Von Pawel-Rammingen et al., 1992</xref>; <xref ref-type="bibr" rid="bib4">Aström et al., 1993</xref>; <xref ref-type="bibr" rid="bib30">Kolitz and Lorsch, 2010</xref>). Yet, the initiator and the elongator represent extreme cases of tRNAs that are used at different stages in translation. The present chimeric tRNAs that emerged in our lab-evolution experiments successfully replaced the deleted tRNA, despite differences in 20 nucleotide positions between the 2 tRNA scaffolds. If tRNA scaffolds are interchangeable in terms of the effect of their function on the fitness, what can explain the high sequence similarity observed among tRNA gene copies of the same family in yeast? It is possible that the sequence of the tRNA scaffold is indeed important under specific conditions that were not examined in this work, or that our measurements were not sensitive enough to detect small selective disadvantages that can act against chimeric tRNAs in nature. Under these scenarios, the high sequence similarity can be explained by purifying selection that maintains sequence identity within tRNA families. Yet, there is also a possibility that the sequence similarity is not due to purifying selection but the result of ‘concerted evolution’, an evolutionary process that maintains sequence identity by frequent recombination events among copies of the same gene family (<xref ref-type="bibr" rid="bib41">Munz et al., 1982</xref>; <xref ref-type="bibr" rid="bib2">Amstutz et al., 1985</xref>; <xref ref-type="bibr" rid="bib62">Teshima and Innan, 2004</xref>). This possibility implies that the high conservation observed within tRNA gene families is not due to functionality, but is rather the result of neutral evolution. At present, it is not possible to determine which of the two possibilities explains best the observed sequence identity.</p><p>If a single point mutation in one of the tRNA<sup>Arg</sup><sub>UCU</sub> copies enables it to function like a tRNA from a different family, what were the evolutionary constraints that have left some families with more members while others with fewer? Is purifying selection acting to deliberately maintain low levels of certain tRNAs? Such selection would render their corresponding codons ‘non-optimal’. To examine potential adaptive functions of tRNA family sizes, we tested the consequences of increasing the sizes of several tRNA families. We found that keeping low copy tRNA families is adaptive, as increasing their copy number can result in a proteotoxic stress due to problems in protein folding.</p><p>Most of the published work on the functionality of codons that correspond to rare tRNAs have so far tested how modified codon usage of specific proteins influences their proper folding (<xref ref-type="bibr" rid="bib16">Crombie et al., 1992</xref>; <xref ref-type="bibr" rid="bib31">Komar et al., 1999</xref>; <xref ref-type="bibr" rid="bib15">Cortazzo et al., 2002</xref>; <xref ref-type="bibr" rid="bib64">Tsai et al., 2008</xref>; <xref ref-type="bibr" rid="bib73">Zhang et al., 2009</xref>; <xref ref-type="bibr" rid="bib75">Zhou et al., 2013</xref>). In contrast, we took a different approach, in which no protein coding gene sequence is modified, but rather the tRNA supply is manipulated. Thus, the effect we generated could be exerted on all genes, and we could indeed detect it as a global proteotoxic stress in the cell. Our observations are consistent with the theory that programmed pauses in the translation process could promote proper folding during translation (<xref ref-type="bibr" rid="bib63">Thanaraj, 1996</xref>; <xref ref-type="bibr" rid="bib33">Kramer et al., 2009</xref>; <xref ref-type="bibr" rid="bib11">Cabrita et al., 2010</xref>; <xref ref-type="bibr" rid="bib70">Wilke and Drummond, 2010</xref>). The overexpression of the rare tRNA could have thus impaired with cotranslation folding. Yet, there could be additional reasons for the observed proteotoxic stress, which are not necessarily mutually exclusive. First, overexpression of tRNA<sup>Arg</sup><sub>CCU</sub> may result in misincorporation of arginine into non-arginine codons. Second, other aminoacyl tRNA synthetases may aminoacylate an incorrect amino acid to the highly expressed tRNA. Misloading will result in the incorporation of a different amino acid instead of arginine. A potential part of the observed proteotoxic stress due to misincorporation still remains to be studied. Yet, such an effect should be relevant not only for the overexpression of the rare tRNA<sup>Arg</sup><sub>CCU</sub> but also for the overexpression of the abundant tRNA<sup>Arg</sup><sub>UCU</sub>. Our results show a sever proteotoxic stress only upon expression of the rare tRNA, thus landing more support to the intriguing hypothesis, that the proteotoxic phenotypes observed are due to converting a slow translating codon, scattered in many genes in the genome, into a fast one. This notion is consistent with and complementary to the picture that emerges from single gene-based analyses.</p><p>When facing the need to adapt, the tRNA pool (that is the supply) provides evolutionary plasticity to the translation machinery. The ability of the tRNA pool to change rapidly can be mainly attributed to its unique architecture in the form of multimember gene families. Only on a much longer evolutionary timescale, will the genome-wide codon usage of genes change so as to further fine tune the translational balance. Notably, the plasticity of the tRNA genes is constrained by the need to maintain proper protein folding (<xref ref-type="bibr" rid="bib20">Drummond and Wilke, 2008</xref>). Thus, the need to accommodate changes in codon usage demands acts together with protein folding constrains to shape the tRNA pool in the living cells.</p></sec><sec id="s4" sec-type="materials|methods"><title>Materials and methods</title><sec id="s4-1"><title>Yeast strains and plasmids</title><p>The following <italic>S. cerevisiae</italic> strains were used in this study: ΔtRNA<sup>Arg</sup><sub>CCU</sub> (based on Y5565, genetic background: <italic>ΔtR(CCU)J::Hyg, MATα, can1Δ::MFA1pr-HIS3 mfα1Δ::MFα1pr-LEU2 lyp1Δ ura3Δ0 leu2Δ0)</italic> (<xref ref-type="bibr" rid="bib7">Bloom-Ackermann et al., In press</xref>) was used for lab-evolution experiments. <italic>MutΔtRNA</italic><sup><italic>Arg</italic></sup><sub><italic>CCU</italic></sub> is based on ΔtRNA<sup>Arg</sup><sub>CCU</sub> and carries a mutation (T→C transition) in <italic>tR(UCU)K</italic> gene plus a <italic>URA3</italic> selection marker. BY384 (<italic>MATa leu2Δ1 lys2Δ202 trp1Δ63 ura3-52 his3Δ200</italic>) was used to generate a complete deletion of the tRNA<sup>Ser</sup><sub>GCU</sub> gene family. BY4741 (<italic>MATa his3Δ1 leu2Δ0 met15Δ0 ura3Δ0</italic>) and BY4742 (<italic>MATα his3Δ1 leu2Δ0 lys2Δ0 ura3Δ0</italic>) were used for examining the effect of increasing tRNA gene copy number.</p><p>Plasmids used in this study to express tRNA genes were pRS316 (CEN, <italic>URA3</italic>), pRS425 (2µ, <italic>LEU2</italic>), and pRS426 (2µ, <italic>URA3</italic>). For the rescue assays of the quadruple serine deletion, the pQF50 (2µ, <italic>URA3</italic>) and pQF150 (2µ, <italic>LEU2</italic>) were used. For the protein quality control assays the pGAL-VHL-mCherry (2µ, <italic>LEU2)</italic> plasmid was used (<xref ref-type="bibr" rid="bib27">Kaganovich et al., 2008</xref>). For additional information on plasmids and primers see <xref ref-type="supplementary-material" rid="SD3-data">Supplementary file 2</xref>.</p></sec><sec id="s4-2"><title>Media</title><p>Cultures were grown at 30°C in either rich medium (1% bacto-yeast extract, 2% bacto-peptone and 2% dextrose [YPD]) or synthetic medium (0.67% yeast nitrogen base with ammonium sulfate and without amino acids and 2% dextrose, containing the appropriate supplements for plasmid selection). Protein quality control assays were performed on synthetic medium supplemented with 2% galactose as a carbon source. All chemicals used to create the media were manufactured by BD. All sugars, nucleic acids and amino acids were manufactured by Sigma-Aldrich.</p></sec><sec id="s4-3"><title>Evolution experiments</title><p>Lab-evolution experiments were carried out by serial dilution. Cells were grown on 1.2 ml of YPD at 30°C until reaching stationary phase and then diluted by a factor of 1:120 into fresh media (6.9 generations per dilution). This procedure was repeated daily until population growth under the applied condition matched the wild type. In all measurements of evolved populations, we used a population sample and not selected clones.</p></sec><sec id="s4-4"><title>Liquid growth measurements</title><p>The cultures were grown at the relevant condition, and optical density (OD)<sub>600</sub> measurements were taken during the growth at 30–45 min intervals until reaching early stationary phase. Qualitative growth comparisons were performed using 96-well plates (Thermo Scientific) in which 2 strains were divided on the plate in a checkerboard manner on the plate to cancel out positional geographical effects. For each strain, a growth curve was obtained by averaging over 48 wells.</p></sec><sec id="s4-5"><title>Growth on 5-fluoro-orotic acid (5-FOA) plates</title><p>Strains were grown for 2 days in a non-selective liquid medium, which contains uracil (YPD), to allow growth of cells that lost the plasmid containing the <italic>URA3</italic> counterselectable marker (<xref ref-type="bibr" rid="bib8">Boeke et al., 1984</xref>). Then, 100 µl were plated on a YPD plate and replicated on the following day on either YPD or standard 5-fluoro-orotic acid (5-FOA) (US Biological) plates to identify potential colonies that lost the plasmid. Following 2 days of incubation at 30°C, growth of the colonies was scored.</p></sec><sec id="s4-6"><title>Measurements for saturation of the protein quality control machinery</title><p>We used a previously published method that allows examination of the protein quality control of the cell (<xref ref-type="bibr" rid="bib27">Kaganovich et al., 2008</xref>). This assay provides an indication for the protein unfolding stress in cells by assessing the load on the protein quality control machinery. In this assay, the cells were introduced with a high copy number plasmid that contains the human gene VHL fused to a fluorescent tag (mCherry). VHL is a naturally unstructured protein and is dependent on two additional proteins (Elongin B and C) for proper folding in human cells. Expressing VHL in yeast cells, which lack VHL’s complex partners, leads to misfolding of the translated proteins. Under normal conditions, the misfolded VHL proteins are handled by the cell’s quality control machinery. When the quality control machinery is not saturated, the fluorescently tagged VHL appears in the cytosol. However, under stress, in which the quality control machinery is fully occupied, misfolded proteins in the cytosol are processed into dedicated inclusions (JUNQ and IPOD) and form punctum structures. Hence, a punctum phenotype of the VHL-mCherry construct is an indication for cells that experience proteotoxic stress and saturation of the protein quality control machinery.</p><p>Wild-type yeast cells harboring the pGAL-VHL-mCherry (CHFP) fusion plasmid (<xref ref-type="bibr" rid="bib27">Kaganovich et al., 2008</xref>) and either an additional empty plasmid or the tRNA overexpression plasmid, were grown overnight on SC+2% raffinose, diluted into SC+2% galactose and grown at 30°C for 6 hr. The cells were visualized using an Olympus IX71 microscope (Olympus) controlled by Delta Vision SoftWoRx 3.5.1 software, with X60 oil lens. Images were captured by a Photometrics Coolsnap HQ camera with excitation at 555/28 nm and emission at 617/73 nm (mCherry). The images were scored using ImageJ image processing and analysis software. The percentage of cells harboring VHL-CHFP foci (Puncta) in the overexpression strains were normalized to the level in a control strain carrying an empty plasmid.</p></sec><sec id="s4-7"><title>Computational identification of anticodon switching events</title><p>To characterize the extent of anticodon switching across the tree of life, we first downloaded all tRNA sequences from the Genomic tRNA Database (<xref ref-type="bibr" rid="bib12">Chan and Lowe, 2009</xref>). The sequences in this database were discovered with the tRNAscan algorithm (<xref ref-type="bibr" rid="bib36">Lowe and Eddy, 1997</xref>), which finds tRNA sequences by scanning genomic DNA. We removed psuedogene tRNAs, which are defined as those tRNAs with a COVE score less than 40.0 (<xref ref-type="bibr" rid="bib36">Lowe and Eddy, 1997</xref>) . For each remaining tRNA sequence, we masked its anticodon triplet as ‘NNN’. We next grouped all tRNA sequences by their species, and then aligned the sequences for each species using Muscle with default settings (<xref ref-type="bibr" rid="bib21">Edgar, 2004</xref>). For each species, we inferred a maximum likelihood phylogeny of its tRNA sequences using RAxML with the GTRCAT model (<xref ref-type="bibr" rid="bib59">Stamatakis, 2006</xref>). We calculated statistical support for tree branches using SH-like approximate likelihood ratio test (<xref ref-type="bibr" rid="bib3">Anisimova and Gascuel, 2006</xref>). We next interrogated each species’ tRNA phylogeny using DendroPy (<xref ref-type="bibr" rid="bib61">Sukumaran and Holder, 2010</xref>). Specifically, we identified those tRNA sequences harboring an anticodon that appears in the genome more than once; for each of these tRNA sequences, we found the shortest distance to another tRNA with the same anticodon (d<sub>same</sub>) and the shortest distance to another tRNA with a different anticodon (d<sub>diff</sub>). We labeled tRNAs as putatively ‘switched’ if d<sub>diff</sub><d<sub>same</sub>.</p></sec><sec id="s4-8"><title>RT-qPCR measurements of HSP genes in strains overexpressing tRNAs</title><p>Cultures were grown in rich medium at 30°C to a cell concentration of 1 × 10<sup>7</sup> cells/ml. Then, RNA was extracted using MasterPure kit (Epicentre-illumina) (EPICENTER Biotechnologies), and used as a template for quantitative RT-PCR using light cycler 480 SYBR I master kit (Roche Applied Science) and the LightCycler 480 system (Roche Applied Science), according to the manufacturer’s instructions.</p></sec><sec id="s4-9"><title>Genomic copies of tRNA<sup>Arg</sup><sub>UCU</sub> mutated during lab-evolution experiments</title><p>In all, 4 independent lab-evolution experiments that started with ΔtRNA<sup>Arg</sup><sub>CCU</sub> as the ancestral strain showed full recovery of the deletion phenotype after 200 generations. In each of the evolved populations a mutation in one of the copies of tRNA<sup>Arg</sup><sub>UCU</sub> was found to change the anticodon from UCU to CCU. The genomic copies of tRNA<sup>Arg</sup><sub>UCU</sub> that were found to carry the mutation were: <italic>tR(UCU)K</italic>, <italic>tR(UCU)G1</italic> and <italic>tR(UCU)D</italic> that was changed in two of the independent cultures.</p></sec><sec id="s4-10"><title>The contribution of different anticodon positions to tRNA switching events</title><p>Out of 4245 anticodon switching events that we detected, the first position in the anticodon was changed in 2540 cases while the second and third were only involved in 1448 and 1330 cases, respectively.</p></sec></sec></body><back><ack id="ack"><title>Acknowledgements</title><p>We thank M Schuldiner, D Kaganovich, and S Leidel for kindly providing plasmids. We also thank the Pilpel lab, and especially H Gingold, for fruitful discussions. We acknowledge T Ast, Y Cohen and R Ackermann for critical reading of the manuscript. We thank the European Research Council (ERC) (YP), the Ben-May Charitable Trust (YP) and the NIH (grant 1P50GM107632) (JDB) for grant support. YP is an incumbent of the Ben May Professorial Chair.</p></ack><sec sec-type="additional-information"><title>Additional information</title><fn-group content-type="competing-interest"><title>Competing interests</title><fn fn-type="conflict" id="conf1"><p>The authors declare that no competing interests exist.</p></fn></fn-group><fn-group content-type="author-contribution"><title>Author contributions</title><fn fn-type="con" id="con1"><p>AHY, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con2"><p>IF, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con3"><p>ZB-A, Design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con4"><p>VH-S, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con5"><p>YC-A, Acquisition of data</p></fn><fn fn-type="con" id="con6"><p>QF, Acquisition of data, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con7"><p>JDB, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con8"><p>OD, Design, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con9"><p>YP, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn></fn-group></sec><sec sec-type="supplementary-material"><title>Additional files</title><supplementary-material id="SD2-data"><object-id pub-id-type="doi">10.7554/eLife.01339.014</object-id><label>Supplementary file 1.</label><caption><title>Usage of arginine codons in <italic>Saccharomyces cerevisiae</italic>.</title><p>(<bold>A</bold>) The AGG codon constitutes approximately 21% of the arginine codons in the yeast genome and approximately 16% of the arginine codons in the yeast transcriptome under standard lab conditions. In comparison, the AGA codon constitutes approximately 47.5% of the arginine codons in the genome and approximately 56% of the arginine codons in the transcriptome (<xref ref-type="bibr" rid="bib22">Gingold et al., 2012</xref>). (<bold>B</bold>) Statistical significance of differences between proportions of discrepancies when masking the anticodon triplet versus control triplets. (a) Six control triplets were masked (<xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1</xref>), and the analysis to identify switched tRNAs was repeated (see ‘Materials and methods’). For each masked triplet, we computed the percentage of tRNA sequences that are not clustered according to their triplet content in each species. (b) The proportion of switched tRNAs was averaged across all species, and (c) the standard error of this average was computed. (d) The distribution of tRNA switching proportions from the original anticodon analysis was compared to the distribution of discrepancies from the control triplet, using a paired <italic>t</italic> test. The p value from this comparison is reported here. (<bold>C</bold>) Various yeast species tend to keep tRNAArgCCU in a single copy. All examined yeast species maintain a single copy of tRNA<sup>Arg</sup><sub>CCU</sub> compared to tRNA<sup>Arg</sup><sub>UCU</sub>, which is mostly found in multiple copies. The copy number of tRNA<sup>Arg</sup><sub>UCU</sub> and tRNA<sup>Arg</sup><sub>CCU</sub> genes is shown together with the codon usage of AGG and AGA codons in the different yeast genomes.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.014">http://dx.doi.org/10.7554/eLife.01339.014</ext-link></p></caption><media mime-subtype="xlsx" mimetype="application" xlink:href="elife01339s002.xlsx"/></supplementary-material><supplementary-material id="SD3-data"><object-id pub-id-type="doi">10.7554/eLife.01339.015</object-id><label>Supplementary file 2.</label><caption><title>Plasmids and primers.</title><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.01339.015">http://dx.doi.org/10.7554/eLife.01339.015</ext-link></p></caption><media mime-subtype="xlsx" mimetype="application" xlink:href="elife01339s003.xlsx"/></supplementary-material><sec sec-type="datasets"><title>Major dataset</title><p>The following previously published datasets were used:</p><p><related-object content-type="generated-dataset" document-id="Dataset ID and/or url" document-id-type="dataset" document-type="data" id="dataro1"><name><surname>Chan</surname><given-names>PP</given-names></name>, <name><surname>Lowe</surname><given-names>TM</given-names></name>, <year>2009</year><x>, </x><source>Genomic tRNA database</source><x>, </x><ext-link ext-link-type="uri" xlink:href="http://gtrnadb.ucsc.edu/">http://gtrnadb.ucsc.edu/</ext-link><x>, </x><comment>We downloaded all the known tRNA sequences that are publicly available from the Genomic tRNA Database.</comment></related-object></p></sec></sec><ref-list><title>References</title><ref id="bib1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Akashi</surname><given-names>H</given-names></name></person-group><year>1994</year><article-title>Synonymous codon usage in <italic>Drosophila melanogaster</italic>: natural selection and translational accuracy</article-title><source>Genetics</source><volume>136</volume><fpage>927</fpage><lpage>935</lpage></element-citation></ref><ref id="bib2"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Amstutz</surname><given-names>H</given-names></name><name><surname>Munz</surname><given-names>P</given-names></name><name><surname>Heyer</surname><given-names>WD</given-names></name><name><surname>Leupoid</surname><given-names>U</given-names></name><name><surname>Kohli</surname><given-names>J</given-names></name></person-group><year>1985</year><article-title>Concerted evolution of tRNA genes: intergenic conversion among three unlinked serine tRNA genes in <italic>S. pombe</italic></article-title><source>Cell</source><volume>40</volume><fpage>879</fpage><lpage>886</lpage><pub-id pub-id-type="doi">10.1016/0092-8674(85)90347-2</pub-id></element-citation></ref><ref id="bib3"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Anisimova</surname><given-names>M</given-names></name><name><surname>Gascuel</surname><given-names>O</given-names></name></person-group><year>2006</year><article-title>Approximate likelihood-ratio test for branches: a fast, accurate, and powerful alternative</article-title><source>Systematic Biology</source><volume>55</volume><fpage>539</fpage><lpage>552</lpage><pub-id pub-id-type="doi">10.1080/10635150600755453</pub-id></element-citation></ref><ref id="bib4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Aström</surname><given-names>SU</given-names></name><name><surname>von Pawel-Rammingen</surname><given-names>U</given-names></name><name><surname>Byström</surname><given-names>AS</given-names></name></person-group><year>1993</year><article-title>The yeast initiator tRNAMet can act as an elongator tRNA(Met) in vivo</article-title><source>Journal of Molecular Biology</source><volume>233</volume><fpage>43</fpage><lpage>58</lpage><pub-id pub-id-type="doi">10.1006/jmbi.1993.1483</pub-id></element-citation></ref><ref id="bib5"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Begley</surname><given-names>U</given-names></name><name><surname>Dyavaiah</surname><given-names>M</given-names></name><name><surname>Patil</surname><given-names>A</given-names></name><name><surname>Rooney</surname><given-names>JP</given-names></name><name><surname>DiRenzo</surname><given-names>D</given-names></name><name><surname>Young</surname><given-names>CM</given-names></name><name><surname>Conklin</surname><given-names>DS</given-names></name><name><surname>Zitomer</surname><given-names>RS</given-names></name><name><surname>Begley</surname><given-names>TJ</given-names></name></person-group><year>2007</year><article-title>Trm9-catalyzed tRNA modifications link translation to the DNA damage response</article-title><source>Molecular Cell</source><volume>28</volume><fpage>860</fpage><lpage>870</lpage><pub-id pub-id-type="doi">10.1016/j.molcel.2007.09.021</pub-id></element-citation></ref><ref id="bib6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bermudez-Santana</surname><given-names>C</given-names></name><name><surname>Attolini</surname><given-names>CS</given-names></name><name><surname>Kirsten</surname><given-names>T</given-names></name><name><surname>Engelhardt</surname><given-names>J</given-names></name><name><surname>Prohaska</surname><given-names>SJ</given-names></name><name><surname>Steigele</surname><given-names>S</given-names></name><name><surname>Stadler</surname><given-names>PF</given-names></name></person-group><year>2010</year><article-title>Genomic organization of eukaryotic tRNAs</article-title><source>BMC Genomics</source><volume>11</volume><fpage>270</fpage><pub-id pub-id-type="doi">10.1186/1471-2164-11-270</pub-id></element-citation></ref><ref id="bib7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bloom-Ackermann</surname><given-names>Z</given-names></name><name><surname>Navon</surname><given-names>S</given-names></name><name><surname>Gingold</surname><given-names>H</given-names></name><name><surname>Towers</surname><given-names>R</given-names></name><name><surname>Pilpel</surname><given-names>Y</given-names></name><name><surname>Dahan</surname><given-names>O</given-names></name></person-group><year>2014</year><article-title>A Comprehensive tRNA deletion library unravels the genetic architecture of the tRNA pool</article-title><source>Public Library of Science Genetics</source><volume>10</volume><fpage>e1004084</fpage><pub-id pub-id-type="doi">10.1371/journal.pgen.1004084</pub-id></element-citation></ref><ref id="bib8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Boeke</surname><given-names>JD</given-names></name><name><surname>LaCroute</surname><given-names>F</given-names></name><name><surname>Fink</surname><given-names>GR</given-names></name></person-group><year>1984</year><article-title>A positive selection for mutants lacking orotidine-5’-phosphate decarboxylase activity in yeast: 5-fluoro-orotic acid resistance</article-title><source>Molecular and General Genetics</source><volume>197</volume><fpage>345</fpage><lpage>346</lpage><pub-id pub-id-type="doi">10.1007/BF00330984</pub-id></element-citation></ref><ref id="bib9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bulmer</surname><given-names>M</given-names></name></person-group><year>1991</year><article-title>The selection-mutation-drift theory of synonymous codon usage</article-title><source>Genetics</source><volume>129</volume><fpage>897</fpage><lpage>907</lpage></element-citation></ref><ref id="bib10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Byström</surname><given-names>AS</given-names></name><name><surname>Fink</surname><given-names>GR</given-names></name></person-group><year>1989</year><article-title>A functional analysis of the repeated methionine initiator tRNA genes (IMT) in yeast</article-title><source>Molecular and General Genetics</source><volume>216</volume><fpage>276</fpage><lpage>286</lpage><pub-id pub-id-type="doi">10.1007/BF00334366</pub-id></element-citation></ref><ref id="bib11"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cabrita</surname><given-names>LD</given-names></name><name><surname>Dobson</surname><given-names>CM</given-names></name><name><surname>Christodoulou</surname><given-names>J</given-names></name></person-group><year>2010</year><article-title>Protein folding on the ribosome</article-title><source>Current Opinion In Structural Biology</source><volume>20</volume><fpage>33</fpage><lpage>45</lpage><pub-id pub-id-type="doi">10.1016/j.sbi.2010.01.005</pub-id></element-citation></ref><ref id="bib12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chan</surname><given-names>PP</given-names></name><name><surname>Lowe</surname><given-names>TM</given-names></name></person-group><year>2009</year><article-title>GtRNAdb: a database of transfer RNA genes detected in genomic sequence</article-title><source>Nucleic Acids Research</source><volume>37</volume><fpage>D93</fpage><lpage>D97</lpage><pub-id pub-id-type="doi">10.1093/nar/gkn787</pub-id></element-citation></ref><ref id="bib13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Clare</surname><given-names>JJ</given-names></name><name><surname>Belcourt</surname><given-names>M</given-names></name><name><surname>Farabaugh</surname><given-names>PJ</given-names></name></person-group><year>1988</year><article-title>Efficient translational frameshifting occurs within a conserved sequence of the overlap between the two genes of a yeast Ty1 transposon</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>85</volume><fpage>6816</fpage><lpage>6820</lpage><pub-id pub-id-type="doi">10.1073/pnas.85.18.6816</pub-id></element-citation></ref><ref id="bib14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cochella</surname><given-names>L</given-names></name><name><surname>Green</surname><given-names>R</given-names></name></person-group><year>2005</year><article-title>An active role for tRNA in decoding beyond codon:anticodon pairing</article-title><source>Science</source><volume>308</volume><fpage>1178</fpage><lpage>1180</lpage><pub-id pub-id-type="doi">10.1126/science.1111408</pub-id></element-citation></ref><ref id="bib15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cortazzo</surname><given-names>P</given-names></name><name><surname>Cerveñansky</surname><given-names>C</given-names></name><name><surname>Marín</surname><given-names>M</given-names></name><name><surname>Reiss</surname><given-names>C</given-names></name><name><surname>Ehrlich</surname><given-names>R</given-names></name><name><surname>Deana</surname><given-names>A</given-names></name></person-group><year>2002</year><article-title>Silent mutations affect in vivo protein folding in <italic>Escherichia coli</italic></article-title><source>Biochemical and Biophysical Research</source><volume>293</volume><fpage>537</fpage><lpage>541</lpage><pub-id pub-id-type="doi">10.1016/S0006-291X(02)00226-7</pub-id></element-citation></ref><ref id="bib16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Crombie</surname><given-names>T</given-names></name><name><surname>Swaffield</surname><given-names>JC</given-names></name><name><surname>Brown</surname><given-names>AJ</given-names></name></person-group><year>1992</year><article-title>Protein folding within the cell is influenced by controlled rates of polypeptide elongation</article-title><source>Journal of Molecular Biology</source><volume>228</volume><fpage>7</fpage><lpage>12</lpage><pub-id pub-id-type="doi">10.1016/0022-2836(92)90486-4</pub-id></element-citation></ref><ref id="bib17"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Doherty</surname><given-names>A</given-names></name><name><surname>McInerney</surname><given-names>JO</given-names></name></person-group><year>2013</year><article-title>Translational selection frequently Overcomes genetic drift in shaping synonymous codon usage patterns in vertebrates</article-title><source>Molecular Biology and Evolution</source><volume>30</volume><fpage>2263</fpage><lpage>2267</lpage><pub-id pub-id-type="doi">10.1093/molbev/mst128</pub-id></element-citation></ref><ref id="bib18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dong</surname><given-names>H</given-names></name><name><surname>Nilsson</surname><given-names>L</given-names></name><name><surname>Kurland</surname><given-names>CG</given-names></name></person-group><year>1996</year><article-title>Co-variation of tRNA abundance and codon usage in <italic>Escherichia coli</italic> at different growth rates</article-title><source>Journal of Molecular Biology</source><volume>260</volume><fpage>649</fpage><lpage>663</lpage><pub-id pub-id-type="doi">10.1006/jmbi.1996.0428</pub-id></element-citation></ref><ref id="bib19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Dos Reis</surname><given-names>M</given-names></name><name><surname>Savva</surname><given-names>R</given-names></name><name><surname>Wernisch</surname><given-names>L</given-names></name></person-group><year>2004</year><article-title>Solving the riddle of codon usage preferences: a test for translational selection</article-title><source>Nucleic Acids Research</source><volume>32</volume><fpage>5036</fpage><lpage>5044</lpage><pub-id pub-id-type="doi">10.1093/nar/gkh834</pub-id></element-citation></ref><ref id="bib20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Drummond</surname><given-names>DA</given-names></name><name><surname>Wilke</surname><given-names>CO</given-names></name></person-group><year>2008</year><article-title>Mistranslation-induced protein misfolding as a dominant constraint on coding-sequence evolution</article-title><source>Cell</source><volume>134</volume><fpage>341</fpage><lpage>352</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2008.05.042</pub-id></element-citation></ref><ref id="bib21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Edgar</surname><given-names>RC</given-names></name></person-group><year>2004</year><article-title>MUSCLE: multiple sequence alignment with high accuracy and high throughput</article-title><source>Nucleic Acids Research</source><volume>32</volume><fpage>1792</fpage><lpage>1797</lpage><pub-id pub-id-type="doi">10.1093/nar/gkh340</pub-id></element-citation></ref><ref id="bib22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gingold</surname><given-names>H</given-names></name><name><surname>Dahan</surname><given-names>O</given-names></name><name><surname>Pilpel</surname><given-names>Y</given-names></name></person-group><year>2012</year><article-title>Dynamic changes in translational efficiency are deduced from codon usage of the transcriptome</article-title><source>Nucleic Acids Research</source><volume>40</volume><fpage>1</fpage><lpage>11</lpage><pub-id pub-id-type="doi">10.1093/nar/gks772</pub-id></element-citation></ref><ref id="bib23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Goodman</surname><given-names>DB</given-names></name><name><surname>Church</surname><given-names>GM</given-names></name><name><surname>Kosuri</surname><given-names>S</given-names></name></person-group><year>2013</year><article-title>Causes and effects of N-Terminal codon bias in bacterial genes</article-title><source>Science</source><volume>342</volume><fpage>475</fpage><lpage>479</lpage><pub-id pub-id-type="doi">10.1126/science.1241934</pub-id></element-citation></ref><ref id="bib24"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Higgs</surname><given-names>PG</given-names></name><name><surname>Ran</surname><given-names>W</given-names></name></person-group><year>2008</year><article-title>Coevolution of codon usage and tRNA genes leads to alternative stable states of biased codon usage</article-title><source>Molecular Biology and Evolution</source><volume>25</volume><fpage>2279</fpage><lpage>2291</lpage><pub-id pub-id-type="doi">10.1093/molbev/msn173</pub-id></element-citation></ref><ref id="bib25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hudson</surname><given-names>NJ</given-names></name><name><surname>Gu</surname><given-names>Q</given-names></name><name><surname>Nagaraj</surname><given-names>SH</given-names></name><name><surname>Ding</surname><given-names>YS</given-names></name><name><surname>Dalrymple</surname><given-names>BP</given-names></name><name><surname>Reverter</surname><given-names>A</given-names></name></person-group><year>2011</year><article-title>Eukaryotic evolutionary transitions are associated with extreme codon bias in functionally-related proteins</article-title><source>Public Library of Science One</source><volume>6</volume><fpage>e25457</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0025457</pub-id></element-citation></ref><ref id="bib26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Itoh</surname><given-names>T</given-names></name><name><surname>Takemoto</surname><given-names>K</given-names></name><name><surname>Mori</surname><given-names>H</given-names></name><name><surname>Gojobori</surname><given-names>T</given-names></name></person-group><year>1999</year><article-title>Evolutionary instability of operon structures disclosed by sequence comparisons of complete microbial genomes</article-title><source>Molecular Biology and Evolution</source><volume>16</volume><fpage>332</fpage><lpage>346</lpage><pub-id pub-id-type="doi">10.1093/oxfordjournals.molbev.a026114</pub-id></element-citation></ref><ref id="bib27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kaganovich</surname><given-names>D</given-names></name><name><surname>Kopito</surname><given-names>R</given-names></name><name><surname>Frydman</surname><given-names>J</given-names></name></person-group><year>2008</year><article-title>Misfolded proteins partition between two distinct quality control compartments</article-title><source>Nature</source><volume>454</volume><fpage>1088</fpage><lpage>1095</lpage><pub-id pub-id-type="doi">10.1038/nature07195</pub-id></element-citation></ref><ref id="bib28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kanaya</surname><given-names>S</given-names></name><name><surname>Yamada</surname><given-names>Y</given-names></name><name><surname>Kudo</surname><given-names>Y</given-names></name><name><surname>Ikemura</surname><given-names>T</given-names></name></person-group><year>1999</year><article-title>Studies of codon usage and tRNA genes of 18 unicellular organisms and quantification of <italic>Bacillus subtilis</italic> tRNAs: gene expression level and species-specific diversity of codon usage based on multivariate analysis</article-title><source>Gene</source><volume>238</volume><fpage>143</fpage><lpage>155</lpage><pub-id pub-id-type="doi">10.1016/S0378-1119(99)00225-5</pub-id></element-citation></ref><ref id="bib29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kawakami</surname><given-names>K</given-names></name><name><surname>Pande</surname><given-names>S</given-names></name><name><surname>Faiola</surname><given-names>B</given-names></name><name><surname>Moore</surname><given-names>DP</given-names></name><name><surname>Boeke</surname><given-names>JD</given-names></name><name><surname>Farabaugh</surname><given-names>PJ</given-names></name><name><surname>Strathern</surname><given-names>JN</given-names></name><name><surname>Nakamura</surname><given-names>Y</given-names></name><name><surname>Garfinkel</surname><given-names>DJ</given-names></name></person-group><year>1993</year><article-title>A rare tRNA-Arg(CCU) that regulates Ty1 element ribosomal frameshifting is essential for Ty1 retrotransposition in <italic>Saccharomyces cerevisiae</italic></article-title><source>Genetics</source><volume>135</volume><fpage>309</fpage><lpage>320</lpage></element-citation></ref><ref id="bib30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kolitz</surname><given-names>SE</given-names></name><name><surname>Lorsch</surname><given-names>JR</given-names></name></person-group><year>2010</year><article-title>Eukaryotic initiator tRNA: finely tuned and ready for action</article-title><source>FEBS Letters</source><volume>584</volume><fpage>396</fpage><lpage>404</lpage><pub-id pub-id-type="doi">10.1016/j.febslet.2009.11.047</pub-id></element-citation></ref><ref id="bib31"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Komar</surname><given-names>AA</given-names></name><name><surname>Lesnik</surname><given-names>T</given-names></name><name><surname>Reiss</surname><given-names>C</given-names></name></person-group><year>1999</year><article-title>Synonymous codon substitutions affect ribosome traffic and protein folding during in vitro translation</article-title><source>FEBS Letters</source><volume>462</volume><fpage>387</fpage><lpage>391</lpage><pub-id pub-id-type="doi">10.1016/S0014-5793(99)01566-5</pub-id></element-citation></ref><ref id="bib32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Konevega</surname><given-names>AL</given-names></name><name><surname>Soboleva</surname><given-names>NG</given-names></name><name><surname>Makhno</surname><given-names>VI</given-names></name><name><surname>Semenkov</surname><given-names>YP</given-names></name><name><surname>Wintermeyer</surname><given-names>W</given-names></name><name><surname>Rodnina</surname><given-names>MV</given-names></name><name><surname>Katunin</surname><given-names>VI</given-names></name></person-group><year>2004</year><article-title>Purine bases at position 37 of tRNA stabilize codon-anticodon interaction in the ribosomal A site by stacking and Mg2+-dependent interactions</article-title><source>RNA</source><volume>10</volume><fpage>90</fpage><lpage>101</lpage><pub-id pub-id-type="doi">10.1261/rna.5142404</pub-id></element-citation></ref><ref id="bib33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kramer</surname><given-names>G</given-names></name><name><surname>Boehringer</surname><given-names>D</given-names></name><name><surname>Ban</surname><given-names>N</given-names></name><name><surname>Bukau</surname><given-names>B</given-names></name></person-group><year>2009</year><article-title>The ribosome as a platform for co-translational processing, folding and targeting of newly synthesized proteins</article-title><source>Nature Structural and Molecular Biology</source><volume>16</volume><fpage>589</fpage><lpage>597</lpage><pub-id pub-id-type="doi">10.1038/nsmb.1614</pub-id></element-citation></ref><ref id="bib34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kudla</surname><given-names>G</given-names></name><name><surname>Murray</surname><given-names>A</given-names></name><name><surname>Tollervey</surname><given-names>D</given-names></name><name><surname>Plotkin</surname><given-names>J</given-names></name></person-group><year>2009</year><article-title>Coding-sequence determinants of gene expression in <italic>Escherichia coli</italic></article-title><source>Science</source><volume>324</volume><fpage>255</fpage><lpage>258</lpage><pub-id pub-id-type="doi">10.1126/science.1170160</pub-id></element-citation></ref><ref id="bib35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lenski</surname><given-names>RE</given-names></name><name><surname>Rose</surname><given-names>MR</given-names></name><name><surname>Simpson</surname><given-names>SC</given-names></name><name><surname>Tadler</surname><given-names>SC</given-names></name></person-group><year>1991</year><article-title>Long-term experimental evolution in <italic>Escherichia coli</italic>. I. Adaptation and divergence during 2,000 generations</article-title><source>The American Naturalist</source><volume>138</volume><fpage>1315</fpage><pub-id pub-id-type="doi">10.1086/285289</pub-id></element-citation></ref><ref id="bib36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lowe</surname><given-names>TM</given-names></name><name><surname>Eddy</surname><given-names>SR</given-names></name></person-group><year>1997</year><article-title>tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence</article-title><source>Nucleic Acids Research</source><volume>25</volume><fpage>955</fpage><lpage>964</lpage></element-citation></ref><ref id="bib37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Man</surname><given-names>O</given-names></name><name><surname>Pilpel</surname><given-names>Y</given-names></name></person-group><year>2007</year><article-title>Differential translation efficiency of orthologous genes is involved in phenotypic divergence of yeast species</article-title><source>Nature Genetics</source><volume>39</volume><fpage>415</fpage><lpage>421</lpage><pub-id pub-id-type="doi">10.1038/ng1967</pub-id></element-citation></ref><ref id="bib38"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McClellan</surname><given-names>AJ</given-names></name><name><surname>Scott</surname><given-names>MD</given-names></name><name><surname>Frydman</surname><given-names>J</given-names></name></person-group><year>2005</year><article-title>Folding and quality control of the VHL tumor suppressor proceed through distinct chaperone pathways</article-title><source>Cell</source><volume>121</volume><fpage>739</fpage><lpage>748</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2005.03.024</pub-id></element-citation></ref><ref id="bib39"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>McFarlane</surname><given-names>RJ</given-names></name><name><surname>Whitehall</surname><given-names>SK</given-names></name></person-group><year>2009</year><article-title>tRNA genes in eukaryotic genome organization and reorganization</article-title><source>Cell Cycle</source><volume>8</volume><fpage>3102</fpage><lpage>3106</lpage><pub-id pub-id-type="doi">10.4161/cc.8.19.9625</pub-id></element-citation></ref><ref id="bib40"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Müller</surname><given-names>EC</given-names></name><name><surname>Wittmann-Liebold</surname><given-names>B</given-names></name></person-group><year>1997</year><article-title>Phylogenetic relationship of organisms obtained by ribosomal protein comparison</article-title><source>Cellular and Molecular Life Sciences</source><volume>53</volume><fpage>34</fpage><lpage>50</lpage><pub-id pub-id-type="doi">10.1007/PL00000578</pub-id></element-citation></ref><ref id="bib41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Munz</surname><given-names>P</given-names></name><name><surname>Amstutz</surname><given-names>H</given-names></name><name><surname>Kohli</surname><given-names>J</given-names></name><name><surname>Leupold</surname><given-names>U</given-names></name></person-group><year>1982</year><article-title>Recombination between dispersed serine tRNA genes in <italic>Schizosaccharomyces pombe</italic></article-title><source>Nature</source><volume>300</volume><fpage>225</fpage><lpage>231</lpage><pub-id pub-id-type="doi">10.1038/300225a0</pub-id></element-citation></ref><ref id="bib42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Olejniczak</surname><given-names>M</given-names></name><name><surname>Dale</surname><given-names>T</given-names></name><name><surname>Fahlman</surname><given-names>RP</given-names></name><name><surname>Uhlenbeck</surname><given-names>OC</given-names></name></person-group><year>2005</year><article-title>Idiosyncratic tuning of tRNAs to achieve uniform ribosome binding</article-title><source>Nature Structural & Molecular Biology</source><volume>12</volume><fpage>788</fpage><lpage>793</lpage><pub-id pub-id-type="doi">10.1038/nsmb978</pub-id></element-citation></ref><ref id="bib43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Parmley</surname><given-names>JL</given-names></name><name><surname>Hurst</surname><given-names>LD</given-names></name></person-group><year>2007</year><article-title>Exonic splicing regulatory elements skew synonymous codon usage near intron-exon boundaries in mammals</article-title><source>Molecular Biology and Evolution</source><volume>24</volume><fpage>1600</fpage><lpage>1603</lpage><pub-id pub-id-type="doi">10.1093/molbev/msm104</pub-id></element-citation></ref><ref id="bib44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pavon-Eternod</surname><given-names>M</given-names></name><name><surname>David</surname><given-names>A</given-names></name><name><surname>Dittmar</surname><given-names>K</given-names></name><name><surname>Berglund</surname><given-names>P</given-names></name><name><surname>Pan</surname><given-names>T</given-names></name><name><surname>Bennink</surname><given-names>JR</given-names></name><name><surname>Yewdell</surname><given-names>JW</given-names></name></person-group><year>2013a</year><article-title>Vaccinia and influenza A viruses select rather than adjust tRNAs to optimize translation</article-title><source>Nucleic Acids Research</source><volume>41</volume><fpage>1914</fpage><lpage>1921</lpage><pub-id pub-id-type="doi">10.1093/nar/gks986</pub-id></element-citation></ref><ref id="bib45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pavon-Eternod</surname><given-names>M</given-names></name><name><surname>Gomes</surname><given-names>S</given-names></name><name><surname>Rosner</surname><given-names>MR</given-names></name><name><surname>Pan</surname><given-names>T</given-names></name></person-group><year>2013b</year><article-title>Overexpression of initiator methionine tRNA leads to global reprogramming of tRNA expression and increased proliferation in human epithelial cells</article-title><source>RNA</source><volume>19</volume><fpage>461</fpage><lpage>466</lpage><pub-id pub-id-type="doi">10.1261/rna.037507.112</pub-id></element-citation></ref><ref id="bib46"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Pechmann</surname><given-names>S</given-names></name><name><surname>Frydman</surname><given-names>J</given-names></name></person-group><year>2012</year><article-title>Evolutionary conservation of codon optimality reveals hidden signatures of cotranslational folding</article-title><source>Nature Structural & Molecular Biology</source><volume>20</volume><fpage>1</fpage><lpage>8</lpage><pub-id pub-id-type="doi">10.1038/nsmb.2466</pub-id></element-citation></ref><ref id="bib47"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Percudani</surname><given-names>R</given-names></name><name><surname>Pavesi</surname><given-names>A</given-names></name><name><surname>Ottonello</surname><given-names>S</given-names></name></person-group><year>1997</year><article-title>Transfer RNA gene redundancy and translational selection in Saccharomyces cerevisiae</article-title><source>Journal of Molecular Biology</source><volume>268</volume><fpage>322</fpage><lpage>330</lpage><pub-id pub-id-type="doi">10.1006/jmbi.1997.0942</pub-id></element-citation></ref><ref id="bib48"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rao</surname><given-names>Y</given-names></name><name><surname>Wu</surname><given-names>G</given-names></name><name><surname>Wang</surname><given-names>Z</given-names></name><name><surname>Chai</surname><given-names>X</given-names></name><name><surname>Nie</surname><given-names>Q</given-names></name><name><surname>Zhang</surname><given-names>X</given-names></name></person-group><year>2011</year><article-title>Mutation bias is the driving force of codon usage in the Gallus gallus genome</article-title><source>DNA Research</source><volume>18</volume><fpage>499</fpage><lpage>512</lpage><pub-id pub-id-type="doi">10.1093/dnares/dsr035</pub-id></element-citation></ref><ref id="bib49"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rawlings</surname><given-names>TA</given-names></name><name><surname>Collins</surname><given-names>TM</given-names></name><name><surname>Bieler</surname><given-names>R</given-names></name></person-group><year>2003</year><article-title>Changing identities: tRNA duplication and remolding within animal mitochondrial genomes</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>100</volume><fpage>15700</fpage><lpage>15705</lpage><pub-id pub-id-type="doi">10.1073/pnas.2535036100</pub-id></element-citation></ref><ref id="bib50"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Rogers</surname><given-names>HH</given-names></name><name><surname>Bergman</surname><given-names>CM</given-names></name><name><surname>Griffiths-Jones</surname><given-names>S</given-names></name></person-group><year>2010</year><article-title>The evolution of tRNA genes in <italic>Drosophila</italic></article-title><source>Genome Biology and Evolution</source><volume>2</volume><fpage>467</fpage><lpage>477</lpage><pub-id pub-id-type="doi">10.1093/gbe/evq034</pub-id></element-citation></ref><ref id="bib51"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Saikia</surname><given-names>M</given-names></name><name><surname>Krokowski</surname><given-names>D</given-names></name><name><surname>Guan</surname><given-names>BJ</given-names></name><name><surname>Ivanov</surname><given-names>P</given-names></name><name><surname>Parisien</surname><given-names>M</given-names></name><name><surname>Hu</surname><given-names>GF</given-names></name><name><surname>Anderson</surname><given-names>P</given-names></name><name><surname>Pan</surname><given-names>T</given-names></name><name><surname>Hatzoglou</surname><given-names>M</given-names></name></person-group><year>2012</year><article-title>Genome-wide identification and quantitative analysis of cleaved tRNA fragments induced by cellular stress</article-title><source>Journal of Biological Chemistry</source><volume>287</volume><fpage>42708</fpage><lpage>42725</lpage><pub-id pub-id-type="doi">10.1074/jbc.M112.371799</pub-id></element-citation></ref><ref id="bib52"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Saks</surname><given-names>ME</given-names></name><name><surname>Conery</surname><given-names>JS</given-names></name></person-group><year>2007</year><article-title>Anticodon-dependent conservation of bacterial tRNA gene sequences</article-title><source>RNA</source><volume>13</volume><fpage>651</fpage><lpage>660</lpage><pub-id pub-id-type="doi">10.1261/rna.345907</pub-id></element-citation></ref><ref id="bib53"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Saks</surname><given-names>ME</given-names></name><name><surname>Sampson</surname><given-names>JR</given-names></name><name><surname>Abelson</surname><given-names>J</given-names></name></person-group><year>1998</year><article-title>Evolution of a transfer RNA gene through a point mutation in the anticodon</article-title><source>Science</source><volume>279</volume><fpage>1665</fpage><lpage>1670</lpage><pub-id pub-id-type="doi">10.1126/science.279.5357.1665</pub-id></element-citation></ref><ref id="bib54"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schmeing</surname><given-names>TM</given-names></name><name><surname>Voorhees</surname><given-names>RM</given-names></name><name><surname>Kelley</surname><given-names>AC</given-names></name><name><surname>Ramakrishnan</surname><given-names>V</given-names></name></person-group><year>2011</year><article-title>How mutations in tRNA distant from the anticodon affect the fidelity of decoding</article-title><source>Nature Structural & Molecular Biology</source><volume>18</volume><fpage>432</fpage><lpage>436</lpage><pub-id pub-id-type="doi">10.1038/nsmb.2003</pub-id></element-citation></ref><ref id="bib55"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schultz</surname><given-names>DW</given-names></name><name><surname>Yarus</surname><given-names>M</given-names></name></person-group><year>1994</year><article-title>tRNA structure and ribosomal function. I. tRNA nucleotide 27-43 mutations enhance first position wobble</article-title><source>Journal of Molecular Biology</source><volume>235</volume><fpage>1381</fpage><lpage>1394</lpage><pub-id pub-id-type="doi">10.1006/jmbi.1994.1095</pub-id></element-citation></ref><ref id="bib56"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sharp</surname><given-names>PM</given-names></name><name><surname>Li</surname><given-names>WH</given-names></name></person-group><year>1986a</year><article-title>An evolutionary perspective on synonymous codon usage in unicellular organisms</article-title><source>Journal of Molecular Evolution</source><volume>24</volume><fpage>28</fpage><lpage>38</lpage><pub-id pub-id-type="doi">10.1007/BF02099948</pub-id></element-citation></ref><ref id="bib57"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sharp</surname><given-names>PM</given-names></name><name><surname>Li</surname><given-names>WH</given-names></name></person-group><year>1986b</year><article-title>Codon usage in regulatory genes in Escherichia coli does not reflect selection for “rare” codons</article-title><source>Nucleic Acids Research</source><volume>14</volume><fpage>7737</fpage><lpage>7749</lpage><pub-id pub-id-type="doi">10.1093/nar/14.19.7737</pub-id></element-citation></ref><ref id="bib58"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sharp</surname><given-names>PM</given-names></name><name><surname>Li</surname><given-names>WH</given-names></name></person-group><year>1987</year><article-title>The codon Adaptation Index–a measure of directional synonymous codon usage bias, and its potential applications</article-title><source>Nucleic Acids Research</source><volume>15</volume><fpage>1281</fpage><lpage>1295</lpage><pub-id pub-id-type="doi">10.1093/nar/15.3.1281</pub-id></element-citation></ref><ref id="bib59"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stamatakis</surname><given-names>A</given-names></name></person-group><year>2006</year><article-title>RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models</article-title><source>Bioinformatics</source><volume>22</volume><fpage>2688</fpage><lpage>2690</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btl446</pub-id></element-citation></ref><ref id="bib60"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Stoletzki</surname><given-names>N</given-names></name><name><surname>Eyre-Walker</surname><given-names>A</given-names></name></person-group><year>2007</year><article-title>Synonymous codon usage in <italic>Escherichia coli</italic>: selection for translational accuracy</article-title><source>Molecular Biology and Evolution</source><volume>24</volume><fpage>374</fpage><lpage>381</lpage><pub-id pub-id-type="doi">10.1093/molbev/msl166</pub-id></element-citation></ref><ref id="bib61"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sukumaran</surname><given-names>J</given-names></name><name><surname>Holder</surname><given-names>MT</given-names></name></person-group><year>2010</year><article-title>DendroPy: a Python library for phylogenetic computing</article-title><source>Bioinformatics</source><volume>26</volume><fpage>1569</fpage><lpage>1571</lpage><pub-id pub-id-type="doi">10.1093/bioinformatics/btq228</pub-id></element-citation></ref><ref id="bib62"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Teshima</surname><given-names>KM</given-names></name><name><surname>Innan</surname><given-names>H</given-names></name></person-group><year>2004</year><article-title>The effect of gene conversion on the divergence between duplicated genes</article-title><source>Genetics</source><volume>166</volume><fpage>1553</fpage><lpage>1560</lpage><pub-id pub-id-type="doi">10.1534/genetics.166.3.1553</pub-id></element-citation></ref><ref id="bib63"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thanaraj</surname><given-names>TA</given-names></name><name><surname>Argos</surname><given-names>P</given-names></name></person-group><year>1996</year><article-title>Ribosome-mediated translational pause and protein domain organization</article-title><source>Protein Science</source><volume>5</volume><fpage>1594</fpage><lpage>1612</lpage><pub-id pub-id-type="doi">10.1002/pro.5560050814</pub-id></element-citation></ref><ref id="bib64"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tsai</surname><given-names>C-J</given-names></name><name><surname>Sauna</surname><given-names>ZE</given-names></name><name><surname>Kimchi-Sarfaty</surname><given-names>C</given-names></name><name><surname>Ambudkar</surname><given-names>SV</given-names></name><name><surname>Gottesman</surname><given-names>MM</given-names></name><name><surname>Nussinov</surname><given-names>R</given-names></name></person-group><year>2008</year><article-title>Synonymous mutations and ribosome stalling can lead to altered folding pathways and distinct minima</article-title><source>Journal of Molecular Biology</source><volume>383</volume><fpage>281</fpage><lpage>291</lpage><pub-id pub-id-type="doi">10.1016/j.jmb.2008.08.012</pub-id></element-citation></ref><ref id="bib65"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tuller</surname><given-names>T</given-names></name><name><surname>Carmi</surname><given-names>A</given-names></name><name><surname>Vestsigian</surname><given-names>K</given-names></name><name><surname>Navon</surname><given-names>S</given-names></name><name><surname>Dorfan</surname><given-names>Y</given-names></name><name><surname>Zaborske</surname><given-names>J</given-names></name><name><surname>Pan</surname><given-names>T</given-names></name><name><surname>Dahan</surname><given-names>O</given-names></name><name><surname>Furman</surname><given-names>I</given-names></name><name><surname>Pilpel</surname><given-names>Y</given-names></name></person-group><year>2010</year><article-title>An evolutionarily conserved mechanism for controlling the efficiency of protein translation</article-title><source>Cell</source><volume>141</volume><fpage>344</fpage><lpage>354</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2010.03.031</pub-id></element-citation></ref><ref id="bib66"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Urrutia</surname><given-names>AO</given-names></name><name><surname>Hurst</surname><given-names>LD</given-names></name></person-group><year>2001</year><article-title>Codon usage bias covaries with expression breadth and the rate of synonymous evolution in humans, but this is not evidence for selection</article-title><source>Genetics</source><volume>159</volume><fpage>1191</fpage><lpage>1199</lpage></element-citation></ref><ref id="bib67"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Varenne</surname><given-names>S</given-names></name><name><surname>Buc</surname><given-names>J</given-names></name><name><surname>Lloubes</surname><given-names>R</given-names></name><name><surname>Lazdunski</surname><given-names>C</given-names></name></person-group><year>1984</year><article-title>Translation is a non-uniform process. Effect of tRNA availability on the rate of elongation of nascent polypeptide chains</article-title><source>Journal of Molecular Biology</source><volume>180</volume><fpage>549</fpage><lpage>576</lpage><pub-id pub-id-type="doi">10.1016/0022-2836(84)90027-5</pub-id></element-citation></ref><ref id="bib68"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Von Pawel-Rammingen</surname><given-names>U</given-names></name><name><surname>Aström</surname><given-names>S</given-names></name><name><surname>Byström</surname><given-names>AS</given-names></name></person-group><year>1992</year><article-title>Mutational analysis of conserved positions potentially important for initiator tRNA function in <italic>Saccharomyces cerevisiae</italic></article-title><source>Molecular and Cellular Biology</source><volume>12</volume><fpage>1432</fpage><lpage>1442</lpage></element-citation></ref><ref id="bib69"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Widmann</surname><given-names>J</given-names></name><name><surname>Harris</surname><given-names>J</given-names></name><name><surname>Lozupone</surname><given-names>C</given-names></name></person-group><year>2010</year><article-title>Stable tRNA-based phylogenies using only 76 nucleotides</article-title><source>RNA</source><fpage>1469</fpage><lpage>1477</lpage><pub-id pub-id-type="doi">10.1261/rna.726010</pub-id></element-citation></ref><ref id="bib70"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wilke</surname><given-names>CO</given-names></name><name><surname>Drummond</surname><given-names>DA</given-names></name></person-group><year>2010</year><article-title>Signatures of protein biophysics in coding sequence evolution</article-title><source>Current Opinion In Structural Biology</source><volume>20</volume><fpage>385</fpage><lpage>389</lpage><pub-id pub-id-type="doi">10.1016/j.sbi.2010.03.004</pub-id></element-citation></ref><ref id="bib71"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Withers</surname><given-names>M</given-names></name><name><surname>Wernisch</surname><given-names>L</given-names></name><name><surname>dos Reis</surname><given-names>M</given-names></name></person-group><year>2006</year><article-title>Archaeology and evolution of transfer RNA genes in the <italic>Escherichia coli</italic> genome</article-title><source>RNA</source><volume>12</volume><fpage>933</fpage><lpage>942</lpage><pub-id pub-id-type="doi">10.1261/rna.2272306</pub-id></element-citation></ref><ref id="bib72"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wolf</surname><given-names>YI</given-names></name><name><surname>Rogozin</surname><given-names>IB</given-names></name><name><surname>Kondrashov</surname><given-names>AS</given-names></name><name><surname>Koonin</surname><given-names>EV</given-names></name></person-group><year>2001</year><article-title>Genome alignment, evolution of prokaryotic genome organization, and prediction of gene function using genomic context</article-title><source>Genome Research</source><volume>11</volume><fpage>356</fpage><lpage>372</lpage><pub-id pub-id-type="doi">10.1101/gr.161901</pub-id></element-citation></ref><ref id="bib73"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname><given-names>G</given-names></name><name><surname>Hubalewska</surname><given-names>M</given-names></name><name><surname>Ignatova</surname><given-names>Z</given-names></name></person-group><year>2009</year><article-title>Transient ribosomal attenuation coordinates protein synthesis and co-translational folding</article-title><source>Nature Structural & Molecular Biology</source><volume>16</volume><fpage>274</fpage><lpage>280</lpage><pub-id pub-id-type="doi">10.1038/nsmb.1554</pub-id></element-citation></ref><ref id="bib74"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname><given-names>T</given-names></name><name><surname>Weems</surname><given-names>M</given-names></name><name><surname>Wilke</surname><given-names>CO</given-names></name></person-group><year>2009</year><article-title>Translationally optimal codons associate with structurally sensitive sites in proteins</article-title><source>Molecular Biology and Evolution</source><volume>26</volume><fpage>1571</fpage><lpage>1580</lpage><pub-id pub-id-type="doi">10.1093/molbev/msp070</pub-id></element-citation></ref><ref id="bib75"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname><given-names>M</given-names></name><name><surname>Guo</surname><given-names>J</given-names></name><name><surname>Cha</surname><given-names>J</given-names></name><name><surname>Chae</surname><given-names>M</given-names></name><name><surname>Chen</surname><given-names>S</given-names></name><name><surname>Barral</surname><given-names>JM</given-names></name><name><surname>Sachs</surname><given-names>MS</given-names></name><name><surname>Liu</surname><given-names>Y</given-names></name></person-group><year>2013</year><article-title>Non-optimal codon usage affects expression, structure and function of clock protein FRQ</article-title><source>Nature</source><volume>495</volume><fpage>111</fpage><lpage>115</lpage><pub-id pub-id-type="doi">10.1038/nature11833</pub-id></element-citation></ref></ref-list></back><sub-article article-type="article-commentary" id="SA1"><front-stub><article-id pub-id-type="doi">10.7554/eLife.01339.016</article-id><title-group><article-title>Decision letter</article-title></title-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Laub</surname><given-names>Michael</given-names></name><role>Reviewing editor</role><aff><institution>Massachusetts Institute of Technology</institution>, <country>United States</country></aff></contrib></contrib-group></front-stub><body><boxed-text><p>eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see <ext-link ext-link-type="uri" xlink:href="http://elife.elifesciences.org/review-process">review process</ext-link>). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.</p></boxed-text><p>Thank you for sending your work entitled “tRNA genes rapidly change in evolution to meet novel translational demands” for consideration at <italic>eLife</italic>. Your article has been favorably evaluated by a Senior editor and 3 reviewers, one of whom is a member of our Board of Reviewing Editors.</p><p>The following individuals responsible for the peer review of your submission have agreed to reveal their identity: Michael Laub (Reviewing editor) and Chris Marx (peer reviewer).</p><p>The Reviewing editor and the other reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.</p><p>Overall, this paper was favorably received by all three reviewers. The general consensus was that the results provide important new insights into the evolution of tRNAs and represents a nice example of the power of experimental evolution to provide novel biological insights. There are, however, several concerns about the writing and interpretation of a few things that we felt warranted revision. More generally, there was a concern that the paper felt like a loose collection of vignettes and better writing to integrate them would serve the reader well. Each of these concerns about the writing must be dealt with through changes in the text. There is also one experiment that is straightforward that would be necessary to bolster one of the primary claims made in the text (item #1 below) and one experiment (item #6 below) that a reviewer suggests (it's not essential, but could really enhance the paper). All of the major issues to address, collated from the three independent reviewers, are listed below:</p><p>1) The authors argue that overproducing low-copy tRNAs causes proteotoxic stress. This conclusion rests primarily on a single assay involving the expression of an mCherry-tagged human protein. I think it's important to assess this by at least one other means; for example, testing whether various HSPs are induced (either individual HSPs or, preferably, a genome-wide transcriptional analysis by arrays or RNA-Seq).</p><p>(As a side note: it was not clear whether the potential proteotoxic stress observed resulted from misincorporation of amino acids per se or defects in co-translational folding that result from changes in translation speed/efficiency without misincorporation. The reviewers did not think experiments to address this were necessary for a revised manuscript, but we felt future experiments ought to try and distinguish between these models. For instance, one could measure rates of misincorporation. Or, if the defect is in cotranslational folding, one could imagine taking a protein in which a “slow” and a “fast” codon are near each other: it should be true (if cotranslational folding is the culprit) that making both codons “fast” is detrimental, but that swapping them (and retaining one slow, one fast) is tolerated. The misincorporation model would, by contrast, predict that the exact position of the “slow” codon matters.)</p><p>2) The authors often refer to “translation efficiency” in the text, but does this mean translation speed, accuracy, or some combination thereof? Can you be more explicit about the intended meaning?</p><p>3) The data are clear and compelling that “chimeric” tRNAs harboring the backbone of one tRNA and the anticodon of another can compete with wild-type under a range of conditions. But this begs the question of why duplicated tRNAs in yeast retain nearly 100% sequence identity - one would expect, given the data presented here, to see much more sequence diversity. Can this be commented on?</p><p>4) The first paragraph introduces codon bias, then states that the existence of codon usage implies that translation has undergone adaptation. This is not the case, because mutational biases are known to cause codon usage bias. The manuscript states that adaptation is “presumably” to fine-tune codon usage to gene expression. One need not presume, as the subject has been studied for decades. The authors should expand this section to give a more accurate picture of the question at hand. Plotkin & Kudla (Nat Rev Genet 2011) provide a solid review.</p><p>5) It has not been established that any tRNA scaffold with an anticodon change suffices to restore function and fitness, only isoaccepting scaffolds.</p><p>One would not expect the scaffold to be entirely interchangeable. Some aa-tRNA synthetases do not bind the anticodon (Leu, Ala and Ser-accepting aaRS), and the anticodon of the corresponding tRNAs can be changed without affecting charging (see Beuning and Musier-Forsyth, Biopolymers 1999). For these species, the scaffold is decisive for proper charging.</p><p>Thus, the claim that the tRNA scaffold does not confer tRNA functionality is incorrect. It also does not follow that the anticodon is the dominant feature. It may be the dominant feature so long as the tRNA scaffold is from a synonymous family.</p><p>6) Not reporting results after starving the cells for arginine seems a missed opportunity.</p><p>7) The levels of tRNAs achieved after multicopy plasmid expression remain unclear. At present, we do not know whether the fold-change in tRNA for each tRNA variant is equivalent, or whether the absolute change is equivalent. We also have no guidance as to which one matters. It would seem that if achieving an increase of x molecules of tRNA means a 10% increase in the abundant tRNA and a 10-fold increase in the rare tRNA, then the results obtained may not be surprising. Please address why adding a putatively constant amount of additional tRNA is the correct experiment.</p><p>8) The paper jumps to assuming codon usage bias emerges from selection before providing evidence that it isn't just random drift. Particularly for a broad audience like <italic>eLife</italic>, I think it is critical to distinguish that the existence of a bias is not enough to claim selection. It is the ways in which there are bias (such as between genes and within them, not just between organisms) that make it clear adaptation plays a role. This idea shows up at the end of the paragraph beginning “The fitness effects of an unmet translational demand…”, but by then adaptation has long since been assumed.</p><p>9) Related to the above point, the paragraph beginning “Furthermore, the forces that direct and maintain low copy tRNA families remain unclear…” is confusing as currently written. I think the point trying to be made is that, although on average selection favors a particular set of codons that are generally labeled as optimal, might there be an adaptive reason for low frequency/sub-optimal codons that – in certain circumstances – are actually selectively advantageous. Please try to clarify this point, as it is so central to the paper. Along these lines, the argument could be strengthened by citing one or more of the available papers that a gene encoded with 100% ‘optimal’ codons can actually be poorly expressed or cause decreased fitness.</p><p>10) An interesting point was made that detectable anticodon switching events occurred in genomes with a greater than average number of tRNA. Might there be ascertainment bias here? If such switches are easiest to detect for clades with many sequenced members, and these are generally fast pathogens with many tRNA (and rRNA copies, etc), then this may make at least some of this signal a simple consequence of where it was easiest to identify the events.</p><p>11) In their Discussion, they say, “our…search…revealed the prevalence of this adaptation mechanism in nature,” where I think they can at best say “It is not established that natural anticodon changes are adaptive, nor that, if they are adaptive, that they are a response to novel translational demands. Our search revealed the prevalence of tRNA anticodon mutations in nature, consistent with a response to novel translational demands as in our experimental evolution results”, or something of the sort.</p></body></sub-article><sub-article article-type="reply" id="SA2"><front-stub><article-id pub-id-type="doi">10.7554/eLife.01339.017</article-id><title-group><article-title>Author response</article-title></title-group></front-stub><body><p><italic>1) The authors argue that overproducing low-copy tRNAs causes proteotoxic stress. This conclusion rests primarily on a single assay involving the expression of an mCherry-tagged human protein. I think it's important to assess this by at least one other means; for example, testing whether various HSPs are induced (either individual HSPs or, preferably, a genome-wide transcriptional analysis by arrays or RNA-Seq)</italic>.</p><p>Following the reviewers’ suggestion, we measured the expression level of six canonical heat shock proteins (HSPs) by RT-qPCR in order to establish the proteotoxic stress in the two tRNA overexpressing strains. Our new results (<xref ref-type="fig" rid="fig4">Figure 4E</xref>) demonstrate that indeed HSP genes are induced 2–8-fold when overexpressing the originally low copy number tRNA gene. In contrast, none of the HSP genes was induced in the cells that overexpress the high copy number tRNA gene. Thus, the RT-PCR measurements support our initial claim that over-expressing a rare tRNA gene causes a proteotoxic stress, ascribing a pro-folding function to rare tRNAs, and suggest that at least part of the fitness reduction observed upon over expression of a rare tRNA is caused by proteotoxic stress.</p><p><italic>(As a side note: it was not clear whether the potential proteotoxic stress observed resulted from misincorporation of amino acids</italic> per se <italic>or defects in co-translational folding that result from changes in translation speed/efficiency without misincorporation. The reviewers did not think experiments to address this were necessary for a revised manuscript, but we felt future experiments ought to try and distinguish between these models. For instance, one could measure rates of misincorporation. Or, if the defect is in cotranslational folding, one could imagine taking a protein in which a “slow” and a “fast” codon are near each other: it should be true (if cotranslational folding is the culprit) that making both codons “fast” is detrimental, but that swapping them (and retaining one slow, one fast) is tolerated. The misincorporation model would, by contrast, predict that the exact position of the “slow” codon matters.</italic>)</p><p>We agree that the proteotoxic stress may have several origins and thus discuss three mechanisms that may lead to it in the text. We also further clarify that the potential part of the observed proteotoxic stress due to mis-incorporation still remains to be studied.</p><p>Currently in our lab, we aim to pursue this line of research by protein sequencing. We use mass spectrometry technology on designated proteins enriched with arginine. We also thank the reviews for the suggested “slow” and “fast” adjacent arginine codons as it is an elegant experiment that can be performed in future by using the same mass-spec method.</p><p><italic>2) The authors often refer to “translation efficiency” in the text, but does this mean translation speed, accuracy, or some combination thereof? Can you be more explicit about the intended meaning</italic>?</p><p>We now clarify in the text the definition of translation efficiency. We note that the translation efficiency affects both production throughput and accuracy.</p><p><italic>3) The data are clear and compelling that “chimeric” tRNAs harboring the backbone of one tRNA and the anticodon of another can compete with wild-type under a range of conditions. But this begs the question of why duplicated tRNAs in yeast retain nearly 100% sequence identity - one would expect, given the data presented here, to see much more sequence diversity. Can this be commented on</italic>?</p><p>We agree that in light of our results the origin of the high sequence identity among tRNA genes is an interesting question. We added the following paragraph to the Discussion as follows:</p><p>“The present chimeric tRNAs that emerged in our lab-evolution experiments successfully replaced the deleted tRNA, despite differences in 20 nucleotide positions between the two tRNA scaffolds…”</p><p><italic>4) The first paragraph introduces codon bias, then states that the existence of codon usage implies that translation has undergone adaptation. This is not the case, because mutational biases are known to cause codon usage bias. The manuscript states that adaptation is “presumably” to fine-tune codon usage to gene expression. One need not presume, as the subject has been studied for decades. The authors should expand this section to give a more accurate picture of the question at hand. Plotkin & Kudla (Nat Rev Genet 2011) provide a solid review</italic>.</p><p>We thank the reviewers for this important point. We extensively edited the first three paragraphs of the Introduction and elaborated on the evolutionary processes of codon usage bias.</p><p><italic>5) It has not been established that any tRNA scaffold with an anticodon change suffices to restore function and fitness, only isoaccepting scaffolds</italic>.</p><p>One would not expect the scaffold to be entirely interchangeable. Some aa-tRNA synthetases do not bind the anticodon (Leu, Ala and Ser-accepting aaRS), and the anticodon of the corresponding tRNAs can be changed without affecting charging (see Beuning and Musier-Forsyth, Biopolymers 1999). For these species, the scaffold is decisive for proper charging.</p><p><italic>Thus, the claim that the tRNA scaffold does not confer tRNA functionality is incorrect. It also does not follow that the anticodon is the dominant feature. It may be the dominant feature so long as the tRNA scaffold is from a synonymous family</italic>.</p><p>We agree that our observations focus on synonymous tRNA families only. Thus, we changed the text accordingly:</p><p>“Hence, the chimeric tRNA provides a direct in vivo indication that the scaffolds of tRNAs, which encode for the same amino-acid, may be interchangeable in terms of their effect on cellular growth under the conditions we tested.”</p><p><italic>6) Not reporting results after starving the cells for arginine seems a missed opportunity</italic>.</p><p>We now preformed arginine-starving assays on the deletion strain (ΔtRNA<sup>Arg</sup><sub>CCU</sub>) compared to the WT (not shown). Yet, we did not observe any significant difference in the growth, comparing this medium to a medium containing arginine. Accordingly, we did not pursue this line of research further.</p><p><italic>7) The levels of tRNAs achieved after multicopy plasmid expression remain unclear. At present, we do not know whether the fold-change in tRNA for each tRNA variant is equivalent, or whether the absolute change is equivalent. We also have no guidance as to which one matters. It would seem that if achieving an increase of x molecules of tRNA means a 10% increase in the abundant tRNA and a 10-fold increase in the rare tRNA, then the results obtained may not be surprising. Please address why adding a putatively constant amount of additional tRNA is the correct experiment</italic>.</p><p>We acknowledge that this is a valid point, yet with the aim to recapitulates the natural evolutionary process of gene dosage changes we used similar increments in gene copy number for both rare and abundant tRNA. We started with addition of a centromeric plasmid (considered as an equivalent to addition of a single copy) and observed minor effects under severe growth conditions only. Therefore, we aimed to shift tRNA copies of both families to the extreme. For that we used high-copy plasmids in which both families were supplemented with the same number of copies. Since expression of each tRNA is determined by its endogenous promoter we could not control for differences in tRNA expression, only for differences in tRNA gene copy number. It is also worth noting that the copy number of the plasmid that we used (2μ plasmid) ranges from 20 to 80 in a single cell. Thus, we believe that the expression level of both rare and abundant tRNAs were greatly influenced in our experiments.</p><p><italic>8) The paper jumps to assuming codon usage bias emerges from selection before providing evidence that it isn't just random drift. Particularly for a broad audience like</italic> eLife<italic>, I think it is critical to distinguish that the existence of a bias is not enough to claim selection. It is the ways in which there are bias (such as between genes and within them, not just between organisms) that make it clear adaptation plays a role. This idea shows up at the end of the paragraph beginning “The fitness effects of an unmet translational demand…”, but by then adaptation has long since been assumed</italic>.</p><p>We thank the reviewers for this important point and we have now edited the first three paragraphs of the Introduction to clarify the role of selection and drift in the phenomenon of codon usage bias.</p><p><italic>9) Related to the above point, the paragraph beginning “Furthermore, the forces that direct and maintain low copy tRNA families remain unclear…” is confusing as currently written. I think the point trying to be made is that, although on average selection favors a particular set of codons that are generally labeled as optimal, might there be an adaptive reason for low frequency/sub-optimal codons that</italic> – <italic>in certain circumstances</italic> – <italic>are actually selectively advantageous. Please try to clarify this point, as it is so central to the paper. Along these lines, the argument could be strengthened by citing one or more of the available papers that a gene encoded with 100% ‘optimal’ codons can actually be poorly expressed or cause decreased fitness</italic>.</p><p>We have revised this paragraph completely in order to express the point in a clearer fashion. We also added published examples, according to the reviewers’ suggestion, that demonstrate how “non-optimal codons” have advantageous functions in translation. The new paragraph now reads as follows:</p><p>“Further, the forces that direct and maintain low copy tRNA families remain unclear…”</p><p><italic>10) An interesting point was made that detectable anticodon switching events occurred in genomes with a greater than average number of tRNA. Might there be ascertainment bias here? If such switches are easiest to detect for clades with many sequenced members, and these are generally fast pathogens with many tRNA (and rRNA copies, etc), then this may make at least some of this signal a simple consequence of where it was easiest to identify the events</italic>.</p><p>Indeed, we observed that genomes with large tRNA pools demonstrate more tRNA switching events. However, thanks to the reviewer comments, we noticed that there are no simple relations between the number of tRNA genes and the number of tRNA switching events per species. We thus took out our original statement about the difference between the species and merely state that this issue warrants future analysis.</p><p>Additionally, it is worth noting that our tRNA switching analysis was performed on individual species and not clades and that all tRNAs are discovered by scanning whole genomic sequences with the tRNA-SCAN software and then scored with the COVE algorithm for the likelihood of being a tRNA gene. Thus, even though there might be a bias in the available genomic data (the species sequenced so far) our screen is able to detect most of the switching event per species.</p><p>We thank the reviewers for raising this important point and clarify the issue in the text of both the Results and Materials and methods sections.</p><p><italic>11) In their Discussion, they say, “our…search…revealed the prevalence of this adaptation mechanism in nature,” where I think they can at best say “It is not established that natural anticodon changes are adaptive, nor that, if they are adaptive, that they are a response to novel translational demands. Our search revealed the prevalence of tRNA anticodon mutations in nature, consistent with a response to novel translational demands as in our experimental evolution results”, or something of the sort</italic>.</p><p>We changed the text according to the reviewers’ comment by modifying the end of the paragraph on as follows:</p><p>“Furthermore, our systematic search for tRNA switching events throughout the tree of life revealed the prevalence of tRNA anticodon mutations in nature. This observation is consistent with the results of our lab-evolution experiment and may be the evolutionary outcome to novel translational demands in the wild.”</p></body></sub-article></article> |