Tri-split tRNA is a transfer RNA made from 3 transcripts that provides insight into the evolution of fragmented tRNAs 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.

tREPs-A New Class of Functional tRNA-Encoded Peptides

We asked if transfer RNA (tRNA) ever got an opportunity of translating its own sequence during evolution, what would have been the function of such tRNA-encoded peptides (tREPs)? If not, could one artificially synthesize tREPs to study the corresponding functional outcomes? Here, we report a novel, first-in-the-class, chemically synthesized tREP-18 molecule originating from the Escherichia coli tRNA sequence showing potent antileishmanial property. As a first step, E. coli tRNAs were computationally translated into peptide sequence equivalents and a database of full-length hypothetical tREPs was created. The tREP sequences were sent into sequence, structure, and energy filters to narrow down potential peptides for experimental validation. Based on the functional predictions, tREPs were screened against antiparasitic targets, leading to the identification of tREP-18 as a potential antiparasitic peptide. The in vitro assay of chemically synthesized tREP-18 on the Ag83 strain of Leishmania donovani showed its potent antileishmanial property (IC50 value of 22.13 nM). The atomic force microscopy and scanning electron microscopy images indicated significant alteration in the cytoskeletal architecture of tREP-18-treated parasites. Also, tREP-18 seems to destabilize the mitochondrial membrane potential of parasites, disrupting their cellular integrity and leading to parasitic death. The cellular assays of the tREP-18 peptide on the BS12 strain, a clinical isolate of post-kala azar dermal leishmaniasis, demonstrated its significant efficacy at an IC50 value of 15 nM. The tREP-18 peptide showed a toxic effect on the amastigote stage of the parasite, showing macrophage pathogen clearance at a concentration of 22.5 nM. This study provides the proof of the concept of making a new class of functional peptides from tRNA sequences. It also opens a huge untapped tRNA-peptide space toward novel discoveries and applications. In the future, it would be interesting to perform tREP edits and redesign tREPs toward specific applications.

PltRNAdb: Plant transfer RNA database

Transfer RNAs (tRNAs) are intermediate-sized non-coding RNAs found in all organisms that help translate messenger RNA into protein. Recently, the number of sequenced plant genomes has increased dramatically. The availability of this extensive data greatly accelerates the study of tRNAs on a large scale. Here, 8,768,261 scaffolds/chromosomes containing 229,093 giga-base pairs representing whole-genome sequences of 256 plant species were analyzed to identify tRNA genes. As a result, 331,242 nuclear, 3,216 chloroplast, and 1,467 mitochondrial tRNA genes were identified. The nuclear tRNA genes include 275,134 tRNAs decoding 20 standard amino acids, 1,325 suppressor tRNAs, 6,273 tRNAs with unknown isotypes, 48,475 predicted pseudogenes, and 37,873 tRNAs with introns. Efforts also extended to the creation of PltRNAdb (https://bioinformatics.um6p.ma/PltRNAdb/index.php), a data source for tRNA genes from 256 plant species. PltRNAdb website allows researchers to search, browse, visualize, BLAST, and download predicted tRNA genes. PltRNAdb will help improve our understanding of plant tRNAs and open the door to discovering the unknown regulatory roles of tRNAs in plant genomes.

The Clinical Significance of Transfer RNAs Present in Extracellular Vesicles

Extracellular vesicles (EVs) are important for intercellular signalling in multi-cellular organisms. However, the role of mature transfer RNAs (tRNAs) and tRNA fragments in EVs has yet to be characterised. This systematic review aimed to identify up-to-date literature on tRNAs present within human EVs and explores their potential clinical significance in health and disease. A comprehensive and systematic literature search was performed, and the study was conducted in accordance with PRISMA guidelines. Electronic databases MEDLINE and EMBASE were searched up until 1 January 2022. From 685 papers, 60 studies were identified for analysis. The majority of papers reviewed focussed on the role of EV tRNAs in cancers (31.7%), with numerous other conditions represented. Blood and cell lines were the most common EV sources, representing 85.9% of protocols used. EV isolation methods included most known methods, precipitation being the most common (49.3%). The proportion of EV tRNAs was highly variable, ranging between 0.04% to >95% depending on tissue source. EV tRNAs are present in a multitude of sources and show promise as disease markers in breast cancer, gastrointestinal cancers, and other diseases. EV tRNA research is an emerging field, with increasing numbers of papers highlighting novel methodologies for tRNA and tRNA fragment discovery.

Distinct Expansion of Group II Introns During Evolution of Prokaryotes and Possible Factors Involved in Its Regulation

Group II introns (G2Is) are ribozymes that have retroelement characteristics in prokaryotes. Although G2Is are suggested to have been an important evolutionary factor in the prokaryote-to-eukaryote transition, comprehensive analyses of these introns among the tens of thousands of prokaryotic genomes currently available are still limited. Here, we developed a bioinformatic pipeline that systematically collects G2Is and applied it to prokaryotic genomes. We found that in bacteria, 25% (447 of 1,790) of the total representative genomes had an average of 5.3 G2Is, and in archaea, 9% (28 of 296) of the total representative genomes had an average of 3.0 G2Is. The greatest number of G2Is per genome was 101 in Arthrospira platensis (phylum Cyanobacteriota). A comprehensive sequence analysis of the intron-encoded protein (IEP) in each G2I sequence was conducted and resulted in the addition of three new IEP classes (U1-U3) to the previous classification. This analysis suggested that about 30% of all IEPs are non-canonical IEPs. The number of G2Is per genome was defined almost at the phylum level, and at least in the following two phyla, Firmicutes, and Cyanobacteriota, the type of IEP was largely associated as a factor in the G2I increase, i.e., there was an explosive increase in G2Is with bacterial C-type IEPs, mainly in the phylum Firmicutes, and in G2Is with CL-type IEPs, mainly in the phylum Cyanobacteriota. We also systematically analyzed the relationship between genomic signatures and the mechanism of these increases in G2Is. This is the first study to systematically characterize G2Is in the prokaryotic phylogenies.

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.

Large-Scale Molecular Evolutionary Analysis Uncovers a Variety of Polynucleotide Kinase Clp1 Family Proteins in the Three Domains of Life

Clp1, a polyribonucleotide 5′-hydroxyl kinase in eukaryotes, is involved in pretRNA splicing and mRNA 3′-end formation. Enzymes similar in amino acid sequence to Clp1, Nol9, and Grc3, are present in some eukaryotes and are involved in prerRNA processing. However, our knowledge of how these Clp1 family proteins evolved and diversified is limited. We conducted a large-scale molecular evolutionary analysis of the Clp1 family proteins in all living organisms for which protein sequences are available in public databases. The phylogenetic distribution and frequencies of the Clp1 family proteins were investigated in complete genomes of Bacteria, Archaea and Eukarya. In total, 3,557 Clp1 family proteins were detected in the three domains of life, Bacteria, Archaea, and Eukarya. Many were from Archaea and Eukarya, but a few were found in restricted, phylogenetically diverse bacterial species. The domain structures of the Clp1 family proteins also differed among the three domains of life. Although the proteins were, on average, 555 amino acids long (range, 196–2,728), 122 large proteins with >1,000 amino acids were detected in eukaryotes. These novel proteins contain the conserved Clp1 polynucleotide kinase domain and various other functional domains. Of these proteins, >80% were from Fungi or Protostomia. The polyribonucleotide kinase activity of Thermus scotoductus Clp1 (Ts-Clp1) was characterized experimentally. Ts-Clp1 preferentially phosphorylates single-stranded RNA oligonucleotides (Km value for ATP, 2.5 µM), or single-stranded DNA at higher enzyme concentrations. We propose a comprehensive assessment of the diversification of the Clp1 family proteins and the molecular evolution of their functional domains.

Impact of intron removal from tRNA genes on Saccharomyces cerevisiae

In eukaryotes and archaea, tRNA genes frequently contain introns, which are removed during maturation. However, biological roles of tRNA introns remain elusive. Here, we constructed a complete set of Saccharomyces cerevisiae strains in which the introns were removed from all the synonymous genes encoding 10 different tRNA species. All the intronless strains were viable, but the tRNA^Phe_GAA and tRNA^Tyr_GUA intronless strains displayed slow growth, cold sensitivity and defective growth under respiratory conditions, indicating physiological importance of certain tRNA introns. Northern analyses revealed that removal of the introns from genes encoding three tRNAs reduced the amounts of the corresponding mature tRNAs, while it did not affect aminoacylation. Unexpectedly, the tRNA^Leu_CAA intronless strain showed reduced 5.8S rRNA levels and abnormal nucleolar morphology. Because pseudouridine (Ψ) occurs at position 34 of the tRNA^Ile_UAU anticodon in an intron-dependent manner, tRNA^Ile_UAU in the intronless strain lost Ψ34. However, in a portion of tRNA^Ile_UAU population, position 34 was converted into 5-carbamoylmethyluridine (ncm⁵U), which could reduce decoding fidelity. In summary, our results demonstrate that, while introns are dispensable for cell viability, some introns have diverse roles, such as ensuring proper growth under various conditions and controlling the appropriate anticodon modifications for accurate pairing with the codon.

Identified Hybrid tRNA Structure Genes in Archaeal Genome

Background

In Archaea, previous studies have revealed the presence of multiple intron-containing tRNAs and split tRNAs. The full unexpurgated analysis of archaeal tRNA genes remains a challenging task in the field of bioinformatics, because of the presence of various types of hidden tRNA genes in archaea. Here, we suggested a computational method that searched for widely separated genes encoding tRNA halves to generate suppressive variants of missing tRNAs.

Objectives

The exploration of tRNA genes from a genome with varying hypotheses, among all three domain of life (eukaryotes, bacteria and archaea), has been rapidly identified in different ways in the field of bioinformatics. Like eukaryotic tRNA genes, it has been established that two separated regions of the coding sequence of a tRNA gene are essential and sufficient for promotion of transcription. Our objective is to find out the two essential regions in the genome sequence which comprises two halves of the hidden tRNAs.

Material and Methods

Considering the existence of split tRNA genes widely separated throughout the genome, we developed our tRNA search algorithm to predict such separated tRNA genes by searching both a conserved terminal 5'- and 3'-motif of tRNA in agreement with the split hypothesis on the basis of cloverleaf prediction and precise insilico determination of bulge-helix-bulge secondary structure at the splice sites.

Results

By a comprehensive search for all kinds of missing tRNA genes, we have constructed hybrid tRNA genes containing one essential region from tDNA (XYZ) and the other from tDNA (ABC), both from same species in the archaea. We have also found, this type of hybrid tRNA genes are identified in the different species of the archaea (XYZ ASN, ARG and MET; ABC ASP,SER, ARG and PRO).These hybrid split tRNA share a common structural motif called bulge-helix-bulge (BHB) a more relaxed bulge-helix loop (BHL), at the leader exon boundary and suggested to be evolutionary interrelated.

Conclusions

Analysis of the complete genome sequences of Metallosphaera sedula DSM 5348, Desulfurococcus kamchatkensis 1221n and Ignicoccus hospitalis KIN4/I in archaea by our algorithm revealed that a number of hybrid tRNAs are constructed from different tDNAs . Asymmetric combination of 5' and 3' tRNA halves may have generated the diversity of tRNA molecules. Our study of hybrid tRNA genes will provide a new molecular basis for upcoming tRNA studies.

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.

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.

The importance of designating type material for uncultured taxa

Naming of uncultured Bacteria and Archaea is often inconsistent with the International Code of Nomenclature of Prokaryotes. The recent practice of proposing names for higher taxa without designation of lower ranks and nomenclature types is one of the most important inconsistencies that needs to be addressed to avoid nomenclatural instability. The Code requires names of higher taxa up to the rank of class to be derived from the type genus name, with a proposal pending to formalise this requirement for the rank of phylum. Designation of nomenclature types is crucial for providing priority to names and ensures their uniqueness and stability. However, only legitimate names proposed for axenic cultures can be used for this purpose. Candidatus names reserved for taxa lacking cultured representatives may be granted this right if recent proposals to use genome sequences as type material are endorsed, thereby allowing the Code to be fully applied to lineages represented by metagenome-assembled genomes (MAGs) or single amplified genomes (SAGs). Genome quality standards need to be considered to ensure unambiguous assignment of type material. Here, we illustrate the recommended practice by proposing nomenclature type material for four major uncultured prokaryotic lineages based on high-quality MAGs in accordance with the Code.

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.

The emerging complexity of the tRNA world: mammalian tRNAs beyond protein synthesis

The discovery of the genetic code and tRNAs as decoders of the code transformed life science. However, after establishing the role of tRNAs in protein synthesis, the field moved to other parts of the RNA world. Now, tRNA research is blooming again, with demonstration of the involvement of tRNAs in various other pathways beyond translation and in adapting translation to environmental cues. These roles are linked to the presence of tRNA sequence variants known as isoacceptors and isodecoders, various tRNA base modifications, the versatility of protein binding partners and tRNA fragmentation events, all of which collectively create an incalculable complexity. This complexity provides a vast repertoire of tRNA species that can serve various functions in cellular homeostasis and in adaptation of cellular functions to changing environments, and it likely arose from the fundamental role of RNAs in early evolution.

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.

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.

A Post-Genomic View of the Ecophysiology, Catabolism and Biotechnological Relevance of Sulphate-Reducing Prokaryotes

Dissimilatory sulphate reduction is the unifying and defining trait of sulphate-reducing prokaryotes (SRP). In their predominant habitats, sulphate-rich marine sediments, SRP have long been recognized to be major players in the carbon and sulphur cycles. Other, more recently appreciated, ecophysiological roles include activity in the deep biosphere, symbiotic relations, syntrophic associations, human microbiome/health and long-distance electron transfer. SRP include a high diversity of organisms, with large nutritional versatility and broad metabolic capacities, including anaerobic degradation of aromatic compounds and hydrocarbons. Elucidation of novel catabolic capacities as well as progress in the understanding of metabolic and regulatory networks, energy metabolism, evolutionary processes and adaptation to changing environmental conditions has greatly benefited from genomics, functional OMICS approaches and advances in genetic accessibility and biochemical studies. Important biotechnological roles of SRP range from (i) wastewater and off gas treatment, (ii) bioremediation of metals and hydrocarbons and (iii) bioelectrochemistry, to undesired impacts such as (iv) souring in oil reservoirs and other environments, and (v) corrosion of iron and concrete. Here we review recent advances in our understanding of SRPs focusing mainly on works published after 2000. The wealth of publications in this period, covering many diverse areas, is a testimony to the large environmental, biogeochemical and technological relevance of these organisms and how much the field has progressed in these years, although many important questions and applications remain to be explored.

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.

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

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

Handling tRNA introns, archaeal way and eukaryotic way

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

tRNA gene diversity in the three domains of life

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

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

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

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.

Class-specific prediction of ncRNAs

Many RNA families, i.e., groups of homologous RNA genes, belong to RNA classes, such as tRNAs, snoRNAs, or microRNAs, that are characterized by common sequence motifs and/or common secondary structure features. The detection of new members of RNA classes, as well as the comprehensive annotation of genomes with members of RNA classes is a challenging task that goes beyond simple homology search. Computational methods addressing this problem typically use a three-tiered approach: In the first step an efficient and sensitive filter is employed. In the second step the candidate set is narrowed down using computationally expensive methods geared towards specificity. In the final step the hits are annotated with class-specific features and scored. Here we review the tools that are currently available for a diverse set of RNA classes.

The impact of genomics on research in diversity and evolution of archaea

Since the definition of archaea as a separate domain of life along with bacteria and eukaryotes, they have become one of the most interesting objects of modern microbiology, molecular biology, and biochemistry. Sequencing and analysis of archaeal genomes were especially important for studies on archaea because of a limited availability of genetic tools for the majority of these microorganisms and problems associated with their cultivation. Fifteen years since the publication of the first genome of an archaeon, more than one hundred complete genome sequences of representatives of different phylogenetic groups have been determined. Analysis of these genomes has expanded our knowledge of biology of archaea, their diversity and evolution, and allowed identification and characterization of new deep phylogenetic lineages of archaea. The development of genome technologies has allowed sequencing the genomes of uncultivated archaea directly from enrichment cultures, metagenomic samples, and even from single cells. Insights have been gained into the evolution of key biochemical processes in archaea, such as cell division and DNA replication, the role of horizontal gene transfer in the evolution of archaea, and new relationships between archaea and eukaryotes have been revealed.

Putative anticodons in mitochondrial tRNA sidearm loops: Pocketknife tRNAs?

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

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

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

Mapping the RNA-Seq trash bin: unusual transcripts in prokaryotic transcriptome sequencing data

Prokaryotic transcripts constitute almost always uninterrupted intervals when mapped back to the genome. Split reads, i.e., RNA-seq reads consisting of parts that only map to discontiguous loci, are thus disregarded in most analysis pipelines. There are, however, some well-known exceptions, in particular, tRNA splicing and circularized small RNAs in Archaea as well as self-splicing introns. Here, we reanalyze a series of published RNA-seq data sets, screening them specifically for non-contiguously mapping reads. We recover most of the known cases together with several novel archaeal ncRNAs associated with circularized products. In Eubacteria, only a handful of interesting candidates were obtained beyond a few previously described group I and group II introns. Most of the atypically mapping reads do not appear to correspond to well-defined, specifically processed products. Whether this diffuse background is, at least in part, an incidental by-product of prokaryotic RNA processing or whether it consists entirely of technical artifacts of reverse transcription or amplification remains unknown.

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.

Diversity and roles of (t)RNA ligases

The discovery of discontiguous tRNA genes triggered studies dissecting the process of tRNA splicing. As a result, we have gained detailed mechanistic knowledge on enzymatic removal of tRNA introns catalyzed by endonuclease and ligase proteins. In addition to the elucidation of tRNA processing, these studies facilitated the discovery of additional functions of RNA ligases such as RNA repair and non-conventional mRNA splicing events. Recently, the identification of a new type of RNA ligases in bacteria, archaea, and humans closed a long-standing gap in the field of tRNA processing. This review summarizes past and recent findings in the field of tRNA splicing with a focus on RNA ligation as it preferentially occurs in archaea and humans. In addition to providing an integrated view of the types and phyletic distribution of RNA ligase proteins known to date, this survey also aims at highlighting known and potential accessory biological functions of RNA ligases.

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.

Unscrambling genetic information at the RNA level

Genomics aims at unraveling the blueprint of life; however, DNA sequence alone does not always reveal the proteins and structural RNAs encoded by the genome. The reason is that genetic information is often encrypted. Recognizing the logic of encryption, and understanding how living cells decode hidden information--at the level of DNA, RNA or protein--is challenging. RNA-level decryption includes topical RNA editing and more 'macroscopic' transcript rearrangements. The latter events involve the four types of introns recognized to date, notably spliceosomal, group I, group II, and archaeal/tRNA splicing. Intricate variants, such as alternative splicing and trans-splicing, have been reported for each intron type, but the biological significance has not always been confirmed. Novel RNA-level unscrambling processes were recently discovered in mitochondria of dinoflagellates and diplonemids, and potentially euglenids. These processes seem not to rely on known introns, and the corresponding molecular mechanisms remain to be elucidated.

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.

tRNA-like elements in Haloferax volcanii

All functional RNAs are generated from precursor molecules by a plethora of processing steps. The generation of mature RNA molecules by processing is an important layer of gene expression regulation catalysed by ribonucleases. Here, we analysed 5S rRNA processing in the halophilic Archaeon Haloferax volcanii. Earlier experiments showed that the 5S rRNA is cleaved at its 5' end by the endonuclease tRNase Z. Interestingly, a tRNA-like structure was identified upstream of the 5S rRNA that might be used as a processing signal. Here, we show that this tRNA-like element is indeed recognised as a processing signal by tRNase Z. Substrates containing mutations in the tRNA-like sequence are no longer processed, whereas a substrate containing a deletion in the 5S rRNA sequence is still cleaved. Therefore, an intact 5S rRNA structure is not required for processing. Further, we used bioinformatics analyses to identify additional sequences in Haloferax containing tRNA-like structures. This search resulted in the identification of all tRNAs, the tRNA-like structure upstream of the 5S RNA and 47 new tRNA-like structural elements. However, the in vitro processing of selected examples showed no cleavage of these newly identified elements. Thus, tRNA-like elements are not a general processing signal in Haloferax.

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.

Unassigned codons, nonsense suppression, and anticodon modifications in the evolution of the genetic code

The origin of the genetic code is a central open problem regarding the early evolution of life. Here, we consider two undeveloped but important aspects of possible scenarios for the evolutionary pathway of the translation machinery: the role of unassigned codons in early stages of the code and the incorporation of tRNA anticodon modifications. As the first codons started to encode amino acids, the translation machinery likely was faced with a large number of unassigned codons. Current molecular scenarios for the evolution of the code usually assume the very rapid assignment of all codons before all 20 amino acids became encoded. We show that the phenomenon of nonsense suppression as observed in current organisms allows for a scenario in which many unassigned codons persisted throughout most of the evolutionary development of the code. In addition, we demonstrate that incorporation of anticodon modifications at a late stage is feasible. The wobble rules allow a set of 20 tRNAs fully lacking anticodon modifications to encode all 20 canonical amino acids. These observations have implications for the biochemical plausibility of early stages in the evolution of the genetic code predating tRNA anticodon modifications and allow for effective translation by a relatively small and simple early tRNA set.

The intron of tRNA-TrpCCA is dispensable for growth and translation of Saccharomyces cerevisiae

A part of eukaryotic tRNA genes harbor an intron at one nucleotide 3' to the anticodon, so that removal of the intron is an essential processing step for tRNA maturation. While some tRNA introns have important roles in modification of certain nucleotides, essentiality of the tRNA intron in eukaryotes has not been tested extensively. This is partly because most of the eukaryotic genomes have multiple genes encoding an isoacceptor tRNA. Here, we examined whether the intron of tRNA-Trp(CCA) genes, six copies of which are scattered on the genome of yeast, Saccharomyces cerevisiae, is essential for growth or translation of the yeast in vivo. We devised a procedure to remove all of the tRNA introns from the yeast genome iteratively with marker cassettes containing both positive and negative markers. Using this procedure, we removed all the introns from the six tRNA-Trp(CCA) genes, and found that the intronless strain grew normally and expressed tRNA-Trp(CCA) in an amount similar to that of the wild-type genes. Neither incorporation of (35)S-labeled amino acids into a TCA-insoluble fraction nor the major protein pattern on SDS-PAGE/2D gel were affected by complete removal of the intron, while expression levels of some proteins were marginally affected. Therefore, the tRNA-Trp(CCA) intron is dispensable for growth and bulk translation of the yeast. This raises the possibility that some mechanism other than selective pressure from translational efficiency maintains the tRNA intron on the yeast genome.

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.

Identification of an entire set of tRNA molecules and characterization of cleavage sites of the intron-containing tRNA precursors in acidothermophilic crenarchaeon Sulfolobus tokodaii strain7

The acidothermophilic crenarchaeon, Sulfolobus tokodaii strain7, was isolated from a hot spring in Beppu, Kyushu, Japan. Whole genomic data of this microorganism indicated that among 46 putative tRNA genes identified, 24 were interrupted tRNA genes containing an intron. A sequence comparison between the cDNA sequences for unspliced and spliced tRNAs indicated that all predicted tRNAs were expressed and all intron portions were spliced in this microorganism. However, the actual cleavage site in the splicing process was not determined for 13 interrupted tRNAs because of the presence of the same nucleotides at both 5' and 3' border regions of each intron. The cleavage sites for all the introns, which were determined by an in vitro cleavage experiment with recombinant splicing endonuclease as well as cDNA sequencing of the spliced tRNAs, indicated that non-canonical BHB structure motifs were also recognized and processed by the splicing machinery in this organism. This is the first report to empirically determine the actual cleavage and splice sites of introns in the whole set of archaeal tRNA genes, and reassigns the exon-intron borders with a novel and more plausible non-canonical BHB structure.

[Virus, phage, transposon and their regulatory small non-coding RNAs]

Many reports have been accumulated describing not a few microRNAs (miRNAs) in eukaryotes target viral genomes, whereas a number of viruses also encode miRNA genes. These small RNAs play important roles on viral infection and their replication. In germ cells, another small RNA, piRNA is reported to repress endogenous transposons. Furthermore, CRISPR RNA target virus/phage genomes in both archaea and bacteria. Therefore, small RNA is deeply involved in a broad range of biological defense systems. This system may be applied not only to control replication of viruses or phages but also provide implication on regulating the growth of microorganisms including pathogenic bacteria.