L. Steven Johnson

BACKGROUND: Profile hidden Markov models (profile-HMMs) are sensitive tools for remote protein homology detection, but the main scoring algorithms, Viterbi or Forward, require considerable time to search large sequence databases. RESULTS: We have designed a series of database filtering steps, HMMERHEAD, that are applied prior to the scoring algorithms, as implemented in the HMMER package, in an effort to reduce search time. Using this heuristic, we obtain a 20-fold decrease in Forward and a 6-fold decrease in Viterbi search time with a minimal loss in sensitivity relative to the unfiltered approaches. We then implemented an iterative profile-HMM search method, JackHMMER, which employs the HMMERHEAD heuristic. Due to our search heuristic, we eliminated the subdatabase creation that is common in current iterative profile-HMM approaches. On our benchmark, JackHMMER detects 14% more remote protein homologs than SAM's iterative method T2K. CONCLUSIONS: Our search heuristic, HMMERHEAD, significantly reduces the time needed to score a profile-HMM against large sequence databases. This search heuristic allowed us to implement an iterative profile-HMM search method, JackHMMER, which detects significantly more remote protein homologs than SAM's T2K and NCBI's PSI-BLAST.

We applied deep sequencing technology to small RNA fractions from cells lytically infected with murine gammaherpesvirus 68 (gammaHV68) in order to define in detail small RNAs generated from a cluster of tRNA-related polycistronic structures located at the left end of the viral genome. We detected 10 new candidate microRNAs (miRNAs), six of which were confirmed by Northern blot analysis, leaving four as provisional. In addition, we determined that previously identified and annotated viral miRNA molecules were not necessarily represented as the most abundant sequence originating from a transcript. Based on these new small RNAs and previously reported gammaHV68 miRNAs, we were able to further describe and annotate the distinctive gammaHV68 tRNA-miRNA structures. We used this deep sequencing data and computational analysis to identify similar structures in the mouse genome and validated that these host structures also give rise to small RNAs. This reveals a possible convergent usage of tRNA/polymerase III (pol III) transcripts to generate small RNAs from both mammalian and viral genomes.

Box C/D RNAs are small, noncoding RNAs that function in RNA modification in eukaryotes and archaea. Here, we report that box C/D RNAs exist in the rare biological form of RNA circles in the hyperthermophilic archaeon Pyrococcus furiosus. Northern analysis of box C/D RNAs reveals two prominent RNA species of different electrophoretic mobilities in total P. furiosus RNA preparations. Together, the results of Northern, ribozyme, RT-PCR, and lariat debranching analyses indicate that the two species are circular and linear RNAs of similar length and abundance. It seems that most, if not all, species of box C/D RNAs exist as circles in P. furiosus. In addition, the circular RNAs are found in complexes with proteins required for box C/D RNA function. Our finding places box C/D RNAs among the extremely few circular RNAs known to exist in nature. Moreover, the unexpected discovery of circular box C/D RNAs points to the existence of a previously unrecognized biogenesis pathway for box C/D RNAs in archaea.