Genetics, Vol. 167, 977-988, June 2004, Copyright © 2004
doi:10.1534/genetics.103.026146

Comparative Evaluation of a New Effective Population Size Estimator Based on Approximate Bayesian Computation

* Laboratoire d'Ecologie Alpine, UMR Centre National de la Recherche Scientifique 5553, Université Joseph Fourier, F38041 BP 53 Cedex 09, Grenoble, France
{dagger} School of Animal and Microbial Sciences, University of Reading, Whiteknights, Reading RG6 6AJ, United Kingdom

1 Corresponding author: LECA-Génomique des Population et Biodiversité, UMR CNRS 5553, Université Joseph Fourier, F-38041 BP 53 Cedex 09, Grenoble, France.
E-mail: dtallmon42{at}yahoo.com

We describe and evaluate a new estimator of the effective population size (Ne), a critical parameter in evolutionary and conservation biology. This new "SummStat" Ne estimator is based upon the use of summary statistics in an approximate Bayesian computation framework to infer Ne. Simulations of a Wright-Fisher population with known Ne show that the SummStat estimator is useful across a realistic range of individuals and loci sampled, generations between samples, and Ne values. We also address the paucity of information about the relative performance of Ne estimators by comparing the SummStat estimator to two recently developed likelihood-based estimators and a traditional moment-based estimator. The SummStat estimator is the least biased of the four estimators compared. In 32 of 36 parameter combinations investigated using initial allele frequencies drawn from a Dirichlet distribution, it has the lowest bias. The relative mean square error (RMSE) of the SummStat estimator was generally intermediate to the others. All of the estimators had RMSE > 1 when small samples (n = 20, five loci) were collected a generation apart. In contrast, when samples were separated by three or more generations and Ne ≤ 50, the SummStat and likelihood-based estimators all had greatly reduced RMSE. Under the conditions simulated, SummStat confidence intervals were more conservative than the likelihood-based estimators and more likely to include true Ne. The greatest strength of the SummStat estimator is its flexible structure. This flexibility allows it to incorporate any potentially informative summary statistic from population genetic data.




This article has been cited by other articles:


Home page
GeneticsHome page
V. C. Sousa, M. Fritz, M. A. Beaumont, and L. Chikhi
Approximate Bayesian Computation Without Summary Statistics: The Case of Admixture
Genetics, April 1, 2009; 181(4): 1507 - 1519.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
P. E. Jorde and N. Ryman
Unbiased Estimator for Genetic Drift and Effective Population Size
Genetics, October 1, 2007; 177(2): 927 - 935.
[Abstract] [Full Text] [PDF]


Home page
BioinformaticsHome page
T. Antao, A. Beja-Pereira, and G. Luikart
MODELER4SIMCOAL2: A user-friendly, extensible modeler of demography and linked loci for coalescent simulations
Bioinformatics, July 15, 2007; 23(14): 1848 - 1850.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
R. S. Waples and M. Yokota
Temporal Estimates of Effective Population Size in Species With Overlapping Generations
Genetics, January 1, 2007; 175(1): 219 - 233.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
M. M. Tanaka, A. R. Francis, F. Luciani, and S. A. Sisson
Using Approximate Bayesian Computation to Estimate Tuberculosis Transmission Parameters From Genotype Data
Genetics, July 1, 2006; 173(3): 1511 - 1520.
[Abstract] [Full Text] [PDF]


Home page
Am J Trop Med HygHome page
I. MORLAIS, R. GIROD, R. HUNT, F. SIMARD, and D. FONTENILLE
POPULATION STRUCTURE OF ANOPHELES ARABIENSIS ON LA REUNION ISLAND, INDIAN OCEAN
Am J Trop Med Hyg, December 1, 2005; 73(6): 1077 - 1082.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
S. C. Leman, Y. Chen, J. E. Stajich, M. A. F. Noor, and M. K. Uyenoyama
Likelihoods From Summary Statistics: Recent Divergence Between Species
Genetics, November 1, 2005; 171(3): 1419 - 1436.
[Abstract] [Full Text] [PDF]


Home page
Phil Trans R Soc BHome page
J. Wang
Estimation of effective population sizes from data on genetic markers
Phil Trans R Soc B, July 29, 2005; 360(1459): 1395 - 1409.
[Abstract] [Full Text] [PDF]


Home page
GeneticsHome page
E. C. Anderson
An Efficient Monte Carlo Method for Estimating Ne From Temporally Spaced Samples Using a Coalescent-Based Likelihood
Genetics, June 1, 2005; 170(2): 955 - 967.
[Abstract] [Full Text] [PDF]