Discovery of permuted and recently split transfer RNAs in Archaea

Emerging functional principles of tRNA-derived small RNAs and other regulatory small RNAs

Recent advancements in small RNA sequencing have unveiled a previously hidden world of regulatory small noncoding RNAs (sncRNAs) that extend beyond the well-studied small interfering RNAs, microRNAs, and piwi-interacting RNAs. This exploration, starting with tRNA-derived small RNAs, has led to the discovery of a diverse universe of sncRNAs derived from various longer structured RNAs such as rRNAs, small nucleolar RNAs, small nuclear RNAs, Y RNAs, and vault RNAs, with exciting uncharted functional possibilities. In this perspective, we discuss the emerging functional principles of sncRNAs beyond the well-known RNAi-like mechanisms, focusing on those that operate independent of linear sequence complementarity but rather function in an aptamer-like fashion. Aptamers use 3D structure for specific interactions with ligands and are modulated by RNA modifications and subcellular environments. Given that aptamer-like sncRNA functions are widespread and present in species lacking RNAi, they may represent an ancient functional principle that predates RNAi. We propose a rethinking of the origin of RNAi and its relationship with these aptamer-like functions in sncRNAs and how these complementary mechanisms shape biological processes. Lastly, the aptamer-like function of sncRNAs highlights the need for caution in using small RNA mimics in research and therapeutics, as their specificity is not restricted solely to linear sequence.

Crystal structure of Nanoarchaeum equitans tyrosyl-tRNA synthetase and its aminoacylation activity toward tRNATyr with an extra guanosine residue at the 5'-terminus

tRNA^Tyr of Nanoarchaeum equitans has a remarkable feature with an extra guanosine residue at the 5'-terminus. However, the N. equitans tRNA^Tyr mutant without extra guanosine at the 5'-end was tyrosylated by tyrosyl-tRNA synthase (TyrRS). We solved the crystal structure of N. equitans TyrRS at 2.80 Å resolution. By comparing the present solved structure with the complex structures TyrRS with tRNA^Tyr of Thermus thermophilus and Methanocaldococcus jannaschii, an arginine substitution mutant of N. equitans TyrRS at Ile200 (I200R), which is the putative closest candidate to the 5'-phosphate of C1 of N. equitans tRNA^Tyr, was prepared. The I200R mutant tyrosylated not only wild-type tRNA^Tyr but also the tRNA without the G-1 residue. Further tyrosylation analysis revealed that the second base of the anticodon (U35), discriminator base (A73), and C1:G72 base pair are strong recognition sites.

tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes

tRNAscan-SE has been widely used for transfer RNA (tRNA) gene prediction for over twenty years, developed just as the first genomes were decoded. With the massive increase in quantity and phylogenetic diversity of genomes, the accurate detection and functional prediction of tRNAs has become more challenging. Utilizing a vastly larger training set, we created nearly one hundred specialized isotype- and clade-specific models, greatly improving tRNAscan-SE’s ability to identify and classify both typical and atypical tRNAs. We employ a new comparative multi-model strategy where predicted tRNAs are scored against a full set of isotype-specific covariance models, allowing functional prediction based on both the anticodon and the highest-scoring isotype model. Comparative model scoring has also enhanced the program's ability to detect tRNA-derived SINEs and other likely pseudogenes. For the first time, tRNAscan-SE also includes fast and highly accurate detection of mitochondrial tRNAs using newly developed models. Overall, tRNA detection sensitivity and specificity is improved for all isotypes, particularly those utilizing specialized models for selenocysteine and the three subtypes of tRNA genes encoding a CAU anticodon. These enhancements will provide researchers with more accurate and detailed tRNA annotation for a wider variety of tRNAs, and may direct attention to tRNAs with novel traits.

Phylogeny and evolution of chloroplast tRNAs in Adoxaceae

Chloroplasts are semiautonomous organelles found in photosynthetic plants. The major functions of chloroplasts include photosynthesis and carbon fixation, which are mainly regulated by its circular genomes. In the highly conserved chloroplast genome, the chloroplast transfer RNA genes (cp tRNA) play important roles in protein translation within chloroplasts. However, the evolution of cp tRNAs remains unclear. Thus, in the present study, we investigated the evolutionary characteristics of chloroplast tRNAs in five Adoxaceae species using 185 tRNA gene sequences. In total, 37 tRNAs encoding 28 anticodons are found in the chloroplast genome in Adoxaceae species. Some consensus sequences are found within the Ψ-stem and anticodon loop of the tRNAs. Some putative novel structures were also identified, including a new stem located in the variable region of tRNATyr in a similar manner to the anticodon stem. Furthermore, phylogenetic and evolutionary analyses indicated that synonymous tRNAs may have evolved from multiple ancestors and frequent tRNA duplications during the evolutionary process may have been primarily caused by positive selection and adaptive evolution. The transition and transversion rates are uneven among different tRNA isotypes. For all tRNAs, the transition rate is greater with a transition/transversion bias of 3.13. Phylogenetic analysis of cp tRNA suggested that the type I introns in different taxa (including eukaryote organisms and cyanobacteria) share the conserved sequences "U-U-x2-C" and "U-x-G-x2-T," thereby indicating the diverse cyanobacterial origins of organelles. This detailed study of cp tRNAs in Adoxaceae may facilitate further investigations of the evolution, phylogeny, structure, and related functions of chloroplast tRNAs.

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.

Dimerization confers increased stability to nucleases in 5' halves from glycine and glutamic acid tRNAs

We have previously shown that 5′ halves from tRNA^Gly_GCC and tRNA^Glu_CUC are the most enriched small RNAs in the extracellular space of human cell lines, and especially in the non-vesicular fraction. Extracellular RNAs are believed to require protection by either encapsulation in vesicles or ribonucleoprotein complex formation. However, deproteinization of non-vesicular tRNA halves does not affect their retention in size-exclusion chromatography. Thus, we considered alternative explanations for their extracellular stability. In-silico analysis of the sequence of these tRNA-derived fragments showed that tRNA^Gly 5′ halves can form homodimers or heterodimers with tRNA^Glu 5′ halves. This capacity is virtually unique to glycine tRNAs. By analyzing synthetic oligonucleotides by size exclusion chromatography, we provide evidence that dimerization is possible in vitro. tRNA halves with single point substitutions preventing dimerization are degraded faster both in controlled nuclease digestion assays and after transfection in cells, showing that dimerization can stabilize tRNA halves against the action of cellular nucleases. Finally, we give evidence supporting dimerization of endogenous tRNA^Gly_GCC 5′ halves inside cells. Considering recent reports have shown that 5′ tRNA halves from Ala and Cys can form tetramers, our results highlight RNA intermolecular structures as a new layer of complexity in the biology of tRNA-derived fragments.

Glycyl-tRNA synthetase from Nanoarchaeum equitans: The first crystal structure of archaeal GlyRS and analysis of its tRNA glycylation

This study reports the X-ray crystallographic structure of the glycyl-tRNA synthetase (GlyRS) of Nanoarchaeum equitans - a hyperthermophilic archaeal species. This is the first archaeal GlyRS crystal structure elucidated. The GlyRS comprises an N-terminal catalytic domain and a C-terminal anticodon-binding domain with a long β-sheet inserted between these domains. An unmodified transcript of the wild-type N. equitans tRNA^Gly was successfully glycylated using GlyRS. Substitution of the discriminator base A73 of tRNA^Gly with any other nucleotide caused a significant decrease in glycylation activity. Mutational analysis of the second base-pair C2G71 of the acceptor stem of tRNA^Gly elucidated the importance of the base-pair, especially G71, as an identity element for recognition by GlyRS. Glycylation assays using tRNA^Gly G71 substitution mutants and a GlyRS mutant where Arg223 is mutated to alanine strengthen the possibility that the carbonyl oxygen at position 6 of G71 would hydrogen-bond with the guanidine nitrogen of Arg223 in N. equitans GlyRS.

Recent Insights Into the Structure, Function, and Evolution of the RNA-Splicing Endonucleases

RNA-splicing endonuclease (EndA) cleaves out introns from archaeal and eukaryotic precursor (pre)-tRNA and is essential for tRNA maturation. In archaeal EndA, the molecular mechanisms underlying complex assembly, substrate recognition, and catalysis have been well understood. Recently, certain studies have reported novel findings including the identification of new subunit types in archaeal EndA structures, providing insights into the mechanism underlying broad substrate specificity. Further, metagenomics analyses have enabled the acquisition of numerous DNA sequences of EndAs and intron-containing pre-tRNAs from various species, providing information regarding the co-evolution of substrate specificity of archaeal EndAs and tRNA genetic diversity, and the evolutionary pathway of archaeal and eukaryotic EndAs. Although the complex structure of the heterothermic form of eukaryotic EndAs is unknown, previous reports regarding their functions indicated that mutations in human EndA cause neurological disorders including pontocerebellar hypoplasia and progressive microcephaly, and yeast EndA significantly cleaves mitochondria-localized mRNA encoding cytochrome b mRNA processing 1 (Cpb1) for mRNA maturation. This mini-review summarizes the aforementioned results, discusses their implications, and offers my personal opinion regarding future directions for the analysis of the structure and function of EndAs.

Insights into RNA-processing pathways and associated RNA-degrading enzymes in Archaea

RNA-processing pathways are at the centre of regulation of gene expression. All RNA transcripts undergo multiple maturation steps in addition to covalent chemical modifications to become functional in the cell. This includes destroying unnecessary or defective cellular RNAs. In Archaea, information on mechanisms by which RNA species reach their mature forms and associated RNA-modifying enzymes are still fragmentary. To date, most archaeal actors and pathways have been proposed in light of information gathered from Bacteria and Eukarya. In this context, this review provides a state of the art overview of archaeal endoribonucleases and exoribonucleases that cleave and trim RNA species and also of the key small archaeal proteins that bind RNAs. Furthermore, synthetic up-to-date views of processing and biogenesis pathways of archaeal transfer and ribosomal RNAs as well as of maturation of stable small non-coding RNAs such as CRISPR RNAs, small C/D and H/ACA box guide RNAs, and other emerging classes of small RNAs are described. Finally, prospective post-transcriptional mechanisms to control archaeal messenger RNA quality and quantity are discussed.

Unique Archaeal Small RNAs

Advances in genome-wide sequence technologies allow for detailed insights into the complexity of RNA landscapes of organisms from all three domains of life. Recent analyses of archaeal transcriptomes identified interaction and regulation networks of noncoding RNAs in this understudied domain. Here, we review current knowledge of small, noncoding RNAs with important functions for the archaeal lifestyle, which often requires adaptation to extreme environments. One focus is RNA metabolism at elevated temperatures in hyperthermophilic archaea, which reveals elevated amounts of RNA-guided RNA modification and virus defense strategies. Genome rearrangement events result in unique fragmentation patterns of noncoding RNA genes that require elaborate maturation pathways to yield functional transcripts. RNA-binding proteins, e.g., L7Ae and LSm, are important for many posttranscriptional control functions of RNA molecules in archaeal cells. We also discuss recent insights into the regulatory potential of their noncoding RNA partners.

Structural Basis for tRNA Mimicry by a Bacterial Y RNA

Noncoding Y RNAs are present in both animal cells and many bacteria. In all species examined, Y RNAs tether the Ro60 protein to an effector protein to perform various cellular functions. Recently, a new Y RNA subfamily was identified in bacteria. Bioinformatic analyses of these YrlA (Y RNA-like A) RNAs predict that the effector-binding domain resembles tRNA. We present the structure of this domain, the overall folding of which is strikingly similar to canonical tRNAs. The tertiary interactions that are responsible for stabilizing tRNA are present in YrlA, making it a close tRNA mimic. However, YrlA lacks a free CCA end and contains a kink in the stem corresponding to the anticodon stem. Since nucleotides in the D and T stems are conserved among YrlAs, they may be an interaction site for an unknown factor. Our experiments identify YrlA RNAs as a new class of tRNA mimics.

The RNA-splicing endonuclease from the euryarchaeaon Methanopyrus kandleri is a heterotetramer with constrained substrate specificity

Four different types (α4, α'2, (αβ)2 and ϵ2) of RNA-splicing endonucleases (EndAs) for RNA processing are known to exist in the Archaea. Only the (αβ)2 and ϵ2 types can cleave non-canonical introns in precursor (pre)-tRNA. Both enzyme types possess an insert associated with a specific loop, allowing broad substrate specificity in the catalytic α units. Here, the hyperthermophilic euryarchaeon Methanopyrus kandleri (MKA) was predicted to harbor an (αβ)2-type EndA lacking the specific loop. To characterize MKA EndA enzymatic activity, we constructed a fusion protein derived from MKA α and β subunits (fMKA EndA). In vitro assessment demonstrated complete removal of the canonical bulge-helix-bulge (BHB) intron structure from MKA pre-tRNAAsn. However, removal of the relaxed BHB structure in MKA pre-tRNAGlu was inefficient compared to crenarchaeal (αβ)2 EndA, and the ability to process the relaxed intron within mini-helix RNA was not detected. fMKA EndA X-ray structure revealed a shape similar to that of other EndA types, with no specific loop. Mapping of EndA types and their specific loops and the tRNA gene diversity among various Archaea suggest that MKA EndA is evolutionarily related to other (αβ)2-type EndAs found in the Thaumarchaeota, Crenarchaeota and Aigarchaeota but uniquely represents constrained substrate specificity.

Novel Genomic and Evolutionary Perspective of Cyanobacterial tRNAs

Transfer RNA (tRNA) plays a central role in protein synthesis and acts as an adaptor molecule between an mRNA and an amino acid. A tRNA has an L-shaped clover leaf-like structure and contains an acceptor arm, D-arm, D-loop, anti-codon arm, anti-codon loop, variable loop, Ψ-arm and Ψ-loop. All of these arms and loops are important in protein translation. Here, we aimed to delineate the genomic architecture of these arms and loops in cyanobacterial tRNA. Studies from tRNA sequences from 61 cyanobacterial species showed that, except for few tRNAs (tRNA^Asn, tRNA^Leu, tRNA^Gln, and tRNA^Met), all contained a G nucleotide at the 1st position in the acceptor arm. tRNA^Leu and tRNA^Met did not contain any conserved nucleotides at the 1st position whereas tRNA^Asn and tRNA^Gln contained a conserved U¹ nucleotide. In several tRNA families, the variable region also contained conserved nucleotides. Except for tRNA^Met and tRNA^Glu, all other tRNAs contained a conserved A nucleotide at the 1st position in the D-loop. The Ψ-loop contained a conserved U¹-U²-C³-x-A⁵-x-U⁷ sequence, except for tRNA^Gly, tRNA^Ala, tRNA^Val, tRNA^Phe, tRNA^Thr, and tRNA^Gln in which the U⁷ nucleotide was not conserved. However, in tRNA^Asp, the U⁷ nucleotide was substituted with a C⁷ nucleotide. Additionally, tRNA^Arg, tRNA^Gly, and tRNA^Lys of cyanobacteria contained a group I intron within the anti-codon loop region. Maximum composite likelihood study on the transition/transversion of cyanobacterial tRNA revealed that the rate of transition was higher than the rate of transversion. An evolutionary tree was constructed to understand the evolution of cyanobacterial tRNA and analyses revealed that cyanobacterial tRNA may have evolved polyphyletically with high rate of gene loss.

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.

Binding Properties of Split tRNA to the C-terminal Domain of Methionyl-tRNA Synthetase of Nanoarchaeum equitans

The C-terminal domain of methionyl-tRNA synthetase (MetRS-C) from Nanoarchaeum equitans is homologous to a tRNA-binding protein consisting of 111 amino acids (Trbp111) from Aquifex aeolicus. The crystal structure of MetRS-C showed that it existed as a homodimer, and that each monomer possessed an oligonucleotide/oligosaccharide-binding fold (OB-fold). Analysis using a quartz crystal microbalance indicated that MetRS-C freshly isolated from N. equitans was bound to tRNA. However, binding of the split 3'-half tRNA species was stronger than that of the 5'-half species. The T-loop and the 3'-end regions of the split 3'-half tRNA were found to be responsible for the binding. The minimum structure for binding to MetRS-C might be a minihelix-like stem-loop with single-stranded 3'-terminus. After successive duplications of such a small hairpin structure with the assistance of a Trbp-like structure, the interaction of the T-loop region of the 3'-half with a Trbp-like structure could have been evolutionarily replaced by RNA-RNA interactions, along with many combinational tertiary interactions, to form the modern tRNA structure.

Evolutionary Insights into RNA trans-Splicing in Vertebrates

Pre-RNA splicing is an essential step in generating mature mRNA. RNA trans-splicing combines two separate pre-mRNA molecules to form a chimeric non-co-linear RNA, which may exert a function distinct from its original molecules. Trans-spliced RNAs may encode novel proteins or serve as noncoding or regulatory RNAs. These novel RNAs not only increase the complexity of the proteome but also provide new regulatory mechanisms for gene expression. An increasing amount of evidence indicates that trans-splicing occurs frequently in both physiological and pathological processes. In addition, mRNA reprogramming based on trans-splicing has been successfully applied in RNA-based therapies for human genetic diseases. Nevertheless, clarifying the extent and evolution of trans-splicing in vertebrates and developing detection methods for trans-splicing remain challenging. In this review, we summarize previous research, highlight recent advances in trans-splicing, and discuss possible splicing mechanisms and functions from an evolutionary viewpoint.

A Comprehensive tRNA Genomic Survey Unravels the Evolutionary History of tRNA Arrays in Prokaryotes

Considering the importance of tRNAs in the translation machinery, scant attention has been paid to tRNA array units defined as genomic regions containing at least 20 tRNA genes with a minimal tRNA gene density of two tRNA genes per kilobase. Our analysis of Acidithiobacillus ferrivorans CF27 and Acidithiobacillus ferrooxidans ATCC 23270(T) genomes showed that both display a tRNA array unit with syntenic conservation which mainly contributed to the tRNA gene redundancy in these two organisms. Our investigations into the occurrence and distribution of tRNA array units revealed that 1) this tRNA organization is limited to few phyla and mainly found in Gram-positive bacteria; and 2) the presence of tRNA arrays favors the redundancy of tRNA genes, in particular those encoding the core tRNA isoacceptors. Finally, comparative array organization revealed that tRNA arrays were acquired through horizontal gene transfer (from Firmicutes or unknown donor), before being subjected to tRNA rearrangements, deletions, and duplications. In Bacilli, the most parsimonious evolutionary history involved two common ancestors and the acquisition of their arrays arose late in evolution, in the genera branches. Functional roles of the array units in organism lifestyle, selective genetic advantage and translation efficiency, as well as the evolutionary advantages of organisms harboring them were proposed. Our study offers new insight into the structural organization and evolution of tRNA arrays in prokaryotic organisms.

AccessFold: predicting RNA-RNA interactions with consideration for competing self-structure

MOTIVATION

There are numerous examples of RNA-RNA complexes, including microRNA-mRNA and small RNA-mRNA duplexes for regulation of translation, guide RNA interactions with target RNA for post-transcriptional modification and small nuclear RNA duplexes for splicing. Predicting the base pairs formed between two interacting sequences remains difficult, at least in part because of the competition between unimolecular and bimolecular structure.

RESULTS

Two algorithms were developed for improved prediction of bimolecular RNA structure that consider the competition between self-structure and bimolecular structure. These algorithms utilize two novel approaches to evaluate accessibility: free energy density minimization and pseudo-energy minimization. Free energy density minimization minimizes the folding free energy change per nucleotide involved in an intermolecular secondary structure. Pseudo-energy minimization (called AccessFold) minimizes the sum of free energy change and a pseudo-free energy penalty for bimolecular pairing of nucleotides that are unlikely to be accessible for bimolecular structure. The pseudo-free energy, derived from unimolecular pairing probabilities, is applied per nucleotide in bimolecular pairs, and this approach is able to predict binding sites that are split by unimolecular structures. A benchmark set of 17 bimolecular RNA structures was assembled to assess structure prediction. Pseudo-energy minimization provides a statistically significant improvement in sensitivity over the method that was found in a benchmark to be the most accurate previously available method, with an improvement from 36.8% to 57.8% in mean sensitivity for base pair prediction.

AVAILABILITY AND IMPLEMENTATION

Pseudo-energy minimization is available for download as AccessFold, under an open-source license and as part of the RNAstructure package, at: http://rna.urmc.rochester.edu/RNAstructure.html

CONTACT

david_mathews@urmc.rochester.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Alternative genetic code for amino acids and transfer RNA revisited

The genetic code is highly conserved among all organisms and its evolution is thought to be strictly limited. However, an increasing number of studies have reported non-standard codes in prokaryotic and eukaryotic genomes. Most of these deviations from the standard code are attributable to tRNA changes relating to, for example, codon/anticodon base pairing and tRNA/aminoacyl-tRNA synthetase recognition. In this review, we focus on tRNA, a key molecule in the translation of the genetic code, and summarize the most recently published information on the evolutionary divergence of the tRNAs. Surprisingly, although higher eukaryotes, such as the nematode (worm), utilize the standard genetic code, newly identified nematode-specific tRNAs (nev-tRNAs) translate nucleotides in a manner that transgresses the code. Furthermore, a variety of additional functions of tRNAs, beyond their translation of the genetic code, have emerged rapidly. We also review these intriguing new aspects of tRNA, which have potential impacts on translational control, RNA silencing, antibiotic resistance, RNA biosynthesis, and transcriptional regulation.

Precise mapping and dynamics of tRNA-derived fragments (tRFs) in the development of Triops cancriformis (tadpole shrimp)

Background

In a deep sequencing analysis of small RNAs prepared from a living fossil, the tadpole shrimp Triops cancriformis, a 32-nt small RNA was specifically detected in the adult stage. A nucleotide sequence comparison between the 32-nt small RNA and predicted tRNA sequences in the draft nuclear genomic DNA showed that the small RNA was derived from tRNA^Gly(GCC). To determine the overall features of the tRNA-derived fragments (tRFs) of T. cancriformis, the small RNA sequences in each of the six developmental stages (egg, 1st − 4th instar larvae, and adult) were compared with the mitochondrial and nuclear tRNA sequences.

Results

We found that the tRFs were derived from mitochondrial and nuclear tRNAs corresponding to 16 and 39 anticodons, respectively. The total read number of nuclear tRFs was approximately 400 times larger than the number of mitochondrial tRFs. Interestingly, the main regions in each parental tRNA from which these tRFs were derived differed, depending on the parental anticodon. Mitochondrial tRF^Ser(GCU)s were abundantly produced from the 5’ half regions of the parental tRNA, whereas mitochondrial tRF^Val(UAC)s were mainly produced from the 3’ end regions. Highly abundant nuclear tRFs, tRF^Gly(GCC)s, tRF^Gly(CCC)s, tRF^Glu(CUC)s, and tRF^Lys(CUU)s were derived from the 5’ half regions of the parental tRNAs. Further analysis of the tRF read counts in the individual developmental stages suggested that the expression of mitochondrial and nuclear tRFs differed during the six stages. Based on these data, we precisely summarized the positions of the tRFs in their parental tRNAs and their expression changes during development.

Conclusions

Our results reveal the entire dynamics of the tRFs from both the nuclear and mitochondrial genomes of T. cancriformis and indicate that the majority of tRFs in the cell are derived from nuclear tRNAs. This study provides the first examples of developmentally expressed mitochondrial tRFs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0245-5) contains supplementary material, which is available to authorized users.

Capture, Unfolding, and Detection of Individual tRNA Molecules Using a Nanopore Device

Transfer RNAs (tRNA) are the most common RNA molecules in cells and have critical roles as both translators of the genetic code and regulators of protein synthesis. As such, numerous methods have focused on studying tRNA abundance and regulation, with the most widely used methods being RNA-seq and microarrays. Though revolutionary to transcriptomics, these assays are limited by an inability to encode tRNA modifications in the requisite cDNA. These modifications are abundant in tRNA and critical to their function. Here, we describe proof-of-concept experiments where individual tRNA molecules are examined as linear strands using a biological nanopore. This method utilizes an enzymatically ligated synthetic DNA adapter to concentrate tRNA at the lipid bilayer of the nanopore device and efficiently denature individual tRNA molecules, as they are pulled through the α-hemolysin (α-HL) nanopore. Additionally, the DNA adapter provides a loading site for ϕ29 DNA polymerase (ϕ29 DNAP), which acts as a brake on the translocating tRNA. This increases the dwell time of adapted tRNA in the nanopore, allowing us to identify the region of the nanopore signal that is produced by the translocating tRNA itself. Using adapter-modified Escherichia coli tRNA^fMet and tRNA^Lys, we show that the nanopore signal during controlled translocation is dependent on the identity of the tRNA. This confirms that adapter-modified tRNA can translocate end-to-end through nanopores and provide the foundation for future work in direct sequencing of individual transfer RNA with a nanopore-based device.

Pyrobaculum yellowstonensis Strain WP30 Respires on Elemental Sulfur and/or Arsenate in Circumneutral Sulfidic Geothermal Sediments of Yellowstone National Park

Thermoproteales (phylum Crenarchaeota) populations are abundant in high-temperature (>70°C) environments of Yellowstone National Park (YNP) and are important in mediating the biogeochemical cycles of sulfur, arsenic, and carbon. The objectives of this study were to determine the specific physiological attributes of the isolate Pyrobaculum yellowstonensis strain WP30, which was obtained from an elemental sulfur sediment (Joseph's Coat Hot Spring [JCHS], 80°C, pH 6.1, 135 μM As) and relate this organism to geochemical processes occurring in situ. Strain WP30 is a chemoorganoheterotroph and requires elemental sulfur and/or arsenate as an electron acceptor. Growth in the presence of elemental sulfur and arsenate resulted in the formation of thioarsenates and polysulfides. The complete genome of this organism was sequenced (1.99 Mb, 58% G+C content), revealing numerous metabolic pathways for the degradation of carbohydrates, amino acids, and lipids. Multiple dimethyl sulfoxide-molybdopterin (DMSO-MPT) oxidoreductase genes, which are implicated in the reduction of sulfur and arsenic, were identified. Pathways for the de novo synthesis of nearly all required cofactors and metabolites were identified. The comparative genomics of P. yellowstonensis and the assembled metagenome sequence from JCHS showed that this organism is highly related (∼95% average nucleotide sequence identity) to in situ populations. The physiological attributes and metabolic capabilities of P. yellowstonensis provide an important foundation for developing an understanding of the distribution and function of these populations in YNP.

Evolution of small guide RNA genes in hyperthermophilic archaea

Profiling the RNA production in hyperthermophilic archaea revealed an abundance of small RNA-guided processes near the upper temperature limit of life. Archaea utilize the base-pairing ability of RNA guide sequences to target ribosomal RNAs, transfer RNAs, messenger RNAs, and viral genomes. Cellular processes that are guided by small RNAs include the modification of RNA molecules, trans-splicing, gene regulation, and RNA and DNA degradation. Here, a brief overview of our knowledge on small guide RNA genes in archaeal genomes is provided and examples of their putative roles in genome evolution are described.

Conservation and losses of non-coding RNAs in avian genomes

Here we present the results of a large-scale bioinformatics annotation of non-coding RNA loci in 48 avian genomes. Our approach uses probabilistic models of hand-curated families from the Rfam database to infer conserved RNA families within each avian genome. We supplement these annotations with predictions from the tRNA annotation tool, tRNAscan-SE and microRNAs from miRBase. We identify 34 lncRNA-associated loci that are conserved between birds and mammals and validate 12 of these in chicken. We report several intriguing cases where a reported mammalian lncRNA, but not its function, is conserved. We also demonstrate extensive conservation of classical ncRNAs (e.g., tRNAs) and more recently discovered ncRNAs (e.g., snoRNAs and miRNAs) in birds. Furthermore, we describe numerous “losses” of several RNA families, and attribute these to either genuine loss, divergence or missing data. In particular, we show that many of these losses are due to the challenges associated with assembling avian microchromosomes. These combined results illustrate the utility of applying homology-based methods for annotating novel vertebrate genomes.

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.

Experimental confirmation of a whole set of tRNA molecules in two archaeal species

Based on the genomic sequences for most archaeal species, only one tRNA gene (isodecoder) is predicted for each triplet codon. This observation promotes analysis of a whole set of tRNA molecules and actual splicing patterns of interrupted tRNA in one organism. The entire genomic sequences of two Creanarchaeota, Aeropyrum pernix and Sulfolobus tokodaii, were determined approximately 15 years ago. In these genome datasets, 47 and 46 tRNA genes were detected, respectively. Among them, 14 and 24 genes, respectively, were predicted to be interrupted tRNA genes. To confirm the actual transcription of these predicted tRNA genes and identify the actual splicing patterns of the predicted interrupted tRNA molecules, RNA samples were prepared from each archaeal species and used to synthesize cDNA molecules with tRNA sequence-specific primers. Comparison of the nucleotide sequences of cDNA clones representing unspliced and spliced forms of interrupted tRNA molecules indicated that some introns were located at positions other than one base 3' from anticodon region and that bulge-helix-bulge structures were detected around the actual splicing sites in each interrupted tRNA molecule. Whole-set analyses of tRNA molecules revealed that the archaeal tRNA splicing mechanism may be essential for efficient splicing of all tRNAs produced from interrupted tRNA genes in these archaea.

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.

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.

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.

Avatar pre-tRNAs help elucidate the properties of tRNA-splicing endonucleases that produce tRNA from permuted genes

Unusual tRNA genes, found in some algae, have their mature terminal 3' portion in front of their 5' portion in the genome. The transcripts from such genes must be cleaved by a pre-tRNA endonuclease to form a functional tRNA. We present a mechanism for the generation of "corrected" tRNAs from such a "permuted" pre-tRNA configuration. We used two avatar (av) or model pre-tRNAs and two splicing endonucleases with distinct mechanisms of recognition of the pre-tRNA. The splicing results are compatible with an evolutionary route in which permuted genes result from a duplication event followed by DNA rearrangement. The model pre-tRNAs permit description of the features that a transcript, derived from a rearranged duplicated gene, must have to give rise to functional tRNA. The two tRNA endonucleases are a eukaryal enzyme that normally acts in a mature domain-dependent mode and an archaeal enzyme that acts in a mature domain-independent mode. Both av pre-tRNAs are able to fold into two conformations: 1 and 2. We find that only conformation 2 can yield a corrected functional tRNA. This result is consistent with contemporary algae representing snapshots of different evolutionary stages, with duplicated genes preceding recombinatorial events generating a permutated gene. In a scenario elucidated by the use of the av pre-tRNAs, algal permuted tRNA genes could have further lost one of two mature domains, eliminating steric problems for the algal tRNA endonuclease, which remains a typical eukaryal enzyme capable of correcting the permuted transcript to a functional tRNA.

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 (α′₂, α₄ and α₂β₂) 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 ε₂ homodimer and has broad substrate specificity like the α₂β₂ 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 ε₂ 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 ε₂ and α₂β₂ EndAs were separately acquired through different evolutionary processes.

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.

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

Class II transfer RNAs (tRNAs), including tRNA^Leu and tRNA^Ser, 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 tRNA^Leu. In vitro aminoacylation assays showed that nev-tRNA^Gly and nev-tRNA^Ile 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-tRNA^Gly 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.

The UCSC Archaeal Genome Browser: 2012 update

The UCSC Archaeal Genome Browser (http://archaea.ucsc.edu) offers a graphical web-based resource for exploration and discovery within archaeal and other selected microbial genomes. By bringing together existing gene annotations, gene expression data, multiple-genome alignments, pre-computed sequence comparisons and other specialized analysis tracks, the genome browser is a powerful aggregator of varied genomic information. The genome browser environment maintains the current look-and-feel of the vertebrate UCSC Genome Browser, but also integrates archaeal and bacterial-specific tracks with a few graphic display enhancements. The browser currently contains 115 archaeal genomes, plus 31 genomes of viruses known to infect archaea. Some of the recently developed or enhanced tracks visualize data from published high-throughput RNA-sequencing studies, the NCBI Conserved Domain Database, sequences from pre-genome sequencing studies, predicted gene boundaries from three different protein gene prediction algorithms, tRNAscan-SE gene predictions with RNA secondary structures and CRISPR locus predictions. We have also developed a companion resource, the Archaeal COG Browser, to provide better search and display of arCOG gene function classifications, including their phylogenetic distribution among available archaeal genomes.

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: α₄, homodimeric: α₂, and heterotetrameric: (αβ)₂] have been identified, each representing different substrate specificity during the tRNA intron cleavage. Here, we discovered a fourth type of archaeal tRNA splicing endonuclease (ε₂) 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 ε₂ endonuclease cleaves both canonical and non-canonical bulge–helix–bulge motifs, similar to that of (αβ)₂ 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 (αβ)₂ endonuclease. Thus, the discovery of ε₂ endonuclease in an archaeon deeply branched within Euryarchaeota represents a new example of the coevolution of tRNA and their processing enzymes.