Composition-based statistics and translated nucleotide searches: Improving the TBLASTN module of BLAST

BACKGROUND: TBLASTN is a mode of operation for BLAST that aligns protein sequences to a nucleotide database translated in all six frames. We present the first description of the modern implementation of TBLASTN, focusing on new techniques that were used to implement composition-based statistics for translated nucleotide searches. Composition-based statistics use the composition of the sequences being aligned to generate more accurate E-values, which allows for a more accurate distinction between true and false matches. Until recently, composition-based statistics were available only for protein-protein searches. They are now available as a command line option for recent versions of TBLASTN and as an option for TBLASTN on the NCBI BLAST web server. RESULTS: We evaluate the statistical and retrieval accuracy of the E-values reported by a baseline version of TBLASTN and by two variants that use different types of composition-based statistics. To test the statistical accuracy of TBLASTN, we ran 1000 searches using scrambled proteins from the mouse genome and a database of human chromosomes. To test retrieval accuracy, we modernize and adapt to translated searches a test set previously used to evaluate the retrieval accuracy of protein-protein searches. We show that composition-based statistics greatly improve the statistical accuracy of TBLASTN, at a small cost to the retrieval accuracy. CONCLUSION: TBLASTN is widely used, as it is common to wish to compare proteins to chromosomes or to libraries of mRNAs. Composition-based statistics improve the statistical accuracy, and therefore the reliability, of TBLASTN results. The algorithms used by TBLASTN are not widely known, and some of the most important are reported here. The data used to test TBLASTN are available for download and may be useful in other studies of translated search algorithms

Analysis of the mitochondrial and nuclear genomes of two basidiomycetes, Coprinus cinereus and Coprinus stercorarius

The mitochondrial and nuclear genomes of Coprinus stercorarius and C. cinereus were compared to assess their evolutionary relatedness and to characterize at the molecular level changes that have occurred since they diverged from a common ancestor. The mitochondrial genome of C. stercorarius (91.1 kb) is approximately twice as large as that of C. cinereus (43.3 kb). The pattern of restriction enzyme recognition sites shows both genomes to be circular, but reveals no clear homologies; furthermore, the order of structural genes is different in each species. The C. stercorarius mitochondrial genome contains a region homologous to a probe derived from the yeast mitochondrial var1 gene, whereas its nuclear genome does not. By contrast, the C. cinereus nuclear, but not mitochondrial, genome contains a region homologous to the var1 probe. Only a small fraction of either the nuclear or mitochondrial genomes, perhaps corresponding to the coding sequences, is capable of forming duplexes in interspecies solution reassociations, as measured by binding to hydroxylapatite. Those sequences capable of reassociating were found to have approximately 15% divergence for the mitochondrial genomes and 7%–15% divergence for the nuclear genomes, depending on the conditions of reassociation.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/46959/1/294_2004_Article_BF00447385.pd

The QDE-3 homologue RecQ-2 co-operates with QDE-3 in DNA repair in Neurospora crassa

The post-transcriptional gene silencing mechanism in Neurospora crassa, called quelling, was shown to involve the products of three genes termed quelling-defective. A homologue to the qde-3 gene encoding a putative RecQ-type DNA helicase was isolated and was named RecQ-2. Characterisation of the RecQ-2 gene has revealed that it is not involved in quelling, but may cooperate with the qde-3 gene product in a pathway that repairs damage to DNA caused by the chemical mutagens methyl methanesulfonate and N-methyl-N′-nitro-N-nitrosoguanidine. These results indicate that the qde-3 RecQ helicase may have a dual role in N. crassa, either acting alone as an essential component of the quelling mechanism or together with the RecQ-2 RecQ helicase, as part of a process to repair DNA lesions during replication

Heterogeneity of mitochondrial DNA from Saccharomyces cerevisiae and genetic information for tRNA.

Mitochondrial DNA from wild-type Saccharomyces cerevisiae and from an "extreme" petite mutant were analyzed by hybridization of several tRNAs on DNA fragments of different buoyant density, obtained by sonication and fractionation on a CsCl gradient. The hybridization patterns show that the genes for tRNAser, tRNAphe, tRNAhis, tRNAval, tRNAileu are present on wild-type mitochondrial DNA, while only genes for tRNAser and tRNAhis are present on petite mitochondrial DNA; moreover the hybridization patterns indicate that these genes are not clustered and suggest that more than one gene might exist for tRNAser and tRNAhis