Treewidth-based algorithms for the small parsimony problem on networks
Background Phylogenetic reconstruction is one of the paramount challenges of contemporary bioinformatics. A subtask of existing tree reconstruction algorithms is modeled by the Small Parsimony problem: given a tree T and an assignment of character-states to its leaves, assign states to the internal...
Uložené v:
| Vydané v: | Algorithms for molecular biology Ročník 17; číslo 1; s. 1 - 31 |
|---|---|
| Hlavní autori: | , |
| Médium: | Journal Article |
| Jazyk: | English |
| Vydavateľské údaje: |
London
BioMed Central
20.08.2022
Springer Nature B.V BMC |
| Predmet: | |
| ISSN: | 1748-7188, 1748-7188 |
| On-line prístup: | Získať plný text |
| Tagy: |
Pridať tag
Žiadne tagy, Buďte prvý, kto otaguje tento záznam!
|
| Shrnutí: | Background
Phylogenetic reconstruction is one of the paramount challenges of contemporary bioinformatics. A subtask of existing tree reconstruction algorithms is modeled by the
Small Parsimony
problem: given a tree
T
and an assignment of character-states to its leaves, assign states to the internal nodes of
T
such as to minimize the
parsimony score
, that is, the number of edges of
T
connecting nodes with different states. While this problem is polynomial-time solvable on trees, the matter is more complicated if
T
contains reticulate events such as hybridizations or recombinations, i.e. when
T
is a network. Indeed, three different versions of the parsimony score on networks have been proposed and each of them is NP-hard to decide. Existing parameterized algorithms focus on combining the number
c
of possible character-states with the number of reticulate events (per biconnected component).
Results
We consider the parameter treewidth
t
of the underlying undirected graph of the input network, presenting dynamic programming algorithms for (slight generalizations of) all three versions of the parsimony problem on size-
n
networks running in times
c
t
n
O
(
1
)
,
(
3
c
)
t
n
O
(
1
)
, and
6
tc
n
O
(
1
)
, respectively. Our algorithms use a formulation of the treewidth that may facilitate formalizing treewidth-based dynamic programming algorithms on phylogenetic networks for other problems.
Conclusions
Our algorithms allow the computation of the three popular parsimony scores, modeling the evolutionary development of a (multistate) character on a given phylogenetic network of low treewidth. Our results subsume and improve previously known algorithm for all three variants. While our results rely on being given a “good” tree-decomposition of the input, encouraging theoretical results as well as practical implementations producing them are publicly available. We present a reformulation of tree decompositions in terms of “agreeing trees” on the same set of nodes. As this formulation may come more natural to researchers and engineers developing algorithms for phylogenetic networks, we hope to render exploiting the input network’s treewidth as parameter more accessible to this audience. |
|---|---|
| Bibliografia: | ObjectType-Article-1 SourceType-Scholarly Journals-1 ObjectType-Feature-2 content type line 14 content type line 23 |
| ISSN: | 1748-7188 1748-7188 |
| DOI: | 10.1186/s13015-022-00216-w |