MBE Advance Access published online on January 8, 2009
Molecular Biology and Evolution, doi:10.1093/molbev/msp003
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Research Article |
Learning to Count: Robust Estimates for Labeled Distances between Molecular Sequences
1 Department of Biomathematics, UCLA, Los Angeles, CA, 90095, USA
2 Department of Statistics, University of Washington, Seattle, WA, 98195, USA
3 Department of Biostatistics, UCLA, Los Angeles, CA, 90095, USA
4 Department of Human Genetics, UCLA, Los Angeles, CA, 90095, USA
* To whom correspondence should be addressed.
Corresponding author: Marc A. Suchard, Departments of Biomathematics, Biostatistics and Human Genetics, 6558 Gonda Building, 695 Charles E. Young Drive, South, Los Angeles, CA 90095-1766, email: msuchard{at}ucla.edu, telephone: (310) 825-7442, fax: (310) 825-8685
Received for publication July 1, 2008. Revision received November 3, 2008. Revision received December 29, 2008. Accepted for publication December 30, 2008.
Researchers routinely estimate distances between molecular sequences using continuous-time Markov chain models. We present a new method, robust counting, that protects against the possibly severe bias arising from model misspecification. We achieve this robustness by generalizing conventional distance estimation to incorporate the empirical distribution of site patterns found in the observed pairwise sequence alignment. Our flexible framework allows for computing distances based only on a subset of possible substitutions. From this, we show how to estimate labeled codon distances, such as expected numbers of synonymous or nonsynonymous substitutions. We present two simulation studies. The first compares the relative bias and variance of conventional and robust labeled nucleotide estimators. In the second simulation, we demonstrate that robust counting furnishes accurate synonymous and nonsynonymous distance estimates based only on easy-to-fit models of nucleotide substitution, bypassing the need for computationally expensive codon models. We conclude with three empirical examples. In the first two examples, we investigate the evolutionary dynamics of the influenza A hemagglutinin gene using labeled codon distances. In the final example, we demonstrate the advantages of using robust synonymous distances to alleviate the effect of convergent evolution on phylogenetic analysis of an HIV transmission network.
Key Words: robust counting • labeled codon distance • empirical distribution • Markov chain substitution model
These authors wish to be considered as joint first authors.