Introduction

The effect of predation on changes in life-history traits is well documented (Roff 1992; Stearns 1992), but we still know relatively little how predation structures the genetic composition of individual prey populations. Life-history theory predicts a shift in resource allocation towards the least vulnerable life-stages (Law 1979; Taylor and Gabriel 1992; Lampert 2006). For example, if small individuals are targeted, selection might favour faster growth to escape the vulnerable size ranges, whereas if large individuals are targeted, earlier maturation at smaller size might be favoured. Selection is based on phenotypes, but if the individual phenotypic differences have genetic basis, then size-selective predation will result in evolutionary change. However, predation can also affect phenotypes directly, through induced defences (Gilbert 1966; Tollrian and Harvell 1999), or by changing the environment, for example resource availability and thereby growth (Claessen et al. 2000; Enberg et al. 2012). Therefore, it may prove difficult to identify evolutionary changes and disentangle them from phenotypically plastic changes.

In this study we take advantage of a semi-natural experiment that occurred when European perch (Perca fluviatilis) was introduced to Lake Myravatn, a much-studied, small (area 0.6 km2) lake in Western Norway. European perch, a planktivorous fish, was in 2008 discovered to have been illegally released (Regmi et al. 2013) to this lake that previously only hosted one piscivorous fish species, northern pike (Esox lucius) (Giske 1986; Kvam and Kleiven 1995; Knudsen et al. 2006). Before the perch introduction, the lake had a characteristic, species-poor zooplankton community where unusually large-sized Daphnia pulex and D. longispina, one calanoid copepod (Eudiaptomus gracilis), and larvae of phantom midge (Chaoborus flavicans) were the dominant species. Chaoborus larvae were the main predators of D. pulex and occurred at high densities (Regmi et al. 2013), exerting high predation pressure on small zooplankton prey (Pastorok 1981; Taylor and Gabriel 1992). The introduction of perch released the predation pressure on juvenile D. pulex by eliminating Chaoborus larvae from the lake, but increased vulnerability of large adults to predation (Taylor and Gabriel 1992). This complete turnover in the nature of size-selective predation, akin to the one theoretically studied by Taylor and Gabriel (1992), gave us the opportunity to study how size-dependent predation affects the life history of D. pulex.

We hypothesized that, first, the perch introduction caused evolution in Daphnia towards faster life history. Life-history changes that could bring about such acceleration include earlier maturation at smaller size, faster growth prior to maturation, and increased reproductive rate coupled to reduced growth after maturation. Second, the offspring quality–size trade-off should shift towards producing a larger number of smaller offspring when predation on the smallest Daphnia is relaxed. To test these hypotheses, we characterized life-history changes that occurred in Daphnia pulex from Lake Myravatn, disentangling possible genetic change from phenotypic change. Importantly, Daphnia clones collected just around the time of the perch introduction were maintained in our laboratory, allowing a direct comparison of pre- and post-introduction clones under common-garden conditions. We measured a comprehensive set of life-history traits: number of instars, growth trajectories before and after maturation, mean age and body size at maturation, probabilistic maturation reaction norm (PMRN) for age and size at maturation, egg size, and fecundity. Because the clones were reared under standardized laboratory conditions over several generations prior to the experiment, differences in these phenotypic traits would represent a genetically-based phenotypic change, i.e. life-history evolution, most likely in response to the shift in main predator from Chaoborus larvae to planktivorous fish, i.e., from negatively to positively size-dependent predation.

Materials and methods

Origin and isolation of clonal lineages

The exact timing of the perch release is unknown. Analysis of age and growth of the introduced perch suggests that it happened between late summer 2006 and spring 2007 (Regmi et al. 2013). The large ecological changes become apparent only in 2008, after the perch had successfully reproduced in summer 2007 (Regmi et al. 2013). Because the clones collected in winter 2006–2007 represent pre-introduction ecological conditions (the number of released perch was probably small), we refer to them as the pre-introduction clones, even though the actual introduction may have happened slightly before. These clones were collected by W. Lampert (Lampert et al. 2010). Seven of the clones came from females collected in December 2006 and six from ephippia (sexually produced dormant embryos), probably produced late in 2006, and collected from the surficial sediments in March 2007. The selected clones represent a random subset of all available clones (11 clones established from adults and 10 clones established from ephippia, see Lampert et al. 2010). Lampert et al. (2010) could not find significant differences in relatedness, genetic differentiation, or juvenile growth between the clones established from females or from ephippia.

Post-introduction clones were established from females collected in January–February 2010, i.e. about 3 years after introduction of perch. Animals were sampled by vertical net hauls (60 µm) from bottom to surface at the deepest part of the lake (18 m). Adult animals were isolated and cultured in separate glasses to create separate lineages (clones). When isolating the clones, they were maintained in aged and filtrated (20 µm and 0.20 µm) lake water at 19 °C and fed the green alga Scenedesmus sp. The clones were genotyped by use of microsatellites by K. Lampert [see Lampert et al. (2010, 2014) for methods for genotyping]. 13 pre-introduction and 14 post-introduction clones, in total 27 clones, were used in the experiment.

Experimental setup

For logistic reasons, the experiment was performed in three separate but procedurally identical batches, with respectively ten, seven and ten clones randomly chosen from the clone pool. Thus, each batch had up to five pre- and five post-introduction clones, with 20–30 individuals per clone, for a total of 630 individuals. Juveniles from third and fourth clutch were used to remove maternal effects (Lampert 1993). Juveniles born within 3-hour intervals were kept in 200 ml glass jars until 45 h old, and then transferred individually into 70 ml glass jars. Animals were maintained in a climate room at 19 °C with a 16:8 light:dark regime throughout the experiment. Water was changed every third day using artificial Daphnia medium (Klüttgen et al. 1994) modified according to Ebert et al. (1998). The animals were fed daily with Scenedesmus sp. grown in batch cultures. Each animal was fed an increasing amount of algal suspension; at 1–2 days, 3 days, 4–6 days, and then throughout the experiment 1, 2, 3 and 4 million algal cells, respectively. Cell concentrations were determined daily with a CASY1 cell counter (Schärfe System GmbH).

For length-at-age measurements, all animals were first photographed at 45 h; we considered juveniles too fragile for handling at an earlier stage. Thereafter, they were photographed daily until maturation, defined as the first instar with eggs visible in the brood chamber. After maturation, animals were photographed after the release of juveniles for the first three clutches. The photos were taken through a Wild stereomicroscope with a Sony DSC-17 camera. Each individual was photographed three times. The length of the animals was measured using the image analysis software ImageJ (Schneider et al. 2012). Measurements were taken from the top of the head to the base of the spine. All three measurements were averaged to obtain a length measure for each animal. When approaching time of maturation, the animals were checked every three hours until all were mature. Juveniles were counted for the first three clutches.

Egg size was measured from the photographs of females bearing their first clutch. Egg size was expressed as egg volume and estimated by assuming spheroid eggs. The volume was then calculated as \(\pi /6\, d^{2} l\), where \(l\) is the length of an egg along its main axis and \(d\) its diameter perpendicular to the main axis. For each female, all eggs that could easily be distinguished were measured (median 4 eggs, range 1–10) and their mean was taken as the egg size of that female. Size at birth was assumed to be equal to the mean egg length over all females.

Statistical analysis

Data were analysed using (generalized) linear mixed models using the R package “lme4” (Bates et al. 2015). Detailed model output is provided in the Supplementary Materials. Egg size-body length and clutch size-body length relationships were treated as allometric relationships and estimated after log-transformation of both variables. In the first case, the null hypothesis is no relationship, corresponding to an allometric exponent of zero. In the second case, the null hypothesis is that clutch size-body size relationship has an allometric exponent equal to three. Assuming that body length–weight relationship is isometric, this corresponds to clutch size being proportional to body weight. Experimental batch was treated as nuisance factor (except for the probabilistic maturation models, where this factor was omitted for convergence reasons), with contrasts set such that the values of other model coefficients apply to the mean level of the three batches (option “contr.sum” in R). “Clone” was treated as a random effect in all analyses. In addition, “individual” was treated as a random effect where repeated measurements from a single individual were used, i.e., when estimating probabilistic maturation reaction norms (Heino et al. 2002). To account for variable duration of the time between instars, log-transformed value of instar duration (in days) was used as an offset variable (Harney et al. 2013). Error distribution in this analysis was assumed to be binomial, whereas normal distribution was used in other analyses. Likelihood ratio testing was used for testing significance of fixed effects.

Results

Age and length at maturation

Our experimental Daphnia pulex reached maturation, here defined as the instar before the release of the first clutch of eggs, after three to five post-45 h instars (Table 1). Because the number of the instar at 45 h of age is unknown, in all that follows we will pragmatically define the instar at 45 h as the first instar of the observation period. Of the 630 tested individuals, 40% matured during the third observed instar, 58% during the fourth one, and the remaining 2% during the fifth one. Individuals from the post-introduction clones had a higher maturation probability during the third observed instar than those collected before (Table 1). Maturation probability was similar between the groups during the fourth observed instar, and all the remaining individuals matured during the fifth instar.

Table 1 Maturation by instar number before and after the fish introduction

The average age and observed instar number at maturation was lower for the post-introduction clones compared to pre-introduction ones, but their length at maturation was marginally larger (Fig. 1). The difference was, however, only significant for age at maturation (number of observed instars, estimate ± SE, d.f. = 1: − 0.23 ± 0.13, P = 0.079; age: − 0.50 ± 0.22 d, P = 0.034; length: 0.025 ± 0.026 mm, P = 0.344). For all three cases, between-clone variability was important: the standard deviation of “clone” as a random effect is of similar magnitude as the effect of the fish introduction (number of instars, SD with 95% confidence limits: 0.32 [0.22…0.40]; age: 0.56 [0.39…0.70] d; length: 0.063 [0.043…0.080] mm).

Fig. 1
figure 1

Instar number (instar at 45 h = 1st observed instar) and age and length at maturation before and after the perch introduction. The boxes show the medians and quartiles while the whiskers show the extremes within 1.5 times the interquartile range. Black circles show raw means across all individuals. The numbers of measured individuals were 299 and 331 before and after the introduction, respectively

Size and growth

Post-introduction clones had larger instar-specific size for all observed instars (Fig. 2a). However, because maturation can occur at any of the instars three to five, the effect of maturation on growth might have confounded these simple estimates. Length differences for the instar at maturation and at clutches one to three were small (Fig. 2b).

Fig. 2
figure 2

Size and growth before and after the perch introduction. a Size-at-age for observed instars one to five (instar at 45 h = 1st observed instar). Only instars up to maturation are included. b Size for instar at maturation (instar prior to release of the first clutch) and instars at clutches one to three. Age was not recorded for clutch two and three. In a, b, the error bars show one standard deviation. Box-and-whisker plots ce show average growth in mm per day for three different intervals, with black circles showing the raw means. In c, zero size at birth is assumed. The number of measured individuals was 630 for size at maturation and the 1st clutch, 625 for 2nd clutch, and 434 for 3rd clutch

Growth rate to maturation was significantly higher after the introduction, both when measured from birth (Fig. 2c: 0.030 ± 0.013 mm/d, d.f. = 1, P = 0.027) and from 45 h (Fig. 2d: 0.029 ± 0.014 mm/d, d.f. = 1, P = 0.044). In contrast, growth rate from maturation to the instar after the release of the first clutch was significantly slower after the introduction (Fig. 2e: − 0.035 ± 0.015 mm/d, d.f. = 1, P = 0.025). This result is not essentially changed if growth rate is square root transformed to obtain residuals that are close to normally distributed (− 0.032 ± 0.014 (mm/d)1/2, d.f. = 1, P = 0.025). No significant differences were found for the growth increment between the release of the first and second clutch (0.007 ± 0.009 mm, d.f. = 1, P = 0.444), nor between the second and third clutch (0.010 ± 0.009 mm, d.f. = 1, P = 0.271).

Variability in growth was large between the clones. For example, the standard deviation of between-clone variability was 0.034 [0.026…0.043] mm/d for growth from 45 h to maturation, and 0.030 [0.019…0.039] (mm/d)1/2 from maturation to the first clutch—similar to the magnitude of the estimated difference in growth between the pre- and post-introduction clones.

Dependence of maturation on age and size

To disentangle whether earlier maturation was caused by faster growth or by lower size threshold for maturation, we estimated maturation conditional to individual age, length, and/or instar number, i.e., a probabilistic maturation reaction norm. From the outset, it is not obvious which of the said variables, or a combination thereof, best can explain maturation. Figure 3 illustrates the dependence of maturation on age, size, and instar, suggesting that maturation is primarily determined by length: both before and after the introduction, maturation probability rapidly increases from very low to very high values at around 1.6–1.7 mm of length. The fast-growing individuals reach this critical size range during the third observed instar, and almost all others during the fourth instar. Before the introduction, a larger proportion of 3rd instar individuals were below this size threshold than afterwards (Fig. 3a, b), in agreement with their lower overall maturation probability (Table 1).

Fig. 3
figure 3

Probability of maturation for all observed instars (top row) and for the third instar (bottom row; instar at 45 h = 1st observed instar) before and after the perch introduction. The grey contour lines enclose 25%, 50%, 75%, and 90% of individuals. Maturation probability (colour gradient) is calculated based on dividing the kernel density estimate of maturing individuals with that of all individuals for the combinations of age and length enclosed by the 75% contour line

The qualitative impression from Fig. 3 is confirmed by models for probability of maturation (Table 2). Length is the single most important explanatory variable. Age is also important and improves the model regardless whether it is measured in terms of time or instar number. Of these two, time improves the model more than instar number (ΔAIC = 5.1). There is no significant difference in maturation of pre- and post-introduction clones (Table 2).

Table 2 Probabilistic models for maturation for the second observed instar and later (instar at 45 h = 1st observed instar)

We focus then on the third observed instar (occurring between age 4 and 5 days), which shows the largest contrast between individuals that “decided” whether to mature or not (Fig. 3). Without yet considering the effects of age or size, the odds of maturation were 7.8 times higher for after-introduction clones than before (estimate ± SE in log(odds): 2.06 ± 1.05, likelihood ratio (LR) test: Χ2 = 3.70, d.f. = 1, P = 0.055). However, considering also the effect of length yields a vastly better model (Χ2 = 263, d.f. = 1, P < 0.001). Adding fish introduction to this model does not improve it (Χ2 = 0.182, d.f. = 1, P = 0.662). The length threshold for 50% maturation probability from the length-only model is 1.67 mm (so-called PMRN midpoint), with the corresponding 25% and 75% maturation thresholds at respectively 1.64 mm and 1.70 mm (i.e., the maturation envelope). There is some indication of a positive effect of within-instar age on maturation (estimate ± SE in log(odds), model with length but without introduction: 1.6 ± 0.93, LR test: Χ2 = 3.17, d.f. = 1, P = 0.075).

Clutch and egg size

Clutch size increased with body length (Fig. 4). For the first clutch, the allometric exponent was not different from three (estimated exponent: 2.79 ± 0.15, d.f. = 1, P = 0.162), suggesting an isometric relationship between fecundity and body weight. There was no significant difference in the clutch size of pre- and post- introduction clones (period-dependent intercept vs. null model, LR test: Χ2 = 0.203, d.f. = 1, P = 0.652). Moreover, while egg size increased with body length (Fig. S1), egg size–body length relationship was similar for pre- and post-introduction clones (LR test: Χ2 = 1.18, d.f. = 1, P = 0.278). Taken together, these results suggest that reproductive effort for the first clutch did not change between the two periods.

Fig. 4
figure 4

The relationship between body length and clutch size before and after the perch introduction. Length refers to the body size for the instar during which the clutch is produced. Regression lines were fitted to log–log transformed data to capture the allometric relationship. The difference in the relationship is significant for the second and third clutch. Body size is rounded to the nearest 0.01 mm, with the area of symbol proportional to the number of observations. The number of measured individuals was 611 for the 1st clutch, 618 for 2nd clutch, and 617 for 3rd clutch

For the second and third clutch, the relationship between clutch size and body size differed between the two periods. While there was no overall difference in clutch size between the periods (period-dependent intercept vs. null model, LR test, 2nd clutch: Χ2 = 1.13, d.f. = 1, P = 0.288; 3rd clutch: Χ2 = 2.06, d.f. = 1, P = 0.151), the model with period-dependent allometric exponent and intercept was significantly better than the null model (2nd clutch: Χ2 = 11.3, d.f. = 2, P = 0.004; 3rd clutch: Χ2 = 11.2, d.f. = 2, P = 0.004). The allometric exponent is less than three for the second and third clutch (estimated exponent without accounting for the introduction, 2nd clutch: 1.91 ± 0.13, P < 0.001; 3rd clutch: 1.61 ± 0.17, P < 0.001). When the allometric exponent is period-dependent, it is lower after the introduction. Thus, these models suggest a length-fecundity relationship with higher fecundity for small- and medium-sized individuals after the introduction, but not for larger ones (Fig. 4b, c). At the mean body length across all individuals, which is 2.17 mm for the second clutch and 2.49 mm for the third one, fecundity is estimated to be respectively 4% and 6% higher after the introduction than before it.

Discussion

The aim of this study was to identify potential genetically-based life-history change in a population of Daphnia pulex as a result of size-dependent predation by newly introduced perch. Because we had maintained pre-introduction clones in the laboratory, without allowing for sexual reproduction, we could use a common-garden experiment to establish whether genetic adaptation had occurred between clones collected at the time of the introduction and 3 years later. Our results reveal genetic changes in growth and clutch size, and a change in maturation that was triggered by the change in growth, thus showing that predation can lead to rapid evolution of faster life history in D. pulex.

Our results show that the post-introduction clones matured earlier than the pre-introduction clones, as predicted by life-history theory (Roff 1992; Stearns 1992). However, the analysis of probabilistic maturation reaction norms (PMRNs, the size- and age-dependent maturation probabilities) did not display a significant change, showing that direct genetic adaptation in maturation was negligible or absent. Instead, growth to maturation was significantly accelerated in the post-introduction clones compared to the pre-introduction clones. This shows that the earlier maturation occurred as an indirect response to genetic change in another trait, in a similar way when maturation plastically responds to environmentally-induced variations in growth (Stearns 1992; Heino et al. 2002): fast-growing individuals reached the size threshold for maturation described by the PMRNs earlier.

There was also a tendency towards a reduction in number of instars to maturation post introduction, but it was not significant. Several studies have shown a decrease in number of instars when Daphnia are exposed to fish (e.g., Beckerman et al. 2010). We did not count instars directly, but deduced their number from the growth curves of each individual based on frequent size measurement. However, as we needed to ensure that the animals survive the handling when photographed, the first pictures, and hence the first size measurements, were taken when the animals were 45 h old; younger animals proved too small and vulnerable for photographing. Instar-specific length and age distributions between the two periods are best aligned by assuming that the instar at 45 h corresponds to the same instar for pre- and post-introduction clones (Fig. S2–S3, see also Fig. 3). The slight shift of the length distributions towards larger sizes in post-introduction clones (Fig. S2) is consistent with an additional pre-45 h instar, but this would have been expected to lead to a parallel shift towards slightly older ages, which is opposite to what we observed (Fig. 2a). We therefore conclude that a systematic change in instar numbers before 45 h is unlikely, but only direct observations could have provided a conclusive answer to this question. At any rate, our main results on growth, maturation, and clutch size are not dependent on instar numbers before maturation.

We could also show a number of evolutionary changes in the adult part of the life cycle. Theoretically, we would have expected changes in energy allocation from growth towards reproduction (Roff 1992; Stearns 1992), possibly in combination with reduced energy acquisition (Enberg et al. 2012). Indeed, growth from maturation to the release of the first clutch was lower for the post-introduction clones compared to the pre-introduction clones. This could represent either reduced energy acquisition or increased allocation of energy to other purposes. We did not measure the former, whereas increased reproductive allocation is an obvious contender for the latter. However, the results did not show an increase in reproductive allocation, at least not in clutch size and egg size, the proxies of energetic costs of reproduction we measured. However, we did observe an increase in the size of the second and the third clutches. Growth for the corresponding time intervals did not change. Thus, we could show the theoretically predicted changes towards lesser adult growth and increased reproduction, albeit only for some of the measured parameters.

Egg size was positively correlated with female size, as observed in earlier studies (Lampert 1993). We could not detect a difference in the egg size between pre- and post-introduction clones, while life-history theory predicts that offspring size should decrease when size-dependent predation pressure during the early life history is reduced (e.g., Jørgensen et al. 2011). Note that we photographed eggs in their mother’s brood pouch, which is a quick but not very precise method; more precise measurements might have revealed the theoretically predicted response.

For logistic reasons, we conducted the experiment in three consecutive but procedurally identical batches. Because we were not interested in the batches per se, we treated batch as a nuisance factor. Nevertheless, there were significant differences among the batches for all traits except clutch size (see the Supplementary Material). We are unable to identify the cause of these differences, but they might be related to the laboratory being relatively new when the experiment was started, manifested as, for example, some variation in the food quality.

Life-history changes in a population are not necessarily manifesting genetically-based phenotypic change. One of the prominent features of Daphnia is their ability to detect and respond to kairomones (chemical signals) released by predators, manifested as changes in morphology, behaviour, and life history (e.g., Spitze 1991; Stibor 1992; Stibor and Lüning 1994; Tollrian and Dodson 1999; Lass and Spaak 2003). The discovery of inducible defences (Gilbert 1966) in zooplankton resulted in a number of experimental studies (for review see Lass and Spaak 2003) aiming to separate changes in phenotype expression from direct selection effects. Importantly, these have demonstrated inducible life-history changes similar to what we have observed in our experiment. Daphnia exposed to water infused with fish kairomones reduce body size and age at maturation, and increase clutch sizes (Stibor 1992; Stibor and Lampert 2000; Beckerman et al. 2010). When exposed to invertebrate kairomones, Daphnia increase body sizes and age at maturation (Spitze 1991, 1992; Beckerman et al. 2010). We obtained Daphnia clones from their natural habitat, where they had been exposed to both direct size-dependent predation and the predator’s chemical impact, during a continuous transition from invertebrate to (predominantly) vertebrate predation. However, maintaining Daphnia as asexually reproducing clones made it possible to take a snapshot of the genetic composition of a population and keep it in laboratory conditions for years. Plastic effects of predator’s chemical influence will be “reset” under common-garden conditions, and hence any differences in clones collected during different size-dependent predator regimes will essentially represent a genetically-based phenotypic change.

The change in size-selection regime in Lake Myravatn has intriguing parallels with the effect of fishing on size-dependent mortality in fishes. In fishes, mortality from predation is typically negatively size-dependent (Lorenzen 1996; Sogard 1997; Perez and Munch 2010), reflecting factors such as gape limitation and other restrictions in handling of large prey (Sogard 1997; Scharf et al. 2000) as well as better predator-evasion abilities among large-sized prey (Miller et al. 1988; Lundvall et al. 1999; Schürch and Taborsky 2005). The situation for Daphnia in lakes without planktivorous fish but with Chaoborus larvae is broadly similar.

Size-dependent mortality from fishing has a pattern that is almost exactly the opposite to natural mortality: fisheries typically target large-sized fish (Fenberg and Roy 2008; Allendorf and Hard 2009), leading fishing to be the main source of mortality for adult fish (Mertz and Myers 1998). Thus, fishing turns around the mortality regime in a way similar to what the perch introduction did in Lake Myravatn, with similar expected life-history consequences. Trends toward earlier age at maturation are commonplace in exploited fish populations (Trippel 1995), and analyses of phenotypic data using the PMRN approach indicate that such trends often have a genetic component (reviewed by Heino et al. 2015). In contrast, while we could demonstrate changes in both juvenile and adult growth, relatively few fisheries studies have addressed growth changes, and even fewer have separated juvenile and adult growth (Enberg et al. 2012; Heino et al. 2015). Similarly, changes in measures of reproductive investment such as clutch size remain little studied in exploited fish (Heino et al. 2015).

Our results with Daphnia are in agreement with the fishing-induced evolution literature in showing earlier maturation as a prominent response to positively size-dependent mortality (Audzijonyte et al. 2013; Heino et al. 2015). However, whereas most fish studies have shown that size threshold for maturation has declined, indicating a genetic change in maturation conducive to earlier maturation, we did not observe this. Instead, earlier maturation in Daphnia was a purely secondary, indirect response to faster juvenile growth, similar to phenotypic plasticity, except that the cause of faster growth here was genetic, not environmental. Evolution of faster growth in response to positively size-dependent mortality is known from theoretical models, but has not been empirically documented in exploited fish (Dunlop et al. 2009; Enberg et al. 2012).

The distinct advantage we have over observational studies on fishing-induced evolution is that we have been able to unequivocally show that organisms exposed to positively size-dependent predation genetically differed from those that had been exposed to negatively size-dependent predation. However, because we lack replication and controls, we share the same impediment as the field-based fishing-induced evolution studies in that we cannot unequivocally identify the driver of the change, that is, change in mortality, as opposed to some other change coinciding with the period of observation.

Adaptive changes in life histories, be they phenotypically plastic or genetic, will help populations to cope with changing environments (Chevin et al. 2010). Clearly, phenotypic plasticity is often advantageous because it allows rapid adjustments over time scales of individual life spans. Nevertheless, plasticity is also costly and therefore has its limits (DeWitt et al. 1998; Ernande and Dieckmann 2004). Large and persistent changes in the environment, beyond what a population has experienced in its recent evolutionary past, are likely to exceed the capacity of plasticity to buffer environmental change. The scenario in which a population’s continued existence is only made possible through adaptive evolution has become to be known as “evolutionary rescue” (Bürger and Lynch 1995; Gomulkiewicz and Holt 1995; Gonzalez et al. 2013; Carlson et al. 2014). Alas, in the case of Daphnia pulex from Lake Myravatn, despite showing significant genetic changes in the life history, the environmental challenge exceeded the population’s adaptive capacity: the population went extinct in about 2010/2011 (unpublished), not long after the post-introduction clones for this study were collected.