Adaptation and acclimation of traits associated with swimming capacity in Lake Whitefish (coregonus clupeaformis) ecotypes
© The Author(s). 2016
Received: 9 March 2016
Accepted: 1 August 2016
Published: 11 August 2016
Improved performance in a given ecological niche can occur through local adaptation, phenotypic plasticity, or a combination of these mechanisms. Evaluating the relative importance of these two mechanisms is needed to better understand the cause of intra specific polymorphism. In this study, we reared populations of Lake Whitefish (Coregonus clupeaformis) representing the’normal’ (benthic form) and the ‘dwarf’ (derived limnetic form) ecotypes in two different conditions (control and swim-training) to test the relative importance of adaptation and acclimation in the differentiation of traits related to swimming capacity. The dwarf whitefish is a more active swimmer than the normal ecotype, and also has a higher capacity for aerobic energy production in the swimming musculature. We hypothesized that dwarf fish would show changes in morphological and physiological traits consistent with reductions in the energetic costs of swimming and maintenance metabolism.
We found differences in traits predicted to decrease the costs of prolonged swimming and standard metabolic rate and allow for a more active lifestyle in dwarf whitefish. Dwarf whitefish evolved a more streamlined body shape, predicted to lead to a decreased drag, and a smaller brain, which may decrease their standard metabolic rate. Contrary to predictions, we also found evidence of acclimation in liver size and metabolic enzyme activities.
Results support the view that local adaptation has contributed to the genetically-based divergence of traits associated with swimming activity. Presence of post-zygotic barriers limiting gene flow between these ecotype pairs may have favoured repeated local adaptation to the limnetic niches.
Local adaptation occurs when individuals within a population evolve in response to selective pressure leading to increased fitness in their local environment relative to a foreign environment . It is notoriously difficult to differentiate between local adaptation and environmentally induced plasticity in wild populations because phenotypic variation could be the product of either mechanism, or a combination of genetically based differentiation and plastic responses. Distinguishing between these two non-exclusive mechanisms is needed to better understand the cause of intra specific phenotypic polymorphism and the interaction of organisms with their environment [2–4]. This is particularly true for traits that commonly respond to a gradual change in the environment, a process known as acclimation (i.e. reversible, “physiological” phenotypic plasticity, or flexibility; ). Because it often increases survival in new environments, the process of acclimation is predicted to favour gene flow [6, 7], and hinder local adaptation . Therefore, it is hypothesized that the extent of local adaptation should be inversely proportional to the level of gene flow between sympatric morphs exploiting different niches.
Northern freshwater fishes inhabiting postglacial environments show remarkable amounts of intra specific polymorphism [8–12] and numerous cases of parallel phenotypic evolution [13–17]. After the retreat of the Pleistocene ice sheets (~12,000 years ago), several species colonised newly formed lakes [18, 19], with free ecological niches allowing for phenotypic diversification [20, 21]. More specifically, sympatric ecotypes often show patterns of phenotypic differentiation associated with the benthic and limnetic niches [10, 11, 22–24]. A major difference between benthic and limnetic environments is the need to continually swim when foraging in the limnetic niche [24–27], as limnetic prey (i.e. zooplankton) show greater variation in abundance and distribution than benthic prey [24, 28]. The energy that is spent on swimming activity takes away from the energy available for other metabolic requirements such as growth, tissue maintenance and reproduction [29–31]. To maintain a balance between energy supply and demand and cope with the requirements for more active swimming in the limnetic niche, fish can adopt a number of strategies [25–27, 29, 30, 32]. These include decreasing the costs of swimming or the standard metabolic rate (i.e. the minimum metabolic rate needed to sustain life when animals are post-absorptive and inactive) [27, 31–33]. For instance, a common feature observed in limnetic morphs is a more hydrodynamic, slender body shape, which will diminish drag during steady swimming [24, 34, 35]. In contrast, benthic morphs often display deeper bodies, which enhance manoeuvrability and burst swimming capacity and increase foraging efficiency on benthic prey [24, 32, 33, 36]. Fish may also decrease the costs of maintenance metabolism by reducing the size of ‘metabolically-expensive’ organs. Indeed, studies of intra specific variation in standard metabolic rate often find positive correlations between standard/basal metabolic rate and the size of metabolically-expensive organs and a negative correlation with growth [31–33]. Brain and liver are two of the most metabolically-costly organs, so decreasing their size is predicted to reduce whole-animal metabolic maintenance cost [37–40]. A reduction of total gill surface area may also decrease whole organism energy expenditure because a larger gill surface area will lead to greater ion loss in freshwater, requiring fish to expend energy pumping ions to maintain osmoregulatory homeostasis . These reductions in overall gill surface area may decrease oxygen uptake and hinder aerobic energy metabolism. However, fish normally use all of their gill area for oxygen uptake only during maximum aerobic swimming , which is not the case during limnetic foraging. Therefore, a decrease in gill surface area is predicted to be selected in freshwater environments when fish are energy limited, but oxygen is not limiting (i.e. normoxic waters). Overall, fish living in the limnetic niche are predicted to have more streamlined bodies, with smaller organ sizes (e.g. liver, brain, gills) compared to their benthic congeners, and thus save energy on maintenance costs and swimming which will allow for more energy to be spent on reproduction and growth (mainly in the white skeletal muscle) while actively foraging [29, 31].
Two Lake Whitefish (Coregonus clupeaformis) ecotypes are found in several postglacial lakes of the St. John River Drainage (Québec, Canada and Maine, USA) . The ‘normal’ ecotype (Atlantic lineage) is specialized to a benthic niche while the repeatedly derived ‘dwarf’ ecotype (Acadian lineage) has colonised the free limnetic niche following secondary contact after the last glacial maximum [lineages were divided ~ 60,000 years ago or about 15,000 generations during the last Pleistocene glaciation and the secondary contact occurred ~ 12,000 years ago or about 3000 generations] [23, 44–47]. Reduced viability and sperm performance, reactivation of transposable elements and aneuploidy in hybrids are post-zygotic barriers that restrict gene flow between these ecotypes [48–52], in addition to ecological processes that may reduce the fitness of hybrids in natural conditions [11, 49]. As predicted by their ecological divergence, the limnetic ‘dwarf’ whitefish are more active swimmers and are younger at maturity (2–3 years vs 5 years or more) than normal fish, whereas normal whitefish grow more quickly, live longer, attain larger sizes, and have a higher condition factor than dwarf fish [23, 53, 54]. Studies examining the mechanisms underlying these differences in whole-animal performance have revealed that wild dwarf fish are more streamlined  and have a higher capacity for oxygen transport and use, which is required to fuel aerobic swimming (e.g. larger hearts, more red muscle, higher muscle mitochondrial content) [56–58]. Some traits related to maintenance and growth also follow predictions, as gill surface area in dwarf fish is smaller than in normal fish, but no differences in brain size are detected in wild caught fish . Differences in gene expression between limnetic and benthic whitefish ecotypes have found extensive divergence in central energy metabolism (i.e. mitochondrial oxidative phosphorylation, citric acid cycle, glycolysis, glycogen metabolism, creatine phosphokinase). This metabolic divergence is predicted to underlie the observed trade-off in life-history traits, wherein enhanced survival via more active swimming is necessary for increased foraging and predator avoidance in the limnetic zone, but leads to energetic costs that translate into slower growth rates and reduced fecundity in dwarf relative to normal whitefish [23, 60–64].
However, the relative roles of acclimation and adaptation in the expression of these traits, including brain size, gill size and body shape is unclear because previous studies have all been performed on wild caught fish where dwarf and normal whitefish may be exposed to different environmental conditions [55, 59]. Moreover, juvenile and adult fish often show phenotypic plasticity in these traits in similar conditions [23, 65–67]. Therefore, the goal of this study is to test if variation in traits is due to acclimation and/or adaptation in controlled environments. We predicted that both mechanisms are used to compensate for the increase in energy allocation to locomotion in dwarf fish, but that adaptation should be a major mechanism underlying phenotypic divergence, considering the reproductive barriers previously documented between ecotypes pairs [11, 48–52]. We also further investigate the energetic strategies used by dwarf fish to inhabit the limnetic niche by measuring additional traits predicted to influence maintenance metabolism (i.e. liver size, liver metabolic capacity and a more detailed study of gill surface area).
Crosses were produced from parents collected from Témiscouata (dwarf whitefish, Acadian lineage 47°36 N, 68°45 W) and Aylmer lakes (normal whitefish, Atlantic lineage, 45°50 N, 71°26 W) in Québec, Canada, caught in November 2011. It should be noted that these two populations are the only populations of whitefish ecotypes for which the breeding sites are known and crosses can be produced from. However, the divergence between these two allopatric ‘dwarf’ and ‘normal’ ecotypes is representative of the differentiation among other wild sympatric ecotype pairs. This is because transcriptomic comparisons among ecotypes show similar patterns of differentiation in these two populations and wild sympatric pairs of ecotypes . As well, dwarf individuals from Témiscouata and normal individuals from Alymer lakes raised in the lab show heritable differentiation in behavioural traits in accordance with the different life history traits expected between limnetic and benthic fishes  and show the same postzygotic isolation mechanisms described above.
After the images were captured, wet brain and liver mass were determined and gills were preserved in 4 % formaldehyde. The first arch of the left gill was laid flat on a microscope slide and photographed. We counted the number of filaments and estimated the average length of filaments, total length of filaments, average space between lamellae, number of lamellae and total hemibranch area using ImageJ . The two latter estimates were our gill area proxies. More precisely, the average filament length was estimated by measuring each fifth filament and divided the total by the number of filaments measured (Fig. 2b). These points were also used to obtain shape and measure of hemibranch area (Fig. 2b). Total filament length was obtained by multiplying the average filament length and number of filaments. Average space between lamellae was estimated by the distance between ten lamellae, and measured five times on different parts of the gill to control for measurement error (Fig. 2c). Finally, the total number of lamellae was calculated by dividing the total filament length by the average space between lamellae and multiplied by ten (since average space between lamella was measured on 10 lamellae).
To study the metabolic capacity of the liver, we measured the activity of COX (cytochrome c oxidase, EC 184.108.40.206, complex IV in the electron transport chain, and found on the inner mitochondrial membrane) and CS (citrate synthase, EC 220.127.116.11, a citric acid cycle enzyme found in the mitochondrial matrix) per gram of liver tissue and calculated enzyme activities for the entire liver by multiplying this by total liver mass. COX and CS activities are indicators of mitochondrial content and should reflect the metabolic capacity of the liver . To do this, we re-weighed frozen liver (g), and immediately added 20 volumes of chilled homogenization buffer (50mmoll −1 hepes, 1mmoll −1 EDTA and 0.1 % Triton X-100; pH7.4) in 4 ml Wheaton glass homogenizers kept on ice. Enzyme activities were measured on whole-cell extracts at 25 °C using non-limiting substrate concentrations outlined by .
We first tested for a relationship between each trait and body mass and used residuals to remove allometric effects when there was a significant correlation (P < 0.05). To test the effect of ecotype (fixed effect, dwarf or normal whitefish) and treatment (fixed effect, control or swim) we ran mixed effects linear models using the ‘nlme’ package in R v 3.1.1 with tanks nested as a random effect and individual fish nested within tanks (Two-way nested ANOVAs). Acclimation was defined as a significant fixed-effect of treatment (e.g. both ecotypes display a similar direction and extent of plasticity) or a significant interaction effect (e.g. only one ecotype is plastic or ecotypes display different directions and/or extents of plasticity) and a genetic control was defined by similar criteria (fixed effect of ‘ecotypes’ or a significant interaction indicating evolutionary variation in plasticity). We also performed a multivariate analysis to determine if traits measured in this study could differentiate our four experimental groups (i.e. ecotype X treatment: dwarf/control - dwarf/swim - normal/control - normal/swim) and determine which traits best differentiated them. To do this we performed a linear discriminant function analysis (DFA) on mean value specific to each ecotypes for each tank. All variables were standardized and the function ‘lda’ in R software v 3.1.1 was used to perform the DFA.
The estimated variance between the two body shape measures for the same individual, measured from two independent pictures taken after repositioning the fish, was on average 2.0 % of the variance between any two different individuals. This indicates that the variation among individuals is 50 times higher than the variation produced between two pictures of the same individual (i.e. measurement error). The differences among individuals explained a significant portion of the shape variation (MANOVA; d.f. = 126; approx. F = 4.43; P < 0.001). Therefore, the shape differentiation observed among experimental groups (if any) should not have been produced by fish manipulation and/or landmark positioning.
We found a significant effect of ecotype for eight traits that are predicted to mediate energetic trade-offs between dwarf and normal whitefish by contributing to reductions in the cost of prolonged swimming (body streamlining) and standard metabolic rate (liver, brain and gill size). Three of these eight traits (all linked to the liver: liver mass, total liver COX activity, total liver CS activity) also showed a significant interaction between ecotype and environment, but no treatment effect was observed. Based on ecological and functional expectations, only body shape and brain mass match our predictions for reductions in swimming and metabolic maintenance costs in the more active dwarf whitefish, while liver mass varied in the opposite direction to that which was predicted. In accordance with the presence of postzygotic barriers between these ecotypes [11, 48–52], our results indicate that adaptation (more than acclimation) is the major mechanism allowing specialisation to the limnetic niche. Furthermore, fish shape was the trait that best differentiated ecotypes in the multivariate analysis, suggesting that the evolution of fish shape is of major importance in the colonisation of a new trophic niche which requires a more active lifestyle. In contrast, traits associated with liver metabolism and gill size do not match our predictions. This could be the result of either genetic or functional constraints (i.e. limitations on how traits could evolve due to functional roles in other performance traits/biochemical pathways that strongly select against change) or stochastic evolutionary processes (e.g. genetic drift). Finally, the three traits associated with liver metabolism showing ecotype–environment interactions suggest that liver energy metabolism is more responsive to swimming activity than the other traits measured in this study.
Local adaptation vs. Acclimation in Lake Whitefish ecotypes
Three traits had significant ecotype–environment interactions and all were associated with liver metabolism (liver mass, total COX activity and total CS activity). In comparison, the genetic and environmental bases of traits related to oxygen uptake, transport and use in the same fish were also tested . Eight traits related to ventricle size and metabolism and white muscle mitochondrial respiration displayed evidence of acclimation for at least one ecotype. In contrast, twice (n = 15) as many traits, including hematocrit, skeletal muscle mitochondrial function and the activities of multiples enzyme displayed a strong genetic basis and limited acclimation. This suggests that acclimation to high swimming activity leads to changes in liver metabolism, ventricle size and mitochondrial respiration, but that traits related to body shape, brain size, gill size, hematocrit and skeletal muscle mitochondrial content are mainly under genetic control. Similar experiments in other salmonids commonly find a significant effect of exercise on morphological and physiological traits (e.g. heart size and mitochondrial enzyme activities, skeletal muscle metabolic enzyme activities and body shape) [34, 65, 83–85]. Here, we chose a water velocity between one and two body lengths per second to emulate the constant aerobic swimming required for the exploitation of the limnetic niche (zooplankton are more aggregate and dispersed than benthic prey in lakes) . This swimming speed should emulate the cruising speeds of many fish species in the wild . As well, our experimental fish did not respond well to higher water velocities, so we did not increase swimming speeds to the levels that were used in other studies with salmonid fishes . Our fish were swim-trained for six months, which is normally sufficient to induce a plastic response [34, 65, 85]. For instance, a significant effect of swim-training on body shape after training at a similar water velocity was found in the benthic fish Salaria fluviatilis after only 28 days . Together, these data suggest that Lake Whitefish acclimation to swim-training is relatively subtle in comparison to evolutionary divergence. This hypothesis is also supported by the first axis of the discriminant function analysis that accounts for 98.2 % of the variation among groups and strongly differentiates ecotypes but not treatments. Lake Whitefish ecotypes are known to have the strongest reproductive barrier among all studied species-pairs of northern temperate fishes undergoing ecological speciation , and our results are consistent with the hypothesis that strong reproductive barriers are linked to the accumulation of adaptive genetic divergence between sympatric ecotypes [87, 88]. Admittedly, our results do not exclude the possibility that other environmental variables (e.g. diet, predation, temperature, dissolved oxygen content), which vary between benthic and limnetic niches could induce a plastic response on these traits, nor can we rule out the presence of developmental plasticity at an earlier life stage.
Limits of the experimental setup
Comparing organisms reared in different controlled environments allow biologists to test for the presence of genetic differentiation and local adaptation . Differentiation observed in a trait of interest among environmental treatments indicates phenotypic plasticity in this trait, differentiation in a trait among ecotypes in the same environment indicates genetic divergence, and differential responses to environments among ecotypes suggests that phenotypic plasticity has evolved differently . To rigorously test for genetic differentiation, second-generation offspring reared in a common garden must be used in the controlled experiment, since long-term phenotypic plasticity produced by prior parental exposure to different environments is possible [3, 89]. In the present study, parents could not be chosen from grandparents reared in a common garden because of the long generation time of the Lake Whitefish (2–3 years for the dwarf and 5 years or more for the normal whitefish) and the difficulty of maintaining enough individuals for a second round of artificial fertilization without losing a substantial amount of genetic diversity. Despite the fact that long term environmental effects cannot be excluded, we clearly ruled out short term acclimation as a major mechanism leading to traits variation. In addition, previous studies comparing reciprocal hybrid crosses (dwarf female X normal male; normal female X dwarf male) found no significant differences in swimming activity among crosses . These data suggest that the parent of origin does not have a significant effect on these performance traits and that, if present, parental effects are equally transmitted by males and females . Finally, the findings from this study extend mainly to the comparison of this particular set of dwarf and normal whitefish.
Decreasing the cost of swimming
When a species begins to forage on a new resource, a first and crucial step in successful foraging is having the ability to reach this new resource, which comes at an energetic cost of transporting the body to the resource’s location [86, 90, 91]. Fish shape influences locomotion, and ultimately foraging efficiency. Thus, it is predicted to be under selection when major changes in diet occur that require a change in swimming mode [25–27, 32]. We found that dwarf whitefish have evolved a more slender body with smaller pectoral fins. Such shape differentiation is generally associated with the limnetic niche in salmoniform swimmers as it decreases drag, and allows fish to expend less energy foraging on dispersed zooplankton [25–28, 32]. In contrast, a stouter body and longer pectoral fins are expected to increase manoeuvrability and are predicted to be beneficial when foraging on benthic organisms [24–27, 32]. Fish shape was the trait that best discriminated our four experimental groups (five times more than any other trait), suggesting the importance of shape evolution during the colonisation of the limnetic niche. Because no effect of treatment was observed on body shape, the use of the benthic niche by a dwarf individual or the limnetic niche by a normal individual would likely result in a decrease in foraging efficiency. Such adaptive phenotypic divergence in foraging efficiency has been observed in the European Whitefish ecotypes as well . As discussed above, a long term environmental effect, such as parental effects, could not be ruled out by our experimental setup. However, similar shape differentiation in wild fish and a signal of selection on shape-QTL (Quantitative Trait Loci) have previously been documented, indicating that shape is genetically based and reflects adaptive divergence .
Diminishing costs of organ maintenance metabolism
With respect to maintenance metabolism, including the costs of tissue maintenance, reproduction and growth, we predicted that dwarf whitefish would show reductions in liver and brain size as both are metabolically expensive organs to maintain [37–39]. Reductions in the size of these organs may allow for an increase in energy available for foraging in the limnetic niche. As predicted, dwarf whitefish had a smaller brain mass than the normal ecotype, suggesting that a reduction in energy expenditure could be used to divert additional energy towards locomotion. However, no parallel trend in brain mass differentiation was observed among wild ecotypes of Lake Whitefish . Nevertheless, a bigger or a smaller brain in dwarf whitefish could be observed when considering a single sympatric pair . It is possible that other ecological conditions induce a plastic response in brain mass in wild Lake Whitefish [66, 93–95]. For example, changes in brain structures were found in fish actively escaping predators [66, 93] and it is thought that the smaller dwarf whitefish suffer from higher predation pressure than normal fish [46, 53, 54, 96]. Further research specifically addressing brain structure differentiation will be needed to unravel this organ’s role in the origin of the dwarf ecotype.
Contrary to our predictions, we found that the liver was larger in dwarf whitefish. We further investigated this finding by testing if the maintenance costs of liver tissue might vary among ecotypes. A higher mitochondrial content is predicted to be associated with a higher standard metabolic rate because of increases in membrane density and transmembrane ion gradient maintenance costs [97, 98]. Indeed, variation in liver mitochondrial content (COX and CS enzyme activity) correlates positively with differences in whole-fish standard metabolic rate among species, while total liver size does not . As predicted, we found that the enzyme activities of COX and CS were higher in normal whitefish per gram of liver indicating a higher mitochondrial content per gram of tissue than for dwarf fish. However, total liver enzyme activity, which is a function of both activity and liver size, remained higher in dwarf than normal fish under control conditions, opposite to predictions. Interestingly, there was an interaction between ecotype and environment for both enzymes such that when fish were subjected to swim-training, differences in overall liver enzyme activity were no longer found. This suggests that potential reduction in liver size in dwarf whitefish is constrained by the many other physiological functions this organ is involved in . For example, the higher food consumption rate of the dwarf whitefish could necessitate a relatively larger liver to support increases in gluconeogenesis, blood filtration and waste management [30, 61].
We also predicted dwarf fish to have a smaller total gill surface area because a reduction in ion loss in these freshwater fish could decrease whole organism energy expenditure [41, 42]. We observed an ecotype effect that revealed fewer gill filaments and smaller space between lamellae in the dwarf whitefish. Fewer gill filaments matches our prediction but smaller space between lamellae does not. Since no significant differentiation was observed in the number of lamellae, the variation in filament number and lamellar spacing negate each other leading to similar gill surface areas. In contrast, a nearly significant ecotype-environment interaction in hemibranch area was also observed. However, this trait does not take into account differences in space between lamellae when estimating gill surface area. The similar total gill surface areas in dwarf and normal fish suggests that differences in underlying gill morphological traits (e.g. lamellar spacing and number of filaments) is likely due to stochastic evolutionary processes . In wild fish, there were differences in hemibranch area, hemibranch perimeter, number of filaments, average length of filaments and total length of filaments, leading to a slightly, but not significantly, larger gill surface area in normal whitefish, but space between lamellae was not measured . Here, we found that space between lamellae may have contributed to this difference between wild ecotypes and that there are no differences in overall gill surface area among lab-reared ecotypes when it is taken into account.
We found differences in a number of morphological and physiological traits predicted to respond to the more active and energy demanding lifestyle of the dwarf fish, associated with occupying the limnetic niche. In combination with Dalziel et al. , a total of 23 traits related to body shape, brain size, liver size, gill surface area, ventricle mass, ventricle and skeletal muscle metabolism, mitochondrial function and hematocrit show a predominant genetic basis for the variation observed in these traits between dwarf and normal Lake Whitefish. The two studies together also revealed that eleven traits related to ventricle mass, mitochondrial function and liver mass show an induced-environmental response for at least one ecotype. Overall, these data suggest that local adaptation (more than acclimation) is the major mechanism underlying the divergence of swimming activity between dwarf and normal whitefish, suggesting that there is a low probability for fish to switch trophic niches within a lifetime in the wild. The well-documented reproductive barriers between these ecotypes and genetically based differences in several behavioural, physiological and morphological traits underlying differences in swimming activity also support this hypothesis [48–55].
We would like to thank J-C. Therrien, M-C. Dion, S. Higgins, and I. Langlois-Parisé from the Laboratoire Regional des Sciences Aquatiques for assistance rearing whitefish, G. Côté, R. Janvier, L. Papillion, and A-M. Dion-Côté for assistance in the laboratory, and A. R. Kusler for valuable suggestions to improve earlier versions of this manuscript. We are also grateful to M. Barluenga and two anonymous reviewers for their constructive inputs on a previous version.
This work was supported by a Natural Science and Engineering Research Council of Canada (NSERC) discovery grant to LB, the Canadian Research Chair in Genomics and Conservation of Aquatic Resources to LB, an NSERC postdoctoral fellowship to ACD and ML, in addition to a Fonds de recherche du Québec – Nature et technologies postdoctoral fellowship to ML. This project is a contribution to the research program of Ressources Aquatique Québec (RAQ).
Availability of data and materials
Supporting data can be made available upon request.
ML wrote the initial manuscript and analyse the data. All authors have contributed to editing and revising the manuscript. ML, ACD and LB designed the project. ML, ACD and NM carried out the experiment and fish sampling. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publications
Ethics approval and consent to participate
All protocols were approved by Université Laval’s animal care committee prior to the sampling (Protocol 82178).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Kawacki TJ, Ebert D. Conceptual issues in local adaptation. Ecol Lett. 2004;7:1225–41.View ArticleGoogle Scholar
- Savolainen O, Lascoux M, Merilä J. Ecological genomics of local adaptation. Nat Rev Genet. 2013;14:807–20.View ArticlePubMedGoogle Scholar
- Willette DA, Allendorf FW, Barber PH, Barshis DJ, Carpenter KE, Crandall ED, Cresko WA, Fernandez-Silva I, Matz MV, Meyer E, Santos MD, Seeb LW, Seeb JE. So, you want to use next-generation sequencing in marine systems? Insight from the Pan-Pacific Advanced Studies Institute. Bull Mar Sci. 2014;90:79–122.View ArticleGoogle Scholar
- Merilä J, Hendry AP. Climate change, adaptation, and phenotypic plasticity: the problem and the evidence. Evol Appl. 2014;7:1–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Piersma T, Drent J. Phenotypic flexibility and the evolution of organismal design. Trends Ecol Evol. 2003;18:228–33.View ArticleGoogle Scholar
- Price TD, Qvarnstrom A, Irwin DE. The role of phenotypic plasticity in driving genetic evolution. Proc Roy Soc B. 2003;270:1433–40.View ArticleGoogle Scholar
- Thibert-Plante X, Hendry AP. The consequences of phenotypic plasticity for ecological speciation. J Evol Biol. 2011;24:326–42.View ArticlePubMedGoogle Scholar
- Robinson BW, Wilson DS. Genetic variation and phenotypic plasticity in a polymorphic population of pumpkinseed sunfish (Lepomis gibbosus). Evol Ecol. 1996;10:631–52.View ArticleGoogle Scholar
- Taylor EB. Species pairs of north temperate freshwater fishes: evolution, taxonomy, and conservation. Rev Fish Biol Fish. 1999;9:299–334.View ArticleGoogle Scholar
- Robinson BW, Parson KJ. Changing times, spaces, and faces: tests and implications of adaptive morphological plasticity in the fish of northern postglacial lakes. Can J Fish Aquat Sci. 2002;59:1819–33.View ArticleGoogle Scholar
- Hendry AP. Ecological speciation! Or the lack thereof ? Can J Fish Aquat Sci. 2009;66:1383–98.View ArticleGoogle Scholar
- Seehausen O, Wagner CE. Speciation in freshwater fishes. Annu Rev Ecol Evol Syst. 2014;45:621–51.View ArticleGoogle Scholar
- Dynes J, Magnan P, Bernatchez L, Rodriguez MA. Genetic and morphological variations between two forms of lacustrine brook charr. J Fish Biol. 1999;54:955–72.View ArticleGoogle Scholar
- Østbye K, Amundsen PA, Bernatchez L, Klemetsen A, Knudsen R, Kristoffersen R, Naesje TF, Hindar K. Parallel evolution of ecomorphological traits in the European Whitefish Coregonus lavaretus (L.) species complex during postglacial times. Mol Ecol. 2006;15:3983–4001.View ArticlePubMedGoogle Scholar
- Laporte M, Magnan P, Angers B. Genetic differentiation between the blue and the yellow phenotypes of walleye (Sander vitreus): An example of parallel evolution. Ecoscience. 2011;18:124–9.View ArticleGoogle Scholar
- Dalziel AC, Vines TM, Schulte PM. Repeated reductions in prolonged swimming capacity following freshwater colonization in threespine sticklebacks. Evolution. 2012;66:1226–39.View ArticlePubMedGoogle Scholar
- Merilä J. Lakes and ponds as model systems to study parallel evolution. J Limnol. 2014;73:33–45.Google Scholar
- Bernatchez L, Wilson CC. Comparative phylogeography of Nearctic and Palearctic fishes. Mol Ecol. 1998;7:431–52.View ArticleGoogle Scholar
- Taberlet P, Fumagalli L, Wustsaucy AG, Cosson JF. Comparative phylogeography and postglacial colonisation routes in Europe. Mol Ecol. 1998;7:453–64.View ArticlePubMedGoogle Scholar
- Skùlason S, Smith TB. Resource polymorphisms in vertebrates. Trends Ecol Evol. 1995;10:366–70.View ArticlePubMedGoogle Scholar
- Smith TB, Skùlason S. Evolutionary significance of resource polymorphisms in fish, amphibians and birds. Annu Rev Ecol Evol Syst. 1996;27:111–33.View ArticleGoogle Scholar
- Schluter D, McPhail JD. Ecological character displacement and speciation in sticklebacks. Am Nat. 1992;140:85–108.View ArticlePubMedGoogle Scholar
- Bernatchez L, Renaut S, Whiteley AR, Campbell D, Derôme N, Jeukens J, Landry L, Lu G, Nolte AW, Østbye K, Rogers SM, St-Cyr J. On the origin of species: Insights from the ecological genomics of whitefish. Phil Trans Roy Soc Lond. 2010;365:1783–800.View ArticleGoogle Scholar
- Willacker JJ, Von Hippel FA, Wilton PR, Walton KM. Classification of threespine stickleback along the benthic–limnetic axis. Biol J Lin Soc. 2010;101:595–608.View ArticleGoogle Scholar
- Webb PW. Locomotor patterns in the evolution of actinopterygian fishes. Am Zool. 1982;22:329–42.View ArticleGoogle Scholar
- Webb PW. Form and function in fish swimming. Sci Am. 1984;251:58–68.View ArticleGoogle Scholar
- Langerhans RB, Reznick DN. Ecology and evolution of swimming performance in fishes: predicting evolution with biomechanics. In: Domenici P, Kapoor BG, editors. Fish locomotion: an etho-ecological perspective. Enfield: Science Publishers; 2010. p. 200–48.View ArticleGoogle Scholar
- Del Giorgio PA, Gasol JM. Biomass distribution in freshwater plankton communities. Am Nat. 1995;146:135–52.View ArticleGoogle Scholar
- Roff DA. The evolution of life history parameters in teleosts. Can J Fish Aquat Sci. 1984;41:989–1000.View ArticleGoogle Scholar
- Trudel M, Tremblay A, Schetagne R, Rasmussen JB. Why are dwarf fish so small? An energetic analysis of polymorphism in lake whitefish (Coregonus clupeaformis). Can J Fish Aquat Sci. 2001;58:394–405.View ArticleGoogle Scholar
- Rosenfeld J, Van Leeuwen T, Richards J, Allen D. Relationship between growth and standard metabolic rate: measurement artefacts and implications for habitat use and life-history adaptation in salmonids. J Anim Ecol. 2015;84:4–20.View ArticlePubMedGoogle Scholar
- Walker JA. An integrative model of evolutionary covariance: a symposium on body shape in fishes. Integr Comp Biol. 2010;50:1051–6.View ArticlePubMedGoogle Scholar
- Konarzewski M, Książek A. Determinants of intra-specific variation in basal metabolic rate. J Comp Physiol B. 2013;183:27–41.View ArticlePubMedGoogle Scholar
- Proulx R, Magnan P. Contribution of phenotypic plasticity and heredity to the trophic polymorphism of lacustrine brook charr (Salvelinus fontinalis M.). Evol Ecol Res. 2004;6:503–22.Google Scholar
- Harrod C, Mallela J, Kahilainen KK. Phenotype-environment correlations in a putative whitefish adaptive radiation. J Anim Ecol. 2010;79:1057–68.View ArticlePubMedGoogle Scholar
- Walker JA. Ecological morphology of lacustrine threespine stickleback, Gasterosteus aculeatus L. (Gasterosteidae) body shape. Biol J Linn Soc. 1997;61:3–50.Google Scholar
- Aiello LC, Wheeler P. The expensive-tissue hypothesis: the brain and the digestive system in human and primate evolution. Curr Anthropol. 1995;36:199–221.View ArticleGoogle Scholar
- Gallagher D, Belmonte D, Deurenberg P, Wang Z, Krasnow N, Pi-Sunyer FX, Heymsfield SB. Organ-tissue mass measurement allows modeling of REE and metabolically active tissue mass. Am J Physiol Endocrinol Metab. 1998;275:249–58.Google Scholar
- Müller MJ, Langemann D, Gehrke I, Later W, Heller M, Glüer CC, Heymsfield SB, Bosy-Westphal A. Effect of constitution on mass of individual organs and their association with metabolic rate in humans—a detailed view on allometric scaling. PLoS One. 2011;6:e22732.View ArticlePubMedPubMed CentralGoogle Scholar
- Boldsen MM, Norin T, Malte H. Temporal repeatability of metabolic rate and the effect of organ mass and enzyme activity on metabolism in European eel (Anguilla anguilla). Comp Biochem Physiol. 2013;165:22–9.View ArticleGoogle Scholar
- Gonzalez RJ, McDonald DG. The relationship between oxygen consumption and ion loss in a freshwater fish. J Exp Biol. 1992;163:317–32.Google Scholar
- Duthie GG, Hughes GM. The effects of reduced gill area and hyperoxia on the oxygen consumption and swimmin speed of rainbow trout. J Exp Biol. 1986;127:349–54.Google Scholar
- Bernatchez L, Dodson JJ. Allopatric origin of sympatric populations of lake whitefish (Coregonus clupeaformis) as revealed by mitochondrial DNA restriction analysis. Evolution. 1990;44:1263–71.View ArticleGoogle Scholar
- Lu G, Bernatchez L. Correlated trophic specialization and genetic divergence in sympatric lake whitefish ecotypes (Coregonus clupeaformis): support for the ecological speciation hypothesis. Evolution. 1999;53:1491–505.View ArticleGoogle Scholar
- Lu G, Basley DJ, Bernatchez L. Contrasting patterns of mitochondrial DNA and microsatellite introgressive hybridization between lineages of lake whitefish (Coregonus clupeaformis): relevance for speciation. Mol Ecol. 2001;10:965–85.View ArticlePubMedGoogle Scholar
- Landry L, Vincent WF, Bernatchez L. Parallelism between limnological features and phenotypic evolution of lake whitefish dwarf ecotypes. J Evol Biol. 2007;20:971–84.View ArticlePubMedGoogle Scholar
- Landry L, Bernatchez L. Role of epibenthic resource opportunities in the parallel evolution of lake whitefish species pairs (Coregonus sp.). J Evol Biol. 2010;23:2602–13.View ArticlePubMedGoogle Scholar
- Lu G, Bernatchez L. Experimental evidence for reduced hybrid viability between dwarf and normal ecotypes of lake whitefisg (Coregonus clupeaformis Mitchill). Proc Roy Soc B. 1998;265:1025–30.View ArticleGoogle Scholar
- Rogers SM, Bernatchez L. The genetic basis of intrinsic and extrinsic postzygotic reproductive isolation jointly promoting speciation in the lake whitefish species complex (Coregonus Clupeaformis). J Evol Biol. 2006;19:1979–94.View ArticlePubMedGoogle Scholar
- Dion-Côté A-M, Renaut S, Normandeau E, Bernatchez L. RNA-seq reveals transcriptomic shock involving transposable elements reactivation in hybrids of young lake whitefish species. Mol Biol Evol. 2014;31:1188–99.View ArticlePubMedGoogle Scholar
- Dion-Côté A-M, Symonová R, Ráb P, Bernatchez L. Reproductive isolation in a nascent species pair is associated with aneuploidy in hybrid offspring. Proc Roy Soc B. 2015;282:20142862.View ArticleGoogle Scholar
- Whiteley AR, Persaud KN, Derôme N, Montgomerie R, Bernatchez L. Reduced sperm performance in backcross hybrids between whitefish species-pairs (Coregonus sp.). Can J Zool. 2009;87:566–72.View ArticleGoogle Scholar
- Rogers MR, Gagnon V, Bernatchez L. Genetically based phenotype-environment association for swimming behavior in Lake Whitefish ecotypes (Coregonus clupeaformis Mitchill). Evolution. 2002;56:2322–9.View ArticlePubMedGoogle Scholar
- Rogers SM, Bernatchez L. The Genetic architecture of ecological speciation and the association with signatures of selection in natural lake whitefish (Coregonus sp. Salmonidae) species pairs. Mol Biol Evol. 2007;24:1423–38.View ArticlePubMedGoogle Scholar
- Laporte M, Rogers SM, Dion-Côté A-M, Normandeau E, Gagnaire P-A, Dalziel AC, Chebib J, Bernatchez L. RAD‐QTL mapping reveals both genome‐level parallelism and different genetic architecture underlying the evolution of body shape in Lake Whitefish (Coregonus clupeaformis) species pairs. G3 (Bethesda). 2015;5:1481–91.View ArticleGoogle Scholar
- Evans ML, Bernatchez L. Oxidative phosphorylation gene transcription in whitefish species pairs reveals patterns of parallel and non-parallel physiological divergence. J Evol Biol. 2012;25:1823–34.View ArticlePubMedGoogle Scholar
- Evans ML, Præbel K, Peruzzi S, Bernatchez L. Parallelism in the oxygen transport system of the lake whitefish: the role of physiological adaptive divergence in ecological speciation. Mol Ecol. 2012;21:4038–50.View ArticlePubMedGoogle Scholar
- Dalziel AC, Martin N, Laporte M, Guderley H, Bernatchez L. Acclimation and adaptation of aerobic exercise physiology in Lake Whitefish ecotype (Coregonus clupeaformis). Evolution. 2015;69:2167–86.View ArticlePubMedGoogle Scholar
- Evans ML, Chapman LJ, Mitrofanov I, Bernatchez L. Variable extent of parallelism in respiratory, circulatory, and neurological traits across lake whitefish species pairs. Ecol Evol. 2013;3:546–57.View ArticlePubMedPubMed CentralGoogle Scholar
- Derôme N, Bernatchez L. The transcriptomics of ecological convergence between two limnetic coregonine fishes (Salmonidae). Mol Biol Evol. 2006;23:2370–8.View ArticlePubMedGoogle Scholar
- St-Cyr J, Derôme N, Bernatchez L. The transcriptomics of life-history trade-offs in whitefish species pairs (Coregonus sp.). Mol Ecol. 2008;17:1850–70.View ArticlePubMedGoogle Scholar
- Jeukens J, Bittner D, Knudsen R, Bernatchez L. Candidate genes and adaptive radiation: Insights from transcriptional adaptation to the limnetic niche among coregonine fishes (Coregonus sp., Salmonidae). Mol Biol Evol. 2009;26:155–66.View ArticlePubMedGoogle Scholar
- Nolte AW, Renaut S, Bernatchez L. Divergence in gene regulation at young life history stages of whitefish (Coregonus sp.) and the emergence of genomic isolation. BMC Evol Biol. 2009;9:59.View ArticlePubMedPubMed CentralGoogle Scholar
- Renaut S, Nolte AW, Bernatchez L. Mining transcriptome sequences towards identifying adaptive Single Nucleotide Polymorphisms in lake whitefish species pairs (Coregonus sp.). Mol Ecol. 2010;19:115–31.View ArticlePubMedGoogle Scholar
- Laporte M, Claude J, Berrebi P, Perret P, Magnan P. Shape plasticity in response to water velocity in the freshwater blenny Salaria fluviatilis. J Fish Biol. 2016;88:1191–203.View ArticlePubMedGoogle Scholar
- Gonda A, Välimäki K, Herczeg G, Merilä J. Brain development and predation: plastic responses depend on evolutionary history. Biol Lett. 2012;8:249–52.View ArticlePubMedGoogle Scholar
- Crispo E, Chapman LJ. Geographic variation in phenotypic plasticity in response to dissolved oxygen in an African cichlid fish. J Evol Biol. 2010;23:2091–103.View ArticlePubMedGoogle Scholar
- Jacobsen MW, Hansen MM, Orlando L, Bekkevold D, Bernatchez L, Willerslev E, Gilbert MT. Mitogenome sequencing reveals shallow evolutionary histories and recent divergence time between morphologically and ecologically distinct European whitefish (Coregonus spp.). Mol Ecol. 2012;21:2727–42.View ArticlePubMedGoogle Scholar
- Bernatchez L, Dodson JJ. Influence of temperature and current speed on the swmiming capacity of lake whitefish (Coregonus clupeaformis) and cisco (C. artedii). Can J Fish Aquat Sci. 1985;42:1522–9.View ArticleGoogle Scholar
- Anttila K, Jaevilehto M, Manttari S. The swimming performance of brown trout and whitefish: the effects of exercise on Ca2+ handling and oxidative capacity of swimming muscles. J Comp Physiol B. 2008;178:465–75.View ArticlePubMedGoogle Scholar
- Claude J. Morphometrics with R. New York: Springer Science + Business Media, L.L.C; 2008.Google Scholar
- Rohlf FJ. tpsDig, digitize landmarks and outlines, version 2.16. Department of Ecology and Evolution, State University of New York at Stony Brook. 2010.Google Scholar
- Zelditch ML, Swiderki DL, Sheets HD, Fink WL. Geometric morphometrics for biologists. San Diego: Elsevier; 2004.Google Scholar
- Albert AY, Sawaya S, Vines TH, Knetch AK, Miller CT, Summers BR, Balabhadra S, Kingsley DM, Schulter D. The genetics of adaptive shape shift in stickleback: pleiotropy and effect size. Evolution. 2008;62:76–85.PubMedGoogle Scholar
- Yezerinac SM, Lougheed SC, Handford P. Measurement error and morphometric studies: Statistical power and observer experience. Syst Biol. 1992;41:471–82.View ArticleGoogle Scholar
- Claude J, Paradis E, Tong H, Auffray J-C. A geometric morphometric assessment of the effects of environment and cladogenesis on the evolution of the turtle shell. Biol J Linn Soc. 2003;79:485–501.View ArticleGoogle Scholar
- Rohlf FJ, Slice D. Extensions of the Procrustes method for the optimal superimposition of landmarks. Syst Zool. 1990;39:40–59.View ArticleGoogle Scholar
- Dryden IL, Mardia KV. Statistical shape analysis. New York: Wiley; 1998.Google Scholar
- Klingenberg CP. MorphoJ: an integrated software package for geometric morphometrics. Mol Ecol Res. 2011;11:353–7.View ArticleGoogle Scholar
- Legendre P, Legendre L. Numerical Ecology, Second English Edition. Amsterdam: Elsevier Science BV; 1998.Google Scholar
- Rasband WS. ImageJ. Bethesda, MD: US National Institutes of Health, http://imagej.nih.gov/ij/. 2011.
- Larsen S, Nielsen J, Hansen CN, Nielsen LB, Wibrand F, Stride N, Schroder HD, Boushel R, Helge JW, Dela F, Hey-Mogensen M. Biomarkers of mitochondrial content in skeletal muscle of healthy young human subjects. J Physiol. 2012;590:3349–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Farrell AP, Johansen JA, Suarez RK. Effects of exercise training on cardiac performance and muscle enzymes in rainbow trout, Oncorhynchus mykiss. Fish Physiol Biochem. 1991;9:303–12.View ArticlePubMedGoogle Scholar
- Martin CI, Johnston IA. The role of myostatin and the calcineurin-signalling pathway in regulating muscle mass in response to exercise training in the rainbow trout Oncorhynchus mykiss. J Exp Biol. 2005;208:2083–90.View ArticlePubMedGoogle Scholar
- Palstra AP, Mes D, Kusters K, Roques JAC, Flik G, Kloet K, Blonk RJW. Forced sustained swimming exercise at optimal speed enhances growth of juvenile yellowtail kingfish (Seriola lalandi). Front Physiol. 2015;5:506.View ArticlePubMedPubMed CentralGoogle Scholar
- Fish FE. Swimming strategies for energy economy. In: Domenici P, Kapoor BG, editors. Fish locomotion: an etho-ecological perspective. Enfield: Science Publishers; 2010. p. 90–122.View ArticleGoogle Scholar
- Schluter D. The Ecology of adaptive Radiation. Oxford: Oxford University Press; 2000.Google Scholar
- Nosil P. Ecological speciation. Oxford: Oxford University Press; 2012.View ArticleGoogle Scholar
- Angers B, Castonguay E, Massicotte R. Environmentally induced phenotypes and DNA methylation: how to deal with unpredictable conditions until the next generation and after. Mol Ecol. 2010;19:1283–95.View ArticlePubMedGoogle Scholar
- Domenici P, Kapoor BG. Fish locomotion: an etho-ecological perspective. Enfield: Science Publishers, UK; 2010.View ArticleGoogle Scholar
- Rice AN, Hale ME. Role of locomotions in feeding. In: Domenici P, Kapoor BG, editors. Fish locomotion: an etho-ecological perspective. Enfield: Science Publishers; 2010. p. 171–99.View ArticleGoogle Scholar
- Lundsgaard-Hansen B, Matthews B, Vonlanthen P, Taverna A, Seehausen O. Adaptive plasticity and genetic divergence in feeding efficiency during parallel adaptive radiation of whitefish (Coregonus spp.). J Evol Biol. 2013;26:483–98.View ArticlePubMedGoogle Scholar
- Bauchot R, Randall JE, Ridet JM, Bauchot ML. Encephalization in tropical teleost fishes and comparison with their mode of life. J für Hirnf. 1989;30:645–69.Google Scholar
- Kotrschal K, Van Staaden MJ, Huber R. Fish brains: evolution and environmental relationships. Rev Fish Biol Fish. 1998;8:373–408.View ArticleGoogle Scholar
- Chapman LJ, Hulen KG. Implications of hypoxia for the brain size and gill morphometry of mormyrid fishes. J Zool. 2001;254:461–72.View ArticleGoogle Scholar
- Brooks JL, Dodson SI. Predation, body size, and composition of plankton. Science. 1965;150:28–35.View ArticlePubMedGoogle Scholar
- Rolfe DFS, Brown GC. Cellular energy utilization and molecular origin of standard metabolic rate in mammals. Physiol Rev. 1997;77:731–58.PubMedGoogle Scholar
- Hulbert AJ, Else PL. Membranes and the setting of energy demand. J Exp Biol. 2005;208:1593–9.View ArticlePubMedGoogle Scholar
- Walker JA. A general model of functional constraints on phenotypic evolution. Am Nat. 2007;170:681–9.View ArticlePubMedGoogle Scholar