Quantifying population structure on short timescales

Quantifying the contribution of the various processes that influence population genetic structure is important, but difficult. One of the reasons is that no single measure appropriately quantifies all aspects of genetic structure. An increasing number of studies is analysing population structure using the statistic D, which measures genetic differentiation, next to G(ST) , which quantifies the standardized variance in allele frequencies among populations. Few studies have evaluated which statistic is most appropriate in particular situations. In this study, we evaluated which index is more suitable in quantifying postglacial divergence between three-spined stickleback (Gasterosteus aculeatus) populations from Western Europe. Population structure on this short timescale (10 000 generations) is probably shaped by colonization history, followed by migration and drift. Using microsatellite markers and anticipating that D and G(ST) might have different capacities to reveal these processes, we evaluated population structure at two levels: (i) between lowland and upland populations, aiming to infer historical processes; and (ii) among upland populations, aiming to quantify contemporary processes. In the first case, only D revealed clear clusters of populations, putatively indicative of population ancestry. In the second case, only G(ST) was indicative for the balance between migration and drift. Simulations of colonization and subsequent divergence in a hierarchical stepping stone model confirmed this discrepancy, which becomes particularly strong for markers with moderate to high mutation rates. We conclude that on short timescales, and across strong clines in population size and connectivity, D is useful to infer colonization history, whereas G(ST) is sensitive to more recent demographic events.