Bayesian Phylogenetic Model Selection Using Reversible Jump Markov Chain Monte Carlo

Abstract

A common problem in molecular phylogenetics is choosing a model of DNA substitution that does a good job of explaining the DNA sequence alignment without introducing superfluous parameters. A number of methods have been used to choose among a small set of candidate substitution models, such as the likelihood ratio test, the Akaike Information Criterion (AIC), the Bayesian Information Criterion (BIC), and Bayes factors. Current implementations of any of these criteria suffer from the limitation that only a small set of models are examined, or that the test does not allow easy comparison of non-nested models. In this article, we expand the pool of candidate substitution models to include all possible time-reversible models. This set includes seven models that have already been described. We show how Bayes factors can be calculated for these models using reversible jump Markov chain Monte Carlo, and apply the method to 16 DNA sequence alignments. For each data set, we compare the model with the best Bayes factor to the best models chosen using AIC and BIC. We find that the best model under any of these criteria is not necessarily the most complicated one; models with an intermediate number of substitution types typically do best. Moreover, almost all of the models that are chosen as best do not constrain a transition rate to be the same as a transversion rate, suggesting that it is the transition/ transversion rate bias that plays the largest role in determining which models are selected. Importantly, the reversible jump Markov chain Monte Carlo algorithm described here allows estimation of phylogeny (and other phylogenetic model parameters) to be performed while accounting for uncertainty in the model of DNA substitution.

Hybridization of powdery mildew strains gives rise to pathogens on novel agricultural crop species

“…Variation in substitution rates across sites was modeled with a discretized (four categories) gamma distribution 39,40 . The chains have been let free to sample all models of the GTR model family using reversible-jump Monte Carlo Markov chain 44 . The divergence time between B.g.…”

Phylogeny of hoplocercine lizards (Squamata: Iguania) with estimates of relative divergence times

“…The individual likelihoods of the various models at each site are then summed, after weighting the models by the probability that they apply to that site, which can be estimated from the data Meade, 2004, 2005). The model of evolution (i.e., number of Q matrices, each conforming to the GTR + C model) was selected using reversible jump Markov Chain Monte Carlo (Huelsenbeck et al, 2004) as implemented in a test version of BayesPhylogenies provided by its authors. One run consisting of 10 7 generations and four Markov chains was conducted in BayesPhylogenies.…”

New isolates of Trichoderma antagonistic to Sclerotinia sclerotiorum

“…A phylogenetic hypothesis was calculated using MrBayes (v3.2.5; Ronquist & Huelsenbeck 2003). The Bayesian analysis was conducted using default program settings except that reversible-jump MCMC (Huelsenbeck et al 2004) was used to optimize base substitution model parameters within the general time reversible (GTR) framework. MCMC sampling was run for 2 million generations, which allowed the standard deviation of split frequencies to fall to 0.006853, strongly indicating stable convergence.…”