Analysis of genomic tRNA sets from Bacteria, Archaea, and Eukarya points to anticodon-codon hydrogen bonds as a major determinant of tRNA compositional variations

Analysis of 100 complete sets of the cytoplasmic elongator tRNA genes from Bacteria, Archaea, and Eukarya pointed to correspondences between types of anticodon and composition of the rest of the tRNA body. The number of the hydrogen bonds formed between the complementary nucleotides in the anticodon-codon duplex appeared as a major quantitative parameter determining covariations in all three domains of life. Our analysis has supported and advanced the "extended anticodon" concept that is based on the argument that the decoding performance of the anticodon is enhanced by selection of a matching anticodon stem-loop sequence, as reported by Yarus in 1982. In addition to the anticodon stem-loop, we have found covariations between the anticodon nucleotides and the composition of the distant regions of their respective tRNAs that include dihydrouridine (D) and thymidyl (T) stem-loops. The majority of the covariable tRNA positions were found at the regions with the increased dynamic potential--such as stem-loop and stem-stem junctions. The consistent occurrences of the covariations on the multigenomic level suggest that the number and pattern of the hydrogen bonds in the anticodon-codon duplex constitute a major factor in the course of translation that is reflected in the fine-tuning of the tRNA composition and structure.     

Conflict between translation initiation and elongation in vertebrate mitochondrial genomes

The strand-biased mutation spectrum in vertebrate mitochondrial genomes results in an AC-rich L-strand and a GT-rich H-strand. Because the L-strand is the sense strand of 12 protein-coding genes out of the 13, the third codon position is overall strongly AC-biased. The wobble site of the anticodon of the 22 mitochondrial tRNAs is either U or G to pair with the most abundant synonymous codon, with only one exception. The wobble site of Met-tRNA is C instead of U, forming the Watson-Crick match with AUG instead of AUA, the latter being much more frequent than the former. This has been attributed to a compromise between translation initiation and elongation; i.e., AUG is not only a methionine codon, but also an initiation codon, and an anticodon matching AUG will increase the initiation rate. However, such an anticodon would impose selection against the use of AUA codons because AUA needs to be wobble-translated. According to this translation conflict hypothesis, AUA should be used relatively less frequently compared to UUA in the UUR codon family. A comprehensive analysis of mitochondrial genomes from a variety of vertebrate species revealed a general deficiency of AUA codons relative to UUA codons. In contrast, urochordate mitochondrial genomes with two tRNA(Met) genes with CAU and UAU anticodons exhibit increased AUA codon usage. Furthermore, six bivalve mitochondrial genomes with both of their tRNA-Met genes with a CAU anticodon have reduced AUA usage relative to three other bivalve mitochondrial genomes with one of their two tRNA-Met genes having a CAU anticodon and the other having a UAU anticodon. We conclude that the translation conflict hypothesis is empirically supported, and our results highlight the fine details of selection in shaping molecular evolution.     

Identification and characterization of a tRNA decoding the rare AUA codon in Haloarcula marismortui

Annotation of the complete genome of the extreme halophilic archaeon Haloarcula marismortui does not include a tRNA for translation of AUA, the rare codon for isoleucine. This is a situation typical for most archaeal genomes sequenced to date. Based on computational analysis, it has been proposed recently that a single intron-containing tRNA gene produces two very similar but functionally different tRNAs by means of alternative splicing; a UGG-decoding tRNA(TrpCCA) and an AUA-decoding tRNA(IleUAU). Through analysis of tRNAs from H. marismortui, we have confirmed the presence of tRNA(TrpCCA), but found no evidence for the presence of tRNA(IleUAU). Instead, we have shown that a tRNA, currently annotated as elongator methionine tRNA and containing CAU as the anticodon, is aminoacylated with isoleucine in vivo and that this tRNA represents the missing isoleucine tRNA. Interestingly, this tRNA carries a base modification of C34 in the anticodon different from the well-known lysidine found in eubacteria, which switches the amino acid identity of the tRNA from methionine to isoleucine and its decoding specificity from AUG to AUA. The methods described in this work for the identification of individual tRNAs present in H. marismortui provide the tools necessary for experimentally confirming the presence of any tRNA in a cell and, thereby, to test computational predictions of tRNA genes.     

A primordial tRNA modification required for the evolution of life?

The evolution of reading frame maintenance must have been an early event, and presumably preceded the emergence of the three domains Archaea, Bacteria and Eukarya. Features evolved early in reading frame maintenance may still exist in present-day organisms. We show that one such feature may be the modified nucleoside 1-methylguanosine (m(1)G37), which prevents frameshifting and is present adjacent to and 3' of the anticodon (position 37) in the same subset of tRNAs from all organisms, including that with the smallest sequenced genome (Mycoplasma genitalium), and organelles. We have identified the genes encoding the enzyme tRNA(m(1)G37)methyltransferase from all three domains. We also show that they are orthologues, and suggest that they originated from a primordial gene. Lack of m(1)G37 severely impairs the growth of a bacterium and a eukaryote to a similar degree. Yeast tRNA(m(1)G37)methyltransferase also synthesizes 1-methylinosine and participates in the formation of the Y-base (yW). Our results suggest that m(1)G37 existed in tRNA before the divergence of the three domains, and that a tRNA(m(1)G37)methyltrans ferase is part of the minimal set of gene products required for life.     

Identification and codon reading properties of 5-cyanomethyl uridine,a new modified nucleoside found in the anticodon wobble position of mutant haloarchaeal isoleucine tRNAs

Most archaea and bacteria use a modified C in the anticodon wobble position of isoleucine tRNA to base pair with A but not with G of the mRNA. This allows the tRNA to read the isoleucine codon AUA without also reading the methionine codon AUG. To understand why a modified C, and not U or modified U, is used to base pair with A, we mutated the C34 in the anticodon of Haloarcula marismortui isoleucine tRNA (tRNA₂^Ile) to U, expressed the mutant tRNA in Haloferax volcanii, and purified and analyzed the tRNA. Ribosome binding experiments show that although the wild-type tRNA₂^Ile binds exclusively to the isoleucine codon AUA, the mutant tRNA binds not only to AUA but also to AUU, another isoleucine codon, and to AUG, a methionine codon. The G34 to U mutant in the anticodon of another H. marismortui isoleucine tRNA species showed similar codon binding properties. Binding of the mutant tRNA to AUG could lead to misreading of the AUG codon and insertion of isoleucine in place of methionine. This result would explain why most archaea and bacteria do not normally use U or a modified U in the anticodon wobble position of isoleucine tRNA for reading the codon AUA. Biochemical and mass spectrometric analyses of the mutant tRNAs have led to the discovery of a new modified nucleoside, 5-cyanomethyl U in the anticodon wobble position of the mutant tRNAs. 5-Cyanomethyl U is present in total tRNAs from euryarchaea but not in crenarchaea, eubacteria, or eukaryotes.     

The in vivo stability, maturation and aminoacylation of anticodon-substituted Escherichia coli initiator methionine tRNAs

We have constructed eight anticodon-modified Escherichia coli initiator methionine (fMet) tRNAs by insertion of synthetic ribotrinucleotides between two fragments ('half molecules') derived from the initiator tRNA. The trinucleotides, namely CAU (the normal anticodon), CAA, CAC, CAG, GAA, GAC, GAG and GAU, were joined to the 5' and 3' tRNA fragments with T4 RNA ligase. The strategy of reconstruction permitted the insertion of radioactive 32P label between nucleotides 36 and 37. tRNAs were microinjected into the cytoplasm of Xenopus laevis oocytes, and the following properties were evaluated: the stability of these eubacterial tRNA variants in the eukaryotic oocytes; the enzymatic modification of the adenosine at position 37 (3' adjacent to the anticodon) and aminoacylation of the chimeric tRNAs by endogenous oocyte aminoacyl-tRNA synthetases. In contrast to other variants, the two RNAs having CAU and GAU anticodons were stable and underwent quantitative modification at A-37. These results show that the enzyme responsible for the modification of A-37 to N-[N-(9-beta-D-ribofuranosylpurine-6-yl)carbamoyl]threonine (t6A) is present in the cytoplasm of oocytes and is very sensitive to the anticodon environment of the tRNA. Also, these same GAU and CAU anticodon-containing tRNAs are fully aminoacylated with the heterologous oocyte aminoacyl-tRNA synthetases in vivo. During the course of this work we developed a generally applicable assay for the aminoacylation of femtomole amounts of labelled tRNAs.     

Suppression of amber codons in vivo as evidence that mutants derived from Escherichia coli initiator tRNA can act at the step of elongation in protein synthesis

The absence of a Watson-Crick base pair at the end of the amino acid acceptor stem is one of the features which distinguishes prokaryotic initiator tRNAs as a class from all other tRNAs. We show that this structural feature prevents Escherichia coli initiator tRNA from acting as an elongator in protein synthesis in vivo. We generated a mutant of E. coli initiator tRNA in which the anticodon sequence is changed from CAU to CUA (the T35A36 mutant). This mutant tRNA has the potential to read the amber termination codon UAG. We then coupled this mutation to others which change the C1.A72 mismatch at the end of the acceptor stem to either a U1:A72 base pair (T1 mutant) or a C1:G72 base pair (G72 mutant). Transformation of E. coli CA274 (HfrC Su- lacZ125am trpEam) with multicopy plasmids carrying the mutant initiator tRNA genes show that mutant tRNAs carrying changes in both the anticodon sequence and the acceptor stem suppress amber codons in vivo, whereas mutant tRNA with changes in the anticodon sequence alone does not. Mutant tRNAs with the above anticodon sequence change are aminoacylated with glutamine in vitro. Measurement of kinetic parameters for aminoacylation by E. coli glutaminyl-tRNA synthetase show that both the nature of the base pair at the end of the acceptor stem and the presence or absence of a base pair at this position can affect aminoacylation kinetics. We discuss the implications of this result on recognition of tRNAs by E. coli glutaminyl-tRNA synthetase.     

Evolution of the mitochondrial genetic code II. Reassignment of codon AUA from isoleucine to methionine

Summary The reassignment of codon AUA from isoleucine to methionine during mitochondrial evolution may be explained by the codon reassignment (capture) hypothesis without assuming direct replacement of isoleucine by methionine in mitochondrial proteins. According to this hypothesis, codon AUA would have disappeared from the reading frames of messenger RNA. AUA codons would have mutated mainly to AUU isoleucine codons because of constraints resulting from elimination of tRNA Ile with anticodon *CAU (in which *C is lysidine). Later, tRNA Met (CAU) would have undergone structural changes enabling it to pair with both AUG and AUA. AUA codons, formed by mutations of other codons, including AUG, would have reappeared and would have been translated as methionine.     

Initiator-Elongator Discrimination in Vertebrate tRNAs for Protein Synthesis

Initiator tRNAs are used exclusively for initiation of protein synthesis and not for the elongation step. We show, in vivo and in vitro, that the primary sequence feature that prevents the human initiator tRNA from acting in the elongation step is the nature of base pairs 50:64 and 51:63 in the TΨC stem of the initiator tRNA. Various considerations suggest that this is due to sequence-dependent perturbation of the sugar phosphate backbone in the TΨC stem of initiator tRNA, which most likely blocks binding of the elongation factor to the tRNA. Because the sequences of all vertebrate initiator tRNAs are identical, our findings with the human initiator tRNA are likely to be valid for all vertebrate systems. We have developed reporter systems that can be used to monitor, in mammalian cells, the activity in elongation of mutant human initiator tRNAs carrying anticodon sequence mutations from CAU to CCU (the C35 mutant) or to CUA (the U35A36 mutant). Combination of the anticodon sequence mutation with mutations in base pairs 50:64 and 51:63 yielded tRNAs that act as elongators in mammalian cells. Further mutation of the A1:U72 base pair, which is conserved in virtually all eukaryotic initiator tRNAs, to G1:C72 in the C35 mutant background yielded tRNAs that were even more active in elongation. In addition, in a rabbit reticulocyte in vitro protein-synthesizing system, a tRNA carrying the TΨC stem and the A1:U72-to-G1:C72 mutations was almost as active in elongation as the elongator methionine tRNA. The combination of mutant initiator tRNA with the CCU anticodon and the reporter system developed here provides the first example of missense suppression in mammalian cells.     

Role of the three consecutive G:C base pairs conserved in the anticodon stem of initiator tRNAs in initiation of protein synthesis in Escherichia coli.

The three consecutive G:C base pairs, G29:C41, G30:C40, and G31:C39, are conserved in the anticodon stem of virtually all initiator tRNAs from eubacteria, eukaryotes, and archaebacteria. We show that these G:C base pairs are important for function of the tRNA in initiation of protein synthesis in vivo. We changed these base pairs individually and in combinations and analyzed the activities of the mutant Escherichia coli initiator tRNAs in initiation in vivo. For assessment of activity of the mutant tRNAs in vivo, mutations in the G:C base pairs were coupled to mutation in the anticodon sequence from CAU to CUA. Mutations in each of the G:C base pairs reduced activity of the mutant tRNA in initiation, with mutation in the second G:C base pair having the most severe effect. The greatly reduced activity of this C30:G40 mutant tRNA is not due to defects in aminoacylation or formulation of the tRNA or defects in base modification of the A37, next to the anticodon, which we had previously shown to be important for activity of the mutant tRNAs in initiation. The anticodon stem mutants are most likely affected specifically at the step of binding to the ribosomal P site. The pattern of cleavages in the anticodon loop of mutant tRNAs by S1 nuclease indicate that the G:C base pairs may be involved directly in interactions of the tRNA with components of the P site on the ribosome rather than indirectly by inducing a particular conformation of the anticodon loop critical for function of the tRNA in initiation.     

Nucleotide sequences of animal mitochondrial tRNAs(Met) possibly recognizing both AUG and AUA codons.

To elucidate the role of modified nucleosides of tRNA in mitochondrial translation systems, especially with regard to their codon recognition, we purified mitochondrial tRNAs(Met) isolated from liver of frog, chicken and rat, and determined their nucleotide sequences. All of these tRNAs(Met) were found to possess 5-formylcytidine in the first letter of the anticodon, which is known to be prerequisite for bovine mt tRNA(Met) to decode AUA codon as well as AUG codon. These tRNA possesses two pseudeuridines in similar positions, and only chicken tRNA(Met) had ribothymidine at the first position of the T-loop, which is always found in the usual tRNAs. Considering that AUA codon is used as five times frequently as AUG codon in these animal mitochondrial genomes, it is deduced that 5-formylcytidine at the wobble position is essential for the recognition of both AUA and AUG codons.     

N2-methylation of guanosine at position 10 in tRNA is catalyzed by a THUMP domain-containing, S-adenosylmethionine-dependent methyltransferase, conserved in Archaea and Eukaryota

In sequenced genomes, genes belonging to the cluster of orthologous group COG1041 are exclusively, and almost ubiquitously, found in Eukaryota and Archaea but never in Bacteria. The corresponding gene products exhibit a characteristic Rossmann fold, S-adenosylmethionine-dependent methyltransferase domain in the C terminus and a predicted RNA-binding THUMP (thiouridine synthases, RNA methyltransferases, and pseudouridine synthases) domain in the N terminus. Recombinant PAB1283 protein from the archaeon Pyrococcus abyssi GE5, a member of COG1041, was purified and shown to behave as a monomeric 39-kDa entity. This protein (EC 2.1.1.32), now renamed (Pab)Trm-G10, which is extremely thermostable, forms a 1:1 complex with tRNA and catalyzes the adenosylmethionine-dependent methylation of the exocyclic amino group (N(2)) of guanosine located at position 10. Depending on the experimental conditions used, as well as the tRNA substrate tested, the enzymatic reaction leads to the formation of either N(2)-monomethyl (m(2)G) or N(2)-dimethylguanosine (m(2)(2)G). Interestingly, (Pab)Trm-G10 exhibits different domain organization and different catalytic site architecture from another, earlier characterized, tRNA-dimethyltransferase from Pyrococcus furiosus ((Pfu)Trm-G26, also known as (Pfu)Trm1, a member of COG1867) that catalyzes an identical two-step dimethylation of guanosine but at position 26 in tRNAs and is also conserved among all sequenced Eukaryota and Archaea. The co-occurrence of these two guanosine dimethyltransferases in both Archaea and Eukaryota but not in Bacteria is a hallmark of distinct tRNAs maturation strategies between these domains of life.     

TFAM detects co-evolution of tRNA identity rules with lateral transfer of histidyl-tRNA synthetase

We present TFAM, an automated, statistical method to classify the identity of tRNAs. TFAM, currently optimized for bacteria, classifies initiator tRNAs and predicts the charging identity of both typical and atypical tRNAs such as suppressors with high confidence. We show statistical evidence for extensive variation in tRNA identity determinants among bacterial genomes due to variation in overall tDNA base content. With TFAM we have detected the first case of eukaryotic-like tRNA identity rules in bacteria. An alpha-proteobacterial clade encompassing Rhizobiales, Caulobacter crescentus and Silicibacter pomeroyi, unlike a sister clade containing the Rickettsiales, Zymomonas mobilis and Gluconobacter oxydans, uses the eukaryotic identity element A73 instead of the highly conserved prokaryotic element C73. We confirm divergence of bacterial histidylation rules by demonstrating perfect covariation of alpha-proteobacterial tRNA(His) acceptor stems and residues in the motif IIb tRNA-binding pocket of their histidyl-tRNA synthetases (HisRS). Phylogenomic analysis supports lateral transfer of a eukaryotic-like HisRS into the alpha-proteobacteria followed by in situ adaptation of the bacterial tDNA(His) and identity rule divergence. Our results demonstrate that TFAM is an effective tool for the bioinformatics, comparative genomics and evolutionary study of tRNA identity.     

Ribosome-initiator tRNA complex as an intermediate in translation initiation in Escherichia coli revealed by use of mutant initiator tRNAs and specialized ribosomes.

For functional studies of mutant Escherichia coli initiator tRNAs in vivo, we previously described a strategy based on the use of tRNA genes carrying an anticodon sequence change from CAU to CUA along with a mutant chloramphenicol acetyltransferase (CAT) gene carrying an initiation codon change from AUG to UAG. Surprisingly, under conditions where the mutant initiator tRNA is optimally active, the CAT gene with the UAG initiation codon produced more CAT protein (3- to 9-fold more depending on the conditions) than the wild-type CAT gene. Here we show that two new mutant CAT genes having GUC and AUC initiation codons also produce more of the CAT protein in the presence of the corresponding mutant initiator tRNAs. These results are most easily understood if assembly of the 30S ribosome-initiator tRNA-mRNA initiation complex in vivo proceeds with the 30S ribosome binding first to the initiator tRNA and then to the mRNA. In cells overproducing the mutant initiator tRNAs, most ribosomes would carry the mutant initiator tRNA and these ribosomes would select the mutant CAT mRNA over the other mRNAs.     

tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes

tRNAscan-SE has been widely used for transfer RNA (tRNA) gene prediction for over twenty years, developed just as the first genomes were decoded. With the massive increase in quantity and phylogenetic diversity of genomes, the accurate detection and functional prediction of tRNAs has become more challenging. Utilizing a vastly larger training set, we created nearly one hundred specialized isotype- and clade-specific models, greatly improving tRNAscan-SE’s ability to identify and classify both typical and atypical tRNAs. We employ a new comparative multi-model strategy where predicted tRNAs are scored against a full set of isotype-specific covariance models, allowing functional prediction based on both the anticodon and the highest-scoring isotype model. Comparative model scoring has also enhanced the program''s ability to detect tRNA-derived SINEs and other likely pseudogenes. For the first time, tRNAscan-SE also includes fast and highly accurate detection of mitochondrial tRNAs using newly developed models. Overall, tRNA detection sensitivity and specificity is improved for all isotypes, particularly those utilizing specialized models for selenocysteine and the three subtypes of tRNA genes encoding a CAU anticodon. These enhancements will provide researchers with more accurate and detailed tRNA annotation for a wider variety of tRNAs, and may direct attention to tRNAs with novel traits.     

The anticodon and discriminator base are major determinants of cysteine tRNA identity in vivo.

Mutants of the Escherichia coli initiator tRNA (tRNA(fMet)) have been used to examine the role of the anticodon and discriminator base in in vivo aminoacylation of tRNAs by cysteinyl-tRNA synthetase. Substitution of the methionine anticodon CAU with the cysteine anticodon GCA was found to allow initiation of protein synthesis by the mutant tRNA from a complementary initiation codon in a reporter protein. Sequencing of the protein revealed that cysteine comprised about half of the amino acid at the N terminus. An additional mutation, converting the discriminator base of tRNA(GCAfMet) from A73 to the base present in tRNA(Cys) (U73), resulted in a 6-fold increase in the amount of protein produced and insertion of greater than or equal to 90% cysteine in response to the complementary initiation codon. Substitution of C73 or G73 at the discriminator position led to insertion of little or no cysteine, indicating the importance of U73 for recognition of the tRNA by cysteinyl-tRNA synthetase. Single base changes in the anticodon of tRNA(GCAfMet) containing U73 from GCA to UCA, GUA, GCC, and GCG (changes underlined) eliminated or dramatically reduced cysteine insertion by the mutant initiator tRNA indicating that all three cysteine anticodon bases are essential for specific aminoacylation of the tRNA with cysteine in vivo.