Studying the evolutionary relationships and phylogenetic trees of 21 groups of tRNA sequences based on complex networks

To find out the evolutionary relationships among different tRNA sequences of 21 amino acids, 22 networks are constructed. One is constructed from whole tRNAs, and the other 21 networks are constructed from the tRNAs which carry the same amino acids. A new method is proposed such that the alignment scores of any two amino acids groups are determined by the average degree and the average clustering coefficient of their networks. The anticodon feature of isolated tRNA and the phylogenetic trees of 21 group networks are discussed. We find that some isolated tRNA sequences in 21 networks still connect with other tRNAs outside their group, which reflects the fact that those tRNAs might evolve by intercrossing among these 21 groups. We also find that most anticodons among the same cluster are only one base different in the same sites when S ≥ 70, and they stay in the same rank in the ladder of evolutionary relationships. Those observations seem to agree on that some tRNAs might mutate from the same ancestor sequences based on point mutation mechanisms.     

tRNA序列的进化网络研究

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

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]     

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.     

Parallel evolution of truncated transfer RNA genes in arachnid mitochondrial genomes

The cloverleaf secondary structure of transfer RNA (tRNA) is highly conserved across all forms of life. Here, we provide sequence data and inferred secondary structures for all tRNA genes from 8 new arachnid mitochondrial genomes, including representatives from 6 orders. These data show remarkable reductions in tRNA gene sequences, indicating that T-arms are missing from many of the 22 tRNAs in the genomes of 4 out of 7 orders of arachnids. Additionally, all opisthothele spiders possess some tRNA genes that lack sequences that could form well-paired aminoacyl acceptor stems. We trace the evolution of T-arm loss onto phylogenies of arachnids and show that a genome-wide propensity to lose sequences that encode canonical cloverleaf structures likely evolved multiple times within arachnids. Mapping of structural characters also shows that certain tRNA genes appear more evolutionarily prone to lose the sequence coding for the T-arm and that once a T-arm is lost, it is not regained. We use tRNA structural data to construct a phylogeny of arachnids and find high bootstrap support for a clade that is not supported in phylogenies that are based on more traditional morphological characters. Together, our data demonstrate variability in structural evolution among different tRNAs as well as evidence for parallel evolution of the loss of sequence coding for tRNA arms within an ancient and diverse group of animals.     

Functional defects of pathogenic human mitochondrial tRNAs related to structural fragility

Aminoacylation of transfer RNAs (tRNAs) is essential for protein synthesis. A growing number of human diseases correlate with point mutations in tRNA genes within the mitochondrial genome. These tRNAs have unique sequences that suggest they have fragile structures. However, the structural significance of pathology-related tRNA mutations and their effects on molecular function have not been explored. Here, opthalmoplegia related mutants of a human mitochondrial tRNA have been investigated. Each mutation replaces either an A-U or G-C pair in the predicted secondary structure with an A-C pair. Aminoacylation of each mutant tRNA was severely attenuated. Moreover, each strongly inhibited aminoacylation of the wild type substrate, suggesting that the effects of these mutations might not be bypassed in the potentially heteroplasmic environment of mitochondria. The function of mutant tRNAs was rescued by single compensatory mutations that restored Watson-Crick base pairing and reintroduced stability into regions of predicted secondary structure, even though the pairs introduced were different from those found in the wild type tRNA. Thus, functional defects caused by a subset of pathogenic mutations may result from the inherent structural fragility of human mitochondrial tRNAs.     

The 51-63 base pair of tRNA confers specificity for binding by EF-Tu

Elongation factor Tu (EF-Tu) exhibits significant specificity for the different elongator tRNA bodies in order to offset its variable affinity to the esterified amino acid. Three X-ray cocrystal structures reveal that while most of the contacts with the protein involve the phosphodiester backbone of tRNA, a single hydrogen bond is observed between the Glu390 and the amino group of a guanine in the 51-63 base pair in the T-stem of tRNA. Here we show that the Glu390Ala mutation of Thermus thermophilus EF-Tu selectively destabilizes binding of those tRNAs containing a guanine at either position 51 or 63 and that mutagenesis of the 51-63 base pair in several tRNAs modulates their binding affinities to EF-Tu. A comparison of Escherichia coli tRNA sequences suggests that this specificity mechanism is conserved across the bacterial domain. While this contact is an important specificity determinant, it is clear that others remain to be identified.     

Disrupted tRNA Gene Diversity and Possible Evolutionary Scenarios

The following unusual tRNAs have recently been discovered in the genomes of Archaea and primitive Eukaryota: multiple-intron-containing tRNAs, which have more than one intron; split tRNAs, which are produced from two pieces of RNA transcribed from separate genes; tri-split tRNAs, which are produced from three separate genes; and permuted tRNA, in which the 5′ and 3′ halves are encoded with permuted orientations within a single gene. All these disrupted tRNA genes can form mature contiguous tRNA, which is aminoacylated after processing by cis or trans splicing. The discovery of such tRNA disruptions has raised the question of when and why these complex tRNA processing pathways emerged during the evolution of life. Many previous reports have noted that tRNA genes contain a single intron in the anticodon loop region, a feature common throughout all three domains of life, suggesting an ancient trait of the last universal common ancestor. In this context, these unique tRNA disruptions recently found only in Archaea and primitive Eukaryota provide new insight into the origin and evolution of tRNA genes, encouraging further research in this field. In this paper, we summarize the phylogeny, structure, and processing machinery of all known types of disrupted tRNAs and discuss possible evolutionary scenarios for these tRNA genes.     

The impact of including tRNA content on the optimality of the genetic code

Statistical and biochemical studies have revealed nonrandom patterns in codon assignments. The canonical genetic code is known to be highly efficient in minimizing the effects of mistranslational errors and point mutations, since it is known that, when an amino acid is converted to another due to error, the biochemical properties of the resulted amino acid are usually very similar to those of the original one. In this study, we have taken into consideration both relative frequencies of amino acids and relative gene copy frequencies of tRNAs in genomic sequences in order to introduce a fitness function which models the mistranslational probabilities more accurately in modern organisms. The relative gene copy frequencies of tRNAs are used as estimates of the tRNA content. We also altered the rule previously used for the calculation of the probabilities of single base mutation occurrences. Our model signifies higher optimality of the genetic code towards load minimization and suggests the presence of a coevolution of tRNA frequency and the genetic code.     

Epistasis can lead to fragmented neutral spaces and contingency in evolution

In evolution, the effects of a single deleterious mutation can sometimes be compensated for by a second mutation which recovers the original phenotype. Such epistatic interactions have implications for the structure of genome space--namely, that networks of genomes encoding the same phenotype may not be connected by single mutational moves. We use the folding of RNA sequences into secondary structures as a model genotype-phenotype map and explore the neutral spaces corresponding to networks of genotypes with the same phenotype. In most of these networks, we find that it is not possible to connect all genotypes to one another by single point mutations. Instead, a network for a phenotypic structure with n bonds typically fragments into at least 2(n) neutral components, often of similar size. While components of the same network generate the same phenotype, they show important variations in their properties, most strikingly in their evolvability and mutational robustness. This heterogeneity implies contingency in the evolutionary process.     

A disease-causing point mutation in human mitochondrial tRNAMet rsults in tRNA misfolding leading to defects in translational initiation and elongation

The mitochondrial tRNA genes are hot spots for mutations that lead to human disease. A single point mutation (T4409C) in the gene for human mitochondrial tRNA(Met) (hmtRNA(Met)) has been found to cause mitochondrial myopathy. This mutation results in the replacement of U8 in hmtRNA(Met) with a C8. The hmtRNA(Met) serves both in translational initiation and elongation in human mitochondria making this tRNA of particular interest in mitochondrial protein synthesis. Here we show that the single 8U-->C mutation leads to a failure of the tRNA to respond conformationally to Mg(2+). This mutation results in a drastic disruption of the structure of the hmtRNA(Met), which significantly reduces its aminoacylation. The small fraction of hmtRNA(Met) that can be aminoacylated is not formylated by the mitochondrial Met-tRNA transformylase preventing its function in initiation, and it is unable to form a stable ternary complex with elongation factor EF-Tu preventing any participation in chain elongation. We have used structural probing and molecular reconstitution experiments to examine the structures formed by the normal and mutated tRNAs. In the presence of Mg(2+), the normal tRNA displays the structural features expected of a tRNA. However, even in the presence of Mg(2+), the mutated tRNA does not form the cloverleaf structure typical of tRNAs. Thus, we believe that this mutation has disrupted a critical Mg(2+)-binding site on the tRNA required for formation of the biologically active structure. This work establishes a foundation for understanding the physiological consequences of the numerous mitochondrial tRNA mutations that result in disease in humans.     

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     

Stochastic context-free grammars for tRNA modeling.

Stochastic context-free grammars (SCFGs) are applied to the problems of folding, aligning and modeling families of tRNA sequences. SCFGs capture the sequences' common primary and secondary structure and generalize the hidden Markov models (HMMs) used in related work on protein and DNA. Results show that after having been trained on as few as 20 tRNA sequences from only two tRNA subfamilies (mitochondrial and cytoplasmic), the model can discern general tRNA from similar-length RNA sequences of other kinds, can find secondary structure of new tRNA sequences, and can produce multiple alignments of large sets of tRNA sequences. Our results suggest potential improvements in the alignments of the D- and T-domains in some mitochondrial tRNAs that cannot be fit into the canonical secondary structure.     

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.     

Recognition of exon-intron boundaries by the Halobacterium volcanii tRNA intron endonuclease.

The intron-containing tRNA(Trp) precursor from Halobacterium volcanii, like many intron-containing archaebacterial precursor tRNAs, can assume a structure in which the two intron endonuclease cleavage sites are localized in two three-nucleotide loops separated by four base pairs. To investigate the role of this structure in cleavage by the halophilic endonuclease, a series of mutant tRNA(Trp) RNAs were prepared and evaluated as substrates. We find that alterations in this structure result in the loss of cleavage at both 5' and 3' sites. Cleavage of a 35-nucleotide model RNA substrate, containing only these features, demonstrates that sequences and structures present at the exon-intron boundaries are sufficient for recognition and cleavage. We have also examined the mechanism used by the halophilic endonuclease to identify the cleavage sites. Addition of a single base, or a base pair in the anticodon stem above the cleavage sites, does not affect the cleavage site selection. The addition of nucleotides between the two cleavage sites significantly decreases cleavage efficiency and has an effect on the cleavage site selection. These results demonstrate that the halophilic endonuclease requires a defined structure at the exon-intron boundaries and does not identify its cleavage sites by a measurement mechanism like that employed by eukaryotic tRNA intron endonucleases.     

Anticodons, frameshifts, and hidden periodicities in tRNA sequences

Fourier analysis of the short-range periodicities for the complete set of sequences coding for tRNA genes in genome of Bacillus subtilis proves that periodicities with periods p = 2, 3, 4, and 6 sites are the inherent properties of tRNAs. The related periodicities should be understood in a broad statistical sense and their identifying needs the elaborate statistical methods. To improve the statistics, the analysis of significant periodicities was performed for the binary R-Y, S-W, and K-M sequences. Generally, such short-range periodicities are produced via biased positioning of particular nucleotides rather than via the tandem multiplication and subsequent modifications of repeats, though the latter mechanism may also be realized. Quasi-coherently piercing long segments of tRNA, the short-range periodicities create the effective long-range structural coupling between the acceptor stem and the anticodon loop and may participate in the mechanisms of molecular recognition. The periodicities with p = 2 and 4 provide the natural ground for the translation with spontaneous or programmed frameshifting and are present in tRNAs decoding the most frameshift-prone codons. The observation of short-range periodicities suggests that the mechanisms of amino-acylation of tRNAs and codon-anticodon pairing are not independent. Their study may also provide the important information related to the origin and evolution of the genetic code.     

Towards a comprehensive picture of alloacceptor tRNA remolding in metazoan mitochondrial genomes

Remolding of tRNAs is a well-documented process in mitochondrial genomes that changes the identity of a tRNA. It involves a duplication of a tRNA gene, a mutation that changes the anticodon and the loss of the ancestral tRNA gene. The net effect is a functional tRNA that is more closely related to tRNAs of a different alloacceptor family than to tRNAs with the same anticodon in related species. Beyond being of interest for understanding mitochondrial tRNA function and evolution, tRNA remolding events can lead to artifacts in the annotation of mitogenomes and thus in studies of mitogenomic evolution. Therefore, it is important to identify and catalog these events. Here we describe novel methods to detect tRNA remolding in large-scale data sets and apply them to survey tRNA remolding throughout animal evolution. We identify several novel remolding events in addition to the ones previously mentioned in the literature. A detailed analysis of these remoldings showed that many of them are derived from ancestral events.     

Missense and nonsense suppressors can correct frameshift mutations

Missense and nonsense suppressor tRNAs, selected for their ability to read a new triplet codon, were observed to suppress one or more frameshift mutations in trpA of Escherichia coli. Two of the suppressible frameshift mutants, trpA8 and trpA46AspPR3, were cloned, sequenced, and found to be of the +1 type, resulting from the insertion of four nucleotides and one nucleotide, respectively. Twenty-two suppressor tRNAs were examined, 20 derived from one of the 3 glycine isoacceptor species, one from lysT, and one from trpT. The sequences of all but four of the mutant tRNAs are known, and two of those four were converted to suppressor tRNAs that were subsequently sequenced. Consideration of the coding specificities and anticodon sequences of the suppressor tRNAs does not suggest a unitary mechanism of frameshift suppression. Rather, the results indicate that different suppressors may shift frame according to different mechanisms. Examination of the suppression windows of the suppressible frameshift mutations indicates that some of the suppressors may work at cognate codons, either in the 0 frame or in the +1 frame, and others may act at noncognate codons (in either frame) by some as-yet-unspecified mechanism. Whatever the mechanisms, it is clear that some +1 frameshifting can occur at non-monotonous sequences. A striking example of a frameshifting missense suppressor is a mutant lysine tRNA that differs from wild-type lysine tRNA by only a single base in the amino acid acceptor stem, a C to U70 transition that results in a G.U base pair. It is suggested that when this mutant lysine tRNA reads its cognate codon, AAA, the presence of the G.U base pair sometimes leads either to a conformational change in the tRNA or to an altered interaction with some component of the translation machinery involved in translocation, resulting in a shift of reading frame. In general, the results indicate that translocation is not simply a function of anticodon loop size, that different frameshifting mechanisms may operate with different tRNAs, and that conformational features, some far removed from the anticodon region, are involved in maintaining fidelity in translocation.