Nøjgaard Nikolai, Geiß Manuela, Merkle Daniel, Stadler Peter F, Wieseke Nicolas, Hellmuth Marc
1Institute of Mathematics and Computer Science, University of Greifswald, Walther-Rathenau-Strasse 47, 17487 Greifswald, Germany.
2Department of Mathematics and Computer Science, University of Southern Denmark, Campusvej 55, 5230 Odense M, Denmark.
Algorithms Mol Biol. 2018 Feb 6;13:2. doi: 10.1186/s13015-018-0121-8. eCollection 2018.
In the absence of horizontal gene transfer it is possible to reconstruct the history of gene families from empirically determined orthology relations, which are equivalent to gene trees. Knowledge of the event labels considerably simplifies the problem of reconciling a gene tree with a species trees , relative to the reconciliation problem without prior knowledge of the event types. It is well-known that optimal reconciliations in the unlabeled case may violate time-consistency and thus are not biologically feasible. Here we investigate the mathematical structure of the event labeled reconciliation problem with horizontal transfer.
We investigate the issue of time-consistency for the event-labeled version of the reconciliation problem, provide a convenient axiomatic framework, and derive a complete characterization of time-consistent reconciliations. This characterization depends on certain weak conditions on the event-labeled gene trees that reflect conditions under which evolutionary events are observable at least in principle. We give an [Formula: see text]-time algorithm to decide whether a time-consistent reconciliation map exists. It does not require the construction of explicit timing maps, but relies entirely on the comparably easy task of checking whether a small auxiliary graph is acyclic. The algorithms are implemented in C++ using the boost graph library and are freely available at https://github.com/Nojgaard/tc-recon.
The combinatorial characterization of time consistency and thus biologically feasible reconciliation is an important step towards the inference of gene family histories with horizontal transfer from orthology data, i.e., without presupposed gene and species trees. The fast algorithm to decide time consistency is useful in a broader context because it constitutes an attractive component for all tools that address tree reconciliation problems.
在不存在水平基因转移的情况下,有可能从经验确定的直系同源关系(等同于基因树)重建基因家族的历史。相对于没有事件类型先验知识的和解问题,事件标签的知识极大地简化了使基因树与物种树和解的问题。众所周知,未标记情况下的最优和解可能违反时间一致性,因此在生物学上不可行。在这里,我们研究具有水平转移的事件标记和解问题的数学结构。
我们研究和解问题的事件标记版本的时间一致性问题,提供一个方便的公理框架,并推导时间一致和解的完整特征。这种特征取决于事件标记基因树上的某些弱条件,这些条件反映了进化事件至少在原则上可观察的条件。我们给出一个[公式:见正文]时间算法来决定是否存在时间一致的和解映射。它不需要构建明确的时间映射,而是完全依赖于检查一个小辅助图是否无环这个相对容易的任务。这些算法使用boost图库用C++实现,可在https://github.com/Nojgaard/tc-recon免费获取。
时间一致性以及因此生物学上可行的和解的组合特征是朝着从直系同源数据推断具有水平转移的基因家族历史迈出的重要一步,即无需预先假定基因树和物种树。用于决定时间一致性的快速算法在更广泛的背景下是有用的,因为它构成了所有解决树和解问题的工具的一个有吸引力的组件。