- Research article
Behavioural and physiological adaptations to low-temperature environments in the common frog, Rana temporaria
BMC Evolutionary Biologyvolume 14, Article number: 110 (2014)
Extreme environments can impose strong ecological and evolutionary pressures at a local level. Ectotherms are particularly sensitive to low-temperature environments, which can result in a reduced activity period, slowed physiological processes and increased exposure to sub-zero temperatures. The aim of this study was to assess the behavioural and physiological responses that facilitate survival in low-temperature environments. In particular, we asked: 1) do high-altitude common frog (Rana temporaria) adults extend the time available for larval growth by breeding at lower temperatures than low-altitude individuals?; and 2) do tadpoles sampled from high-altitude sites differ physiologically from those from low-altitude sites, in terms of routine metabolic rate (RMR) and freeze tolerance? Breeding date was assessed as the first day of spawn observation and local temperature recorded for five, paired high- and low-altitude R. temporaria breeding sites in Scotland. Spawn was collected and tadpoles raised in a common laboratory environment, where RMR was measured as oxygen consumed using a closed respiratory tube system. Freeze tolerance was measured as survival following slow cooling to the point when all container water had frozen.
We found that breeding did not occur below 5°C at any site and there was no significant relationship between breeding temperature and altitude, leading to a delay in spawning of five days for every 100 m increase in altitude. The relationship between altitude and RMR varied by mountain but was lower for individuals sampled from high- than low-altitude sites within the three mountains with the highest high-altitude sites (≥900 m). In contrast, individuals sampled from low-altitudes survived freezing significantly better than those from high-altitudes, across all mountains.
Our results suggest that adults at high-altitude do not show behavioural adaptations in terms of breeding at lower temperatures. However, tadpoles appear to have the potential to adapt physiologically to surviving at high-altitude via reduced RMR but without an increase in freeze tolerance. Therefore, survival at high-altitude may be facilitated by physiological mechanisms that permit faster growth rates, allowing completion of larval development within a shorter time period, alleviating the need for adaptations that extend the time available for larval growth.
Stressful environments (environments outside the optimum conditions for a particular species) can impose strong ecological and evolutionary pressures at a local level [1, 2]. Population persistence depends on the ability of individuals to respond to environmental stress through adaptive, plastic or behavioural mechanisms that maximise fitness . Extremes of pH (common frog; ), water availability (wild mustard ), and temperature (redband trout; ) have been observed to drive adaptive population divergence. High-latitudes and altitudes experience low temperatures that can result in shorter activity periods and longer periods of freezing [3, 6, 7]. Plastic and adaptive responses to low temperature environments have been widely recorded (for a review see ) and can result in cryptic divergence between populations inhabiting different temperature regimes (counter-gradient variation; ). Temperature is often the major abiotic factor that influences physiological mechanisms in ectotherms [10, 11] and growth slows in response to cold environments . Reduced activity periods in low-temperature environments, in combination with low-temperature driven growth-rate reductions, can result in lower sizes at important life-history events such as metamorphosis and reproduction . Smaller sizes can translate to lower fitness when weight is positively correlated with survival or reproductive success [13, 14].
Assessing the mechanisms that facilitate survival in challenging environments is important for understanding how populations respond to ecological and evolutionary pressures, particularly in a globally changing climate . Potential responses to maximise size at important life history stages in low-temperature environments include altering metabolic rate (e.g. to allow more resources to be allocated to growth; [6, 16]), developmental period (e.g. delaying sexual maturity; [17–19]), or temperature activity range (e.g. breeding at lower temperatures; ). Populations that inhabit high-altitude environments experience lower temperatures and shorter activity periods than their low-altitude neighbours and offer an excellent opportunity to assess how survival is facilitated in environments where growth is constrained . Amphibians are a particularly good model for studying physiological and behavioural responses to growth constraints, as size at metamorphosis is positively correlated with survival in the subsequent terrestrial life-history stages [22, 23].
Variation in metabolic rates between individuals is a common occurrence in nature , but the effects on fitness are still relatively unknown . Resting metabolic rate, here defined as the energetic cost of self-maintenance , has been linked to multiple physiological and behavioural traits including predator avoidance, foraging behaviour, swimming performance and growth [13, 24]. Growth imposes a significant physiological cost and can result in a trade-off with other physiological mechanisms [13, 27, 28], especially when resources are limited . As the vast majority of energy expenditure in ectotherms is maintenance costs (80-85%), small differences in resting metabolic rate can result in large differences in energy available for growth . An increased growth rate can result in a larger size at important life history events and has been linked to a reduced resting metabolic rate in sagebrush lizards at high-altitude , Sydney rock oysters from growth rate-selected stock  and snapping turtles . Attempts to assess the physiological trade-offs facilitating higher growth rates in larval common frog (Rana temporaria) at high latitudes have found no link to reduced metabolic rates . However, as temperature is not linearly related to latitude in Sweden , where these experiments were conducted, these results may mask the true nature of the temperature-metabolic rate relationship. Therefore, further research in a system with a linear temperature change is required to elucidate the relationship between resting metabolic rate, growth rate and temperature.
Another potential response to maximise size at important life-history events, is to increase the time available for growth prior to metamorphosis or reproduction by extending development over multiple growth periods . The concept of delayed development, or diapause, has been commonly observed in insects, often in terms of cohort splitting where different cohorts within a population complete development at different times of the year, or even in different years . In amphibians, the period immediately prior to metamorphic climax is accompanied by a loss of weight , but a lower weight decreases the chances of adult survival . Therefore, overwintering at a higher weight, but still at the larval stage, and metamorphosing the following year has the potential to increase survival, and has been recorded in a number of temperate amphibian species [19, 32]. However, low winter temperatures at high-altitude can lead to prolonged periods of freezing . Therefore, in order to survive, overwintering larval amphibians must be able to respond to freezing temperatures via freeze avoidance (i.e. inhabiting environments that buffer individuals from freezing temperatures) or freeze tolerance (survival of extensive freezing of body fluids; [3, 34, 35]). Freeze tolerance depends on the ability to restrict ice formation to extra-cellular areas, which is mediated by accumulation of low molecular weight carbohydrates in the blood [3, 34]. The ability to tolerate freezing has been linked to glucose accumulation in the blood, via release of liver glycogen, in the frogs R. sylvatica, R. lessonae and R. esculenta [3, 34], and with glycerol accumulation in Hyla versicolor . However, all previous studies have focussed on freeze tolerance in adult amphibians and the potential for freeze survival in the larval stage has, to the best of our knowledge, never been studied. The ability of a tadpole to survive freezing would extend the time available for growth to, and thus size at, metamorphosis in amphibians breeding in temperate climates.
A third alternative response to larval growth constraints would be for adults to adapt behaviourally rather than amphibian larvae adapting physiologically. Adults have the potential to expand the growing season for larvae by breeding earlier in the year . In temperate amphibians, breeding is closely linked to temperature  and frequently occurs immediately after winter dormancy (e.g. Bufo bufo, R. chinensis, R. sylvatica and R. temporaria; [38, 39]). By adults becoming active and breeding at lower temperatures, larvae would have longer to grow and develop prior to winter dormancy. The longer time available for growth would allow larvae to reach a larger size at metamorphosis and thus have an increased chance of survival as adults .
The common frog (R. temporaria) is the most widespread amphibian in Europe and occurs from zero to 2742 metres above sea level within its range, and to over a thousand metres on the mountains of Scotland [40–42]. It is an explosive breeder, with communal spawning taking place immediately after winter dormancy ; a 5°C temperature threshold is generally considered to initiate activity and spawning . R. temporaria larvae show increased growth rates in response to low temperatures experienced at high-latitudes and altitudes throughout its range [16, 45]. We have previously shown that local adaptation to high-altitude environments occurs even in the face of high gene flow, suggesting that temperature exerts a strong selective pressure . However, there are also reports of R. temporaria overwintering as tadpoles in Scotland, although it is currently unclear whether this response is particularly linked to low-temperature environments . The mountains of Scotland offer an excellent opportunity to study the responses that facilitate survival in low-temperature environments, as there is continuous habitat along altitudinal gradients, with temperature decreasing linearly by 0.65°C for every 100 m gain in altitude [40, 46, 47]. Individuals from high-altitude sites in Scotland experience substantially lower temperatures than their low-altitude counterparts, with an average mean annual temperature reduction of 4.5°C at high- compared to low-altitude breeding sites .
The overall aim of this study was to assess the physiological and behavioural responses of common frogs in Scotland that facilitate survival in low-temperature environments. In particular, this study answers the questions: 1) do high-altitude adults extend the time available for larval growth by breeding at lower temperatures than low-altitude individuals?; and 2) do tadpoles sampled from high-altitude sites differ physiologically from those from low-altitude sites, in terms of routine metabolic rate and freeze tolerance?
Adult spawning behaviour in relation to altitude
Temperature on the day of egg mass observation was, on average, 7.5 ± 2.1°C (Table 1) and did not vary predictably with altitude (r2 = -0.03, p = 0.41). Likewise, no significant regression with altitude was found for the average temperature in the week prior to egg mass observation (mean temperature = 4.8 ± 0.9°C; r2 = -0.06, p = 0.49). Degree days prior to egg mass collection was highly variable across sites (24.5 ± 19.1; Table 1) but also did not show a significant relationship with altitude (r2 = 0.26, p = 0.09). The date of egg mass collection was on average 30 days later at high- compared to low-altitude sites (Table 1) and Julian spawning day showed a significant positive relationship with altitude (r2 = 0.80, p < 0.01): individuals spawned 5 days later for every 100 m gain in altitude (Figure 1). The daily mean temperature at all sites had exceeded the threshold value of 5°C in the week prior to egg mass collection (Figure 2).
Larval physiology in relation to altitude
Routine metabolic rate
Between eight and 20 individuals per site were measured for RMR, due to varying levels of mortality (Mean = 16 ± 5; Table 2). Mean RMR per site varied between 0.02 ml O2 g-1 h-1 (LOMHIGH) and 0.10 ml O2 g-1 h-1 (DUBLOW), with an overall average of 0.07 ± 0.02 ml O2 g-1 h-1 (Table 2). Mountain, altitude, and their interaction were found to be significant in predicting RMR. A Tukey’s HSD test showed a significant difference between high- and low-altitude RMR in individuals from three of the mountains: DUB (diff = 0.03, p < 0.01), MNT (diff = 0.02, p = 0.03) and LOM (diff = -0.07, p < 0.01) (Table 2). The difference between high- and low-altitude RMR was not significant for IME (diff = -0.01, p = 0.49) and LAW (diff = 0.01, p = 0.94). The direction of the relationship varied between mountains, with individuals from DUB, LAW and MNT showing a trend for lower RMR at high- compared to low-altitude, whereas individuals from IME and LOM had higher RMR at high-altitude (Figure 3). The post hoc power analysis of the ANOVA used in the Tukey’s HSD revealed an effect size of 0.29 giving an achieved power of 0.69 to determine a significant difference between the means of RMR by altitude within each mountain. The power to confidently conclude that no significant interactions have been missed (a type II error), is generally set at 0.8 .
Ten individuals per site were tested for freeze tolerance, except for LOMHIGH (seven individuals) and DUBLOW (zero individuals) (Table 2), due to variable tadpole mortality prior to the experiment. Between 10% (LOMLOW, MNTLOW and LAWLOW) and 80% (IMEHIGH) mortality was observed post-freezing across sites (Mean survival = 0.60 ± 0.28; Table 2). Out of 50 tadpoles tested for freeze survival from low-altitude sites, 41 survived (82%), compared with 12/37 (32%) tadpoles from high-altitude sites (Figure 4). Altitude and weight, but not their interaction, significantly changed the log likelihood when removed from the GLMM and were thus included in the final model. There was no significant effect on the log likelihood of the model when mountain was removed from the model and thus mountain was not included in the final model (Additional file 1). The results of the GLMM using the final model (RMR ~ altitude + weight) showed that individuals from low-altitude sites had significantly higher survival than those from high-altitude sites (z = 4.20, p < 0.01).
Adult spawning behaviour in relation to altitude
No significant relationship was found between spawning date and either the temperatures recorded on the day of, or in the week prior to, egg mass observation that would have suggested that breeding occurs at lower temperatures at high altitude sites. Degree days also did not show a significant relationship with altitude but they were highly variable across sites (Table 1), suggesting that degree days are not an accurate predictor of spawning activity in R. temporaria. All sites had exceeded the 5°C temperature threshold generally thought to initiate activity and breeding in R. temporaria  in the week prior to spawning (Figure 2). Our results therefore support 5°C as the activity threshold for R. temporaria regardless of altitude of breeding site, and demonstrate that high-altitude individuals experience a longer period of low-temperatures and delayed spawning compared to low-altitude individuals. The date of spawning was 5 days later for every 100 m gain in altitude (Figure 1). As mean annual temperature decreases by 0.65°C for every 100 m increase in altitude in this system , we can infer that spawning is on average one day earlier for every 0.65°C increase in mean annual temperature. However, further information regarding the date of spawning between 223 m and 720 m (currently not available for this system) is needed to confirm this pattern. It is interesting to note that for the mountains with the highest high-altitude sites (>900 m: DUB, LAW and MNT), the average time difference between spawning at low-and high-altitude is only 24 days, compared with 38 days observed difference between spawning at low- and high-altitude in IME and LOM (high-altitude sites at 703 m and 720 m, respectively; Tables 1 and 3). This could suggest that high-altitude individuals at the highest breeding sites are breeding at a lower temperature. However, this is not reflected as a significant difference between high- and low-altitude spawning temperatures in these mountains (Table 1). The longer period of low temperatures at high- compared to low-altitudes prior to spawning supports the hypothesis that high-altitude individuals experience a shorter annual activity period but does not support the hypothesis that breeding occurs at lower temperatures at high altitudes to provide a longer developmental period, facilitating survival at high-altitude within R. temporaria.
Breeding at a lower temperature would allow individuals to spawn earlier in the year, thus providing their offspring with a longer period in which to develop prior to metamorphosis , such as has been found for amphibian species at higher latitudes compared to their lower latitude counterparts  but, to the best of our knowledge, has not yet been observed within-species. However, phenological studies that have quantified within-species breeding temperature have typically used local weather station data [49, 50] and so the spatial scale might not be not fine enough to represent local conditions in mountain areas, where temperatures can vary rapidly over short geographical distances [21, 47, 51]. Therefore, such studies are liable to miss within-species differences in spawning temperature in relation to altitude. When we quantified temperature at a local level, our results support the absence of a shift in breeding temperature in low-temperature environments, contributing to knowledge of the fundamental niche of R. temporaria. However, we did find that R. temporaria can substantially alter the date they breed in response to the environmental conditions experienced, even between geographically close breeding sites, demonstrating a high degree of plasticity in terms of breeding timing.
Larval physiology in relation to altitude
Routine metabolic rate
There was a significant interaction between mountain and altitude, which complicated interpretation of RMR. Posthoc tests indicated no significant differences in RMR between tadpoles sampled from high and low sites on IME or LAW. In regards to the other three mountains, there were significant differences but not always in the same direction: for DUB and MNT tadpoles sampled from high elevation showed a decreased RMR relative to those from low elevation, but the largest difference due to elevation was at LOM, where tadpoles from high elevation showed a substantial increase in RMR compared to those from low elevation (Figure 3; Table 2). However, due to the moderate power of the model (power = 0.69), it is possible that some significant relationships were missed. It has been suggested that a lower resting metabolic rate can allow more energy to be allocated to growth in resource-limited environments  and a link between lower RMR and increased growth rate has been found in the southern toad , Sydney rock oyster  and snapping turtle . Furthermore, Sears  found that both increased growth rates and reduced RMR were positively correlated with altitude in sagebrush lizards. The three mountains in our study system where RMR was lower in individuals from high-altitudes (although the difference was very small for individuals from LAW) have been shown to be locally adapted to temperature parameters, with larval period decreasing and growth rate increasing at high-altitude . Therefore, the lower RMR of individuals from high-altitude from DUB, MNT and LAW, is in line with the increased growth rates observed at these sites. Lindgren and Laurila  did not find a link between growth rates and RMR in R. temporaria along a latitudinal gradient in Sweden. Therefore, this is the first tentative evidence of reduced RMR being linked to increased growth rate as an adaptation to low-temperatures in an amphibian. However, the significantly higher RMR at high- than low-altitude in LOM (Figure 3), and the non-significant difference between LAWHIGH and LAWLOW in terms of RMR, could suggest that local conditions other than altitude are also important in driving divergence in RMR. For instance, the lowest RMR were observed in individuals from IMELOW and LOMLOW, where spawning occurred at the lowest temperatures of all the sites (apart from LAWHIGH), and this could suggest that low RMR is beneficial in tolerating low temperatures at all altitudes, particularly early in the season when temperatures can fall after spawning. It has also been speculated that increased growth rate can be a result of increased time spent foraging, facilitated by lower predator presence in low-temperature environments, and is unrelated to RMR [52–54]. Furthermore, RMR is influenced by temperature and there was up to a 3°C difference in temperature whilst measuring RMR of individuals from different sites, as well as a difference between the temperature at which the tadpoles were raised (15°C) and the temperature at which RMR was measured (19-22°C). Therefore, further research is needed to assess whether site-specific differences in RMR remain constant through time and whether other mountains with breeding sites of above 900 m also show reduced metabolic rate at high- vs. low-altitude.
Just over half of the tadpoles that were frozen survived (61%; Table 2) and altitude was significant in predicting freeze survival: individuals sampled from low-altitude survived freezing significantly better than those from high-altitude (z = 4.20, p < 0.01). Voituron et al.  suggested that R. temporaria adults were freeze intolerant, as 100% mortality was observed after eight hours of complete bodily fluid freezing. However, Pasanen and Karhapää  found that R. temporaria adults could survive 24 hours in a sub-zero environment but died within three days (the actual period an individual was frozen was not measured in that study). Our results suggest that tadpoles of R. temporaria are also capable of surviving short periods of freezing. This is the first time, to the best of our knowledge, that larval freeze tolerance has been demonstrated in any amphibian.
Although the results presented here suggest that tadpoles are capable of surviving being frozen, the finding of a greater survival of low- compared to high-altitude individuals appears counterintuitive, given the longer period of sub-zero temperatures at high- than low-altitude in Scotland  and the sub-zero average winter temperature of -2.2°C at high-altitude, compared to 1.0°C at low-altitude at the breeding sites used in this study . Indeed an increase in freeze tolerance with altitude has been found in plants (Arabidopsis thaliana ) and insects [18, 57], but freeze tolerance and altitude has not been explicitly linked in herptiles. However, evolution of freeze tolerance in frogs, lizards and turtles has been linked to ecological pressures relating to winter temperatures experienced [35, 58, 59]. In general, freeze tolerant species are those that terrestrially overwinter in sub-zero temperatures, as opposed to avoiding freezing in deep water bodies . Therefore, a higher freeze tolerance would suggest that low-altitude tadpoles are more often exposed to freezing temperatures, whereas high-altitude tadpoles may avoid such exposure altogether by inhabiting deep water bodies or metamorphosing prior to winter. Formation of deep water pools is inhibited at high-altitudes in Scotland due to the rocky, exposed landscape . Therefore, freeze exposure is more likely avoided in high-altitude tadpoles by metamorphosing within a single active season, facilitated by a faster growth rate in conjunction with a lower RMR (, this study). Therefore, our results potentially suggest that of the tadpoles found overwintering in Scotland , low-altitude individuals are more likely to overwinter as tadpoles than high-altitude individuals. However, it is possible that long periods of snow cover at high-altitude  could insulate tadpoles from freezing temperatures even in shallow water bodies. Therefore, further field research is needed to assess whether fewer, if any, high- than low-altitude individuals overwinter as tadpoles.
R. temporaria adults do not show behavioural adaptations in terms of breeding at lower temperatures, instead they delay spawning based on the temperature experienced. However, R. temporaria larvae appear to have the potential to physiologically adapt to low-temperature environments, although these relationships, and the results we obtained for routine metabolic rate and freeze tolerance, are not always straightforward to interpret. Therefore, our results suggest that survival at high-altitude may be facilitated by physiological mechanisms that permit faster growth rates, allowing completion of larval development within a shorter period of time, alleviating the need for adaptations that extend the time available for larval growth. How individuals respond to environmental temperature at a local level is an important step in relating ecological and evolutionary pressures to phenotypes.
Adult spawning behaviour in relation to altitude
Paired high- (above 700 m; R. temporaria occur to over 1000 m in Scotland ) and low-altitude (below 300 m) sites from five mountains within west central Scotland were selected for study (Table 3). The study sites have low neutral genetic population structuring (high between-site gene flow) and show local adaptation of larval traits in relation to temperature [40, 45]. Air temperatures were recorded every two hours at each site between March 2010 and October 2011 using Thermocron i-buttons (Dallas Semiconductor/Maxim, London) and downloaded to a laptop every six months using a USB i-button adapter (Dallas Semiconductor/Maxim, London) and the software, Thermodata viewer (Thermodata pty Ltd., Melbourne) (as per [40, 45]). Sites were visited from early February 2011 and date of spawning was recorded as the day egg masses were first observed at each site. The daily mean temperature on the day of egg mass observation was calculated for all sites. Although the majority of egg masses were at or below Gosner stage 10 on this date, and thus likely to have been laid no more than 100 hours previously , it is possible that spawning activity started prior to this. Therefore, the daily mean temperatures and the overall average for the week prior to egg mass observation were also calculated. In addition, degree days to egg observation were calculated for each site (a statistic commonly used to predict flowering date in plants; ), using the following approach: From the 1st January 2011, degrees above the threshold for development (set at 5°C) were calculated per day using the formula: ((daily maximum temperature + daily minimum temperature)/2)-threshold value. The resulting values were summed to give the total degree days .
Linear regression models were used to assess whether the temperature and date at which individuals spawned had a significant relationship with altitude (continuous variable; m) using: 1) the date parameter, Julian day; and 2) the temperature parameters: mean temperature on the date spawn were observed, the mean temperature for the week prior to spawn observation, and degree days over the threshold 5°C.
Larval physiology in relation to altitude
Ten R. temporaria egg masses were collected soon after laying (Gosner stage 10 or below; ) from each of the sites monitored for adult breeding phenology. Egg masses were collected during the 2011 breeding season, transferred to the laboratory and maintained in individual sterilised water tanks until hatching at 10°C. At hatching (Gosner stage 22), a randomly selected subset of ten tadpoles were removed from each egg mass and placed in groups of five in two individual 1.3 L plastic baskets with a 0.1 cm mesh. Baskets were placed in large tanks in a common 15°C treatment room. Water quality was maintained using a flow-through system and tadpoles were fed ad libitum with a 1:2 mixture of finely ground dried fish and rabbit food (for further details see ).
Routine metabolic rate
Tadpoles were allowed to develop until hind leg toe differentiation became apparent, in the early stages of metamorphosis (Gosner stages 36–39). Twenty individuals per site (one individual per basket*two baskets per family*ten families per site) were transferred to individual containers and allowed to acclimatise to laboratory conditions in 100% oxygenated water for an hour prior to commencement of experimental procedures. After this period, tadpoles were moved into 8 ml respiration tubes filled with 100% oxygenated water and the lids immediately sealed. Tubes were placed in a dim, quiet location to reduce disturbance during the experiment. Respiration tubes remained closed for one hour. At the end of this period, the lid was removed and the oxygen saturation of the water was measured using an oxygen meter and probe (Strathkelvin Instruments, UK). The same procedure was carried out using a control tube containing no tadpole, to account for any oxygen consumption caused by microbial action. The oxygen meter was calibrated prior to each use using 100% and 0% oxygenated distilled water as standards. Distilled water was fully oxygenated using an aquatic bubbler and fully deoxygenated by adding sodium sulphite . The oxygen probe was maintained at a constant temperature to avoid biases caused by thermal fluctuation using a flowing water bath. The temperature of the water bath was monitored throughout using a submerged thermometer. Once the experiment was completed, tadpoles were immediately blotted dry to remove excess water, weighed using an electric balance (to the nearest 0.1 g) to account for size differences among individuals in metabolic rate calculations, and returned to 100% oxygenated water. Although tadpole activity levels were very low whilst sealed in the respiration tubes, some short bursts of spontaneous activity were observed. Therefore, the metabolic rate estimates are considered as routine, which includes resting metabolic rate plus any extra energy expenditure due to spontaneous activity and stress . Percentage of oxygen consumed was calculated by subtracting the oxygen saturation of the control tube (i.e. the oxygen used by microbial activity) from the oxygen saturation of each respiration tube (i.e. the total oxygen consumed by both tadpole and bacterial activity). Percentage saturation was converted to ml l-1 using standard conversion tables based on water temperature during the experiment (water temperature varied between 19-22°C depending on the date of the experiment). Routine metabolic rate (RMR) was calculated for each individual as millilitres of oxygen consumed per gram weight per hour (ml O2 g-1 h-1).
Ten individuals per site (one individual per family*ten families per site; different individuals to those used in the RMR experiment) were moved to individual containers at Gosner stage 36–39 (hind leg toe differentiation) and deprived of food for 48 hours. Tadpoles were sealed within individual containers containing 80 ml of water and cooled to 4°C for 24 hours to cause inactivity. Containers were then gradually cooled, over a period of six hours, just until all the water in the container became completely frozen. Following this, tanks were gradually warmed (over a period of 14 hours) to 15°C and this temperature was maintained for one hour. All tadpoles were assessed for normal swimming behaviour at this point and the number of individuals that were still alive and exhibiting normal behaviour were recorded as the measure of freeze survival.
To evaluate whether RMR or freeze survival varied by altitude (fixed factor of interest; considered as a categorical variable of low or high) a generalised linear mixed model approach (GLMM) was used, as implemented in R v2.12.1  using the lme4 package . Mountain was included as a random factor in all models and family was nested within altitude in the RMR model, but multiple individuals from the same egg mass were not used in the freeze tolerance experiments thus family was not included in the freeze tolerance model. Weight was also included as a linear covariate in the freeze survival models, but was already accounted for in the measurement of RMR (ml O2 g-1 h-1). For RMR the model was assessed under a normal distribution and for freeze survival a binomial distribution was used. Each model parameter and interaction was sequentially removed from the models and a likelihood ratio test used to evaluate parameter significance. Only parameters that significantly changed the log likelihood when removed from the model were included in the final model. Significant differences due to the model parameters that showed interactions were evaluated using a Tukey’s HSD test , which works in conjunction with an ANOVA of the GLMM. A post-hoc power analysis of the ANOVA was carried out for RMR using G*Power v3.1.3 ; effect size was calculated using the formula in Krebs .
Animal ethics statement
All the protocols used in this study were approved by the U.K. Home Office: Project License 60/4110.
Availability of supporting data
The data sets supporting the results of this article are available in the Dryad repository, doi:10.5061/dryad.ks2j1, https://datadryad.org/resource/doi:10.5061/dryad.ks2jl .
Hangartner SB, Laurila A, Räsänen K: Adaptive divergence of the moor frog (Rana arvalis) along an acidification gradient. BMC Evol Biol. 2011, 11: 366-378. 10.1186/1471-2148-11-366.
Narum SR, Campbell NR, Kozfkay CC, Meyer KA: Adaptation of redband trout in desert and montane environments. Mol Ecol. 2010, 19: 4622-4637. 10.1111/j.1365-294X.2010.04839.x.
Voituron Y, Joly P, Eugène M, Barré H: Freezing tolerance of the European water frogs: the good, the bad, and the ugly. Am J Physiol Regul Integr Comp Physiol. 2005, 288: 1563-1570. 10.1152/ajpregu.00711.2004.
Hangartner SB, Laurila A, Räsänen K: Adaptive divergence in moor frog (Rana arvalis) populations along an acidification gradient: Inferences from QST-FST correlations. Evolution. 2012, 66: 867-881. 10.1111/j.1558-5646.2011.01472.x.
Lee C-R, Mitchell-Olds T: Quantifying effects of environmental and geographical factors on patterns of genetic differentiation. Mol Ecol. 2011, 20: 4631-4642. 10.1111/j.1365-294X.2011.05310.x.
Sears MW: Resting metabolic expenditure as a potential source of variation in growth rates of the sagebrush lizard. Comp Biochem Physiol Part A. 2005, 140: 171-177.
Laugen AT, Laurila A, Räsänen K, Merilä J: Latitudinal countergradient variation in the common frog (Rana temporaria) development rates – evidence for local adaptation. J Evol Biol. 2003, 16: 996-1005. 10.1046/j.1420-9101.2003.00560.x.
Bullock T: Compensation for temperature in the metabolism and activity of poikilotherms. Biol Rev. 1955, 30: 311-342.
Conover DO, Present TMC: Countergradient variation in growth rate: compensation for length of the growing season among Atlantic silversides from different latitudes. Oecologia. 1990, 83: 316-324.
Orizaola G, Quintela M, Laurila A: Climatic adaptation in an isolated and genetically impoverished amphibian population. Ecography. 2010, 33: 730-737. 10.1111/j.1600-0587.2009.06033.x.
Sinsch U: Thermal Influences on the habitat preference and the diurnal activity in three European Rana species. Oecologia. 1984, 64: 125-131. 10.1007/BF00377554.
Conover DO, Schultz ET: Phenotypic similarity and the evolutionary significance of countergradient variation. Trends Ecol Evol. 1995, 10: 248-252. 10.1016/S0169-5347(00)89081-3.
Lee W-S, Monaghan P, Metcalfe NB: The trade-off between growth rate and locomotor performance varies with perceived time until breeding. J Exp Biol. 2010, 213: 3289-3298. 10.1242/jeb.043083.
Räsänen K, Laurila A, Merilä J: Maternal investment in egg size: environment- and population-specific effects on offspring performance. Oecologia. 2005, 142: 546-553. 10.1007/s00442-004-1762-5.
Tanino KK, Storey KB: Temperature Adaptation in a Changing Climate: Nature at Risk. 2012, 256-Google ebook: http://books.google.com/books?hl=en&lr=&id=vzuDsMEmUKcC&pgis=1 CABI,
Lindgren B, Laurila A: Physiological variation along a geographical gradient: is growth rate correlated with routine metabolic rate in Rana temporaria tadpoles?. Biol J Linn Soc. 2009, 98: 217-224. 10.1111/j.1095-8312.2009.01255.x.
Block W, Somme L: Low temperature adaptations in beetles from the Sub-Antarctic Island of South Georgia. Polar Biol. 1983, 2: 109-114. 10.1007/BF00303176.
Lencioni V: Survival strategies of freshwater insects in cold environments. J Limnol. 2004, 63: 45-55.
Walsh PT, Downie JR, Monaghan P: Larval over-wintering: plasticity in the timing of life-history events in the common frog. J Zool. 2008, 276: 394-401. 10.1111/j.1469-7998.2008.00503.x.
Moore JA: Temperature tolerance and rates of development in the eggs of Amphibia. Ecology. 1939, 20: 459-478. 10.2307/1930439.
Körner C: The use of "altitude" in ecological research. Trends Ecol Evol. 2007, 22: 569-574. 10.1016/j.tree.2007.09.006.
Altwegg R, Reyer H-U: Patterns of natural selection on size at metamorphosis in water frogs. Evolution. 2003, 57: 872-882. 10.1111/j.0014-3820.2003.tb00298.x.
Lesbarrères D, Schmeller DS, Primmer CR, Merilä J: Genetic variability predicts common frog (Rana temporaria) size at metamorphosis in the wild. Heredity. 2007, 99: 41-46. 10.1038/sj.hdy.6800961.
Killen SS, Marras S, McKenzie DJ: Fuel, fasting, fear: routine metabolic rate and food deprivation exert synergistic effects on risk-taking in individual juvenile European sea bass. J Anim Ecol. 2011, 80: 1024-1033. 10.1111/j.1365-2656.2011.01844.x.
Reid D, Armstrong JD, Metcalfe NB: The performance advantage of a high resting metabolic rate in juvenile salmon is habitat dependent. J Anim Ecol. 2012, 81: 868-875. 10.1111/j.1365-2656.2012.01969.x.
Burton T, Killen SS, Armstrong JD, Metcalfe NB: What causes intraspecific variation in resting metabolic rate and what are its ecological consequences?. Proc R Soc B. 2011, 278: 3465-3473. 10.1098/rspb.2011.1778.
Arendt JD: Size-fecundity relationships, growth trajectories, and the temperature-size rule for ectotherms. Evolution. 2010, 65: 43-51.
Alvarez D, Metcalfe NB: The tradeoff between catch-up growth and escape speed: variation between habitats in the cost of compensation. Oikos. 2007, 116: 1144-1151. 10.1111/j.0030-1299.2007.15861.x.
Beck CW, Congdon JD: Effects of age and size at metamorphosis on performance and metabolic rates of Southern toad, Bufo terrestris, metamorphs. Funct Ecol. 2000, 14: 32-38. 10.1046/j.1365-2435.2000.00386.x.
Bayne B: Relations between variable rates of growth, metabolic costs and growth efficiencies in individual Sydney rock oysters (Saccostrea commercialis). J Exp Mar Bio Ecol. 2000, 251: 185-203. 10.1016/S0022-0981(00)00211-2.
Steyermark AC: A high standard metabolic rate constrains juvenile growth. Zoology. 2002, 105: 147-151. 10.1078/0944-2006-00055.
Kuan S-H, Lin YK: Bigger or faster? Spring and summer tadpole cohorts use different life-history strategies. J Zool. 2011, 285: 165-171. 10.1111/j.1469-7998.2011.00836.x.
Trivedi MR, Browne MK, Berry PM, Dawson TP, Morecroft MD: Projecting climate change impacts on mountain snow cover in central Scotland from historical patterns. Arctic Antarct Alp Res. 2007, 39: 488-499. 10.1657/1523-0430(06-006)[TRIVEDI]2.0.CO;2.
Costanzo JP, Lee RE, Lortz PH: Glucose concentration regulates freeze tolerance in the wood frog Rana sylvatica. J Exp Biol. 1993, 181: 245-255.
Storey KB, Storey JM: Natural freeze tolerance in ectothermic vertebrates. Annu Rev Physiol. 1992, 54: 619-637. 10.1146/annurev.ph.54.030192.003155.
Layne JR, Jones AL: Freeze tolerance in the gray treefrog: cryoprotectant mobilization and organ dehydration. J Exp Zool. 2001, 290: 1-5. 10.1002/jez.1030.
Beebee TJC: Amphibian breeding and climate. Nature. 1995, 374: 219-220. 10.1038/374219a0.
Duellman W, Trueb L: Biology of Amphibians. 1986, Baltimore: The John Hopkins University Press, 670-
Wells KD: The Ecology and Behavior of Amphibians. Volume 2007. 2007, 1400-Google eBook: http://books.google.com/books?id=eDKEKy5JJbIC&pgis=1 University of Chicago Press,
Muir AP, Thomas R, Biek R, Mable BK: Using genetic variation to infer associations with climate in the common frog, Rana temporaria. Mol Ecol. 2013, 22: 3737-3751. 10.1111/mec.12334.
Gasc JP, Maurin H, Oliviera ME, Sofianidou TS, Veith M, Zuiderwijk ACM, Cabela A, Crnobrnja-isailovic D, Dolmen K, Grossenbacher K, Haffner P, Lescure J, Martens H, Martinez-Rica JP: Atlas of Amphibians and Reptiles in Europe. 1997, Paris, France: Societas Europaea Herpetologica & Museum National d’Histoire, 496-
Vences M, Grossenbacher K, Puente M, Palanca A, Vieites DR: The Cambales fairy tale: elevational limits of Rana temporaria (Amphibia: Ranidae) and other European amphibians revisited. Folia Zool. 2003, 52: 189-202.
Inns H: Britain’s Reptiles and Amphibians. 2009, UK: WildGuides, 164-
Laurila A, Merilä J, Pakkasmaa S: Influence of seasonal time constraints on growth and development of common frog tadpoles: a photoperiod experiment. Oikos. 2001, 95: 451-460. 10.1034/j.1600-0706.2001.950310.x.
Muir AP, Biek R, Thomas R, Mable BK: Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria). Mol Ecol. 2014, 23: 561-574. 10.1111/mec.12624.
Thompson DBA, Brown A: Biodiversity in montane Britain: habitat variation, vegetation diversity and some objectives for conservation. Biodivers Conserv. 1992, 1: 179-208. 10.1007/BF00695915.
Briggs D, Smithson P, Addison K, Atkinson K: Fundamentals of the Physical Environment. 1997, London and New York: Routledge, 557-2
Crawley MJ: Statistics: An Introduction Using R. 2005, Chichester, UK: Wiley, 327-
Blaustein AR, Belden LK, Olson DH, Green DM, Root TL, Kiesecker JM: Amphibian breeding and climate change. Conserv Biol. 2001, 15: 1804-1809. 10.1046/j.1523-1739.2001.00307.x.
Phillimore AB, Hadfield JD, Jones OR, Smithers RJ: Differences in spawning date between populations of common frog reveal local adaptation. PNAS. 2010, 107: 8292-8297. 10.1073/pnas.0913792107.
Corn PS: Amphibian breeding and climate change: importance of snow in the mountains. Conserv Biol. 2003, 17: 622-625. 10.1046/j.1523-1739.2003.02111.x.
Urban MC: Microgeographic adaptations of spotted salamander morphological defenses in response to a predaceous salamander and beetle. Oikos. 2010, 119: 646-658. 10.1111/j.1600-0706.2009.17970.x.
Torres-Dowdall J, Handelsman CA, Reznick DN, Ghalambor CK: Local adaptation and the evolution of phenotypic plasticity in Trinidadian guppies (Poecilia reticulata). Evolution. 2012, 66: 3432-3443. 10.1111/j.1558-5646.2012.01694.x.
Dahl E, Backström T, Winberg S, Laurila A: Is growth hormone expression correlated with variation in growth rate along a latitudinal gradient in Rana temporaria?. J Zool. 2011, 285: 85-92. 10.1111/j.1469-7998.2011.00817.x.
Pasanen S, Karhapää M: Can boreal common frog (Rana temporaria L.) survive in frost?. Ann Zool Fennici. 1997, 34: 247-250.
Zhen Y, Ungerer MC: Clinal variation in freezing tolerance among natural accessions of Arabidopsis thaliana. New Phytol. 2008, 177: 419-427.
Somme L, Zachariassen KE: Adaptations to low temperature in high altitude insects from Mount Kenya. Ecol Entomol. 1981, 6: 199-204. 10.1111/j.1365-2311.1981.tb00606.x.
Claussen DL, Townsley M, Bausch R: Supercooling and freeze-tolerance in the European wall lizard, Podarcis muralis, with a revisional history of the discovery of freeze-tolerance in vertebrates. J Comp Physiol B. 1990, 160: 137-143. 10.1007/BF00300945.
Costanzo JP, Claussen DL: Natural freeze tolerance in the terrestrial turtle, Terrapene carolina. J Exp Zool. 1990, 254: 228-232. 10.1002/jez.1402540215.
Voituron Y, Barré H, Ramløv H, Douady CJ: Freeze tolerance evolution among anurans: frequency and timing of appearance. Cryobiology. 2009, 58: 241-247. 10.1016/j.cryobiol.2009.01.001.
Kernan M, Hughes M, Helliwell RC: Chemical variation and catchment characteristics in high altitude lochs in Scotland, U.K. Water Air Soil Pollut Focus. 2002, 2: 61-73.
Herreid CF, Kinney S: Temperature and development of the wood frog, Rana sylvatica, in Alaska. Ecology. 1967, 48: 579-590. 10.2307/1936502.
Sork VL, Davis FW, Westfall R, Flint A, Ikegami M, Wang H, Grivet D: Gene movement and genetic association with regional climate gradients in California valley oak (Quercus lobata Née) in the face of climate change. Mol Ecol. 2010, 19: 3806-3823. 10.1111/j.1365-294X.2010.04726.x.
Gosner KL: A simplified table for staging anuran embryos and larvae with notes on identification. Herpetologica. 1960, 16: 183-190.
Lewis WM: Morphological adaptations of cyprinodontoids for inhabiting oxygen deficient waters. Copeia. 1970, 2: 319-326.
R Development Core Team: R: A Language and Environment for Statistical Computing. 2013, Austria: R Found Stat Comput Vienna, http://www.R–project.org/,
Bates D, Maechler M, Bolker B: lme4: linear mixed-effects models using S4 classes. R Packag. 2011, Version 0.999375–42-
Tukey JW: Some selected quick and easy methods of statistical analysis. Trans N Y Acad Sci. 1953, 16: 88-97. 10.1111/j.2164-0947.1953.tb01326.x.
Erdfelder E, Faul F, Buchner A: GPOWER: a general power analysis program. Behav Res Methods Instrum Comput. 1996, 28: 1-11. 10.3758/BF03203630.
Krebs CJ: Ecological Methodology. 1989, New York: Harper and Row, 654-
Muir AP, Biek R, Mable BK: Data from: Behavioural and physiological adaptations to low-temperature environments in the common frog, Rana temporaria. BMC Evol Biol. 2014, https://datadryad.org/resource/doi:10.5061/dryad.ks2jl,
We thank Rose Hanley-Nickolls, Romaine Furmston-Evans, David Fettes and Martin Muir for assistance with fieldwork; and Aileen Adam and Elizabeth Kilbride for laboratory support. Thanks to the landowners that permitted access to sites. Fieldwork was supported by grants from the Royal Geographic Society, the Glasgow Natural History Society and the Scottish Mountaineering Trust. Permission for sampling from protected areas was granted by Scottish Natural Heritage. This study was supported by PhD CASE studentship funding from the Biotechnology and Biological Sciences Research Council, in partnership with the Royal Zoological Society of Scotland.
The authors declare that they have no competing interests.
This research forms part of APM’s PhD thesis work on the population genetics of R. temporaria in Scotland with BKM and RB; she was responsible for all aspects of the work, from experimental design, sampling, experimental procedures, analyses and writing. BKM contributed to experimental design, advice on analyses, and editorial content. RB contributed to editorial content. All authors read and approved the final manuscript.