Statistical power of phylo-HMM for evolutionarily conserved element detection.
BACKGROUND
An important goal of comparative genomics is the identification of functional elements through conservation analysis. Phylo-HMM was recently introduced to detect conserved elements based on multiple genome alignments, but the method has not been rigorously evaluated.


RESULTS
We report here a simulation study to investigate the power of phylo-HMM. We show that the power of the phylo-HMM approach depends on many factors, the most important being the number of species-specific genomes used and evolutionary distances between pairs of species. This finding is consistent with results reported by other groups for simpler comparative genomics models. In addition, the conservation ratio of conserved elements and the expected length of the conserved elements are also major factors. In contrast, the influence of the topology and the nucleotide substitution model are relatively minor factors.


CONCLUSIONS
Our results provide for general guidelines on how to select the number of genomes and their evolutionary distance in comparative genomics studies, as well as the level of power we can expect under different parameter settings.
Fan
Xiaodan
X
Department of Statistics, Harvard University, Boston, MA, USA. xfan@fas.harvard.edu
Zhu
Jun
J
Schadt
Eric E
EE
Liu
Jun S
JS
D000465:Algorithms; D001483:Base Sequence; D002874:Chromosome Mapping / Q000379:methods*; D017124:Conserved Sequence / Q000235:genetics*; D003627:Data Interpretation, Statistical; D019143:Evolution, Molecular; D008390:Markov Chains; D008969:Molecular Sequence Data; D010802:Phylogeny; D016415:Sequence Alignment / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods*