Cargando…

Bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation

BACKGROUND: Bayesian phylogenetic inference holds promise as an alternative to maximum likelihood, particularly for large molecular-sequence data sets. We have investigated the performance of Bayesian inference with empirical and simulated protein-sequence data under conditions of relative branch-le...

Descripción completa

Detalles Bibliográficos
Autores principales: Mar, Jessica C, Harlow, Timothy J, Ragan, Mark A
Formato: Texto
Lenguaje:English
Publicado: BioMed Central 2005
Materias:
Acceso en línea:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC549035/
https://www.ncbi.nlm.nih.gov/pubmed/15676079
http://dx.doi.org/10.1186/1471-2148-5-8
_version_ 1782122379456020480
author Mar, Jessica C
Harlow, Timothy J
Ragan, Mark A
author_facet Mar, Jessica C
Harlow, Timothy J
Ragan, Mark A
author_sort Mar, Jessica C
collection PubMed
description BACKGROUND: Bayesian phylogenetic inference holds promise as an alternative to maximum likelihood, particularly for large molecular-sequence data sets. We have investigated the performance of Bayesian inference with empirical and simulated protein-sequence data under conditions of relative branch-length differences and model violation. RESULTS: With empirical protein-sequence data, Bayesian posterior probabilities provide more-generous estimates of subtree reliability than does the nonparametric bootstrap combined with maximum likelihood inference, reaching 100% posterior probability at bootstrap proportions around 80%. With simulated 7-taxon protein-sequence datasets, Bayesian posterior probabilities are somewhat more generous than bootstrap proportions, but do not saturate. Compared with likelihood, Bayesian phylogenetic inference can be as or more robust to relative branch-length differences for datasets of this size, particularly when among-sites rate variation is modeled using a gamma distribution. When the (known) correct model was used to infer trees, Bayesian inference recovered the (known) correct tree in 100% of instances in which one or two branches were up to 20-fold longer than the others. At ratios more extreme than 20-fold, topological accuracy of reconstruction degraded only slowly when only one branch was of relatively greater length, but more rapidly when there were two such branches. Under an incorrect model of sequence change, inaccurate trees were sometimes observed at less extreme branch-length ratios, and (particularly for trees with single long branches) such trees tended to be more inaccurate. The effect of model violation on accuracy of reconstruction for trees with two long branches was more variable, but gamma-corrected Bayesian inference nonetheless yielded more-accurate trees than did either maximum likelihood or uncorrected Bayesian inference across the range of conditions we examined. Assuming an exponential Bayesian prior on branch lengths did not improve, and under certain extreme conditions significantly diminished, performance. The two topology-comparison metrics we employed, edit distance and Robinson-Foulds symmetric distance, yielded different but highly complementary measures of performance. CONCLUSIONS: Our results demonstrate that Bayesian inference can be relatively robust against biologically reasonable levels of relative branch-length differences and model violation, and thus may provide a promising alternative to maximum likelihood for inference of phylogenetic trees from protein-sequence data.
format Text
id pubmed-549035
institution National Center for Biotechnology Information
language English
publishDate 2005
publisher BioMed Central
record_format MEDLINE/PubMed
spelling pubmed-5490352005-02-18 Bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation Mar, Jessica C Harlow, Timothy J Ragan, Mark A BMC Evol Biol Research Article BACKGROUND: Bayesian phylogenetic inference holds promise as an alternative to maximum likelihood, particularly for large molecular-sequence data sets. We have investigated the performance of Bayesian inference with empirical and simulated protein-sequence data under conditions of relative branch-length differences and model violation. RESULTS: With empirical protein-sequence data, Bayesian posterior probabilities provide more-generous estimates of subtree reliability than does the nonparametric bootstrap combined with maximum likelihood inference, reaching 100% posterior probability at bootstrap proportions around 80%. With simulated 7-taxon protein-sequence datasets, Bayesian posterior probabilities are somewhat more generous than bootstrap proportions, but do not saturate. Compared with likelihood, Bayesian phylogenetic inference can be as or more robust to relative branch-length differences for datasets of this size, particularly when among-sites rate variation is modeled using a gamma distribution. When the (known) correct model was used to infer trees, Bayesian inference recovered the (known) correct tree in 100% of instances in which one or two branches were up to 20-fold longer than the others. At ratios more extreme than 20-fold, topological accuracy of reconstruction degraded only slowly when only one branch was of relatively greater length, but more rapidly when there were two such branches. Under an incorrect model of sequence change, inaccurate trees were sometimes observed at less extreme branch-length ratios, and (particularly for trees with single long branches) such trees tended to be more inaccurate. The effect of model violation on accuracy of reconstruction for trees with two long branches was more variable, but gamma-corrected Bayesian inference nonetheless yielded more-accurate trees than did either maximum likelihood or uncorrected Bayesian inference across the range of conditions we examined. Assuming an exponential Bayesian prior on branch lengths did not improve, and under certain extreme conditions significantly diminished, performance. The two topology-comparison metrics we employed, edit distance and Robinson-Foulds symmetric distance, yielded different but highly complementary measures of performance. CONCLUSIONS: Our results demonstrate that Bayesian inference can be relatively robust against biologically reasonable levels of relative branch-length differences and model violation, and thus may provide a promising alternative to maximum likelihood for inference of phylogenetic trees from protein-sequence data. BioMed Central 2005-01-28 /pmc/articles/PMC549035/ /pubmed/15676079 http://dx.doi.org/10.1186/1471-2148-5-8 Text en Copyright © 2005 Mar et al; licensee BioMed Central Ltd. http://creativecommons.org/licenses/by/2.0 This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( (http://creativecommons.org/licenses/by/2.0) ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
spellingShingle Research Article
Mar, Jessica C
Harlow, Timothy J
Ragan, Mark A
Bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation
title Bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation
title_full Bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation
title_fullStr Bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation
title_full_unstemmed Bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation
title_short Bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation
title_sort bayesian and maximum likelihood phylogenetic analyses of protein sequence data under relative branch-length differences and model violation
topic Research Article
url https://www.ncbi.nlm.nih.gov/pmc/articles/PMC549035/
https://www.ncbi.nlm.nih.gov/pubmed/15676079
http://dx.doi.org/10.1186/1471-2148-5-8
work_keys_str_mv AT marjessicac bayesianandmaximumlikelihoodphylogeneticanalysesofproteinsequencedataunderrelativebranchlengthdifferencesandmodelviolation
AT harlowtimothyj bayesianandmaximumlikelihoodphylogeneticanalysesofproteinsequencedataunderrelativebranchlengthdifferencesandmodelviolation
AT raganmarka bayesianandmaximumlikelihoodphylogeneticanalysesofproteinsequencedataunderrelativebranchlengthdifferencesandmodelviolation