A Continuous Method for Gene Flow

Michal Palczewski

doi:10.1534/genetics.113.150904

Return

A Continuous Method for Gene Flow

2013-07-01

Michal Palczewski *

Peter Beerli

DOI：10.1534/genetics.113.150904

Original

Original request for help

Save

Abstract

En 中文

Most modern population genetics inference methods are based on the coalescence framework. Methods that allow estimating parameters of structured populations commonly insert migration events into the genealogies. For these methods the calculation of the coalescence probability density of a genealogy requires a product over all time periods between events. Data sets that contain populations with high rates of gene flow among them require an enormous number of calculations. A new method, transition probability-structured coalescence (TPSC), replaces the discrete migration events with probability statements. Because the speed of calculation is independent of the amount of gene flow, this method allows calculating the coalescence densities efficiently. The current implementation of TPSC uses an approximation simplifying the interaction among lineages. Simulations and coverage comparisons of TPSC vs. MIGRATE show that TPSC allows estimation of high migration rates more precisely, but because of the approximation the estimation of low migration rates is biased. The implementation of TPSC into programs that calculate quantities on phylogenetic tree structures is straightforward, so the TPSC approach will facilitate more general inferences in many computer programs.

Keywords:

EFFECTIVE POPULATION-SIZE

MAXIMUM-LIKELIHOOD

MIGRATION RATES

COALESCENT

PARAMETERS

INFERENCE

SUBPOPULATION

EVOLUTION

MODELS

NUMBER