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

The 3 31 Nucleotide Minihelix tRNA Evolution Theorem and the Origin of Life

There are no theorems (proven theories) in the biological sciences. We propose that the 3 31 nt minihelix tRNA evolution theorem be universally accepted as one. The 3 31 nt minihelix theorem completely describes the evolution of type I and type II tRNAs from ordered precursors (RNA repeats and inverted repeats). Despite the diversification of tRNAome sequences, statistical tests overwhelmingly support the theorem. Furthermore, the theorem relates the dominant pathway for the origin of life on Earth, specifically, how tRNAomes and the genetic code may have coevolved. Alternate models for tRNA evolution (i.e., 2 minihelix, convergent and accretion models) are falsified. In the context of the pre-life world, tRNA was a molecule that, via mutation, could modify anticodon sequences and teach itself to code. Based on the tRNA sequence, we relate the clearest history to date of the chemical evolution of life. From analysis of tRNA evolution, ribozyme-mediated RNA ligation was a primary driving force in the evolution of complexity during the pre-life-to-life transition. TRNA formed the core for the evolution of living systems on Earth.

Comparisons between small ribosomal RNA and theoretical minimal RNA ring secondary structures confirm phylogenetic and structural accretion histories

Ribosomal RNAs are complex structures that presumably evolved by tRNA accretions. Statistical properties of tRNA secondary structures correlate with genetic code integration orders of their cognate amino acids. Ribosomal RNA secondary structures resemble those of tRNAs with recent cognates. Hence, rRNAs presumably evolved from ancestral tRNAs. Here, analyses compare secondary structure subcomponents of small ribosomal RNA subunits with secondary structures of theoretical minimal RNA rings, presumed proto-tRNAs. Two independent methods determined different accretion orders of rRNA structural subelements: (a) classical comparative homology and phylogenetic reconstruction, and (b) a structural hypothesis assuming an inverted onion ring growth where the three-dimensional ribosome's core is most ancient and peripheral elements most recent. Comparisons between (a) and (b) accretions orders with RNA ring secondary structure scales show that recent rRNA subelements are: 1. more like RNA rings with recent cognates, indicating ongoing coevolution between tRNA and rRNA secondary structures; 2. less similar to theoretical minimal RNA rings with ancient cognates. Our method fits (a) and (b) in all examined organisms, more with (a) than (b). Results stress the need to integrate independent methods. Theoretical minimal RNA rings are potential evolutionary references for any sequence-based evolutionary analyses, independent of the focal data from that study.

Accretion history of large ribosomal subunits deduced from theoretical minimal RNA rings is congruent with histories derived from phylogenetic and structural methods

Accretions of tRNAs presumably formed the large complex ribosomal RNA structures. Similarities of tRNA secondary structures with rRNA secondary structures increase with the integration order of their cognate amino acid in the genetic code, indicating tRNA evolution towards rRNA-like structures. Here analyses rank secondary structure subelements of three large ribosomal RNAs (Prokaryota: Archaea: Thermus thermophilus; Bacteria: Escherichia coli; Eukaryota: Saccharomyces cerevisiae) in relation to their similarities with secondary structures formed by presumed proto-tRNAs, represented by 25 theoretical minimal RNA rings. These ranks are compared to those derived from two independent methods (ranks provide a relative evolutionary age to the rRNA substructure), (a) cladistic phylogenetic analyses and (b) 3D-crystallography where core subelements are presumed ancient and peripheral ones recent. Comparisons of rRNA secondary structure subelements with RNA ring secondary structures show congruence between ranks deduced by this method and both (a) and (b) (more with (a) than (b)), especially for RNA rings with predicted ancient cognate amino acid. Reconstruction of accretion histories of large rRNAs will gain from adequately integrating information from independent methods. Theoretical minimal RNA rings, sequences deterministically designed in silico according to specific coding constraints, might produce adequate scales for prebiotic and early life molecular evolution.

An RNA Ring was Not the Progenitor of the tRNA Molecule

I analyzed the model that suggests that an RNA ring might have been the progenitor of the tRNA molecule (Demongeot and Moreira in J Theor Biol 249:314-324, 2007; Demongeot and Seligmann in J Mol Evol 1-23, 2019a; Demongeot and Norris in Life 9(2):51, 2019). In particular, I analyze three ways in which this precursor, especially in its RNA hairpin form, could have evolved into the complete tRNA molecule. These three modalities are based on multiple duplication events, and therefore, appear to be less parsimonious than that which assumes that this molecule originated through one duplication of a single hairpin structure. The conclusion is, therefore, that the latter model appears to be preferable with respect to that of the RNA ring, also because there are many independent observations and some of a historical nature that would corroborate it in an extraordinary way.

How amino acids and peptides shaped the RNA world

The “RNA world” hypothesis is seen as one of the main contenders for a viable theory on the origin of life. Relatively small RNAs have catalytic power, RNA is everywhere in present-day life, the ribosome is seen as a ribozyme, and rRNA and tRNA are crucial for modern protein synthesis. However, this view is incomplete at best. The modern protein-RNA ribosome most probably is not a distorted form of a “pure RNA ribosome” evolution started out with. Though the oldest center of the ribosome seems “RNA only”, we cannot conclude from this that it ever functioned in an environment without amino acids and/or peptides. Very small RNAs (versatile and stable due to basepairing) and amino acids, as well as dipeptides, coevolved. Remember, it is the amino group of aminoacylated tRNA that attacks peptidyl-tRNA, destroying the bond between peptide and tRNA. This activity of the amino acid part of aminoacyl-tRNA illustrates the centrality of amino acids in life. With the rise of the “RNA world” view of early life, the pendulum seems to have swung too much towards the ribozymatic part of early biochemistry. The necessary presence and activity of amino acids and peptides is in need of highlighting. In this article, we try to bring the role of the peptide component of early life back into focus. We argue that an RNA world completely independent of amino acids never existed.

The stereochemical basis of the genetic code and the (mostly) autotrophic origin of life

Spark-tube experiments and analysis of meteorite contents have led to the widespread notion that abiotic organic molecules were the first life components. However, there is a contradiction between the abundance of simple molecules, such as the amino acids glycine and alanine, observed in these studies, and the minimal functional complexity that even the least sophisticated living system should require. I will argue that although simple abiotic molecules must have primed proto-metabolic pathways, only Darwinian evolving systems could have generated life. This condition may have been initially fulfilled by both replicating RNAs and autocatalytic reaction chains, such as the reductive citric acid cycle. The interactions between nucleotides and biotic amino acids, which conferred new functionalities to the former, also resulted in the progressive stereochemical recognition of the latter by cognate anticodons. At this point only large enough amino acids would be recognized by the primordial RNA adaptors and could polymerize forming the first peptides. The gene duplication of RNA adaptors was a crucial event. By removing one of the anticodons from the acceptor stem the new RNA adaptor liberated itself from the stereochemical constraint and could be acylated by smaller amino acids. The emergence of messenger RNA and codon capture followed.

Mitochondrial swinger replication: DNA replication systematically exchanging nucleotides and short 16S ribosomal DNA swinger inserts

Assuming systematic exchanges between nucleotides (swinger RNAs) resolves genomic 'parenthood' of some orphan mitochondrial transcripts. Twenty-three different systematic nucleotide exchanges (bijective transformations) exist. Similarities between transcription and replication suggest occurrence of swinger DNA. GenBank searches for swinger DNA matching the 23 swinger versions of human and mouse mitogenomes detect only vertebrate mitochondrial swinger DNA for swinger type AT+CG (from five different studies, 149 sequences) matching three human and mouse mitochondrial genes: 12S and 16S ribosomal RNAs, and cytochrome oxidase subunit I. Exchange A<->T+C<->G conserves self-hybridization properties, putatively explaining swinger biases for rDNA, against protein coding genes. Twenty percent of the regular human mitochondrial 16S rDNA consists of short swinger repeats (from 13 exchanges). Swinger repeats could originate from recombinations between regular and swinger DNA: duplicated mitochondrial genes of the parthenogenetic gecko Heteronotia binoei include fewer short A<->T+C<->G swinger repeats than non-duplicated mitochondrial genomes of that species. Presumably, rare recombinations between female and male mitochondrial genes (and in parthenogenetic situations between duplicated genes), favors reverse-mutations of swinger repeat insertions, probably because most inserts affect negatively ribosomal function. Results show that swinger DNA exists, and indicate that swinger polymerization contributes to the genesis of genetic material and polymorphism.

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.

Identification of highly-disrupted tRNA genes in nuclear genome of the red alga, Cyanidioschyzon merolae 10D

The limited locations of tRNA introns are crucial for eukaryal tRNA-splicing endonuclease recognition. However, our analysis of the nuclear genome of an early-diverged red alga, Cyanidioschyzon merolae, demonstrated the first evidence of nuclear-encoded tRNA genes that contain ectopic and/or multiple introns. Some genes exhibited both intronic and permuted structures in which the 3′-half of the tRNA coding sequence lies upstream of the 5′-half, and an intron is inserted into either half. These highly disrupted tRNA genes, which account for 63% of all nuclear tRNA genes, are expressed via the orderly and sequential processing of bulge-helix-bulge (BHB) motifs at intron-exon junctions and termini of permuted tRNA precursors, probably by a C. merolae tRNA-splicing endonuclease with an unidentified subunit architecture. The results revealed a considerable diversity in eukaryal tRNA intron properties and endonuclease architectures, which will help to elucidate the acquisition mechanism of the BHB-mediated disrupted tRNA genes.

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.

Discovery of permuted and recently split transfer RNAs in Archaea

Background

As in eukaryotes, precursor transfer RNAs in Archaea often contain introns that are removed in tRNA maturation. Two unrelated archaeal species display unique pre-tRNA processing complexity in the form of split tRNA genes, in which two to three segments of tRNAs are transcribed from different loci, then trans-spliced to form a mature tRNA. Another rare type of pre-tRNA, found only in eukaryotic algae, is permuted, where the 3' half is encoded upstream of the 5' half, and must be processed to be functional.

Results

Using an improved version of the gene-finding program tRNAscan-SE, comparative analyses and experimental verifications, we have now identified four novel trans-spliced tRNA genes, each in a different species of the Desulfurococcales branch of the Archaea: tRNA^Asp(GUC) in Aeropyrum pernix and Thermosphaera aggregans, and tRNA^Lys(CUU) in Staphylothermus hellenicus and Staphylothermus marinus. Each of these includes features surprisingly similar to previously studied split tRNAs, yet comparative genomic context analysis and phylogenetic distribution suggest several independent, relatively recent splitting events. Additionally, we identified the first examples of permuted tRNA genes in Archaea: tRNA^iMet(CAU) and tRNA^Tyr(GUA) in Thermofilum pendens, which appear to be permuted in the same arrangement seen previously in red alga.

Conclusions

Our findings illustrate that split tRNAs are sporadically spread across a major branch of the Archaea, and that permuted tRNAs are a new shared characteristic between archaeal and eukaryotic species. The split tRNA discoveries also provide new clues to their evolutionary history, supporting hypotheses for recent acquisition via viral or other mobile elements.

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.

RESULTS

The aa-binding sites of arginine, isoleucine and tyrosine contain both their cognate triplets, anticodons and codons. We have noticed that these cases might be associated with palindrome-dinucleotides. For example, one-base shift to the left brings arginine codons CGN, with CG at 1-2 positions, to the respective anticodons NCG, with CG at 2-3 positions. Formally, the concomitant presence of codons and anticodons is also expected in the reverse situation, with codons containing palindrome-dinucleotides at their 2-3 positions, and anticodons exhibiting them at 1-2 positions. A closer analysis reveals that, surprisingly, RNA binding sites for Arg, Ile and Tyr "prefer" (exactly as in the actual genetic code) the anticodon(2-3)/codon(1-2) tetramers to their anticodon(1-2)/codon(2-3) counterparts, despite the seemingly perfect symmetry of the latter. However, since in vitro selection of aa-specific RNA aptamers apparently had nothing to do with translation, this striking preference provides a new strong support to the notion of the genetic code emerging before translation, in response to catalytic (and possibly other) needs of ancient RNA life. Consistently with the pre-translation origin of the code, we propose here a new model of tRNA origin by the gradual, Fibonacci process-like, elongation of a tRNA molecule from a primordial coding triplet and 5'DCCA3' quadruplet (D is a base-determinator) to the eventual 76 base-long cloverleaf-shaped molecule.

CONCLUSION

Taken together, our findings necessarily imply that primordial tRNAs, tRNA aminoacylating ribozymes, and (later) the translation machinery in general have been co-evolving to ''fit'' the (likely already defined) genetic code, rather than the opposite way around. Coding triplets in this primal pre-translational code were likely similar to the anticodons, with second and third nucleotides being more important than the less specific first one. Later, when the code was expanding in co-evolution with the translation apparatus, the importance of 2-3 nucleotides of coding triplets "transferred" to the 1-2 nucleotides of their complements, thus distinguishing anticodons from codons. This evolutionary primacy of anticodons in genetic coding makes the hypothesis of primal stereo-chemical affinity between amino acids and cognate triplets, the hypothesis of coding coenzyme handles for amino acids, the hypothesis of tRNA-like genomic 3' tags suggesting that tRNAs originated in replication, and the hypothesis of ancient ribozymes-mediated operational code of tRNA aminoacylation not mutually contradicting but rather co-existing in harmony.

Permuted tRNA genes in the nuclear and nucleomorph genomes of photosynthetic eukaryotes

Transfer RNA (tRNA) is a central genetic element in the decoding of genome information for all of Earth's life forms. Nevertheless, there are a great number of missing tRNAs that have been left without examination, especially in microbial genomes. Two tRNA gene families remarkable in their structure and expression mechanism have been reported: split and permuted tRNAs. Split tRNAs in archaea are encoded on the genome as two or three fragmented genes and then processed into single tRNA molecules. Permuted tRNAs are organized with the 5' and 3' halves of the gene positioned in reverse on the genome and hitherto have been found only in one tiny red alga. Here we reveal a wide-ranging distribution of permuted tRNA genes in the genomes of photosynthetic eukaryotes. This includes in the smallest eukaryotic genome known to date, the nucleomorph genome of the chlorarachniophyte alga Bigelowiella natans. Comparison between cDNA and genomic DNA sequences of two nucleomorph-encoded tRNA(Ser) genes confirms that precursors are circularized and processed into mature tRNA molecules in vivo. In the tRNA(Ser)(AGA), adenine at the wobble position of the codon is likely modified to inosine to expand capacity of the codon recognition. We also show the presence of permuted tRNAs in the ultrasmall free-living green algae Ostreococcus and Micromonas, which are closely related to the B. natans nucleomorph. Conserved intron/leader sequence structures in the intron-containing and permuted tRNAs suggest the ancient origin of the splicing machinery in the common ancestor of eukaryotes and archaea. Meanwhile, a wide but patchy distribution of permuted tRNA genes in the photosynthetic eukaryotes implies that extant permuted tRNAs might have emerged multiple times. Taken together, our data demonstrate that permuted tRNA is an evolutionarily conserved and fundamental element in tiny eukaryotic 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.