Kosuke Fujishima

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 tRNAGlys 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.

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 tRNAiMet, tRNAThr, tRNAIle, tRNAGly, tRNAGln, tRNAGlu, tRNAAsp, tRNALys, tRNAArg and tRNALeu but also phylum-specific combinations in tRNAPro, tRNAAla, and tRNATrp. 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.

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.

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.

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 tRNASer and tRNALeu 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.

Proteome-wide prediction of novel DNA/RNA-binding proteins using amino acid composition and periodicity in the hyperthermophilic archaeon Pyrococcus furiosus.

Abstract Proteins play a critical role in complex biological systems, yet about half of the proteins in publicly available databases are annotated as functionally unknown. Proteome-wide functional classification using bioinformatics approaches thus is becoming an important method for revealing unknown protein functions. Using the hyperthermophilic archaeon Pyrococcus furiosus as a model species, we used the support vector machine (SVM) method to discriminate DNA/RNA-binding proteins from proteins with other functions, using amino acid composition and periodicities as feature vectors. We defined this value as the composition score (CO) and periodicity score (PD). The P. furiosus proteins were classified into three classes (I–III) on the basis of the two-dimensional correlation analysis of CO score and PD score. As a result, approximately 87% of the functionally known proteins categorized as class I proteins (CO score + PD score > 0.6) were found to be DNA/RNA-binding proteins. Applying the two-dimensional correlation analysis to the 994 hypothetical proteins in P. furiosus, a total of 151 proteins were predicted to be novel DNA/RNA-binding protein candidates. DNA/RNA-binding activities of randomly chosen hypothetical proteins were experimentally verified. Six out of seven candidate proteins in class I possessed DNA/RNA-binding activities, supporting the efficacy of our method.

A novel three-unit tRNA splicing endonuclease found in ultrasmall Archaea possesses broad substrate specificity

tRNA splicing endonucleases, essential enzymes found in Archaea and Eukaryotes, are involved in the processing of pre-tRNA molecules. In Archaea, three types of splicing endonuclease [homotetrameric: α4, homodimeric: α2, and heterotetrameric: (αβ)2] have been identified, each representing different substrate specificity during the tRNA intron cleavage. Here, we discovered a fourth type of archaeal tRNA splicing endonuclease (ε2) in the genome of the acidophilic archaeon Candidatus Micrarchaeum acidiphilum, referred to as ARMAN-2 and its closely related species, ARMAN-1. The enzyme consists of two duplicated catalytic units and one structural unit encoded on a single gene, representing a novel three-unit architecture. Homodimeric formation was confirmed by cross-linking assay, and site-directed mutagenesis determined that the conserved L10-pocket interaction between catalytic and structural unit is necessary for the assembly. A tRNA splicing assay reveal that ε2 endonuclease cleaves both canonical and non-canonical bulge–helix–bulge motifs, similar to that of (αβ)2 endonuclease. Unlike other ARMAN and Euryarchaeota, tRNAs found in ARMAN-2 are highly disrupted by introns at various positions, which again resemble the properties of archaeal species with (αβ)2 endonuclease. Thus, the discovery of ε2 endonuclease in an archaeon deeply branched within Euryarchaeota represents a new example of the coevolution of tRNA and their processing enzymes.

Characterization of RNase HII substrate recognition using RNase HII-argonaute chimaeric enzymes from Pyrococcus furiosus

RNase H (ribonuclease H) is an endonuclease that cleaves the RNA strand of RNA-DNA duplexes. It has been reported that the three-dimensional structure of RNase H is similar to that of the PIWI domain of the Pyrococcus furiosus Ago (argonaute) protein, although the two enzymes share almost no similarity in their amino acid sequences. Eukaryotic Ago proteins are key components of the RNA-induced silencing complex and are involved in microRNA or siRNA (small interfering RNA) recognition. In contrast, prokaryotic Ago proteins show greater affinity for RNA-DNA hybrids than for RNA-RNA hybrids. Interestingly, we found that wild-type Pf-RNase HII (P. furiosus, RNase HII) digests RNA-RNA duplexes in the presence of Mn2+ ions. To characterize the substrate specificity of Pf-RNase HII, we aligned the amino acid sequences of Pf-RNase HII and Pf-Ago, based on their protein secondary structures. We found that one of the conserved secondary structural regions (the fourth beta-sheet and the fifth alpha-helix of Pf-RNase HII) contains family-specific amino acid residues. Using a series of Pf-RNase HII-Pf-Ago chimaeric mutants of the region, we discovered that residues Asp110, Arg113 and Phe114 are responsible for the dsRNA (double-stranded RNA) digestion activity of Pf-RNase HII. On the basis of the reported three-dimensional structure of Ph-RNase HII from Pyrococcus horikoshii, we built a three-dimensional structural model of RNase HII complexed with its substrate, which suggests that these amino acids are located in the region that discriminates DNA from RNA in the non-substrate strand of the duplexes.

X-ray structure of the fourth type of archaeal tRNA splicing endonuclease: insights into the evolution of a novel three-unit composition and a unique loop involved in broad substrate specificity.

Cleavage of introns from precursor transfer RNAs (tRNAs) by tRNA splicing endonuclease (EndA) is essential for tRNA maturation in Archaea and Eukarya. In the past, archaeal EndAs were classified into three types (α′2, α4 and α2β2) according to subunit composition. Recently, we have identified a fourth type of archaeal EndA from an uncultivated archaeon Candidatus Micrarchaeum acidiphilum, referred to as ARMAN-2, which is deeply branched within Euryarchaea. The ARMAN-2 EndA forms an ε2 homodimer and has broad substrate specificity like the α2β2 type EndAs found in Crenarchaea and Nanoarchaea. However, the precise architecture of ARMAN-2 EndA was unknown. Here, we report the crystal structure of the ε2 homodimer of ARMAN-2 EndA. The structure reveals that the ε protomer is separated into three novel units (αN, α and βC) fused by two distinct linkers, although the overall structure of ARMAN-2 EndA is similar to those of the other three types of archaeal EndAs. Structural comparison and mutational analyses reveal that an ARMAN-2 type-specific loop (ASL) is involved in the broad substrate specificity and that K161 in the ASL functions as the RNA recognition site. These findings suggest that the broad substrate specificities of ε2 and α2β2 EndAs were separately acquired through different evolutionary processes.

Nematode-specific tRNAs that decode an alternative genetic code for leucine

Class II transfer RNAs (tRNAs), including tRNALeu and tRNASer, have an additional stem and loop structure, the long variable arm (V-arm). Here, we describe Class II tRNAs with a unique anticodon corresponding to neither leucine nor serine. Because these tRNAs are specifically conserved among the nematodes, we have called them ‘nematode-specific V-arm-containing tRNAs’ (nev-tRNAs). The expression of nev-tRNA genes in Caenorhabditis elegans was confirmed experimentally. A comparative sequence analysis suggested that the nev-tRNAs derived phylogenetically from tRNALeu. In vitro aminoacylation assays showed that nev-tRNAGly and nev-tRNAIle are only charged with leucine, which is inconsistent with their anticodons. Furthermore, the deletion and mutation of crucial determinants for leucylation in nev-tRNA led to a marked loss of activity. An in vitro translation analysis showed that nev-tRNAGly decodes GGG as leucine instead of the universal glycine code, indicating that nev-tRNAs can be incorporated into ribosomes and participate in protein biosynthesis. Our findings provide the first example of unexpected tRNAs that do not consistently obey the general translation rules for higher eukaryotes.