Parallel and nonparallel ecological,morphological and genetic divergence in lake–stream stickleback from a single catchment

Parallel phenotypic evolution in similar environments has been well studied in evolutionary biology; however, comparatively little is known about the influence of determinism and historical contingency on the nature, extent and generality of this divergence. Taking advantage of a novel system containing multiple lake–stream stickleback populations, we examined the extent of ecological, morphological and genetic divergence between three‐spined stickleback present in parapatric environments. Consistent with other lake–stream studies, we found a shift towards a deeper body and shorter gill rakers in stream fish. Morphological shifts were concurrent with changes in diet, indicated by both stable isotope and stomach contents analysis. Performing a multivariate test for shared and unique components of evolutionary response to the distance gradient from the lake, we found a strong signature of parallel adaptation. Nonparallel divergence was also present, attributable mainly to differences between river locations. We additionally found evidence of genetic substructuring across five lake–stream transitions, indicating that some level of reproductive isolation occurs between populations in these habitats. Strong correlations between pairwise measures of morphological, ecological and genetic distance between lake and stream populations supports the hypothesis that divergent natural selection between habitats drives adaptive divergence and reproductive isolation. Lake–stream stickleback divergence in Lough Neagh provides evidence for the deterministic role of selection and supports the hypothesis that parallel selection in similar environments may initiate parallel speciation.     

Analyzing nested variation in the body form of Lepomid sunfishes

Phylogenetic hierarchies are often composed of younger diverging lineages nested within older diverging lineages. Comparing phenotypic variation among several hierarchical levels can be used to test hypotheses about selection, phenotypic evolution and speciation. Such hierarchical comparisons have only been performed in threespine stickleback, and so here we use a hierarchical pattern of divergences between near-shore littoral and off-shore pelagic habitats to test for selection on the evolution of body form in Lepomis sunfish in lakes. We compare variation in external body form between fish from littoral and pelagic habitats at three levels: among ecomorphs within individual lake populations (intrapopulation), among populations of the same species in different lakes (interpopulation), and between bluegill and pumpkinseed sunfish species (interspecifically). Using geometric morphometric methods, we first demonstrate that interpopulation variation in mean body form of pumpkinseed sunfish varies with the presence of pelagic habitat. We then incorporate these results with existing data in order to test the similarity of phenotypic divergence between littoral and pelagic habitats at different hierarchical levels. Parallel relationships between certain body form traits (head length, caudal length and pectoral length) and habitat occur at all three levels suggesting that selection persistently acts at all levels to diversify these traits and so may contribute to species formation. For other traits (caudal depth and pectoral altitude), divergence between habitats is inconsistent at different hierarchical levels. Thus, nested biological variation in Lepomid body form reflects a history of deterministic selection and historical contingency, and also identifies traits that likely have likely influenced fitness and so serve important functions.     

Divergent selection maintains adaptive differentiation despite high gene flow between sympatric rainbow smelt ecotypes (Osmerus mordax Mitchill)

In this study, we investigate the relative role of historical factors and evolutionary forces in promoting population differentiation in a new case of sympatric dwarf and normal ecotypes of the rainbow smelt (Osmerus mordax Mitchill) in Lac Saint-Jean (Québec, Canada). Our first objective was to test the hypothesis that the evolution of sympatric smelt ecotypes in Lac Saint-Jean has been contingent upon the secondary contact between two evolutionary lineages in postglacial times. Secondly, the QST method was applied to test the null hypothesis that the extent of phenotypic differences relative to that of neutral marker variation would be similar in comparisons involving populations within and among ecotypes. Thirdly, we applied a quantitative-genetic method as an exploratory assessment as to whether the amount of gene flow observed between populations could affect divergence in adaptive traits under specific conditions. This study revealed a unique situation of dwarf and normal smelt ecotypes that are, respectively, characterized by selmiparous and iteroparous life histories and the occurrence in each of two genetically distinct populations that synchronously use the same spawning habitat in two tributaries. Historical contingency has apparently played little role in the origin of these populations. In contrast, an important role of divergent natural selection in driving their phenotypic divergence was suggested. While divergent selection has apparently been strong enough to maintain phenotypic differentiation in the face of migration, this study suggests that gene flow has been sufficiently important to modulate the extent of adaptive differentiation being achieved between ecotypes, unless the extent of stabilizing selection acting on smelt ecotypes is much more pronounced than usually reported in natural populations.     

Layers of contingency shroud pervasive ecological divergence in a local radiation of land snails

The predictability of evolution depends on the roles that selection and historical contingency play in determining its outcomes, but the relative importance of these evolutionary mechanisms has attracted considerable debate. One view is that historical events have such a profound impact on the genetic structure of populations that patterns of phenotypic evolution are essentially unpredictable. The opposing view is that selection is so powerful that evolutionary change is primarily deterministic, and thus highly predictable. By controlling for the effects of phylogeny, geographic location and habitat, this study examined the relative roles of contingency and determinism in a local radiation of land snails, genus Rhagada, in a continental archipelago. Informed by previous studies on a single island, which revealed a strong association between low‐spired shells and rocky habitats, 28 population pairs were sampled in directly adjoining rocky and spinifex plain habitats. When considered in their respective pairs, the effect of habitat was remarkably consistent, with lower‐spired shells observed in the rocky habitat in 24 of the comparisons. However, when analyzed outside the context of those pairs, the association was obscured by broad variation in shell shape within habitat types and among lineages. These results reveal the complex nature of a morphological radiation; while the pattern of ecological divergence is highly predictable at the scale that selection acts, deterministic evolution is largely obscured by phylogenetic and population history.     

Environmental factors associated with genetic and phenotypic divergence among sympatric populations of Arctic charr (Salvelinus alpinus)

The mechanisms by which phenotypic and genetic divergence may occur among sympatric, conspecific populations have been widely discussed but are still not well understood. Possible mechanisms include assortative mating based on morphology or variation in the reproductive behaviour of phenotypes, and both have been suggested to be relevant to the differentiation of salmonid populations in post-glacial lakes. Here, we studied Arctic charr (Salvelinus alpinus) in Windermere, where putative populations are defined by spatial and temporal variation in spawning. Genetic differentiation was assessed based on nine microsatellite loci, and phenotypic variation was assessed from morphometric characters. We test hypotheses about the relative role of morphology, spawning season and spawning habitat in the evolution of genetic divergence among these populations. Distinct from other lake systems, we find that both morphological and genetic differentiation are restricted primarily to one of two interconnecting basins, that genetic and morphological differentiation are decoupled in this lake and that both phenotype and environment have changed over the last 20 years. The implication is that breeding habitat plays a primary role in isolating populations that differentiate by drift and that phenotypically plastic changes, potentially related to foraging specializations, have either become secondarily decoupled from the genetically defined populations or were never fundamental in driving the evolution of genetic diversity in the Windermere system.     

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.     

Spatial phenotypic and genetic structure of threespine stickleback (Gasterosteus aculeatus) in a heterogeneous natural system,Lake Mývatn,Iceland

Eco‐evolutionary responses of natural populations to spatial environmental variation strongly depend on the relative strength of environmental differences/natural selection and dispersal/gene flow. In absence of geographic barriers, as often is the case in lake ecosystems, gene flow is expected to constrain adaptive divergence between environments – favoring phenotypic plasticity or high trait variability. However, if divergent natural selection is sufficiently strong, adaptive divergence can occur in face of gene flow. The extent of divergence is most often studied between two contrasting environments, whereas potential for multimodal divergence is little explored. We investigated phenotypic (body size, defensive structures, and feeding morphology) and genetic (microsatellites) structure in threespine stickleback (Gasterosteus aculeatus) across five habitat types and two basins (North and South) within the geologically young and highly heterogeneous Lake Mývatn, North East Iceland. We found that (1) North basin stickleback were, on average, larger and had relatively longer spines than South basin stickleback, whereas (2) feeding morphology (gill raker number and gill raker gap width) differed among three of five habitat types, and (3) there was only subtle genetic differentiation across the lake. Overall, our results indicate predator and prey mediated phenotypic divergence across multiple habitats in the lake, in face of gene flow.     

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.     

PHENOTYPE-DEPENDENT NATIVE HABITAT PREFERENCE FACILITATES DIVERGENCE BETWEEN PARAPATRIC LAKE AND STREAM STICKLEBACK

Adaptive divergence between adjoining populations reflects a balance between the diversifying effect of divergent selection and the potentially homogenizing effect of gene flow. In most models of migration-selection balance, gene flow is assumed to reflect individuals' inherent capacity to disperse, without regard to the match between individuals' phenotypes and the available habitats. However, habitat preferences can reduce dispersal between contrasting habitats, thereby alleviating migration load and facilitating adaptive divergence. We tested whether habitat preferences contribute to adaptive divergence in a classic example of migration-selection balance: parapatric lake and stream populations of three-spine stickleback ( Gasterosteus aculeatus ). Using a mark-transplant-recapture experiment on morphologically divergent parapatric populations, we showed that 90% of lake and stream stickleback returned to their native habitat, reducing migration between habitats by 76%. Furthermore, we found that dispersal into a nonnative habitat was phenotype dependent. Stream fish moving into the lake were morphologically more lake-like than those returning to the stream (and the converse for lake fish entering the stream). The strong native habitat preference documented here increases the extent of adaptive divergence between populations two- to fivefold relative to expectations with random movement. These results illustrate the potential importance of adaptive habitat choice in driving parapatric divergence.     

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.     

Contemporary ecotypic divergence during a recent range expansion was facilitated by adaptive introgression

Although rapid phenotypic evolution during range expansion associated with colonization of contrasting habitats has been documented in several taxa, the evolutionary mechanisms that underlie such phenotypic divergence have less often been investigated. A strong candidate for rapid ecotype formation within an invaded range is the three‐spine stickleback in the Lake Geneva region of central Europe. Since its introduction only about 140 years ago, it has undergone a significant expansion of its range and its niche, now forming phenotypically differentiated parapatric ecotypes that occupy either the pelagic zone of the large lake or small inlet streams, respectively. By comparing museum collections from different times with contemporary population samples, we here reconstruct the evolution of parapatric phenotypic divergence through time. Using genetic data from modern samples, we infer the underlying invasion history. We find that parapatric habitat‐dependent phenotypic divergence between the lake and stream was already present in the first half of the twentieth century, but the magnitude of differentiation increased through time, particularly in antipredator defence traits. This suggests that divergent selection between the habitats occurred and was stable through much of the time since colonization. Recently, increased phenotypic differentiation in antipredator defence traits likely results from habitat‐dependent selection on alleles that arrived through introgression from a distantly related lineage from outside the Lake Geneva region. This illustrates how hybridization can quickly promote phenotypic divergence in a system where adaptation from standing genetic variation was constrained.     

Quick divergence but slow convergence during ecotype formation in lake and stream stickleback pairs of variable age

When genetic constraints restrict phenotypic evolution, diversification can be predicted to evolve along so‐called lines of least resistance. To address the importance of such constraints and their resolution, studies of parallel phenotypic divergence that differ in their age are valuable. Here, we investigate the parapatric evolution of six lake and stream threespine stickleback systems from Iceland and Switzerland, ranging in age from a few decades to several millennia. Using phenotypic data, we test for parallelism in ecotypic divergence between parapatric lake and stream populations and compare the observed patterns to an ancestral‐like marine population. We find strong and consistent phenotypic divergence, both among lake and stream populations and between our freshwater populations and the marine population. Interestingly, ecotypic divergence in low‐dimensional phenotype space (i.e. single traits) is rapid and seems to be often completed within 100 years. Yet, the dimensionality of ecotypic divergence was highest in our oldest systems and only there parallel evolution of unrelated ecotypes was strong enough to overwrite phylogenetic contingency. Moreover, the dimensionality of divergence in different systems varies between trait complexes, suggesting different constraints and evolutionary pathways to their resolution among freshwater systems.     

Drift effects on the multivariate floral phenotype of Calceolaria polyrhiza during a post‐glacial expansion in Patagonia

Quaternary environmental changes substantially impacted the landscape and promoted rapid evolutionary changes in many species; however, analyses of adaptive phenotypic variation in plants have usually neglected the underlying historical context. Here, we associate phylogeography and phenotypic evolution by analysing the divergence of Calceolaria polyrhiza multivariate floral phenotype after a Pleistocene post‐glacial expansion in Patagonia. Phenotypic matrix ( P ) properties (size, shape, orientation and phenotypic integration) of six refugium and six recent populations from two different phylogroups were compared following different approaches. We found that P ‐matrix shape and orientation remained stable despite the strong phylogeographic footprint of post‐glacial expansion. However, average proportional reductions in matrix size supported the expectation that drift had a significant effect on the floral phenotype in the northern phylogroup. When phylogeographic history was not included in the analyses, the results overestimated phenotypic differences, whereas under explicit phylogeographic control, drift appeared as the best explanation for matrix differences. In general, recent populations showed a larger phenotypic divergence among them, but a lower overall phenotypic variation than refugium populations. Random Skewers analyses indicated a lower potential response to selection in recently colonized populations than in refugium populations. We discuss that the combination of phylogeographic analyses with geographical distribution of functional phenotypic (genotypic) variation is critical not only to understand how historical effects influence adaptive evolution, but also to improve field comparisons in evolutionary ecology studies.     

Spatial predictors of genomic and phenotypic variation differ in a lowland Middle American bird (Icterus gularis)

Spatial patterns of intraspecific variation are shaped by geographical distance among populations, historical changes in gene flow and interactions with local environments. Although these factors are not mutually exclusive and operate on both genomic and phenotypic variation, it is unclear how they affect these two axes of variation. We address this question by exploring the predictors of genomic and phenotypic divergence in Icterus gularis, a broadly distributed Middle American bird that exhibits marked geographical variation in body size across its range. We combined a comprehensive single nucleotide polymorphism and phenotypic data set to test whether genome‐wide genetic and phenotypic differentiation are best explained by (i) isolation by distance, (ii) isolation by history or (iii) isolation by environment. We find that the pronounced genetic and phenotypic variation in I. gularis are only partially correlated and differ regarding spatial predictors. Whereas genomic variation is largely explained by historical barriers to gene flow, phenotypic diversity can be best predicted by contemporary environmental heterogeneity. Our genomic analyses reveal strong phylogeographical structure coinciding with the Chivela Pass at the Isthmus of Tehuantepec that was formed during the Pleistocene, when populations were isolated in north–south refugia. In contrast, we found a strong association between body size and environmental variables, such as temperature and precipitation. The relationship between body size and local climate is consistent with a pattern produced by either natural selection or environmental plasticity. Overall, these results provide empirical evidence for why phenotypic and genomic data are often in conflict in taxonomic and phylogeographical studies.     

Introgressive hybridization erodes morphological divergence between lentic and lotic habitats in an endangered minnow

Introgressive hybridization may erode phenotypic divergence along environmental gradients, collapsing locally adapted populations into a hybrid swarm. Alternatively, introgression may promote phenotypic divergence by providing variation on which natural selection can act. In freshwater fishes, water flow often selects for divergent morphological traits in lake versus stream habitats. We tested the effects of introgression on lake–stream morphological divergence in the minnow Owens Tui Chub (Siphateles bicolor snyderi), which has been rendered endangered by introgession from the introduced Lahontan Tui Chub (Siphateles bicolor obesa). Using geometric morphometric analysis of 457 individual Tui Chub from thirteen populations, we found that both native and introgressing parent taxa exhibited divergent body and caudal fin shapes in lake versus stream habitats, but their trajectories of divergence were distinct. In contrast, introgressed populations exhibited intermediate body and caudal fin shapes that were not differentiated by habitat type, indicating that introgression has eroded phenotypic divergence along the lentic–lotic gradient throughout the historic range of the Owens Tui Chub. Individuals within hybrid populations were less morphologically variable than those within parent populations, suggesting hybrid adaptation to selective agents other than water flow or loss of variance by drift.     

Parallel and nonparallel aspects of ecological, phenotypic, and genetic divergence across replicate population pairs of lake and stream stickleback

Parallel (or convergent) evolution provides strong evidence for a deterministic role of natural selection: similar phenotypes evolve when independent populations colonize similar environments. In reality, however, independent populations in similar environments always show some differences: some nonparallel evolution is present. It is therefore important to explicitly quantify the parallel and nonparallel aspects of trait variation, and to investigate the ecological and genetic explanations for each. We performed such an analysis for threespine stickleback (Gasterosteus aculeatus) populations inhabiting lake and stream habitats in six independent watersheds. Morphological traits differed in the degree to which lake-stream divergence was parallel across watersheds. Some aspects of this variation were correlated with ecological variables related to diet, presumably reflecting the strength and specifics of divergent selection. Furthermore, a genetic scan revealed some markers that diverged between lakes and streams in many of the watersheds and some that diverged in only a few watersheds. Moreover, some of the lake-stream divergence in genetic markers was associated within some of the lake-stream divergence in morphological traits. Our results suggest that parallel evolution, and deviations from it, are primarily the result of natural selection, which corresponds in only some respects to the dichotomous habitat classifications frequently used in such studies.     

Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks

Comparisons of neutral marker and quantitative trait divergence can provide important insights into the relative roles of natural selection and neutral genetic drift in population differentiation. We investigated phenotypic and genetic differentiation among Fennoscandian threespine stickleback (Gasterosteus aculeatus) populations, and found that the highest degree of differentiation occurred between sea and freshwater habitats. Within habitats, morphological divergence was highest among the different freshwater populations. Pairwise phenotypic and neutral genetic distances among populations were positively correlated, suggesting that genetic drift may have contributed to the morphological differentiation among habitats. On the other hand, the degree of phenotypic differentiation (PST) clearly surpassed the neutral expectation set by FST, suggesting a predominant role for natural selection over genetic drift as an explanation for the observed differentiation. However, separate PST/FST comparisons by habitats revealed that body shape divergence between lake and marine populations, and even among marine populations, can be strongly influenced by natural selection. On the other hand, genetic drift can play an important role in the differentiation among lake populations.     

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild‐caught individuals. Here, we investigated the contribution of habitat‐specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake–river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake‐like pond set‐up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set‐up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set‐up were based on a plastic response. Finally, gene expression and bacterial communities of wild‐caught individuals and individuals acclimatized to lake‐like pond conditions showed shifts underlying adaptive phenotypic plasticity.     

Local adaptation in the rock pocket mouse (Chaetodipus intermedius): natural selection and phylogenetic history of populations

Elucidating the causes of population divergence is a central goal of evolutionary biology. Rock pocket mice, Chaeotdipus intermedius, are an ideal system in which to study intraspecific phenotypic divergence because of the extensive color variation observed within this species. Here, we investigate whether phenotypic variation in color is correlated with local environmental conditions or with phylogenetic history. First, we quantified variation in pelage color (n=107 mice) and habitat color (n=51 rocks) using a spectrophotometer, and showed that there was a correlation between pelage color and habitat color across 14 sampled populations (R2=0.43). Analyses of mtDNA sequences from these same individuals revealed strong population structure in this species across its range, where most variation (63%) was partitioned between five geographic regions. Using Mantel tests, we show that there is no correlation between color variation and mtDNA phylogeny, suggesting that pelage coloration has evolved rapidly. At a finer geographical scale, high levels of gene flow between neighboring melanic and light populations suggest the selection acting on color must be quite strong to maintain habitat-specific phenotypic distributions. Finally, we raise the possibility that, in some cases, migration between populations of pocket mice inhabiting different lava flows may be responsible for similar melanic phenotypes in different populations. Together, the results suggest that color variation can evolve very rapidly over small geographic scales and that gene flow can both hinder and promote local adaptation.     

Evidence for morphological and adaptive genetic divergence between lake and stream habitats in European minnows (Phoxinus phoxinus, Cyprinidae)

Natural selection drives local adaptation, potentially even at small temporal and spatial scales. As a result, adaptive genetic and phenotypic divergence can occur among populations living in different habitats. We investigated patterns of differentiation between contrasting lake and stream habitats in the cyprinid fish European minnow (Phoxinus phoxinus) at both the morphological and genomic levels using geometric morphometrics and AFLP markers, respectively. We also used a spatial correlative approach to identify AFLP loci associated with environmental variables representing potential selective forces responsible for adaptation to divergent habitats. Our results identified different morphologies between lakes and streams, with lake fish presenting a deeper body and caudal peduncle compared to stream fish. Body shape variation conformed to a priori predictions concerning biomechanics and swimming performance in lakes vs. streams. Moreover, morphological differentiation was found to be associated with several environmental variables, which could impose selection on body and caudal peduncle shape. We found adaptive genetic divergence between these contrasting habitats in the form of 'outlier' loci (2.9%) whose genetic divergence exceeded neutral expectations. We also detected additional loci (6.6%) not associated with habitat type (lake vs. stream), but contributing to genetic divergence between populations. Specific environmental variables related to trophic dynamics, landscape topography and geography were associated with several neutral and outlier loci. These results provide new insights into the morphological divergence and genetic basis of adaptation to differentiated habitats.