Abstract:Background: Phylogeny estimation from aligned haplotype sequences has attracted more and more attention in the recent years due to its importance in analysis of many fine-scale genetic data. Its application fields range from medical research, to drug discovery, to epidemiology, to population dynamics. The literature on molecular phylogenetics proposes a number of criteria for selecting a phylogeny from among plausible alternatives. Usually, such criteria can be expressed by means of objective functions, and the phylogenies that optimize them are referred to as optimal. One of the most important estimation criteria is the parsimony which states that the optimal phylogeny T∗for a set H of n haplotype sequences over a common set of variable loci is the one that satisfies the following requirements: (i) it has the shortest length and (ii) it is such that, for each pair of distinct haplotypes hi,hj∈H, the sum of the edge weights belonging to the path from hi to hj in T∗ is not smaller than the observed number of changes between hi and hj. Finding the most parsimonious phylogeny for H involves solving an optimization problem, called the Most Parsimonious Phylogeny Estimation Problem (MPPEP), which is NP-hard in many of its versions. Results: In this article we investigate a recent version of the MPPEP that arises when input data consist of single nucleotide polymorphism haplotypes extracted from a population of individuals on a common genomic region. Specifically, we explore the prospects for improving on the implicit enumeration strategy of implicit enumeration strategy used in previous work using a novel problem formulation and a series of strengthening valid inequalities and preliminary symmetry breaking constraints to more precisely bound the solution space and accelerate implicit enumeration of possible optimal phylogenies. We present the basic formulation and then introduce a series of provable valid constraints to reduce the solution space. We then prove that these constraints can often lead to significant reductions in the gap between the optimal solution and its non-integral linear programming bound relative to the prior art as well as often substantially faster processing of moderately hard problem instances. Conclusion: We provide an indication of the conditions under which such an optimal enumeration approach is likely to be feasible, suggesting that these strategies are usable for relatively large numbers of taxa, although with stricter limits on numbers of variable sites. The work thus provides methodology suitable for provably optimal solution of some harder instances that resist all prior approaches.

POY version 5: phylogenetic analysis using dynamic homologies under multiple optimality criteria

POY version 4: phylogenetic analysis using dynamic homologies

VEHoP: A Versatile, Easy-to-use, and Homology-based Phylogenomic pipeline accommodating diverse sequences

A mixed integer linear programming model to reconstruct phylogenies from single nucleotide polymorphism haplotypes under the maximum parsimony criterion

BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis

An Alignment-free Method for Phylogeny Estimation using Maximum Likelihood

PASTA with many application-aware optimization criteria for alignment based phylogeny inference

Dollo-CDP: a polynomial-time algorithm for the clade-constrained large Dollo parsimony problem

DendroPy 5: a mature Python library for phylogenetic computing

ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R

HyPhy: hypothesis testing using phylogenies

ProtParCon: A Framework for Processing Molecular Data and Identifying Parallel and Convergent Amino Acid Replacements

Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion

PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies

Fast and exact gap-affine partial order alignment with POASTA

A MOSAIC of methods: Improving ortholog detection through integration of algorithmic diversity

A multiobjective memetic algorithm for PPI network alignment

Co-Phylog: an Assembly-Free Phylogenomic Approach for Closely Related Organisms

Differentiable Phylogenetics via Hyperbolic Embeddings with Dodonaphy

Hapsolutely: a user-friendly tool integrating haplotype phasing, network construction, and haploweb calculation

diverse-seq: an application for alignment-free selecting and clustering biological sequences