Evolutionary patterns in the sequence and structure of transfer RNA: early origins of archaea and viruses

Transfer RNAs (tRNAs) are ancient molecules that are central to translation. Since they probably carry evolutionary signatures that were left behind when the living world diversified, we reconstructed phylogenies directly from the sequence and structure of tRNA using well-established phylogenetic methods. The trees placed tRNAs with long variable arms charging Sec, Tyr, Ser, and Leu consistently at the base of the rooted phylogenies, but failed to reveal groupings that would indicate clear evolutionary links to organismal origin or molecular functions. In order to uncover evolutionary patterns in the trees, we forced tRNAs into monophyletic groups using constraint analyses to generate timelines of organismal diversification and test competing evolutionary hypotheses. Remarkably, organismal timelines showed Archaea was the most ancestral superkingdom, followed by viruses, then superkingdoms Eukarya and Bacteria, in that order, supporting conclusions from recent phylogenomic studies of protein architecture. Strikingly, constraint analyses showed that the origin of viruses was not only ancient, but was linked to Archaea. Our findings have important implications. They support the notion that the archaeal lineage was very ancient, resulted in the first organismal divide, and predated diversification of tRNA function and specificity. Results are also consistent with the concept that viruses contributed to the development of the DNA replication machinery during the early diversification of the living world.     

Evolutionary changes in the genetic code

The genetic code has been influenced by directional mutation pressure affecting the base composition of DNA, sometimes in the direction of increased GC content and at other times, in the direction of AT. Such pressure led to changes in species-specific usages of codons and tRNA anticodons, and also in amino acid assignments of codons in mitochondria and in several intact organisms. These code changes are probably recent evolutionary events. The genetic code is not 'frozen', but instead it is still evolving.     

A new noncanonical nuclear genetic code: translation of UAA into glutamate

Deviant genetic codes reported in ciliates share the same feature: one (UGA) or two (UAR) of the three canonical stop codons are translated into one particular amino acid. In many genera, such as Oxytricha, Paramecium, and Tetrahymena, UAR codons are translated into glutamine. UGA is translated into cysteine in Euplotes or into tryptophan in Colpoda inflata and Blepharisma americanum. Here, we show that three peritrich species (Vorticella microstoma, Opisthonecta henneguyi, and Opisthonecta matiensis) translate UAA into glutamate and that at least UAA in O. matiensis is decoded through a mutant suppressor-like tRNA. This kind of genetic code has never been reported for any living organism. Phylogenetic analysis with alpha-tubulin sequences corroborates that peritrichs, peniculines (Paramecium), and hymenostomates (Tetrahymena) form a monophyletic group (class Oligohymenophorea). The differential translation (glu/gln) of UAR codons, the monophyly of the Oligohymenophorea, and the common evolutionary origin of glutamate and glutamine suggest that deviant genetic codes of present-day oligohymenophoreans could have the same origin.     

Evidence from glycine transfer RNA of a frozen accident at the dawn of the genetic code

Background

Transfer RNA (tRNA) is the means by which the cell translates DNA sequence into protein according to the rules of the genetic code. A credible proposition is that tRNA was formed from the duplication of an RNA hairpin half the length of the contemporary tRNA molecule, with the point at which the hairpins were joined marked by the canonical intron insertion position found today within tRNA genes. If these hairpins possessed a 3'-CCA terminus with different combinations of stem nucleotides (the ancestral operational RNA code), specific aminoacylation and perhaps participation in some form of noncoded protein synthesis might have occurred. However, the identity of the first tRNA and the initial steps in the origin of the genetic code remain elusive.

Results

Here we show evidence that glycine tRNA was the first tRNA, as revealed by a vestigial imprint in the anticodon loop sequences of contemporary descendents. This provides a plausible mechanism for the missing first step in the origin of the genetic code. In 448 of 466 glycine tRNA gene sequences from bacteria, archaea and eukaryote cytoplasm analyzed, CCA occurs immediately upstream of the canonical intron insertion position, suggesting the first anticodon (NCC for glycine) has been captured from the 3'-terminal CCA of one of the interacting hairpins as a result of an ancestral ligation.

Conclusion

That this imprint (including the second and third nucleotides of the glycine tRNA anticodon) has been retained through billions of years of evolution suggests Crick's 'frozen accident' hypothesis has validity for at least this very first step at the dawn of the genetic code.

Reviewers

This article was reviewed by Dr Eugene V. Koonin, Dr Rob Knight and Dr David H Ardell.     

Evolutionary deviations from the universal genetic code in ciliates

The review surveys the information, including the most recent data, on the evolution of genetic code in ciliates, which is among the few codes deviating from the universal one. We discuss the cases of recurrent, independently arising deviations from the assignments of standard codons of polypeptide chain termination in the mitochondrial and nuclear genomes of ciliates and some other protozoans. Possible molecular mechanisms are considered, which underlie deviations from standard termination code to coding glutamine (codon UAA and UAG) and cystein or tryptophane (codon UGA) in the nuclear genome. Critical analysis of the main hypotheses on the evolution of secondary deviations from the universal code in ciliates is presented.     

The translation of DNA primary base sequence into three-dimensional structure

A procedure is outlined to obtain a reliable computer-generatedrepresentation of the DNA duplex from its primary sequence ofbase pairs. The calculations are based on the potential energiesof interaction of adjacent side groups. The methods are, however,completely general and can be adapted to any set of base sequencedependent conformational rules. Static representations of theDNA are compared with the distributions of conformations obtainedfrom Monte Carlo simulation studies. Direct matrix generatorcalculations of the average (equilibrium) extension and orientationof various sequences and numerical estimates of the flexibilityof the chains as a whole are also reported. The methods areapplied to three short fragments of kinetoplast DNA from Crithidiafasciculata which exhibit dramatically different behavior onnon-denaturing poly-acrylamide gels. Received on August 17, 1987; accepted on December 19, 1987     

The fidelity of the translation of the genetic code.

Aminoacyl-tRNA synthetases play a central role in maintaining accuracy during the translation of the genetic code. To achieve this challenging task they have to discriminate against amino acids that are very closely related not only in structure but also in chemical nature. A 'double-sieve' editing model was proposed in the late seventies to explain how two closely related amino acids may be discriminated. However, a clear understanding of this mechanism required structural information on synthetases that are faced with such a problem of amino acid discrimination. The first structural basis for the editing model came recently from the crystal structure of isoleucyl-tRNA synthetase, a class I synthetase, which has to discriminate against valine. The structure showed the presence of two catalytic sites in the same enzyme, one for activation, a coarse sieve which binds both isoleucine and valine, and another for editing, a fine sieve which binds only valine and rejects isoleucine. Another structure of the enzyme in complex with tRNA showed that the tRNA is responsible for the translocation of the misactivated amino-acid substrate from the catalytic site to the editing site. These studies were mainly focused on class I synthetases and the situation was not clear about how class II enzymes discriminate against similar amino acids. The recent structural and enzymatic studies on threonyl-tRNA synthetase, a class II enzyme, reveal how this challenging task is achieved by using a unique zinc ion in the active site as well as by employing a separate domain for specific editing activity. These studies led us to propose a model which emphasizes the mirror symmetrical approach of the two classes of enzymes and highlights that tRNA is the key player in the evolution of these class of enzymes.     

On origin of genetic code and tRNA before translation

Background  

Synthesis of proteins is based on the genetic code - a nearly universal assignment of codons to amino acids (aas). A major challenge to the understanding of the origins of this assignment is the archetypal "key-lock vs. frozen accident" dilemma. Here we re-examine this dilemma in light of 1) the fundamental veto on "foresight evolution", 2) modular structures of tRNAs and aminoacyl-tRNA synthetases, and 3) the updated library of aa-binding sites in RNA aptamers successfully selected in vitro for eight amino acids.     

Evolution of the genetic code: partial optimization of a random code for robustness to translation error in a rugged fitness landscape

Background  

The standard genetic code table has a distinctly non-random structure, with similar amino acids often encoded by codons series that differ by a single nucleotide substitution, typically, in the third or the first position of the codon. It has been repeatedly argued that this structure of the code results from selective optimization for robustness to translation errors such that translational misreading has the minimal adverse effect. Indeed, it has been shown in several studies that the standard code is more robust than a substantial majority of random codes. However, it remains unclear how much evolution the standard code underwent, what is the level of optimization, and what is the likely starting point.     

No accident: genetic codes freeze in error-correcting patterns of the standard genetic code

The standard genetic code poses a challenge in understanding the evolution of information processing at a fundamental level of biological organization. Genetic codes are generally coadapted with, or 'frozen' by, the protein-coding genes that they translate, and so cannot easily change by natural selection. Yet the standard code has a significantly non-random pattern that corrects common errors in the transmission of information in protein-coding genes. Because of the freezing effect and for other reasons, this pattern has been proposed not to be due to selection but rather to be incidental to other evolutionary forces or even entirely accidental. We present results from a deterministic population genetic model of code-message coevolution. We explicitly represent the freezing effect of genes on genetic codes and the perturbative effect of changes in genetic codes on genes. We incorporate characteristic patterns of mutation and translational error, namely, transition bias and positional asymmetry, respectively. Repeated selection over small successive changes produces genetic codes that are substantially, but not optimally, error correcting. In particular, our model reproduces the error-correcting patterns of the standard genetic code. Aspects of our model and results may be applicable to the general problem of adaptation to error in other natural information-processing systems.     

The evolutionary change of the genetic code as restricted by the anticodon and identity of transfer RNA

The discovery of non-universal genetic codes in several mitochondria and nuclear systems during the past ten years has necessitated a reconsideration of the concept that the genetic code is universal and frozen, as was once believed. Here, the flexibility of the relationship between codons and amino acids is discussed on the basis of the distribution of non-universal genetic codes in various organisms insofar as has been observed to date. Judging from the result of recent investigations into tRNA identity, it would appear that the non-participation of the anticodon in recognition by aminoacyl-tRNA synthetase has significantly influenced the variability of codons.     

Block structure and stability of the genetic code

It is known that different codons may be unified into larger groups related to the hierarchical structure, approximate hidden symmetries, and evolutionary origin of the universal genetic code. Using a simplified evolutionary motivated two-letter version of genetic code, the general principles of the most stable coding are discussed. By the complete enumeration in such a reduced code it is strictly proved that the maximum stability with respect to point mutations and shifts in the reading frame needs the fixation of the middle letters within codons in groups with different physico-chemical properties, thus, explaining a key feature of the universal genetic code. The translational stability of the genetic code is studied by the mapping of code onto de Bruijn graph providing both the compact visual representation of mutual relationships between different codons as well as between codons and protein coding DNA sequence and a powerful tool for the investigation of stability of protein coding. Then, the results are extended to four-letter codes. As is shown, the universal genetic code obeys mainly the principles of optimal coding. These results demonstrate the hierarchical character of optimization of universal genetic code with strictly optimal coding being evolved at the earliest stages of molecular evolution. Finally, the universal genetic code is compared with the other natural variants of genetic codes.     

A harmonic structure of the genetic code

In this paper is presented a new, very harmonic structure of the genetic code (GC) within a system of "4 x 5" (and/or of "5 x 4") of amino acids (AAs) in two variants. In first variant, the five rows within the system start with one polar charged amino acid (AA) each, making first column, consisting from five polar charged AAs (D, R, K, H, E). Five polar non-charged AAs (N, P, Y, W, Q) follow, then five non-polar AAs as last column (A, L, F, V, I) and, finally, five polar or non-polar AAs, in a combination, as first to last column (A as non-polar; S, T as polar, and G, P as ambivalent AAs). A second variant is subsequent to this one-"4 x 5" system with five nitrogen AAs (K, R, P, H, W), five oxygen (D, E, Y, S, T), five solely carbon (A, L, F, V, I) and five "combined" AAs (G with hydrogen as side chain; C and M with carbon and sulfur; N and Q with carbon, oxygen and nitrogen). A strict balance of atom and nucleon number as well as molecule mass follows the classification in both system variants.