Identification of BHB splicing motifs in intron-containing tRNAs from 18 archaea: evolutionary implications

Variation of tRNA modifications with and without intron dependency

tRNAs have recently gained attention for their novel regulatory roles in translation and for their diverse functions beyond translation. One of the most remarkable aspects of tRNA biogenesis is the incorporation of various chemical modifications, ranging from simple base or ribose methylation to more complex hypermodifications such as formation of queuosine and wybutosine. Some tRNAs are transcribed as intron-containing pre-tRNAs. While the majority of these modifications occur independently of introns, some are catalyzed in an intron-inhibitory manner, and in certain cases, they occur in an intron-dependent manner. This review focuses on pre-tRNA modification, including intron-containing pre-tRNA, in both intron-inhibitory and intron-dependent fashions. Any perturbations in the modification and processing of tRNAs may lead to a range of diseases and disorders, highlighting the importance of understanding these mechanisms in molecular biology and medicine.

Asgard archaea modulate potential methanogenesis substrates in wetland soil

The roles of Asgard archaea in eukaryogenesis and marine biogeochemical cycles are well studied, yet their contributions in soil ecosystems remain unknown. Of particular interest are Asgard archaeal contributions to methane cycling in wetland soils. To investigate this, we reconstructed two complete genomes for soil-associated Atabeyarchaeia, a new Asgard lineage, and a complete genome of Freyarchaeia, and predicted their metabolism in situ. Metatranscriptomics reveals expression of genes for [NiFe]-hydrogenases, pyruvate oxidation and carbon fixation via the Wood-Ljungdahl pathway. Also expressed are genes encoding enzymes for amino acid metabolism, anaerobic aldehyde oxidation, hydrogen peroxide detoxification and carbohydrate breakdown to acetate and formate. Overall, soil-associated Asgard archaea are predicted to include non-methanogenic acetogens, highlighting their potential role in carbon cycling in terrestrial environments.

Recognition and cleavage mechanism of intron-containing pre-tRNA by human TSEN endonuclease complex

Removal of introns from transfer RNA precursors (pre-tRNAs) occurs in all living organisms. This is a vital phase in the maturation and functionality of tRNA. Here we present a 3.2 Å-resolution cryo-EM structure of an active human tRNA splicing endonuclease complex bound to an intron-containing pre-tRNA. TSEN54, along with the unique regions of TSEN34 and TSEN2, cooperatively recognizes the mature body of pre-tRNA and guides the anticodon-intron stem to the correct position for splicing. We capture the moment when the endonucleases are poised for cleavage, illuminating the molecular mechanism for both 3' and 5' cleavage reactions. Two insertion loops from TSEN54 and TSEN2 cover the 3' and 5' splice sites, respectively, trapping the scissile phosphate in the center of the catalytic triad of residues. Our findings reveal the molecular mechanism for eukaryotic pre-tRNA recognition and cleavage, as well as the evolutionary relationship between archaeal and eukaryotic TSENs.

New insights into RNA processing by the eukaryotic tRNA splicing endonuclease

Through its role in intron cleavage, tRNA splicing endonuclease (TSEN) plays a critical function in the maturation of intron-containing pre-tRNAs. The catalytic mechanism and core requirement for this process is conserved between archaea and eukaryotes, but for decades, it has been known that eukaryotic TSENs have evolved additional modes of RNA recognition, which have remained poorly understood. Recent research identified new roles for eukaryotic TSEN, including processing or degradation of additional RNA substrates, and determined the first structures of pre-tRNA-bound human TSEN complexes. These recent discoveries have changed our understanding of how the eukaryotic TSEN targets and recognizes substrates. Here, we review these recent discoveries, their implications, and the new questions raised by these findings.

Intron-Dependent or Independent Pseudouridylation of Precursor tRNA Containing Atypical Introns in Cyanidioschyzon merolae

Eukaryotic precursor tRNAs (pre-tRNAs) often have an intron between positions 37 and 38 of the anticodon loop. However, atypical introns are found in some eukaryotes and archaea. In an early-diverged red alga Cyanidioschyzon merolae, the tRNA^Ile(UAU) gene contains three intron coding regions, located in the D-, anticodon, and T-arms. In this study, we focused on the relationship between the intron removal and formation of pseudouridine (Ψ), one of the most universally modified nucleosides. It had been reported that yeast Pus1 is a multiple-site-specific enzyme that synthesizes Ψ34 and Ψ36 in tRNA^Ile(UAU) in an intron-dependent manner. Unexpectedly, our biochemical experiments showed that the C. merolae ortholog of Pus1 pseudouridylated an intronless tRNA^Ile(UAU) and that the modification position was determined to be 55 which is the target of Pus4 but not Pus1 in yeast. Furthermore, unlike yeast Pus1, cmPus1 mediates Ψ modification at positions 34, 36, and/or 55 only in some specific intron-containing pre-tRNA^Ile(UAU) variants. cmPus4 was confirmed to be a single-site-specific enzyme that only converts U55 to Ψ, in a similar manner to yeast Pus4. cmPus4 did not catalyze the pseudouridine formation in pre-tRNAs containing an intron in the T-arm.

Transfer RNA processing - from a structural and disease perspective

Transfer RNAs (tRNAs) are highly structured non-coding RNAs which play key roles in translation and cellular homeostasis. tRNAs are initially transcribed as precursor molecules and mature by tightly controlled, multistep processes that involve the removal of flanking and intervening sequences, over 100 base modifications, addition of non-templated nucleotides and aminoacylation. These molecular events are intertwined with the nucleocytoplasmic shuttling of tRNAs to make them available at translating ribosomes. Defects in tRNA processing are linked to the development of neurodegenerative disorders. Here, we summarize structural aspects of tRNA processing steps with a special emphasis on intron-containing tRNA splicing involving tRNA splicing endonuclease and ligase. Their role in neurological pathologies will be discussed. Identification of novel RNA substrates of the tRNA splicing machinery has uncovered functions unrelated to tRNA processing. Future structural and biochemical studies will unravel their mechanistic underpinnings and deepen our understanding of neurological diseases.

Eukaryotic tRNA splicing - one goal, two strategies, many players

Transfer RNAs (tRNAs) are transcribed as precursor molecules that undergo several maturation steps before becoming functional for protein synthesis. One such processing mechanism is the enzyme-catalysed splicing of intron-containing pre-tRNAs. Eukaryotic tRNA splicing is an essential process since intron-containing tRNAs cannot fulfil their canonical function at the ribosome. Splicing of pre-tRNAs occurs in two steps: The introns are first excised by a tRNA-splicing endonuclease and the exons are subsequently sealed by an RNA ligase. An intriguing complexity has emerged from newly identified tRNA splicing factors and their interplay with other RNA processing pathways during the past few years. This review summarises our current understanding of eukaryotic tRNA splicing and the underlying enzyme machinery. We highlight recent structural advances and how they have shaped our mechanistic understanding of tRNA splicing in eukaryotic cells. A special focus lies on biochemically distinct strategies for exon-exon ligation in fungi versus metazoans.

The tRNA discriminator base defines the mutual orthogonality of two distinct pyrrolysyl-tRNA synthetase/tRNAPyl pairs in the same organism

Site-specific incorporation of distinct non-canonical amino acids into proteins via genetic code expansion requires mutually orthogonal aminoacyl-tRNA synthetase/tRNA pairs. Pyrrolysyl-tRNA synthetase (PylRS)/tRNAPyl pairs are ideal for genetic code expansion and have been extensively engineered for developing mutually orthogonal pairs. Here, we identify two novel wild-type PylRS/tRNAPyl pairs simultaneously present in the deep-rooted extremely halophilic euryarchaeal methanogen Candidatus Methanohalarchaeum thermophilum HMET1, and show that both pairs are functional in the model halophilic archaeon Haloferax volcanii. These pairs consist of two different PylRS enzymes and two distinct tRNAs with dissimilar discriminator bases. Surprisingly, these two PylRS/tRNAPyl pairs display mutual orthogonality enabled by two unique features, the A73 discriminator base of tRNAPyl2 and a shorter motif 2 loop in PylRS2. In vivo translation experiments show that tRNAPyl2 charging by PylRS2 is defined by the enzyme's shortened motif 2 loop. Finally, we demonstrate that the two HMET1 PylRS/tRNAPyl pairs can simultaneously decode UAG and UAA codons for incorporation of two distinct noncanonical amino acids into protein. This example of a single base change in a tRNA leading to additional coding capacity suggests that the growth of the genetic code is not yet limited by the number of identity elements fitting into the tRNA structure.

Archaeal tRNA-Splicing Endonuclease as an Effector for RNA Recombination and Novel Trans-Splicing Pathways in Eukaryotes

We have characterized a homodimeric tRNA endonuclease from the euryarchaeota Ferroplasma acidarmanus (FERAC), a facultative anaerobe which can grow at temperatures ranging from 35 to 42 °C. This enzyme, contrary to the eukaryal tRNA endonucleases and the homotetrameric Methanocaldococcus jannaschii (METJA) homologs, is able to cleave minimal BHB (bulge–helix–bulge) substrates at 30 °C. The expression of this enzyme in Schizosaccharomyces pombe (SCHPO) enables the use of its properties as effectors by inserting BHB motif introns into hairpin loops normally seen in mRNA transcripts. In addition, the FERAC endonuclease can create proteins with new functionalities through the recombination of protein domains.

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.

Splicing Endonuclease Is an Important Player in rRNA and tRNA Maturation in Archaea

In all three domains of life, tRNA genes contain introns that must be removed to yield functional tRNA. In archaea and eukarya, the first step of this process is catalyzed by a splicing endonuclease. The consensus structure recognized by the splicing endonuclease is a bulge-helix-bulge (BHB) motif which is also found in rRNA precursors. So far, a systematic analysis to identify all biological substrates of the splicing endonuclease has not been carried out. In this study, we employed CRISPRi to repress expression of the splicing endonuclease in the archaeon Haloferax volcanii to identify all substrates of this enzyme. Expression of the splicing endonuclease was reduced to 1% of its normal level, resulting in a significant extension of lag phase in H. volcanii growth. In the repression strain, 41 genes were down-regulated and 102 were up-regulated. As an additional approach in identifying new substrates of the splicing endonuclease, we isolated and sequenced circular RNAs, which identified excised introns removed from tRNA and rRNA precursors as well as from the 5' UTR of the gene HVO_1309. In vitro processing assays showed that the BHB sites in the 5' UTR of HVO_1309 and in a 16S rRNA-like precursor are processed by the recombinant splicing endonuclease. The splicing endonuclease is therefore an important player in RNA maturation in archaea.

Comprehensive analysis of the pre-ribosomal RNA maturation pathway in a methanoarchaeon exposes the conserved circularization and linearization mode in archaea

The ribosomal RNA (rRNA) genes are generally organized as an operon and cotranscribed into a polycistronic precursor; therefore, processing and maturation of pre-rRNAs are essential for ribosome biogenesis. However, rRNA maturation pathways of archaea, particularly of methanoarchaea, are scarcely known. Here, we thoroughly elucidated the maturation pathway of the rRNA operon (16S-tRNA^Ala-23S-tRNA^Cys-5S) in Methanolobus psychrophilus, one representative of methanoarchaea. Enzymatic assay demonstrated that EndA, a tRNA splicing endoribonuclease, cleaved bulge-helix-bulge (BHB) motifs buried in the processing stems of pre-16S and pre-23S rRNAs. Northern blot and quantitative PCR detected splicing-coupled circularization of pre-16S and pre-23S rRNAs, which accounted for 2% and 12% of the corresponding rRNAs, respectively. Importantly, endoribonuclease Nob1 was determined to linearize circular pre-16S rRNA at the mature 3' end so to expose the anti-Shine-Dalgarno sequence, while circular pre-23S rRNA was linearized at the mature 5' end by an unknown endoribonuclease. The resultant 5' and 3' extension in linearized pre-16S and pre-23S rRNAs were finally matured through 5'-3' and 3'-5' exoribonucleolytic trimming, respectively. Additionally, a novel processing pathway of endoribonucleolysis coupled with exoribonucleolysis was identified for the pre-5S rRNA maturation in this methanogen, which could be also conserved in most methanogenic euryarchaea. Based on evaluating the phylogenetic conservation of the key elements that are involved in circularization and linearization of pre-rRNA maturation, we predict that the rRNA maturation mode revealed here could be prevalent among archaea.

The function of KptA/Tpt1 gene - a minor review

Rapid response of uni- and multicellular organisms to environmental changes and their own growth is achieved through a series of molecular mechanisms, often involving modification of macromolecules, including nucleic acids, proteins and lipids. The ADP-ribosylation process has ability to modify these different macromolecules in cells, and is closely related to the biological processes, such as DNA replication, transcription, signal transduction, cell division, stress, microbial aging and pathogenesis. In addition, tRNA plays an essential role in the regulation of gene expression, as effector molecules, no-load tRNA affects the overall gene expression level of cells under some nutritional stress. KptA/Tpt1 is an essential phosphotransferase in the process of pre-tRNA splicing, releasing mature tRNA and participating in ADP-ribose. The objective of this review is concluding the gene structure, the evolution history and the function of KptA/Tpt1 from prokaryote to eukaryote organisms. At the same time, the results of promoter elements analysis were also shown in the present study. Moreover, the problems in the function of KptA/Tpt1 that have not been clarified at the present time are summarised, and some suggestions to solve those problems are given. This review presents no only a summary of clear function of KptA/Tpt1 in the process of tRNA splicing and ADP-ribosylation of organisms, but also gives some proposals to clarify unclear problems of it in the future.

Molecular determinants of metazoan tricRNA biogenesis

Mature tRNAs are generated by multiple post-transcriptional processing steps, which can include intron removal. Recently, we discovered a new class of circular non-coding RNAs in metazoans, called tRNA intronic circular (tric)RNAs. To investigate the mechanism of tricRNA biogenesis, we generated constructs that replace native introns of human and fruit fly tRNA genes with the Broccoli fluorescent RNA aptamer. Using these reporters, we identified cis-acting elements required for tricRNA formation in vivo. Disrupting a conserved base pair in the anticodon-intron helix dramatically reduces tricRNA levels. Although the integrity of this base pair is necessary for proper splicing, it is not sufficient. In contrast, strengthening weak bases in the helix also interferes with splicing and tricRNA production. Furthermore, we identified trans-acting factors important for tricRNA biogenesis, including several known tRNA processing enzymes such as the RtcB ligase and components of the TSEN endonuclease complex. Depletion of these factors inhibits Drosophila tRNA intron circularization. Notably, RtcB is missing from fungal genomes and these organisms normally produce linear tRNA introns. Here, we show that in the presence of ectopic RtcB, yeast lacking the tRNA ligase Rlg1/Trl1 are converted into producing tricRNAs. In summary, our work characterizes the major players in eukaryotic tricRNA biogenesis.

Mutations in the anticodon stem of tRNA cause accumulation and Met22-dependent decay of pre-tRNA in yeast

During tRNA maturation in yeast, aberrant pre-tRNAs are targeted for 3'-5' degradation by the nuclear surveillance pathway, and aberrant mature tRNAs are targeted for 5'-3' degradation by the rapid tRNA decay (RTD) pathway. RTD is catalyzed by the 5'-3' exonucleases Xrn1 and Rat1, which act on tRNAs with an exposed 5' end due to the lack of certain body modifications or the presence of destabilizing mutations in the acceptor stem, T-stem, or tRNA fold. RTD is inhibited by mutation of MET22, likely due to accumulation of the Met22 substrate adenosine 3',5' bis-phosphate, which inhibits 5'-3' exonucleases. Here we provide evidence for a new tRNA quality control pathway in which intron-containing pre-tRNAs with destabilizing mutations in the anticodon stem are targeted for Met22-dependent pre-tRNA decay (MPD). Multiple SUP4_οc anticodon stem variants that are subject to MPD each perturb the bulge-helix-bulge structure formed by the anticodon stem-loop and intron, which is important for splicing, resulting in substantial accumulation of end-matured unspliced pre-tRNA as well as pre-tRNA decay. Mutations that restore exon-intron structure commensurately reduce pre-tRNA accumulation and MPD. The MPD pathway can contribute substantially to decay of anticodon stem variants, since pre-tRNA decay is largely suppressed by removal of the intron or by restoration of exon-intron structure, each also resulting in increased tRNA levels. The MPD pathway is general as it extends to variants of tRNA^Tyr(GUA) and tRNA^Ser(CGA) These results demonstrate that the integrity of the anticodon stem-loop and the efficiency of tRNA splicing are monitored by a quality control pathway.

TA, GT and AC are significantly under-represented in open reading frames of prokaryotic and eukaryotic protein-coding genes

Genomes can be considered a combination of 16 dinucleotides. Analysing the relative abundance of different dinucleotides may reveal important features of genome evolution. In present study, we conducted extensive surveys on the relative abundances of dinucleotides in various genomic components of 28 bacterial, 20 archaean, 19 fungal, 24 plant and 29 animal species. We found that TA, GT and AC are significantly under-represented in open reading frames of all organisms and in intergenic regions and introns of most organisms. Specific dinucleotides are of greatly varied usage at different codon positions. The significantly low representations of TA, GT and AC are considered the evolutionary consequences of preventing formation of pre-mature stop codons and of reducing intron-splicing options in candidate primary mRNA sequences. These data suggest that a reduction of TA and GT occurred on both strands of the DNA sequence at an early stage of de novo gene birth. Interestingly, GT and AC are also significantly under-represented in current prokaryotic genomes, suggesting that ancient prokaryotic protein-coding genes might have contained introns. The greatly varied usages of specific dinucleotides at different codon positions are considered evolutionary accommodations to compensate the unavailability of specific codons and to avoid formation of pre-mature stop codons. This is the first report presenting data of dinucleotide relative abundance to indicate the possible existence of spliceosomal introns in ancient prokaryotic genes and to hypothesize early steps of de novo gene birth.

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.

Genome-Wide Analysis Reveals Ancestral Lack of Seventeen Different tRNAs and Clade-Specific Loss of tRNA-CNNs in Archaea

Transfer RNA (tRNA) is a category of RNAs that specifically decode messenger RNAs (mRNAs) into proteins by recognizing a set of 61 codons commonly adopted by different life domains. The composition and abundance of tRNAs play critical roles in shaping codon usage and pairing bias, which subsequently modulate mRNA translation efficiency and accuracy. Over the past few decades, effort has been concentrated on evaluating the specificity and redundancy of different tRNA families. However, the mechanism and processes underlying tRNA evolution have only rarely been investigated. In this study, by surveying tRNA genes in 167 completely sequenced genomes, we systematically investigated the composition and evolution of tRNAs in Archaea from a phylogenetic perspective. Our data revealed that archaeal genomes are compact in both tRNA types and copy number. Generally, no more than 44 different types of tRNA are present in archaeal genomes to decode the 61 canonical codons, and most of them have only one gene copy per genome. Among them, tRNA-Met was significantly overrepresented, with an average of three copies per genome. In contrast, the tRNA-UAU and 16 tRNAs with A-starting anticodons (tRNA-ANNs) were rarely detected in all archaeal genomes. The conspicuous absence of these tRNAs across the archaeal phylogeny suggests they might have not been evolved in the common ancestor of Archaea, rather than have lost independently from different clades. Furthermore, widespread absence of tRNA-CNNs in the Methanococcales and Methanobacteriales genomes indicates convergent loss of these tRNAs in the two clades. This clade-specific tRNA loss may be attributing to the reductive evolution of their genomes. Our data suggest that the current tRNA profiles in Archaea are contributed not only by the ancestral tRNA composition, but also by differential maintenance and loss of redundant tRNAs.

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.

Current Bacterial Gene Encoding Capsule Biosynthesis Protein CapI Contains Nucleotides Derived from Exonization

Since the proposition of introns-early hypothesis, although many studies have shown that most eukaryotic ancestors possessed intron-rich genomes, evidence of intron existence in genomes of ancestral bacteria has still been absent. While not a single intron has been found in all protein-coding genes of current bacteria, analyses on bacterial genes horizontally transferred into eukaryotes at ancient time may provide evidence of intron existence in bacterial ancestors. In this study, a bacterial gene encoding capsule biosynthesis protein CapI was found in the genome of sea anemone, Nematostella vectensis. This horizontally transferred gene contains a phase 1 intron of 40 base pairs. The nucleotides of this intron have high sequence identity with those encoding amino acids in current bacterial CapI gene, indicating that the intron and the amino acid-coding nucleotides are originated from the same ancestor sequence. Moreover, 5′-splice site of this intron is located in a GT-poor region associated with a closely following AG-rich region, suggesting that deletion mutation at 5′-splice site has been employed to remove this intron and the intron-like amino acid-coding nucleotides in current bacterial CapI gene are derived from exonization. These data suggest that bacterial CapI gene contained intron(s) at ancient time. This is the first report providing the result of sequence analysis to suggest possible existence of spliceosomal introns in ancestral bacterial genes. The methodology employed in this study may be used to identify more such evidence that would aid in settlement of the dispute between introns-early and introns-late theories.

Coevolution Theory of the Genetic Code at Age Forty: Pathway to Translation and Synthetic Life

The origins of the components of genetic coding are examined in the present study. Genetic information arose from replicator induction by metabolite in accordance with the metabolic expansion law. Messenger RNA and transfer RNA stemmed from a template for binding the aminoacyl-RNA synthetase ribozymes employed to synthesize peptide prosthetic groups on RNAs in the Peptidated RNA World. Coevolution of the genetic code with amino acid biosynthesis generated tRNA paralogs that identify a last universal common ancestor (LUCA) of extant life close to Methanopyrus, which in turn points to archaeal tRNA introns as the most primitive introns and the anticodon usage of Methanopyrus as an ancient mode of wobble. The prediction of the coevolution theory of the genetic code that the code should be a mutable code has led to the isolation of optional and mandatory synthetic life forms with altered protein alphabets.

The distribution, diversity, and importance of 16S rRNA gene introns in the order Thermoproteales

Background

Intron sequences are common in 16S rRNA genes of specific thermophilic lineages of Archaea, specifically the Thermoproteales (phylum Crenarchaeota). Environmental sequencing (16S rRNA gene and metagenome) from geothermal habitats in Yellowstone National Park (YNP) has expanded the available datasets for investigating 16S rRNA gene introns. The objectives of this study were to characterize and curate archaeal 16S rRNA gene introns from high-temperature habitats, evaluate the conservation and distribution of archaeal 16S rRNA introns in geothermal systems, and determine which “universal” archaeal 16S rRNA gene primers are impacted by the presence of intron sequences.

Results

Several new introns were identified and their insertion loci were constrained to thirteen locations across the 16S rRNA gene. Many of these introns encode homing endonucleases, although some introns were short or partial sequences. Pyrobaculum, Thermoproteus, and Caldivirga 16S rRNA genes contained the most abundant and diverse intron sequences. Phylogenetic analysis of introns revealed that sequences within the same locus are distributed biogeographically. The most diverse set of introns were observed in a high-temperature, circumneutral (pH 6) sulfur sediment environment, which also contained the greatest diversity of different Thermoproteales phylotypes.

Conclusions

The widespread presence of introns in the Thermoproteales indicates a high probability of misalignments using different “universal” 16S rRNA primers employed in environmental microbial community analysis.

Reviewers

This article was reviewed by Dr. Eugene Koonin and Dr. W. Ford Doolittle.

Electronic supplementary material

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

Cutting, dicing, healing and sealing: the molecular surgery of tRNA

All organisms encode transfer RNAs (tRNAs) that are synthesized as precursor molecules bearing extra sequences at their 5' and 3' ends; some tRNAs also contain introns, which are removed by splicing. Despite commonality in what the ultimate goal is (i.e., producing a mature tRNA), mechanistically, tRNA splicing differs between Bacteria and Archaea or Eukarya. The number and position of tRNA introns varies between organisms and even between different tRNAs within the same organism, suggesting a degree of plasticity in both the evolution and persistence of modern tRNA splicing systems. Here we will review recent findings that not only highlight nuances in splicing pathways but also provide potential reasons for the maintenance of introns in tRNA. Recently, connections between defects in the components of the tRNA splicing machinery and medically relevant phenotypes in humans have been reported. These differences will be discussed in terms of the importance of splicing for tRNA function and in a broader context on how tRNA splicing defects can often have unpredictable consequences.

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.

Eukaryotic tRNAs fingerprint invertebrates vis-à-vis vertebrates

During translation, aminoacyl-tRNA synthetases recognize the identities of the tRNAs to charge them with their respective amino acids. The conserved identities of 58,244 eukaryotic tRNAs of 24 invertebrates and 45 vertebrates in genomic tRNA database were analyzed and their novel features extracted. The internal promoter sequences, namely, A-Box and B-Box, were investigated and evidence gathered that the intervention of optional nucleotides at 17a and 17b correlated with the optimal length of the A-Box. The presence of canonical transcription terminator sequences at the immediate vicinity of tRNA genes was ventured. Even though non-canonical introns had been reported in red alga, green alga, and nucleomorph so far, fairly motivating evidence of their existence emerged in tRNA genes of other eukaryotes. Non-canonical introns were seen to interfere with the internal promoters in two cases, questioning their transcription fidelity. In a first of its kind, phylogenetic constructs based on tRNA molecules delineated and built the trees of the vast and diverse invertebrates and vertebrates. Finally, two tRNA models representing the invertebrates and the vertebrates were drawn, by isolating the dominant consensus in the positional fluctuations of nucleotide compositions.

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.

Origin of spliceosomal introns and alternative splicing

In this work we review the current knowledge on the prehistory, origins, and evolution of spliceosomal introns. First, we briefly outline the major features of the different types of introns, with particular emphasis on the nonspliceosomal self-splicing group II introns, which are widely thought to be the ancestors of spliceosomal introns. Next, we discuss the main scenarios proposed for the origin and proliferation of spliceosomal introns, an event intimately linked to eukaryogenesis. We then summarize the evidence that suggests that the last eukaryotic common ancestor (LECA) had remarkably high intron densities and many associated characteristics resembling modern intron-rich genomes. From this intron-rich LECA, the different eukaryotic lineages have taken very distinct evolutionary paths leading to profoundly diverged modern genome structures. Finally, we discuss the origins of alternative splicing and the qualitative differences in alternative splicing forms and functions across lineages.

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.

RNA-Seq analyses reveal the order of tRNA processing events and the maturation of C/D box and CRISPR RNAs in the hyperthermophile Methanopyrus kandleri

The methanogenic archaeon Methanopyrus kandleri grows near the upper temperature limit for life. Genome analyses revealed strategies to adapt to these harsh conditions and elucidated a unique transfer RNA (tRNA) C-to-U editing mechanism at base 8 for 30 different tRNA species. Here, RNA-Seq deep sequencing methodology was combined with computational analyses to characterize the small RNome of this hyperthermophilic organism and to obtain insights into the RNA metabolism at extreme temperatures. A large number of 132 small RNAs were identified that guide RNA modifications, which are expected to stabilize structured RNA molecules. The C/D box guide RNAs were shown to exist as circular RNA molecules. In addition, clustered regularly interspaced short palindromic repeats RNA processing and potential regulatory RNAs were identified. Finally, the identification of tRNA precursors before and after the unique C8-to-U8 editing activity enabled the determination of the order of tRNA processing events with termini truncation preceding intron removal. This order of tRNA maturation follows the compartmentalized tRNA processing order found in Eukaryotes and suggests its conservation during evolution.

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.

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.

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.

Fidelity of tRNA 5'-maturation: a possible basis for the functional dependence of archaeal and eukaryal RNase P on multiple protein cofactors

RNase P, which catalyzes tRNA 5′-maturation, typically comprises a catalytic RNase P RNA (RPR) and a varying number of RNase P proteins (RPPs): 1 in bacteria, at least 4 in archaea and 9 in eukarya. The four archaeal RPPs have eukaryotic homologs and function as heterodimers (POP5•RPP30 and RPP21•RPP29). By studying the archaeal Methanocaldococcus jannaschii RPR's cis cleavage of precursor tRNA^Gln (pre-tRNA^Gln), which lacks certain consensus structures/sequences needed for substrate recognition, we demonstrate that RPP21•RPP29 and POP5•RPP30 can rescue the RPR's mis-cleavage tendency independently by 4-fold and together by 25-fold, suggesting that they operate by distinct mechanisms. This synergistic and preferential shift toward correct cleavage results from the ability of archaeal RPPs to selectively increase the RPR's apparent rate of correct cleavage by 11 140-fold, compared to only 480-fold for mis-cleavage. Moreover, POP5•RPP30, like the bacterial RPP, helps normalize the RPR's rates of cleavage of non-consensus and consensus pre-tRNAs. We also show that archaeal and eukaryal RNase P, compared to their bacterial relatives, exhibit higher fidelity of 5′-maturation of pre-tRNA^Gln and some of its mutant derivatives. Our results suggest that protein-rich RNase P variants might have evolved to support flexibility in substrate recognition while catalyzing efficient, high-fidelity 5′-processing.

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.

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.

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.

Cleavage of intron from the standard or non-standard position of the precursor tRNA by the splicing endonuclease of Aeropyrum pernix, a hyper-thermophilic Crenarchaeon, involves a novel RNA recognition site in the Crenarchaea specific loop

In Crenarchaea, several tRNA genes are predicted to express precursor-tRNAs (pre-tRNAs) with canonical or non-canonical introns at various positions. We initially focused on the tRNA(Thr) species of hyperthermophilic crenarchaeon, Aeropyrum pernix (APE) and found that in the living APE cells three tRNA(Thr) species were transcribed and subsequently matured to functional tRNAs. During maturation, introns in two of them were cleaved from standard and non-standard positions. Biochemical studies revealed that the APE splicing endonuclease (APE-EndA) removed both types of introns, including the non-canonical introns, without any nucleotide modification. To clarify the underlying reasons for broad substrate specificity of APE-EndA, we determined the crystal structure of wild-type APE-EndA and subsequently compared its structure with that of Archaeaoglobus fulgidus (AFU)-EndA, which has narrow substrate specificity. Remarkably, structural comparison revealed that APE-EndA possesses a Crenarchaea specific loop (CSL). Introduction of CSL into AFU-EndA enhanced its intron-cleaving activity irrespective of the position or motif of the intron. Thus, our biochemical and crystallographic analyses of the chimera-EndA demonstrated that the CSL is responsible for the broad substrate specificity of APE-EndA. Furthermore, mutagenesis studies revealed that Lys44 in CSL functions as the RNA recognition site.

Discovery of permuted and recently split transfer RNAs in Archaea

Background

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

Results

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

Conclusions

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