Accuracy and power of statistical methods for detecting adaptive evolution in protein coding sequences and for identifying positively selected sites.
The parsimony method of Suzuki and Gojobori (1999) and the maximum likelihood method developed from the work of Nielsen and Yang (1998) are two widely used methods for detecting positive selection in homologous protein coding sequences. Both methods consider an excess of nonsynonymous (replacement) substitutions as evidence for positive selection. Previously published simulation studies comparing the performance of the two methods show contradictory results. Here we conduct a more thorough simulation study to cover and extend the parameter space used in previous studies. We also reanalyzed an HLA data set that was previously proposed to cause problems when analyzed using the maximum likelihood method. Our new simulations and a reanalysis of the HLA data demonstrate that the maximum likelihood method has good power and accuracy in detecting positive selection over a wide range of parameter values. Previous studies reporting poor performance of the method appear to be due to numerical problems in the optimization algorithms and did not reflect the true performance of the method. The parsimony method has a very low rate of false positives but very little power for detecting positive selection or identifying positively selected sites.
Wong
Wendy S W
WS
Department of Biological Statistics and Computational Biology, Cornell University, Ithaca, New York 14850, USA. sww8@cornell.edu
Yang
Ziheng
Z
Goldman
Nick
N
Nielsen
Rasmus
R
D000465:Algorithms; D000596:Amino Acids / Q000235:genetics*; D000818:Animals; D003062:Codon / Q000235:genetics*; D003198:Computer Simulation; D019143:Evolution, Molecular; D006680:HLA Antigens / Q000737:chemistry / Q000235:genetics; D006801:Humans; D008957:Models, Genetic; D015233:Models, Statistical; D016366:Open Reading Frames / Q000235:genetics*; D010802:Phylogeny; D012641:Selection, Genetic
D000596:Amino Acids; D003062:Codon; D006680:HLA Antigens