Bioinformatics Advance Access originally published online on April 15, 2009
Bioinformatics 2009 25(11):1370-1376; doi:10.1093/bioinformatics/btp244
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Many-core algorithms for statistical phylogenetics
1Department of Biomathematics, 2Department of Biostatistics, 3Department of Human Genetics, University of California, Los Angeles, CA 90095, USA and 4Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3JT, UK
*To whom correspondence should be addressed.
| Abstract |
|---|
Motivation: Statistical phylogenetics is computationally intensive, resulting in considerable attention meted on techniques for parallelization. Codon-based models allow for independent rates of synonymous and replacement substitutions and have the potential to more adequately model the process of protein-coding sequence evolution with a resulting increase in phylogenetic accuracy. Unfortunately, due to the high number of codon states, computational burden has largely thwarted phylogenetic reconstruction under codon models, particularly at the genomic-scale. Here, we describe novel algorithms and methods for evaluating phylogenies under arbitrary molecular evolutionary models on graphics processing units (GPUs), making use of the large number of processing cores to efficiently parallelize calculations even for large state-size models.
Results: We implement the approach in an existing Bayesian framework and apply the algorithms to estimating the phylogeny of 62 complete mitochondrial genomes of carnivores under a 60-state codon model. We see a near 90-fold speed increase over an optimized CPU-based computation and a >140-fold increase over the currently available implementation, making this the first practical use of codon models for phylogenetic inference over whole mitochondrial or microorganism genomes.
Availability and implementation: Source code provided in BEAGLE: Broad-platform Evolutionary Analysis General Likelihood Evaluator, a cross-platform/processor library for phylogenetic likelihood computation (http://beagle-lib.googlecode.com/). We employ a BEAGLE-implementation using the Bayesian phylogenetics framework BEAST (http://beast.bio.ed.ac.uk/).
Contact: msuchard{at}ucla.edu; a.rambaut{at}ed.ac.uk
Associate Editor: Martin Bishop
Received on March 16, 2009; revised on April 2, 2009; accepted on April 4, 2009
This article has been cited by other articles:
![]() |
E. W. Bloomquist and M. A. Suchard Unifying Vertical and Nonvertical Evolution: A Stochastic ARG-based Framework Syst Biol, November 9, 2009; (2009) syp076v1. [Abstract] [Full Text] [PDF] |
||||
