Abstract:The combination of a time reversible Markov process with a 'hidden' mixture of gamma distributed relative site rates plus invariant sites have become the most favoured options for likelihood and other probabilistic models of nucleotide evolution (e.g., tr4gi which approximates a gamma with four rate classes). However, these models assume a homogeneous and stationary distribution of nucleotide (character or base) frequencies. Here, we explore the potential benefits and pitfalls of allowing each rate category (rate class) of a 4gi mixture model to have its own base frequencies. This is achieved by starting each of the five rate classes, at the tree's root, with its own free choice of nucleotide frequencies to create a 4gi5rf model or a 5rf model in shorthand We assess the practical identifiability of this approach with a BEAST 2 implementation, aiming to determine if it can accurately estimate credibility intervals and expected values for a wide range of plausible parameter values. Practical identifiability, as distinguished from mathematical identifiability, gauges the model's ability to identify parameters in real-world scenarios, as opposed to theoretically with infinite data. One of the most common types of phylogenetic data is mitochondrial DNA (mtDNA) protein coding sequence. It is often assumed current models analyse robustly such data and that higher likelihood/posterior probability models do better. However, this abstract shows that vertebrate mtDNA remains a very difficult type of data to fully model, and that dramatically higher likelihoods do not mean a model is measurably more accurate with respect to recovering key parameters of biological interest (e.g., monophyletic groups, their support and their ages). The 4gi5rf model considerably improves marginal likelihoods and seems to reverse some apparent errors exacerbated by the 4gi model, while introducing others. Problems appear to be linked to non-stationary DNA repair processes that alter the mutation/substitution spectra across lineages and time. We also show such problems are not unique to mtDNA and are encountered in analysing nuclear sequences. Non-stationarity of DNA repair processes mutation/substitution spectra thus pose an active challenge to obtaining reliable inferences of relationships and divergence times near the root of placental mammals, for example. An open source implementation is available under the LGPL 3.0 license in the beastbooster package for BEAST 2, available from https://github.com/rbouckaert/beastbooster.

Estimating the rate of molecular evolution: incorporating non-contemporaneous sequences into maximum likelihood phylogenies

Estimating evolutionary rates using time-structured data: a general comparison of phylogenetic methods

Bayesian Inference of Evolutionary Histories under Time-Dependent Substitution Rates

Dating the time of viral subtype divergence

Bayesian inference of phylogenetic distances: revisiting the eigenvalue approach

Tree and Rate Estimation by Local Evaluation of Heterochronous Nucleotide Data

Molecular evolution of dengue virus: a Bayesian approach using 1581 whole-genome sequences from January 1944 to July 2022

Model-based inference of punctuated molecular evolution

Examining the molecular clock hypothesis for the contemporary evolution of the rabies virus

Scalable Bayesian divergence time estimation with ratio transformations

Mind your Ps: A probabilistic model to aid the interpretation of molecular epidemiology data

Maximum likelihood estimators for scaled mutation rates in an equilibrium mutation-drift model

Estimating the rate of adaptive molecular evolution in the presence of slightly deleterious mutations and population size change

Combined Sum of Squares Penalties for Molecular Divergence Time Estimation

A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences

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

Estimating Diversity Through Time using Molecular Phylogenies: Old and Species-Poor Frog Families are the Remnants of a Diverse Past

Estimation of demography and mutation rates from one million haploid genomes

Rate variation and recurrent sequence errors in pandemic-scale phylogenetics

PhyloAcc-GT: A Bayesian method for inferring patterns of substitution rate shifts on targeted lineages accounting for gene tree discordance

Relaxed Phylogenetics and Dating with Confidence