Testing for local adaptation in brown trout using reciprocal transplants
© Stelkens et al.; licensee BioMed Central Ltd. 2012
Received: 27 April 2012
Accepted: 15 December 2012
Published: 18 December 2012
Skip to main content
© Stelkens et al.; licensee BioMed Central Ltd. 2012
Received: 27 April 2012
Accepted: 15 December 2012
Published: 18 December 2012
Local adaptation can drive the divergence of populations but identification of the traits under selection remains a major challenge in evolutionary biology. Reciprocal transplant experiments are ideal tests of local adaptation, yet rarely used for higher vertebrates because of the mobility and potential invasiveness of non-native organisms. Here, we reciprocally transplanted 2500 brown trout (Salmo trutta) embryos from five populations to investigate local adaptation in early life history traits. Embryos were bred in a full-factorial design and raised in natural riverbeds until emergence. Customized egg capsules were used to simulate the natural redd environment and allowed tracking the fate of every individual until retrieval. We predicted that 1) within sites, native populations would outperform non-natives, and 2) across sites, populations would show higher performance at ‘home’ compared to ‘away’ sites.
There was no evidence for local adaptation but we found large differences in survival and hatching rates between sites, indicative of considerable variation in habitat quality. Survival was generally high across all populations (55% ± 3%), but ranged from 4% to 89% between sites. Average hatching rate was 25% ± 3% across populations ranging from 0% to 62% between sites.
This study provides rare empirical data on variation in early life history traits in a population network of a salmonid, and large-scale breeding and transplantation experiments like ours provide powerful tests for local adaptation. Despite the recently reported genetic and morphological differences between the populations in our study area, local adaptation at the embryo level is small, non-existent, or confined to ecological conditions that our experiment could not capture.
Natural selection can lead to local adaptation, allowing individuals to have higher relative fitness in their native habitat compared to non-native individuals [1, 2]. As a consequence, selection against immigrants can reduce gene flow between populations [3, 4], which can lead to divergence without physical isolation [5–7]. But evidence for local adaptation is hard to come by. It usually requires a sound knowledge of the populations’ ecology, the prevailing selection agents, and ideally the underlying genetic basis of phenotypic adaptations. An alternative way, not requiring much preexisting knowledge, is reciprocal transplant experiments [8–11]. Reciprocal transplants can also help to disentangle local adaptation (in the sense of genetically based phenotypic divergence) from phenotypic plasticity (in the sense of differential reactions norms).
Local adaptation in salmonid fishes seems to occur with high frequency  and a recent review showed that in 55-70% of all comparisons, native populations outperformed non-natives . Even in populations separated by only small geographic distances, local adaptation has been suggested to contribute to adaptive divergence [13–17]. While these studies suggest the possibility of local adaptation, they could only weakly infer it since no reciprocal transplants were performed. Experiments testing for local adaptation in salmonids in a natural setting are generally rare for practical reasons (dispersal, invasiveness of non-natives, recapture difficulties). Thus, most experiments are carried out in common garden settings in hatcheries or in the laboratory [15, 16]. The conclusiveness of these experiments is however limited, since the key ecological factors driving local adaptation may not be replicated in the common garden environment and domestication selection may promote genotypes performing well under laboratory or hatchery conditions but not in the wild. To our knowledge, only a few reciprocal transplant experiments in salmonids have investigated fitness (or some proxy of fitness) and identified traits under selection references in 8, [18–21], and only three studies used the populations’ natural habitat as an egg rearing environment [22–24] even though selection in the wild might be substantial during egg incubation and fry emergence mortality may approach 90%; . Considering the growing concern in recent years about the adaptability of temperature-dependent life history traits such as embryonic developmental rates and timing of hatching e.g. , local adaptation of salmonids during the embryonic stage has not been sufficiently understood.
We made two predictions with regard to survival: 1) within sites: native populations outperform non-natives; 2) across sites: populations show higher performance at ‘home’ compared to ‘away’ sites. Generally, the first criterion is a better diagnostic for the presence of local adaptation because it is not confounded with differences in intrinsic habitat quality . The first criterion directly tests for the presence of divergent selection within habitats, acting on genetic differences in relative fitness, where non-adapted immigrant genotypes are selected against. The second criterion would only be reliable if all test habitats were of the same quality, which is not usually the case (some habitats may generally cause lower or higher fitness in all populations).
Location- and population-specific survival rates
4.5 ± 0.2
4.7 ± 0.3
4.9 ± 0.2
5.0 ± 0.3
5.1 ± 0.2
Likelihood ratio tests comparing generalized linear mixed models
P, L, P×L
P, L, P×L
P, L, P×L
P, L, P×L
P, L, P×L
P, L, P×L
In agreement with the above, none of the five separate GLMMs within location (used to test prediction 1) found a native population to outperform non-native populations within location. Similarly, the five GLMMs within population (used to test prediction 2) did not find populations to generally perform better at their ‘home’ than at ‘away’ sites.
We conclude that i) there is no evidence for local adaptation in survival overall, and ii) that locations differed considerably in their quality for egg incubation, causing responses in embryo survival independent of their population origin.
To compensate for different rates of development caused by temperature differences between streams, we estimated the time when 50% hatching should have occurred for retrieval, taking into account local temperature regimes. However, at two sites hatching had not occurred yet (in Dorfbach and Biberenbach embryos were still alive but had not hatched yet; Table 1 and Figure 2b), and hatching rate at the Kiese was based on only 19 survivors. We therefore excluded these locations and used the two remaining streams (Giesse and Worble) for the analysis of hatching timing. Mean hatching rate across these two sites was 38.8% (Figure 2b). Population egg size was not related to the timing of hatching (F 1,5 = 1.37, p = 0.33).
The full GLMM explained significant amounts of variance in hatching time (Table 2b). No overdispersion was found. Again, Population × Location did not result in a significantly worse or better fit, i.e. there is no evidence for local differentiation in the timing of hatching. Also, the vertical position of the egg within the capsule did not have a significant effect on the timing of hatching of the embryos. Removal of Population, Location, and Capsule on the other hand decreased the goodness of fit significantly and were retained. Tukey-Kramer comparisons between locations (based on the GLMM without the interaction term and Position) showed that the timing of hatching differed significantly between the two burial sites (Z = −2.94, P = 0.003), but only between two of the five populations (Dorfbach vs. Biberenbach: Z = −3.69, P = 0.002, and Worble vs. Dorfbach: Z = 3.34, P = 0.008; the comparison between Kiese and Dorfbach was marginally significant: Z = 2.72, P = 0.052).
We tested for local adaptation in early life-history traits in brown trout using a 5 × 5 population-by-habitat transplant experiment. Embryos were raised in the natural streambed at ‘home’ and at four ‘away’ sites until 50% hatching was expected to have occurred. We tested for both the ‘local vs. foreign’ criterion (native populations show higher fitness than non-natives) and the ‘home vs. away’ criterion (populations have higher fitness in their own habitat than in other habitats) .
We found neither of the two predictions confirmed. The response in embryo survival to environmental conditions was largely independent of population origin, i.e. we found no evidence for local adaptation at the embryonic stage for the trait examined.
There are several possible explanations for the absence of local adaptation. For one, the environments we chose (and the populations sampled) may have not been sufficiently different from each other. Speaking against this, however, our results revealed large differences in survival between sites, ranging from 4% to almost 90%. Since we applied a ‘unique environment’ approach , we only assumed that every habitat was distinct and different in its ecological parameters, and no link was established between the prevailing ecological conditions and the phenotypes they select for. Yet, the survival data suggests that sites indeed differed in habitat quality parameters that are crucial to brown trout, such as chemophysical (e.g. temperature, oxygen availability, rate of discharge) and ecological characteristics (e.g. substrate grain size, degree of shading) [27, 28]. Alternatively, the microclimates at the burial sites may not have been representative of the habitat conditions typically encountered by the native population and the egg containers we used may not sufficiently simulate the conditions in a natural redd. Furthermore, if mothers exhibited preferences for particular characteristics of their redd, the use of artificial nests could miss aspects of offspring adaptation to the specific incubation conditions locally adapted mothers usually create for them.
In one case (Kiese), conditions seemed unsuitable for egg incubation altogether. Within-stream replication of burial sites would have helped to distinguish between a genuine absence of local adaptation and microclimatic effects. Arguing for the importance of microclimatic conditions during egg incubation, we found strong effects of the individual egg capsules, i.e. the positioning even within one square meter of riverbed made a difference to survival. Also the vertical positioning within the capsule was important, such that eggs closer to the surface had a higher chance of survival than eggs positioned deeper in the riverbed.
There are many environmental threats to an embryo developing in the riverbed. Hatching has been shown to be an inducible defense mechanism in salmonid fishes, such that populations inhabiting risky habitats may hatch earlier or later to escape suboptimal conditions [29, 30]. Variation in the timing of hatching can be caused, for example, by desiccation, oxygen shortage, pathogens, or predation. Alternatively, selection may target female spawning time, which can also affect the emergence time of the fry. Whatever the selective mechanism, if populations are genetically adapted to their home site, and if sites differ in selection regimes, differences in the timing of hatching between populations within sites can be expected. Indeed, timing of hatching differed significantly between locations, but all five populations responded similarly to the conditions prevailing at these sites, indicating that, in this population network, the timing of hatching depends largely on the local environmental conditions rather than on population differences.
There are more scenarios that could explain the absence of adaptation in survival and the timing of hatching. For one, maternal investment can potentially confound effects of local adaptation [2, 8]. Egg size provides a proxy for the quality of maternal investment (larger eggs contain more yolk) and has been shown to affect growth and survival . Although we could not statistically control for within-population variation in maternal effects since maternity was not tracked until retrieval, we used large numbers of randomly selected dams per population, which should help to obtain reliable estimate of average population fitness. Although populations differed in average egg size, we found no evidence for a direct association between mean egg size and survival or hatching, suggesting that differential maternal investment is not responsible for the absence of local adaptation here.
Both, high levels of gene flow and temporally varying selection can prevent local adaptation [32, 33]. However, high gene flow is unlikely in this system. A recent study on brown trout in an area overlapping with our study area found substantial genotypic and phenotypic differences in locomotory and trophic morphology between geographically close populations (2-40 km), arguing for the presence of local adaptation or phenotypically plastic response . Another recent study on brown trout in Switzerland found evidence for adaptation of brown trout to altitude . Generally, brown trout are known to have spatially restricted ranges and genetic differentiation on small geographic scales has been repeatedly reported [35–37]. One possibility is, that population differences in survival due to local adaptation only become evident at later developmental stages in this system.
Although cantonal stocking policies in this area stipulate that streams can only be stocked with offspring from native individuals (further details in ), domestication selection can reduce the efficacy of divergent selection and local adaptation. Individuals spending their embryonic and juvenile development in hatcheries are not exposed to natural selection at their site of origin, and the alleles beneficial for survival in the wild may never come to fixation in continuously stocked populations. Hansen et al. argued that repeated introduction and admixture of wild and hatchery brown trout in Denmark have led to reduced local adaptation, and hatchery strains of salmonids are often subject to domestication selection  with the result that their fitness in the wild is decreased relative to wild fish [40–42]. Some of the adults, used to generate embryos in this experiment, may have been hatchery-raised individuals released into the streams in previous years. Additional file 1: Table S1 shows the stocking effort in the five study streams over the last 30 years. Interestingly, the Worble was stocked with by far the lowest numbers of hatchery-raised fish. At the same time, embryos from the Worble were the only ones performing better at their ‘home’ site than at three other locations. We agree with Fraser et al. and Meier et al. that more information is needed on how hatchery fish might affect adaptation at local and regional scales.
We conclude that there is a great need to evaluate the prevalence and spatial scales of local adaptation in the wild, and to establish the link between molecular variation and variation in fitness in local adaptation research. More integrative approaches and collaboration between disciplines will help to understand the molecular basis of local adaptation, and how it may vary with developmental stage, selection gradients, geographic scale, and demography.
The five test streams were part of a hatchery’s yearly spawner collection: every autumn, the state hatchery in Reutigen collects adult breeders from different rivers and streams of the Berne canton, and uses them for supportive breeding, i.e. males and females are stripped of their gametes in the hatchery, embryos are produced by artificial insemination, and the fry are released into their streams of origin in the following spring or summer. Geographical distance between the test streams was maximized (within the sample of streams available) to increase the likelihood that ecologically and genetically distinct populations were sampled. A recent study on brown trout in an area overlapping with the area here, showed that populations can differ significantly in both genotype and phenotype on a scale of 2-40 km . Another criterion for stream selection was that the beds of the stream had to be accessible and the substrate suitable for the burial of egg capsules (see below for further details). Egg burial took place in locations where natural spawning sites had been observed in previous years by hatchery staff (pers. comm. Ulrich Gutmann, Fisheries Inspectorate Bern; see GPS coordinates in Figure 1).
In November 2010, sexually mature S. trutta males and females were collected within a 200 m stretch up - and downstream from the five determined egg burial sites (Figure 1). Animals were collected by electrofishing and transported to the hatchery where they were held until further use. Note that for analysis purposes, individuals collected at the same site are assumed to belong to the same ‘population’ although they may represent a sample rather than a biological population.
On December 1st 2010, fish were anesthetized with clove oil. Before fertilization, ten eggs from each female were photographed in a Petri dish on graph paper. Egg diameter was measured digitally in ImageJ 1.44 to the nearest 0.01 mm. Three measurements per egg were taken from five eggs of each female and averaged to obtain estimates for egg size (Table 1).
For burial in the natural riverbeds, a total of 2,500 fertilized eggs (500 eggs per population) were distributed into 250 custom-made egg capsules constructed from stainless steel mesh wire, bolts and washers (Ross-Gillespie et al. in prep). Each capsule consisted of ten equally sized, vertically arranged compartments, each for one egg. The mesh wire allowed for good through flow of water in the stream while keeping embryos separate and identifiable. Eggs were distributed systematically, such that each capsule contained two eggs per population, filling each vertical position within the capsule with equal frequency across populations (i.e. yielding 50 replicate capsules with five different distribution schemes). On December 2nd 2010, capsules were transported in chilled standardized water to the five burial sites (Figure 1) and inserted into the gravel using a customized metal tool (similar to a pottiputki tree planting tube). A total of 50 capsules were buried at each site in rows of five capsules, keeping a distance of ca. 10 cm between capsules and rows. Each of the five burial sites received ten replicate capsules of each of the five distribution schemes, making sure that each population was represented in equal numbers in every capsule and every vertical position at each site. In summary, we allocated 100 eggs of each population to 1) their native spawning site, and 2) to the native spawning sites of the four other populations. Buried eggs remained in the stream until their retrieval between February and March 2011.
Digital data loggers were installed at each site (iLog, ESCORT data loggers Inc., USA), recording water temperature hourly until the end of the experiment. On February 6th 2011, data from each logger were downloaded and loggers were placed back in the streams. Eggs were left untouched. The average stream temperature until February 6th was then used to predict the time until 50% hatching, using the method provided in Elliott and Hurley equation 1a; , with the aim to retrieve all embryos at approximately the same degree day (to be able to compare survival and hatching rates between sites). After including the remaining temperature data (from February 6th to retrieval), capsules from all sites were retrieved at 493 ± 26 degree days (between 69 and 80 actual days) after fertilization (Table 1). Eggs were visually inspected in the field directly after removal from the riverbed and scored for survival and hatching. Individuals alive at the time of retrieval were considered survived (0 = dead, 1 = alive, both unhatched and hatched individuals). Individuals that had fully emerged from their chorion were considered hatched (0 = alive and unhatched embryos, 1 = hatched embryos). The specific causes of death could not be assessed in this study, but viral, bacterial and fungal pathogens as well as oxygen deprivation or pollution can be considered potential agents. In case of death, alevins were mostly visible inside the egg, and often in an early eyed stage.
Overall, 64 embryos (2.6% of 2,500 eggs) were missing at the time of retrieval, i.e. the compartment within the egg capsule was found empty. In two cases (Kiese and Giesse), an entire capsule was not found (Table 1). All missing embryos were excluded from further analysis.
Because embryos could not be constantly monitored in the streams and initial variation in fertilization success between populations can potentially distort measures of population-specific survival, fertilization success was monitored in the laboratory under sterile and temperature-controlled conditions as follows. On December 1st 2010, a subset of 240 eggs of each population was distributed singly into the wells of 24-well cell culture plates, which contained 2 ml per well of sterile, standardized, aerated water. Culture plates were transported in a chilled cooler box to a climate chamber at the University of Lausanne (Switzerland) where the embryos were kept undisturbed at constant 6.5°C. Every 3 days eggs were checked under the stereoscope for embryos from three weeks to seven weeks after fertilization had occurred (Dec 21st 2010 - Jan 18th 2011). If no developing embryo was visible by the end of this period, eggs were scored as not fertilized. Survival rates were recorded once hatching (or death) had occurred.
The conclusiveness of reciprocal transplant experiments increases with the dimension of the population-by-habitat matrix (i.e. the number of populations and habitats tested). Generally, more than two populations are needed to differentiate effects of local adaptation from the effects of drift or migration, causing similar population-by-habitat interactions , and only replication on the population level allows for identifying fitness differences caused by divergent selection. The larger the transplanting matrix, however, the less likely it is that native populations outperform non-natives in every case, and the prediction must be reformulated such that, across all population-by-habitat combinations, mean population fitness should be higher for natives than for non-natives . We thus tested for differences in performance between populations and burial sites applying a generalized linear mixed model (GLMM) using the Laplace approximation with Population (5 levels), Location (5 levels), and the interaction Population × Location as fixed effects. We also tested for the effects of the position of the embryo in the egg capsule (Position ranging from 1 to 10), and capsule number (50 capsules per location), included as random effects. We used a binomial fit for the binary response variables survival and hatching. To evaluate the significance of fixed and random effects, alternative models without the variable of interest were compared to the full model using likelihood ratio tests. If an alternative model had a significantly better fit, this model was subsequently compared against further reduced models. To test for overdispersion, all variables were treated as fixed effects (converting to a GLM for this purpose), and the residual deviance of the model was divided by the number of degrees of freedom (which should approach 1 in case of no overdispersion ). Tukey-Kramer posthoc tests were applied to test for significance of pairwise comparisons between populations and locations using the multcomp package .
To specifically test our first prediction (native populations outperform non-natives within sites) we ran another five GLMMs on survival, separately per location, using Population as fixed effect and Capsule and Position as random effects. To test our second prediction (populations outperform others at ‘home’ but not at ‘away’ sites) we ran five more GLMMs on survival, separately per population, using Location as fixed effect and Capsule and Position as random effects.
ANOVA was used to test for population differences in egg size and for associations between mean population egg size and survival or the timing of hatching. All statistical analyses were performed in R (version 2.9.2; http://www.r-project.org) using the lme4 package .
We thank G Brazzola, JP Danko, M dos Santos, E Clark, M Escher, F Glauser, J Guthruf, R Hari, R Nicolet, T Kawecki, A Ross-Gillespie, L Wilkins, the staff of the Fisheries Inspectorate Bern, especially C Küng, and U Gutmann, M Schmid, and H Walther from the Reutigen hatchery for discussion, support and/or assistance in the field, and two reviewers for valuable comments on the manuscript. This study complied with the relevant ethical regulations imposed by the university, canton and country in which it was carried out. The Swiss National Science Foundation, the Foundation Maison de la Rivière,and a Marie Curie Intra-European Fellowship to RBS provided funding for this study.
This article is published under license to BioMed Central Ltd. 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.