Computing the Rooted Triplet Distance Between Phylogenetic Networks

The rooted triplet distance measures the structural dissimilarity of two phylogenetic trees or networks by counting the number of rooted trees with exactly three leaf labels that occur as embedded subtrees in one, but not both of them. Suppose that N-1 = (V-1, E-1) and N-2 = (V-2, E-2) are rooted phylogenetic networks over a common leaf label set of size., that N-i has level k(i) and maximum in-degree d(i) for i. {1, 2}, and that the networks' out-degrees are unbounded. Denote n = max(vertical bar V-1 vertical bar, vertical bar V-2 vertical bar|), m = max(vertical bar E-1 vertical bar, vertical bar E-2 vertical bar), k = max(k(1), k(2)), and d = max(d(1), d(2)). Previous work has shown how to compute the rooted triplet distance between N-1 and N-2 in O(. log.) time in the special case k = 1. For k > 1, no efficient algorithms are known; a trivial approach leads to a running time of Omega(n(7) lambda(3)) and the only existing non-trivial algorithm imposes restrictions on the networks' in- and out-degrees (in particular, it does not work when non-binary nodes are allowed). In this paper, we develop two new algorithms that have no such restrictions. Their running times are O(n(2)m +lambda(3)) and O(m + k(3)d(3) lambda + lambda(3)), respectively. We also provide implementations of our algorithms and evaluate their performance in practice. This is the first publicly available software for computing the rooted triplet distance between unrestricted networks of arbitrary levels.