Statistical properties of the methods for detecting positively selected amino acid sites.
Parsimony and Bayesian methods have been developed for detecting positively selected amino acid sites. It has been reported that the parsimony method is generally conservative. In contrast, the Bayesian method is known to identify more positively selected sites than the parsimony method, especially when the number of sequences analyzed is small, although the interpretation of results obtained from the former method is controversial. Here I show that the likelihood-ratio test (LRT) of the Bayesian method corresponds to the parsimony method with window analysis, by analyzing the nucleotide sequences encoding the trans-activator (tax) gene of human T-cell lymphotropic virus type I (HTLV-I). It is also indicated that in the parsimony method, the test of selective neutrality using the binomial probability tends to be conservative, but the Monte Carlo simulation is useful for solving this problem. In addition, in the Bayesian method, the bootstrap method appears to produce similar results to the LRT. This information may be useful for improving the methods for detecting positively selected amino acid sites.
Suzuki
Yoshiyuki
Y
Center for Information Biology and DNA Data Bank of Japan, National Institute of Genetics, Mishima-shi, Shizuoka-ken 411-8540, Japan. yossuzuk@lab.nig.ac.jp
D000595:Amino Acid Sequence; D000596:Amino Acids / Q000235:genetics*; D001483:Base Sequence; D001499:Bayes Theorem; D001665:Binding Sites; D003062:Codon; D003198:Computer Simulation; D016208:Databases, Factual; D016356:Gene Products, tax / Q000737:chemistry*; D015368:Human T-lymphotropic virus 1 / Q000235:genetics*; D006801:Humans; D016013:Likelihood Functions; D008969:Molecular Sequence Data; D009010:Monte Carlo Method; D010802:Phylogeny; D012641:Selection, Genetic; D017386:Sequence Homology, Amino Acid; D013223:Statistics as Topic
D000596:Amino Acids; D003062:Codon; D016356:Gene Products, tax