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.     

Ecosystem size matters: the dimensionality of intralacustrine diversification in Icelandic stickleback is predicted by lake size

Cases of evolutionary diversification can be characterized along a continuum from weak to strong genetic and phenotypic differentiation. Several factors may facilitate or constrain the differentiation process. Comparative analyses of replicates of the same taxon at different stages of differentiation can be useful to identify these factors. We estimated the number of distinct phenotypic groups in three‐spine stickleback populations from nine lakes in Iceland and in one marine population. Using the inferred number of phenotypic groups in each lake, genetic divergence from the marine population, and physical lake and landscape variables, we tested whether ecosystem size, approximated by lake size and depth, or isolation from the ancestral marine gene pool predicts the occurrence and the extent of phenotypic and genetic diversification within lakes. We find intralacustrine phenotypic diversification to be the rule rather than the exception, occurring in all but the youngest lake population and being manifest in ecologically important phenotypic traits. Neutral genetic data further indicate nonrandom mating in four of nine studied lakes, and restricted gene flow between sympatric phenotypic groups in two. Although neither the phenotypic variation nor the number of intralacustrine phenotypic groups was associated with any of our environmental variables, the number of phenotypic traits that were differentiated was significantly positively related to lake size, and evidence for restricted gene flow between sympatric phenotypic groups was only found in the largest lakes where trait specific phenotypic differentiation was highest.     

ADAPTIVE RADIATION AND ECOLOGICAL DIVERSIFICATION OF SULAWESI'S ANCIENT LAKE SHRIMPS

Ancient lakes are natural laboratories for the study of adaptive radiation. Recently, two lake systems on the Indonesian island Sulawesi have emerged as promising new model systems. A species flock of atyid freshwater shrimps in the Malili lake system comprises 15 colorful endemic taxa. Mitochondrial DNA data suggest two independent colonizations by riverine ancestors. Only one colonization event led to subsequent radiation into 14 species, while the second clade comprises just one species. The vast majority of species (n= 12) are habitat specialists, which are confined to the larger Malili clade and include all taxa with species‐specific color patterns and a restricted distribution within the five connected Malili lakes. Morphological, genetic, and ecological data are consistent with the existence of an adaptive radiation in the Malili lakes, involving the habitat‐specific diversification of trophic morphology. In addition to testing criteria for the recognition of an adaptive radiation, an ancestral state reconstruction reveals an equal probability for either a riverine generalist or a lacustrine specialist as ancestor of the large Malili clade, which is interpreted as indicative of an early stage of habitat specialization within this radiation. Finally, our results suggest that species diversification may have been primarily driven by ecological specialization and allopatric speciation.     

Speciation patterns and processes in the zooplankton of the ancient lakes of Sulawesi Island,Indonesia

Although studies of ancient lake fauna have provided important insights about speciation patterns and processes of organisms in heterogeneous benthic environments, evolutionary forces responsible for speciation in the relatively homogenous planktonic environment remain largely unexplored. In this study, we investigate possible mechanisms of speciation in zooplankton using the freshwater diaptomids of the ancient lakes of Sulawesi, Indonesia, as a model system. We integrate phylogenetic and population genetic analyses of mitochondrial and nuclear genes with morphological and genome size data. Overall, our results support the conclusion that colonization order and local adaptation are dominant at the large, island scale, whereas at local and intralacustrine scales, speciation processes are regulated by gene flow among genetically differentiated and locally adapted populations. In the Malili lakes, the diaptomid populations are homogenous at nuclear loci, but show two highly divergent mitochondrial clades that are geographically restricted to single lakes despite the interconnectivity of the lake systems. Our study, based on coalescent simulations and population genetic analyses, indicates that unidirectional hybridization allows gene flow across the nuclear genome, but prevents the introgression of mitochondria into downstream populations. We suggest that hybridization and introgression between young lineages is a significant evolutionary force in freshwater plankton.     

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.     

Ecosystem size predicts eco‐morphological variability in a postglacial diversification

Identifying the processes by which new phenotypes and species emerge has been a long‐standing effort in evolutionary biology. Young adaptive radiations provide a model to study patterns of morphological and ecological diversification in environmental context. Here, we use the recent radiation (ca. 12k years old) of the freshwater fish Arctic charr (Salvelinus alpinus) to identify abiotic and biotic environmental factors associated with adaptive morphological variation. Arctic charr are exceptionally diverse, and in postglacial lakes there is strong evidence of repeated parallel evolution of similar morphologies associated with foraging. We measured head depth (a trait reflecting general eco‐morphology and foraging ecology) of 1,091 individuals across 30 lake populations to test whether fish morphological variation was associated with lake bathymetry and/or ecological parameters. Across populations, we found a significant relationship between the variation in head depth of the charr and abiotic environmental characteristics: positively with ecosystem size (i.e., lake volume, surface area, depth) and negatively with the amount of littoral zone. In addition, extremely robust‐headed phenotypes tended to be associated with larger and deeper lakes. We identified no influence of co‐existing biotic community on Arctic charr trophic morphology. This study evidences the role of the extrinsic environment as a facilitator of rapid eco‐morphological diversification.     

Novel trophic niches drive variable progress towards ecological speciation within an adaptive radiation of pupfishes

Adaptive radiation is recognized by a rapid burst of phenotypic, ecological and species diversification. However, it is unknown whether different species within an adaptive radiation evolve reproductive isolation at different rates. We compared patterns of genetic differentiation between nascent species within an adaptive radiation of Cyprinodon pupfishes using genotyping by sequencing. Similar to classic adaptive radiations, this clade exhibits rapid morphological diversification rates and two species are novel trophic specialists, a scale‐eater and hard‐shelled prey specialist (durophage), yet the radiation is <10 000 years old. Both specialists and an abundant generalist species all coexist in the benthic zone of lakes on San Salvador Island, Bahamas. Based on 13 912 single‐nucleotide polymorphisms (SNPs), we found consistent differences in genetic differentiation between each specialist species and the generalist across seven lakes. The scale‐eater showed the greatest genetic differentiation and clustered by species across lakes, whereas durophage populations often clustered with sympatric generalist populations, consistent with parallel speciation across lakes. However, we found strong evidence of admixture between durophage populations in different lakes, supporting a single origin of this species and genome‐wide introgression with sympatric generalist populations. We conclude that the scale‐eater is further along the speciation‐with‐gene‐flow continuum than the durophage and suggest that different adaptive landscapes underlying these two niche environments drive variable progress towards speciation within the same habitat. Our previous measurements of fitness surfaces in these lakes support this conclusion: the scale‐eating fitness peak may be more distant than the durophage peak on the complex adaptive landscape driving adaptive radiation.     

Ecological speciation in postglacial European whitefish: rapid adaptive radiations into the littoral,pelagic, and profundal lake habitats

Understanding how a monophyletic lineage of a species diverges into several adaptive forms has received increased attention in recent years, but the underlying mechanisms in this process are still under debate. Postglacial fishes are excellent model organisms for exploring this process, especially the initial stages of ecological speciation, as postglacial lakes represent replicated discrete environments with variation in available niches. Here, we combine data of niche utilization, trophic morphology, and 17 microsatellite loci to investigate the diversification process of three sympatric European whitefish morphs from three northern Fennoscandian lakes. The morphological divergence in the gill raker number among the whitefish morphs was related to the utilization of different trophic niches and was associated with reproductive isolation within and across lakes. The intralacustrine comparison of whitefish morphs showed that these systems represent two levels of adaptive divergence: (1) a consistent littoral–pelagic resource axis; and (2) a more variable littoral–profundal resource axis. The results also indicate that the profundal whitefish morph has diverged repeatedly from the ancestral littoral whitefish morph in sympatry in two different watercourses. In contrast, all the analyses performed revealed clustering of the pelagic whitefish morphs across lakes suggesting parallel postglacial immigration with the littoral whitefish morph into each lake. Finally, the analyses strongly suggested that the trophic adaptive trait, number of gill rakers, was under diversifying selection in the different whitefish morphs. Together, the results support a complex evolutionary scenario where ecological speciation acts, but where both allopatric (colonization history) and sympatric (within watercourse divergence) processes are involved.     

Ancient lakes revisited: from the ecology to the genetics of speciation

Explosive speciation in ancient lakes has fascinated biologists for centuries and has inspired classical work on the tempo and modes of speciation. Considerable attention has been directed towards the extrinsic forces of speciation—the geological, geographical and ecological peculiarities of ancient lakes. Recently, there has been a resurgence of interest in the intrinsic nature of these radiations, the biological characteristics conducive to speciation. While new species are thought to arise mainly by the gradual enhancement of reproductive isolation among geographically isolated populations, ancient lakes provide little evidence for a predominant role of geography in speciation. Recent phylogenetic work provides strong evidence that multiple colonization waves were followed by parallel intralacustrine radiations that proceeded at relatively rapid rates despite long‐term gene flow through hybridization and introgression. Several studies suggest that hybridization itself might act as a key evolutionary mechanism by triggering major genomic reorganization/revolution and enabling the colonization of new ecological niches in ancient lakes. These studies propose that hybridization is not only of little impediment to diversification but could act as an important force in facilitating habitat transitions, promoting postcolonization adaptations and accelerating diversification. Emerging ecological genomic approaches are beginning to shed light on the long‐standing evolutionary dilemma of speciation in the face of gene flow. We propose an integrative programme for future studies on speciation in ancient lakes.     

Role of ancient lakes in genetic and phenotypic diversification of freshwater snails

Endemic organisms of ancient lakes have been studied as models to understand processes of speciation and adaptive radiation. However, it remains unclear how ancient lakes play roles in genetic and phenotypic diversity of freshwater mollusks. In the present study, we focus on viviparid freshwater snails in the ancient lakes of East and Southeast Asia (Japan and China) to address this question. Using molecular phylogenetic analyses based on mitochondrial (COI, 16S) and nuclear genes (18S, 28S, H3), we show that patterns of species diversification in viviparid lineages. Colonization to ancient lakes occurred independently in China and Japan at least four times, with subsequent diversification into more than two species within each lake group. Morphological analyses of fossil related viviparids suggest parallel phenotypic evolution occurred in the different lakes and ages. Each lake contained a single lineage which was phenotypically diversified relative to those from other sites. Using genome‐wide SNPs obtained by MIG‐seq, we also examined the genetic structure of three Japanese viviparids, including two endemic species of ancient Lake Biwa. The results suggest that these two species diversified from the population of the third species living in wetlands surrounding the lake. These findings suggest that rapid diversification of lineages and phenotypic divergence can occur in ancient lakes compared to other habitats. Formation of large lakes probably promotes speciation and phenotypic divergence as a result of adaptation into different microhabitats. High numbers of ancient lakes could be a driver of species diversity in Asian viviparid snails.     

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.     

Genomic signatures of divergent selection and speciation patterns in a ‘natural experiment’, the young parallel radiations of Nicaraguan crater lake cichlid fishes

Divergent selection is the main driving force in sympatric ecological speciation and may also play a strong role in divergence between allopatric populations. Characterizing the genome‐wide impact of divergent selection often constitutes a first step in unravelling the genetic bases underlying adaptation and ecological speciation. The Midas cichlid fish (Amphilophus citrinellus) species complex in Nicaragua is a powerful system for studying evolutionary processes. Independent colonizations of isolated young crater lakes by Midas cichlid populations from the older and great lakes of Nicaragua resulted in the repeated evolution of adaptive radiations by intralacustrine sympatric speciation. In this study we performed genome scans on two repeated radiations of crater lake species and their great lake source populations (1030 polymorphic AFLPs, n ～ 30 individuals per species). We detected regions under divergent selection (0.3% in the crater lake Xiloá flock and 1.7% in the older crater lake Apoyo radiation) that might be responsible for the sympatric diversifications. We find no evidence that the same genomic regions have been involved in the repeated evolution of parallel adaptations across crater lake flocks. However, there is some genetic parallelism apparent (seven out of 51 crater lake to great lake outlier loci are shared; 13.7%) that is associated with the allopatric divergence of both crater lake flocks. Interestingly, our results suggest that the number of outlier loci involved in sympatric and allopatric divergence increases over time. A phylogeny based on the AFLP data clearly supports the monophyly of both crater lake species flocks and indicates a parallel branching order with a primary split along the limnetic‐benthic axis in both radiations.     

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.     

Fitness differences between parapatric lake and stream stickleback revealed by a field transplant

Molecular comparisons of populations diverging into ecologically different environments often reveal strong differentiation in localized genomic regions, with the remainder of the genome being weakly differentiated. This pattern of heterogeneous genomic divergence, however, is rarely connected to direct measurements of fitness differences among populations. We here do so by performing a field enclosure experiment in threespine stickleback fish residing in a lake and in three replicate adjoining streams, and displaying weak yet heterogeneous genomic divergence between these habitats. Tracking survival over 29 weeks, we consistently find that lake genotypes transplanted into the streams suffer greatly reduced viability relative to local stream genotypes and that the performance of F1 hybrid genotypes is intermediate. This observed selection against migrants and hybrids combines to a total reduction in gene flow from the lake into streams of around 80%. Overall, our study identifies a strong reproductive barrier between parapatric stickleback populations, and cautions against inferring weak fitness differences between populations exhibiting weak overall genomic differentiation.     

Ecological opportunity alters the timing and shape of adaptive radiation

The uneven distribution of diversity is a conspicuous phenomenon across the tree of life. Ecological opportunity is a prominent catalyst of adaptive radiation and therefore may alter patterns of diversification. We evaluated the distribution of shifts in diversification rates across the cichlid phylogeny and the distribution of major clades across phylogenetic space. We also tested if ecological opportunity influenced these patterns. Colonization‐associated ecological opportunity altered the tempo and mode of diversification during the adaptive radiation of cichlid fishes. Clades that arose following colonization events diversified faster than other clades. Speciation rate shifts were nonrandomly distributed across the phylogeny such that they were disproportionally concentrated around nodes that corresponded with colonization events (i.e., of continents, river basins, or lakes). Young clades tend to expand faster than older clades; however, colonization‐associated ecological opportunity accentuated this pattern. There was an interaction between clade age and ecological opportunity that explained the trajectory of clades through phylogenetic space over time. Our results indicate that ecological opportunities afforded by continental and ecosystem‐scale colonization events explain the dramatic speciation rate heterogeneity and phylogenetic imbalance that arose during the evolutionary history of cichlid fishes.     

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.     

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.     

Effects of population,family, and diet on craniofacial morphology of Icelandic Arctic charr (Salvelinus alpinus)

We evaluated hypotheses of intralacustrine diversification and plastic responses to two diet environments in Icelandic Arctic charr (Salvelinus alpinus). Full‐sib families of progeny of wild polymorphic charr from two lakes where morphs vary in their degree of phenotypic and ecological divergence were split, with half of the offspring reared on a benthic and half on a limnetic type of diet to estimate family norms of reaction. We focused on variation in craniofacial traits because they are probably functionally related to diet and complement a previous study of body shape in these charr. A hierarchical analysis of phenotypic variation between lakes, pairs of morphs within each lake, and two families within each morph found that phenotypic variation partitioned between families relative to morphs was reduced in the more ecologically diversified population, which is consistent with adaptive diversification. The effect size of plastic responses between lake populations was similar, suggesting little difference in the degree of canalization in contrast to a previous analysis of body form plasticity. Thus, the role that plastic morphological responses play in the adaptive diversification of morphs and different lake populations of Arctic charr may depend on the trait. © 2013 The Linnean Society of London     

Demography and genome divergence of lake and stream populations of an East African cichlid fish

Disentangling the processes and mechanisms underlying adaptive diversification is facilitated by the comparative study of replicate population pairs that have diverged along a similar environmental gradient. Such a setting is realized in a cichlid fish from southern Lake Tanganyika, Astatotilapia burtoni, which occurs within the lake proper as well as in various affluent rivers. Previously, we demonstrated that independent lake and stream populations show similar adaptations to the two habitat regimes. However, little is known about the evolutionary and demographic history of the A. burtoni populations in question and the patterns of genome divergence among them. Here, we apply restriction site‐associated DNA sequencing (RADseq) to examine the evolutionary history, the population structure and genomic differentiation of lake and stream populations in A. burtoni. A phylogenetic reconstruction based on genome‐wide molecular data largely resolved the evolutionary relationships among populations, allowing us to re‐evaluate the independence of replicate lake–stream population clusters. Further, we detected a strong pattern of isolation by distance, with baseline genomic divergence increasing with geographic distance and decreasing with the level of gene flow between lake and stream populations. Genome divergence patterns were heterogeneous and inconsistent among lake‐stream population clusters, which is explained by differences in divergence times, levels of gene flow and local selection regimes. In line with the latter, we only detected consistent outlier loci when the most divergent lake–stream population pair was excluded. Several of the thus identified candidate genes have inferred functions in immune and neuronal systems and show differences in gene expression between lake and stream populations.