  The increasing availability of population-level allele frequency data across
one or more related populations necessitates the development of methods that
can efficiently estimate population genetics parameters, such as the strength
of selection acting on the population(s), from such data. Existing methods for
this problem in the setting of the Wright-Fisher diffusion model are primarily
likelihood-based, and rely on numerical approximation for likelihood
computation and on bootstrapping for assessment of variability in the resulting
estimates, requiring extensive computation. Recent work (Jenkins and Spano,
2015) has provided a method for obtaining exact samples from general
Wright-Fisher diffusion processes, enabling the development of methods for
Bayesian estimation in this setting. We develop and implement a Bayesian method
for estimating the strength of selection based on the Wright-Fisher diffusion
for data sampled at a single time point. The method utilizes the work of
Jenkins and Spano (2015) to develop a Markov chain Monte Carlo algorithm to
draw samples from the joint posterior distribution of the selection coefficient
and the allele frequencies. We demonstrate that when assumptions about the
initial allele frequencies are accurate the method performs well for both
simulated data and for an empirical data set on hypoxia in flies (Zhou et al.
2011), where we find evidence for strong positive selection in a region of
chromosome 2L previously identified by Ronen et al. (2013). We discuss possible
extensions of our method to the more general settings commonly encountered in
practice, highlighting the advantages of Bayesian approaches to inference in
this setting.
