tRNA序列的进化网络研究

根据tRNA序列的反密码子,把3420 条tRNA 序列分成了21组,其中包括20 组氨基酸及1 组Stop。通过tRNA序列的相似度构建了1 组整体网络和21 组子网络,并计算了它们在不同相似度下的平均度、平均聚类系数以及平均最短路径。通过分析、比较和讨论网络中的三个重要参数,进一步说明点突变是tRNA序列进化的重要机制,并反映了它们的进化近似符合中性理论,且在同一组氨基酸和Stop 内的tRNA序列在进化历史上的同源关系更密切;同时表明了tRNA 序列在进化过程中具有自相似性。     

Comparing two evolutionary mechanisms of modern tRNAs

All modern tRNA gene families have a high similarity in their primary structure, and share the same cloverleaf secondary structure and an inverted L tertiary structure, which provide the clues for the study of their origin and evolution. There are two important mechanisms of the tRNA sequences evolution. One is point mutation, another is complementary duplication method. Both of them are supported by some specific examples. To find out the superior one of the two mechanisms or find out the most suitable mechanism for modern tRNAs evolution, we constructed two types of networks, parallel and antiparallel networks, based on the two mechanisms respectively, and then compared the degree distribution and clustering coefficient of networks constructed by the tRNA sequences of the single anticodon group, single isoaccepting group, and the whole tRNAs group of the two types of networks. The result of the comparison seems consistent with the idea that modern tRNA sequences evolved primarily by the mechanism of complementary method, and point mutation is an important and indispensable auxiliary mechanism during the evolutionary event.     

Relationships Among Isoacceptor tRNAs Seems to Support the Coevolution Theory of the Origin of the Genetic Code

A new method for looking at relationships between nucleotide sequences has been used to analyze divergence both within and between the families of isoaccepting tRNA sets. A dendrogram of the relationships between 21 tRNA sets with different amino acid specificities is presented as the result of the analysis. Methionine initiator tRNAs are included as a separate set. The dendrogram has been interpreted with respect to the final stage of the evolutionary pathway with the development of highly specific tRNAs from ambiguous molecular adaptors. The location of the sets on the dendrogram was therefore analyzed in relation to hypotheses on the origin of the genetic code: the coevolution theory, the physicochemical hypothesis, and the hypothesis of ambiguity reduction of the genetic code. Pairs of 16 sets of isoacceptor tRNAs, whose amino acids are in biosynthetic relationships, occupied contiguous positions on the dendrogram, thus supporting the coevolution theory of the genetic code. Received: 4 May 1998 / Accepted: 11 July 1998     

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.     

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.     

The non-monophyletic origin of the tRNA molecule

The hypothesis that the tRNA molecule may have originated from the assembly of two similar RNA hairpin structures is utilised to understand the evolutionary period in which this molecule originated. Consistent with the exon theory of genes is the observation that the introns in tRNA genes are found almost exclusively in the anticodon loop and "stitched together" the two halves of the molecule, which originally may have been simply two hairpin structures and which can still be observed in the three-dimensional structure of tRNAs. This theory therefore considers these hairpin structures as minigenes on which complex protein synthesis may have been achieved. This in turn leads to the belief that the organisation of the genetic code may have been determined by use of the hairpin structures but not the complete tRNA molecule. In view of this, it can be conjectured that tRNA molecules might have been assembled only after the establishment of the main phyletic lines. If this is all true, then the origin of the tRNA molecule might have been non-monophyletic, i.e. a tRNA specific for a certain amino acid might have been assembled in different phyletic lines with a second and different hairpin structure. This leads to the belief that tRNAs specific for different amino acids but belonging to the same phyletic line might have been more similar to one another than to tRNAs specific for the same amino acid but belonging to different phyletic lines. This prediction seems to be supported by phylogenetic analysis making major use of the bootstrap technique performed on the tRNA sequences and by analysis already existing in the literature which supports the non-monophyletic origin of the tRNA molecule. The main conclusion of this paper is that if the tRNA molecule was assembled in the main phyletic lines this would imply a still rapidly evolving translation apparatus which, in turn, seems to imply that the last universal common ancestor was a progenote.     

The phylogeny of trna molecules and the origin of the genetic code

The evolutionary relationships between transfer RNA (tRNA) molecules are analyzed by parsimony algorithms. The position of the topologies expected on the basis of the hypotheses made to explain the origin of the genetic code, on the frequency distribution of all the possible tree topologies of the evolutionary relationships between tRNAs seems to lead to the following conclusion: The hypothesis (Wong, J. X, Proc. Natl. Acad. Sci. USA, 1975, 72: 1909–1912) that sees the genetic code as a map of the biosynthetic relationships between amino acids seems to occupy a statistically significant position on these frequency distributions, thus reflecting a significant part of the tRNA phylogeny.     

Formal Proof that the Split Genes of tRNAs of <Emphasis Type="Italic">Nanoarchaeum equitans</Emphasis> Are an Ancestral Character

A proof is given that the genes of the tRNA molecule of Nanoarchaeum equitans split into the 5′ and 3′ halves are an ancestral trait. First, the existence of a natural succession of evolutionary stages will be proven, formed in the order of the three gene structures of tRNAs known today: (i) the split genes of tRNAs, (ii) the genes of tRNAs with introns, and (iii) the genes of tRNAs continuously codifying for the tRNA molecule. This succession of evolutionary stages identifies the split genes of tRNAs as a pleisiomorphic character. The proof that this succession of evolutionary stages is, moreover, true is performed by proving that all the possible remaining five successions of evolutionary stages are false. Indeed, the succession of evolutionary stages considering split genes as a derived character turns out to be false in that the increase in complexity inherent to this succession cannot be justified by the split genes of tRNAs because these could not have conferred any selective advantage justifying this increase in complexity. Furthermore, genetic drift is unable to explain the evolution of split genes of tRNAs because of the enormous genetic effective size of the population observed in these organisms. The remaining four successions of evolutionary stages are also false because: (i) they are not natural successions of evolutionary stages, (ii) the absolute observed frequencies of these evolutionary stages are such as to exclude categorically that they might be natural successions of evolutionary stages, and also (iii) two of these are falsified by the fact that they do not place the evolutionary stage of genes of tRNAs with introns in a close evolutionary relationship with that of the split genes of tRNAs which can, instead, be proven to have a close evolutionary link. Therefore, there remains only the succession of evolutionary stages considering the split genes of tRNAs codifying for the 5′ and 3′ halves, as a pleisiomorphic character, as the only succession compatible with all the arguments presented in this article and as the one that actually operated during the evolution of the tRNA molecule. This proof has two very important implications. One regards how the tRNA molecule originated; considering how tRNA originated as the union of two hairpin-like structures, the split genes of tRNAs might be the transition stage through which the evolution of this molecule passed. The other regards when the genes of tRNAs originated, reaching the conclusion that the origin of these genes is polyphyletic, i.e. not monophyletic and hence contrary to the assumptions of the current paradigm.     

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.     

Comparative and evolutionary analysis in natural diploid and tetraploid weather loach Misgurnus anguillicaudatus based on cytochrome b sequence data in central China

To obtain the phylogenetic relationship between diploid and tetraploid Misgurnus anguillicaudatus, the mitochondrial cyt b gene in the diploid and tetraploid weather loach were isolated and sequenced. The DNA sequences were analyzed using MEGA 3.0 software to determine the phylogenetic relationship. Forty-five variable sites among cyt b gene sequences and 18 amino acid substitutions occurred within the diploid and tetraploid loaches as deduced from the nucleotide sequences analysis of the cyt b gene. The nucleotide pairwise distance between diploid and tetraploid loach ranged from 0.001 to 0.025. Phylogenetic analysis revealed evolutionary relationships between diploid and tetraploid loach. Our results indicated a significant difference between diploid and tetraploid loach about the cyt b gene. AMOVA analysis indicated that there were no significant genetic variations within diploid loaches (Fst = 0.2529, P > 0.05) and within tetraploid loaches (Fst = 0.0564, P > 0.05), neither. However, significant genetic differences were found between diploid and tetraploid loaches (Fst = 0.7634, P < 0.05). Thus, it is concluded that no reproductive isolation was found within the same cytotypes of different localities, but there was reproductive isolation between these two cytotypes. The diploid loach existed before the tetraploid loach in nature. The present study is the first to describe the phylogenetic relationships of natural polyploidy weather loach using mtDNA cyt b gene.     

Phylogeny and evolution of chloroplast tRNAs in Adoxaceae

Chloroplasts are semiautonomous organelles found in photosynthetic plants. The major functions of chloroplasts include photosynthesis and carbon fixation, which are mainly regulated by its circular genomes. In the highly conserved chloroplast genome, the chloroplast transfer RNA genes (cp tRNA) play important roles in protein translation within chloroplasts. However, the evolution of cp tRNAs remains unclear. Thus, in the present study, we investigated the evolutionary characteristics of chloroplast tRNAs in five Adoxaceae species using 185 tRNA gene sequences. In total, 37 tRNAs encoding 28 anticodons are found in the chloroplast genome in Adoxaceae species. Some consensus sequences are found within the Ψ‐stem and anticodon loop of the tRNAs. Some putative novel structures were also identified, including a new stem located in the variable region of tRNA^Tyr in a similar manner to the anticodon stem. Furthermore, phylogenetic and evolutionary analyses indicated that synonymous tRNAs may have evolved from multiple ancestors and frequent tRNA duplications during the evolutionary process may have been primarily caused by positive selection and adaptive evolution. The transition and transversion rates are uneven among different tRNA isotypes. For all tRNAs, the transition rate is greater with a transition/transversion bias of 3.13. Phylogenetic analysis of cp tRNA suggested that the type I introns in different taxa (including eukaryote organisms and cyanobacteria) share the conserved sequences “U‐U‐x2‐C” and “U‐x‐G‐x2‐T,” thereby indicating the diverse cyanobacterial origins of organelles. This detailed study of cp tRNAs in Adoxaceae may facilitate further investigations of the evolution, phylogeny, structure, and related functions of chloroplast tRNAs.     

The RNA origin of transfer RNA aminoacylation and beyond

Aminoacylation of tRNA is an essential event in the translation system. Although in the modern system protein enzymes play the sole role in tRNA aminoacylation, in the primitive translation system RNA molecules could have catalysed aminoacylation onto tRNA or tRNA-like molecules. Even though such RNA enzymes so far are not identified from known organisms, in vitro selection has generated such RNA catalysts from a pool of random RNA sequences. Among them, a set of RNA sequences, referred to as flexizymes (Fxs), discovered in our laboratory are able to charge amino acids onto tRNAs. Significantly, Fxs allow us to charge a wide variety of amino acids, including those that are non-proteinogenic, onto tRNAs bearing any desired anticodons, and thus enable us to reprogramme the genetic code at our will. This article summarizes the evolutionary history of Fxs and also the most recent advances in manipulating a translation system by integration with Fxs.     

tRNA Creation by Hairpin Duplication

Many studies have suggested that the modern cloverleaf structure of tRNA may have arisen through duplication of a primordial hairpin, but the timing of this duplication event has been unclear. Here we measure the level of sequence identity between the two halves of each of a large sample of tRNAs and compare this level to that of chimeric tRNAs constructed either within or between groups defined by phylogeny and/or specificity. We find that actual tRNAs have significantly more matches between the two halves than do random sequences that can form the tRNA structure, but there is no difference in the average level of matching between the two halves of an individual tRNA and the average level of matching between the two halves of the chimeric tRNAs in any of the sets we constructed. These results support the hypothesis that the modern tRNA cloverleaf arose from a single hairpin duplication prior to the divergence of modern tRNA specificities and the three domains of life. [Reviewing Editor: Dr. Niles Lehman]     

Transfer RNA Identity Change in Anticodon Variants of E. coli tRNAPhe in Vivo

The anticodon sequence is a major recognition element for most aminoacyl-tRNA synthetases. We investigated the in vivo effects of changing the anticodon on the aminoacylation specificity in the example of E. coli tRNA^Phe. Constructing different anticodon mutants of E. coli tRNA^Phe by site-directed mutagenesis, we isolated 22 anticodon mutant tRNA^Phe; the anticodons corresponded to 16 amino acids and an opal stop codon. To examine whether the mutant tRNAs had changed their amino acid acceptor specificity in vivo, we tested the viability of E. coli strains containing these tRNA^Phe genes in a medium which permitted tRNA induction. Fourteen mutant tRNA genes did not affect host viability. However, eight mutant tRNA genes were toxic to the host and prevented growth, presumably because the anticodon mutants led to translational errors. Many mutant tRNAs which did not affect host viability were not aminoacylated in vivo. Three mutant tRNAs containing anticodon sequences corresponding to lysine (UUU), methionine (CAU) and threonine (UGU) were charged with the amino acid corresponding to their anticodon, but not with phenylalanine. These three tRNAs and tRNA^Phe are located in the same cluster in a sequence similarity dendrogram of total E. coli tRNAs. The results support the idea that such tRNAs arising from in vivo evolution are derived by anticodon change from the same ancestor 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.     

Molecular basis for the differential interaction of plant mitochondrial VDAC proteins with tRNAs

In plants, the voltage-dependent anion-selective channel (VDAC) is a major component of a pathway involved in transfer RNA (tRNA) translocation through the mitochondrial outer membrane. However, the way in which VDAC proteins interact with tRNAs is still unknown. Potato mitochondria contain two major mitochondrial VDAC proteins, VDAC34 and VDAC36. These two proteins, composed of a N-terminal α-helix and of 19 β-strands forming a β-barrel structure, share 75% sequence identity. Here, using both northwestern and gel shift experiments, we report that these two proteins interact differentially with nucleic acids. VDAC34 binds more efficiently with tRNAs or other nucleic acids than VDAC36. To further identify specific features and critical amino acids required for tRNA binding, 21 VDAC34 mutants were constructed and analyzed by northwestern. This allowed us to show that the β-barrel structure of VDAC34 and the first 50 amino acids that contain the α-helix are essential for RNA binding. Altogether the work shows that during evolution, plant mitochondrial VDAC proteins have diverged so as to interact differentially with nucleic acids, and this may reflect their involvement in various specialized biological functions.     

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.     

The phylogeny of trnas seems to confirm the predictions of the coevolution theory of the origin of the genetic code

An extensive analysis of the evolutionary relationships existing between transfer RNAs, performed using parsimony algorithms, is presented. After building up an estimate of the tRNA ancestral sequences, these sequences are then compared using certain methods. The results seem to suggest that the coevolution hypothesis (Wong, J.T., 1975, Proc. Natl. Acad. Sci. USA 72, 1909–1912) that sees the genetic code as a map of the biosynthetic relationships between amino acids is further supported by these results, as compared to the hypotheses that see the physicochemical properties of amino acids as the main adaptative theme that led to the structuring of the genetic code.     

Evolution of tRNA recognition systems and tRNA gene sequences

The aminoacylation of tRNAs by the aminoacyl-tRNA synthetases recapitulates the genetic code by dictating the association between amino acids and tRNA anticodons. The sequences of tRNAs were analyzed to investigate the nature of primordial recognition systems and to make inferences about the evolution of tRNA gene sequences and the evolution of the genetic code. Evidence is presented that primordial synthetases recognized acceptor stem nucleotides prior to the establishment of the three major phylogenetic lineages. However, acceptor stem sequences probably did not achieve a level of sequence diversity sufficient to faithfully specify the anticodon assignments of all 20 amino acids. This putative bottleneck in the evolution of the genetic code may have been alleviated by the advent of anticodon recognition. A phylogenetic analysis of tRNA gene sequences from the deep Archaea revealed groups that are united by sequence motifs which are located within a region of the tRNA that is involved in determining its tertiary structure. An association between the third anticodon nucleotide (N36) and these sequence motifs suggests that a tRNA-like structure existed close to the time that amino acid-anticodon assignments were being established. The sequence analysis also revealed that tRNA genes may evolve by anticodon mutations that recruit tRNAs from one isoaccepting group to another. Thus tRNA gene evolution may not always be monophyletic with respect to each isoaccepting group.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 Correspondence to: M.E. Saks     

Undetected antisense tRNAs in mitochondrial genomes?

Background  

The hypothesis that both mitochondrial (mt) complementary DNA strands of tRNA genes code for tRNAs (sense-antisense coding) is explored. This could explain why mt tRNA mutations are 6.5 times more frequently pathogenic than in other mt sequences. Antisense tRNA expression is plausible because tRNA punctuation signals mt sense RNA maturation: both sense and antisense tRNAs form secondary structures potentially signalling processing. Sense RNA maturation processes by default 11 antisense tRNAs neighbouring sense genes. If antisense tRNAs are expressed, processed antisense tRNAs should have adapted more for translational activity than unprocessed ones. Four tRNA properties are examined: antisense tRNA 5′ and 3′ end processing by sense RNA maturation and its accuracy, cloverleaf stability and misacylation potential.