Nucleic Acids Research Advance Access published online on March 15, 2008
Nucleic Acids Research, doi:10.1093/nar/gkn105
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Computational Biology |
A space-efficient and accurate method for mapping and aligning cDNA sequences onto genomic sequence
1Department of Intelligence Science and Technology, Graduate School of Informatics, Kyoto University, Yoshida Honmachi, Sakyo-ku, Kyoto 606-8501 and 2National Institute of Advanced Industrial Science and Technology, Computational Biology Research Center, 2-42 Aomi, Koto-ku, Tokyo 135-0064, Japan
*To whom correspondence should be addressed. Tel: +81 75 753 9109; Fax: +81 75 753 9110; Email: o.gotoh{at}i.kyoto-u.ac.jp
Received November 22, 2007. Revised January 30, 2008. Accepted February 25, 2008.
The mapping and alignment of transcripts (cDNA, expressed sequence tag or amino acid sequences) onto a genomic sequence is a fundamental step for genome annotation, including gene finding and analyses of transcriptional activity, alternative splicing and nucleotide polymorphisms. As DNA sequence data of genomes and transcripts are accumulating at an unprecedented rate, steady improvement in accuracy, speed and space requirement in the computational tools for mapping/alignment is desired. We devised a multi-phase heuristic algorithm and implemented it in the development of the stand-alone computer program Spaln (space-efficient spliced alignment). Spaln is reasonably fast and space efficient; it requires <1 Gb of memory to map and align >120 000 Unigene sequences onto the unmasked whole human genome with a conventional computer, finishing the job in <6 h. With artificially introduced noise of various levels, Spaln significantly outperforms other leading alignment programs currently available with respect to the accuracy of mapped exon–intron structures. This performance is achieved without extensive learning procedures to adjust parameter values to a particular organism. According to the handiness and accuracy, Spaln may be used for studies on a wide area of genome analyses.