Abstract:Reconstructing the evolutionary history of different groups of organisms provides insight into how life originated and diversified on Earth. Phylogenetic trees are commonly used to estimate this evolutionary history. Within Bayesian phylogenetics a major step in estimating a tree is in choosing an appropriate model of character evolution. While the most common character data used is molecular sequence data, morphological data remains a vital source of information. The use of morphological characters allows for the incorporation fossil taxa, and despite advances in molecular sequencing, continues to play a significant role in neontology. Moreover, it is the main data source that allows us to unite extinct and extant taxa directly under the same generating process. We therefore require suitable models of morphological character evolution, the most common being the Mk Lewis model. While it is frequently used in both palaeobiology and neontology, it is not known whether the simple Mk substitution model, or any extensions to it, provide a sufficiently good description of the process of morphological evolution. In this study we investigate the impact of different morphological models on empirical tetrapod data sets. Specifically, we compare unpartitioned Mk models with those where characters are partitioned by the number of observed states, both with and without allowing for rate variation across sites and accounting for ascertainment bias. We show that the choice of substitution model has an impact on both topology and branch lengths, highlighting the importance of model choice. Through simulations, we validate the use of the model adequacy approach, posterior predictive simulations, for choosing an appropriate model. Additionally, we compare the performance of model adequacy with Bayesian model selection. We demonstrate how model selection approaches based on marginal likelihoods are not appropriate for choosing between models with partition schemes that vary in character state space (i.e., that vary in Q-matrix state size). Using posterior predictive simulations, we found that current variations of the Mk model are often performing adequately in capturing the evolutionary dynamics that generated our data. We do not find any preference for a particular model extension across multiple data sets, indicating that there is no 'one size fits all' when it comes to morphological data and that careful consideration should be given to choosing models of discrete character evolution. By using suitable models of character evolution, we can increase our confidence in our phylogenetic estimates, which should in turn allow us to gain more accurate insights into the evolutionary history of both extinct and extant taxa.

Automatic generation of evolutionary hypotheses using mixed Gaussian phylogenetic models

Bayesian inference of mixed Gaussian phylogenetic models

Efficient Bayesian Inference of General Gaussian Models on Large Phylogenetic Trees

Inferring ancestral states without assuming neutrality or gradualism using a stable model of continuous character evolution

A discrete character evolution model for phylogenetic comparative biology with Γ-distributed rate heterogeneity among branches of the tree

Infinite Mixture Models for Improved Modeling of Across-Site Evolutionary Variation

Assessing the Adequacy of Morphological Models using Posterior Predictive Simulations

Multi-Response Phylogenetic Mixed Models: Concepts and Application

Fast fitting of phylogenetic mixed effects models

Evaluation of the Models Handling Heterotachy in Phylogenetic Inference

Multivariate Trait Evolution: Models for the Evolution of the Quantitative Genetic G-Matrix on Phylogenies.

Probabilistic Graphical Model Representation in Phylogenetics

Automatic Discovery of Optimal Discrete Character Models

An independent base composition of each rate class for improved likelihood-based phylogeny estimation: The 5rf model

EREM: Parameter Estimation and Ancestral Reconstruction by Expectation-Maximization Algorithm for a Probabilistic Model of Genomic Binary Characters Evolution

Bayesian inference of phylogenetic distances: revisiting the eigenvalue approach

Assessing the impact of character evolution models on phylogenetic and macroevolutionary inferences from fossil data

A Unifying Model for the Analysis of Phenotypic, Genetic and Geographic Data

Assessing the Adequacy of Morphological Models used in Palaeobiology

Expectation-Maximization enables Phylogenetic Dating under a Categorical Rate Model

A Bayesian Supertree Model for Genome-Wide Species Tree Reconstruction