Contemporary divergence in early life history in grayling (Thymallus thymallus)
© Thomassen et al; licensee BioMed Central Ltd. 2011
Received: 19 August 2011
Accepted: 13 December 2011
Published: 13 December 2011
Skip to main content
© Thomassen et al; licensee BioMed Central Ltd. 2011
Received: 19 August 2011
Accepted: 13 December 2011
Published: 13 December 2011
Following colonization of new habitats and subsequent selection, adaptation to environmental conditions might be expected to be rapid. In a mountain lake in Norway, Lesjaskogsvatnet, more than 20 distinct spawning demes of grayling have been established since the lake was colonized, some 20-25 generations ago. The demes spawn in tributaries consistently exhibiting either colder or warmer temperature conditions during spawning in spring and subsequent early development during early summer. In order to explore the degree of temperature-related divergence in early development, a multi-temperature common-garden experiment was performed on embryos from four different demes experiencing different spring temperatures.
Early developmental characters were measured to test if individuals from the four demes respond differently to the treatment temperatures. There was clear evidence of among-deme differences (genotype - environment interactions) in larval growth and yolk-to-body-size conversion efficiency. Under the cold treatment regime, larval growth rates were highest for individuals belonging to cold streams. Individuals from warm streams had the highest yolk-consumption rate under cold conditions. As a consequence, yolk-to-body-mass conversion efficiency was highest for cold-deme individuals under cold conditions. As we observed response parallelism between individuals from demes belonging to similar thermal groups for these traits, some of the differentiation seems likely to result from local adaptation
The observed differences in length at age during early larval development most likely have a genetic component, even though both directional and random processes are likely to have influenced evolutionary change in the demes under study.
Large changes in environmental conditions put populations at risk of becoming maladapted, thereby leading to negative population growth and potential extinction . Such large changes in environmental conditions may be due to anthropogenic activity , but may also be experienced as individuals invade a novel environment. For ectotherms in general, temperature is considered to be by far the most important external factor, controlling much of the variation seen in embryonic ontogenetic rates [3–5]. Following an environmental perturbation a population may be maladapted to the new conditions. If populations cannot rapidly move towards the new evolutionary optimum they may risk extinction .
Modification of traits may happen through purely plastic processes or through evolutionary change [7–9]. Phenotypic plasticity is defined as the ability of a single genotype to develop different phenotypes under different environmental conditions . Norms of reaction (i.e., the variation in trait value across environmental gradients) are commonly used to illustrate the phenotypic plasticity of genotypes (reviewed in ). When a perturbation to the environment is outside of the normal range of environmental fluctuations, it is expected that evolutionary change will be required to reach the new fitness optimum . In many situations, for example where there is substantial spatial habitat heterogeneity, the norms of reaction may be the target of selection . Divergence in the slope and/or elevation of reaction norms may develop among populations occupying habitats that differ in their environmental characteristics. Since many traits depend on conditions experienced earlier in life, the reaction norm concept has been extended to include developmental processes . These developmental reaction norms may be the target of selection, where individuals that most optimally modify their developmental trajectories under varying environmental conditions are selected.
A change in the population reaction norm requires genetic change, either through stochastic processes like genetic drift, or through adaptive evolution. Today it is realized that ecological and evolutionary dynamics do not necessarily operate on vastly different timescales [13–15]. A mounting number of studies provide evidence of evolutionary change over contemporary time scales [16, 17].
Grayling, a spring spawning salmonid fish, colonized the lake Lesjaskogsvatnet through a man-made connection in the upper reaches of the river Gudbrandsdalslågen that was opened in the late 1880s (20-25 generations ago). Subsequent closing of the connection made further immigration of grayling difficult . Genetic evidence suggests that Lesjaskogsvatnet was colonized by a low number of individuals from Gudbrandsdalslågen . Since the colonization event, the grayling have established more than 20 demes spawning in separate tributaries to the lake . It is likely that grayling return to their natal tributary to spawn. Lacustrine populations of grayling prefer to move to running water to spawn and are known to return to their natal stream with high fidelity .
During spring, the south-facing slopes of Lesjaskogsvatnet receive considerably more sunlight and are less steep than the shaded north-facing slope resulting in more rapid spring warming. In addition, small tributaries tend to warm more rapidly than large ones. As a result of these differences in aspect coupled with size differences among tributaries, tributaries can generally be classified as either large-and-cold (LC) or small-and-warm (SW) . Generally, spring ice-melt commences earlier in SW tributaries than in LC tributaries, and SW tributaries have a higher temperature than LC tributaries through spring and early summer . Typically, mean daily June- July water temperatures differ by 1-1.5 °C among LC and SW streams, which add up to substantial differences in accumulated temperature sums over the stream phase.
Haugen and Vøllestad  conducted a common-garden study of several early growth and developmental traits of grayling in the same area, but on an inter-lake scale. They investigated differences between the Lesjaskogsvatnet population and two other populations recently derived from Lesjaskogsvatnet grayling. In their study, they discovered apparent adaptations in many of the traits. Based on the same data, evidence for divergence in time to hatching was found between two of the populations . These results demonstrate that in the absence of gene flow local adaptation may progress significantly over less than a 20-generation time span. Here, we use a common environment approach at three different temperature conditions to investigate how temperature reaction norms for early life history traits vary for grayling from two cold and two warm streams where gene flow among streams is possible.
In high latitude mountain lakes like Lesjaskogsvatnet the growth season is short. For north-temperate fishes the first winter is usually a period of energy deficit , potentially selecting for faster growth. The selection intensity is expected to be strongest for the demes with the shortest growth season (i.e. LC demes). Individuals from these demes are expected to compensate for the shorter growth season with a higher general capacity for growth .
As grayling occupy divergent habitats during development it is possible that evolutionary divergence could result. However, one should bear in mind that this divergence is observed under conditions of continuing gene flow and potential meta-population dynamics [34, 35]. The grayling could have responded to the differing environmental conditions by evolutionary changes; alternatively, the grayling could have a plastic response common to individuals of all demes . If the latter is the case, reaction norms for all demes should appear more or less identical. If reaction norms differ this would be indicative of genetic effects.
In this study, we test if there is evidence for evolution in early life-history traits in larval grayling from four different spawning demes that are experiencing different environmental selection pressures. To contrast the development of larval grayling from different demes a common-environment experiment was conducted. In the experiment, individuals from four tributaries (two LC, two SW) were subjected to three experimental temperatures and studied from fertilization to swim up, the stage when the juveniles emerge from the gravel and become free-living .
Mean initial egg size (measured as egg diameter, mm) differed significantly among demes (one-way ANOVA: F3, 116 = 7.13, P = 0.0002, R2= 0.123). The largest and the smallest eggs were found in the warm stream demes (mean ± 2se; Hyrjon: 3.88 ± 0.06; Valåe: 3.94 ± 0.07; Steinbekken: 4.03 ± 0.09; Sandbekken: 3.79 ± 0.09).
Effect test for variation in larval length over time.
Sum of squares
Replicate [Deme, Treatment]
Effect test for variation in larval yolk sac size over time.
Sum of sq
Replicate [Deme, Treatment]
There was little support for thermal group parallelism in the yolk consumption pattern as fitting the fully factorial model with the deme effect substituted by group increased AIC by 74 units (full model with deme effect; AIC = 670, full model with temperature group effect; AIC = 744).
Effect test for variation in conversion efficiency.
Sum of sq
Replicate [Deme, Treatment]
The results presented in this study demonstrate genetic differences in growth-related traits between grayling from four demes within Lesjaskogsvatnet. Grayling from different demes responded differently to changes in developmental temperature, with respect to growth in length and yolk-to-body-size transition efficiency. This appeared as differentiated temperature-related reaction norms for these traits. A major part of the among-deme variation in growth-related traits is likely to be genetically based since the results were obtained from a common-environment experiment. Parallelism of developmental phenotypes between demes belonging to similar thermal groups indicates that the differentiation is consistent with local adaptation.
We found significant among-deme variation in growth, estimated as the change in length over time, and in how yolk was allocated to growth. As expected from previous work ( and references therein), these traits were also significantly influenced by temperature. However, the different demes responded differently to the treatment temperatures. Haugen and Vøllestad  also discovered similar differences when testing three grayling populations having the same ancestry, but experiencing vastly different environmental conditions.
The differences in the growth rate estimates for LC deme individuals relative to the SW individuals may, to some extent, be interpreted as indicative of local adaptation. Larvae from the LC demes grow faster than SW deme individuals under experimental conditions comparable to what is observed in the LC deme habitats (i.e., the cold treatment temperature, Figure 2). LC deme individuals also have more efficient yolk-to-body size transition efficiencies under cold and medium temperature conditions (Figure 3). All grayling demes originate from fish spawning downstream in the river Gudbrandsdalslågen. Overall, temperatures in Gudbrandsdalslågen is comparable to the temperatures in the SW streams, indicating that only the LC deme individuals have been exposed to temperatures leading to strong selection. The full picture emerging from the data, however, does not fully correspond with predictions for the "bigger-is-better"-hypothesis or counter-gradient variation hypothesis [39, 40]: at the medium and warm treatment temperatures, the rank order of demes appears not to follow a pattern that can be ascribed to conditions in their natal tributaries. Overall, there is strong evidence for genetically based differences in growth rate among demes. However, since the optimal strategy is unknown it is premature to decide whether the result is due to local adaptation or not. Further, there is a general trend in all demes for growth rates to increase with temperature, much in line with results from comparable studies .
Jensen et al.  found significant divergence in traits such as larvae length and growth rate of brown trout juveniles in a common-garden experiment, and suggested that this was due to local adaptation. However, as for the present study, fitness was not estimated so conclusions in this aspect are premature. In our study, differences in traits related to growth and energy utilization do not unambiguously correspond with predictions for local adaptation. Any lack of adaptation might be caused by evolutionary constraints, resulting from the history of the Lesjaskogsvatnet grayling. It is possible that the demes under investigation have not had equal opportunities for environmental adaptation and that equilibria have not been reached [34, 35]. When investigating potential genetic effects in divergence and adaptation studies one should never overlook the potential importance of random genetic drift. Through the colonization process, the grayling of Lesjaskogsvatnet have gone through various severe bottlenecks [34, 35]. Normally, bottlenecks are expected to reduce genetic variability and hence impede response to selection (e.g. [42, 43]). However, because founder populations may contain a high fraction of the heterozygosity and additive genetic variance in quantitative characters that existed in a large source population, prospects for adaptation may not necessarily have to be impeded . Potentially, genetic drift could complicate the interpretation of results since the outcome of evolution by genetic drift cannot be predicted (e.g. ). Further, in young systems such as in the Lesjaskogsvatnet grayling system gene flow and drift may vary in strength, leading to large among year variation in strength of selection . Gene flow most often is considered to hinder local adaptation [46, 47], however, this is not necessarily always the case [48, 49]. Gene flow among inbred demes may counteract the effect of drift and increase the evolvability of these demes by introducing genetic variation, and as such enhance short-term adaptation . However, in most situations where you have heterogeneous habitats and genotype-environment interaction effects, gene flow is predicted to constrain local adaptation [49, 50].
The different traits studied here may be influenced by maternal effects. If such effects are present, and if they are physiological side effects and not consequences of adaptive evolution, they have the potential to confound our results. Potential maternal effects might be egg provisioning of energy and other important nutrients . This cannot be directly accounted for in a study like ours, where the maternal environment is not controlled. One complicating factor is the time delay between maturation of gonads for LC and SW demes, potentially leading to quality differences in the eggs. In Haugen and Vøllestad's  study on grayling from the same area, egg size was found to have a significant contribution to the total female effect for three out of four early-life size traits and also on early survival rate. Our study design did not allow for family-level resolution. Egg size - a main determinant of maternal effects - did differ among the four demes. However, these differences cannot explain the parallelism in growth rates and energy use as both the largest and the smallest eggs were observed in the SW demes.
We have shown that several early life history traits have diverged among the four grayling demes studied here, and some of the divergence is consistent with adaptive predictions. Taken together with several related studies from the same area [29, 30] this indicates that these traits can evolve rapidly and in spite of conditions that are generally assumed to constrain adaptation. We have shown that most of the demes studied have low genetic variability and have experienced bottlenecks and founder effects [34, 35]. Gene flow among demes is strong and highly variable . Additionally, some demes are small indicating that genetic drift may be strong. These conditions can potentially lead to strong genetic drift, a process that can constrain adaptive evolution [43, 46, 47]. A recent theoretical study does, however, show that dispersal, mutation and sexual reproduction may accelerate local adaptation in growing populations , especially when populations are small. However, these same processes may reduce local adaptation in the long run.
In summary, we find evidence for significant among-deme variation in important early life-history traits in grayling from the Lesjaskogsvatnet. These differences have evolved during 20-25 generations in the face of gene flow, recent bottlenecks, and the fact that the grayling live in sympatry during most of their life. Some of these results are consistent with adaptation expectations as inferred from signs of parallelism, but we cannot rule out that other evolutionary or non-evolutionary processes have influenced the results.
Information on parental fish.
226.8 ± 130.3
287.6 ± 56.4
337.6 ± 89.3
333.25 ± 28.2
187.2 ± 85.6
268.6 ± 11.5
263.9 ± 100.7
305.6 ± 35.2
158.8 ± 90.4
254.8 ± 30.6
244.0 ± 61.2
299.3 ± 22.1
236.4 ± 106.8
284.4 ± 51.7
298.7 ± 82.2
327.0 ± 41.9
The eggs were pooled by deme using an equal volume (100 ml) of eggs from each female. The pooled batch of eggs from each deme was then split into a number of batches equal to the number of males from that deme. To avoid sperm competition and so maximise the number of families, each batch of pooled eggs was subsequently fertilized with sperm from one male. Following fertilization, the batches of eggs from each deme were mixed together again before they were partitioned into three treatment groups. Each treatment group was split into two replicates. Three separate experimental tanks containing water of three different temperatures (5.83 ± 0.43°C, 8.14 ± 0.29 °C, 10.02 ± 0.28 °C) were used. These temperatures were chosen to represent lower, medium and upper temperatures experienced by developing grayling larvae in nature , and thus we would be able to investigate norms of reaction covering the range of temperatures that grayling may be expected to tolerate. Mean summer (June - July) temperatures in the four streams investigated here differ strongly, with the two small and warm streams being approximately 1-1.5 ºC higher than the large and cold streams (Sandbekken 8.44 ± 0.52 (n = 5); Steinbekken 8.81 ± 0.60 (n = 4); Hyrjon 7.40 ± 0.94 (n = 8); Valåe 7.28 ± 0.69 (n = 0.69)). This adds to a large temperature-sum difference among streams during egg and larvae development.
Inflowing water was comprised of activated-charcoal filtered tap water. Temperatures were registered using temperature loggers (HOBO). Fertilized eggs were placed in porous containers suspended in the large treatment tanks as described in detail previously . Each of the three tanks contained two replicate containers from each deme.
Animal sampling and experimentation were performed in compliance with the recommendations of National Animal Research Authority (permission ID 2008/7368.5) and under the supervision of authorized investigators.
44 - 12
59 - 22
89 - 57
34 - 16
51 - 20
82 - 48
39 - 19
60 - 21
86 - 55
35 - 16
52 - 18
81 - 51
42 - 18
61 - 32
85 - 16
36 - 15
52 - 24
62 - 29
50 - 17
62 - 43
89 - 63
40 - 16
60 - 27
80 - 56
175 - 66
242 - 118
349 - 191
145 - 63
215 - 89
305 - 184
All samples were visually inspected to check for hatched individuals. Numbers of hatched individuals per sample was recorded. However, precise estimates of time at hatching could not be obtained because of large among-sample variation. Photographs of eggs and larvae were obtained using a Leica DC300 digital camera mounted on a Leica MZ8 stereo microscope connected to a computer and captured using IrfanView (version 3.99) http://www.irfanview.com/. One hundred and twenty eggs (5 eggs from each deme in each treatment for the first two days after fertilization) were photographed and initial egg size was estimated as the mean egg diameter measured in UTHSCSA Image Tool (version 3.0; http://ddsdx.uthscsa.edu/dig/itdesc.html). Out of the 2142 individuals sampled (see Table 5), 717 had hatched and these where measured for length and yolk sac area. All larvae were photographed twice; usually once with the left side facing the camera, and once with the right side facing the camera to allow repeated estimates of traits to be made. The mean of the two measurements was used in the analysis. Based on a random subsample of 100 individuals there was a mean deviation of 1.06% between the two length measurements done on each individual. The corresponding number for yolk sac area estimates was 7.77%. On certain occasions larvae were fixed in awkward positions and were too brittle to straighten out. If so, photographs in three positions (i.e. side facing up, back facing up and/or yolk sac facing up) were obtained in order to get the measurements needed. Length was measured from the tip of the cranium to the visible end of the notochord. ImageJ version 1.38 http://rsbweb.nih.gov/ij/download.html was used to obtain a measurement of the yolk sac area. Yolk sac area was estimated by fitting an ellipse to a polygon drawn along the edge of the yolk sac in order to estimate the longest and shortest possible axis of the fitted ellipse. Axis lengths obtained from the two photographs of each individual were averaged. The resulting estimates were used to estimate each individual's yolk sac area using the formula for the area of an ellipse (πab; a and b being one-half of the ellipse's major and minor axes, respectively).
Six individuals were completely excluded from the statistical analysis due to obvious disfigurement or severe damage during handling making any measurements uncertain. In four individuals the yolk sac was damaged or completely removed due to handling, making measurements only on length possible. One individual had an intact yolk sac, but a damaged body, making only measurement of the yolk sac area possible.
A one-way ANOVA on egg diameter with deme as effect was carried out, and a Tukey-Kramer HSD (honestly significant difference) test was used to compare deme means.
Here, "Y" is the response variable, ßs correspond to fixed parameters under estimation and b(ij)k corresponds to the random effect of replicate (k) nested under treatment level i and deme j. Variables are provided in italics where "°d" indicates degree-days (ln-transformed), ε ijk represent the between individual variance (i.e., the residual variance) which is assumed to be normally distributed within replicates, and for a given link function. An asterisk indicates an interaction term. If non-linear trends were found in LME residual plots, new models including a squared term for the degree-days effect were fitted. In order to check for possible first-order temporal autocorrelation violations, Durbin-Watson tests were conducted for the fitted LME models.
In order to explore differentiation in the yolk-to-body-size conversion efficiency (see ), we fitted LME models where body length was predicted from yolk sac area and degree-days and with rearing temperature and deme as fixed effects. Replicates were nested under deme and temperature regime and were modelled as random effects. The full model structure was then compared with a comparable model where the deme structure was exchanged with a temperature-group (SW vs LC) structure, using Akaike Information Criterion (AIC)-values for comparison . Due to collinearity between the predictor variables yolk-sac area and degree-days, information about variance inflation (assessed from the general variance inflation factor, GVIF) and the corresponding tolerance is presented together with the chosen model .
The statistical analyses of the data were carried out using the statistical software JMP 5.0 (SAS Institute Inc., Cary, NC, USA) and R version 2.10.1 . LME models were fitted in R using the lmer-procedure in the LME4 library. Information about variance inflation due to collinearity (assessed from the general variance inflation factor, GVIF) was retrieved using the CAR library.
The study was financed by the Norwegian Research Council. We thank Erica Leder, Kai Rune Båtstad and Håvard Kauserud for important help during fieldwork.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.