Partitioning the effects of isolation by distance, environment, and physical barriers on genomic divergence between parapatric threespine stickleback

Abstract

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation-by-distance (IBD) and isolation-by-environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake-stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake-stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ∼1.9 km of IBD. Remarkably, of our three lake-stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.