Microgeographical diversification of threespine stickleback: body shape–habitat correlations in a small, ecologically diverse Alaskan drainage

Adaptive radiations are a major source of evolutionary diversity in nature, and understanding how they originate and how organisms diversify during the early stages of adaptive radiation is a major problem in evolutionary biology. The relationship between habitat type and body shape variation was investigated in a postglacial radiation of threespine stickleback in the upper Fish Creek drainage of Cook Inlet, Alaska. Although small, the upper Fish Creek drainage includes ecologically diverse lakes and streams in close proximity to one another that harbour abundant stickleback. Specimens from ancestral anadromous and derived resident freshwater populations differed substantially and could be distinguished by body shape alone, suggesting that the initial stages of adaptation contribute disproportionately to evolutionary divergence. Body shape divergence among resident freshwater populations was also considerable, and phenotypic distances among samples from freshwater populations were associated with habitat type but not geographical distance. As expected, stream stickleback from slow-moving, structurally complex environments tended to have the deepest bodies, stickleback from lakes with a mostly benthic habitat were similar but less extreme, and stickleback from lakes with a mostly limnetic habitat were the most shallow-bodied, elongate fish. Beyond adapting rapidly to conditions in freshwater environments, stickleback can diversify rapidly over small geographical scales in freshwater systems despite opportunities for gene flow. This study highlights the importance of ecological heterogeneity over small geographical scales for evolutionary diversification during the early stages of adaptive radiation, and lays the foundation for future research on this ecologically diverse, postglacial system.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 98 , 139–151.     

Natural selection drives patterns of lake-stream divergence in stickleback foraging morphology

To what extent are patterns of biological diversification determined by natural selection? We addressed this question by exploring divergence in foraging morphology of threespine stickleback fish inhabiting lake and stream habitats within eight independent watersheds. We found that lake fish generally displayed more developed gill structures and had more streamlined bodies than did stream fish. Diet analysis revealed that these morphological differences were associated with limnetic vs. benthic foraging modes, and that the extent of morphological divergence within watersheds reflected differences in prey resources utilized by lake and stream fish. We also found that patterns of divergence were unrelated to patterns of phenotypic trait (co)variance within populations (i.e. the ‘line of least resistance’). Instead, phenotypic (co)variances were more likely to have been shaped by adaptation to lake vs. stream habitats. Our study thus implicates natural selection as a strong deterministic force driving morphological diversification in lake–stream stickleback. The strength of this inference was obtained by complementing a standard analysis of parallel divergence in means between discrete habitat categories (lake vs. stream) with quantitative estimates of selective forces and information on trait (co)variances.     

Does plasticity enhance or dampen phenotypic parallelism? A test with three lake–stream stickleback pairs

Parallel (and convergent) phenotypic variation is most often studied in the wild, where it is difficult to disentangle genetic vs. environmentally induced effects. As a result, the potential contributions of phenotypic plasticity to parallelism (and nonparallelism) are rarely evaluated in a formal sense. Phenotypic parallelism could be enhanced by plasticity that causes stronger parallelism across populations in the wild than would be expected from genetic differences alone. Phenotypic parallelism could be dampened if site‐specific plasticity induced differences between otherwise genetically parallel populations. We used a common‐garden study of three independent lake–stream stickleback population pairs to evaluate the extent to which adaptive divergence has a genetic or plastic basis, and to investigate the enhancing vs. dampening effects of plasticity on phenotypic parallelism. We found that lake–stream differences in most traits had a genetic basis, but that several traits also showed contributions from plasticity. Moreover, plasticity was much more prevalent in one watershed than in the other two. In most cases, plasticity enhanced phenotypic parallelism, whereas in a few cases, plasticity had a dampening effect. Genetic and plastic contributions to divergence seem to play a complimentary, likely adaptive, role in phenotypic parallelism of lake–stream stickleback. These findings highlight the value of formally comparing wild‐caught and laboratory‐reared individuals in the study of phenotypic parallelism.     

Molecular investigation of genetic assimilation during the rapid adaptive radiations of East African cichlid fishes

Adaptive radiations are characterized by adaptive diversification intertwined with rapid speciation within a lineage resulting in many ecologically specialized, phenotypically diverse species. It has been proposed that adaptive radiations can originate from ancestral lineages with pronounced phenotypic plasticity in adaptive traits, facilitating ecologically driven phenotypic diversification that is ultimately fixed through genetic assimilation of gene regulatory regions. This study aimed to investigate how phenotypic plasticity is reflected in gene expression patterns in the trophic apparatus of several lineages of East African cichlid fishes, and whether the observed patterns support genetic assimilation. This investigation used a split brood experimental design to compare adaptive plasticity in species from within and outside of adaptive radiations. The plastic response was induced in the crushing pharyngeal jaws through feeding individuals either a hard or soft diet. We find that nonradiating, basal lineages show higher levels of adaptive morphological plasticity than the derived, radiated lineages, suggesting that these differences have become partially genetically fixed during the formation of the adaptive radiations. Two candidate genes that may have undergone genetic assimilation, gif and alas1, were identified, in addition to alterations in the wiring of LPJ patterning networks. Taken together, our results suggest that genetic assimilation may have dampened the inducibility of plasticity related genes during the adaptive radiations of East African cichlids, flattening the reaction norms and canalizing their feeding phenotypes, driving adaptation to progressively more narrow ecological niches.     

Ancestral Plasticity and Allometry in Threespine Stickleback Fish Reveal Phenotypes Associated with Derived, Freshwater Ecotypes

For over a century, evolutionary biologists have debated whether and how phenotypic plasticity impacts the processes of adaptation and diversification. The empirical tests required to resolve these issues have proven elusive, mainly because it requires documentation of ancestral reaction norms, a difficult prospect where many ancestors are either extinct or have evolved. The threespine stickleback radiation is not limited in this regard, making it an ideal system in which to address general questions regarding the role of plasticity in adaptive evolution. As retreating ice sheets have exposed new habitats, oceanic stickleback founded innumerable freshwater populations, many of which have evolved parallel adaptations to their new environments. Because the founding oceanic population is extant, we can directly evaluate whether specific patterns of ancestral phenotypic expression in the context of novel environments (plasticity), or over ontogeny, predisposed the repeated evolution of "benthic" and "limnetic" ecotypes in shallow and deep lakes, respectively. Consistent with this hypothesis, we found that oceanic stickleback raised in a complex habitat and fed a macroinvertebrate diet expressed traits resembling derived, benthic fish. Alternatively, when reared in a simple environment on a diet of zooplankton, oceanic stickleback developed phenotypes resembling derived, limnetic fish. As fish in both treatments grew, their body depths increased allometrically, as did the size of their mouths, while their eyes became relatively smaller. Allometric trajectories were subtly but significantly impacted by rearing environment. Thus, both environmental and allometric influences on development, along with their interactive effects, produced variation in phenotypes consistent with derived benthic and limnetic fish, which may have predisposed the repeated genetic accommodation of this specific suite of traits. We also found significant shape differences between marine and anadromous stickleback, which has implications for evaluating the ancestral state of stickleback traits.     

Ontogenetic and evolutionary effects of predation and competition on nine-spined stickleback (Pungitius pungitius) body size

1.?Individual- and population-level variation in body size and growth often correlates with many fitness traits. Predation and food availability are expected to affect body size and growth as important agents of both natural selection and phenotypic plasticity. How differences in predation and food availability affect body size/growth during ontogeny in populations adapted to different predation and competition regimes is rarely studied. 2.?Nine-spined stickleback (Pungitius pungitius) populations originating from habitats with varying levels of predation and competition are known to be locally adapted to their respective habitats in terms of body size and growth. Here, we studied how different levels of perceived predation risk and competition during ontogeny affect the reaction norms of body size and growth in (i) marine and pond populations adapted to different levels of predation and competition and (ii) different sexes. We reared nine-spined stickleback in a factorial experiment under two levels of perceived predation risk (present/absent) and competition (high/low food supply). 3.?We found divergence in the reaction norms at two levels: (i) predation-adapted marine stickleback had stronger reactions to predatory cues than intraspecific competition-adapted pond stickleback, the latter being more sensitive to available food than the marine fish and (ii) females reacting more strongly to the treatments than males. 4.?The repeated, habitat-dependent nature of the differences suggests that natural selection is the agent behind the observed patterns. Our results suggest that genetic adaptation to certain environmental factors also involves an increase in the range of expressible phenotypic plasticity. We found support for this phenomenon at two levels: (i) across populations driven by habitat type and (ii) within populations driven by sex.     

DISENTANGLING THE ROLE OF PHENOTYPIC PLASTICITY AND GENETIC DIVERGENCE IN CONTEMPORARY ECOTYPE FORMATION DURING A BIOLOGICAL INVASION

The occurrence of contemporary ecotype formation through adaptive divergence of populations within the range of an invasive species typically requires standing genetic variation but can be facilitated by phenotypic plasticity. The relative contributions of both of these to adaptive trait differentiation have rarely been simultaneously quantified in recently diverging vertebrate populations. Here we study a case of intraspecific divergence into distinct lake and stream ecotypes of threespine stickleback that evolved in the past 140 years within the invasive range in Switzerland. Using a controlled laboratory experiment with full‐sib crosses and treatments mimicking a key feature of ecotypic niche divergence, we test if the phenotypic divergence that we observe in the wild results from phenotypic plasticity or divergent genetic predisposition. Our experimental groups show qualitatively similar phenotypic divergence as those observed among wild adults. The relative contribution of plasticity and divergent genetic predisposition differs among the traits studied, with traits related to the biomechanics of feeding showing a stronger genetic predisposition, whereas traits related to locomotion are mainly plastic. These results implicate that phenotypic plasticity and standing genetic variation interacted during contemporary ecotype formation in this case.     

Differences in rheotactic responses contribute to divergent habitat use between parapatric lake and stream threespine stickleback

Migration among populations is widely thought to undermine adaptive divergence, assuming gene flow arises from random movement of individuals. If individuals instead differ in dispersal behavior, phenotype‐dependent dispersal can reduce the effective rate of gene flow or even facilitate divergence. For example, parapatric populations of lake and stream stickleback tend to actively avoid dispersing into the adjoining habitat. However, the behavioral basis of this nonrandom dispersal was previously unknown. Here, we show that lake and stream stickleback exhibit divergent rheotactic responses (behavioral response to currents). During the breeding season, wild‐caught inlet stream stickleback were better than lake fish at maintaining position in currents, faced upstream more, and spent more time in low‐current areas. As a result, stream fish expended significantly less energy in currents than did lake fish. These divergent rheotactic responses likely contribute to divergent habitat use by lake and stream stickleback. Although rheotactic differences were absent in nonbreeding fish, divergent behavior of breeding‐season fish may suffice for assortative mating by breeding location. The resulting reproductive isolation between lake and stream fish may explain the fine‐scale evolutionary differentiation in parapatric stickleback populations.     

Morphological anti‐predator defences in the nine‐spined stickleback: constitutive,induced or both?

Environmental differences among populations are expected to lead to local adaptation, while spatial or temporal environmental variation within a population will favour evolution of phenotypic plasticity. As plasticity itself can be under selection, locally adapted populations can vary in levels of plasticity. Nine‐spined stickleback (Pungitius pungitius) originating from isolated ponds (low piscine predation risk, high competition) vs. lake and marine populations (high piscine predation risk, low competition) are known to be morphologically adapted to their respective environments. However, nothing is known about their ability to express phenotypic plasticity in morphology in response to perceived predation risk or food availability/competition. We studied predator‐induced phenotypic plasticity in body shape and armour of marine and pond nine‐spined stickleback in a factorial common garden experiment with two predator treatments (present vs. absent) and two feeding regimes (low vs. high). The predation treatment did not induce any morphological shifts in fish from either habitat or food regime. However, strong habitat‐dependent differences between populations as well as strong sexual dimorphism in both body shape and armour were found. The lack of predator‐induced plasticity in development of the defence traits (viz. body armour and body depth) suggests that morphological anti‐predator traits in nine‐spined stickleback are strictly constitutive, rather than inducible. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ??, ??–??.     

Adaptive divergence and the balance between selection and gene flow: lake and stream stickleback in the Misty system

We investigated the interplay between natural selection and gene flow in the adaptive divergence of threespine stickleback (Gasterosteus aculeatus) that reside parapatrically in lakes and streams. Within the Misty Lake system (Vancouver Island, British Columbia), stickleback from the inlet stream (flowing into the lake) have fewer gill rakers and deeper bodies than stickleback from the lake--differences thought to facilitate foraging (benthic macroinvertebrates in the stream vs. zooplankton in the open water of the lake). Common-garden experiments demonstrated that these differences have a genetic basis. Reciprocal transplant enclosure experiments showed that lake and inlet stickleback grow best in their home environments (although differences were subtle and often not significant). Release-recapture experiments in the inlet showed that lake fish are less well-suited than inlet fish for life in the stream (higher mortality or emigration in lake fish). Morphological divergence in the wild and under common rearing was greater between the lake and the inlet than between the lake and the outlet. Genetic divergence (mitochondrial DNA and microsatellites) was greatest between the lake and the upper inlet (1.8 km upstream from the lake), intermediate between the lake and the lower inlet (0.9 km upstream), and least between the lake and the outlet stream (1.2 km downstream). Relative levels of gene flow estimated from genetic data showed the inverse pattern. The negative association between morphological divergence and gene flow is consistent with the expectation that gene flow can constrain adaptation. Estimated absolute levels of gene flow also implied a constraint on adaptation in the outlet but not the inlet. Our results suggest that natural selection promotes the adaptive divergence of lake and stream stickleback. but that the magnitude of divergence can be constrained by gene flow.     

Strong and parallel salinity‐induced phenotypic plasticity in one generation of threespine stickleback

Phenotypic plasticity is a major factor contributing to variation of organisms in nature, yet its evolutionary significance is insufficiently understood. One example system where plasticity might have played an important role in an adaptive radiation is the threespine stickleback (Gasterosteus aculeatus), a fish that has diversified after invading freshwater lakes repeatedly from the marine habitat. The parallel phenotypic changes that occurred in this radiation were extremely rapid. This study evaluates phenotypic plasticity in stickleback body shape in response to salinity in fish stemming from a wild freshwater population. Using a split‐clutch design, we detected surprisingly large phenotypically plastic changes in body shape after one generation. Fish raised in salt water developed shallower bodies and longer jaws, and these changes were consistent and parallel across families. Although this work highlights the effect of phenotypic plasticity, we also find indications that constraints may play a role in biasing the direction of possible phenotypic change. The slopes of the allometric relationship of individual linear traits did not change across treatments, indicating that plastic change does not affect the covariation of traits with overall size. We conclude that stickleback have a large capacity for plastic phenotypic change in response to salinity and that plasticity and evolutionary constraints have likely contributed to the phenotypic diversification of these fish.     

Gene expression stasis and plasticity following migration into a foreign environment

Selection against migrants is key to maintaining genetic differences between populations linked by dispersal. However, migrants may mitigate fitness costs by proactively choosing among available habitats, or by phenotypic plasticity. We previously reported that a reciprocal transplant of lake and stream stickleback (Gasterosteus aculeatus) found little support for divergent selection. Here, we revisit that experiment to test whether phenotypic plasticity in gene expression may have helped migrants adjust to unfamiliar habitats. We measured gene expression profiles in stickleback via TagSeq and tested whether migrants between lake and stream habitats exhibited a plastic response to their new environment that allowed them to converge on the expression profile of adapted natives. We report extensive gene expression differences between genetically divergent lake and stream stickleback, despite gene flow . But for many genes, expression was highly plastic. Fish transplanted into the adjoining habitat partially converged on the expression profile typical of natives from their new habitat. This suggests that expression plasticity may soften the impact of migration. Nonetheless, lake and stream fish differed in survival rates and parasite infection rates in our study, implying that expression plasticity is not fast or extensive enough to fully homogenize fish performance.     

Comparative Analysis of Japanese Three-Spined Stickleback Clades Reveals the Pacific Ocean Lineage Has Adapted to Freshwater Environments while the Japan Sea Has Not

Divergent selection and adaptive divergence can increase phenotypic diversification amongst populations and lineages. Yet adaptive divergence between different environments, habitats or niches does not occur in all lineages. For example, the colonization of freshwater environments by ancestral marine species has triggered adaptive radiation and phenotypic diversification in some taxa but not in others. Studying closely related lineages differing in their ability to diversify is an excellent means of understanding the factors promoting and constraining adaptive evolution. A well-known example of the evolution of increased phenotypic diversification following freshwater colonization is the three-spined stickleback. Two closely related stickleback lineages, the Pacific Ocean and the Japan Sea occur in Japan. However, Japanese freshwater stickleback populations are derived from the Pacific Ocean lineage only, suggesting the Japan Sea lineage is unable to colonize freshwater. Using stable isotope data and trophic morphology, we first show higher rates of phenotypic and ecological diversification between marine and freshwater populations within the Pacific Ocean lineage, confirming adaptive divergence has occurred between the two lineages and within the Pacific Ocean lineage but not in the Japan Sea lineage. We further identified consistent divergence in diet and foraging behaviour between marine forms from each lineage, confirming Pacific Ocean marine sticklebacks, from which all Japanese freshwater populations are derived, are better adapted to freshwater environments than Japan Sea sticklebacks. We suggest adaptive divergence between ancestral marine populations may have played a role in constraining phenotypic diversification and adaptive evolution in Japanese sticklebacks.     

Adaptive divergence vs. environmental plasticity: tracing local genetic adaptation of metamorphosis traits in salamanders

In order to assess the significance of local adaptation relative to environmental plasticity on the evolution of life history traits, we analysed the possible genetic basis of differences between pond- and stream-breeding fire salamanders (Salamandra salamandra) in Germany. These salamanders typically deposit their larvae in small streams, where they grow until they are sufficiently large to metamorphose. However, some populations in Western Germany use ponds as larval habitat. Because habitat quality of streams differs from that of ponds one expects life history differences in the pond animals, which may result either from a plastic response or through genetic differentiation (i.e. local adaptation). Using a phylogeographical analysis of mitochondrial D-loop sequences, we show that both stream and pond populations in Western Germany are derived from a single lineage that recolonized following the last glaciation. This finding suggests that pond breeding originated very recently. Our studies of habitat quality and metamorphic behaviour of larvae in natural ponds and streams disclosed that pond larvae experience a significantly reduced food supply and greater risk of drying than do stream larvae. Pond larvae metamorphose earlier at the cost of reduced mass. Common-environment experiments with pond and stream larvae show that metamorphic behaviour of pond larvae under limited-food conditions is determined genetically and is not simply a plastic response to the differing habitat conditions. These results show that phenotypic plasticity is less important than local adaptation in explaining differences in ecological diversification within this species and suggests the possibility of rapid evolution of genetic adaptations when new habitats are exploited.     

Adaptive phenotypic plasticity contributes to divergence between lake and river populations of an East African cichlid fish

Adaptive phenotypic plasticity and fixed genotypic differences have long been considered opposing strategies in adaptation. More recently, these mechanisms have been proposed to act complementarily and under certain conditions jointly facilitate evolution, speciation, and even adaptive radiations. Here, we investigate the relative contributions of adaptive phenotypic plasticity vs. local adaptation to fitness, using an emerging model system to study early phases of adaptive divergence, the generalist cichlid fish species Astatotilapia burtoni. We tested direct fitness consequences of morphological divergence between lake and river populations in nature by performing two transplant experiments in Lake Tanganyika. In the first experiment, we used wild‐caught juvenile lake and river individuals, while in the second experiment, we used F1 crosses between lake and river fish bred in a common garden setup. By tracking the survival and growth of translocated individuals in enclosures in the lake over several weeks, we revealed local adaptation evidenced by faster growth of the wild‐caught resident population in the first experiment. On the other hand, we did not find difference in growth between different types of F1 crosses in the second experiment, suggesting a substantial contribution of adaptive phenotypic plasticity to increased immigrant fitness. Our findings highlight the value of formally comparing fitness of wild‐caught and common garden‐reared individuals and emphasize the necessity of considering adaptive phenotypic plasticity in the study of adaptive divergence.     

Modifying effects of phenotypic plasticity on interactions among natural selection, adaptation and gene flow

Divergent natural selection, adaptive divergence and gene flow may interact in a number of ways. Recent studies have focused on the balance between selection and gene flow in natural populations, and empirical work has shown that gene flow can constrain adaptive divergence, and that divergent selection can constrain gene flow. A caveat is that phenotypic diversification may be under the direct influence of environmental factors (i.e. it may be due to phenotypic plasticity), in addition to partial genetic influence. In this case, phenotypic divergence may occur between populations despite high gene flow that imposes a constraint on genetic divergence. Plasticity may dampen the effects of natural selection by allowing individuals to rapidly adapt phenotypically to new conditions, thus slowing adaptive genetic divergence. On the other hand, plasticity may promote future adaptive divergence by allowing populations to persist in novel environments. Plasticity may promote gene flow between selective regimes by allowing dispersers to adapt to alternate conditions, or high gene flow may result in the selection for increased plasticity. Here I expand frameworks for understanding relationships among selection, adaptation and gene flow to include the effects of phenotypic plasticity in natural populations, and highlight its importance in evolutionary diversification.     

Allopatric and sympatric diversification within roach (Rutilus rutilus) of large pre‐alpine lakes

Intraspecific differentiation in response to divergent natural selection between environments is a common phenomenon in some northern freshwater fishes, especially salmonids and stickleback. Understanding why these taxa diversify and undergo adaptive radiations while most other fish species in the same environments do not, remains an open question. The possibility for intraspecific diversification has rarely been evaluated for most northern freshwater fish species. Here, we assess the potential for intraspecific differentiation between and within lake populations of roach (Rutilus rutilus)—a widespread and abundant cyprinid species—in lakes in which salmonids have evolved endemic adaptive radiations. Based on more than 3,000 polymorphic RADseq markers, we detected low but significant genetic differentiation between roach populations of two ultraoligotrophic lakes and between these and populations from other lakes. This, together with differentiation in head morphology and stable isotope signatures, suggests evolutionary and ecological differentiation among some of our studied populations. Next, we tested for intralacustrine diversification of roach within Lake Brienz, the most pristine lake surveyed in this study. We found significant phenotypic evidence for ecological intralacustrine differentiation between roach caught over a muddy substrate and those caught over a rocky substrate. However, evidence for intralacustrine genetic differentiation is at best subtle and phenotypic changes may therefore be mostly plastic. Overall, our findings suggest roach can differ between ecologically distinct lakes, but the extent of intralacustrine ecological differentiation is weak, which contrasts with the strong differentiation among endemic species of whitefish in the same lakes.     

Phenotypic plasticity in gene expression contributes to divergence of locally adapted populations of Fundulus heteroclitus

We examine the interaction between phenotypic plasticity and evolutionary adaptation using muscle gene expression levels among populations of the fish Fundulus heteroclitus acclimated to three temperatures. Our analysis reveals shared patterns of phenotypic plasticity due to thermal acclimation as well as non‐neutral patterns of variation among populations adapted to different thermal environments. For the majority of significant differences in gene expression levels, phenotypic plasticity and adaptation operate on different suites of genes. The subset of genes that demonstrate both adaptive differences and phenotypic plasticity, however, exhibit countergradient variation of expression. Thus, expression differences among populations counteract environmental effects, reducing the phenotypic differentiation between populations. Finally, gene‐by‐environment interactions among genes with non‐neutral patterns of expression suggest that the penetrance of adaptive variation depends on the environmental conditions experienced by the individual.     

The effect of top‐predator presence and phenotype on aquatic microbial communities

The presence of predators can impact a variety of organisms within the ecosystem, including microorganisms. Because the effects of fish predators and their phenotypic differences on microbial communities have not received much attention, we tested how the presence/absence, genotype, and plasticity of the predatory three‐spine stickleback (Gasterosteus aculeatus) influence aquatic microbes in outdoor mesocosms. We reared lake and stream stickleback genotypes on contrasting food resources to adulthood, and then added them to aquatic mesocosm ecosystems to assess their impact on the planktonic bacterial community. We also investigated whether the effects of fish persisted following the removal of adults, and the subsequent addition of a homogenous juvenile fish population. The presence of adult stickleback increased the number of bacterial OTUs and altered the size structure of the microbial community, whereas their phenotype affected bacterial community composition. Some of these effects were detectable after adult fish were removed from the mesocosms, and after juvenile fish were placed in the tanks, most of these effects disappeared. Our results suggest that fish can have strong short‐term effects on microbial communities that are partially mediated by phenotypic variation of fish.     

Microhabitat contributes to microgeographic divergence in threespine stickleback

Since the New Synthesis, most migration-selection balance theory has predicted that there should be negligible differentiation over small spatial scales (relative to dispersal), because gene flow should erode any effect of divergent selection. Nevertheless, there are classic examples of microgeographic divergence, which theory suggests can arise under specific conditions: exceptionally strong selection, phenotypic plasticity in philopatric individuals, or nonrandom dispersal. Here, we present evidence of microgeographic morphological variation within lake and stream populations of threespine stickleback (Gasterosteus aculeatus). It seems reasonable to assume that a given lake or stream population of fish is well-mixed. However, we found this assumption to be untenable. We examined trap-to-trap variation in 34 morphological traits measured on stickleback from 16 lakes and 16 streams. Most traits varied appreciably among traps within populations. Both between-trap distance and microhabitat characteristics such as depth and substrate explained some of the within-population morphological variance. Microhabitat was also associated with genotype at particular loci but there was no genetic isolation by distance, implying that heritable habitat preferences may contribute to microgeographic variation. Our study adds to growing evidence that microgeographic divergence can occur across small spatial scales within individuals’ daily dispersal neighborhood where gene flow is expected to be strong.