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.     

Avoidance of antisense, antiterminator tRNA anticodons in vertebrate mitochondria

Protein synthesis (translation) stops at stop codons, codons not complemented by tRNA anticodons. tRNAs matching stops, antitermination (Ter) tRNAs, prevent translational termination, producing dysfunctional proteins. Genomes avoid tRNAs with anticodons whose complement (the anticodon of the ‘antisense’ tRNA) matches stops. This suggests that antisense tRNAs, which also form cloverleaves, are occasionally expressed. Mitochondrial antisense tRNA expression is plausible, because both DNA strands are transcribed as single RNAs, and tRNA structures signal RNA maturation. Results describe potential antisense Ter tRNAs in mammalian mitochondrial genomes detected by tRNAscan-SE, and evidence for adaptations preventing translational antitermination: genomes possessing Ter tRNAs use less corresponding stop codons; antisense Ter tRNAs form weaker cloverleaves than homologuous non-Ter antisense tRNAs; and genomic stop codon usages decrease with stabilities of codon-anticodon interactions and of Ter tRNA cloverleaves. This suggests that antisense tRNAs frequently function in translation. Results suggest that opposite strand coding is exceptional in modern genes, yet might be frequent for mitochondrial tRNAs. This adds antisense tRNA templating to other mitochondrial tRNA functions: sense tRNA templating, formation and regulation of secondary (light strand DNA) replication origins. Antitermination probably affects mitochondrial degenerative diseases and ageing: pathogenic mutations are twice as frequent in tRNAs with antisense Ter anticodons than in other tRNAs, and species lacking mitochondrial antisense Ter tRNAs have longer mean maximal lifespans than those possessing antisense Ter tRNAs.     

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

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

In vitro hyperprocessing of Drosophila tRNAs by the catalytic RNA of RNase P the cloverleaf structure of tRNA is not always stable?

We have previously reported that the catalytic RNA subunit of RNase P of Escherichia coli (M1 RNA) cleaves Drosophila initiator methionine tRNA (tRNA(Met)i) within the mature tRNA sequence to produce specific fragments. This cleavage was dependent on the occurrence of an altered conformation of the tRNA substrate. We call this further cleavage hyperprocessing. In the present paper, to search for another tRNA that can be hyperprocessed in vitro, we used total mature tRNAs from Drosophila as substrates for the in vitro M1 RNA reaction. We found that some tRNAs can be hyperprocessed by M1 RNA and that two such tRNAs are an alanine tRNA and a histidine tRNA. Using mutant substrates of these tRNAs, we also show that the hyperprocessing by M1 RNA is dependent on the occurrence of altered conformations of these tRNAs. The altered conformations were very similar to that of tRNA(Met)i. We show here that M1 RNA can be used as a powerful tool to detect the alternative conformation of tRNAs. The relationship between these hyperprocessing reactions and stability of the tRNA structure will also be discussed.     

Mature methyl-deficient tRNA isolated from a mammary adenocarcinoma

A transplantable rat tumor, mammary adenocarcinoma 13762, accumulates tRNA which can be methylated in vitro by mammalian tRNA (adenine-1) methyltransferase. This unusual ability of the tumor RNA to serve as substrate for a homologous tRNA methylating enzyme is correlated with unusually low levels of the A₅₈-specific adenine-1 methyltransferase. The nature of the methyl-accepting RNA has been examined by separating tumor tRNA on two-dimensional polyacrylamide gels. Comparisons of ethidium bromide-stained gels of tumor vs. liver tRNA show no significant quantitative differences and no accumulation of novel tRNAs or precursor tRNAs in adenocarcinoma RNA. Two-dimensional separations of tumor RNA after in vitro [¹⁴C]methylation using purified adenine-1 methyltransferase indicate that about 25% of the tRNA species are strongly methyl-accepting RNAs. Identification of six of the tRNAs separated on two-dimensional gels has been carried out by hybridization of cloned tRNA genes to Northern blots. Three of these, tRNA^Lys3, tRNA^Gln and tRNA^Met_i, are among the adenocarcinoma methyl-accepting RNAs. The other three RNAs, all of which are leucine-specific tRNAs, show no methyl-accepting properties. Our results suggest that low levels of a tRNA methyltransferase in the adenocarcinoma cause selected species of tRNA to escape the normal A₅₈ methylation, resulting in the appearance of several mature tRNAs which are deficient in 1-methyladenine. The methyl-accepting tRNAs from the tumor appear as ethidium bromide-stained spots of similar intensity to those seen for RNA from rat liver; therefore, methyladenine deficiency does not seem to impair processing of these tRNAs.     

Identification of nuclear encoded precursor tRNAs within the mitochondrion of Trypanosoma brucei.

RNAs that function in mitochondria are typically encoded by the mitochondrial DNA. However, the mitochondrial tRNAs of Trypanosoma brucei are encoded by the nuclear DNA and therefore must be imported into the mitochondrion. It is becoming evident that RNA import into mitochondria is phylogenetically widespread and is essential for cellular processes, but virtually nothing is known about the mechanism of RNA import. We have identified and characterized mitochondrial precursor tRNAs in T. brucei. The identification of mitochondrially located precursor tRNAs clearly indicates that mitochondrial tRNAs are imported as precursors. The mitochondrial precursor tRNAs hybridize to cloned nuclear tRNA genes, label with [alpha-32P]CTP using yeast tRNA nucleotidyltransferase and in isolated mitochondria via an endogenous nucleotidyltransferase-like activity, and are processed to mature tRNAs by Escherichia coli and yeast mitochondrial RNase P. We show that T. brucei mitochondrial extract contains an RNase P activity capable of processing a prokaryotic tRNA precursor as well as the T. brucei tRNA precursors. Precursors for tRNA(Asn) and tRNA(Leu) were detected on Northern blots of mitochondrial RNA, and the 5' ends of these RNAs were characterized by primer extension analysis. The structure of the precursor tRNAs and the significance of nuclear encoded precursor tRNAs within the mitochondrion are discussed.     

In silico screening of archaeal tRNA-encoding genes having multiple introns with bulge-helix-bulge splicing motifs

In archaeal species, several transfer RNA genes have been reported to contain endogenous introns. Although most of the introns are located at anticodon loop regions between nucleotide positions 37 and 38, a number of introns at noncanonical sites and six cases of tRNA genes containing two introns have also been documented. However, these tRNA genes are often missed by tRNAscan-SE, the software most widely used for the annotation of tRNA genes. We previously developed SPLITS, a computational tool to identify tRNA genes containing one intron at a noncanonical position on the basis of its discriminative splicing motif, but the software was limited in the detection of tRNA genes with multiple introns at noncanonical sites. In this study, we initially updated the system as SPLITSX in order to correctly predict known tRNA genes as well as novel ones with multiple introns. By a comprehensive search for tRNA genes in 29 archaeal genomes using SPLITSX, we listed 43 novel candidates that contain introns at noncanonical sites. As a result, 15 contained two introns and three contained three introns within the respective putative tRNA genes. Moreover, the candidates completely complemented all the codons of two archaeal species of uncultured methanogenic archaeon, RC-I and Thermofilum pendens Hrk 5, with novel candidates that were not detectable by tRNAscan-SE alone.     

Human tyrosine tRNA is also internally cleavable by E. coli ribonuclease P RNA ribozyme in vitro

Transfer RNA is an essential molecule for biological system, and each tRNA molecule commonly has a cloverleaf structure. Previously, we experimentally showed that some Drosophila tRNA (tRNA(Ala), tRNA(His), and tRNA(iMet)) molecules fit to form another, non-cloverleaf, structure in which the 3'-half of the tRNA molecules forms an alternative hairpin, and that the tRNA molecules are internally cleaved by the catalytic RNA of bacterial ribonuclease P (RNase P). Until now, the hyperprocessing reaction of tRNA has only been reported with Drosophila tRNAs. This time, we applied the hyperprocessing reaction to one of human tRNAs, human tyrosine tRNA, and we showed that this tRNA was also hyperprocessed by E. coli RNase P RNA. This tRNA is the first example for hyperprocessed non-Drosophila tRNAs. The results suggest that the hyperprocessing reaction can be a useful tool detect destablized tRNA molecules from any species.     

Pathogenic mutations in antisense mitochondrial tRNAs

Pathogenic mutations in mitochondrial tRNAs are 6.5 times more frequent than in other mitochondrial genes. This suggests that tRNA mutations perturb more than one function. A potential additional tRNA gene function is that of templating for antisense tRNAs. Pathogenic mutations weaken cloverleaf secondary structures of sense tRNAs. Analyses here show similar effects for most antisense tRNAs, especially after adjusting for associations between sense and antisense cloverleaf stabilities. These results imply translational activity by antisense tRNAs. For sense tRNAs Ala and Ser UCN, pathogenicity associates as much with sense as with antisense cloverleaf formation. For tRNA Pro, pathogenicity seems associated only with antisense, not sense tRNA cloverleaf formation. Translational activity by antisense tRNAs is expected for the 11 antisense tRNAs processed by regular sense RNA maturation, those recognized by their cognate amino acid’s tRNA synthetase, and those forming relatively stable cloverleaves as compared to their sense counterpart. Most antisense tRNAs probably function routinely in translation and extend the tRNA pool (extension hypothesis); others do not (avoidance hypothesis). The greater the expected translational activity of an antisense tRNA, the more pathogenic mutations weaken its cloverleaf secondary structure. Some evidence for RNA interference, a more classical role for antisense tRNAs, exists only for tRNA Ser UCN. Mutation pathogenicity probably frequently results from a mixture of effects due to sense and antisense tRNA translational activity for many mitochondrial tRNAs. Genomic studies should routinely explore for translational activity by antisense tRNAs.     

RNA sequence analysis using covariance models.

We describe a general approach to several RNA sequence analysis problems using probabilistic models that flexibly describe the secondary structure and primary sequence consensus of an RNA sequence family. We call these models 'covariance models'. A covariance model of tRNA sequences is an extremely sensitive and discriminative tool for searching for additional tRNAs and tRNA-related sequences in sequence databases. A model can be built automatically from an existing sequence alignment. We also describe an algorithm for learning a model and hence a consensus secondary structure from initially unaligned example sequences and no prior structural information. Models trained on unaligned tRNA examples correctly predict tRNA secondary structure and produce high-quality multiple alignments. The approach may be applied to any family of small RNA sequences.     

Undetected antisense tRNAs in mitochondrial genomes?

Background  

The hypothesis that both mitochondrial (mt) complementary DNA strands of tRNA genes code for tRNAs (sense-antisense coding) is explored. This could explain why mt tRNA mutations are 6.5 times more frequently pathogenic than in other mt sequences. Antisense tRNA expression is plausible because tRNA punctuation signals mt sense RNA maturation: both sense and antisense tRNAs form secondary structures potentially signalling processing. Sense RNA maturation processes by default 11 antisense tRNAs neighbouring sense genes. If antisense tRNAs are expressed, processed antisense tRNAs should have adapted more for translational activity than unprocessed ones. Four tRNA properties are examined: antisense tRNA 5′ and 3′ end processing by sense RNA maturation and its accuracy, cloverleaf stability and misacylation potential.     

Conservation and Losses of Non-Coding RNAs in Avian Genomes

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

Mitochondrial import of only one of three nuclear-encoded glutamine tRNAs in Tetrahymena thermophila.

The mitochondrial genome of Tetrahymena does not appear to encode enough tRNAs to perform mitochondrial protein synthesis. It has therefore been proposed that nuclear-encoded tRNAs are imported into the mitochondria. T.thermophila has three major glutamine tRNAs: tRNA(Gln)(UUG), tRNA(Gln)(UUA) and tRNA(Gln)(CUA). Each of these tRNAs functions in cytosolic translation. However, due to differences between the Tetrahymena nuclear and mitochondrial genetic codes, only tRNA(Gln)(UUG) has the capacity to function in mitochondrial translation as well. Here we show that approximately 10-20% of the cellular complement of tRNA(Gln)(UUG) is present in mitochondrial RNA fractions, compared with 1% or less for the other two glutamine tRNAs. Furthermore, this glutamine tRNA is encoded only by a family of nuclear genes, the sequences of several of which are presented. Finally, when marked versions of tRNA(Gln)(UUG) and tRNA(Gln)(UUA) flanked by identical sequences are expressed in the macronucleus, only the former undergoes mitochondrial import; thus sequences within tRNA(Gln)(UUG) direct import. Because tRNA(Gln)(UUG) is a constituent of mitochondrial RNA fractions and is encoded only by nuclear genes, and because ectopically expressed tRNA(Gln)(UUG) fractionates with mitochondria like its endogenous counterpart, we conclude that it is an imported tRNA in T.thermophila.     

Human Tyrosine tRNA Is Also Internally Cleavable by E. coli Ribonuclease P RNA Ribozyme in Vitro

Transfer RNA is an essential molecule for biological system, and each tRNA molecule commonly has a cloverleaf structure. Previously, we experimentally showed that some Drosophila tRNA (tRNA^Ala, tRNA^His, and tRNA_i ^Met) molecules fit to form another, non-cloverleaf, structure in which the 3'-half of the tRNA molecules forms an alternative hairpin, and that the tRNA molecules are internally cleaved by the catalytic RNA of bacterial ribonuclease P (RNase P). Until now, the hyperprocessing reaction of tRNA has only been reported with Drosophila tRNAs. This time, we applied the hyperprocessing reaction to one of human tRNAs, human tyrosine tRNA, and we showed that this tRNA was also hyperprocessed by E. coli RNase P RNA. This tRNA is the first example for hyperprocessed non-Drosophila tRNAs. The results suggest that the hyperprocessing reaction can be a useful tool to detect destablized tRNA molecules from any species.     

tRNA accumulation and suppression of the bldA phenotype during development in Streptomyces coelicolor

Streptomyces coelicolor undergoes distinct morphological changes as it grows on solid media where spores differentiate into vegetative and aerial mycelium that is followed by the production of spores. Deletion of bldA, encoding the rare tRNA(Leu) UAA, blocks development at the stage of vegetative mycelium formation. From previous data it appears that tRNA(Leu) UAA accumulates relatively late during growth while two other tRNAs do not. Here, we studied the expression of 17 different tRNAs including bldA tRNA, and the RNA subunit of the tRNA processing endoribonuclease RNase P. Our results showed that all selected tRNAs and RNase P RNA increased with time during development. However, accumulation of bldA tRNA and another rare tRNA(Leu) isoacceptor started at an earlier stage compared with the other tRNAs. We also introduced the bldA tRNA anticodon (UAA) into other tRNAs and introduced these into a bldA deletion strain. In particular, one such mutant tRNA derived from the tRNA(Leu) CAA isoacceptor suppressed the bldA phenotype. Thus, the bldA tRNA scaffold is not critical for function as a regulator of S. coelicolor cell differentiation. Further substitution experiments, in which the 5'- and 3'-flanking regions of the suppressor tRNA were changed, indicated that these regions were important for the suppression.     

Differential association of transfer RNAs with the genomes of murine, feline and primate retroviruses

The tRNAs that are bound to the genomic RNAs of several murine, feline, and primate retroviruses have been identified. Transfer RNAs were divided into those loosely bound and those tightly bound by stepwise thermal dissociation of the 70 S RNA. They were then identified and semiquantitated by aminoacylation. Proline tRNA is the most tenaciously bound tRNA in several strains of murine leukemia virus, two strains of feline leukemia virus, and the primate viruses simian sarcoma, baboon endogenous, and gibbon ape lymphoma. In the feline xenotropic virus, RD-114, tRNAGly is enriched in the most tightly bound fraction. In Mason-Pfizer monkey virus, as in the murine mammary tumor virus, tRNALys is the tRNA most tenaciously bound to its genomic RNA. Besides the most tightly associated tRNA, one or more different tRNAs are found in relatively large amounts in association with the 70 S RNA. (For convenience, we refer to the largest RNA ccomplex (50-70 S) isolated from any of the retroviruses studies as '70 S' RNA.) These tRNAs can be distinguished from the most tightly bound tRNA by the fact that they can be dissociated at lower temperatures. However, they occur in the same relative abundance as the tightly bound tRNA.     

Another cut for lysine tRNA: application of the hyperprocessing reaction reveals another stabilization strategy in metazoan lysine tRNAs

Recently, we revealed that the cloverleaf structure of some eukaryotic tRNAs is not always stable in vitro, and the denatured structures of these tRNAs are sometimes detected in bacterial RNase P reactions. We have designated the unusual internal cleavage reaction of these tRNAs as hyperprocessing. We have developed this hyperprocessing strategy as a useful tool for examining the stability of the tRNA cloverleaf structure. There are some common features in such unstable, hyperprocessible tRNAs, and the criteria for the hyperprocessing reaction of tRNA are extracted. Metazoan initiator methionine tRNAs and lysine tRNAs commonly fit the criteria, and are predicted to be hyperprocessible. The RNase P reactions of two metazoan lysine tRNAs from Homo sapiens and Caenorhabditis elegans, which fit the criteria, resulted in resistance to the internal cleavage reaction, while one bacterial lysine tRNA from Acholeplasma laidlawii, which also fits the criteria, was internally cleaved by the RNase P. The results showed that the metazoan lysine tRNAs examined are very stable without base modifications even under in vitro conditions. We also examined the 3'-half short construct of the human lysine tRNA, and the results showed that this RNA was internally cleaved by the enzyme. The results indicated that the human lysine tRNA has the ability to be hyperprocessed but is structurally stabilized in spite of lacking base modifications. A comparative study suggested, moreover, that the acceptor-stem bases should take part in the stabilization of metazoan lysine tRNAs. Our data strongly suggest that the cloverleaf shape of other metazoan lysine tRNAs should also be stabilized by means of similar strategies to in the case of human tRNA(Lys3).