Abstract
The balanced minimal evolution (BME) method of creating phylogenetic trees can be formulated as a linear programming problem, minimizing an inner product over the vertices of the BME polytope. In this paper we undertake the project of describing the facets of this polytope. We classify and identify the combinatorial structure and geometry (facet inequalities) of all the facets in dimensions up to five, and classify even more facets in all dimensions. A full set of facet inequalities would allow a full implementation of the simplex method for finding the BME tree—although there are reasons to think this an unreachable goal. However, our results provide the crucial first steps for a more likely-to-be-successful program: finding efficient relaxations of the BME polytope.
Similar content being viewed by others
References
Billera LJ, Sarangarajan A (1996) All 0–1 polytopes are traveling salesman polytopes. Combinatorica 16(2):175–188. doi:10.1007/BF01844844
Catanzaro D, Labbé M, Pesenti R, Salazar-González JJ (2012) The balanced minimum evolution problem. INFORMS J Comput 24(2):276–294. doi:10.1287/ijoc.1110.0455
Day WHE (1987) Computational complexity of inferring phylogenies from dissimilarity matrices. Bull Math Biol 49(4):461–467. doi:10.1016/S0092-8240(87)80007-1
Desper R, Gascuel O (2002) Fast and accurate phylogeny reconstruction algorithms based on the minimum-evolution principle. J Comp Biol 9(5):687–705
Eickmeyer K, Huggins P, Pachter L, Yoshida R (2008) On the optimality of the neighbor-joining algorithm. Algorithms Mol Biol 3:5. doi:10.1186/1748-7188-3-5
Fiorini S, Joret G (2012) Approximating the balanced minimum evolution problem. Oper Res Lett 40(1):31–35. doi:10.1016/j.orl.2011.10.003
Forcey S (2014) Dear NSA: Long-term security depends on freedom. Notices of the AMS 61(1):7–8
Gascuel O, Steel M (2006) Neighbor-joining revealed. Mol Biol Evol 23:1997–2000
Gawrilow E, Joswig M (2000) Polymake: a framework for analyzing convex polytopes. In: Kalai G, Ziegler GM (eds) Polytopes—combinatorics and computation. Birkhäuser, Basel, pp 43–74
Haws DC, Hodge TL, Yoshida R (2011) Optimality of the neighbor joining algorithm and faces of the balanced minimum evolution polytope. Bull Math Biol 73(11):2627–2648. doi:10.1007/s11538-011-9640-x
Huggins P (2008) Polytopes in computational biology. PhD Dissertation, UC Berkeley
Pauplin Y (2000) Direct calculation of a tree length using a distance matrix. J Mol Evol 51(1):41–47
Saitou N, Nei M (1987) The neighbor joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4:406–425
Acknowledgments
We would like to thank the referees, whose suggestions helped to improve the readability of this paper. The first author would like to thank the organizers and participants in the Working Group for geometric approaches to phylogenetic tree reconstructions, at the NSF/CBMS Conference on Mathematical Phylogeny held at Winthrop University in June–July 2014. Especially helpful were conversations with Ruriko Yoshida, Terrell Hodge and Matt Macauley. The first author would also like to thank the American Mathematical Society and the Mathematical Sciences Program of the National Security Agency for supporting this research through Grant H98230-14-0121. (This manuscript is submitted for publication with the understanding that the United States Government is authorized to reproduce and distribute reprints.) The first author’s specific position on the NSA is published in Forcey (2014). Suffice it to say here that he appreciates NSA funding for open research and education, but encourages reformers of the NSA who are working to ensure that protections of civil liberties keep pace with intelligence capabilities.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Forcey, S., Keefe, L. & Sands, W. Facets of the balanced minimal evolution polytope. J. Math. Biol. 73, 447–468 (2016). https://doi.org/10.1007/s00285-015-0957-1
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00285-015-0957-1