CodonTest: modeling amino acid substitution preferences in coding sequences.
Codon models of evolution have facilitated the interpretation of selective forces operating on genomes. These models, however, assume a single rate of non-synonymous substitution irrespective of the nature of amino acids being exchanged. Recent developments have shown that models which allow for amino acid pairs to have independent rates of substitution offer improved fit over single rate models. However, these approaches have been limited by the necessity for large alignments in their estimation. An alternative approach is to assume that substitution rates between amino acid pairs can be subdivided into rate classes, dependent on the information content of the alignment. However, given the combinatorially large number of such models, an efficient model search strategy is needed. Here we develop a Genetic Algorithm (GA) method for the estimation of such models. A GA is used to assign amino acid substitution pairs to a series of rate classes, where is estimated from the alignment. Other parameters of the phylogenetic Markov model, including substitution rates, character frequencies and branch lengths are estimated using standard maximum likelihood optimization procedures. We apply the GA to empirical alignments and show improved model fit over existing models of codon evolution. Our results suggest that current models are poor approximations of protein evolution and thus gene and organism specific multi-rate models that incorporate amino acid substitution biases are preferred. We further anticipate that the clustering of amino acid substitution rates into classes will be biologically informative, such that genes with similar functions exhibit similar clustering, and hence this clustering will be useful for the evolutionary fingerprinting of genes.
Delport
Wayne
W
Department of Pathology, University of California, San Diego, La Jolla, California, United States of America.
Scheffler
Konrad
K
Botha
Gordon
G
Gravenor
Mike B
MB
Muse
Spencer V
SV
Kosakovsky Pond
Sergei L
SL
D000465:Algorithms; D019943:Amino Acid Substitution / Q000235:genetics*; D003062:Codon; D003198:Computer Simulation; D004259:DNA-Directed DNA Polymerase / Q000235:genetics; D019143:Evolution, Molecular; D015497:HIV-1 / Q000235:genetics; D006388:Hemagglutinins / Q000235:genetics; D006801:Humans; D008390:Markov Chains; D008957:Models, Genetic; D016415:Sequence Alignment
D003062:Codon; D006388:Hemagglutinins; D004259:DNA-Directed DNA Polymerase
J Mol Evol. 1985;22(2):160-74
3934395
Mol Biol Evol. 2009 May;26(5):1155-61
19233963
Bioinformatics. 1998;14(9):817-8
9918953
Bioinformatics. 2005 Mar 1;21(5):676-9
15509596
Mol Biol Evol. 1994 Sep;11(5):715-24
7968485
Comput Appl Biosci. 1992 Jun;8(3):275-82
1633570
Mol Biol Evol. 1994 Sep;11(5):725-36
7968486
Bioinformatics. 2005 May 15;21(10):2531-3
15713735
J Virol. 2008 Jan;82(2):596-601
17942553
J Mol Evol. 1996 Apr;42(4):459-68
8642615
Genetics. 2004 Aug;167(4):2027-43
15342538
J Mol Evol. 1981;17(6):368-76
7288891
PLoS Pathog. 2007 Jul;3(7):e94
17616974
Nucleic Acids Res. 2006 Jan 1;34(Database issue):D327-31
16381879
Mol Biol Evol. 2009 Feb;26(2):255-71
18922761
Mol Biol Evol. 2007 Jan;24(1):159-70
17038448
Science. 2008 Apr 18;320(5874):340-6
18420927
PLoS One. 2010 Jul 21;5(7):e11587
20657773
Mol Biol Evol. 2005 Dec;22(12):2375-85
16107593
Mol Phylogenet Evol. 2007 Feb;42(2):298-307
16942891
Mol Biol Evol. 2008 Jul;25(7):1253-6
18397919
BMC Bioinformatics. 2006 Mar 16;7:148
16542458
Mol Biol Evol. 2009 Jul;26(7):1663-76
19383983
PLoS Comput Biol. 2007 Nov;3(11):e231
18039027
Proteins. 1997 Jul;28(3):405-20
9223186
Mol Biol Evol. 2008 Jul;25(7):1307-20
18367465
J Mol Evol. 2002 Jul;55(1):65-73
12165843
J Mol Evol. 2000 Apr;50(4):348-58
10795826
Mol Biol Evol. 1987 Jul;4(4):406-25
3447015
Mol Biol Evol. 2005 May;22(5):1208-22
15703242
Mol Biol Evol. 2007 Jul;24(7):1464-79
17400572
Mol Biol Evol. 2000 Jun;17(6):890-6
10833195
Mol Biol Evol. 2006 Oct;23(10):1891-901
16818476
Nat Med. 2004 Mar;10(3):282-9
14770175
J Mol Evol. 2005 Mar;60(3):315-26
15871042
Brief Bioinform. 2009 Jan;10(1):97-109
18971241
PLoS Comput Biol. 2009 Nov;5(11):e1000581
19956739
J Virol. 2008 Jul;82(13):6434-46
18434400
Mol Biol Evol. 2001 May;18(5):691-9
11319253
Mol Biol Evol. 2007 Feb;24(2):388-97
17110464
Syst Biol. 2003 Oct;52(5):696-704
14530136
Genetics. 2000 May;155(1):431-49
10790415
PLoS One. 2007 Jun 06;2(6):e503
17551583
Philos Trans R Soc Lond B Biol Sci. 2008 Dec 27;363(1512):3941-53
18852098
PLoS One. 2010 Jul 30;5(7):e11230
20689581
Proc Natl Acad Sci U S A. 2008 Sep 9;105(36):13480-5
18768804
J Theor Biol. 1996 Nov 21;183(2):195-205
8977877
Mol Biol Evol. 2010 Mar;27(3):520-36
19864470
Mol Biol Evol. 2005 Mar;22(3):478-85
15509724
Mol Biol Evol. 1993 Nov;10(6):1396-401
8277861