The whole alignment and nothing but the alignment: the problem of spurious alignment flanks.
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.
Frith
Martin C
MC
Computational Biology Research Center, Institute for Advanced Industrial Science and Technology, Tokyo 135-0064, Japan.
Park
Yonil
Y
Sheetlin
Sergey L
SL
Spouge
John L
JL
D000818:Animals; D019295:Computational Biology; D003627:Data Interpretation, Statistical; D023281:Genomics; D006801:Humans; D011336:Probability; D016415:Sequence Alignment / Q000379:methods*
Proc Natl Acad Sci U S A. 1992 Nov 15;89(22):10915-9
1438297
PLoS Genet. 2007 Feb 2;3(2):e22
17274688
J Comput Biol. 1999 Spring;6(1):91-112
10223667
Bioinformatics. 1999 Dec;15(12):1012-9
10745991
Nucleic Acids Res. 2007 Jan;35(Database issue):D61-5
17130148
Trends Genet. 2000 Jun;16(6):276-7
10827456
Proc Natl Acad Sci U S A. 1988 Apr;85(8):2444-8
3162770
Pac Symp Biocomput. 2002;:115-26
11928468
Proc Natl Acad Sci U S A. 2003 Sep 30;100(20):11484-9
14500911
Genome Res. 2003 Jan;13(1):103-7
12529312
PLoS Comput Biol. 2008 May 16;4(5):e1000077
18464927
PLoS Biol. 2004 Sep;2(9):E273
15361937
Ann Stat. 2009 Dec 1;37(6A):3697
20148197
Proc Natl Acad Sci U S A. 2008 Mar 4;105(9):3209-14
18305160
Genome Res. 2002 Jun;12(6):996-1006
12045153
Protein Eng. 1995 Oct;8(10):999-1009
8771180
J Comput Biol. 2002;9(2):243-60
12015880
Genome Res. 2008 Feb;18(2):298-309
18073381
J Mol Biol. 1982 Dec 15;162(3):705-8
7166760
Nucleic Acids Res. 1997 Sep 1;25(17):3389-402
9254694
J Mol Biol. 1981 Mar 25;147(1):195-7
7265238
Proc Natl Acad Sci U S A. 1991 Oct 15;88(20):8880-4
1924347
Protein Eng. 1996 Feb;9(2):127-32
9005433