Cubic time algorithms of amalgamating gene trees and building evolutionary scenarios

<p>Abstract</p> <p>Background</p> <p>A long recognized problem is the inference of the supertree <it>S</it> that amalgamates a given set {<it>G</it>_<it>j</it>} of trees <it>G</it>_<it>j</it>, with leaves in each <it>G</it>_<it>j</it> being assigned homologous elements.</p> <p>We ground on an approach to find the tree <it>S</it> by minimizing the total cost of mappings <it>α</it>_<it>j</it> of individual gene trees <it>G</it>_<it>j</it> into <it>S</it>. Traditionally, this cost is defined basically as a sum of duplications and gaps in each <it>α</it>_<it>j</it>. The classical problem is to minimize the total cost, where <it>S</it> runs over the set of all trees that contain an exhaustive non-redundant set of species from all input <it>G</it>_<it>j</it>.</p> <p>Results</p> <p>We suggest a reformulation of the classical <it>NP</it>-hard problem of building a supertree in terms of the global minimization of the same cost functional but only over species trees <it>S</it> that consist of clades belonging to a fixed set <it>P</it> (e.g., an exhaustive set of clades in all <it>G</it>_<it>j</it>). We developed a deterministic solving algorithm with a low degree polynomial (typically cubic) time complexity with respect to the size of input data.</p> <p>We define an extensive set of elementary evolutionary events and suggest an original definition of mapping <it>β</it> of tree <it>G</it> into tree <it>S</it>. We introduce the cost functional <it>c</it>(<it>G</it>, <it>S</it>, <it>f</it> ) and define the mapping <it>β</it> as the global minimum of this functional with respect to the variable <it>f</it>, in which sense it is a generalization of classical mapping <it>α</it>.</p> <p>We suggest a reformulation of the classical <it>NP</it>-hard mapping (reconciliation) problem by introducing time slices into the species tree <it>S</it> and present a cubic time solving algorithm to compute the mapping <it>β</it>. We introduce two novel definitions of the evolutionary scenario based on mapping <it>β</it> or a random process of gene evolution along a species tree.</p> <p>Conclusions</p> <p>Developed algorithms are mathematically proved, which justifies the following statements. The supertree building algorithm finds exactly the global minimum of the total cost if only gene duplications and losses are allowed and the given sets of gene trees satisfies a certain condition. The mapping algorithm finds exactly the minimal mapping <it>β</it>, the minimal total cost and the evolutionary scenario as a minimum over all possible distributions of elementary evolutionary events along the edges of tree <it>S</it>.</p> <p>The algorithms and their effective software implementations provide useful tools in many biological studies. They facilitate processing of voluminous tree data in acceptable time still largely avoiding heuristics. Performance of the tools is tested with artificial and prokaryotic tree data.</p> <p>Reviewers</p> <p>This article was reviewed by Prof. Anthony Almudevar, Prof. Alexander Bolshoy (nominated by Prof. Peter Olofsson), and Prof. Marek Kimmel.</p>