Skip Navigation


MBE Advance Access originally published online on April 8, 2008
Molecular Biology and Evolution 2008 25(7):1253-1256; doi:10.1093/molbev/msn083
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
25/7/1253    most recent
msn083v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Posada, D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Posada, D.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2008. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oxfordjournals.org

Letters

jModelTest: Phylogenetic Model Averaging

David Posada

Departamento de Genética, Bioquímica e Inmunología, Facultad de Biología, Universidad de Vigo, Vigo, Spain

E-mail: dposada{at}uvigo.es.


   Abstract
TOP
 Abstract
Introduction
 Model Selection with jModelTest
 Model Selection Uncertainty
 Conclusions
 Acknowledgements
 References
 
jModelTest is a new program for the statistical selection of models of nucleotide substitution based on "Phyml" (Guindon and Gascuel 2003Go. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 52:696–704.). It implements 5 different selection strategies, including "hierarchical and dynamical likelihood ratio tests," the "Akaike information criterion," the "Bayesian information criterion," and a "decision-theoretic performance-based" approach. This program also calculates the relative importance and model-averaged estimates of substitution parameters, including a model-averaged estimate of the phylogeny. jModelTest is written in Java and runs under Mac OSX, Windows, and Unix systems with a Java Runtime Environment installed. The program, including documentation, can be freely downloaded from the software section at http://darwin.uvigo.es.

Key Words: model selection • likelihood ratio tests • AIC • BIC • performance-based selection • statistical phylogenetics


   Introduction
TOP
 Abstract
 Introduction
Model Selection with jModelTest
 Model Selection Uncertainty
 Conclusions
 Acknowledgements
 References
 
Models of nucleotide substitution allow for the calculation of probabilities of change between nucleotides along the branches of a phylogenetic tree. The use of a particular substitution model may change the outcome of the phylogenetic analysis (e.g., Buckley 2002Go; Buckley and Cunningham 2002Go; Lemmon and Moriarty 2004Go), and statistical model selection has become an essential step for the estimation of phylogenies from DNA sequence alignments. In-depth reviews about model selection in phylogenetics are available elsewhere (Johnson and Omland 2003Go; Posada and Buckley 2004Go; Sullivan and Joyce 2005Go). Indeed, the performance of different model selection strategies has been the subject of active research (Posada 2001Go; Posada and Crandall 2001Go; Pol 2004Go; Abdo et al. 2005Go; Alfaro and Huelsenbeck 2006Go).

Several programs already exist for the statistical selection of models of nucleotide substitution (e.g., Nylander 2004Go; Keane et al. 2006Go). Among these, Modeltest (Posada and Crandall 1998Go) has been one of the most popular. This note describes a new program called jModelTest that supersedes Modeltest in several aspects. jModelTest allows for the definition of restricted sets of candidate models (table 1), implements customizable "hierarchical likelihood ratio tests" (hLRTs) (Frati et al. 1997Go; Huelsenbeck and Crandall 1997Go; Sullivan et al. 1997Go) and "dynamic likelihood ratio tests" (dLRTs) (Posada and Crandall 2001Go), provides a rank of models according to the "Akaike Information Criterion" (AIC) (Akaike 1973Go), to the "Bayesian Information Criterion" (BIC) (Schwarz 1978Go) or to a "decision-theoretic performance-based" approach (DT) (Minin et al. 2003Go) (table 2), calculates the relative importance of every parameter, and computes model-averaged estimates of these, including a model-averaged estimate of the tree topology (Posada and Buckley 2004Go).


View this table:
[in this window]
[in a new window]

  		Table 1 Substitution Models Available in jModelTest

 

View this table:
[in this window]
[in a new window]

  		Table 2 Model Selection Strategies Implemented in jModelTest

 

   Model Selection with jModelTest
TOP
 Abstract
 Introduction
 Model Selection with jModelTest
Model Selection Uncertainty
 Conclusions
 Acknowledgements
 References
 
jModelTest is essentially a front-end to a computational pipeline that takes advantage of existing programs for running different tasks. Basically, this pipeline (fig. 1) includes:
  • "ReadSeq" (Gilbert 2007Go): for conversion among different DNA sequence alignment formats.
  • "Phyml" (Guindon and Gascuel 2003Go): for the likelihood calculations, including estimates of model parameters and trees.
  • "Ted" (D. Posada): to compute Euclidean distances between trees for performance-based model selection.
  • "Consense" (from the PHYLIP package) (Felsenstein 2005Go): to calculate weighted and strict consensus trees representing model-averaged phylogenies.


Figure 1
View larger version (31K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		FIG. 1.— jModelTest pipeline. Alignments are loaded using the ReadSeq library (Gilbert 2007Go). Likelihood calculations, including estimates of model parameters and trees, are carried out with Phyml (Guindon and Gascuel 2003Go). A custom program called Ted (D. Posada) is used to compute Euclidean distances between trees for performance-based model selection (DT), whereas Consense (Felsenstein 2005Go) is used to calculate weighted and strict consensus trees representing model-averaged phylogenies.

 
Likelihood Calculations
Likelihood calculations, including model parameters and tree estimates, are carried out with Phyml (Guindon and Gascuel 2003Go). The tree topology used in these calculations can be the same across models (fixed) or optimized for each one. Fixed tree topologies can be estimated with the BIONJ algorithm (Gascuel 1997Go) upon JC distances (Jukes and Cantor 1969Go) or user-defined. Alternatively, a BIONJ or an ML tree can be estimated under each model. In all cases, branch lengths are estimated and counted as parameters.

Custom Set of Models
Currently, there are 11 different nucleotide substitution schemes implemented in jModelTest, which combined with equal or unequal base frequencies (+F), a proportion of invariable sites (+I), and rate variation among sites (+G), result in 88 distinct models (table 1). The program offers the possibility of defining to a reasonable extent which models are included in the candidate set.

Sequential Likelihood Ratio Tests
A series of likelihood ratio tests (LRTs) can be implemented under a particular hierarchy (hLRTs), in which the user can specify their order, and whether parameters are added (forward selection) or removed (backward selection). Alternatively, the order of the LRTs can be set dynamically (dLRTs) (Posada and Crandall 2001Go), by comparing the current model with the one that is one hypothesis away and provides the largest increase (under forward selection) or smallest decrease (under backward selection) in likelihood. The hLRTs and dLRTs will be available only if the likelihood scores were calculated upon a fixed topology, due to the nesting requirement of the {chi}2 approximation.

Information Criteria
The program implements 3 different information criteria: the AIC (Akaike 1973Go), the BIC (Schwarz 1978Go), and a performance-based approach based on decision theory (DT) (Minin et al. 2003Go). Under the AIC framework, there is also the possibility of using a corrected version for small samples (AICc) (Sugiura 1978Go; Hurvich and Tsai 1989Go), instead of the standard AIC. In this case, sample size has to be specified, which by default is approximated as the number of sites in the alignment (note that the sample size of an alignment is presently an unknown quantity).


   Model Selection Uncertainty
TOP
 Abstract
 Introduction
 Model Selection with jModelTest
 Model Selection Uncertainty
Conclusions
 Acknowledgements
 References
 
The AIC, BIC, and DT methods assign a score to each model in the candidate set, therefore providing an objective function to rank them. Using the differences in scores, the program can calculate a measure of model support called AIC or BIC weights (Burnham and Anderson 2003Go). For the DT scores, this calculation is not as straightforward, and right now a very gross approach is used instead, where the DT weights are the rescaled reciprocal DT scores. Confidence intervals (CIs) can be defined according to the cumulative weights, including a specified fraction of the models. When the CI includes only partially a given model, this model is included (yes/no) in the CI with a probability equal to the fraction included.

Parameter Importance and Model-Averaged Estimates
The program can also calculate the relative importance of every parameter of the substitution model and model-averaged estimates of these, using all the models in the candidate set, or a fraction included in a particular CI (see Posada and Buckley 2004Go).

Model-Averaged Phylogenies
jModelTest is able to compute an average estimate of the tree topology by building a consensus of the maximum likelihood (ML) trees for every model in the candidate set, weighting them with their model weights (AIC, BIC, or DT) (fig. 2). Indeed, this option is only available when the tree topology has been optimized for every model. The consensus tree is constructed using the Consense program from the PHYLIP package (Felsenstein 2005Go).


Figure 2
View larger version (6K):
[in this window]
[in a new window]
[Download PowerPoint slide]
  		FIG. 2.— Model-averaged tree of HIV-1 pol sequences. The topology shown is the consensus of 88 ML tree topologies, one for every model, weighted according to the AIC weights. The numbers on the branches represent uncertainty due to model selection. In this case, clades (AJ), (AJC), and (HG) are supported by the best and fourth best AIC models (GTR + G, AIC weight = 0.83; TIM3 + G, AIC weight = 0.01; respectively) and others, but not by the second or third best AIC models (GTR + I + G, AIC weight = 0.15; GTR + I, AIC weight = 0.01; respectively).

 
Software Platform and Availability
jModelTest is written in Java and can be started in any operating system with a Java Runtime Environment (see http://www.java.com). However, jModelTest uses other programs for different tasks, and these have been compiled for Mac OSX, Windows XP, and Linux. The package, including installation instructions, documentation, executables, and example data, is distributed free of charge for academic use from the software section at http://darwin.uvigo.es.


   Conclusions
TOP
 Abstract
 Introduction
 Model Selection with jModelTest
 Model Selection Uncertainty
 Conclusions
Acknowledgements
 References
 
Model selection is an important issue in statistical phylogenetics, around which some questions still remain open (Kelchner and Thomas 2007Go). jModelTest addresses some of these, providing an increased flexibility for the user to explore the data and the role of the substitution model on the estimation of phylogenetic trees.


   Acknowledgements
TOP
 Abstract
 Introduction
 Model Selection with jModelTest
 Model Selection Uncertainty
 Conclusions
 Acknowledgements
References
 
I want to thank a number of users of Modeltest that had made numerous comments and suggestions through the years. Special thanks to Stephane Guindon for his generous help with Phyml and to John Huelsenbeck for suggesting the stochastic calculation of CIs. I want to acknowledge Sudhir Kumar for inviting me to present the latest advances in Modeltest at the 2006 SMBE annual meeting, which finally prompted the completion of jModelTest.


   Footnotes
 
Sudhir Kumar, Associate Editor


   References
TOP
 Abstract
 Introduction
 Model Selection with jModelTest
 Model Selection Uncertainty
 Conclusions
 Acknowledgements
 References
 

    Abdo Z, Minin VN, Joyce P, Sullivan J. Accounting for uncertainty in the tree topology has little effect on the decision-theoretic approach to model selection in phylogeny estimation. Mol Biol Evol (2005) 22:691–703.[Abstract/Free Full Text]

    Akaike H. Information theory and an extension of the maximum likelihood principle. In: Second International Symposium on Information Theory—Petrov BN, Csaki F, eds. (1973) Budapest (Hungary): Akademiai Kiado. 267–281.

    Alfaro ME, Huelsenbeck JP. Comparative performance of Bayesian and AIC-based measures of phylogenetic model uncertainty. Syst Biol (2006) 55:89–96.[CrossRef][Web of Science][Medline]

    Buckley TR. Model misspecification and probabilistic tests of topology: evidence from empirical data sets. Syst Biol (2002) 51:509–523.[CrossRef][Web of Science][Medline]

    Buckley TR, Cunningham CW. The effects of nucleotide substitution model assumptions on estimates of nonparametric bootstrap support. Mol Biol Evol (2002) 19:394–405.[Abstract/Free Full Text]

    Burnham KP, Anderson DR. Model selection and multimodel inference. a practical information-theoretic approach (2003) New York: Springer.

    Felsenstein J. Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol (1981) 17:368–376.[CrossRef][Web of Science][Medline]

    Felsenstein J. PHYLIP (phylogeny inference package) (2005) Seattle (WA): Department of Genome Sciences. University of Washington.

    Frati F, Simon C, Sullivan J, Swofford DL. Evolution of the mitochondrial cytochrome oxidase II gene in Collembola. J Mol Evol (1997) 44:145–158.[CrossRef][Web of Science][Medline]

    Gascuel O. BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol Biol Evol (1997) 14:685–695.[Abstract]

    Gilbert D. ReadSeq (2007) Bloomington (IN): Indiana University.

    Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol (2003) 52:696–704.[Abstract/Free Full Text]

    Hasegawa M, Kishino K, Yano T. Dating the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol (1985) 22:160–174.[CrossRef][Web of Science][Medline]

    Huelsenbeck JP, Crandall KA. Phylogeny estimation and hypothesis testing using maximum likelihood. Annu Rev Ecol Syst (1997) 28:437–466.[CrossRef][Web of Science]

    Hurvich CM, Tsai C-L. Regression and time series model selection in small samples. Biometrika (1989) 76:297–307.[Abstract/Free Full Text]

    Johnson JB, Omland KS. Model selection in ecology and evolution. Trends Ecol Evol (2003) 19:101–108.[CrossRef]

    Jukes TH, Cantor CR. Evolution of protein molecules. In: Mammalian protein metabolism—Munro HM, ed. (1969) New York: Academic Press. 21–132.

    Keane TM, Creevey CJ, Pentony MM, Naughton TJ, McLnerney JO. Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol Biol (2006) 6:29.[CrossRef][Medline]

    Kelchner SA, Thomas MA. Model use in phylogenetics: nine key questions. Trends Ecol Evol (2007) 22:87–94.[CrossRef][Medline]

    Kimura M. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J Mol Evol (1980) 16:111–120.[CrossRef][Web of Science][Medline]

    Kimura M. Estimation of evolutionary distances between homologous nucleotide sequences. Proc Natl Acad Sci USA (1981) 78:454–458.[Abstract/Free Full Text]

    Lemmon AR, Moriarty EC. The importance of proper model assumption in Bayesian phylogenetics. Syst Biol (2004) 53:265–277.[Abstract/Free Full Text]

    Minin V, Abdo Z, Joyce P, Sullivan J. Performance-based selection of likelihood models for phylogeny estimation. Syst Biol (2003) 52:674–683.[Abstract/Free Full Text]

    Nylander JA. MrAIC [Internet] (2004) [cited 2008 April 23]. Available from: http://www.abc.se/~nylander/. program distributed by the author.

    Pol D. Empirical problems of the hierarchical likelihood ratio test for model selection. Syst Biol (2004) 53:949–962.[Abstract/Free Full Text]

    Posada D. The effect of branch length variation on the selection of models of molecular evolution. J Mol Evol (2001) 52:434–444.[Web of Science][Medline]

    Posada D. Using Modeltest and PAUP* to select a model of nucleotide substitution. In: Current Protocols in Bioinformatics—Baxevanis AD, Davison DB, Page RDM, Petsko GA, Stein LD, Stormo GD, eds. (2003) New York: John Wiley & Sons. 6.5.1–6.5.14.

    Posada D, Buckley TR. Model selection and model averaging in phylogenetics: advantages of Akaike information criterion and Bayesian approaches over likelihood ratio tests. Syst Biol (2004) 53:793–808.[Abstract/Free Full Text]

    Posada D, Crandall KA. Modeltest: testing the model of DNA substitution. Bioinformatics (1998) 14:817–818.[Abstract/Free Full Text]

    Posada D, Crandall KA. Selecting the best-fit model of nucleotide substitution. Syst Biol (2001) 50:580–601.[Abstract/Free Full Text]

    Schwarz G. Estimating the dimension of a model. Ann Stat (1978) 6:461–464.[CrossRef]

    Sugiura N. Further analysis of the data by Akaike's information criterion and the finite corrections. Commun Stat Theory Methods (1978) A7:13–26.[CrossRef]

    Sullivan J, Joyce P. Model selection in phylogenetics. Annu Rev Ecol Evol Syst (2005) 36:445–466.[CrossRef]

    Sullivan J, Markert JA, Kilpatrick CW. Phylogeography and molecular systematics of the Peromyscus aztecus species group (Rodentia: Muridae) inferred using parsimony and likelihood. Syst Biol (1997) 46:426–440.[CrossRef][Web of Science][Medline]

    Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol (1993) 10:512–526.[Abstract]

    Tavaré S. Some probabilistic and statistical problems in the analysis of DNA sequences. In: Some mathematical questions in biology—DNA sequence analysis—Miura RM, ed. (1986) Providence (RI): American Mathematical Society. 57–86.

    Zharkikh A. Estimation of evolutionary distances between nucleotide sequences. J Mol Evol (1994) 39:315–329.[CrossRef][Web of Science][Medline]

Accepted for publication April 2, 2008.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
 Proc R Soc BHome page
J. B. Searle, P. Kotlik, R. V. Rambau, S. Markova, J. S. Herman, and A. D. McDevitt
The Celtic fringe of Britain: insights from small mammal phylogeography
Proc R Soc B, December 22, 2009; 276(1677): 4287 - 4294.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
H. Ikeda, N. Fujii, and H. Setoguchi
Molecular Evolution of Phytochromes in Cardamine nipponica (Brassicaceae) Suggests the Involvement of PHYE in Local Adaptation
Genetics, June 1, 2009; 182(2): 603 - 614.
[Abstract] [Full Text] [PDF]

Home page
 Am J Trop Med HygHome page
J. A. Farfan-Ale, M. A. Lorono-Pino, J. E. Garcia-Rejon, E. Hovav, A. M. Powers, M. Lin, K. S. Dorman, K. B. Platt, L. C. Bartholomay, V. Soto, et al.
Detection of RNA from a Novel West Nile-like Virus and High Prevalence of an Insect-specific Flavivirus in Mosquitoes in the Yucatan Peninsula of Mexico
Am J Trop Med Hyg, January 1, 2009; 80(1): 85 - 95.
[Abstract] [Full Text] [PDF]

Home page
 Syst BiolHome page
A. Dornburg, F. Santini, and M. E. Alfaro
The Influence of Model Averaging on Clade Posteriors: An Example Using the Triggerfishes (Family Balistidae)
Syst Biol, December 1, 2008; 57(6): 905 - 919.
[Abstract] [Full Text] [PDF]

Home page
 BioinformaticsHome page
S. Ramirez-Flandes and O. Ulloa
Bosque: integrated phylogenetic analysis software
Bioinformatics, November 1, 2008; 24(21): 2539 - 2541.
[Abstract] [Full Text] [PDF]

This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
25/7/1253    most recent
msn083v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Posada, D.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Posada, D.
Social Bookmarking
Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?