A new method to uncover signatures of divergent and stabilizing selection in quantitative traits

While it is well understood that the pace of evolution depends on the interplay between natural selection, random genetic drift, mutation, and gene flow, it is not always easy to disentangle the relative roles of these factors with data from natural populations. One popular approach to infer whether the observed degree of population differentiation has been influenced by local adaptation is the comparison of neutral marker gene differentiation (as reflected in FST) and quantitative trait divergence (as reflected in QST). However, this method may lead to compromised statistical power, because FST and QST are summary statistics which neglect information on specific pairs of populations, and because current multivariate tests of neutrality involve an averaging procedure over the traits. Further, most FST-QST comparisons actually replace QST by its expectation over the evolutionary process and are thus theoretically flawed. To overcome these caveats, we derived the statistical distribution of population means generated by random genetic drift and used the probability density of this distribution to test whether the observed pattern could be generated by drift alone. We show that our method can differentiate between genetic drift and selection as a cause of population differentiation even in cases with FST=QST and demonstrate with simulated data that it disentangles drift from selection more accurately than conventional FST-QST tests especially when data sets are small.