Diversity and similarity in the tRNA world: Overall view and case study on malaria-related tRNAs

Transfer RNAs (tRNAs) are ancient macromolecules that have evolved under various environmental pressures as adaptors in translation in all forms of life but also towards alternative structures and functions. The present knowledge on both “canonical” and “deviating” signature motifs retrieved from vertical and horizontal sequence comparisons is briefly reviewed. Novel characteristics, proper to tRNAs from a given translation system, are revealed by a case study on the nuclear and organellar tRNA sets from malaria-related organisms. Unprecedented distinctive features for Plasmodium falciparum apicoplastic tRNAs appear, which provide novel routes to be explored towards anti-malarial drugs. The ongoing high-throughput sequencing programs are expected to allow for further horizontal comparisons and to reveal other signatures of either full or restricted sets of tRNAs.     

A comparison of mitochondrial tRNAs in five vertebrates

Sequences of several vertebrate mitochondrial tRNAs were aligned and compared. The comparisons were made in pairs of the tRNAs for an identical amino acid. There are 22 genes for different tRNAs in each vertebrate mitochondrial DNA. The closest similarities were between rat and mouse, the next were between mammals, and the widest difference was between human or rat and Xenopus laevis. However, there were very wide variations between different amino acids in each set of comparisons. The time lapse for each percent of difference greatly increased with evolutionary separation. Most of the nucleotide substitutions appeared to be neutral in character.Based on a presentation made at a workshop- Aminoacyl-tRNA Synthetases and the Evolution of the Genetic Code-held at Berkeley, CA, July 17–20, 1994     

Characterization of Small RNAs Derived from tRNAs,rRNAs and snoRNAs and Their Response to Heat Stress in Wheat Seedlings

Small RNAs (sRNAs) derived from non-coding RNAs (ncRNAs), such as tRNAs, rRNAs and snoRNAs, have been identified in various organisms. Several observations have indicated that cleavage of tRNAs and rRNAs is induced by various stresses. To clarify whether sRNAs in wheat derived from tRNAs (stRNAs), rRNAs (srRNAs) and snoRNAs (sdRNAs) are produced specifically in association with heat stress responses, we carried out a bioinformatic analysis of sRNA libraries from wheat seedlings and performed comparisons between control and high-temperature-treated samples to measure the differential abundance of stRNAs, srRNAs and sdRNAs. We found that the production of sRNAs from tRNAs, 5.8S rRNAs, and 28S rRNAs was more specific than that from 5S rRNAs and 18S rRNAs, and more than 95% of the stRNAs were processed asymmetrically from the 3’ or 5’ ends of mature tRNAs. We identified 333 stRNAs and 8,822 srRNAs that were responsive to heat stress. Moreover, the expression of stRNAs derived from tRNA-Val-CAC, tRNA-Thr-UGU, tRNA-Tyr-GUA and tRNA-Ser-UGA was not only up-regulated under heat stress but also induced by osmotic stress, suggesting that the increased cleavage of tRNAs might be a mechanism that developed in wheat seedlings to help them cope with adverse environmental conditions.     

The Presence in tRNA Molecule Sequences of the Double Hairpin,an Evolutionary Stage Through Which the Origin of this Molecule is Thought to have Passed

We analyse 6,810 tRNAs, calculating the free energy of the corresponding double hairpin and ‘cigar’ secondary structures, for which we find a high thermodynamic and statistical significance. We also analyse these tRNAs for similarity and complementarity of their 5′ and 3′ halves or segments of them in intra-and inter-molecular comparisons. We find very clear signs that the two halves of tRNAs had an evident evolutionary relationship, although it is not totally clear whether this was a relationship of homology or complementarity between the 5′ and 3′ halves of tRNAs, even if there is strong evidence in favour of the homology hypothesis. Overall, these data favour models for the origin of the tRNA molecule postulating that a duplication event involving a hairpin structure as a precursor was involved in the origin of this molecule. Moreover, we interpret these results and favour the hypothesis that sees the assembly of two hairpin structures sharing a homology relationship as the intermediate evolutionary stage preceding the appearance of the cloverleaf structure of tRNA.     

Evolution of transfer RNA

Evolution by gene duplication and subsequent divergence is indicated by similarities common to 43 different transfer RNAs. Pairwise comparisons of these tRNAs reveal additional similarity, greatest for certain pairs of tRNAs for the same amino acid in the same organism, and also occurring in certain pairs of tRNAs for different amino acids in the same organism. Although tRNAs functionally interact with several other molecules, there have been surprisingly few restrictions on the divergence of their primary structures. This divergence has proceeded so far that clear phylogenetic separations are absent in most cases: it it impossible to construct a coherent phylogeny for most of the 43. Selection and stochastic processes have both been active in the evolution of tRNA. Selection has favored moderate change more than expected and has reduced radical change below that expected from stochastic processes alone. Two obvious effects of selection are nine invariant loci, another five that are always purines and five others that are always pyrimidines, in the tRNAs involved in protein synthesis. In addition to these constraints in the primary nucleotide sequence, the method of “identical site equivalents”, introduced here, demonstrates that further constraints exist equivalent to about 12 additional invariant loci. These “invisible” restraints reflect disperse chemical forces maintaining the tertiary structure and reducing evolutionary divergence to an extent quantitatively comparable to that of the nine observable invariant loci. The average divergence (49·4%) for pairs of tRNAs for different amino acids involved in protein synthesis represents an equilibrium between natural selection and stochastic processes. These tRNAs have had time to diverge nearly to the 75% maximum expected from stochastic process alone; this is shown by comparing the two glycine tRNAs involved in peptidoglycan synthesis with tRNAs for different amino acids participating in polypeptide synthesis. The rates of nucleotide replacements in genes coding for the tRNAs and the cytochromes c are about the same: 2 × 10 ^?10 replacements per nucleotide site per year.     

Nucleotide sequences of the mitochondrial genes trnS(TGA) encoding tRNA(TGASer) in Oenothera berteriana and Arabidopsis thaliana

The genes encoding tRNA(TGASer) have been investigated in the mitochondrial (mt) genomes of Oenothera berteriana and Arabidopsis thaliana. Sequence analysis shows four nucleotide (nt) differences between the two dicots, but only two differences between each dicot and the available monocot sequences. Similarity comparisons identify these genes as encoding a native mt tRNA(TGASer), with less than 77% of the nt identical to the corresponding chloroplast tRNAs.     

Experimental RNomics and genomic comparative analysis reveal a large group of species-specific small non-message RNAs in the silkworm Bombyx mori

Accumulating evidences show that small non-protein coding RNAs (ncRNAs) play important roles in development, stress response and other cellular processes. The silkworm is an important model for studies on insect genetics and control of lepidopterous pests. Here, we have performed the first systematic identification and analysis of intermediate size ncRNAs (50-500 nt) in the silkworm. We identified 189 novel ncRNAs, including 141 snoRNAs, six snRNAs, three tRNAs, one SRP and 38 unclassified ncRNAs. Forty ncRNAs showed significantly altered expression during silkworm development or across specific stage transitions. Genomic comparisons revealed that 123 of these ncRNAs are potentially silkworm-specific. Analysis of the genomic organization of the ncRNA loci showed that 32.62% of the novel snoRNA loci are intergenic, and that all the intronic snoRNAs follow the pattern of one-snoRNA-per-intron. Target site analysis predicted a total of 95 2'-O-methylation and pseudouridylation modification sites of rRNAs, snRNAs and tRNAs. Together, these findings provide new clues for future functional study of ncRNA during insect development and evolution.     

The role of tRNA and ribosome competition in coupling the expression of different mRNAs in Saccharomyces cerevisiae

Protein synthesis translates information from messenger RNAs into functional proteomes. Because of the finite nature of the resources required by the translational machinery, both the overall protein synthesis activity of a cell and activity on individual mRNAs are controlled by the allocation of limiting resources. Upon introduction of heterologous sequences into an organism-for example for the purposes of bioprocessing or synthetic biology-limiting resources may also become overstretched, thus negatively affecting both endogenous and heterologous gene expression. In this study, we present a mean-field model of translation in Saccharomyces cerevisiae for the investigation of two particular translational resources, namely ribosomes and aminoacylated tRNAs. We firstly use comparisons of experiments with heterologous sequences and simulations of the same conditions to calibrate our model, and then analyse the behaviour of the translational system in yeast upon introduction of different types of heterologous sequences. Our main findings are that: competition for ribosomes, rather than tRNAs, limits global translation in this organism; that tRNA aminoacylation levels exert, at most, weak control over translational activity; and that decoding speeds and codon adaptation exert strong control over local (mRNA specific) translation rates.     

Avoidance of antisense, antiterminator tRNA anticodons in vertebrate mitochondria

Protein synthesis (translation) stops at stop codons, codons not complemented by tRNA anticodons. tRNAs matching stops, antitermination (Ter) tRNAs, prevent translational termination, producing dysfunctional proteins. Genomes avoid tRNAs with anticodons whose complement (the anticodon of the ‘antisense’ tRNA) matches stops. This suggests that antisense tRNAs, which also form cloverleaves, are occasionally expressed. Mitochondrial antisense tRNA expression is plausible, because both DNA strands are transcribed as single RNAs, and tRNA structures signal RNA maturation. Results describe potential antisense Ter tRNAs in mammalian mitochondrial genomes detected by tRNAscan-SE, and evidence for adaptations preventing translational antitermination: genomes possessing Ter tRNAs use less corresponding stop codons; antisense Ter tRNAs form weaker cloverleaves than homologuous non-Ter antisense tRNAs; and genomic stop codon usages decrease with stabilities of codon-anticodon interactions and of Ter tRNA cloverleaves. This suggests that antisense tRNAs frequently function in translation. Results suggest that opposite strand coding is exceptional in modern genes, yet might be frequent for mitochondrial tRNAs. This adds antisense tRNA templating to other mitochondrial tRNA functions: sense tRNA templating, formation and regulation of secondary (light strand DNA) replication origins. Antitermination probably affects mitochondrial degenerative diseases and ageing: pathogenic mutations are twice as frequent in tRNAs with antisense Ter anticodons than in other tRNAs, and species lacking mitochondrial antisense Ter tRNAs have longer mean maximal lifespans than those possessing antisense Ter tRNAs.     

Evolution ofE.coli tRNATrp

Earlier studies (1) have shown there are direct correlations between the hydrophobicity ranking of most amino acids and their anticodonic nucleotides. However, four anticodonic assignments, i.e. those for Trp, Tyr, Ile and the XGA anticodons for Ser, did not correlate. It was our proposal that this failure to correlate was due to the fact that these assignments were made late, relative to the bulk of the assignments, in evolution through the mutation of existing tRNAs. We have shown (2) thatE. coli tRNA^{Ile 1} and tRNA^{Ile 2} were likely derived from tRNA^{Val 1} and tRNA^Lys respectively andE. coli tRNA^Tyr was possibly derived fromE. coli 5s rRNA or a common precursor with 5s rRNA (3). The fact that quite high homologies were observed in these comparisons is consistent with the late evolution of the tRNAs in question. We now examine the evolution ofE. coli tRNA^Trp by comparing its homology with otherE. coli tRNAs. The data suggest a possible evolutionary relationship withE. coli tRNA^Gly or tRNA^Arg. The data support the idea of the late assignment of anticodons to Trp.Deceased.     

ISOACCEPTOR tRNAs FOR GLUTAMATE, GLUTAMINE, ASPARTATE AND ASPARAGINE IN CALF BRAIN

Abstract— Glutamyl, glutaminyl. aspartyl and asparaginyl tRNAs of calf brain were analysed by reverse phase chromatography for isoacceptor tRNAs. The radioactivity profiles revealed two peaks for gluta-mate. three for glutamine, two for aspartate and one for asparagine. Comparison of brain tRNAs with tRNAs from other sources showed that glutamate and aspartate tRNAs of brain closely resembled a majority of other tRNAs in the number and relative abundance of isoacceptors. Glutamine and asparagine tRNAs from different sources exhibited more marked differences.     

Evolution of Pulmonate Gastropod Mitochondrial Genomes: Comparisons of Gene Organizations of Euhadra, Cepaea and Albinaria and Implications of Unusual Trna Secondary Structures

Complete gene organizations of the mitochondrial genomes of three pulmonate gastropods, Euhadra herklotsi, Cepaea nemoralis and Albinaria coerulea, permit comparisons of their gene organizations. Euhadra and Cepaea are classified in the same superfamily, Helicoidea, yet they show several differences in the order of tRNA and protein coding genes. Albinaria is distantly related to the other two genera but shares the same gene order in one part of its mitochondrial genome with Euhadra and in another part with Cepaea. Despite their small size (14.1-14.5 kbp), these snail mtDNAs encode 13 protein genes, two rRNA genes and at least 22 tRNA genes. These genomes exhibit several unusual or unique features compared to other published metazoan mitochondrial genomes, including those of other molluscs. Several tRNAs predicted from the DNA sequences possess bizarre structures lacking either the T stem or the D stem, similar to the situation seen in nematode mt-tRNAs. The acceptor stems of many tRNAs show a considerable number of mismatched basepairs, indicating that the RNA editing process recently demonstrated in Euhadra is widespread in the pulmonate gastropods. Strong selection acting on mitochondrial genomes of these animals would have resulted in frequent occurrence of the mismatched basepairs in regions of overlapping genes.     

Pathogenic mutations in antisense mitochondrial tRNAs

Pathogenic mutations in mitochondrial tRNAs are 6.5 times more frequent than in other mitochondrial genes. This suggests that tRNA mutations perturb more than one function. A potential additional tRNA gene function is that of templating for antisense tRNAs. Pathogenic mutations weaken cloverleaf secondary structures of sense tRNAs. Analyses here show similar effects for most antisense tRNAs, especially after adjusting for associations between sense and antisense cloverleaf stabilities. These results imply translational activity by antisense tRNAs. For sense tRNAs Ala and Ser UCN, pathogenicity associates as much with sense as with antisense cloverleaf formation. For tRNA Pro, pathogenicity seems associated only with antisense, not sense tRNA cloverleaf formation. Translational activity by antisense tRNAs is expected for the 11 antisense tRNAs processed by regular sense RNA maturation, those recognized by their cognate amino acid’s tRNA synthetase, and those forming relatively stable cloverleaves as compared to their sense counterpart. Most antisense tRNAs probably function routinely in translation and extend the tRNA pool (extension hypothesis); others do not (avoidance hypothesis). The greater the expected translational activity of an antisense tRNA, the more pathogenic mutations weaken its cloverleaf secondary structure. Some evidence for RNA interference, a more classical role for antisense tRNAs, exists only for tRNA Ser UCN. Mutation pathogenicity probably frequently results from a mixture of effects due to sense and antisense tRNA translational activity for many mitochondrial tRNAs. Genomic studies should routinely explore for translational activity by antisense tRNAs.     

Differential localization of nuclear-encoded tRNAs between the cytosol and mitochondrion in Leishmania tarentolae

All mitochondrial tRNAs of the kinetoplastid protozoan Leishmania tarentolae are encoded in the nucleus and are imported from the cytosol into the mitochondrion. We previously reported the partitioning of five tRNAs and found that all were shared between the two compartments to different extents. To increase our knowledge of the tRNAs of this organism, and to attempt to understand the signals involved in their subcellular localization, a method to RT-PCR amplify new tRNAs was developed. Various tRNAs were 3' polyadenylated and reverse transcribed with a sequence-tagged primer. The cDNA was tagged by ligation to an anchor oligonucleotide, and the resulting double-tagged cDNA was amplified by PCR. Four new tRNAs were obtained, bringing to 20 the total number of L. tarentolae tRNAs identified to date. The subcellular localization of 17 tRNAs was quantitatively analyzed by two-dimensional gel electrophoresis and northern hybridization. In general, the previously suggested operational classification of tRNAs into three groups (mainly cytosolic, mainly mitochondrial, and shared between the two compartments) is still valid, but the relative abundance of each tRNA in the cytosol or mitochondrion varied greatly as did the level of expression.     

Comparison of dissimilarity patterns of E coli, yeast and mammalian tRNAs.

A number of experimental approaches have been developed for identification of recognition (identity) sites in tRNAs. Along with them a theoretical methodology has been proposed by McClain et al that is based on concomitant analysis of all tRNA sequences from a given species. This approach allows an evaluation of nucleotide combinations present in isoacceptor tRNAs specific for the given amino acid, and not present in equivalent positions in cloverleaf structure in other tRNAs of the same organism. These elements predicted from computer analysis of the databank could be tested experimentally for their participation in forming recognition sites. The correlation between theoretical predictions and experimental data appeared promising. The aim of the present work consisted of introducing further improvements into McClain's procedure by: i), introducing into analysis a variable region in tRNAs which had not been previously considered; to accomplish this, 'normalization' of variable nucleotides was suggested, based on primary and tertiary structures of tRNAs; ii), developing a new procedure for comparison of patterns for synonymous and non-synonymous tRNAs from different organisms; iii), analysis of 3- and 4-positional contacts between tRNAs and enzymes in addition to a formerly used 2-positional model. A systematic application of McClain's procedure to mammalian, yeast and E coli tRNAs led to the following results: i), imitancy patterns for non-synonymous tRNAs of any amino acid specificity and from any organisms analysed so far overlap by no more than 30%, providing a structural basis for discrimination with high fidelity between cognate and non-cognate tRNAs; ii), the predicted identity sites are non-randomly distributed within tRNA molecules; the dominant role is ascribed to only two regions--anticodon and amino acid stem which are located far apart from one another at extremes of all tRNA molecules; iii), the imitancy patterns for synonymous tRNAs in lower (yeast) and higher (mammalian) eukaryotes are similar but not identical; iv), distribution of predicted identity sites in the cloverleaf structure in prokaryotes and eukaryotes is essentially different: in eubacterial tRNAs the major role in recognition plays anticodon and/or amino acid acceptor stem, whereas in eukaryotic (both unicellular and multicellular) tRNAs the remaining part of the molecules is also involved in recognition; v), the imitancy patterns of synonymous tRNAs from prokaryotes and eukaryotes are dissimilar, this observation leads to the prediction that the tRNA identity sites for the same amino acid in prokaryotes and eukaryotes may differ.