Inferring species trees directly from biallelic genetic markers: bypassing gene trees in a full coalescent analysis.

ABSTRACT: The multispecies coalescent provides an elegant theoretical framework for estimating species trees and species demographics from genetic markers. However, practical applications of the multispecies coalescent model are limited by the need to integrate or sample over all gene trees possible for each genetic marker. Here we describe a polynomial-time algorithm that computes the likelihood of a species tree directly from the markers under a finite-sites model of mutation effectively integrating over all possible gene trees. The method applies to independent (unlinked) biallelic markers such as well-spaced single nucleotide polymorphisms, and we have implemented it in SNAPP, a Markov chain Monte Carlo sampler for inferring species trees, divergence dates, and population sizes. We report result

SUBMITTER: Bryant D 

PROVIDER: S-EPMC3408069 | biostudies-literature | 2012 Aug

REPOSITORIES: biostudies-literature