Sampling strategies for frequency spectrum-based population genomic inference

John D. Robinson, Alec J. Coffman, Michael J. Hickerson, Ryan N. Gutenkunst

Research output: Contribution to journalArticlepeer-review

55 Scopus citations


Background: The allele frequency spectrum (AFS) consists of counts of the number of single nucleotide polymorphism (SNP) loci with derived variants present at each given frequency in a sample. Multiple approaches have recently been developed for parameter estimation and calculation of model likelihoods based on the joint AFS from two or more populations. We conducted a simulation study of one of these approaches, implemented in the Python module δaδi, to compare parameter estimation and model selection accuracy given different sample sizes under one- and two-population models. Results: Our simulations included a variety of demographic models and two parameterizations that differed in the timing of events (divergence or size change). Using a number of SNPs reasonably obtained through next-generation sequencing approaches (10,000 - 50,000), accurate parameter estimates and model selection were possible for models with more ancient demographic events, even given relatively small numbers of sampled individuals. However, for recent events, larger numbers of individuals were required to achieve accuracy and precision in parameter estimates similar to that seen for models with older divergence or population size changes. We quantify i) the uncertainty in model selection, using tools from information theory, and ii) the accuracy and precision of parameter estimates, using the root mean squared error, as a function of the timing of demographic events, sample sizes used in the analysis, and complexity of the simulated models. Conclusions: Here, we illustrate the utility of the genome-wide AFS for estimating demographic history and provide recommendations to guide sampling in population genomics studies that seek to draw inference from the AFS. Our results indicate that larger samples of individuals (and thus larger AFS) provide greater power for model selection and parameter estimation for more recent demographic events.

Original languageEnglish (US)
Article number254
JournalBMC Evolutionary Biology
Issue number1
StatePublished - Dec 4 2014


  • Allele frequency spectrum
  • Demographic history
  • Model selection
  • Parameter uncertainty
  • SNP

ASJC Scopus subject areas

  • Ecology, Evolution, Behavior and Systematics


Dive into the research topics of 'Sampling strategies for frequency spectrum-based population genomic inference'. Together they form a unique fingerprint.

Cite this