Sequence evidence in the archaeal genomes that tRNAs emerged through the combination of ancestral genes as 5' and 3' tRNA halves

Footprints of a Singular 22-Nucleotide RNA Ring at the Origin of Life

(1) Background: Previous experimental observations and theoretical hypotheses have been providing insight into a hypothetical world where an RNA hairpin or ring may have debuted as the primary informational and functional molecule. We propose a model revisiting the architecture of RNA-peptide interactions at the origin of life through the evolutionary dynamics of RNA populations. (2) Methods: By performing a step-by-step computation of the smallest possible hairpin/ring RNA sequences compatible with building up a variety of peptides of the primitive network, we inferred the sequence of a singular docosameric RNA molecule, we call the ALPHA sequence. Then, we searched for any relics of the peptides made from ALPHA in sequences deposited in the different public databases. (3) Results: Sequence matching between ALPHA and sequences from organisms among the earliest forms of life on Earth were found at high statistical relevance. We hypothesize that the frequency of appearance of relics from ALPHA sequence in present genomes has a functional necessity. (4) Conclusions: Given the fitness of ALPHA as a supportive sequence of the framework of all existing theories, and the evolution of Archaea and giant viruses, it is anticipated that the unique properties of this singular archetypal ALPHA sequence should prove useful as a model matrix for future applications, ranging from synthetic biology to DNA computing.

Why Is AUG the Start Codon?: Theoretical Minimal RNA Rings: Maximizing Coded Information Biases 1st Codon for the Universal Initiation Codon AUG

The rational design of theoretical minimal RNA rings predetermines AUG as the universal start codon. This design maximizes coded amino acid diversity over minimal sequence length, defining in silico theoretical minimal RNA rings, candidate ancestral genes. RNA rings code for 21 amino acids and a stop codon after three consecutive translation rounds, and form a degradation-delaying stem-loop hairpin. Twenty-five RNA rings match these constraints, ten start with the universal initiation codon AUG. No first codon bias exists among remaining RNA rings. RNA ring design predetermines AUG as initiation codon. This is the only explanation yet for AUG as start codon. RNA ring design determines additional RNA ring gene- and tRNA-like properties described previously, because it presumably mimics constraints on life's primordial RNAs.

The Uroboros Theory of Life's Origin: 22-Nucleotide Theoretical Minimal RNA Rings Reflect Evolution of Genetic Code and tRNA-rRNA Translation Machineries

Theoretical minimal RNA rings attempt to mimick life's primitive RNAs. At most 25 22-nucleotide-long RNA rings code once for each biotic amino acid, a start and a stop codon and form a stem-loop hairpin, resembling consensus tRNAs. We calculated, for each RNA ring's 22 potential splicing positions, similarities of predicted secondary structures with tRNA vs. rRNA secondary structures. Assuming rRNAs partly derived from tRNA accretions, we predict positive associations between relative secondary structure similarities with rRNAs over tRNAs and genetic code integration orders of RNA ring anticodon cognate amino acids. Analyses consider for each secondary structure all nucleotide triplets as potential anticodon. Anticodons for ancient, chemically inert cognate amino acids are most frequent in the 25 RNA rings. For RNA rings with primordial cognate amino acids according to tRNA-homology-derived anticodons, tRNA-homology and coding sequences coincide, these are separate for predicted cognate amino acids that presumably integrated late the genetic code. RNA ring secondary structure similarity with rRNA over tRNA secondary structures associates best with genetic code integration orders of anticodon cognate amino acids when assuming split anticodons (one and two nucleotides at the spliced RNA ring 5' and 3' extremities, respectively), and at predicted anticodon location in the spliced RNA ring's midst. Results confirm RNA ring homologies with tRNAs and CDs, ancestral status of tRNA half genes split at anticodons, the tRNA-rRNA axis of RNA evolution, and that single theoretical minimal RNA rings potentially produce near-complete proto-tRNA sets. Hence genetic code pre-existence determines 25 short circular gene- and tRNA-like RNAs. Accounting for each potential splicing position, each RNA ring potentially translates most amino acids, realistically mimicks evolution of the tRNA-rRNA translation machinery. These RNA rings 'of creation' remind the uroboros' (snake biting its tail) symbolism for creative regeneration.

tRNA diversification among uncultured archeon clones

Whole genome sequences (DNA sequences) of four uncultured archeon clones (1B6:CR626858.1, 4B7:CR626856.1, 22i07:JQ768096.1 and 19c08:JQ768095.1) were collected from NCBI BioSample database for the construction of digital data on tRNA. tRNAscan-SE 2.0 and ENDMEMO tools were used to identify and sketch tRNA structure as well as calculate Guanine-Cytosine (GC) percentage respectively. Eight true/functional tRNAs were identified from above 4 sequences which showed cove score greater than 20% with no variable loop. The tRNAs from the uncultured archeon clones were classified as Ala, Arg, Ile, Thr, Pro and Val type tRNA with cove score ranging from 34.22%-79.03%. The range of GC content was found 42.89%-56.91%; while tRNA contributed GC content ranging from 52%-64.86% to the total GC content in these sequences. The data fabricated in this study could be very useful for studying the diversity of tRNA among prokaryotes.

Evolutionary Analysis of HIV-1 Pol Proteins Reveals Representative Residues for Viral Subtype Differentiation

RNA viruses have been used as model systems to understand the patterns and processes of molecular evolution because they have high mutation rates and are genetically diverse. Human immunodeficiency virus 1 (HIV-1), the etiological agent of acquired immune deficiency syndrome, is highly genetically diverse, and is classified into several groups and subtypes. However, it has been difficult to use its diverse sequences to establish the overall phylogenetic relationships of different strains or the trends in sequence conservation with the construction of phylogenetic trees. Our aims were to systematically characterize HIV-1 subtype evolution and to identify the regions responsible for HIV-1 subtype differentiation at the amino acid level in the Pol protein, which is often used to classify the HIV-1 subtypes. In this study, we systematically characterized the mutation sites in 2,052 Pol proteins from HIV-1 group M (144 subtype A; 1,528 subtype B; 380 subtype C), using sequence similarity networks. We also used spectral clustering to group the sequences based on the network graph structures. A stepwise analysis of the cluster hierarchies allowed us to estimate a possible evolutionary pathway for the Pol proteins. The subtype A sequences also clustered according to when and where the viruses were isolated, whereas both the subtype B and C sequences remained as single clusters. Because the Pol protein has several functional domains, we identified the regions that are discriminative by comparing the structures of the domain-based networks. Our results suggest that sequence changes in the RNase H domain and the reverse transcriptase (RT) connection domain are responsible for the subtype classification. By analyzing the different amino acid compositions at each site in both domain sequences, we found that a few specific amino acid residues (i.e., M357 in the RT connection domain and Q480, Y483, and L491 in the RNase H domain) represent the differences among the subtypes. These residues were located on the surface of the RT structure and in the vicinity of the amino acid sites responsible for RT enzymatic activity or function.

Intercompartmental Piecewise Gene Transfer

Gene relocation from the residual genomes of organelles to the nuclear genome still continues, although as a scaled down evolutionary phenomenon, limited in occurrence mostly to protists (sensu lato) and land plants. During this process, the structural integrity of transferred genes is usually preserved. However, the relocation of mitochondrial genes that code for respiratory chain and ribosomal proteins is sometimes associated with their fragmentation into two complementary genes. Herein, this review compiles cases of piecewise gene transfer from the mitochondria to the nucleus, and discusses hypothesized mechanistic links between the fission and relocation of those genes.

Disrupted tRNA Genes and tRNA Fragments: A Perspective on tRNA Gene Evolution

Transfer RNAs (tRNAs) are small non-coding RNAs with lengths of approximately 70-100 nt. They are directly involved in protein synthesis by carrying amino acids to the ribosome. In this sense, tRNAs are key molecules that connect the RNA world and the protein world. Thus, study of the evolution of tRNA molecules may reveal the processes that led to the establishment of the central dogma: genetic information flows from DNA to RNA to protein. Thanks to the development of DNA sequencers in this century, we have determined a huge number of nucleotide sequences from complete genomes as well as from transcriptomes in many species. Recent analyses of these large data sets have shown that particular tRNA genes, especially in Archaea, are disrupted in unique ways: some tRNA genes contain multiple introns and some are split genes. Even tRNA molecules themselves are fragmented post-transcriptionally in many species. These fragmented small RNAs are known as tRNA-derived fragments (tRFs). In this review, I summarize the progress of research into the disrupted tRNA genes and the tRFs, and propose a possible model for the molecular evolution of tRNAs based on the concept of the combination of fragmented tRNA halves.

Handling tRNA introns, archaeal way and eukaryotic way

Introns are found in various tRNA genes in all the three kingdoms of life. Especially, archaeal and eukaryotic genomes are good sources of tRNA introns that are removed by proteinaceous splicing machinery. Most intron-containing tRNA genes both in archaea and eukaryotes possess an intron at a so-called canonical position, one nucleotide 3′ to their anticodon, while recent bioinformatics have revealed unusual types of tRNA introns and their derivatives especially in archaeal genomes. Gain and loss of tRNA introns during various stages of evolution are obvious both in archaea and eukaryotes from analyses of comparative genomics. The splicing of tRNA molecules has been studied extensively from biochemical and cell biological points of view, and such analyses of eukaryotic systems provided interesting findings in the past years. Here, I summarize recent progresses in the analyses of tRNA introns and the splicing process, and try to clarify new and old questions to be solved in the next stages.

tRNA gene diversity in the three domains of life

Transfer RNA (tRNA) is widely known for its key role in decoding mRNA into protein. Despite their necessity and relatively short nucleotide sequences, a large diversity of gene structures and RNA secondary structures of pre-tRNAs and mature tRNAs have recently been discovered in the three domains of life. Growing evidences of disrupted tRNA genes in the genomes of Archaea reveals unique gene structures such as, intron-containing tRNA, split tRNA, and permuted tRNA. Coding sequence for these tRNAs are either separated with introns, fragmented, or permuted at the genome level. Although evolutionary scenario behind the tRNA gene disruption is still unclear, diversity of tRNA structure seems to be co-evolved with their processing enzyme, so-called RNA splicing endonuclease. Metazoan mitochondrial tRNAs (mtRNAs) are known for their unique lack of either one or two arms from the typical tRNA cloverleaf structure, while still maintaining functionality. Recently identified nematode-specific V-arm containing tRNAs (nev-tRNAs) possess long variable arms that are specific to eukaryotic class II tRNA^Ser and tRNA^Leu but also decode class I tRNA codons. Moreover, many tRNA-like sequences have been found in the genomes of different organisms and viruses. Thus, this review is aimed to cover the latest knowledge on tRNA gene diversity and further recapitulate the evolutionary and biological aspects that caused such uniqueness.

Circularly permuted tRNA genes: their expression and implications for their physiological relevance and development

A number of genome analyses and searches using programs that focus on the RNA-specific bulge-helix-bulge (BHB) motif have uncovered a wide variety of disrupted tRNA genes. The results of these analyses have shown that genetic information encoding functional RNAs is described in the genome cryptically and is retrieved using various strategies. One such strategy is represented by circularly permuted tRNA genes, in which the sequences encoding the 5′-half and 3′-half of the specific tRNA are separated and inverted on the genome. Biochemical analyses have defined a processing pathway in which the termini of tRNA precursors (pre-tRNAs) are ligated to form a characteristic circular RNA intermediate, which is then cleaved at the acceptor-stem to generate the typical cloverleaf structure with functional termini. The sequences adjacent to the processing site located between the 3′-half and the 5′-half of pre-tRNAs potentially form a BHB motif, which is the dominant recognition site for the tRNA-intron splicing endonuclease, suggesting that circularization of pre-tRNAs depends on the splicing machinery. Some permuted tRNAs contain a BHB-mediated intron in their 5′- or 3′-half, meaning that removal of an intron, as well as swapping of the 5′- and 3′-halves, are required during maturation of their pre-tRNAs. To date, 34 permuted tRNA genes have been identified from six species of unicellular algae and one archaeon. Although their physiological significance and mechanism of development remain unclear, the splicing system of BHB motifs seems to have played a key role in the formation of permuted tRNA genes. In this review, current knowledge of circularly permuted tRNA genes is presented and some unanswered questions regarding these species are discussed.

Putative anticodons in mitochondrial tRNA sidearm loops: Pocketknife tRNAs?

The hypothesis that tRNA sidearm loops bear anticodons assumes crossovers between anticodon and sidearms, or translation by expressed aminoacylated tRNA halves forming single stem-loops. Only the latter might require ribosomal adaptations. Drosophila mitochondrial codon usages coevolve with sidearm numbers bearing matching putative anticodons (comparing different codon families in one genome, macroevolution) and when comparing different genomes for single codon families (microevolution). Coevolution between Drosophila and yeast mitochondrial antisense tRNAs and codon usages partly confounds microevolutionary patterns for putative sidearm anticodons. Some tRNA sidearm loops have more than seven nucleotides, putative expanded anticodons potentially matching quadruplet codons (tetracodons, codons expanded by a fourth silent position, forming tetragenes (predicted by alignment analyses of Drosophila mitochondrial genomes)). Tetracodon numbers coevolve with expanded tRNA sidearm loops. Sidearm coevolution with amino acid usages and tetragenes occurs for putative anticodons in 5' and 3' sidearms loops (D and TΨC loops, respectively), are stronger for the D-loop. Results slightly favour isolated stem-loops upon crossover hypotheses. An alternative hypothesis, that patterns observed for sidearm 'anticodons' do not imply translational activity, but recognition signals for tRNA synthetases that aminoacylate tRNAs, is incompatible with tetracodon/tetra-anticodon coevolution. Hence analyses strengthen translational hypotheses for tRNA sidearm anticodons, tetragenes, and antisense tRNAs.

Genome-wide analysis reveals origin of transfer RNA genes from tRNA halves

Transfer RNAs (tRNAs) play an important role linking mitochondrial RNA and amino acids during protein biogenesis. Four types of tRNA genes have been identified in living organisms. However, the evolutionary origin of tRNAs remains largely unknown. In this article, we conduct a deep sequence analysis of diverse genomes that cover all three domains of life to unveil the evolutionary history of tRNA genes from tRNA halves. tRNA half homologs were detected in diverse organisms, and some of them were expressed in mouse tissues. Continuous tRNA genes have a conserved pattern similar to indels, which is, more closely flanking regions have higher single nucleotide substitution rates, whereas tRNA half homologs do not have this pattern. In addition, tRNAs tend to break into tRNA halves when tissues are incubated in vitro, the tendency of tRNA to break into tRNA halves may be a "side-effect" of tRNA genes evolving from tRNA halves. These results suggest that modern tRNAs originated from tRNA halves through a repeat element-mediated mechanism. These findings provide insight into the evolutionary origin of tRNA genes.

Does the genetic code have a eukaryotic origin?

In the RNA world, RNA is assumed to be the dominant macromolecule performing most, if not all, core “house-keeping” functions. The ribo-cell hypothesis suggests that the genetic code and the translation machinery may both be born of the RNA world, and the introduction of DNA to ribo-cells may take over the informational role of RNA gradually, such as a mature set of genetic code and mechanism enabling stable inheritance of sequence and its variation. In this context, we modeled the genetic code in two content variables—GC and purine contents—of protein-coding sequences and measured the purine content sensitivities for each codon when the sensitivity (% usage) is plotted as a function of GC content variation. The analysis leads to a new pattern—the symmetric pattern—where the sensitivity of purine content variation shows diagonally symmetry in the codon table more significantly in the two GC content invariable quarters in addition to the two existing patterns where the table is divided into either four GC content sensitivity quarters or two amino acid diversity halves. The most insensitive codon sets are GUN (valine) and CAN (CAR for asparagine and CAY for aspartic acid) and the most biased amino acid is valine (always over-estimated) followed by alanine (always under-estimated). The unique position of valine and its codons suggests its key roles in the final recruitment of the complete codon set of the canonical table. The distinct choice may only be attributable to sequence signatures or signals of splice sites for spliceosomal introns shared by all extant eukaryotes.

Identification and functional characterization of tRNA-derived RNA fragments (tRFs) in respiratory syncytial virus infection

The discovery of small noncoding RNAs (sncRNAs) with regulatory functions is a recent breakthrough in biology. Among sncRNAs, microRNA (miRNA), derived from host or virus, has emerged as elements with high importance in control of viral replication and host responses. However, the expression pattern and functional aspects of other types of sncRNAs, following viral infection, are unexplored. In order to define expression patterns of sncRNAs, as well as to discover novel regulatory sncRNAs in response to viral infection, we applied deep sequencing to cells infected with human respiratory syncytial virus (RSV), the most common cause of bronchiolitis and pneumonia in babies. RSV infection leads to abundant production of transfer RNA (tRNA)-derived RNA Fragments (tRFs) that are ~30 nucleotides (nts) long and correspond to the 5'-half of mature tRNAs. At least one tRF, which is derived from tRNA-Glu-CTC, represses target mRNA in the cytoplasm and promotes RSV replication. This demonstrates that this tRF is not a random by-product of tRNA degradation but a functional molecule. The biogenesis of this tRF is also specific, as it is mediated by the endonuclease angiogenin (ANG), not by other nucleases. In summary, our study presents novel information on the induction of a functional tRF by viral infection.

The RNA world hypothesis: the worst theory of the early evolution of life (except for all the others)(a)

The problems associated with the RNA world hypothesis are well known. In the following I discuss some of these difficulties, some of the alternative hypotheses that have been proposed, and some of the problems with these alternative models. From a biosynthetic - as well as, arguably, evolutionary - perspective, DNA is a modified RNA, and so the chicken-and-egg dilemma of "which came first?" boils down to a choice between RNA and protein. This is not just a question of cause and effect, but also one of statistical likelihood, as the chance of two such different types of macromolecule arising simultaneously would appear unlikely. The RNA world hypothesis is an example of a 'top down' (or should it be 'present back'?) approach to early evolution: how can we simplify modern biological systems to give a plausible evolutionary pathway that preserves continuity of function? The discovery that RNA possesses catalytic ability provides a potential solution: a single macromolecule could have originally carried out both replication and catalysis. RNA - which constitutes the genome of RNA viruses, and catalyzes peptide synthesis on the ribosome - could have been both the chicken and the egg! However, the following objections have been raised to the RNA world hypothesis: (i) RNA is too complex a molecule to have arisen prebiotically; (ii) RNA is inherently unstable; (iii) catalysis is a relatively rare property of long RNA sequences only; and (iv) the catalytic repertoire of RNA is too limited. I will offer some possible responses to these objections in the light of work by our and other labs. Finally, I will critically discuss an alternative theory to the RNA world hypothesis known as 'proteins first', which holds that proteins either preceded RNA in evolution, or - at the very least - that proteins and RNA coevolved. I will argue that, while theoretically possible, such a hypothesis is probably unprovable, and that the RNA world hypothesis, although far from perfect or complete, is the best we currently have to help understand the backstory to contemporary biology.

Genomic heterogeneity in a natural archaeal population suggests a model of tRNA gene disruption

Understanding the mechanistic basis of the disruption of tRNA genes, as manifested in the intron-containing and split tRNAs found in Archaea, will provide considerable insight into the evolution of the tRNA molecule. However, the evolutionary processes underlying these disruptions have not yet been identified. Previously, a composite genome of the deep-branching archaeon Caldiarchaeum subterraneum was reconstructed from a community genomic library prepared from a C. subterraneum–dominated microbial mat. Here, exploration of tRNA genes from the library reveals that there are at least three types of heterogeneity at the tRNA^Thr(GGU) gene locus in the Caldiarchaeum population. All three involve intronic gain and splitting of the tRNA gene. Of two fosmid clones found that encode tRNA^Thr(GGU), one (tRNA^Thr-I) contains a single intron, whereas another (tRNA^Thr-II) contains two introns. Notably, in the clone possessing tRNA^Thr-II, a 5′ fragment of the tRNA^Thr-I (tRNA^Thr-F) gene was observed 1.8-kb upstream of tRNA^Thr-II. The composite genome contains both tRNA^Thr-II and tRNA^Thr-F, although the loci are >500 kb apart. Given that the 1.8-kb sequence flanked by tRNA^Thr-F and tRNA^Thr-II is predicted to encode a DNA recombinase and occurs in six regions of the composite genome, it may be a transposable element. Furthermore, its dinucleotide composition is most similar to that of the pNOB8-type plasmid, which is known to integrate into archaeal tRNA genes. Based on these results, we propose that the gain of the tRNA intron and the scattering of the tRNA fragment occurred within a short time frame via the integration and recombination of a mobile genetic element.

Mapping hidden potential identity elements by computing the average discriminating power of individual tRNA positions

The recently published discrete mathematical method, extended consensus partition (ECP), identifies nucleotide types at each position that are strictly absent from a given sequence set, while occur in other sets. These are defined as discriminating elements (DEs). In this study using the ECP approach, we mapped potential hidden identity elements that discriminate the 20 different tRNA identities. We filtered the tDNA data set for the obligatory presence of well-established tRNA features, and then separately for each identity set, the presence of already experimentally identified strictly present identity elements. The analysis was performed on the three kingdoms of life. We determined the number of DE, e.g. the number of sets discriminated by the given position, for each tRNA position of each tRNA identity set. Then, from the positional DE numbers obtained from the 380 pairwise comparisons of the 20 identity sets, we calculated the average excluding value (AEV) for each tRNA position. The AEV provides a measure on the overall discriminating power of each position. Using a statistical analysis, we show that positional AEVs correlate with the number of already identified identity elements. Positions having high AEV but lacking published identity elements predict hitherto undiscovered tRNA identity elements.

A model of proto-anti-codon RNA enzymes requiring L-amino acid homochirality

All living organisms encode the 20 natural amino acid units of polypeptides using a universal scheme of triplet nucleotide "codons". Disparate features of this codon scheme are potentially informative of early molecular evolution: (i) the absence of any codons for D-amino acids; (ii) the odd combination of alternate codon patterns for some amino acids; (iii) the confinement of synonymous positions to a codon's third nucleotide; (iv) the use of 20 specific amino acids rather than a number closer to the full coding potential of 64; and (v) the evolutionary relationship of patterns in stop codons to amino acid codons. Here I propose a model for an ancestral proto-anti-codon RNA (pacRNA) auto-aminoacylation system and show that pacRNAs would naturally manifest features of the codon table. I show that pacRNAs could implement all the steps for auto-aminoacylation: amino acid coordination, intermediate activation of the amino acid by the 5'-end of the pacRNA, and 3'-aminoacylation of the pacRNA. The anti-codon cradles of pacRNAs would have been able to recognize and coordinate only a small number of L-amino acids via hydrogen bonding. A need for proper spatial coordination would have limited the number of chargeable amino acids for all anti-codon sequences, in addition to making some anti-codon sequences unsuitable. Thus, the pacRNA model implies that the idiosyncrasies of the anti-codon table and L-amino acid homochirality co-evolved during a single evolutionary period. These results further imply that early life consisted of an aminoacylated RNA world with a richer enzymatic potential than ribonucleotides alone.

Coadaptation of isoacceptor tRNA genes and codon usage bias for translation efficiency in Aedes aegypti and Anopheles gambiae

The transfer RNAs (tRNAs) are essential components of translational machinery. We determined that tRNA isoacceptors (tRNAs with different anticodons but incorporating the same amino acid in protein synthesis) show differential copy number abundance, genomic distribution patterns and sequence evolution between Aedes aegypti and Anopheles gambiae mosquitoes. The tRNA-Ala genes are present in unusually high copy number in the Ae. aegypti genome but not in An. gambiae. Many of the tRNA-Ala genes of Ae. aegypti are flanked by a highly conserved sequence that is not observed in An. gambiae. The relative abundance of tRNA isoacceptor genes is correlated with preferred (or optimal) and nonpreferred (or rare) codons for ∼2-4% of the predicted protein coding genes in both species. The majority (∼74-85%) of these genes are related to pathways involved with translation, energy metabolism and carbohydrate metabolism. Our results suggest that these genes and the related pathways may be under translational selection in these mosquitoes.

Large-scale tRNA intron transposition in the archaeal order Thermoproteales represents a novel mechanism of intron gain

Recently, diverse arrangements of transfer RNA (tRNA) genes have been found in the domain Archaea, in which the tRNA is interrupted by a maximum of three introns or is even fragmented into two or three genes. Whereas most of the eukaryotic tRNA introns are inserted strictly at the canonical nucleotide position (37/38), archaeal intron-containing tRNAs have a wide diversity of small tRNA introns, which differ in their numbers and locations. This feature is especially pronounced in the archaeal order Thermoproteales. In this study, we performed a comprehensive sequence comparison of 286 tRNA introns and their genes in seven Thermoproteales species to clarify how these introns have emerged and diversified during tRNA gene evolution. We identified 46 intron groups containing sets of highly similar sequences (>70%) and showed that 16 of them contain sequences from evolutionarily distinct tRNA genes. The phylogeny of these 16 intron groups indicates that transposition events have occurred at least seven times throughout the evolution of Thermoproteales. These findings suggest that frequent intron transposition occurs among the tRNA genes of Thermoproteales. Further computational analysis revealed limited insertion positions and corresponding amino acid types of tRNA genes. This has arisen because the bulge-helix-bulge splicing motif is required at the newly transposed position if the pre-tRNA is to be correctly processed. These results clearly demonstrate a newly identified mechanism that facilitates the late gain of short introns at various noncanonical positions in archaeal tRNAs.

Genomic organization of eukaryotic tRNAs

BACKGROUND

Surprisingly little is known about the organization and distribution of tRNA genes and tRNA-related sequences on a genome-wide scale. While tRNA gene complements are usually reported in passing as part of genome annotation efforts, and peculiar features such as the tandem arrangements of tRNA gene in Entamoeba histolytica have been described in some detail, systematic comparative studies are rare and mostly restricted to bacteria. We therefore set out to survey the genomic arrangement of tRNA genes and pseudogenes in a wide range of eukaryotes to identify common patterns and taxon-specific peculiarities.

RESULTS

In line with previous reports, we find that tRNA complements evolve rapidly and tRNA gene and pseudogene locations are subject to rapid turnover. At phylum level, the distributions of the number of tRNA genes and pseudogenes numbers are very broad, with standard deviations on the order of the mean. Even among closely related species we observe dramatic changes in local organization. For instance, 65% and 87% of the tRNA genes and pseudogenes are located in genomic clusters in zebrafish and stickleback, resp., while such arrangements are relatively rare in the other three sequenced teleost fish genomes. Among basal metazoa, Trichoplax adherens has hardly any duplicated tRNA gene, while the sea anemone Nematostella vectensis boasts more than 17000 tRNA genes and pseudogenes. Dramatic variations are observed even within the eutherian mammals. Higher primates, for instance, have 616 +/- 120 tRNA genes and pseudogenes of which 17% to 36% are arranged in clusters, while the genome of the bushbaby Otolemur garnetti has 45225 tRNA genes and pseudogenes of which only 5.6% appear in clusters. In contrast, the distribution is surprisingly uniform across plant genomes. Consistent with this variability, syntenic conservation of tRNA genes and pseudogenes is also poor in general, with turn-over rates comparable to those of unconstrained sequence elements. Despite this large variation in abundance in Eukarya we observe a significant correlation between the number of tRNA genes, tRNA pseudogenes, and genome size.

CONCLUSIONS

The genomic organization of tRNA genes and pseudogenes shows complex lineage-specific patterns characterized by an extensive variability that is in striking contrast to the extreme levels of sequence-conservation of the tRNAs themselves. The comprehensive analysis of the genomic organization of tRNA genes and pseudogenes in Eukarya provides a basis for further studies into the interplay of tRNA gene arrangements and genome organization in general.

Transfer RNA processing in archaea: unusual pathways and enzymes

Transfer RNA (tRNA) molecules are highly conserved in length, sequence and structure in order to be functional in the ribosome. However, mostly in archaea, the short genes encoding tRNAs can be found disrupted, fragmented, with permutations or with non-functional mutations of conserved nucleotides. Here, we give an overview of recently discovered tRNA maturation pathways that require intricate processing steps to finally generate the standard tRNA from these unusual tRNA genes.

Micro-RNAs: viral genome and robustness of gene expression in the host

For comparing RNA rings or hairpins with reference or random ring sequences, circular versions of distances and distributions like those of Hamming and Gumbel are needed. We define these circular versions and we apply these new tools to the comparison of RNA relics (such as micro-RNAs and tRNAs) with viral genomes that have coevolved with them. Then we show how robust are the regulation networks incorporating in their boundary micro-RNAs as sources or new feedback loops involving ubiquitous proteins like p53 (which is a micro-RNA transcription factor) or oligopeptides regulating protein translation. Eventually, we propose a new coevolution game between viral and host genomes.

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.

Tri-split tRNA is a transfer RNA made from 3 transcripts that provides insight into the evolution of fragmented tRNAs in archaea

Transfer RNA (tRNA) is essential for decoding the genome sequence into proteins. In Archaea, previous studies have revealed unique multiple intron-containing tRNAs and tRNAs that are encoded on 2 separate genes, so-called split tRNAs. Here, we discovered 10 fragmented tRNA genes in the complete genome of the hyperthermoacidophilic Archaeon Caldivirga maquilingensis that are individually transcribed and further trans-spliced to generate all of the missing tRNAs encoding glycine, alanine, and glutamate. Notably, the 3 mature tRNA(Gly)'s with synonymous codons are created from 1 constitutive 3' half transcript and 4 alternatively switching transcripts, representing tRNA made from a total of 3 transcripts named a "tri-split tRNA." Expression and nucleotide sequences of 10 split tRNA genes and their joined tRNA products were experimentally verified. The intervening sequences of split tRNA have high identity to tRNA intron sequences located at the same positions in intron-containing tRNAs in related Thermoproteales species. This suggests that an evolutionary relationship between intron-containing and split tRNAs exists. Our findings demonstrate the first example of split tRNA genes in a free-living organism and a unique tri-split tRNA gene that provides further insight into the evolution of fragmented tRNAs.

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.

Comprehensive analysis of archaeal tRNA genes reveals rapid increase of tRNA introns in the order thermoproteales

The analysis of archaeal tRNA genes is becoming more important to evaluate the origin and evolution of tRNA molecule. Even with the recent accumulation of complete genomes of numerous archaeal species, several tRNA genes are still required for a full complement of the codon table. We conducted comprehensive screening of tRNA genes from 47 archaeal genomes by using a combination of different types of tRNA prediction programs and extracted a total of 2,143 reliable tRNA gene candidates including 437 intron-containing tRNA genes, which covered more than 99.9% of the codon tables in Archaea. Previously, the content of intron-containing tRNA genes in Archaea was estimated to be approximately 15% of the whole tRNA genes, and most of the introns were known to be located at canonical positions (nucleotide position between 37 and 38) of precursor tRNA (pre-tRNA). Surprisingly, we observed marked enrichment of tRNA introns in five species of the archaeal order Thermoproteales; about 70% of tRNA gene candidates were found to be intron-containing tRNA genes, half of which contained multiple introns, and the introns were located at various noncanonical positions. Sequence similarity analysis revealed that approximately half of the tRNA introns found at Thermoproteales-specific intron locations were highly conserved among several tRNA genes. Intriguingly, identical tRNA intron sequences were found within different types of tRNA genes that completely lacked exon sequence similarity, suggesting that the tRNA introns in Thermoproteales could have been gained via intron insertion events at a later stage of tRNA evolution. Moreover, although the CCA sequence at the 3' terminal of pre-tRNA is added by a CCA-adding enzyme after gene transcription in Archaea, most of the tRNA genes containing highly conserved introns already encode the CCA sequence at their 3' terminal. Based on these results, we propose possible models explaining the rapid increase of tRNA introns as a result of intron insertion events via retrotransposition of pre-tRNAs. The sequences and secondary structures of the tRNA genes and their bulge-helix-bulge motifs were registered in SPLITSdb (http://splits.iab.keio.ac.jp/splitsdb/), a novel and comprehensive database for archaeal tRNA genes.

Permuted tRNA genes of Cyanidioschyzon merolae, the origin of the tRNA molecule and the root of the Eukarya domain

An evolutionary analysis is conducted on the permuted tRNA genes of Cyanidioschyzon merolae, in which the 5' half of the tRNA molecule is codified at the 3' end of the gene and its 3' half is codified at the 5' end. This analysis has shown that permuted genes cannot be considered as derived traits but seem to possess characteristics that suggest they are ancestral traits, i.e. they originated when tRNA molecule genes originated for the first time. In particular, if the hypothesis that permuted genes are a derived trait were true, then we should not have been able to observe that the most frequent class of permuted genes is that of the anticodon loop type, for the simple reason that this class would derive by random permutation from a class of non-permuted tRNA genes, which instead is the rarest. This would not explain the high frequency with which permuted tRNA genes with perfectly separate 5' and 3' halves were observed. Clearly the mechanism that produced this class of permuted genes would envisage the existence, in an advanced stage of evolution, of minigenes codifying for the 5' and 3' halves of tRNAs which were assembled in a permuted way at the origin of the tRNA molecule, thus producing a high frequency of permuted genes of the class here referred. Therefore, this evidence supports the hypothesis that the genes of the tRNA molecule were assembled by minigenes codifying for hairpin-like RNA molecules, as suggested by one model for the origin of tRNA [Di Giulio, M., 1992. On the origin of the transfer RNA molecule. J. Theor. Biol. 159, 199-214; Di Giulio, M., 1999. The non-monophyletic origin of tRNA molecule. J. Theor. Biol. 197, 403-414]. Moreover, the late assembly of the permuted genes of C. merolae, as well as their ancestrality, strengthens the hypothesis of the polyphyletic origins of these genes. Finally, on the basis of the uniqueness and the ancestrality of these permuted genes, I suggest that the root of the Eukarya domain is in the super-ensemble of the Plantae and that the Rhodophyta to which C. merolae belongs are the first line of divergence.