title: The whole alignment and nothing but the alignment: the problem of spurious alignment flanks.
abstract: Pairwise sequence alignment is a ubiquitous tool for inferring the evolution and function of DNA, RNA and protein sequences. It is therefore essential to identify alignments arising by chance alone, i.e. spurious alignments. On one hand, if an entire alignment is spurious, statistical techniques for identifying and eliminating it are well known. On the other hand, if only a part of the alignment is spurious, elimination is much more problematic. In practice, even the sizes and frequencies of spurious subalignments remain unknown. This article shows that some common scoring schemes tend to overextend alignments and generate spurious alignment flanks up to hundreds of base pairs/amino acids in length. In the UCSC genome database, e.g. spurious flanks probably comprise >18% of the human-fugu genome alignment. To evaluate the possibility that chance alone generated a particular flank on a particular pairwise alignment, we provide a simple 'overalignment' P-value. The overalignment P-value can identify spurious alignment flanks, thereby eliminating potentially misleading inferences about evolution and function. Moreover, by explicitly demonstrating the tradeoff between over- and under-alignment, our methods guide the rational choice of scoring schemes for various alignment tasks.
authors.affiliation: Computational Biology Research Center, Institute for Advanced Industrial Science and Technology, Tokyo 135-0064, Japan.
authors.forename: Martin C
authors.initials: MC
authors.lastname: Frith
authors.forename: Yonil
authors.initials: Y
authors.lastname: Park
authors.forename: Sergey L
authors.initials: SL
authors.lastname: Sheetlin
authors.forename: John L
authors.initials: JL
authors.lastname: Spouge
mesh_terms: D000818:Animals; D019295:Computational Biology; D003627:Data Interpretation, Statistical; D023281:Genomics; D006801:Humans; D011336:Probability; D016415:Sequence Alignment / Q000379:methods*
references.citation: Proc Natl Acad Sci U S A. 1992 Nov 15;89(22):10915-9
references.pmid: 1438297
references.citation: PLoS Genet. 2007 Feb 2;3(2):e22
references.pmid: 17274688
references.citation: J Comput Biol. 1999 Spring;6(1):91-112
references.pmid: 10223667
references.citation: Bioinformatics. 1999 Dec;15(12):1012-9
references.pmid: 10745991
references.citation: Nucleic Acids Res. 2007 Jan;35(Database issue):D61-5
references.pmid: 17130148
references.citation: Trends Genet. 2000 Jun;16(6):276-7
references.pmid: 10827456
references.citation: Proc Natl Acad Sci U S A. 1988 Apr;85(8):2444-8
references.pmid: 3162770
references.citation: Pac Symp Biocomput. 2002;:115-26
references.pmid: 11928468
references.citation: Proc Natl Acad Sci U S A. 2003 Sep 30;100(20):11484-9
references.pmid: 14500911
references.citation: Genome Res. 2003 Jan;13(1):103-7
references.pmid: 12529312
references.citation: PLoS Comput Biol. 2008 May 16;4(5):e1000077
references.pmid: 18464927
references.citation: PLoS Biol. 2004 Sep;2(9):E273
references.pmid: 15361937
references.citation: Ann Stat. 2009 Dec 1;37(6A):3697
references.pmid: 20148197
references.citation: Proc Natl Acad Sci U S A. 2008 Mar 4;105(9):3209-14
references.pmid: 18305160
references.citation: Genome Res. 2002 Jun;12(6):996-1006
references.pmid: 12045153
references.citation: Protein Eng. 1995 Oct;8(10):999-1009
references.pmid: 8771180
references.citation: J Comput Biol. 2002;9(2):243-60
references.pmid: 12015880
references.citation: Genome Res. 2008 Feb;18(2):298-309
references.pmid: 18073381
references.citation: J Mol Biol. 1982 Dec 15;162(3):705-8
references.pmid: 7166760
references.citation: Nucleic Acids Res. 1997 Sep 1;25(17):3389-402
references.pmid: 9254694
references.citation: J Mol Biol. 1981 Mar 25;147(1):195-7
references.pmid: 7265238
references.citation: Proc Natl Acad Sci U S A. 1991 Oct 15;88(20):8880-4
references.pmid: 1924347
references.citation: Protein Eng. 1996 Feb;9(2):127-32
references.pmid: 9005433