Phenotypic divergence among threespine stickleback that differ in nuptial coloration

By studying systems in their earliest stages of differentiation, we can learn about the evolutionary forces acting within and among populations and how those forces could contribute to reproductive isolation. Such an understanding would help us to better discern and predict how selection leads to the maintenance of multiple morphs within a species, rather than speciation. The postglacial adaptive radiation of the threespine stickleback (Gasterosteus aculeatus) is one of the best‐studied cases of evolutionary diversification and rapid, repeated speciation. Following deglaciation, marine stickleback have continually invaded freshwater habitats across the northern hemisphere and established resident populations that diverged innumerable times from their oceanic ancestors. Independent freshwater colonization events have yielded broadly parallel patterns of morphological differences in freshwater and marine stickleback. However, there is also much phenotypic diversity within and among freshwater populations. We studied a lesser‐known freshwater “species pair” found in southwest Washington, where male stickleback in numerous locations have lost the ancestral red sexual signal and instead develop black nuptial coloration. We measured phenotypic variation in a suite of traits across sites where red and black stickleback do not overlap in distribution and at one site where they historically co‐occurred. We found substantial phenotypic divergence between red and black morphs in noncolor traits including shape and lateral plating, and additionally find evidence that supports the hypothesis of sensory drive as the mechanism responsible for the evolutionary switch in color from red to black. A newly described third “mixed” morph in Connor Creek, Washington, differs in head shape and size from the red and black morphs, and we suggest that their characteristics are most consistent with hybridization between anadromous and freshwater stickleback. These results lay the foundation for future investigation of the underlying genetic basis of this phenotypic divergence as well as the evolutionary processes that may drive, maintain, or limit divergence among morphs.     

Evolutionary diversification of opercle shape in Cook Inlet threespine stickleback

We investigated the evolution of a large facial bone, the opercle (OP), in lake populations of the threespine stickleback that were founded by anadromous ancestors, in Cook Inlet, Alaska. Recent studies characterized OP variation among marine and lake populations and mapped a quantitative trait locus with a large influence on OP shape. Using populations from diverse environments and independent evolutionary histories, we examined divergence of OP shape from that of the anadromous ancestor. We report preliminary evidence for divergence between benthic and generalist lake ecotypes, necessitating further investigation. Furthermore, rapid divergence of OP shape has occurred in a lake population that was founded by anadromous stickleback in the 1980s, which is consistent with divergence of other phenotypic traits and with OP diversification in other lake populations. By contrast, there has been limited evolution of OP shape in a second lake population that may have experienced a genetic bottleneck early in its history and lacks genetic variation for OP divergence. Taken together, the results obtained from these two populations are consistent with studies of other stickleback phenotypic traits that implicate ancestral variation in postglacial adaptive radiation of threespine stickleback in fresh water.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 832–844.     

Genetic relationships among marine and freshwater populations of the European three-spined stickleback (Gasterosteus aculeatus) revealed by microsatellites

To assess the population genetic structure of the three-spined stickleback, Gasterosteus aculeatus, variability at 18 microsatellite loci was examined in 1724 individuals from 74 locations covering most of the species distribution range in Europe. The results revealed high overall degree of differentiation (F(ST) = 0.21) but contrasting level of divergence and genetic variability between habitat types. Marine populations were genetically relatively uniform even across great geographical distances as compared to substantial differentiation among freshwater populations. Analysis of molecular variance indicated low but significant (2.7%) variation in allele frequencies between geographical regions, but a negligible effect of habitat type (0.2%). The phylogenetic pattern was not explained by habitat type, but a weak signal of populations clustering according to geographical or water system origin was found. The results support the view that three-spined stickleback marine ancestors colonized northern European fresh waters during the postglacial marine submergence c. 10,000 years ago, whereas in the Mediterranean region colonization probably dates back to the Pleistocene. The independent origins of river and lake populations indicate that they originate from multiple colonizations rather than sharing common ancestry. In the continuous marine environment, the low degree of differentiation among populations can be explained by gene flow among subpopulations and large effective population size buffering divergence in neutral markers. In contrast, among postglacially established freshwater populations differentiation appears to be driven by genetic drift and isolation. The stepwise mutations appear to have contributed to the population differentiation in the southern part of the three-spined stickleback distribution range.     

Resource Competition and Coevolution in Sticklebacks

Threespine stickleback in young postglacial lakes provide a compelling example of coevolution between species that compete for resources. Coexisting pairs of stickleback species are highly divergent in habitat, diet, and body size and shape, whereas stickleback occurring alone in lakes are intermediate. We used experiments in ponds to test mechanisms of divergence between coexisting species. The results support the hypothesis of coevolution by resource competition between stickleback, but we found evidence that interactions with natural enemies also contribute to divergence. Natural selection arising from these interactions selects against intermediate phenotypes, included hybrids, and thus has contributed to the origin and persistence of stickleback species.     

Intraguild predation leads to genetically based character shifts in the threespine stickleback

Intraguild predation is a common ecological interaction that occurs when a species preys upon another species with which it competes. The interaction is potentially a mechanism of divergence between intraguild prey (IG‐prey) populations, but it is unknown if cases of character shifts in IG‐prey are an environmental or evolutionary response. We investigated the genetic basis and inducibility of character shifts in threespine stickleback from lakes with and without prickly sculpin, a benthic intraguild predator (IG‐predator). Wild populations of stickleback sympatric with sculpin repeatedly show greater defensive armor and water column height preference. We laboratory‐raised stickleback from lakes with and without sculpin, as well as marine stickleback, and found that differences between populations in armor, body shape, and behavior persisted in a common garden. Within the common garden, we raised stickleback half‐families from multiple populations in the presence and absence of sculpin. Although the presence of sculpin induced trait changes in the marine stickleback, we did not observe an induced response in the freshwater stickleback. Behavioral and morphological trait differences between freshwater populations thus have a genetic basis and suggest an evolutionary response to intraguild predation.     

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.     

Evolutionary divergence and possible incipient speciation in post-glacial populations of a cosmopolitan aquatic plant

Habitat configuration is expected to have a major influence on genetic exchange and evolutionary divergence among populations. Aquatic organisms occur in two fundamentally different habitat types, the sea and freshwater lakes, making them excellent models to study the contrasting effects of continuity vs. isolation on genetic divergence. We compared the divergence in post-glacial populations of a cosmopolitan aquatic plant, the pondweed Potamogeton pectinatus that simultaneously occurs in freshwater lakes and coastal marine sites. Relative levels of gene flow were inferred in 12 lake and 14 Baltic Sea populations in northern Germany using nine highly polymorphic microsatellite markers developed for P. pectinatus. We found highly significant isolation-by-distance in both habitat types (P < 0.001). Genetic differentiation increased approximately 2.5-times faster among freshwater populations compared with those from the Baltic Sea. As different levels of genetic drift or population history cannot explain these differences, higher population connectivity in the sea relative to freshwater populations is the most likely source of contrasting evolutionary divergence. These findings are consistent with the notion that freshwater angiosperms are more conducive to allopatric speciation than their life-history counterparts in the sea, the relative species poor seagrasses. Surprisingly, population pairs from different habitat types revealed almost maximal genetic divergence expected for complete reproductive isolation, regardless of their respective geographical distance. Hence, the barrier to gene flow between lake and sea habitat types cannot be due to dispersal limitation. We may thus have identified a case of rapid incipient speciation in post-glacial populations of a widespread aquatic plant.     

Contemporary and historical evolutionary processes interact to shape patterns of within-lake phenotypic divergences in polyphenic pumpkinseed sunfish, Lepomis gibbosus

Historical and contemporary evolutionary processes can both contribute to patterns of phenotypic variation among populations of a species. Recent studies are revealing how interactions between historical and contemporary processes better explain observed patterns of phenotypic divergence than either process alone. Here, we investigate the roles of evolutionary history and adaptation to current environmental conditions in structuring phenotypic variation among polyphenic populations of sunfish inhabiting 12 postglacial lakes in eastern North America. The pumpkinseed sunfish polyphenism includes sympatric ecomorphs specialized for littoral or pelagic lake habitats. First, we use population genetic methods to test the evolutionary independence of within-lake phenotypic divergences of ecomorphs and to describe patterns of genetic structure among lake populations that clustered into three geographical groupings. We then used multivariate analysis of covariance (MANCOVA) to partition body shape variation (quantified with geometric morphometrics) among the effects of evolutionary history (reflecting phenotypic variation among genetic clusters), the shared phenotypic response of all populations to alternate habitats within lakes (reflecting adaptation to contemporary conditions), and unique phenotypic responses to habitats within lakes nested within genetic clusters. All effects had a significant influence on body form, but the effects of history and the interaction between history and contemporary habitat were larger than contemporary processes in structuring phenotypic variation. This highlights how divergence can be better understood against a known backdrop of evolutionary history.     

Historical contingency and ecological determinism interact to prime speciation in sticklebacks, Gasterosteus

Historical contingency and determinism are often cast as opposing paradigms under which evolutionary diversification operates. It may be, however, that both factors act together to promote evolutionary divergence, although there are few examples of such interaction in nature. We tested phylogenetic predictions of an explicit historical model of divergence (double invasions of freshwater by marine ancestors) in sympatric species of three-spined sticklebacks (Gasterosteus aculeatus) where determinism has been implicated as an important factor driving evolutionary novelty. Microsatellite DNA variation at six loci revealed relatively low genetic variation in freshwater populations, supporting the hypothesis that they were derived by colonization of freshwater by more diverse marine ancestors. Phylogenetic and genetic distance analyses suggested that pairs of sympatric species have evolved multiple times, further implicating determinism as a factor in speciation. Our data also supported predictions based on the hypothesis that the evolution of sympatric species was contingent upon 'double invasions' of postglacial lakes by ancestral marine sticklebacks. Sympatric sticklebacks, therefore, provide an example of adaptive radiation by determinism contingent upon historical conditions promoting unique ecological interactions, and illustrate how contingency and determinism may interact to generate geographical variation in species diversity     

Environmental correlates of adaptive diversification in postglacial freshwater fishes

Determining how environmental conditions contribute to divergence among populations and drive speciation is fundamental to resolving mechanisms and understanding outcomes in evolutionary biology. Postglacial freshwater fish species in the Northern Hemisphere are ideal biological systems to explore the effects of environment on diversification in morphology, ecology, and genetics (ecomorph divergences) within lakes. To date, various environmental factors have been implicated in the presence of multiple ecomorphs within particular lakes or regions. However, concerted evidence for generalizable patterns in environmental variables associated with speciation across geographical regions and across species and genera has been lacking. Here, we aimed to identify key biotic and abiotic factors associated with ecological divergence of postglacial freshwater fish species into multiple sympatric ecomorphs, focusing on species in the well-studied, widespread, and co-distributed genera Gasterosteus, Salvelinus, and Coregonus (stickleback, charr, and whitefish, respectively). We found that the presence of multiple sympatric ecomorphs tended to be associated with increasing lake surface area, maximum depth, and nutrient availability. In addition, predation, competition, and prey availability were suggested to play a role in divergence into multiple ecomorphs, but the effects of biotic factors require further study. Although we identified several environmental factors correlated with the presence of multiple ecomorphs, there were substantial data gaps across species and regions. An improved understanding of these systems may provide insight into both generalizable environmental factors involved in speciation in other systems, and potential ecological and evolutionary responses of species complexes when these variables are altered by environmental change.     

Repeated and predictable patterns of ecotypic differentiation during a biological invasion: lake–stream divergence in parapatric Swiss stickleback

The relative importance of ecological selection and geographical isolation in promoting and constraining genetic and phenotypic differentiation among populations is not always obvious. Interacting with divergent selection, restricted opportunity for gene flow may in some cases be as much a cause as a consequence of adaptation, with the latter being a hallmark of ecological speciation. Ecological speciation is well studied in parts of the native range of the three‐spined stickleback. Here, we study this process in a recently invaded part of its range. Switzerland was colonized within the past 140 years from at least three different colonization events involving different stickleback lineages. They now occupy diverse habitats, ranging from small streams to the pelagic zone of large lakes. We use replicated systems of parapatric lake and stream populations, some of which trace their origins to different invasive lineages, to ask (i) whether phenotypic divergence occurred among populations inhabiting distinct habitats, (ii) whether trajectories of phenotypic divergence follow predictable parallel patterns and (iii) whether gene flow constrains divergent adaptation or vice versa. We find consistent phenotypic divergence between populations occupying distinct habitats. This involves parallel evolution in several traits with known ecological relevance in independent evolutionary lineages. Adaptive divergence supersedes homogenizing gene flow even at a small spatial scale. We find evidence that adaptive phenotypic divergence places constraints on gene flow over and above that imposed by geographical distance, signalling the early onset of ecological speciation.     

Phylogeography and adaptation genetics of stickleback from the Haida Gwaii archipelago revealed using genome‐wide single nucleotide polymorphism genotyping

Threespine stickleback populations are model systems for studying adaptive evolution and the underlying genetics. In lakes on the Haida Gwaii archipelago (off western Canada), stickleback have undergone a remarkable local radiation and show phenotypic diversity matching that seen throughout the species distribution. To provide a historical context for this radiation, we surveyed genetic variation at >1000 single nucleotide polymorphism (SNP) loci in stickleback from over 100 populations. SNPs included markers evenly distributed throughout genome and candidate SNPs tagging adaptive genomic regions. Based on evenly distributed SNPs, the phylogeographic pattern differs substantially from the disjunct pattern previously observed between two highly divergent mtDNA lineages. The SNP tree instead shows extensive within watershed population clustering and different watersheds separated by short branches deep in the tree. These data are consistent with separate colonizations of most watersheds, despite underlying genetic connections between some independent drainages. This supports previous suppositions that morphological diversity observed between watersheds has been shaped independently, with populations exhibiting complete loss of lateral plates and giant size each occurring in several distinct clades. Throughout the archipelago, we see repeated selection of SNPs tagging candidate freshwater adaptive variants at several genomic regions differentiated between marine–freshwater populations on a global scale (e.g. EDA, Na/K ATPase). In estuarine sites, both marine and freshwater allelic variants were commonly detected. We also found typically marine alleles present in a few freshwater lakes, especially those with completely plated morphology. These results provide a general model for postglacial colonization of freshwater habitat by sticklebacks and illustrate the tremendous potential of genome‐wide SNP data sets hold for resolving patterns and processes underlying recent adaptive divergences.     

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.     

A benthic predatory fish does not cause selection on armour traits in three‐spined stickleback Gasterosteus aculeatus (Gasterosteiformes: Gasterosteidae)

Predation can promote divergence between prey populations and contribute to ecological speciation. In theory, predators can also constrain prey population divergence. In coastal British Columbia, Canada, Gasterosteus aculeatus (three‐spined stickleback) species pairs only occur in lakes with a single species of predatory fish: Oncorhynchus clarkii (the cutthroat trout). Similar lakes containing additional predatory fish species (Cottus asper, prickly sculpins; Oncorhynchus mykiss, rainbow trout) contain only single species of morphologically intermediate stickleback, suggesting that these predators prevent the coexistence of stickleback species pairs. We conducted a mesocosm experiment to investigate how prickly sculpins might constrain divergence, by quantifying their impact on survival and natural selection on antipredator (armour) traits in F₂ stickleback from a cross between ecologically divergent populations. We tested three hypotheses: (1) sculpin predation on sticklebacks reduces survival in a way that could result in their exclusion from certain niches; (2) sculpins compete with stickleback; (3) sculpins respond to prey vulnerabilities in similar ways to cutthroat trout, tending to constrain rather than to enhance divergence. We found that sculpins significantly reduce stickleback survival, that their presence per se does not reduce growth in stickleback, and that predation did not result in selection on any of the armour traits measured, or on gill raker length, which is an important trophic trait. These results tend to refute hypotheses (2) and (3), while supporting hypothesis (1). © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 104 , 877–885.     

Habitat‐specific trends in ontogeny of body shape in stickleback from coastal archipelago: Potential for rapid shifts in colonizing populations

We investigated ontogenetic trends in body shape of 54 freshwater (48 lake, seven stream) and six anadromous populations of threespine stickleback (Gasterosteus aculeatus L.) from the Haida Gwaii archipelago off the west coast of Canada. Multivariate analysis of covariance on the partial warp scores generated from 12 homologous landmarks on 1,958 digital images of subadult and adult male stickleback indicated that there was considerable variability of population ontogenetic slopes. We used discriminant function analysis to quantify body shape and determined that anadromous stickleback, which are ancestral to the freshwater populations, have a strongly negative ontogenetic slope (?5.62; increased streamlining with increased size). All freshwater populations exhibit a more positive slope (91% differed significantly from the marine slope), with the differences being most accentuated in populations from ponds and streams. In pristine lakes, ontogenetic slope could be predicted by lake volume as well as multivariate measures of habitat. Evidence from field transplant experiments of one of the intact populations indicates a rapid change (5 years) from allometric to isometric growth, equivalent to about half of the total slope variation among intact populations on the archipelago. We interpret this shift as developmental plasticity and suggest this may comprise the precursor for selection of optimal body shapes in these stickleback populations. J. Morphol., 2011. © 2011 Wiley‐Liss, Inc.     

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.     

Adaptation and constraint in a stickleback radiation

The evolution of threespine sticklebacks in freshwater lakes constitutes a well‐studied example of a phenotypic radiation that has produced numerous instances of parallel evolution, but the exact selective agents that drive these changes are not yet fully understood. We present a comparative study across 74 freshwater populations of threespine stickleback in Norway to test whether evolutionary changes in stickleback morphology are consistent with adaptations to physical parameters such as lake depth, lake area, lake perimeter and shoreline complexity, variables thought to reflect different habitats and feeding niches. Only weak indications of adaptation were found. Instead, populations seem to have diversified in phenotypic directions consistent with allometric scaling relationships. This indicates that evolutionary constraints may have played a role in structuring phenotypic variation across freshwater populations of stickleback. We also tested whether the number of lateral plates evolved in response to lake calcium levels, but found no evidence for this hypothesis.     

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.     

Variation of telencephalon morphology of the threespine stickleback (Gasterosteus aculeatus) in relation to inferred ecology

We tested the hypothesis that increased telencephalon size has evolved in threespine stickleback fish (Gasterosteus aculeatus) from structurally complex habitats using field‐caught samples from one sea‐run (ancestral) and 18 ecologically diverse freshwater (descendant) populations. Freshwater habitats ranged from shallow, structurally complex lakes with benthic‐foraging stickleback (benthics), to deeper, structurally simple lakes in which stickleback depend more heavily on plankton for prey (generalists). Contrary to our expectations, benthics had smaller telencephala than generalists, but the shape of the telencephalon of the sea‐run and benthic populations were more convex laterally. Convex telencephalon shape may indicate enlargement of the dorsolateral region, which is homologous with the tetrapod hippocampus. Telencephalon morphology is also sexually dimorphic, with larger, less convex telencephala in males. Freshwater stickleback from structurally complex habitats have retained the ancestral telencephalon morphology, but populations that feed more in open habitats on plankton have evolved larger, laterally concave telencephala.     

Genetic and morphometric divergence in threespine stickleback in the Chignik catchment,Alaska

Divergent selection pressures induced by different environmental conditions typically lead to variation in life history, behavior, and morphology. When populations are locally adapted to their current environment, selection may limit movement into novel sites, leading to neutral and adaptive genetic divergence in allopatric populations. Subsequently, divergence can be reinforced by development of pre‐ or postzygotic barriers to gene flow. The threespine stickleback, Gasterosteus aculeatus, is a primarily marine fish that has invaded freshwater repeatedly in postglacial times. After invasion, the established freshwater populations typically show rapid diversification of several traits as they become reproductively isolated from their ancestral marine population. In this study, we examine the genetic and morphometric differentiation between sticklebacks living in an open system comprising a brackish water lagoon, two freshwater lakes, and connecting rivers. By applying a set of microsatellite markers, we disentangled the genetic relationship of the individuals across the diverse environments and identified two genetic populations: one associated with brackish and the other with the freshwater environments. The “brackish” sticklebacks were larger and had a different body shape than those in freshwater. However, we found evidence for upstream migration from the brackish lagoon into the freshwater environments, as fish that were genetically and morphometrically similar to the lagoon fish were found in all freshwater sampling sites. Regardless, few F1‐hybrids were identified, and it therefore appears that some pre‐ and/or postzygotic barriers to gene flow rather than geographic distance are causing the divergence in this system.