Variable extent of parallelism in respiratory,circulatory, and neurological traits across lake whitefish species pairs

Parallel adaptive radiation events provide a powerful framework for investigations of ecology's contribution to phenotypic diversification. Ecologically driven divergence has been invoked to explain the repeated evolution of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis) species in multiple lakes in eastern North America. Nevertheless, links between most putatively adaptive traits and ecological variation remain poorly defined within and among whitefish species pairs. Here, we examine four species pairs for variation in gill, heart, and brain size; three traits predicted to show strong phenotypic responses to ecological divergence. In each of the species pairs, normals exhibited larger body size standardized gills compared to dwarfs – a pattern that is suggestive of a common ecological driver of gill size divergence. Within lakes, the seasonal hypoxia experienced in the benthic environment is a likely factor leading to the requirement for larger gills in normals. Interestingly, the morphological pathways used to achieve larger gills varied between species pairs from Québec and Maine, which may imply subtle non‐parallelism in gill size divergence related to differences in genetic background. There was also a non‐significant trend toward larger hearts in dwarfs, the more active species of the two, whereas brain size varied exclusively among the lake populations. Taken together, our results suggest that the diversification of whitefish has been driven by parallel and non‐parallel ecological conditions across lakes. Furthermore, the phenotypic response to ecological variation may depend on genetic background of each population.     

Parallelism in the oxygen transport system of the lake whitefish: the role of physiological divergence in ecological speciation

In North America, populations of lake whitefish (Coregonus clupeaformis) have evolved sympatric 'dwarf' and 'normal' ecotypes that are associated with distinct trophic niches within lakes. Trophic specialization should place diverging physiological demands on individuals, and thus, genes and phenotypes associated with energy production represent ideal candidates for studies of adaptation. Here, we test for the parallel divergence of traits involved in oxygen transport in dwarf and normal lake whitefish from Québec, Canada and Maine, USA. We observed significant differences in red blood cell morphology between the ecotypes. Specifically, dwarfs exhibited larger nuclei and a higher nucleus area/total cell area than normal whitefish in all of the lakes examined. In addition, isoelectric focusing gels revealed variation in the haemoglobin protein components found in whitefish. Dwarf and normal whitefish exhibited a similar number of protein components, but the composition of these components differed, with dwarf whitefish bearing a greater proportion of cathodic components compared to the normals. Furthermore, dwarf whitefish showed significant haemoglobin gene upregulation in the brain compared with the levels shown in normals. Together, our results indicate that metabolic traits involved in oxygen transport differ between the whitefish ecotypes and the strong parallel patterns of divergence observed across lakes implicates ecologically driven selection pressures. We discuss the function of these traits in relation to the differing trophic niches occupied by the whitefish and the potential contributions of trait plasticity and genetic divergence to energetic adaptation.     

Adaptive plasticity and genetic divergence in feeding efficiency during parallel adaptive radiation of whitefish (Coregonus spp.)

Parallel phenotypic divergence in replicated adaptive radiations could either result from parallel genetic divergence in response to similar divergent selection regimes or from equivalent phenotypically plastic response to the repeated occurrence of contrasting environments. In post‐glacial fish, replicated divergence in phenotypes along the benthic‐limnetic habitat axis is commonly observed. Here, we use two benthic‐limnetic species pairs of whitefish from two Swiss lakes, raised in a common garden design, with reciprocal food treatments in one species pair, to experimentally measure whether feeding efficiency on benthic prey has a genetic basis or whether it underlies phenotypic plasticity (or both). To do so, we offered experimental fish mosquito larvae, partially burried in sand, and measured multiple feeding efficiency variables. Our results reveal both, genetic divergence as well as phenotypically plastic divergence in feeding efficiency, with the phenotypically benthic species raised on benthic food being the most efficient forager on benthic prey. This indicates that both, divergent natural selection on genetically heritable traits and adaptive phenotypic plasticity, are likely important mechanisms driving phenotypic divergence in adaptive radiation.     

CORRELATED TROPHIC SPECIALIZATION AND GENETIC DIVERGENCE IN SYMPATRIC LAKE WHITEFISH ECOTYPES (COREGONUS CLUPEAFORMIS): SUPPORT FOR THE ECOLOGICAL SPECIATION HYPOTHESIS

There is ample empirical evidence that phenotypic diversification in an adaptive radiation is the outcome of divergent natural selection related to differential resource use. In contrast, the role of ecological forces in favoring and maintaining reproductive isolation in nature remains poorly understood. If the same forces driving phenotypic divergence are also responsible for speciation, one would predict a correlation between the extent of trophic specialization (reflecting variable intensity of divergent natural selection) and that of reproductive isolation being reached in a given environment. We tested this hypothesis by comparing the extent of morphological and genetic differentiation between sympatric dwarf and normal whitefish ecotypes (Coregonus sp.) from six lakes of the St. John River basin (eastern Canada and northern Maine). Eight meristic variables, 19 morphometric variables, and six microsatellite loci were used to quantify morphological and genetic differentiation, respectively. Dwarf and normal ecotypes in each lake differed primarily by traits related to trophic specialization, but the extent of differentiation varied among lakes. Significant but variable genetic divergence between ecotypes within lakes was also observed. A negative correlation was observed between the extent of gene flow between ecotypes within a lake and that of their morphological differentiation in trophic-related traits. The extent of reproductive isolation reached between dwarf and normal whitefish ecotypes appears to be driven by the potential for occupying distinct trophic niches and, thus, by the same selective forces driving tropic specialization in each lake. These results therefore support the hypothesis of ecological speciation.     

Diversifying selection drives parallel evolution of gill raker number and body size along the speciation continuum of European whitefish

Adaptive radiation is the evolution of ecological and phenotypical diversity. It arises via ecological opportunity that promotes the exploration of underutilized or novel niches mediating specialization and reproductive isolation. The assumed precondition for rapid local adaptation is diversifying natural selection, but random genetic drift could also be a major driver of this process. We used 27 populations of European whitefish (Coregonus lavaretus) from nine lakes distributed in three neighboring subarctic watercourses in northern Fennoscandia as a model to test the importance of random drift versus diversifying natural selection for parallel evolution of adaptive phenotypic traits. We contrasted variation for two key adaptive phenotypic traits correlated with resource utilization of polymorphic fish; the number of gill rakers and the total length of fish, with the posterior distribution of neutral genetic differentiation from 13 microsatellite loci, to test whether the observed phenotypic divergence could be achieved by random genetic drift alone. Our results show that both traits have been under diversifying selection and that the evolution of these morphs has been driven by isolation through habitat adaptations. We conclude that diversifying selection acting on gill raker number and body size has played a significant role in the ongoing adaptive radiation of European whitefish morphs in this region.     

Model‐based demographic inference of introgression history in European whitefish species pairs'

Parallel phenotypic differentiation is generally attributed to parallel adaptive divergence as an evolutionary response to similar environmental contrasts. Such parallelism may actually originate from several evolutionary scenarios ranging from repeated parallel divergence caused by divergent selection to a unique divergence event followed by gene flow. Reconstructing the evolutionary history underlying parallel phenotypic differentiation is thus fundamental to understand the relative contribution of demography and selection on genomic divergence during speciation. In this study, we investigate the divergence history of replicate European whitefish (Coregonus lavaretus), limnetic and benthic species pairs from two lakes in Norway and two lakes in Switzerland. Demographic models accounting for semi‐permeability and linked selection were fitted to the unfolded joint allele frequency spectrum built from genome‐wide SNPs and compared to each other in each species pair. We found strong support for a model of asymmetrical post‐glacial secondary contact between glacial lineages in all four lakes. Moreover, our results suggest that heterogeneous genomic differentiation has been shaped by the joint action of linked selection accelerating lineage sorting during allopatry, and heterogeneous migration eroding divergence at different rates along the genome following secondary contact. Our analyses reveal how the interplay between demography, selection and historical contingency has influenced the levels of diversity observed in previous whitefish phylogeographic studies. This study thus provides new insights into the historical demographic and selective processes that shaped the divergence associated with ecological speciation in European whitefish.     

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.     

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.     

Parallel evolution of lake whitefish dwarf ecotypes in association with limnological features of their adaptive landscape

Sympatric fish populations observed in many north temperate lakes are among the best models to study the processes of population divergence and adaptive radiation. Despite considerable research on such systems, little is known about the associations between ecological conditions and the extent of ecotypic divergence. In this study, we examined the biotic and abiotic properties of postglacial lakes in which lake whitefish, Coregonus clupeaformis, occur as a derived dwarf ecotype in sympatry with an ancestral normal ecotype. We compared 19 limnological variables between two groups of lakes known from previous studies to harbour sympatric dwarf and normal ecotypes with high and low levels of phenotypic and genetic differentiation respectively. We found clear environmental differences between the two lake groups. Namely, oxygen was the most discriminant variable, where lakes harbouring the most divergent populations were characterized by the greatest hypolimnetic oxygen depletion. These lakes also had lower zooplankton densities and a narrower distribution of zooplantonic prey length. These results suggest that the highest differentiation between sympatric ecotypes occurs in lakes with reduced habitat and prey availability that could increase competition for resources. This in turns supports the hypothesis that parallelism in the extent of phenotypic divergence among sympatric whitefish ecotypes is associated with parallelism in adaptive landscape in terms of differences in limnological characteristics, as well as availability and structure of the zooplanktonic community.     

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.     

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     

Parallelism and historical contingency during rapid ecotype divergence in an isopod

Recent studies on parallel evolution have focused on the relative role of selection and historical contingency during adaptive divergence. Here, we study geographically separate and genetically independent lake populations of a freshwater isopod (Asellus aquaticus) in southern Sweden. In two of these lakes, a novel habitat was rapidly colonized by isopods from a source habitat. Rapid phenotypic changes in pigmentation, size and sexual behaviour have occurred, presumably in response to different predatory regimes. We partitioned the phenotypic variation arising from habitat ('selection': 81–94%), lake ('history': 0.1–6%) and lake × habitat interaction ('unique diversification': 0.4–13%) for several traits. There was a limited role for historical contingency but a strong signature of selection. We also found higher phenotypic variation in the source populations. Phenotype sorting during colonization and strong divergent selection might have contributed to these rapid changes. Consequently, phenotypic divergence was only weakly influenced by historical contingency.     

Empirical evaluation of phenotype–environment correlation and trait utility with allopatric and sympatric whitefish, Coregonus lavaretus (L.), populations in subarctic lakes

Adaptive phenotypic divergence of sympatric morphs in a single species may have significant evolutionary consequences. In the present study, phenotypic impacts of predator on zooplankton prey populations were compared in two northern Finnish lakes; one with an allopatric whitefish, Coregonus lavaretus (L.), population and the other with three sympatric whitefish populations. First, we examined whether there were phenotypic associations with specific niches in allopatric and sympatric whitefish. Second, trait utility (i.e. number of gillrakers) of allopatric and sympatric whitefish in utilizing a pelagic resource was explored by comparing predator avoidance of prey, prey size in environment, and prey size in predator diet. The allopatric living large sparsely rakered (LSR) whitefish morph, was a generalist using both pelagic and benthic niches. In contrast, sympatric living whitefish morphs were specialized: LSR whitefish was a littoral benthivore, small sparsely rakered whitefish was a profundal benthivore and densely rakered (DR) whitefish was a pelagic planktivore. In the lake with allopatric whitefish, zooplankton prey did not migrate vertically to avoid predation whereas, in the lake with sympatric whitefish, all important prey taxa migrated significantly. Trait utility was observed as significantly smaller size of prey in environment and predator diet in the lake with DR whitefish than in the lake with only LSR whitefish.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 92 , 561–572.     

Relationships between lake ecology and morphological characters in Icelandic Arctic charr,Salvelinus alpinus

The common occurrence of parallel phenotypic patterns suggests that a strong relationship exists between ecological dynamics and micro‐evolution. Comparative studies from a large number of populations under varying sets of ecological drivers could contribute to a better understanding of this relationship. We used data on morphology of arctic charr (Salvelinus alpinus) and ecological factors from 35 Icelandic lakes to test the hypothesis that morphological patterns among monomorphic charr populations from different lakes are related to interlake variation in ecological characteristics. There is extensive phenotypic diversity among populations of Icelandic charr, and populations are easily distinguished based on overall body morphology. The results obtained in the present study showed that the morphological diversity of charr was related to large‐scale diversity in lake ecology. Variation in charr morphology was related to water origin (e.g. spring fed versus run‐off), bedrock age, and fish community structure. The present study shows how various ecological factors can shape the biological diversity that we observe. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103 , 761–771.     

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.     

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.