Evolutionary patterns in the sequence and structure of transfer RNA: a window into early translation and the genetic code

Transfer RNA (tRNA) molecules play vital roles during protein synthesis. Their acceptor arms are aminoacylated with specific amino acid residues while their anticodons delimit codon specificity. The history of these two functions has been generally linked in evolutionary studies of the genetic code. However, these functions could have been differentially recruited as evolutionary signatures were left embedded in tRNA molecules. Here we built phylogenies derived from the sequence and structure of tRNA, we forced taxa into monophyletic groups using constraint analyses, tested competing evolutionary hypotheses, and generated timelines of amino acid charging and codon discovery. Charging of Sec, Tyr, Ser and Leu appeared ancient, while specificities related to Asn, Met, and Arg were derived. The timelines also uncovered an early role of the second and then first codon bases, identified codons for Ala and Pro as the most ancient, and revealed important evolutionary take-overs related to the loss of the long variable arm in tRNA. The lack of correlation between ancestries of amino acid charging and encoding indicated that the separate discoveries of these functions reflected independent histories of recruitment. These histories were probably curbed by co-options and important take-overs during early diversification of the living world.     

The Origin and Evolution of tRNA Inferred from Phylogenetic Analysis of Structure

The evolutionary history of the two structural and functional domains of tRNA is controversial but harbors the secrets of early translation and the genetic code. To explore the origin and evolution of tRNA, we reconstructed phylogenetic trees directly from molecular structure. Forty-two structural characters describing the geometry of 571 tRNAs and three statistical parameters describing thermodynamic and mechanical features of molecules quantitatively were used to derive phylogenetic trees of molecules and molecular substructures. Trees of molecules failed to group tRNA according to amino acid specificity and did not reveal the tripartite nature of life, probably due to loss of phylogenetic signal or because tRNA diversification predated organismal diversification. Trees of substructures derived from both structural and statistical characters support the origin of tRNA in the acceptor arm and the hypothesis that the top half domain composed of acceptor and pseudouridine (TΨC) arms is more ancient than the bottom half domain composed of dihydrouridine (DHU) and anticodon arms. This constitutes the cornerstone of the genomic tag hypothesis that postulates tRNAs were ancient telomeres in the RNA world. The trees of substructures suggest a model for the evolution of the major functional and structural components of tRNA. In this model, short RNA hairpins with stems homologous to the acceptor arm of present day tRNAs were extended with regions homologous to TΨC and anticodon arms. The DHU arm was then incorporated into the resulting three-stemmed structure to form a proto-cloverleaf structure. The variable region was the last structural addition to the molecular repertoire of evolving tRNA substructures. Electronic Supplementary Material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Stable tRNA-based phylogenies using only 76 nucleotides

tRNAs are among the most ancient, highly conserved sequences on earth, but are often thought to be poor phylogenetic markers because they are short, often subject to horizontal gene transfer, and easily change specificity. Here we use an algorithm now commonly used in microbial ecology, UniFrac, to cluster 175 genomes spanning all three domains of life based on the phylogenetic relationships among their complete tRNA pools. We find that the overall pattern of similarities and differences in the tRNA pools recaptures universal phylogeny to a remarkable extent, and that the resulting tree is similar to the distribution of bootstrapped rRNA trees from the same genomes. In contrast, the trees derived from tRNAs of identical specificity or of individual isoacceptors generally produced trees of lower quality. However, some tRNA isoacceptors were very good predictors of the overall pattern of organismal evolution. These results show that UniFrac can extract meaningful biological patterns from even phylogenies with high level of statistical inaccuracy and horizontal gene transfer, and that, overall, the pattern of tRNA evolution tracks universal phylogeny and provides a background against which we can test hypotheses about the evolution of individual isoacceptors.     

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.     

Origin and Evolution of Protein Fold Designs Inferred from Phylogenomic Analysis of CATH Domain Structures in Proteomes

The spatial arrangements of secondary structures in proteins, irrespective of their connectivity, depict the overall shape and organization of protein domains. These features have been used in the CATH and SCOP classifications to hierarchically partition fold space and define the architectural make up of proteins. Here we use phylogenomic methods and a census of CATH structures in hundreds of genomes to study the origin and diversification of protein architectures (A) and their associated topologies (T) and superfamilies (H). Phylogenies that describe the evolution of domain structures and proteomes were reconstructed from the structural census and used to generate timelines of domain discovery. Phylogenies of CATH domains at T and H levels of structural abstraction and associated chronologies revealed patterns of reductive evolution, the early rise of Archaea, three epochs in the evolution of the protein world, and patterns of structural sharing between superkingdoms. Phylogenies of proteomes confirmed the early appearance of Archaea. While these findings are in agreement with previous phylogenomic studies based on the SCOP classification, phylogenies unveiled sharing patterns between Archaea and Eukarya that are recent and can explain the canonical bacterial rooting typically recovered from sequence analysis. Phylogenies of CATH domains at A level uncovered general patterns of architectural origin and diversification. The tree of A structures showed that ancient structural designs such as the 3-layer (αβα) sandwich (3.40) or the orthogonal bundle (1.10) are comparatively simpler in their makeup and are involved in basic cellular functions. In contrast, modern structural designs such as prisms, propellers, 2-solenoid, super-roll, clam, trefoil and box are not widely distributed and were probably adopted to perform specialized functions. Our timelines therefore uncover a universal tendency towards protein structural complexity that is remarkable.     

Gene trees, species trees, and morphology converge on a similar phylogeny of living gars (Actinopterygii: Holostei: Lepisosteidae), an ancient clade of ray-finned fishes

Extant gars represent the remaining members of a formerly diverse assemblage of ancient ray-finned fishes and have been the subject of multiple phylogenetic analyses using morphological data. Here, we present the first hypothesis of phylogenetic relationships among living gar species based on molecular data, through the examination of gene tree heterogeneity and coalescent species tree analyses of a portion of one mitochondrial (COI) and seven nuclear (ENC1, myh6, plagl2, S7 ribosomal protein intron 1, sreb2, tbr1, and zic1) genes. Individual gene trees displayed varying degrees of resolution with regards to species-level relationships, and the gene trees inferred from COI and the S7 intron were the only two that were completely resolved. Coalescent species tree analyses of nuclear genes resulted in a well-resolved and strongly supported phylogenetic tree of living gar species, for which Bayesian posterior node support was further improved by the inclusion of the mitochondrial gene. Species-level relationships among gars inferred from our molecular data set were highly congruent with previously published morphological phylogenies, with the exception of the placement of two species, Lepisosteus osseus and L. platostomus. Re-examination of the character coding used by previous authors provided partial resolution of this topological discordance, resulting in broad concordance in the phylogenies inferred from individual genes, the coalescent species tree analysis, and morphology. The completely resolved phylogeny inferred from the molecular data set with strong Bayesian posterior support at all nodes provided insights into the potential for introgressive hybridization and patterns of allopatric speciation in the evolutionary history of living gars, as well as a solid foundation for future examinations of functional diversification and evolutionary stasis in a "living fossil" lineage.     

Global phylogeny determined by the combination of protein domains in proteomes

The majority of proteins consist of multiple domains that are either repeated or combined in defined order. In this study, we survey the combination of protein domains defined at fold and fold superfamily levels in 185 genomes belonging to organisms that have been fully sequenced and introduce a method that reconstructs rooted phylogenomic trees from the content and arrangement of domains in proteins at a genomic level. We find that the majority of domain combinations were unique to Archaea, Bacteria, or Eukarya, suggesting most combinations originated after life had diversified. Domain repeat and domain repeat within multidomain proteins increased notably in eukaryotes, mainly at the expense of single-domain and domain-pair proteins. This increase was mostly confined to Metazoa. We also find an unbalanced sharing of domain combinations which suggests that Eukarya is more closely related to Bacteria than to Archaea, an observation that challenges the widely assumed eukaryote-archaebacterial sisterhood relationship. The occurrence and abundance of the molecular repertoire (interactome) of domain combinations was used to generate phylogenomic trees. These global interactome-based phylogenies described organismal histories satisfactorily, revealing the tripartite nature of life, and supporting controversial evolutionary patterns, such as the Coelomata hypothesis, the grouping of plants and animals, and the Gram-positive origin of bacteria. Results suggest strongly that the process of domain combination is not random but curved by evolution, rejecting the null hypothesis of domain modules combining in the absence of natural selection or an optimality criterion.     

About the last common ancestor, the universal life-tree and lateral gene transfer: a reappraisal

An organismal tree rooted in the bacterial branch and derived from a hyperthermophilic last common ancestor (LCA) is still widely assumed to represent the path followed by evolution from the most primeval cells to the three domains recognized among contemporary organisms: Bacteria, Archaea and Eucarya. In the past few years, however, more and more discrepancies between this pattern and individual protein trees have been brought to light. There has been an overall tendency to attribute these incongruities to widespread lateral gene transfer. However, recent developments, a reappraisal of earlier evidence and considerations of our own lead us to a quite different view. It would appear (i) that the role of lateral gene transfer was overemphasized in recent discussions of molecular phylogenies; (ii) that the LCA was probably a non-thermophilic protoeukaryote from which both Archaea and Bacteria emerged by reductive evolution but not as sister groups, in keeping with a current evolutionary scheme for the biosynthesis of membrane lipids; and (iii) that thermophilic Archaea may have been the first branch to diverge from the ancestral line.     

Proteome Evolution and the Metabolic Origins of Translation and Cellular Life

The origin of life has puzzled molecular scientists for over half a century. Yet fundamental questions remain unanswered, including which came first, the metabolic machinery or the encoding nucleic acids. In this study we take a protein-centric view and explore the ancestral origins of proteins. Protein domain structures in proteomes are highly conserved and embody molecular functions and interactions that are needed for cellular and organismal processes. Here we use domain structure to study the evolution of molecular function in the protein world. Timelines describing the age and function of protein domains at fold, fold superfamily, and fold family levels of structural complexity were derived from a structural phylogenomic census in hundreds of fully sequenced genomes. These timelines unfold congruent hourglass patterns in rates of appearance of domain structures and functions, functional diversity, and hierarchical complexity, and revealed a gradual build up of protein repertoires associated with metabolism, translation and DNA, in that order. The most ancient domain architectures were hydrolase enzymes and the first translation domains had catalytic functions for the aminoacylation and the molecular switch-driven transport of RNA. Remarkably, the most ancient domains had metabolic roles, did not interact with RNA, and preceded the gradual build-up of translation. In fact, the first translation domains had also a metabolic origin and were only later followed by specialized translation machinery. Our results explain how the generation of structure in the protein world and the concurrent crystallization of translation and diversified cellular life created further opportunities for proteomic diversification.     

Universal Sharing Patterns in Proteomes and Evolution of Protein Fold Architecture and Life

Protein evolution is imprinted in both the sequence and the structure of evolutionary building blocks known as protein domains. These domains share a common ancestry and can be unified into a comparatively small set of folding architectures, the protein folds. We have traced the distribution of protein folds between and within proteomes belonging to Eukarya, Archaea, and Bacteria along the branches of a universal phylogeny of protein architecture. This tree was reconstructed from global fold-usage statistics derived from a structural census of proteomes. We found that folds shared by the three organismal domains were placed almost exclusively at the base of the rooted tree and that there were marked heterogeneities in fold distribution and clear evolutionary patterns related to protein architecture and organismal diversification. These include a relative timing for the emergence of prokaryotes, congruent episodes of architectural loss and diversification in Archaea and Bacteria, and a late and quite massive rise of architectural novelties in Eukarya perhaps linked to multicellularity.Reviewing Editor : Dr. David Pollock     

BOOM AND BUST: ANCIENT AND RECENT DIVERSIFICATION IN BICHIRS (POLYPTERIDAE: ACTINOPTERYGII), A RELICTUAL LINEAGE OF RAY‐FINNED FISHES

Understanding the history that underlies patterns of species richness across the Tree of Life requires an investigation of the mechanisms that not only generate young species‐rich clades, but also those that maintain species‐poor lineages over long stretches of evolutionary time. However, diversification dynamics that underlie ancient species‐poor lineages are often hidden due to a lack of fossil evidence. Using information from the fossil record and time calibrated molecular phylogenies, we investigate the history of lineage diversification in Polypteridae, which is the sister lineage of all other ray‐finned fishes (Actinopterygii). Despite originating at least 390 million years (Myr) ago, molecular timetrees support a Neogene origin for the living polypterid species. Our analyses demonstrate polypterids are exceptionally species depauperate with a stem lineage duration that exceeds 380 million years (Ma) and is significantly longer than the stem lineage durations observed in other ray‐finned fish lineages. Analyses of the fossil record show an early Late Cretaceous (100.5–83.6 Ma) peak in polypterid genus richness, followed by 60 Ma of low richness. The Neogene species radiation and evidence for high‐diversity intervals in the geological past suggest a “boom and bust” pattern of diversification that contrasts with common perceptions of relative evolutionary stasis in so‐called “living fossils.”     

Evolved RNA Secondary Structure and the Rooting of the Universal Tree of Life

The origin and diversification of RNA secondary structure were traced using cladistic methods. Structural components were coded as polarized and ordered multi-state characters, following a model of character state transformation outlined by considerations in statistical mechanics. Several classes of functional RNA were analyzed, including ribosomal RNA (rRNA). Considerable phylogenetic signal was present in their secondary structure. The intrinsically rooted phylogenies reconstructed from evolved RNA structure depicted those derived from nucleic acid sequence at all taxonomical levels, and grouped organisms in concordance with traditional classification, especially in the archaeal and eukaryal domains. Natural selection appears therefore to operate early in the information flow that originates in sequence and ends in an adapted phenotype. When examining the hierarchical classification of the living world, phylogenetic analysis of secondary structure of the small and large rRNA subunits reconstructed a universal tree of life that branched in three monophyletic groups corresponding to Eucarya, Archaea, and Bacteria, and was rooted in the eukaryotic branch. Ribosomal characters involved in the translational cycle could be easily traced and showed that transfer RNA (tRNA) binding domains in the large rRNA subunit evolved concurrently with the rest of the rRNA molecule. Results suggest it is equally parsimonious to consider that ancestral unicellular eukaryotes or prokaryotes gave rise to all extant life forms and provide a rare insight into the early evolution of nucleic acid and protein biosynthesis. Received: 13 September 2000 / Accepted: 27 August 2001     

Giant viruses coexisted with the cellular ancestors and represent a distinct supergroup along with superkingdoms Archaea, Bacteria and Eukarya

ABSTRACT: BACKGROUND: The discovery of giant viruses with genome and physical size comparable to cellular organisms, remnants of protein translation machinery and virus-specific parasites (virophages) have raised intriguing questions about their origin. Evidence advocates for their inclusion into global phylogenomic studies and their consideration as a distinct and ancient form of life. RESULTS: Here we reconstruct phylogenies describing the evolution of proteomes and protein domain structures of cellular organisms and double-stranded DNA viruses with medium-to-very-large proteomes (giant viruses). Trees of proteomes define viruses as a 'fourth supergroup' along with superkingdoms Archaea, Bacteria, and Eukarya. Trees of domains indicate they have evolved via massive and primordial reductive evolutionary processes. The distribution of domain structures suggests giant viruses harbor a significant number of protein domains including those with no cellular representation. The genomic and structural diversity embedded in the viral proteomes is comparable to the cellular proteomes of organisms with parasitic lifestyles. Since viral domains are widespread among cellular species, we propose that viruses mediate gene transfer between cells and crucially enhance biodiversity. CONCLUSIONS: Results call for a change in the way viruses are perceived. They likely represent a distinct form of life that either predated or coexisted with the last universal common ancestor (LUCA) and constitute a very crucial part of our planet's biosphere.     

Investigating the timing of origin and evolutionary processes shaping regional species diversity: Insights from simulated data and neotropical butterfly diversification rates

Different diversification scenarios have been proposed to explain the origin of extant biodiversity. However, most existing meta‐analyses of time‐calibrated phylogenies rely on approaches that do not quantitatively test alternative diversification processes. Here, I highlight the shortcomings of using species divergence ranks, which is a method widely used in meta‐analyses. Divergence ranks consist of categorizing cladogenetic events to certain periods of time, typically to either Pleistocene or to pre‐Pleistocene ages. This approach has been claimed to shed light on the origin of most extant species and the timing and dynamics of diversification in any biogeographical region. However, interpretations drawn from such method often confound two fundamental questions in macroevolutionary studies, tempo (timing of evolutionary rate shifts) and mode (“how” and “why” of speciation). By using simulated phylogenies under four diversification scenarios, constant‐rate, diversity‐dependence, high extinction, and high speciation rates in the Pleistocene, I showed that interpretations based on species divergence ranks might have been seriously misleading. Future meta‐analyses of dated phylogenies need to be aware of the impacts of incomplete taxonomic sampling, tree topology, and divergence time uncertainties, as well as they might be benefited by including quantitative tests of alternative diversification models that acknowledge extinction and diversity dependence.     

Evolution of the WANCY region in amniote mitochondrial DNA

In most vertebrate mitochondrial genomes, the site for initiation of light-strand replication, OL, is found within a cluster of five transfer RNA (tRNA) genes (tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), and tRNA(Tyr)). This region and part of the adjacent cytochrome c oxydase subunit I (COI) gene were sequenced for two crocodilian, two turtle, and one snake species and for Sphenodon punctatus; part of the adjacent nicotinamide adenine dinucleotide dehydrogenase subunit 2 (ND2) gene was also sequenced for the crocodilian and turtle species. All had the typical vertebrate gene order. The turtles and the snake have a lengthy noncoding sequence between the tRNA(Asn) and tRNA(Cys) genes that we assumed to be homologous to the mammalian OL. The crocodilians and Sphenodon lack such a sequence, a condition they share with birds. Most proposed phylogenies for the amniotes require that OL at this position was lost at least twice during their diversification or was evolved independently more than once. Within the five tRNA genes, frequencies of substitutions are much higher in loops than in stems. Many loops vary dramatically in size among the species; in the most extreme case, the D-arm of the Sphenodon tRNA(Cys) is a "D-arm replacement" loop of seven nucleotides. Frequency of transitions in stems is relatively uniform across tRNAs, but frequency of transversions varies greatly. Mismatches in stems are infrequent, and their relative frequency in a specific tRNA is unrelated to the frequency of substitution in the corresponding gene. Several features of mammalian mitochondrial tRNAs are conserved in WANCY tRNAs throughout amniotes. The inferred initiation codon for COI is GTG in crocodilians, turtles, and the snake, a condition they share with fishes, certain amphibians, and birds. TTG appears to be the initiation codon for COI in Sphenodon; if correct, this would be a novel initiation codon for vertebrate mitochondrial DNA. Phylogenetic analyses of the inferred amino acid sequences of ND2 and COI support the sister-group relationship of birds and crocodilians and suggest that mammals are an early derived lineage within the amniotes.      

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.     

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.