Statistically-Consistent k-mer Methods for Phylogenetic Tree Reconstruction
Abstract
Frequencies of -mers in sequences are sometimes used as a basis for inferring phylogenetic trees without first obtaining a multiple sequence alignment. We show that a standard approach of using the squared-Euclidean distance between -mer vectors to approximate a tree metric can be statistically inconsistent. To remedy this, we derive model-based distance corrections for orthologous sequences without gaps, which lead to consistent tree inference. The identifiability of model parameters from -mer frequencies is also studied. Finally, we report simulations showing the corrected distance out-performs many other -mer methods, even when sequences are generated with an insertion and deletion process. These results have implications for multiple sequence alignment as well, since -mer methods are usually the first step in constructing a guide tree for such algorithms.
Cite
@article{arxiv.1511.01956,
title = {Statistically-Consistent k-mer Methods for Phylogenetic Tree Reconstruction},
author = {Elizabeth S. Allman and John A. Rhodes and Seth Sullivant},
journal= {arXiv preprint arXiv:1511.01956},
year = {2016}
}
Comments
25 pages, 9 figures figure added, to appear, JCB