The evolutionary consequences of temporal variation in selection remain hotly debated. We explored these consequences by studying threespine stickleback in a set of bar-built estuaries along the central California coast. In most years, heavy rains induce water flow strong enough to break through isolating sand bars, connecting streams to the ocean. New sand bars typically re-form within a few weeks or months, thereby re-isolating populations within the estuaries. These breaching events cause severe and often extremely rapid changes in abiotic and biotic conditions, including shifts in predator abundance. We investigated whether this strong temporal environmental variation can maintain within-population variation while eroding adaptive divergence among populations that would be caused by spatial variation in selection. We used neutral genetic markers to explore population structure, and then analyzed how stickleback armor traits, the associated genes Eda and Pitx1, and elemental composition (%P) varies within and among populations. Despite strong gene flow, we detected evidence for divergence in stickleback defensive traits and Eda genotypes associated with predation regime. However, this among-population variation was lower than that observed among other stickleback populations exposed to divergent predator regimes. In addition, within-population variation was very high as compared to populations from environmentally stable locations. Elemental composition was strongly associated with armor traits, Eda genotype, and the presence of predators; thus suggesting that spatiotemporal variation in armor traits generates corresponding variation in elemental phenotypes. We conclude that gene flow, and especially temporal environmental variation, can maintain high levels of within-population variation while reducing, but not eliminating, among-population variation driven by spatial environmental variation.

JEB-2017-00599_phenoPhenotypic dataJEB-2017-00599_NeutralNeutral markersJEB-2017-00599_EDA-PITX1EDA-Pitx1 markersJEB-2017-00599_Morpho2D morphometrics data tablesJEB-2017-00599_stoichStoichiometry data