Abstract:Reconstructing the tree of life is an essential task in evolutionary biology. It demands accurate phylogenetic inference for both extant and extinct organisms, the latter being almost entirely dependent on morphological data. While parsimony methods have traditionally dominated the field of morphological phylogenetics, a rapidly growing number of studies are now employing probabilistic methods (maximum likelihood and Bayesian inference). The present-day toolkit of probabilistic methods offers varied software with distinct algorithms and assumptions for reaching global optimality. However, benchmark performance assessments of different software packages for the analyses of morphological data, particularly in the era of big data, are still lacking. Here, we test the performance of four major probabilistic software under variable taxonomic sampling and missing data conditions: the Bayesian inference-based programs MrBayes and RevBayes, and the maximum likelihood-based IQ-TREE and RAxML. We evaluated software performance by calculating the distance between inferred and true trees using a variety of metrics, including Robinson-Foulds (RF), Matching Splits (MS), and Kuhner-Felsenstein (KF) distances. Our results show that increased taxonomic sampling improves accuracy, precision, and resolution of reconstructed topologies across all tested probabilistic software applications and all levels of missing data. Under the RF metric, Bayesian inference applications were the most consistent, accurate, and robust to variation in taxonomic sampling in all tested conditions, especially at high levels of missing data, with little difference in performance between the two tested programs. The MS metric favored more resolved topologies that were generally produced by IQ-TREE. Adding more taxa dramatically reduced performance disparities between programs. Importantly, our results suggest that the RF metric penalizes incorrectly resolved nodes (false positives) more severely than the MS metric, which instead tends to penalize polytomies. If false positives are to be avoided in systematics, Bayesian inference should be preferred over maximum likelihood for the analysis of morphological data.

IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era

IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies

W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis

RAxML-NG: a fast, scalable and user-friendly tool for maximum likelihood phylogenetic inference

DecentTree: scalable Neighbour-Joining for the genomic era

Terrace Aware Phylogenomic Inference from Supermatrices

MrBayes 3.2: Efficient Bayesian Phylogenetic Inference and Model Choice Across a Large Model Space

RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies

PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies

MixtureFinder: Estimating DNA mixture models for phylogenetic analyses

The Influence of the Number of Tree Searches on Maximum Likelihood Inference in Phylogenomics.

The SplitsTree App: interactive analysis and visualization using phylogenetic trees and networks

BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis

From Easy to Hopeless-Predicting the Difficulty of Phylogenetic Analyses

Sciphy: A Bayesian phylogenetic framework using sequential genetic lineage tracing data.

Evaluating the Performance of Probabilistic Algorithms for Phylogenetic Analysis of Big Morphological Datasets: A Simulation Study

MIKE: an ultrafast, assembly- and alignment-free approach for phylogenetic tree construction

TreeFormer: A transformer-based tree rearrangement operation for phylogenetic reconstruction

GeneSupport Maximum Gene-Support Tree Approach to Species Phylogeny Inference

DAMBE7: New and Improved Tools for Data Analysis in Molecular Biology and Evolution

Phylo2Vec: a vector representation for binary trees