Do branch lengths help to locate a tree in a phylogenetic network?

TL;DR

This paper investigates whether incorporating branch length information in phylogenetic networks improves the accuracy of locating specific trees within the network, revealing computational challenges and efficient solutions for certain cases.

Contribution

It introduces new versions of the tree containment problem that include branch lengths and provides algorithms for special cases where the problem is computationally feasible.

Findings

Branch length information can improve tree localization within networks.

The general problem is computationally hard.

Efficient algorithms exist for networks of limited complexity.

Abstract

Phylogenetic networks are increasingly used in evolutionary biology to represent the history of species that have undergone reticulate events such as horizontal gene transfer, hybrid speciation and recombination. One of the most fundamental questions that arise in this context is whether the evolution of a gene with one copy in all species can be explained by a given network. In mathematical terms, this is often translated in the following way: is a given phylogenetic tree contained in a given phylogenetic network? Recently this tree containment problem has been widely investigated from a computational perspective, but most studies have only focused on the topology of the phylo- genies, ignoring a piece of information that, in the case of phylogenetic trees, is routinely inferred by evolutionary analyses: branch lengths. These measure the amount of change (e.g., nucleotide substitutions) that has occurred along each branch of the phylogeny. Here, we study a number of versions of the tree containment problem that explicitly account for branch lengths. We show that, although length information has the potential to locate more precisely a tree within a network, the problem is computationally hard in its most general form. On a positive note, for a number of special cases of biological relevance, we provide algorithms that solve this problem efficiently. This includes the case of networks of limited complexity, for which it is possible to recover, among the trees contained by the network with the same topology as the input tree, the closest one in terms of branch lengths.