Correcting for sequencing error in maximum likelihood phylogeny inference.
Accurate phylogenies are critical to taxonomy as well as studies of speciation processes and other evolutionary patterns. Accurate branch lengths in phylogenies are critical for dating and rate measurements. Such accuracy may be jeopardized by unacknowledged sequencing error. We use simulated data to test a correction for DNA sequencing error in maximum likelihood phylogeny inference. Over a wide range of data polymorphism and true error rate, we found that correcting for sequencing error improves recovery of the branch lengths, even if the assumed error rate is up to twice the true error rate. Low error rates have little effect on recovery of the topology. When error is high, correction improves topological inference; however, when error is extremely high, using an assumed error rate greater than the true error rate leads to poor recovery of both topology and branch lengths. The error correction approach tested here was proposed in 2004 but has not been widely used, perhaps because researchers do not want to commit to an estimate of the error rate. This study shows that correction with an approximate error rate is generally preferable to ignoring the issue.
Department of Genome Sciences, University of Washington, Seattle, Washington 98195-5065 mkkuhner@uw.com.
Mary K
MK
Kuhner
Department of Genome Sciences, University of Washington, Seattle, Washington 98195-5065.
James
J
McGill
D000818:Animals; D019295:Computational Biology / Q000379:methods*; D006801:Humans; D016013:Likelihood Functions; D008957:Models, Genetic; D015233:Models, Statistical; D010802:Phylogeny; D011336:Probability; D017422:Sequence Analysis, DNA; D012984:Software
maximum likelihood; phylogeny inference; sequencing error