Although general models of the processes involved in insect survival of freezing exist, there have been few studies directly investigating physiological processes during cooling, freezing and thawing, without which these models remain hypothetical. Here, we use open-flow respirometry to investigate the metabolism of the freeze-tolerant sub-Antarctic caterpillar Pringleophaga marioni Viette (Lepidoptera: Tineidae) during cooling,freezing and thawing and to compare animals exposed to non-lethal(–5.8°C) and lethal (–6.0°C, after which caterpillars are moribund for several days, and –18°C, after which caterpillars are completely unresponsive) freezing stress. We found a large decrease in metabolic rate (that is not associated with freezing) at–0.6±0.1°C and calculated a Q10 of 2.14×103 at this breakpoint. This breakpoint is coincident with the critical thermal minimum (CTmin) and is hypothesised to be a metabolic manifestation of the latter, possibly a failure of the Na+/K+-ATPase pump. This provides a plausible link between processes at the cellular level and observations of the action of the CTmin at tissue and whole-organism levels. Caterpillars froze at –4.6±0.1°C and had detectable metabolism when frozen. Post-thaw, metabolic rates were lower than pre-freezing measurements. Post-thaw metabolic rates did not differ between temperatures that did and did not kill the caterpillars, which suggests that mortality may be a result of a breakdown in processes at the organismal,rather than cellular, level of organisation.

Although the majority of insects are killed by freezing and must avoid internal ice formation to survive low temperatures, many insect species from a variety of habitats and taxa survive sub-zero temperatures by withstanding the formation of ice in their bodies and are termed `freeze tolerant'(Lee, 1991; Ramløv, 2000; Sinclair et al., 2003a). Although the biochemical correlates of this freezing survival are reasonably well known (Lee, 1991; Sømme, 1999; Ramløv, 2000; Sinclair et al., 2003b), the actual processes that occur during freezing and thawing have not been investigated to the same extent (but see Sinclair and Wharton, 1997; Joanisse and Storey, 1998; Neufeld and Leader, 1998; Kristiansen and Zachariassen,2001 for studies that do examine the freezing or thawing process explicitly). The general model of the response of freeze-tolerant insects to low temperatures is that as the temperature decreases, the insect will pass through the critical thermal minimum (CTmin –resulting in chill-coma) and eventually reach the temperature of crystallisation (Tc – the temperature at which the body fluids freeze). If the insect continues to cool after freezing, it will eventually reach a lower lethal temperature, after which it will not survive when thawed (Lee, 1991; Ramløv, 2000).

Survival of internal ice formation is generally attributed to a process of osmotic dehydration. Ice forming in the haemolymph creates concentrated fluid pockets that dehydrate cells and restrict ice to the extracellular spaces(Zachariassen, 1985; Sinclair and Wharton, 1997; Ramløv, 2000), although intracellular ice formation probably occurs in some species(Salt, 1962; Lee et al., 1993; Wharton and Ferns, 1995; Davis and Lee, 2001). During thawing, therefore, water must be redistributed and ion gradients re-established. In the woodfly Xylophaga cincta, re-establishing the ion balance takes several hours(Kristiansen and Zachariassen,2001). Nonetheless, there are currently no explicit data about the cause of lower lethal temperatures in freeze-tolerant insects. Although mechanical stress (for example, ice crystals piercing cell membranes or lethal recrystallisation of ice) has been proposed as a cause of mortality in frozen insects (Salt, 1961), the relationship between water-to-ice conversion and temperature means that the overall proportion of water converted to ice will increase continuously with decreasing temperatures (Zachariassen,1985), and most hypotheses of causes of mortality tend to hinge on an osmotic stress threshold being crossed. These hypotheses include the inability of individuals to re-establish ion gradients in the nervous system and transgression of a minimum cell volume resulting in lethal Ca2+concentrations (Meryman, 1971; Kristiansen and Zachariassen,2001).

Although the hypothesised causes of lower lethal temperature act at the sub-organismal level, it is reasonable to assume that the results are manifest at higher levels of organisation. However, aside from the progression of ice formation (see Ramløv,2000 for a summary), there is little information about processes at the organismal level during freezing and thawing of insects. Kozhantsikov(1938) suggested that cold-hardy insects are those that are able to maintain a detectable level of metabolism below 0°. This idea was challenged on methodological grounds by Scholander et al. (1953), who pointed out that Kozhantsikov heated his animals to above 0°C to perform measurements. Scholander et al.(1953) and Salt(1958) both showed that the rates of change of metabolic rates with temperature (i.e. Q10) in frozen and supercooled insects were extremely high at sub-zero temperatures and that the rate of change in metabolic rate increases logarithmically as the 0° threshold is crossed. In later work, Zachariassen et al.(1979) directly addressed the manifestation of freezing injuries by investigating the post-thaw metabolism of freeze-tolerant Eleodes blanchardi beetles exposed to lethal and sub-lethal temperatures. They concluded that metabolic rate was elevated after freezing. Moreover, whole-organism metabolic rate did not differ between those beetles that survived and those that were killed by freezing, suggesting that(at least initially) mortality does not occur at the cellular level. Most recently, Block et al. (1998)investigated the metabolic rates of freeze-tolerant larvae and adults of the perimylopid beetles Hydromedion sparsutum and Perimylops antarcticus on sub-Antarctic South Georgia. They found that, while metabolic rates of both species did not differ pre- and post-chilling(exposure to non-freezing sub-zero temperatures), metabolic rates of both larvae and adults of P. antarcticus were significantly higher after a brief freezing event but did not change after freezing in H. sparsutum. They concluded that the elevated metabolic rate of the former species is indicative of a rapid repair of slight freeze-induced injuries that explains its greater cold tolerance but did not investigate the effects of freezing injury or mortality within the frozen beetles.

All five of the preceding studies were conducted using closed-system respirometry, which integrates instantaneous rate values over time periods that are often prolonged. Therefore, this method generally does not reveal activity on the part of the study animal nor does it provide the fine resolution required to investigate short-term temporal variation in metabolic rate (Lighton, 1991; Addo-Bediako et al., 2002). By contrast, open-flow respirometry has the advantage that it enables the investigation of rapid, real-time metabolic responses to thermal and other stressors. To date, only Irwin and Lee(2002) have used open-flow respirometry to compare the metabolic rate of frozen and supercooled insects. They examined the larvae of the goldenrod gallfly Eurosta solidaginis, concluding that the metabolic rate of frozen larvae was lower than that of supercooled larvae, although they did not present metabolic data from during the freezing, thawing or post-thawing processes. Here, we use open-flow respirometry to explore the metabolic changes that occur during cooling, freezing and thawing, using the sub-Antarctic caterpillar Pringleophaga marioni as a model freeze-tolerant organism, and compare metabolic patterns in individuals exposed to lethal and non-lethal freezing stress.

Study animals

Pringleophaga marioni Viette (Lepidoptera: Tineidae) is a key decomposer on Marion Island (46°54′ S 37°45′ E; Crafford et al., 1986). The long-lived (several years) larvae occur in most of the habitats on the island,from sea level to ∼1000 m, where they are likely to experience multiple freeze–thaw events during their lifetime(Crafford, 1990; Crafford and Chown, 1992). Klok and Chown (1997) studied the environmental physiology of P. marioni larvae and determined that caterpillars have a critical thermal minimum (CTmin) of–0.6±0.2°C and that they are freeze tolerant, with a mean temperature of crystallisation (Tc) of–5.0±0.2°C and mortality at temperatures below approximately–6°C to –10°C. Recently, Sinclair and Chown(2003) have shown that P. marioni larvae do not exhibit a rapid cold-hardening response but that cold hardiness may be enhanced by prior desiccation or heat stress, suggesting that heat shock proteins may play some role in the freezing survival of this species.

Pringleophaga marioni larvae were collected from abandoned albatross nests in the vicinity of the Research Station on Marion Island during the April 2002 relief voyage. Caterpillars were placed in groups of five in 250 ml plastic jars filled with nest material and were returned under refrigeration to the laboratory in Stellenbosch within 10 days, where they were housed in the dark in a refrigerator (6.0±0.4°C) until they were used for experiments. Animals housed in this manner fed and moulted normally, and several individuals pupated.

Freezing and respirometry experiments

Caterpillars (0.202±0.016 g, mean mass ± s.e.m., N=29) were frozen individually in a 3 cm3 plastic respirometry cuvette. A 44-SWG type-T thermocouple was affixed to the bottom of the cuvette's interior. The cuvette was then placed inside a plastic bag in the bath of a Grant LTD-6(Grant Instruments, Cambridge, UK) cooling bath controlled by a Grant PZ-1 temperature programmer. Synthetic air (21% O2, balance N2) was scrubbed of CO2 with soda lime and of water with silica gel and Drierite. The air flowed at 75 ml min–1(Sidetrack mass flow controller) through a narrow-diameter copper coil(∼1.5 m total length) in the cooling bath before passing through the cuvette to a LiCOR Li6262 CO2/H2O infra-red gas analyser(IRGA). The respirometry system was controlled by DATACAN V software (Sable Systems, Henderson, NV, USA), which also provided automatic baselining and recorded the temperature inside the cuvette via a Sable Systems TC1000 thermocouple thermometer. All initial analyses of CO2, water and temperature traces were also performed with DATACAN V software. All measurements were corrected to standard temperature and pressure and expressed in ml CO2 h–1.

At the start of each experiment, the caterpillars were removed from the moist nest material, surface dried, weighed and quickly transferred to the cuvette and placed in the water bath at 1.6°C to equilibrate for one hour. Caterpillars were then cooled at 0.1 deg. min–1 to the test temperature, where they were held for 15 h (5 h for –18°C treatment,which still resulted in ∼10 h frozen and always resulted in mortality)before being rewarmed at 0.1 deg. min–1 to 1.6°C, where they were held for a further 2–5 h to enable the recovery from freezing to be recorded by the respirometer. Preliminary experiments allowed us to choose test temperatures that, under the experimental conditions, would be survived by all the caterpillars (–5.8°C, N=9), would result in mortality of all caterpillars with recovery to a `moribund' state(some uncoordinated movement; –6.0°C, N=8) and would kill all caterpillars outright without any apparent recovery of movement(–18°C, N=5). After freezing, the caterpillars were weighed again and removed to individual Petri dishes filled with nest material. Recovery (coordinated movement and a righting response) was monitored for seven days, after which all caterpillars were preserved in ethanol.

Data analyses

In the CO2 traces, a clear breakpoint in metabolism (b)was visible (Fig. 1), and the temperature at the start of this breakpoint was recorded. Freezing was evident in the temperature traces (owing to the sensitivity of the fine thermocouple),and the Tc – the temperature immediately before the exotherm – was recorded. Freezing was also clearly discernible in the CO2 (Fig. 1) and water traces. Where CO2 production was averaged for metabolic rate determination, temperature was averaged over the same period of recording. Temperatures of crystallisation and breakpoint distributions were not significantly different from normal, and means ± s.e.m. are presented. All caterpillars whose temperature dataset could yield breakpoint and Tc data were included in these analyses, including those from preliminary experiments to determine methods and test temperatures.

Fig. 1.

Representative data plots of temperature (grey, upper lines) and CO2 production (black, lower lines) from Pringleophaga marioni caterpillars exposed to freezing treatments at (A)–5.8°C (in which all caterpillars survive), (B) –6°C(which results in caterpillars being left moribund and then dead) and (C)–18°C (after which no caterpillars show muscle tone or response to stimuli). Breakpoint (b) and temperature of crystallisation(Tc) are marked. Small cycles in the CO2 trace shown in B are a result of minor temperature fluctuations at the IRGA. Whenever the range of these fluctuations exceeded the mean, the data were discarded. The break in the middle graph represents a change in data files (we were restricted by the memory that the acquisition software was able to allot to data gathering), and the slight difference in trace afterwards (equivalent to a maximum of 2.1 μW) is a result of re-baselining mid-data-acquisition. The masses of the caterpillars used in the presented data were 0.2342 g (A),0.1977 g (B) and 0.3826 g (C), accounting for the variation in the scale of the right-hand y-axes.

Fig. 1.

Representative data plots of temperature (grey, upper lines) and CO2 production (black, lower lines) from Pringleophaga marioni caterpillars exposed to freezing treatments at (A)–5.8°C (in which all caterpillars survive), (B) –6°C(which results in caterpillars being left moribund and then dead) and (C)–18°C (after which no caterpillars show muscle tone or response to stimuli). Breakpoint (b) and temperature of crystallisation(Tc) are marked. Small cycles in the CO2 trace shown in B are a result of minor temperature fluctuations at the IRGA. Whenever the range of these fluctuations exceeded the mean, the data were discarded. The break in the middle graph represents a change in data files (we were restricted by the memory that the acquisition software was able to allot to data gathering), and the slight difference in trace afterwards (equivalent to a maximum of 2.1 μW) is a result of re-baselining mid-data-acquisition. The masses of the caterpillars used in the presented data were 0.2342 g (A),0.1977 g (B) and 0.3826 g (C), accounting for the variation in the scale of the right-hand y-axes.

After baselining, mean CO2 production (ml CO2h–1) was calculated for a subset of data for each of four stages during the freeze–thaw cycle: (1) the equilibration period before freezing (24.6±1.2 min, N=16), (2) while the animal was frozen(382.3±43.2 min, N=4 per treatment), (3) during the immediate post-thaw period (the 34.5±2.9 min after cuvette temperature reached 1.6°, N=4 per treatment) and (4) >1 h post-thaw(105.4±20.2 min, N=4 per treatment)(Fig. 1). CO2production was converted to μW, assuming a respiratory quotient of 0.72(Withers, 1992). Minute variations in the temperature of the IRGA chamber meant that some metabolic rate measurements (while the caterpillars were frozen) were characterised by regular (∼30-min period) oscillations, in synchrony with small temperature changes in the room. After we had taken steps to better regulate the IRGA temperature, recordings where the range of metabolic rate (frozen) was greater than the mean metabolic rate (frozen) were not used for any metabolic rate comparisons, resulting in N=4 for each treatment for the majority of metabolic rate comparisons (although the data files were still useable for examination of b and Tc). Metabolic rate was log-transformed and compared between treatments and stages of the freeze–thaw cycle by a repeated-measures analysis of covariance (ANCOVA)using the general linear model (GLM) procedure on Statistica 6.1 (Statsoft Inc., Tulsa, OK, USA), using starting body mass of the caterpillar as the covariate. For ease of display, metabolic rate data are presented graphically as untransformed mass-specific data, although conclusions were drawn from the procedures outlined above.

Q10 was calculated from the slope of a regression of log10metabolic rate on temperature using the equation Q10=10(10×slope)(Cossins and Bowler, 1987). Q10 was calculated at five points during the freeze–thaw cycle: (1) during cooling before the breakpoint, (2) across the breakpoint,(3) during cooling after the caterpillar had frozen (–18°C treatment only), (4) during rewarming between –15°C and –10°C(–18°C treatment only) and (5) during thawing (–5°C to–2°C). Prior to cooling, metabolic rate was highly variable, and relationships where the slope was non-significant (F-test, P>0.05) were excluded from the analyses. Q10 was compared between treatments using Kruskal–Wallis Rank analysis of variance (ANOVA) and between cooling and thawing using Wilcoxon's paired test(Sokal and Rohlf, 1995).

Water loss

Water loss rate measured by the IRGA was baselined and then integrated for the 30-min period immediately prior to temperature decrease and for 240 min(or as long as the trace allowed if less than this) to give total mg H2O lost during these periods. For the post-thaw water loss, time was standardised to 200 min to allow comparison between treatments. Water loss was not normally distributed, despite transformation, and was compared between treatments using a generalized linear model (GLZ) in Statistica 6.1 with starting body mass as a covariate.

Mass loss during the course of the experiment was assumed to be entirely water (although several caterpillars did produce faecal pellets and/or vomit in the cuvette, the non-water mass of these was assumed to be negligible compared with the water lost). An ANCOVA with starting body mass as the covariate was used to compare total water lost between treatments.

Metabolic rate declined as the temperature decreased from 1.6°C to<0°C (Fig. 1). Almost all (16 of 17) animals showed a clear breakpoint in CO2 production where the slightly irregular CO2 trace (indicative of asynchronous spiracular control; Wasserthal,1996) decreased precipitously and disproportionately to the rate of temperature decline (point b in Fig. 1). The mean breakpoint temperature (± s.e.m.) was–0.6±0.1°C (N=16, range –1.2°C to 0.4°C), and the decline in metabolic rate occurred over a period of 14.4±1.1 min (N=16). The freezing event was recorded simultaneously by the thermocouple (the latent heat of crystallisation) and by the water and CO2 channels. The latter was presumably a result of expulsion of water vapour and gases from the tracheal system with the expansion of water in the haemocoel [assuming a tracheal CO2partial pressure (PCO2) of 5 kPa(Levy and Schneiderman, 1966),the volume of gas released is 9.5±1.2% (N=15) of the volume of the caterpillars, which is close to estimates of tracheal volume of insects,which usually vary between 3% and 10%(Bridges et al., 1980; Schmitz and Perry, 1999)]. The mean Tc was –4.6±0.1°C (N=29,range –5.7°C to –3.2°C). Upon rewarming, metabolic rate increased with temperature and remained steady after thawing in all three temperature treatments (Fig. 1). All animals in the –5.8°C treatment survived and were active and coordinated upon removal from the cuvette. Three of eight animals exposed to –6°C showed limited response to stimulus upon removal from the cuvette but died within 5 days (of the remainder, one survived and was excluded from analysis, while the rest never responded to stimulus). There was no consistent qualitative or quantitative difference in any of the variables measured between moribund and dead caterpillars exposed to –6°C, so the results were considered together. None of the–18°C animals showed any response and all completely lacked muscle tone when removed from the cuvette.

In all treatments, metabolic rates before freezing were significantly higher than after thawing and, when frozen, caterpillars had significantly lower metabolic rates (F3,21=18.364, P<0.0001; Fig. 2). Metabolic rates measured in the first 30 min post-thaw did not differ significantly from those measured longer after thawing and did not differ between treatments(F2,7=3.624, P=0.08; Fig. 2).

Fig. 2.

Mass-corrected metabolic rate of Pringleophaga marioni larvae at different stages of freezing and thawing. Treatments: grey bars,–5.8°C; stippled bars, –6°C; hatched bars,–18°C. Means ± s.e.m. are presented. Although mass-specific data are presented here, analyses were performed as a repeated-measures analysis of covariance, with mass as a covariate. Bars with the same letter are not significantly different from one another (Tukey's HSD, P<0.05).

Fig. 2.

Mass-corrected metabolic rate of Pringleophaga marioni larvae at different stages of freezing and thawing. Treatments: grey bars,–5.8°C; stippled bars, –6°C; hatched bars,–18°C. Means ± s.e.m. are presented. Although mass-specific data are presented here, analyses were performed as a repeated-measures analysis of covariance, with mass as a covariate. Bars with the same letter are not significantly different from one another (Tukey's HSD, P<0.05).

There was no significant difference (Kruskal–Wallis Rank ANOVA, P>0.1) between treatments in pre-freezing Q10(comparison only possible between –6.0°C and –18°C groups due to noisy data; median test, χ2=0.667, d.f.=2, P=0.72), so these results were pooled. Q10 before the breakpoint was highly variable (Table 1) but had a median of 2.20 (range 1.18–2.98, N=10)whereas across the metabolic breakpoint it ranged from 4×102to 1×105 (N=15). During cooling after the animal had frozen, Q10 was low, as it was when the caterpillars were being warmed between –15°C and –10°C. During thawing,Q10 was significantly higher in the –6.0°C group(Kruskal–Wallis H=7.05, d.f.=2, N=11, P=0.0294).

Table 1.

Q10 measured at different points during the freezing and thawing of Pringleophaga marioni larvae frozen to different temperatures (-5.8°C, -6.0°C and -18°C)

Point during freeze/thawTreatmentMedian Q10Upper quartileLower quartile
Before breakpoint All 2.20 2.79 1.60 
Across breakpoint All 2.14×103 4.24×103 1.26×103 
Cooling when frozen -18°C 1.49 1.69 1.33 
Warming -15°C to -10°C -18°C 1.83 2.18 1.70 
Thawing -5°C to -2°C -5.8°C 8.65 16.72 7.11 
 -6°C 48.20* 75.14 23.00 
 -18°C 7.37 8.16 5.87 
Point during freeze/thawTreatmentMedian Q10Upper quartileLower quartile
Before breakpoint All 2.20 2.79 1.60 
Across breakpoint All 2.14×103 4.24×103 1.26×103 
Cooling when frozen -18°C 1.49 1.69 1.33 
Warming -15°C to -10°C -18°C 1.83 2.18 1.70 
Thawing -5°C to -2°C -5.8°C 8.65 16.72 7.11 
 -6°C 48.20* 75.14 23.00 
 -18°C 7.37 8.16 5.87 

A significant difference (P<0.05) in post-thaw Q10between the -6.0°C treatment and the other treatments is indicated with an asterisk.

Flow-through gas analysis indicated that water loss rate did not differ between groups before freezing (Wald χ2=1.42, d.f.=2, P=0.49; Fig. 3A). The–18°C group lost considerably more body mass than the other two groups (Wald χ2=44.35, d.f.=2, P<0.0001; Fig. 3B), and the IRGA results suggest that this difference is attributable to the post-thaw period(F2,18=11.313, P<0.001; Fig. 3A).

Fig. 3.

Water lost by Pringleophaga marioni caterpillars after freezing,measured by (A) infra-red gas analyzer before freezing and after thawing during the same run (note that quantity of water lost is standardized to different time periods) and (B) by total mass lost during the run. Grey bars,–5.8°C; stippled bars, –6°C; hatched bars,–18°C. (A) IRGA data are medians ± minimum and maximum. Sample sizes: –5.8°C (N=4); –6°C (N=6);–18°C (N=4). (B) Mass loss data were analysed as an ANCOVA with starting body mass as a covariate, although mass-specific means ± s.e.m. are presented. Sample sizes:–5.8°C (N=9); –6°C (N=8);–18°C (N=5).

Fig. 3.

Water lost by Pringleophaga marioni caterpillars after freezing,measured by (A) infra-red gas analyzer before freezing and after thawing during the same run (note that quantity of water lost is standardized to different time periods) and (B) by total mass lost during the run. Grey bars,–5.8°C; stippled bars, –6°C; hatched bars,–18°C. (A) IRGA data are medians ± minimum and maximum. Sample sizes: –5.8°C (N=4); –6°C (N=6);–18°C (N=4). (B) Mass loss data were analysed as an ANCOVA with starting body mass as a covariate, although mass-specific means ± s.e.m. are presented. Sample sizes:–5.8°C (N=9); –6°C (N=8);–18°C (N=5).

The metabolic breakpoint and critical thermal minimum

There is a steep decline in metabolic rate as Pringleophaga marioni larvae are cooled (the breakpoint; Fig. 1), evidenced by an increase in Q10 of several orders of magnitude at the breakpoint(Table 1). A threshold in metabolism around 0°C that results in an exponential change in metabolic rate was reported by both Salt(1958) and Scholander et al.(1953), who interpreted this decline to be a consequence of freezing. One advantage of flow-through respirometry is that observations can be made in real time, so that Q10 can be calculated from the actual rate of change of metabolic rate over a temperature change. We have shown that there is a rapid change in metabolic rate at a threshold temperature that is clearly not associated with the freezing of the animal (Fig. 1). Rather, the start of this decline(–0.6±0.1°C) closely matches the CTminmeasured previously for this species (–0.6±0.2°C; Klok and Chown, 1997). This strongly suggests that the observed breakpoint is a metabolic manifestation of the CTmin.

Hosler et al. (2000) found a steady decrease in the resting potential of flight muscle neurons of both honey bees and Drosophila with decreasing temperature. They proposed that the chill-coma temperature was the temperature at which the Na+/K+-ATPase pump could no longer maintain nerve cell polarisation to a level where action potential could be produced. Since the contribution of transmembrane ion pumps to total metabolism is estimated to be>55% of basal metabolic rate at normal temperatures(Zachariassen, 1996; Hulbert and Else, 2000), a threshold temperature for pump activity is a plausible explanation for the drop in metabolism we observed.

Because of its ecological significance, the CTmin (or the closely comparable chill-coma temperature) has received attention in both field and laboratory studies of several species of insects(Chown, 2001; Gibert and Huey, 2001; Kelty and Lee, 2001; Hoffmann et al., 2003). For example, Gibert and Huey(2001) found that the chill-coma temperature of Drosophila differed between species and populations from different latitudes (suggestive of natural selection) and could be modified within those populations by developmental temperature(indicative of developmental plasticity). Klok and Chown(2003) found that the CTmin of sub-Antarctic weevils (Ectemnorhinus)responded to acclimation, with higher acclimation temperatures resulting in significantly higher CTmin values. Activity of the Na+/K+-ATPase pump may be modulated to reduce energy expenditure in response to heat shock and anoxia(Hochachka et al., 1996; Wu et al., 2002), and there is a variety of isomeric sodium pumps produced by insects under different environmental conditions, stages of development and tissues(Emery et al., 1998). This variation could be adequate both as a substrate for natural selection and possibly as a mechanistic explanation of acclimation responses in the CTmin.

An alternative mechanism explaining variation in CTminhas been proposed by Pörtner(2001, 2002), who suggested that the CTmin is a consequence of a mismatch between oxygen demand and delivery at low temperatures. However, because of the very large capacity for oxygen delivery in tracheated insects(Chapman, 1998), it seems unlikely that delivery failure sets limits to low temperature activity in these animals (Chapman, 1998;C. J. Klok, B. J. Sinclair and S. L. Chown, manuscript submitted for publication). Nevertheless, an oxygen-limited process (for example, ATP production by mitochondria), rather than a failure in the kinetics of the Na+/K+-ATPase pump at the membrane, would result in the same depolarisation, so the working hypothesis presented here is not inconsistent with Pörtner's(2001, 2002) assertions. Further investigation at the subcellular level is necessary to understand the processes leading to the CTmin in insects.

Metabolic processes during freezing and thawing

The metabolic rate of frozen caterpillars (89.2±24.9 μW g–1) was low but not zero(Fig. 2) and, during cooling after freezing, showed a Q10 comparable with that measured pre-freezing (pre-breakpoint) (Table 1). This suggests that the osmotically dehydrated cells are still functional and that, although the tracheal volume is presumably greatly reduced by the expansion of ice in the haemocoel, the tracheal system is still open and continues to provide an atmospheric interface for the cells. Assuming the animal's habitat allows access to atmospheric oxygen [see Scholander et al. (1953) for examples of chironomids frozen into lake ice; see Conradi-Larsen and Sømme(1973) for an example of terrestrial beetles locked in ice; see Lighton(1998) for a discussion of the hypoxic nature of below-ground habitats], this potentially allows the caterpillars to respire aerobically when frozen. An alternative explanation for this phenomenon is that the conversion of water into ice in freeze-tolerant insects is a slow process, continuing towards an asymptote over a period of hours (Lee and Lewis,1985; Ramløv and Westh,1993) and continuing to expand and expel CO2 from the tracheal system. While this is probably a partial explanation, the rate of CO2 does not attenuate to zero with increasing time frozen(Fig. 1), suggesting that there is a metabolic contribution to CO2 production.

The generally accepted model of the processes of freezing and thawing in freeze-tolerant insects requires that ice formation begins in the haemocoel and is confined to the extracellular spaces, resulting in osmotic dehydration of cells (Zachariassen, 1985;see Sinclair and Wharton(1997) for a demonstration of this). This osmotic dehydration, and subsequent passive movement of ions,means that water and ions must be redistributed and gradients re-established upon thawing (Kristiansen and Zachariassen, 2001). It stands to reason from this model that the freezing process is passive, but the apparent expulsion of CO2 from the tracheal system during freezing did not allow us to test this part of the model.

Redistribution of ions between intra- and extracellular compartments after thawing is probably an active process taking a few hours(Kristiansen and Zachariassen,2001). There is a decrease in metabolic rate in P. marioni in the four hours after thawing relative to before freezing(Fig. 2), suggesting that the initial redistribution of ions is not energetically expensive and that any necessary tissue repair does not begin immediately. Joanisse and Storey(1998) suggest that it would be advantageous for freeze-tolerant animals to depress post-freezing metabolism to reduce oxidative damage caused by reactive species accumulated while the animal is in the frozen state. The metabolic rate of Eleodes blanchardi beetles is elevated 20 h after thawing(Zachariassen et al., 1979),as is the metabolic rate of frozen Perimylops antarcticus beetles from South Georgia (Block et al.,1998), and it seems that this elevation may be associated with the repair of cellular or tissue damage sustained during the freeze–thaw process. Assuming that some cellular or tissue damage is sustained by P. marioni, metabolic rate might be expected to be elevated longer after thawing. Measurements of free radicals and antioxidants (along with longer-term monitoring of post-thaw metabolic rate) will be necessary to test this hypothesis.

Causes of low temperature mortality in frozen P. marioni

Low temperature mortality in frozen cells that would otherwise survive freezing has been hypothesised to be due to a threshold ice content resulting in a lethal minimum cell volume (Meryman,1971; Zachariassen et al.,1979; Zachariassen,1985). However, Zachariassen et al.(1979) did not find any metabolic differences between beetles that were exposed to lethal or sub-lethal freezing and concluded that mortality due to freezing might occur at the organismal, rather than the cellular, level as a result of neuronal membrane impermeability to sodium ions. Recent work by Yi and Lee(2003) on the gallfly Eurosta solidaginis examined the levels of injury sustained by different tissues at the lethal temperature of –80°C. Although nerve tissue was not tested, there is differential mortality between tissues, with gut tissue most tolerant to freezing and several tissues (intertegumentary muscle, the distal segment of Malpighian tubules and haemocytes) least tolerant, suggesting that damage to a specific tissue group such as nervous tissue could account for mortality associated with freezing. We compared metabolic rate and water loss in caterpillars that survived exposure to–5.8°C (which results in no mortality) to caterpillars exposed to–6.0°C (which results in moribund larvae that exhibit uncoordinated movement but die after a few days) and –18°C (after which larvae never show any response to stimuli). In common with Zachariassen et al.(1979), there was no difference in post-freezing metabolism between the three groups(Fig. 2), suggesting that, even in the dead caterpillars, most cells remained intact and capable of metabolism. The major difference between the three groups was that the–18°C group lost substantially more water than the other groups(Fig. 3), and this water loss is concentrated in the post-thaw period. This post-thaw water loss is likely to be a direct consequence of a lack of muscle tone due to nervous system disruption: thawed, dead caterpillars often vomit, and it is also likely that the water from this, as well as from the open anus and spiracles, provides an outlet for water loss from the caterpillar's body(Wharton, 1985). Immediate water loss does not, however, provide an explanation for mortality in the–6.0°C group. The notable feature of this group is their elevated Q10 during thawing (Table 1), which may suggest that the cellular processes associated with recovery from freezing have been disrupted.

Conclusions

In this paper, we have used open-flow respirometry to explicitly examine the metabolic processes associated with cooling, freezing and thawing in a freeze-tolerant insect. We have demonstrated a pronounced decline in metabolic rate that is not associated with freezing but is coincident with the CTmin. This finding provides a plausible link between the CTmin and mechanistic changes at the cellular and whole-organism levels. It also indicates ways in which natural selection might alter mechanisms at the cellular level, so giving rise to environmental variation in the CTmin. We have also shown that frozen caterpillars continue to respire when frozen, but we found no evidence that either freezing or thawing were active metabolic events. Moreover, there were no significant differences in metabolic rate between caterpillars that were and were not killed by freezing, suggesting that mortality from freezing occurs at an organismal, rather than cellular, level.

We thank Sue Jackson for insightful discussion, Annetjie Russell for statistical advice, Elrike Marais and two anonymous referees for helpful comments on earlier drafts of this manuscript, and Jan Crafford, Sarette Slabber and Erika Nortje for assistance collecting caterpillars. The Department of Environmental Affairs and Tourism provided logistic support during the 2002 Relief Voyage. This work was supported by the New Zealand Foundation for Research, Science and Technology, the United States Agency for International Development, and the National Research Foundation (GUN 2053804)and is a contribution to the SCAR RiSCC programme.

Addo-Bediako, A., Chown, S. L. and Gaston, K. J.(
2002
). Metabolic cold adaptation in insects: a large-scale perspective.
Funct. Ecol.
16
,
332
-338.
Block, W., Worland, M. R. and Bale, J. S.(
1998
). Respiratory responses to chilling and freezing in two sub-Antarctic insects.
Cryobiology
37
,
163
-166.
Bridges, C. R., Kestler, P. and Scheid, P.(
1980
). Tracheal volume in the pupa of the saturniid moth Hyalophora cecropia determined with inert gases.
Respir. Physiol.
40
,
281
-291.
Chapman, R. F. (
1998
).
The Insects. Structure and Function
. Cambridge: Cambridge University Press.
Chown, S. L. (
2001
). Physiological variation in insects: hierarchical levels and implications.
J. Insect Physiol.
47
,
649
-660.
Conradi-Larsen, E.-M. and Sømme, L.(
1973
). Anaerobiosis in the overwintering beetle Pelophila borealis.
Nature
245
,
388
-390.
Cossins, A. R. and Bowler, K. (
1987
).
Temperature Biology of Animals
. London: Chapman and Hall.
Crafford, J. E. (
1990
).
Patterns of energy flow in populations of the dominant insect consumers on Marion Island
. Ph.D. Thesis. Pretoria: University of Pretoria.
Crafford, J. E. and Chown, S. L. (
1992
). Microhabitat temperatures at Marion Island.
S. Afr. J. Antarct. Res.
22
,
51
-58.
Crafford, J. E., Scholtz, C. H. and Chown, S. L.(
1986
). The insects of sub-Antarctic Marion and Prince Edward Islands; with a bibliography of entomology of the Kerguelen Biogeographical Province.
S. Afr. J. Antarct. Res.
16
,
42
-84.
Davis, D. J. and Lee, R. E. (
2001
). Intracellular freezing, viability, and composition of fat body cells from freeze-intolerant larvae of Sarcophaga crassipalpis.
Arch. Insect Biochem. Physiol.
48
,
199
-205.
Emery, A. M., Billingsley, P. F., Ready, P. D. and Djamgoz, M. B. A. (
1998
). Insect Na+/K+-ATPase.
J. Insect Physiol.
44
,
197
-209.
Gibert, P. and Huey, R. B. (
2001
). Chill-coma temperature in Drosophila: effects of developmental temperature,latitude, and phylogeny.
Physiol. Biochem. Zool.
74
,
429
-434.
Hochachka, P. W., Buck, L. T., Doll, C. J. and Land, S. C.(
1996
). Unifying theory of hypoxia tolerance: Molecular metabolic defense and rescue mechanisms for surviving oxygen lack.
Proc. Natl. Acad. Sci. USA
93
,
9493
-9498.
Hoffmann, A. A., Sorensen, J. G. and Loeschcke, V.(
2003
). Adaptation of Drosophila to temperature extremes: bringing together quantitative and molecular approaches.
J. Therm. Biol.
28
,
175
-216.
Hosler, J. S., Burns, J. E. and Esch, H. E.(
2000
). Flight muscle resting potential and species-specific differences in chill-coma.
J. Insect Physiol.
46
,
621
-627.
Hulbert, A. J. and Else, P. L. (
2000
). Mechanisms underlying the cost of living in animals.
Annu. Rev. Physiol.
62
,
207
-235.
Irwin, J. T. and Lee, R. E. (
2002
). Energy and water conservation in frozen vs. supercooled larvae of the goldenrod gall fly, Eurosta solidaginis (Fitch) (Diptera: Tephritidae).
J. Exp. Zool.
292
,
345
-350.
Joanisse, D. R. and Storey, K. B. (
1998
). Oxidative stress and antioxidants in stress and recovery of cold-hardy insects.
Insect Biochem. Molec. Biol.
28
,
23
-30.
Kelty, J. D. and Lee, R. E., Jr (
2001
). Rapid cold-hardening of Drosophila melanogaster (Diptera: Drosophilidae)during ecologically based thermoperiodic cycles.
J. Exp. Biol.
204
,
1659
-1666.
Klok, C. J. and Chown, S. L. (
1997
). Critical thermal limits, temperature tolerance and water balance of a sub-Antarctic caterpillar, Pringleophaga marioni (Lepidoptera: Tineidae).
J. Insect Physiol.
43
,
685
-694.
Klok, C. J. and Chown, S. L. (
2003
). Resistance to temperature extremes in sub-Antarctic weevils: interspecific variation,population differentiation and acclimation.
Biol. J. Linn. Soc.
78
,
401
-414.
Kozhantsikov, I. W. (
1938
). Physiological conditions of cold-hardiness in insects.
Bull. Entomol. Res.
29
,
253
-262.
Kristiansen, E. and Zachariassen, K. E. (
2001
). Effect of freezing on the transmembrane distribution of ions in freeze-tolerant larvae of the wood fly Xylophagus cinctus (Diptera,Xylophagidae).
J. Insect Physiol.
47
,
585
-592.
Lee, R. E., Jr (
1991
). Principles of insect low temperature tolerance. In
Insects at Low Temperature
(ed. R. E. Lee, Jr and D. L. Denlinger), pp.
17
-46. New York: Chapman and Hall.
Lee, R. E., Jr and Lewis, E. A. (
1985
). Effect of temperature and duration of exposure on tissue ice formation in the gall fly, Eurosta solidaginis (Diptera, Tephritidae).
Cryo-Lett.
7
,
25
-34.
Lee, R. E., Jr, McGrath, J. J., Morason, R. T. and Taddeo, R. M. (
1993
). Survival of intracellular freezing, lipid coalescence and osmotic fragility in fat body cells of the freeze-tolerant gall fly Eurosta solidaginis.
J. Insect Physiol.
39
,
445
-450.
Levy, R. I. and Schneiderman, H. A. (
1966
). Discontinuous respiration in insects – II. The direct measurement and significance of changes in tracheal gas composition during the respiratory cycle of silkworm pupae.
J. Insect Physiol.
12
,
83
-104.
Lighton, J. R. B. (
1991
). Measurement on insects. In
Concise Encyclopedia on Biological and Biomedical Measurement Systems
(ed. P. A. Payne), pp.
201
-208. Oxford: Pergamon.
Lighton, J. R. B. (
1998
). Notes from underground: towards ultimate hypotheses of cyclic, discontinuous gas-exchange in tracheate arthropods.
Am. Zool.
38
,
483
-491.
Meryman, H. T. (
1971
). Osmotic Stress as a mechanism of freezing injury.
Cryobiology
8
,
489
-500.
Neufeld, D. S. and Leader, J. P. (
1998
). Cold inhibition of cell volume regulation during the freezing of insect Malpighian tubules.
J. Exp. Biol.
201
,
2195
-2204.
Pörtner, H. O. (
2001
). Climate change and temperature-dependent biogeography: oxygen limitation of thermal tolerance in animals.
Naturwissenschaften
88
,
137
-146.
Pörtner, H. O. (
2002
). Climate variations and the physiological basis of temperature dependent biogeography: systemic to molecular hierarchy of thermal tolerance in animals.
Comp. Biochem. Physiol. A
132
,
739
-761.
Ramløv, H. (
2000
). Aspects of natural cold tolerance in ectothermic animals.
Hum. Reprod.
15
, Suppl. 5,
25
-46.
Ramløv, H. and Westh, P. (
1993
). Ice formation in the freeze-tolerant alpine weta Hemideina maori Hutton(Orthoptera: Stenopelmatidae).
Cryo-Lett.
14
,
169
-176.
Salt, R. W. (
1958
). Relationship of respiration rate to temperature in a supercooled insect.
Can. J. Zool.
36
,
265
-268.
Salt, R. W. (
1961
). Principles of insect cold hardiness.
Annu. Rev. Entomol.
6
,
55
-74.
Salt, R. W. (
1962
). Intracellular freezing in insects.
Nature
193
,
1207
-1208.
Schmitz, A. and Perry, S. F. (
1999
). Stereological determination of tracheal volume and diffusing capacity of the tracheal walls in the stick insect Carausius morosus (Phasmatodea,Lonchodidae).
Physiol. Biochem. Zool.
72
,
205
-218.
Scholander, P. F., Flagg, W., Hock, R. J. and Irving, L.(
1953
). Studies on the physiology of frozen plants and animals in the Arctic.
J. Cell. Comp. Physiol.
42
, Suppl. 1,
1
-56.
Sinclair, B. J. and Chown, S. L. (
2003
). Rapid responses to high temperature and desiccation but not to low temperature in the freeze tolerant sub-Antarctic caterpillar Pringleophaga marioni(Lepidoptera, Tineidae).
J. Insect Physiol.
49
,
45
-52.
Sinclair, B. J. and Wharton, D. A. (
1997
). Avoidance of intracellular freezing by the New Zealand alpine weta Hemideina maori (Orthoptera: Stenopelmatidae).
J. Insect Physiol.
43
,
621
-625.
Sinclair, B. J., Addo-Bediako, A. and Chown, S. L.(
2003a
). Climatic variability and the evolution of insect freeze tolerance.
Biol. Rev.
78
,
181
-195.
Sinclair, B. J., Vernon, P., Klok, C. J. and Chown, S. L.(
2003b
). Insects at low temperatures: an ecological perspective.
Trends Ecol. Evol.
18
,
257
-262.
Sokal, R. and Rohlf, F. (
1995
).
Biometry: the Principles and Practice of Statistics in Biological Research
. New York: W. H. Freeman.
Sømme, L. (
1999
). The physiology of cold hardiness in terrestrial arthropods.
Eur. J. Entomol.
96
,
1
-10.
Wasserthal, L. T. (
1996
). Interaction of circulation and tracheal ventilation in holometabolous insects.
Adv. Insect Physiol.
26
,
297
-351.
Wharton, D. A. and Ferns, D. J. (
1995
). Survival of intracellular freezing by the Antarctic nematode Panagrolaimus davidi.
J. Exp. Biol.
198
,
1381
-1387.
Wharton, G. W. (
1985
). Water balance of insects. In
Comprehensive Insect Physiology, Biochemistry and Pharmacology
, vol.
4
(ed. G. A. Kerkut and L. I. Gilbert), pp.
565
-601. Oxford:Pergamon.
Withers, P. C. (
1992
).
Comparative Animal Physiology
. Orlando: Saunders.
Wu, B. S., Lee, J. K., Thompson, K. M., Walker, V. K., Moyes, C. D. and Robertson, R. M. (
2002
). Anoxia induces thermotolerance in the locust flight system.
J. Exp. Biol.
205
,
815
-827.
Yi, S. X. and Lee, R. E. (
2003
). Detecting freeze injury and seasonal cold-hardening of cells and tissues in the gall fly larvae, Eurosta solidaginis (Diptera: Tephritidae) using fluorescent vital dyes.
J. Insect Physiol.
49
,
999
-1004.
Zachariassen, K. E. (
1985
). Physiology of cold tolerance in insects.
Physiol. Rev.
65
,
799
-832.
Zachariassen, K. E. (
1996
). The water conserving physiological compromise of desert insects.
Eur. J. Entomol.
93
,
359
-367.
Zachariassen, K. E., Hammel, H. T. and Schmidek, W.(
1979
). Studies on freezing injuries in Eleodes blanchardi beetles.
Comp. Biochem. Physiol. A
63
,
199
-202.