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

The discovery of separate 5' and 3' halves of transfer RNA (tRNA) molecules-so-called split tRNA-in the archaeal parasite Nanoarchaeum equitans made us wonder whether ancestral tRNA was encoded on 1 or 2 genes. We performed a comprehensive phylogenetic analysis of tRNAs in 45 archaeal species to explore the relationship between the three types of tRNAs (nonintronic, intronic and split). We classified 1953 mature tRNA sequences into 22 clusters. All split tRNAs have shown phylogenetic relationships with other tRNAs possessing the same anticodon. We also mimicked split tRNA by artificially separating the tRNA sequences of 7 primitive archaeal species at the anticodon and analyzed the sequence similarity and diversity of the 5' and 3' tRNA halves. Network analysis revealed specific characteristics of and topological differences between the 5' and 3' tRNA halves: the 5' half sequences were categorized into 6 distinct groups with a sequence similarity of >80%, while the 3' half sequences were categorized into 9 groups with a higher sequence similarity of >88%, suggesting different evolutionary backgrounds of the 2 halves. Furthermore, the combinations of 5' and 3' halves corresponded with the variation of amino acids in the codon table. We found not only universally conserved combinations of 5'-3' tRNA halves in tRNA(iMet), tRNA(Thr), tRNA(Ile), tRNA(Gly), tRNA(Gln), tRNA(Glu), tRNA(Asp), tRNA(Lys), tRNA(Arg) and tRNA(Leu) but also phylum-specific combinations in tRNA(Pro), tRNA(Ala), and tRNA(Trp). Our results support the idea that tRNA emerged through the combination of separate genes and explain the sequence diversity that arose during archaeal tRNA evolution.     

The non-monophyletic origin of the tRNA molecule and the origin of genes only after the evolutionary stage of the last universal common ancestor (LUCA)

A model has been proposed suggesting that the tRNA molecule must have originated by direct duplication of an RNA hairpin structure [Di Giulio, M., 1992. On the origin of the transfer RNA molecule. J. Theor. Biol. 159, 199-214]. A non-monophyletic origin of this molecule has also been theorized [Di Giulio, M., 1999. The non-monophyletic origin of tRNA molecule. J. Theor. Biol. 197, 403-414]. In other words, the tRNA genes evolved only after the evolutionary stage of the last universal common ancestor (LUCA) through the assembly of two minigenes codifying for different RNA hairpin structures, which is what the exon theory of genes suggests when it is applied to the model of tRNA origin. Recent observations strongly corroborate this theorization because it has been found that some tRNA genes are completely separate in two minigenes codifying for the 5' and 3' halves of this molecule [Randau, L., et al., 2005a. Nanoarchaeum equitans creates functional tRNAs from separate genes for their 5'- and 3'-halves. Nature 433, 537-541]. In this paper it is shown that these tRNA genes codifying for the 5' and 3' halves of this molecule are the ancestral form from which the tRNA genes continuously codifying for the complete tRNA molecule are thought to have evolved. This, together with the very existence of completely separate tRNA genes codifying for their 5' and 3' halves, proves a non-monophyletic origin for tRNA genes, as a monophyletic origin would exclude the existence of these genes which have, on the contrary, been observed. Here the polyphyletic origin of genes codifying for proteins is also suggested and discussed. Moreover, a hypothesis is advanced to suggest that the LUCA might have had a fragmented genome made up of RNA and the possibility that 'Paleokaryotes' may exist is outlined. Finally, the characteristic of the indivisibility of homology that these polyphyletic origins seem to remove at the sequence level is discussed.     

SPLITS: a new program for predicting split and intron-containing tRNA genes at the genome level

In the archaea, some tRNA precursors contain intron(s) not only in the anticodon loop region but also in diverse sites of the gene (intron-containing tRNA or cis-spliced tRNA). The parasite Nanoarchaeum equitans, a member of the Nanoarchaeota kingdom, creates functional tRNA from separate genes, one encoding the 5'-half and the other the 3'-half (split tRNA or trans-spliced tRNA). Although recent genome projects have revealed a huge amount of nucleotide sequence data in the archaea, a comprehensive methodology for intron-containing and split tRNA searching is yet to be established. We therefore developed SPLITS, which is aimed at searching for any type of tRNA gene and is especially focused on intron-containing tRNAs or split tRNAs at the genome level. SPLITS initially predicts the bulge-helix-bulge splicing motif (a well-known, required structure in archaeal pre-tRNA introns) to determine and remove the intronic regions of tRNA genes. The intron-removed DNA sequences are automatically queried to tRNAscan-SE. SPLITS can predict known tRNAs with single introns located at unconventional sites on the genes (100%), tRNAs with double introns (85.7%), and known split tRNAs (100%). Our program will be very useful for identifying novel tRNA genes after completion of genome projects. The SPLITS source code is freely downloadable at http://splits.iab.keio.ac.jp/.     

The tree of life might be rooted in the branch leading to Nanoarchaeota

It is suggested that the tree of life might be rooted in the domain of the Archaea, in the branch leading to the phylum of Nanoarchaeota. This hypothesis seems to be corroborated by the uniqueness and ancestrality of some traits possessed by Nanoarchaeum equitans, such as split genes separately codifying for the 5' and 3' halves of the tRNA molecule. These half genes are the oldest ancestral form of gene we have ever seen. This, along with the absence of operons from the genome of N. equitans, would seem to indicate that this genome is a molecular fossil of the evolutionary stage which the ancestral genomes had reached when the first lines of divergence were established. Moreover, the late appearance of DNA coinciding with the rooting of the universal phylogenetic tree would make the genome of N. equitans a witness to this fundamental event.     

Evolutionary and functional genomics of the Archaea

In the past two years, archaeal genomics has achieved several breakthroughs. On the evolutionary front the most exciting development was the sequencing and analysis of the genome of Nanoarchaeum equitans, a tiny parasitic organism that has only approximately 540 genes. The genome of Nanoarchaeum shows signs of extreme rearrangement including the virtual absence of conserved operons and the presence of several split genes. Nanoarchaeum is distantly related to other archaea, and it has been proposed to represent a deep archaeal branch that is distinct from Euryarchaeota and Crenarchaeota. This would imply that many features of its gene repertoire and genome organization might be ancestral. However, additional genome analysis has provided a more conservative suggestion - that Nanoarchaeum is a highly derived euryarchaeon. Also there have been substantial developments in functional genomics, including the discovery of the elusive aminoacyl-tRNA synthetase that is involved in both the biosynthesis of cysteine and its incorporation into proteins in methanogens, and the first experimental validation of the predicted archaeal exosome.     

Assembly of the mitochondrial membrane system. Sequences of yeast mitochondrial tRNA genes

Two cytoplasmic "petite" (rho-) clones of Saccharomyces cerevisiae have been selected for the retention of the aspartic acid tRNA gene. The two clones, designated DS200/A102 and DS200/A5, have tandemly repeated segments of mitochondrial DNA (mtDNA) with unit lengths of 1,000 and 6,400 base pairs, respectively. The DS200/A102 genome has a single tRNA gene with a 3'-CUG-5' anticodon capable of recognizing the 5'-GAC-3' and 5'-GAU-3' codons for aspartic acid. The mtDNA segment of DS200/A102 has been determined to represent the wild type sequence from 5.3 to 6.8 map units. The genome of DS200/A5 is more complex encompassing the region of wild type mtDNA from 3.5 to 12.7 units. A continuous sequence has been obtained from 3.5 to 8.6 units. In addition to the aspartic acid tRNA, this region codes for the tRNAUGCAla,tRNAUCUArg, tRNAACGArg, tRNAGCUSer,tRNAUCCGly and tRNAUUULys. The DNA sequence of the DS200/A5 genome has allowed us to deduce the secondary structures of the seven tRNAs and to assign precise map positions for their genes. All the tRNAs except tRNA GUCAsp exhibit most of the invariant features of prokaryotic and eukaryotic tRNAs. The aspartic acid tRNA has unusual D and T psi C loops. The structure of this tRNA is similar to the mitochondrial initiator tRNA of Neurospora crassa (Heckman, J.E., Hecker, L.I., Shwartzbach, S.D., Barnett, W.E., Baumstark, B., and RajBhandary, U.L. Cell 13, 83-95).     

RNA processing in the minimal organism Nanoarchaeum equitans

ABSTRACT: BACKGROUND: The minimal genome of the tiny, hyperthermophilic archaeon Nanoarchaeum equitans contains several fragmented genes and revealed unusual RNA processing pathways. These include the maturation of tRNA molecules via the trans-splicing of tRNA halves and genomic rearrangements to compensate for the absence of RNase P. RESULTS: Here, the RNA processing events in the N. equitans cell are analyzed using RNA-Seq deep sequencing methodology. All tRNA half precursor and tRNA termini were determined and support the tRNA trans-splicing model. The processing of CRISPR RNAs from two CRISPR clusters was verified. Twenty-seven C/D box small RNAs (sRNAs) and a H/ACA box sRNA were identified. The C/D box sRNAs were found to flank split genes, to form dicistronic tRNA-sRNA precursors and to be encoded within the tRNAMet intron. CONCLUSIONS: The presented data provide an overview of the production and usage of small RNAs in a cell that has to survive with a highly reduced genome. N. equitans lost many essential metabolic pathways but maintains highly active CRISPR/Cas and rRNA modification systems that appear to play an important role in genome fragmentation.     

A Comparison Among the Models Proposed to Explain the Origin of the tRNA Molecule: A Synthesis

A comparison is made among all the models proposed to explain the origin of the tRNA molecule. The conclusion reached is that, for the model predicting that the tRNA molecule originated after the assembly of two hairpin-like structures, molecular fossils have been found in the half-genes of the tRNAs of Nanoarchaeum equitans. These might be the witnesses of the transition stage predicted by the model through which the evolution of the tRNA molecule passed, thus providing considerable corroboration for this model.     

tRNA 3' end maturation in archaea has eukaryotic features: the RNase Z from Haloferax volcanii

Here, we report the first characterization and partial purification of an archaeal tRNA 3' processing activity, the RNase Z from Haloferax volcanii. The activity identified here is an endonuclease, which cleaves tRNA precursors 3' to the discriminator. Thus tRNA 3' processing in archaea resembles the eukaryotic 3' processing pathway. The archaeal RNase Z has a KCl optimum at 5mM, which is in contrast to the intracellular KCl concentration being as high as 4M KCl.The archaeal RNase Z does process 5' extended and intron-containing pretRNAs but with a much lower efficiency than 5' matured, intronless pretRNAs. At least in vitro there is thus no defined order for 5' and 3' processing and splicing. A heterologous precursor tRNA is cleaved efficiently by the archaeal RNase Z. Experiments with precursors containing mutated tRNAs revealed that removal of the anticodon arm reduces cleavage efficiency only slightly, while removal of D and T arm reduces processing effciency drastically, even down to complete inhibition. Comparison with its nuclear and mitochondrial homologs revealed that the substrate specificity of the archaeal RNase Z is narrower than that of the nuclear RNase Z but broader than that of the mitochondrial RNase Z.     

3'-processing of yeast tRNATrp precedes 5'-processing

Previous analyses of eukaryotic pre-tRNAs processing have reported that 5'-cleavage by RNase P precedes 3'-maturation. Here we report that in contrast to all other yeast tRNAs analyzed to date, tRNA(Trp) undergoes 3'-maturation prior to 5'-cleavage. Despite its unusual processing pathway, pre-tRNA(Trp) resembles other pre-tRNAs, showing dependence on the essential Lsm proteins for normal processing and efficient association with the yeast La homolog, Lhp1p. tRNA(Trp) is also unusual in not requiring Lhp1p for 3' processing and stability. However, other Lhp1p-independent tRNAs, tRNA(2)(Lys) and tRNA(1)(Ile), follow the normal pathway of 5'-processing prior to 3-processing.     

End processing precedes mitochondrial importation and editing of tRNAs in Leishmania tarentolae

All mitochondrial tRNAs in Leishmania tarentolae are encoded in the nuclear genome and imported into the mitochondrion from the cytosol. One imported tRNA (tRNA(Trp)) is edited by a C to U modification at the first position of the anticodon. To determine the in vivo substrates for mitochondrial tRNA importation as well as tRNA editing, we examined the subcellular localization and extent of 5'- and 3'-end maturation of tRNA(Trp)(CCA), tRNA(Ile)(UAU), tRNA(Gln)(CUG), tRNA(Lys)(UUU), and tRNA(Val)(CAC). Nuclear, cytosolic, and mitochondrial fractions were obtained with little cross-contamination, as determined by Northern analysis of specific marker RNAs. tRNA(Gln) was mainly cytosolic in localization; tRNA(Ile) and tRNA(Lys) were mainly mitochondrial; and tRNA(Trp) and tRNA(Val) were shared between the two compartments. 5'- and 3'-extended precursors of all five tRNAs were present only in the nuclear fraction, suggesting that the mature tRNAs represent the in vivo substrates for importation into the mitochondrion. Consistent with this model, T7-transcribed mature tRNA(Ile) underwent importation in vitro into isolated mitochondria more efficiently than 5'-extended precursor tRNA(Ile). 5'-Extended precursor tRNA(Trp) was found to be unedited, which is consistent with a mitochondrial localization of this editing reaction. T7-transcribed unedited tRNA(Trp) was imported in vitro more efficiently than edited tRNA(Trp), suggesting the presence of importation determinants in the anticodon.     

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.     

Transfer RNAs of potato (Solanum tuberosum) mitochondria have different genetic origins.

Total transfer RNAs were extracted from highly purified potato mitochondria. From quantitative measurements, the in vivo tRNA concentration in mitochondria was estimated to be in the range of 60 microM. Total potato mitochondrial tRNAs were fractionated by two-dimensional polyacrylamide gel electrophoresis. Thirty one individual tRNAs, which could read all sense codons, were identified by aminoacylation, sequencing or hybridization to specific oligonucleotides. The tRNA population that we have characterized comprises 15 typically mitochondrial, 5 'chloroplast-like' and 11 nuclear-encoded species. One tRNA(Ala), 2 tRNAs(Arg), 1 tRNA(Ile), 5 tRNAs(Leu) and 2 tRNAs(Thr) were shown to be coded for by nuclear DNA. A second, mitochondrial-encoded, tRNA(Ile) was also found. Five 'chloroplast-like' tRNAs, tRNA(Trp), tRNA(Asn), tRNA(His), tRNA(Ser)(GGA) and tRNA(Met)m, presumably transcribed from promiscuous chloroplast DNA sequences inserted in the mitochondrial genome, were identified, but, in contrast to wheat (1), potato mitochondria do not seem to contain 'chloroplast-like' tRNA(Cys) and tRNA(Phe). The two identified tRNAs(Val), as well as the tRNA(Gly), were found to be coded for by the mitochondrial genome, which again contrasts with the situation in wheat, where the mitochondrial genome apparently contains no tRNA(Val) or tRNA(Gly) gene (2).     

Gene transfers from nanoarchaeota to an ancestor of diplomonads and parabasalids

Rare evolutionary events, such as lateral gene transfers and gene fusions, may be useful to pinpoint, and correlate the timing of, key branches across the tree of life. For example, the shared possession of a transferred gene indicates a phylogenetic relationship among organismal lineages by virtue of their shared common ancestral recipient. Here, we present phylogenetic analyses of prolyl-tRNA and alanyl-tRNA synthetase genes that indicate lateral gene transfer events to an ancestor of the diplomonads and parabasalids from lineages more closely related to the newly discovered archaeal hyperthermophile Nanoarchaeum equitans (Nanoarchaeota) than to Crenarchaeota or Euryarchaeota. The support for this scenario is strong from all applied phylogenetic methods for the alanyl-tRNA sequences, whereas the phylogenetic analyses of the prolyl-tRNA sequences show some disagreements between methods, indicating that the donor lineage cannot be identified with a high degree of certainty. However, in both trees, the diplomonads and parabasalids branch together within the Archaea, strongly suggesting that these two groups of unicellular eukaryotes, often regarded as the two earliest independent offshoots of the eukaryotic lineage, share a common ancestor to the exclusion of the eukaryotic root. Unfortunately, the phylogenetic analyses of these two aminoacyl-tRNA synthetase genes are inconclusive regarding the position of the diplomonad/parabasalid group within the eukaryotes. Our results also show that the lineage leading to Nanoarchaeota branched off from Euryarchaeota and Crenarchaeota before the divergence of diplomonads and parabasalids, that this unexplored archaeal diversity, currently only represented by the hyperthermophilic organism Nanoarchaeum equitans, may include members living in close proximity to mesophilic eukaryotes, and that the presence of split genes in the Nanoarchaeum genome is a derived feature.     

A role for tRNA(His) guanylyltransferase (Thg1)-like proteins from Dictyostelium discoideum in mitochondrial 5'-tRNA editing

Genes with sequence similarity to the yeast tRNA(His) guanylyltransferase (Thg1) gene have been identified in all three domains of life, and Thg1 family enzymes are implicated in diverse processes, ranging from tRNA(His) maturation to 5'-end repair of tRNAs. All of these activities take advantage of the ability of Thg1 family enzymes to catalyze 3'-5' nucleotide addition reactions. Although many Thg1-containing organisms have a single Thg1-related gene, certain eukaryotic microbes possess multiple genes with sequence similarity to Thg1. Here we investigate the activities of four Thg1-like proteins (TLPs) encoded by the genome of the slime mold, Dictyostelium discoideum (a member of the eukaryotic supergroup Amoebozoa). We show that one of the four TLPs is a bona fide Thg1 ortholog, a cytoplasmic G(-1) addition enzyme likely to be responsible for tRNA(His) maturation in D. discoideum. Two other D. discoideum TLPs exhibit biochemical activities consistent with a role for these enzymes in mitochondrial 5'-tRNA editing, based on their ability to efficiently repair the 5' ends of mitochondrial tRNA editing substrates. Although 5'-tRNA editing was discovered nearly two decades ago, the identity of the protein(s) that catalyze this activity has remained elusive. This article provides the first identification of any purified protein that appears to play a role in the 5'-tRNA editing reaction. Moreover, the presence of multiple Thg1 family members in D. discoideum suggests that gene duplication and divergence during evolution has resulted in paralogous proteins that use 3'-5' nucleotide addition reactions for diverse biological functions in the same organism.     

The in vivo stability, maturation and aminoacylation of anticodon-substituted Escherichia coli initiator methionine tRNAs

We have constructed eight anticodon-modified Escherichia coli initiator methionine (fMet) tRNAs by insertion of synthetic ribotrinucleotides between two fragments ('half molecules') derived from the initiator tRNA. The trinucleotides, namely CAU (the normal anticodon), CAA, CAC, CAG, GAA, GAC, GAG and GAU, were joined to the 5' and 3' tRNA fragments with T4 RNA ligase. The strategy of reconstruction permitted the insertion of radioactive 32P label between nucleotides 36 and 37. tRNAs were microinjected into the cytoplasm of Xenopus laevis oocytes, and the following properties were evaluated: the stability of these eubacterial tRNA variants in the eukaryotic oocytes; the enzymatic modification of the adenosine at position 37 (3' adjacent to the anticodon) and aminoacylation of the chimeric tRNAs by endogenous oocyte aminoacyl-tRNA synthetases. In contrast to other variants, the two RNAs having CAU and GAU anticodons were stable and underwent quantitative modification at A-37. These results show that the enzyme responsible for the modification of A-37 to N-[N-(9-beta-D-ribofuranosylpurine-6-yl)carbamoyl]threonine (t6A) is present in the cytoplasm of oocytes and is very sensitive to the anticodon environment of the tRNA. Also, these same GAU and CAU anticodon-containing tRNAs are fully aminoacylated with the heterologous oocyte aminoacyl-tRNA synthetases in vivo. During the course of this work we developed a generally applicable assay for the aminoacylation of femtomole amounts of labelled tRNAs.     

Endonucleolytic processing of CCA-less tRNA precursors by RNase Z in Bacillus subtilis

In contrast to Escherichia coli, where the 3' ends of tRNAs are primarily generated by exoribonucleases, maturation of the 3' end of tRNAs is catalysed by an endoribonuclease, known as RNase Z (or 3' tRNase), in many eukaryotic and archaeal systems. RNase Z cleaves tRNA precursors 3' to the discriminator base. Here we show that this activity, previously unsuspected in bacteria, is encoded by the yqjK gene of Bacillus subtilis. Decreased yqjK expression leads to an accumulation of a population of B.subtilis tRNAs in vivo, none of which have a CCA motif encoded in their genes, and YqjK cleaves tRNA precursors with the same specificity as plant RNase Z in vitro. We have thus renamed the gene rnz. A CCA motif downstream of the discriminator base inhibits RNase Z activity in vitro, with most of the inhibition due to the first C residue. Lastly, tRNAs with long 5' extensions are poor substrates for cleavage, suggesting that for some tRNAs, processing of the 5' end by RNase P may have to precede RNase Z cleavage.