Estimation of 2N(e)s from temporal allele frequency data

Research output: Contribution to journalJournal articleResearchpeer-review

We develop a new method for estimating effective population sizes, Ne, and selection coefficients, s, from time-series data of allele frequencies sampled from a single diallelic locus. The method is based on calculating transition probabilities, using a numerical solution of the diffusion process, and assuming independent binomial sampling from this diffusion process at each time point. We apply the method in two example applications. First, we estimate selection coefficients acting on the CCR5-{Delta}32 mutation on the basis of published samples of contemporary and ancient human DNA. We show that the data are compatible with the assumption of s = 0, although moderate amounts of selection acting on this mutation cannot be excluded. In our second example, we estimate the selection coefficient acting on a mutation segregating in an experimental phage population. We show that the selection coefficient acting on this mutation is ~0.43.
Original languageEnglish
JournalGenetics
Volume179
Pages (from-to)497-502
Number of pages6
ISSN0016-6731
DOIs
Publication statusPublished - 2008

ID: 9855393