Coevolution of codon usage and tRNA genes leads to alternative stable states of biased codon usage

The typical number of tRNA genes in bacterial genomes is around 50, but this number varies from under 30 to over 120. We argue that tRNA gene copy numbers evolve in response to translational selection. In rapidly multiplying organisms, the time spent in translation is a limiting factor in cell division; hence, it pays to duplicate tRNA genes, thereby increasing the concentration of tRNA molecules in the cell and speeding up translation. In slowly multiplying organisms, translation time is not a limiting factor, so the overall translational cost is minimized by reducing the tRNAs to only one copy of each required gene. Translational selection also causes a preference for codons that are most rapidly translated by the current tRNAs; hence, codon usage and tRNA gene content will coevolve to a state where each is adapted to the other. We show that there is often more than one stable coevolved state. This explains why different combinations of tRNAs and codon bias can exist for different amino acids in the same organism. We analyze a set of 80 complete bacterial genomes and show that the theory predicts many of the trends that are seen in these data.     

tRNA gene copy number variation in humans

The human tRNAome consists of more than 500 interspersed tRNA genes comprising 51 anticodon families of largely unequal copy number. We examined tRNA gene copy number variation (tgCNV) in six individuals; two kindreds of two parents and a child, using high coverage whole genome sequence data. Such differences may be important because translation of some mRNAs is sensitive to the relative amounts of tRNAs and because tRNA competition determines translational efficiency vs. fidelity and production of native vs. misfolded proteins. We identified several tRNA gene clusters with CNV, which in some cases were part of larger iterations. In addition there was an isolated tRNALysCUU gene that was absent as a homozygous deletion in one of the parents. When assessed by semiquantitative PCR in 98 DNA samples representing a wide variety of ethnicities, this allele was found deleted in hetero- or homozygosity in all groups at ~ 50% frequency. This is the first report of copy number variation of human tRNA genes. We conclude that tgCNV exists at significant levels among individual humans and discuss the results in terms of genetic diversity and prior genome wide association studies (GWAS) that suggest the importance of the ratio of tRNALys isoacceptors in Type-2 diabetes.     

Altered tRNA characteristics and 3' maturation in bacterial symbionts with reduced genomes

Translational efficiency is controlled by tRNAs and other genome-encoded mechanisms. In organelles, translational processes are dramatically altered because of genome shrinkage and horizontal acquisition of gene products. The influence of genome reduction on translation in endosymbionts is largely unknown. Here, we investigate whether divergent lineages of Buchnera aphidicola, the reduced-genome bacterial endosymbiont of aphids, possess altered translational features compared with their free-living relative, Escherichia coli. Our RNAseq data support the hypothesis that translation is less optimal in Buchnera than in E. coli. We observed a specific, convergent, pattern of tRNA loss in Buchnera and other endosymbionts that have undergone genome shrinkage. Furthermore, many modified nucleoside pathways that are important for E. coli translation are lost in Buchnera. Additionally, Buchnera’s A + T compositional bias has resulted in reduced tRNA thermostability, and may have altered aminoacyl-tRNA synthetase recognition sites. Buchnera tRNA genes are shorter than those of E. coli, as the majority no longer has a genome-encoded 3'' CCA; however, all the expressed, shortened tRNAs undergo 3′ CCA maturation. Moreover, expression of tRNA isoacceptors was not correlated with the usage of corresponding codons. Overall, our data suggest that endosymbiont genome evolution alters tRNA characteristics that are known to influence translational efficiency in their free-living relative.     

Over expression of a tRNA(Leu) isoacceptor changes charging pattern of leucine tRNAs and reveals new codon reading

During mRNA translation, synonymous codons for one amino acid are often read by different isoaccepting tRNAs. The theory of selective tRNA charging predicts greatly varying percentages of aminoacylation among isoacceptors in cells starved for their common amino acid. It also predicts major changes in tRNA charging patterns upon concentration changes of single isoacceptors, which suggests a novel type of translational control of gene expression. We therefore tested the theory by measuring with Northern blots the charging of Leu-tRNAs in Escherichia coli under Leu limitation in response to over expression of tRNA(GAG)(Leu). As predicted, the charged level of tRNA(GAG)(Leu) increased and the charged levels of four other Leu isoacceptors decreased or remained unchanged, but the charged level of tRNA(UAG)(Leu) increased unexpectedly. To remove this apparent inconsistency between theory and experiment we postulated a previously unknown common codon for tRNA(GAG)(Leu) and tRNA(UAG)(Leu). Subsequently, we demonstrated that the tRNA(GAG)(Leu) codon CUU is, in fact, read also by tRNA(UAG)(Leu), due to a uridine-5-oxyacetic acid modification.     

Modification of guanine to queuine in transfer RNAs during development and aging

The degree of modification of guanine to queuine in the four queuine-containing tRNAs (Q-tRNAs) has been studied from rats of various age groups, and bacterial cells in different growth phases by measuring the amount of G-tRNA present in these tRNA preparations by tRNA-guanine transferase. In very young (one-week old) animals, only a small amount of G to Q modification was observed. However, this modification was essentially complete in the tRNAs of nine-month old animals, thereafter, the amount of Q decreased steadily. Studies of tRNAs from leukemic lymphocytes and bacterial cells indicated that the degree of G to Q modification was related to the metabolic state of the cell. The possible role of the Q-deficient isoacceptors in translation control is discussed.     

Tissue-specific differences in human transfer RNA expression

Over 450 transfer RNA (tRNA) genes have been annotated in the human genome. Reliable quantitation of tRNA levels in human samples using microarray methods presents a technical challenge. We have developed a microarray method to quantify tRNAs based on a fluorescent dye-labeling technique. The first-generation tRNA microarray consists of 42 probes for nuclear encoded tRNAs and 21 probes for mitochondrial encoded tRNAs. These probes cover tRNAs for all 20 amino acids and 11 isoacceptor families. Using this array, we report that the amounts of tRNA within the total cellular RNA vary widely among eight different human tissues. The brain expresses higher overall levels of nuclear encoded tRNAs than every tissue examined but one and higher levels of mitochondrial encoded tRNAs than every tissue examined. We found tissue-specific differences in the expression of individual tRNA species, and tRNAs decoding amino acids with similar chemical properties exhibited coordinated expression in distinct tissue types. Relative tRNA abundance exhibits a statistically significant correlation to the codon usage of a collection of highly expressed, tissue-specific genes in a subset of tissues or tRNA isoacceptors. Our findings demonstrate the existence of tissue-specific expression of tRNA species that strongly implicates a role for tRNA heterogeneity in regulating translation and possibly additional processes in vertebrate organisms.     

tRNA over-expression in breast cancer and functional consequences

Increased proliferation and elevated levels of protein synthesis are characteristics of transformed and tumor cells. Though components of the translation machinery are often misregulated in cancers, what role tRNA plays in cancer cells has not been explored. We compare genome-wide tRNA expression in cancer-derived versus non-cancer-derived breast cell lines, as well as tRNA expression in breast tumors versus normal breast tissues. In cancer-derived versus non-cancer-derived cell lines, nuclear-encoded tRNAs increase by up to 3-fold and mitochondrial-encoded tRNAs increase by up to 5-fold. In tumors versus normal breast tissues, both nuclear- and mitochondrial-encoded tRNAs increase up to 10-fold. This tRNA over-expression is selective and coordinates with the properties of cognate amino acids. Nuclear- and mitochondrial-encoded tRNAs exhibit distinct expression patterns, indicating that tRNAs can be used as biomarkers for breast cancer. We also performed association analysis for codon usage-tRNA expression for the cell lines. tRNA isoacceptor expression levels are not geared towards optimal translation of house-keeping or cell line specific genes. Instead, tRNA isoacceptor expression levels may favor the translation of cancer-related genes having regulatory roles. Our results suggest a functional consequence of tRNA over-expression in tumor cells. tRNA isoacceptor over-expression may increase the translational efficiency of genes relevant to cancer development and progression.     

Selective charging of tRNA isoacceptors induced by amino-acid starvation

Aminoacylated (charged) transfer RNA isoacceptors read different messenger RNA codons for the same amino acid. The concentration of an isoacceptor and its charged fraction are principal determinants of the translation rate of its codons. A recent theoretical model predicts that amino-acid starvation results in 'selective charging' where the charging levels of some tRNA isoacceptors will be low and those of others will remain high. Here, we developed a microarray for the analysis of charged fractions of tRNAs and measured charging for all Escherichia coli tRNAs before and during leucine, threonine or arginine starvation. Before starvation, most tRNAs were fully charged. During starvation, the isoacceptors in the leucine, threonine or arginine families showed selective charging when cells were starved for their cognate amino acid, directly confirming the theoretical prediction. Codons read by isoacceptors that retain high charging can be used for efficient translation of genes that are essential during amino-acid starvation. Selective charging can explain anomalous patterns of codon usage in the genes for different families of proteins.     

Alterations in tRNA isoaccepting species during erythroid differentiation of the Friend leukemia cell.

The chromatographic profiles of isoaccepting tRNAs were analyzed at five time points during the 96 hr, dimethylsulfoxide induced, erythroid-like differentiation of Friend leukemia cells. Sixty-four isoaccepting species of tRNA for 16 amino acids were resolved by RPC-5 chromatography. The relative amounts of tRNAphe, tRNAile, and tRNAval species were maintained by the cells during differentiation; whereas the relative amounts of some of the isoacceptor tRNAs for the other 13 amino acids changed significantly. Fluctuations in amounts of isoacceptors occurred between 36 and 72 hr after addition of dimethysulfoxide, corresponding to globin mRNA appearance and hemoglobin synthesis, respectively. In most cases, thepredominant tRNA isoacceptors of uninduced cells were retained throughout differentiation. Notable exceptions were tRNA species for threonine, proline, and methionine. Some of the isoacceptors occurring in relatively smaller amounts were not expressed at all times. These changes possibly reflect the cell's functional adaptation of tRNA in differentiation for hemoglobin synthesis.     

Evolution of +1 Programmed Frameshifting Signals and Frameshift-Regulating tRNAs in the Order Saccharomycetales

Programmed translational frameshifting is a ubiquitous but rare mechanism of gene expression in which mRNA sequences cause the translational machinery to shift reading frames with extreme efficiency, up to at least 50%. The mRNA sequences responsible are deceptively simple; the sequence CUU-AGG-C causes about 40% frameshifting when inserted into an mRNA in the yeast Saccharomyces cerevisiae. The high efficiency of this site depends on a set of S. cerevisiae tRNA isoacceptors that perturb the mechanism of translation to cause the programmed translational error. The simplicity of the system might suggest that it could evolve frequently and perhaps be lost as easily. We have investigated the history of programmed +1 frameshifting in fungi. We find that frameshifting has persisted in two structural genes in budding yeasts, ABP140 and EST3 for about 150 million years. Further, the tRNAs that stimulate the event are equally old. Species that diverged from the lineage earlier both do not employ frameshifting and have a different complement of tRNAs predicted to be inimical to frameshifting. The stability of the coevolution of protein coding genes and tRNAs suggests that frameshifting has been selected for during the divergence of these species. [Reviewing Editor: Dr. Niles Lehman]     

Synonymous codon preferences in bacteriophage T4: a distinctive use of transfer RNAs from T4 and from its host Escherichia coli.

Codon usage data of bacteriophage T4 genes were compiled and synonymous codon preferences were investigated in comparison with tRNA availabilities in an infected cell. Since the genome of T4 is highly AT rich and its codon usage pattern is significantly different from that of its host Escherichia coli, certain codons of T4 genes need to be translated by appropriate host transfer RNAs present in minor amounts. To avoid this predicament, T4 phage seems to direct the synthesis of its own tRNA molecules and these phage tRNAs are suggested to supplement the host tRNA population with isoacceptors that are normally present in minor amounts. A positive correlation was found in that the frequency of E. coli optimal codons in T4 genes increases as the number of protein monomers per phage particle increases. A negative correlation was also found between the number of protein monomers per phage and the frequency of "T4 optimal codons", which are defined as those codons that are efficiently recognized by T4 tRNAs. From these observations it was proposed that tRNAs from the host are predominantly used for translation of highly expressed T4 genes while tRNAs from T4 tend to be used for translation of weakly expressed T4 genes. This distinctive tRNA-usage in T4 may be an optimization of translational efficiency, and an adjustment of T4-encoded tRNAs to the synonymous codon preferences, which are largely influenced by the high genomic AT-content, would have occurred during evolution.     

Changes in transfer ribonucleic acids of Bacillus subtilis during different growth phases

The transfer ribonucleic acids (tRNAs) of B. subtilis at different growth phases are examined for changes in the composition and the methylation of minor constituents. The composition of the tRNAs indicates about equal amounts of adenosine and uridine, and of guanosine and cytidine. About 3-4 residues are present as modified bases in the average tRNA molecule. The net composition of tRNAs appears to remain unaltered during different growth phases. In vitro methylation of tRNAs indicates lack of methyl groups in both exponentially growing cells and spores. In vivo methylation studies show tRNA methylation occurs during the stationary phase in the absence of net tRNA synthesis. Thus, both in vitro and in vivo methylation indicates that the tRNAs in exponentially growing cells do not contain their full complement of modified bases. More complete modification is noted in tRNAs from stationary cells or spores. Hence, tRNA modifications in general are preserved with fidelity even in the dormant spore but the possibility is left open that specific modifications of selected isoacceptors of tRNAs may occur.     

Scp160p is required for translational efficiency of codon-optimized mRNAs in yeast

The budding yeast multi-K homology domain RNA-binding protein Scp160p binds to >1000 messenger RNAs (mRNAs) and polyribosomes, and its mammalian homolog vigilin binds transfer RNAs (tRNAs) and translation elongation factor EF1alpha. Despite its implication in translation, studies on Scp160p''s molecular function are lacking to date. We applied translational profiling approaches and demonstrate that the association of a specific subset of mRNAs with ribosomes or heavy polysomes depends on Scp160p. Interaction of Scp160p with these mRNAs requires the conserved K homology domains 13 and 14. Transfer RNA pairing index analysis of Scp160p target mRNAs indicates a high degree of consecutive use of iso-decoding codons. As shown for one target mRNA encoding the glycoprotein Pry3p, Scp160p depletion results in translational downregulation but increased association with polysomes, suggesting that it is required for efficient translation elongation. Depletion of Scp160p also decreased the relative abundance of ribosome-associated tRNAs whose codons show low potential for autocorrelation on mRNAs. Conversely, tRNAs with highly autocorrelated codons in mRNAs are less impaired. Our data indicate that Scp160p might increase the efficiency of tRNA recharge, or prevent diffusion of discharged tRNAs, both of which were also proposed to be the likely basis for the translational fitness effect of tRNA pairing.     

Eukaryotic tRNA sequences present conserved and amino acid-specific structural signatures

Metazoan organisms have many tRNA genes responsible for decoding amino acids. The set of all tRNA genes can be grouped in sets of common amino acids and isoacceptor tRNAs that are aminoacylated by corresponding aminoacyl-tRNA synthetases. Analysis of tRNA alignments shows that, despite the high number of tRNA genes, specific tRNA sequence motifs are highly conserved across multicellular eukaryotes. The conservation often extends throughout the isoacceptors and isodecoders with, in some cases, two sets of conserved isodecoders. This study is focused on non-Watson–Crick base pairs in the helical stems, especially GoU pairs. Each of the four helical stems may contain one or more conserved GoU pairs. Some are amino acid specific and could represent identity elements for the cognate aminoacyl tRNA synthetases. Other GoU pairs are found in more than a single amino acid and could be critical for native folding of the tRNAs. Interestingly, some GoU pairs are anticodon-specific, and others are found in phylogenetically-specific clades. Although the distribution of conservation likely reflects a balance between accommodating isotype-specific functions as well as those shared by all tRNAs essential for ribosomal translation, such conservations may indicate the existence of specialized tRNAs for specific translation targets, cellular conditions, or alternative functions.     

Dihydrouridine-deficient tRNAs in Saccharomyces cerevisiae.

A mutation in Saccharomyces cerevisiae, designated mia, is responsible for the production of isoaccepting tRNA molecules with reduced extents of nucleoside modifications. The mia isoacceptors of tRNAPhe and one of the mutant isoacceptors of tRNATyr were highly purified for nucleoside composition analyses. The data indicate that the mutant isoacceptors are lacking some of the dihydrouridine moieties. This is consistent with our previous hypothesis that the mutant isoacceptors were accumulated due to a defect in a modification process [Lo, R.Y.C. and Bell, J.B. (1981) Current Genetics 3, 73-82). Data from in vitro poly-U translation experiments also support the previous results, suggesting in vivo biological activity of these mutant tRNAs.     

tRNA thiolation links translation to stress responses in Saccharomyces cerevisiae

Although tRNA modifications have been well catalogued, the precise functions of many modifications and their roles in mediating gene expression are still being elucidated. Whereas tRNA modifications were long assumed to be constitutive, it is now apparent that the modification status of tRNAs changes in response to different environmental conditions. The URM1 pathway is required for thiolation of the cytoplasmic tRNAs tGlu^UUC, tGln^UUG, and tLys^UUU in Saccharomyces cerevisiae. We demonstrate that URM1 pathway mutants have impaired translation, which results in increased basal activation of the Hsf1-mediated heat shock response; we also find that tRNA thiolation levels in wild-type cells decrease when cells are grown at elevated temperature. We show that defects in tRNA thiolation can be conditionally advantageous, conferring resistance to endoplasmic reticulum stress. URM1 pathway proteins are unstable and hence are more sensitive to changes in the translational capacity of cells, which is decreased in cells experiencing stresses. We propose a model in which a stress-induced decrease in translation results in decreased levels of URM1 pathway components, which results in decreased tRNA thiolation levels, which further serves to decrease translation. This mechanism ensures that tRNA thiolation and translation are tightly coupled and coregulated according to need.