A program for the identification of tRNA-like structures in DNA sequence data.

A computer algorithm has been developed which identifies tRNA genes and tRNA-like structures in DNA sequences. The program searches the sequence string for specific base positions that correspond to the invariant and semi-invariant bases found in tRNAs. The tRNA nature of the sequence is confirmed by the presence of complementary base pairing at the tRNA's calculated 5' and 3' ends (which in situ constitutes the amino-acyl stem region). The program achieves greater than 96% accuracy when run against known tRNA sequences in the Genbank database. The program is modular and is readily modified to allow searching either a file or database. The program is written in "C" and operates on a D.E.C. Vax 750. The utility of the algorithm is demonstrated by the identification of a distinctive tRNA structure in an intron of a published bovine hemoglobin gene.     

A microcomputer program for the identification of tRNA genes

A microcomputer program which locates tRNA genes within longDNA sequences is described. The search is performed either byidentifying tRNA-like secondary structures or by locating eukaryoticRNA polymerase III promoter consensus sequences. The programis also useful in finding inverted repeats allowing the formationof stem-loop secondary structures in tRNA. The program has beendeveloped in BASIC and 6502 Assembler and runs on the AppleII plus and He microcomputers. The execution is quite fast;all the operations are carried out in 1–90 s, dependingon the required task and on the sequence length. Received on March 1, 1985; accepted on April 25, 1985     

ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences

A computer program, ARAGORN, identifies tRNA and tmRNA genes. The program employs heuristic algorithms to predict tRNA secondary structure, based on homology with recognized tRNA consensus sequences and ability to form a base-paired cloverleaf. tmRNA genes are identified using a modified version of the BRUCE program. ARAGORN achieves a detection sensitivity of 99% from a set of 1290 eubacterial, eukaryotic and archaeal tRNA genes and detects all complete tmRNA sequences in the tmRNA database, improving on the performance of the BRUCE program. Recently discovered tmRNA genes in the chloroplasts of two species from the ‘green’ algae lineage are detected. The output of the program reports the proposed tRNA secondary structure and, for tmRNA genes, the secondary structure of the tRNA domain, the tmRNA gene sequence, the tag peptide and a list of organisms with matching tmRNA peptide tags.     

Mitochondrial genome sequence of <Emphasis Type="Italic">Unionicola foili</Emphasis> (Acari: Unionicolidae): a unique gene order with implications for phylogenetic inference

The mitochondrial genome of Unionicola foili is circular, 14,738 bp in length, and contains several notable features. The sequence and annotation revealed a unique gene order, continuing a pattern of highly rearranged mitochondrial genomes in the Trombidiformes. U. foili mitochondrial tRNA sequences predict non-canonical secondary structures for these molecules, and our annotation suggests an in-frame fusion between the nad4L and nad5 genes in this genome. The unique gene order and unusual tRNA structures could serve as idiosyncratic characters and have the potential to be phylogenetically informative.     

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.     

DNAMAT: an efficient graphic matrix sequence homology algorithm and its application to structural analysis

We present a fast algorithm to produce a graphic matrix representationof sequence homology. The algorithm is based on lexicographicalordering of fragments. It preserves most of the options of asimple naive algorithm with a significant increase in speed.This algorithm was the basis for a program, called DNAMAT, thathas been extensively tested during the last three years at theWeizmann Institute of Science and has proven to be very useful.In addition we suggest a way to extend our approach to analysea series of related DNA or RNA sequences, in order to determinecertain common structural features. The analysis is done by‘summing’ a set of dot-matrices to produce an overallmatrix that displays structural elements common to most of thesequences. We give an example of this procedure by analysingtRNA sequences. Received on June 26, 1986; accepted on September 28, 1986     

Nucleotide sequence of leucine transfer RNA 1 from Candida (Torulopsis) utilis.

Two species of 32P-labelled leucine tRNA were highly purified from Candida (Torulopsis) utilis by successive column chromatographies. The purified major species of leucine tRNA 1 was completely digested with ribonuclease T1 [EC 3.1.4.8] and with pancreatic ribonuclease A [EC 3.1.4.22]. The resulting fragments were fractionated, and their nucleotide sequences were determined according to Barrell (1). The results of analyses of the two ribonuclease digests were consistent with each other, and indicated that this tRNA is composed of 85 nucleotide residues, including 14 modified nucleotides. A tentative total sequence has been derived on the basis of several features in the cloverleaf structure for tRNA.     

Transfer RNA splicing in Saccharomyces cerevisiae: defining the substrates.

The primary sequences of all the tRNA precursors which contain intervening sequences and which accumulate in the Saccharomyces cerevisiae rnal mutant are presented. A combination of DNA and RNA sequence analysis has led to elucidation of the primary sequence of four hitherto uncharacterized precursors. The location of the intervening sequence has in all cases been unambiguously determined by analysis of the intermediates in the splicing reaction. Secondary structures based upon the tRNA cloverleaf are shown for all the tRNA precursors and discussed with respect to common recognition by the yeast splicing endonuclease.     

BIOSITE: a program for the interactive comparison of aligned homologous protein sequences

A program, BIOSITE, providing for the interactive visual comparisonof aligned homologous amino-acid sequences is presented, includingan example of its application. The program allows for two typesof comparison sequence to be generated: an ‘identity’sequence and a ‘difference’ sequence. These maybe used on subsets of sequences and in further comparisons toidentify candidate sites involved in a distinct functional property.The program should prove a useful tool for biologists engagedin understanding sequence—function relationships.     

Rodent type 2 Alu family,rat identifier sequence,rabbit C family,and bovine or goat 73-bp repeat may have evolved from tRNA genes

Summary Close structural resemblances between several mammalian highly or moderately repetitive families and some specific tRNAs were detected. The rodent type 2 Alu family, rat identifier (ID) sequences, rabbit C family, and bovine or goat 73-bp repeat are most homologous with lysine tRNA₅, phenylalanine tRNA, glycine tRNA, and glycine tRNA, respectively. The homologies extend to secondary structures, and the homologous nucleotides are located on nearly the same secondary structures. The repetitive families mentioned have a common structural organization, with a tRNA-like sequence devoid of an aminoacyl stem region. These features suggest that these repetitive families may be generated by nonhomologous recombination between a tRNA gene and a tRNA-unrelated block.     

Identifying potential tRNA genes in genomic DNA sequences.

We have developed an algorithm that automatically and reproducibly identifies potential tRNA genes in genomic DNA sequences, and we present a general strategy for testing the sensitivity of such algorithms. This algorithm is useful for the flagging and characterization of long genomic sequences that have not been experimentally analyzed for identification of functional regions, and for the scanning of nucleotide sequence databases for errors in the sequences and the functional assignments associated with them. In an exhaustive scan of the GenBank database, 97.5% of the 744 known tRNA genes were correctly identified (true-positives), and 42 previously unidentified sequences were predicted to be tRNAs. A detailed analysis of these latter predictions reveals that 16 of the 42 are very similar to known tRNA genes, and we predict that they do, in fact, code for tRNA, yielding a false-positive rate for the algorithm of 0.003%. The new algorithm and testing strategy are a considerable improvement over any previously described strategies for recognizing tRNA genes, and they allow detections of genes (including introns) embedded in long genomic sequences.     

The nucleotide sequence of a precursor to the glycine- and threonine-specific transfer ribonucleic acids of Escherichia coli.

The primary nucleotide sequence of an Escherichia coli tRNA precursor molecule has been determined. This precursor RNA, specified by the transducing phage lambdah80dglyTsuA36 thrT tyrT, accumulates in a mutant strain temperature-sensitive for RNase P activity. The 170-nucleotide precursor RNA is processed by E. coli extracts to form mature tRNA Gly 2 suA36 and tRNA Thr ACU/C. The sequence of the precursor is pG-U-U-C-C-A-G-G-A-U-G-C-G-G-G-C-A-U-C-G-U-A-U-A-A-U-G-G-C-U-A-U-U-A-C-C-U-C-A-G-C-C-U-N-C-U-A-A-G-C-U-G-A-U-G-A-U-G-C-G-G-G-T-psi-C-G-A-U-U-C-C-C-G-C-U-G-C-C-C-G-C-U-C-C-A-A-G-A-U-G-U-G-C-U-G-A-U-A-U-A-G-C-U-C-A-G-D-D-G-G-D-A-G-A-G-C-G-C-A-C-C-C-U-U-G-G-U-mt6A-A-G-G-G-U-G-A-G-m7G-U-C-G-G-C-A-G-T-psi-C-G-A-A-U-C-U-G-C-C-U-A-U-C-A-G-C-A-C-C-A-C-U-UOH(tRNA sequences are italicized). It contains the entire primary nucleotide sequences of tRNA Gly2 suA36 and tRNA Thr ACU/C, including the common 3'-terminal sequence, CCA. Nineteen additional nucleotides are present, with 10 at the 5' end, 3 at the 3' end, and the remaining 6 in the inter-tRNA spacer region. RNase P cleaves the precursor specifically at the 5' ends of the mature tRNA sequences.     

The nucleotide sequences of several tRNA genes from rat mitochondria: common features and relatedness to homologous species.

We have determined the nucleotide sequences of thirteen rat mt tRNA genes. The features of the primary and secondary structures of these tRNAs show that those for Gln, Ser, and f-Met resemble, while those for Lys, Cys, and Trp depart strikingly from the universal type. The remainder are slightly abnormal. Among many mammalian mt DNA sequences, those of mt tRNA genes are highly conserved, thus suggesting for those genes an additional, perhaps regulatory, function. A simple evolutionary relationship between the tRNAs of animal mitochondria and those of eukaryotic cytoplasm, of lower eukaryotic mitochondria or of prokaryotes, is not evident owing to the extreme divergence of the tRNA sequences in the two groups. However, a slightly higher homology does exist between a few animal mt tRNAs and those from prokaryotes or from lower eukaryotic mitochondria.     

A phenylalanine tRNA gene from Neurospora crassa: conservation of secondary structure involving an intervening sequence.

We have isolated and sequenced a tRNAPhe gene from Neurospora crassa. Hybridization analyses suggest that trnaPhe is the only tRNA encoded on the cloned 5 kb DNA fragment. The tRNAPhe gene contains an intervening sequence 16 nucleotides in length located one nucleotide 3' to the anticodon position. The tRNAPhe coding region of Neurospora and yeast are 91% conserved, whereas their intervening sequences are only 50% identical. The pattern of sequence conservation is consistent with a proposed secondary structure for the tRNA precursor in which the anticodon is base paired with the middle of the intervening sequence and the splice points are located in adjacent single-stranded loops. The DNA sequence following the tRNAPhe coding region is similar to sequences following other genes transcribed by RNA polymerase III in that it is AT-rich and includes a tract of A residues in the coding strand. In contrast, the sequence preceding the Neurospora tRNAPhe coding region does not resemble sequences preceding other sequenced tRNA genes.     

A FastA based compilation of higher plant mitochondrial tRNA genes.

A new version of the compilation of higher plant mitochondrial tRNA genes (http://www.ebi.ac.uk/service ) has been obtained by means of the FastA program for similarity searching in nucleotide sequence Databases. This approach improves the previous collection, which was based on literature data analysis. The current compilation contains 158 sequences with an increase of 43 units. In this paper, some interesting features of the new entries are briefly presented.