The many applications of acid urea polyacrylamide gel electrophoresis to studies of tRNAs and aminoacyl-tRNA synthetases

Here we describe the many applications of acid urea polyacrylamide gel electrophoresis (acid urea PAGE) followed by Northern blot analysis to studies of tRNAs and aminoacyl-tRNA synthetases. Acid urea PAGE allows the electrophoretic separation of different forms of a tRNA, discriminated by changes in bulk, charge, and/or conformation that are brought about by aminoacylation, formylation, or modification of a tRNA. Among the examples described are (i) analysis of the effect of mutations in the Escherichia coli initiator tRNA on its aminoacylation and formylation; (ii) evidence of orthogonality of suppressor tRNAs in mammalian cells and yeast; (iii) analysis of aminoacylation specificity of an archaeal prolyl-tRNA synthetase that can aminoacylate archaeal tRNA(Pro) with cysteine, but does not aminoacylate archaeal tRNA(Cys) with cysteine; (iv) identification and characterization of the AUA-decoding minor tRNA(Ile) in archaea; and (v) evidence that the archaeal minor tRNA(Ile) contains a modified base in the wobble position different from lysidine found in the corresponding eubacterial tRNA.     

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.     

The complete set of tRNA species in Nanoarchaeum equitans

The archaeal parasite Nanoarchaeum equitans was found to generate five tRNA species via a unique process requiring the assembly of seperate 5' and 3' tRNA halves [Randau, L., Munch, R., Hohn, M.J., Jahn, D. and Soll, D. (2005) Nanoarchaeum equitans creates functional tRNAs from separate genes for their 5'- and 3'-halves. Nature 433, 537-541]. Biochemical evidence was missing for one of the computationally-predicted, joined tRNAs designated as tRNA(Trp). Our RT-PCR and sequencing results identify this tRNA as tRNA(Lys) (CUU) joined at the alternative position between bases 30 and 31. We show that the intron-containing tRNA(Trp) was misidentified in the initial Nanoarchaeum equitans genome annotation [E. Waters et al. (2003) The genome of Nanoarchaeum equitans: insights into early archaeal evolution and derived parasitism. Proc. Natl. Acad. Sci. USA 100, 12984-12988]. Along with a previously unidentified joined tRNA(Gln) (UUG), Nanoarchaeum equitans exhibits 44 tRNAs and is enabled to read all 61 sense codons. Features unique to this set of tRNA molecules are discussed.     

Sequence evidence in the archaeal genomes that tRNAs emerged through the combination of ancestral genes as 5' and 3' tRNA halves

The discovery of separate 5' and 3' halves of transfer RNA (tRNA) molecules-so-called split tRNA-in the archaeal parasite Nanoarchaeum equitans made us wonder whether ancestral tRNA was encoded on 1 or 2 genes. We performed a comprehensive phylogenetic analysis of tRNAs in 45 archaeal species to explore the relationship between the three types of tRNAs (nonintronic, intronic and split). We classified 1953 mature tRNA sequences into 22 clusters. All split tRNAs have shown phylogenetic relationships with other tRNAs possessing the same anticodon. We also mimicked split tRNA by artificially separating the tRNA sequences of 7 primitive archaeal species at the anticodon and analyzed the sequence similarity and diversity of the 5' and 3' tRNA halves. Network analysis revealed specific characteristics of and topological differences between the 5' and 3' tRNA halves: the 5' half sequences were categorized into 6 distinct groups with a sequence similarity of >80%, while the 3' half sequences were categorized into 9 groups with a higher sequence similarity of >88%, suggesting different evolutionary backgrounds of the 2 halves. Furthermore, the combinations of 5' and 3' halves corresponded with the variation of amino acids in the codon table. We found not only universally conserved combinations of 5'-3' tRNA halves in tRNA(iMet), tRNA(Thr), tRNA(Ile), tRNA(Gly), tRNA(Gln), tRNA(Glu), tRNA(Asp), tRNA(Lys), tRNA(Arg) and tRNA(Leu) but also phylum-specific combinations in tRNA(Pro), tRNA(Ala), and tRNA(Trp). Our results support the idea that tRNA emerged through the combination of separate genes and explain the sequence diversity that arose during archaeal tRNA evolution.     

Structural compensation in an archaeal selenocysteine transfer RNA.

A new type of structural compensation between the lengths of two perpendicularly oriented RNA double helices was found in the archaeal selenocysteine tRNA from Methanococcus jannascii. This tRNA contains only four base-pairs in the T-stem, one base-pair less than in all other cytosolic tRNAs. Our analysis shows that such a T-stem in an otherwise normal tRNA cannot guarantee the formation of the normal interactions between the D and T-loops. The absence of these interactions would affect the juxtaposition of the two tRNA helical domains, potentially damaging the tRNA function. In addition to the short T-stem, this tRNA possesses another unprecedented feature, a very long D-stem consisting of seven base-pairs. Taken as such, a seven base-pair D-stem will also disrupt the normal interaction between the D and T-loops. On the other hand, the presence of the universal nucleotides in both the D and T-loops suggests that these loops probably interact with each other in the same way as in other tRNAs. Here, we demonstrate that the short T-stem and the long D-stem can naturally compensate each other, thus providing the normal D/T interactions. Molecular modeling has helped suggest a detailed scheme of mutual compensation between these two unique structural aspects of the archaeal selenocysteine tRNA. In the light of this analysis, other structural and functional characteristics of the selenocysteine tRNAs are discussed.     

Distinct tRNA modifications in the thermo-acidophilic archaeon, Thermoplasma acidophilum

Thermoplasma acidophilum is a thermo-acidophilic archaeon. We purified tRNA^Leu (UAG) from T. acidophilum using a solid-phase DNA probe method and determined the RNA sequence after determining via nucleoside analysis and m⁷G-specific aniline cleavage because it has been reported that T. acidophilum tRNA contains m⁷G, which is generally not found in archaeal tRNAs. RNA sequencing and liquid chromatography–mass spectrometry revealed that the m⁷G modification exists at a novel position 49. Furthermore, we found several distinct modifications, which have not previously been found in archaeal tRNA, such as 4-thiouridine9, archaeosine13 and 5-carbamoylmethyuridine34. The related tRNA modification enzymes and their genes are discussed.     

Cytosolic yeast tRNA(His) is covalently modified when imported into mitochondria of Trypanosoma brucei.

The mitochondrial genome of Trypanosoma brucei does not encode any tRNAs. Instead, mitochondrial tRNAs are synthesized in the nucleus and subsequently imported into mitochondria. The great majority of mitochondrial tRNAs have cytosolic counterparts showing identical primary sequences. The only difference found between mitochondrial and cytosolic isotypes of the tRNAs are mitochondria-specific nucleotide modifications which appear to be a common feature of imported tRNAs in trypanosomes. In this study, a mutated yeast cytosolic tRNAHis was expressed in trypanosomes and its import phenotype was analyzed by cell fractionation and nuclease treatment of intact mitochondria. Furthermore, cytosolic and mitochondrial isotypes of the yeast tRNA(His) were specifically labeled and analyzed by limited alkaline hydrolysis. These experiments revealed the presence of mitochondria-specific nucleotide modifications in the yeast tRNA(His). The positions of the modifications were determined by direct enzymatic sequencing of the tRNA(His) and shown to correspond to the ultimate and penultimate nucleotides before the anticodon, the same relative positions which are modified in the mitochondrial isotype of trypanosomal tRNA(Tyr). The results demonstrate that covalent modification of tRNAs; in trypanosomal mitochondria can be used, in analogy to processing of precursor proteins during mitochondrial protein import, as a marker for import of both endogenous and heterologous tRNAs.     

Another cut for lysine tRNA: application of the hyperprocessing reaction reveals another stabilization strategy in metazoan lysine tRNAs

Recently, we revealed that the cloverleaf structure of some eukaryotic tRNAs is not always stable in vitro, and the denatured structures of these tRNAs are sometimes detected in bacterial RNase P reactions. We have designated the unusual internal cleavage reaction of these tRNAs as hyperprocessing. We have developed this hyperprocessing strategy as a useful tool for examining the stability of the tRNA cloverleaf structure. There are some common features in such unstable, hyperprocessible tRNAs, and the criteria for the hyperprocessing reaction of tRNA are extracted. Metazoan initiator methionine tRNAs and lysine tRNAs commonly fit the criteria, and are predicted to be hyperprocessible. The RNase P reactions of two metazoan lysine tRNAs from Homo sapiens and Caenorhabditis elegans, which fit the criteria, resulted in resistance to the internal cleavage reaction, while one bacterial lysine tRNA from Acholeplasma laidlawii, which also fits the criteria, was internally cleaved by the RNase P. The results showed that the metazoan lysine tRNAs examined are very stable without base modifications even under in vitro conditions. We also examined the 3'-half short construct of the human lysine tRNA, and the results showed that this RNA was internally cleaved by the enzyme. The results indicated that the human lysine tRNA has the ability to be hyperprocessed but is structurally stabilized in spite of lacking base modifications. A comparative study suggested, moreover, that the acceptor-stem bases should take part in the stabilization of metazoan lysine tRNAs. Our data strongly suggest that the cloverleaf shape of other metazoan lysine tRNAs should also be stabilized by means of similar strategies to in the case of human tRNA(Lys3).     

Modified view of tRNA: stability amid sequence diversity

In a recent issue of Molecular Cell, a report from the Phizicky lab challenges two commonly held notions regarding tRNAs: (1) tRNAs are stable molecules shielded from the machinery that cause other RNAs to have short half-lives, and (2) the many modifications found on tRNAs serve only subtle biological functions. Alexandrov et al. (2006) show that tRNA modifications can play an important role in tRNA stability, as without particular modifications, tRNA half-life can be reduced from hours to minutes by a new turnover pathway.     

28 The study of tRNA modifications by molecular dynamics

Modified nucleosides are prevalent in tRNA. Experimental studies reveal that modifications play an important role in tuning tRNA activity. In this study, molecular dynamics (MD) simulations were used to investigate how modifications alter tRNA structure and dynamics. The X-ray crystal structures of tRNA-Asp, tRNA-Phe, and tRNA-iMet, both with and without modifications, were used as initial structures for 333-ns time-scale MD trajectories with AMBER. For each tRNA molecule, three independent trajectory calculations were performed. Force field parameters were built using the RESP procedure of Cieplak et al. for 17 nonstandard tRNA residues. The global root-mean-square deviations (RMSDs) of atomic positions show that modifications only introduce significant rigidity to tRNA-Phe’s global structure. Interestingly, regional RMSDs of anticodon stem-loop suggest that modified tRNA has more rigid structure compared to the unmodified tRNA in this domain. The anticodon RMSDs of the modified tRNAs, however, are higher than those of corresponding unmodified tRNAs. These findings suggest that rigidity of the anticodon arm is essential for tRNA translocation in the ribosome complex, and, on the other hand, flexibility of anticodon might be critical for anticodon–codon recognition. We also measure the angle between the 3D L-shaped arms of tRNA; backbone atoms of acceptor stem and TψC stem loop are selected to indicate one vector, and backbone atoms of anticodon stem and D stem loop are selected to indicate the other vector. By measuring the angle between two vectors, we find that the initiator tRNA has a narrower range of hinge motion compared to tRNA-Asp and tRNA-Phe, which are elongator tRNA. This suggests that elongator tRNAs, which might require significant flexibility in this hinge to transition from the A–to-P site in the ribosome, have evolved to specifically accommodate this need.     

SPLITS: a new program for predicting split and intron-containing tRNA genes at the genome level

In the archaea, some tRNA precursors contain intron(s) not only in the anticodon loop region but also in diverse sites of the gene (intron-containing tRNA or cis-spliced tRNA). The parasite Nanoarchaeum equitans, a member of the Nanoarchaeota kingdom, creates functional tRNA from separate genes, one encoding the 5'-half and the other the 3'-half (split tRNA or trans-spliced tRNA). Although recent genome projects have revealed a huge amount of nucleotide sequence data in the archaea, a comprehensive methodology for intron-containing and split tRNA searching is yet to be established. We therefore developed SPLITS, which is aimed at searching for any type of tRNA gene and is especially focused on intron-containing tRNAs or split tRNAs at the genome level. SPLITS initially predicts the bulge-helix-bulge splicing motif (a well-known, required structure in archaeal pre-tRNA introns) to determine and remove the intronic regions of tRNA genes. The intron-removed DNA sequences are automatically queried to tRNAscan-SE. SPLITS can predict known tRNAs with single introns located at unconventional sites on the genes (100%), tRNAs with double introns (85.7%), and known split tRNAs (100%). Our program will be very useful for identifying novel tRNA genes after completion of genome projects. The SPLITS source code is freely downloadable at http://splits.iab.keio.ac.jp/.     

Diverse cell stresses induce unique patterns of tRNA up- and down-regulation: tRNA-seq for quantifying changes in tRNA copy number

Emerging evidence points to roles for tRNA modifications and tRNA abundance in cellular stress responses. While isolated instances of stress-induced tRNA degradation have been reported, we sought to assess the effects of stress on tRNA levels at a systems level. To this end, we developed a next-generation sequencing method that exploits the paucity of ribonucleoside modifications at the 3′-end of tRNAs to quantify changes in all cellular tRNA molecules. Application of this tRNA-seq method to Saccharomyces cerevisiae identified all 76 expressed unique tRNA species out of 295 coded in the yeast genome, including all isoacceptor variants, with highly precise relative (fold-change) quantification of tRNAs. In studies of stress-induced changes in tRNA levels, we found that oxidation (H₂O₂) and alkylation (methylmethane sulfonate, MMS) stresses induced nearly identical patterns of up- and down-regulation for 58 tRNAs. However, 18 tRNAs showed opposing changes for the stresses, which parallels our observation of signature reprogramming of tRNA modifications caused by H₂O₂ and MMS. Further, stress-induced degradation was limited to only a small proportion of a few tRNA species. With tRNA-seq applicable to any organism, these results suggest that translational control of stress response involves a contribution from tRNA abundance.     

Archaeal Pus10 proteins can produce both pseudouridine 54 and 55 in tRNA

Pus10, a recently identified pseudouridine (Ψ) synthase, does not belong to any of the five commonly identified families of Ψ synthases. Pyrococcus furiosus Pus10 has been shown to produce Ψ55 in tRNAs. However, in vitro studies have identified another mechanism for tRNA Ψ55 production in Archaea, which uses Cbf5 and other core proteins of the H/ACA ribonucleoprotein complex, in a guide RNA-independent manner. Pus10 homologs have been observed in nearly all sequenced archaeal genomes and in some higher eukaryotes, but not in yeast and bacteria. This coincides with the presence of Ψ54 in the tRNAs of Archaea and higher eukaryotes and its absence in yeast and bacteria. No tRNA Ψ54 synthase has been reported so far. Here, using recombinant Methanocaldococcus jannaschii and P. furiosus Pus10, we show that these proteins can function as synthase for both tRNA Ψ54 and Ψ55. The two modifications seem to occur independently. Salt concentration dependent variations in these activities of both proteins are observed. The Ψ54 synthase activity of M. jannaschii protein is robust, while the same activity of P. furiosus protein is weak. Probable reasons for these differences are discussed. Furthermore, unlike bacterial TruB and yeast Pus4, archaeal Pus10 does not require a U54•A58 reverse Hoogstein base pair and pyrimidine at position 56 to convert tRNA U55 to Ψ55. The homology of eukaryal Pus10 with archaeal Pus10 suggests that the former may also have a tRNA Ψ54 synthase activity.     

Identifying the ligated amino acid of archaeal tRNAs based on positions outside the anticodon

Proper recognition of tRNAs by their aminoacyl-tRNA synthetase is essential for translation accuracy. Following evidence that the enzymes can recognize the correct tRNA even when anticodon information is masked, we search for additional nucleotide positions within the tRNA molecule that potentially contain information for amino acid identification. Analyzing 3936 sequences of tRNA genes from 86 archaeal species, we show that the tRNAs’ cognate amino acids can be identified by the information embedded in the tRNAs’ nucleotide positions without relying on the anticodon information. We present a small set of six to 10 informative positions along the tRNA, which allow for amino acid identification accuracy of 90.6% to 97.4%, respectively. We inspected tRNAs for each of the 20 amino acid types for such informative positions and found that tRNA genes for some amino acids are distinguishable from others by as few as one or two positions. The informative nucleotide positions are in agreement with nucleotide positions that were experimentally shown to affect the loaded amino acid identity. Interestingly, the knowledge gained from the tRNA genes of one archaeal phylum does not extrapolate well to another phylum. Furthermore, each species has a unique ensemble of nucleotides in the informative tRNA positions, and the similarity between the sets of positions of two distinct species reflects their evolutionary distance. Hence, we term this set of informative positions a “tRNA cipher.” It is tempting to suggest that the diverging code identified here might also serve the aminoacyl tRNA synthetase in the task of tRNA recognition.