Probabilistic phylogenetic inference with insertions and deletions.
A fundamental task in sequence analysis is to calculate the probability of a multiple alignment given a phylogenetic tree relating the sequences and an evolutionary model describing how sequences change over time. However, the most widely used phylogenetic models only account for residue substitution events. We describe a probabilistic model of a multiple sequence alignment that accounts for insertion and deletion events in addition to substitutions, given a phylogenetic tree, using a rate matrix augmented by the gap character. Starting from a continuous Markov process, we construct a non-reversible generative (birth-death) evolutionary model for insertions and deletions. The model assumes that insertion and deletion events occur one residue at a time. We apply this model to phylogenetic tree inference by extending the program dnaml in phylip. Using standard benchmarking methods on simulated data and a new "concordance test" benchmark on real ribosomal RNA alignments, we show that the extended program dnamlepsilon improves accuracy relative to the usual approach of ignoring gaps, while retaining the computational efficiency of the Felsenstein peeling algorithm.
Rivas
Elena
E
Janelia Farm Research Campus, Howard Hughes Medical Institute, Ashburn, Virginia, United States of America. rivase@janelia.hhmi.org
Eddy
Sean R
SR
D000465:Algorithms; D019295:Computational Biology; D003198:Computer Simulation; D004247:DNA / Q000235:genetics; D019143:Evolution, Molecular; D008390:Markov Chains; D008957:Models, Genetic; D015233:Models, Statistical; D016254:Mutagenesis, Insertional; D010802:Phylogeny; D016012:Poisson Distribution; D012335:RNA, Ribosomal / Q000235:genetics; D016415:Sequence Alignment / Q000706:statistics & numerical data*; D017384:Sequence Deletion; D012984:Software
D012335:RNA, Ribosomal; D004247:DNA
BMC Evol Biol. 2004 Sep 17;4:33
15377393
J Mol Evol. 1991 Aug;33(2):114-24
1920447
BMC Bioinformatics. 2005 Apr 01;6:83
15804354
Mol Biol Evol. 2004 Mar;21(3):419-27
14660689
Proteins. 2003 Aug 15;52(3):446-53
12866055
Mol Biol Evol. 2001 May;18(5):691-9
11319253
BMC Bioinformatics. 2006 Jun 08;7:292
16762073
J Comput Biol. 2004;11(2-3):413-28
15285899
Biometrics. 1995 Mar;51(1):100-13
7766767
Methods Enzymol. 1996;266:198-212
8743686
J Mol Evol. 1980 Dec;16(2):111-20
7463489
Bioinformatics. 2007 Feb 1;23(3):289-97
17110370
J Mol Evol. 1999 Jul;49(1):11-22
10368429
Mol Biol Evol. 1996 Jan;13(1):93-104
8583911
Mol Biol Evol. 2004 Mar;21(3):529-40
14694074
BMC Bioinformatics. 2004 Oct 26;5:166
15507142
J Mol Evol. 1985;22(2):160-74
3934395
Bioinformatics. 1999 Jun;15(6):446-54
10383470
Bioinformatics. 2003;19 Suppl 1:i147-57
12855451
Syst Biol. 2005 Aug;54(4):548-61
16085574
J Mol Evol. 1981;17(6):368-76
7288891
Bull Math Biol. 1989;51(1):79-94
2706403
BMC Bioinformatics. 2005 Mar 21;6:63
15780137
BMC Bioinformatics. 2005 Sep 27;6:236
16188037
J Comput Biol. 2003;10(6):869-89
14980015
Bioinformatics. 2004 Aug 12;20(12):1850-60
14988105
Comput Biol Chem. 2004 Jul;28(3):219-26
15261152
Genome Res. 2005 Aug;15(8):1034-50
16024819
Mol Biol Evol. 1996 Jan;13(1):105-14
8583885
Bioinformatics. 2001 Aug;17(8):754-5
11524383
Bioinformatics. 2001 Sep;17(9):803-20
11590097
J Mol Biol. 2004 Aug 6;341(2):617-31
15276848
Syst Biol. 2006 Oct;55(5):756-68
17060197
Bioinformatics. 2003 Jan 22;19(2):219-27
12538242
Biometrics. 1999 Mar;55(1):1-12
11318142
Bioinformatics. 2004 Aug 4;20 Suppl 1:i94-100
15262786
Mol Biol Evol. 1998 Dec;15(12):1600-11
9866196
J Mol Evol. 1996 Sep;43(3):304-11
8703097
Bioinformatics. 2003 Mar 1;19(4):490-9
12611804
Proteins. 2001 Oct 1;45(1):102-4
11536366
J Mol Evol. 1994 Jul;39(1):105-11
8064867
Bioinformatics. 2004 Sep 22;20(14):2175-80
15044240
J Mol Evol. 1992 Jan;34(1):3-16
1556741
J Mol Biol. 2000 Sep 8;302(1):265-79
10964574
Syst Biol. 2003 Oct;52(5):696-704
14530136
Proc Natl Acad Sci U S A. 1955 Jun 15;41(6):387-91
16589686
Genetics. 2000 May;155(1):431-49
10790415
Proc Natl Acad Sci U S A. 2004 Sep 28;101(39):13994-4001
15292512
Comput Appl Biosci. 1997 Oct;13(5):555-6
9367129
Genetics. 1995 Feb;139(2):993-1005
7713447
J Theor Biol. 2004 May 7;228(1):97-106
15064085
Genome Res. 2004 Mar;14(3):442-50
14993209
Proc Natl Acad Sci U S A. 2003 Dec 9;100(25):14960-5
14657378
Bioinformatics. 1998;14(2):157-63
9545448
Pac Symp Biocomput. 2001;:179-90
11262938
J Bioinform Comput Biol. 2006 Jun;4(3):721-44
16960972
Bioinformatics. 2005 Feb 15;21(4):456-63
15608047
Mol Biol Evol. 1994 May;11(3):459-68
8015439
J Mol Biol. 1996 Oct 25;263(2):196-208
8913301
Mol Biol Evol. 2001 Apr;18(4):481-90
11264399
Genome Res. 2004 Dec;14(12):2412-23
15574820
J Mol Biol. 2003 Oct 17;333(2):453-60
14529629
BMC Bioinformatics. 2002;3:2
11869452
Bioinformatics. 1998;14(8):726-33
9789098
J Mol Biol. 1986 Jul 20;190(2):159-65
3641921