Pacific hagfish, Eptatretus stoutii, can recover from 36 h of anoxia at 10°C. Such anoxia tolerance demands the mobilization of anaerobic fuels and the removal of metabolic wastes – processes that require a functional heart. The purpose of this study was to measure the metabolic response of the excised, cannulated hagfish heart to anoxia using direct calorimetry. These experiments were coupled with measurements of cardiac pH and metabolite concentrations, at multiple time points, to monitor acid–base balance and anaerobic ATP production. We also exposed hagfish to anoxia to compare the in vitro responses of the excised hearts with the in vivo responses. The calorimetry results revealed a significant reduction in the rate of metabolic heat production over the first hour of anoxia exposure, and a recovery over the subsequent 6 h. This response is likely attributable to a rapid anoxia-induced depression of aerobic ATP-production pathways followed by an upregulation of anaerobic ATP-production pathways such that the ATP production rate was restored to that measured in normoxia. Glycogen-depletion measurements suggest that metabolic processes were initially supported by glycolysis but that an alternative fuel source was used to support the sustained rates of ATP production. The maintenance of intracellular pH during anoxia indicates a remarkable ability of the myocytes to buffer/regulate protons and thus protect cardiac function. Altogether, these results illustrate that the low metabolic demand of the hagfish heart allows for near-routine levels of cardiac metabolism to be supported anaerobically. This is probably a significant contributor to the hagfish's exceptional anoxia tolerance.

Hagfish are benthic fishes that spend considerable periods of time buried in the sediment and the putrefying carcasses of aquatic animals where they feed (Lesser et al., 1997; Martini, 1998). These environments are routinely hypoxic and even anoxic, so it is unsurprising that hagfish display an exceptional tolerance to hypoxia and anoxia. Hansen and Sidell (1983) demonstrated that anoxia and cyanide poisoning did not have an effect on the in situ function of the heart from the Atlantic hagfish Myxine glutinosa. Similarly, the Pacific hagfish Eptatretus stoutii has been shown to tolerate 36 h of anoxia at 10°C with a regular heartbeat and full recovery (Cox et al., 2010,, 2011). Long-term anoxic survival like this demands the maintenance of certain organ functions within the animal, and one such critical organ is the heart. The maintenance of cardiac function under anoxic conditions allows for the mobilization of metabolic fuel stores (e.g. glycogen) from storage tissues (e.g. liver) and the removal of waste from metabolically active tissues (Stecyk et al., 2008). These two functions are critical to the anoxic survival of the animal and indeed, Cox et al. (2010) found that these are well maintained in anoxic hagfish: cardiac power output, a proxy for energy usage, decreased by only 25% in hagfish exposed to 36 h of anoxia. The maintenance of anoxic cardiac function in hagfish is thought to be associated with a relatively low cardiac metabolic requirement (Forster, 1991), such that the hagfish heart's routine power output – the lowest ever measured in a fish – is within the heart's maximum glycolytic potential (Farrell and Stecyk, 2007). The importance of glycolytic pathways to cardiac function in hagfish was demonstrated by Hansen and Sidell (1983). Using an in situ preparation on Atlantic hagfish, they found that despite inhibition of cardiac mitochondrial function with cyanide or sodium azide, cardiac output was maintained for at least 3 h. However, inhibiting glycolytic function with iodoacetate resulted in a significant decrease in cardiac output. The authors also found that the activity ratio of the enzymes pyruvate kinase (PK) to cytochrome c oxidase (CO) were 5.6-fold higher in the hagfish ventricle compared with that of cod. This was interpreted as the hagfish ventricle having a ‘more anaerobically geared’ metabolism (Hansen and Sidell, 1983). Together, these results suggest that anaerobic glycolysis plays an important role in maintaining cardiac function, and thus survival, in the anoxic hagfish.

Because the glycolytic capacity of the hagfish heart is thought to be capable of supporting routine cardiac function, and because this cardiac function is well maintained in anoxic environments, we hypothesized that the metabolic rate of the hagfish heart would be similarly maintained in anoxic environments. We tested this hypothesis by performing direct calorimetry on the excised, cannulated hearts of Pacific hagfish, E. stoutii, which are capable of beating for at least 24 h post-excision under anoxic conditions (Wilson, 2014). Direct calorimetry – the measurement of heat produced by, in this case, an excised heart – was used because it is the only technique that enables the continuous measurement of metabolic rate under anoxic conditions. In addition to the calorimetry measurements, we ran additional parallel anoxic exposures on excised, cannulated hearts and sampled them at three time points throughout the exposures to measure a variety of cardiac metabolites: lactate and glycogen levels were measured to estimate reliance on anaerobic ATP-production pathways, and intracellular pH (pHi) was measured to determine how the heart's ability to maintain acid–base balance was influenced by anoxic exposure.

In addition to the excised heart experiments, we exposed two groups of hagfish to 36 h of normoxia and anoxia. Hagfish heart tissue was analysed for the same metabolites described above and blood was analysed for a number of haematological parameters. This allowed us to compare the response of the in vitro heart with that of the in vivo heart, which was important because the in vitro heart would experience an almost immediate lack of O2 upon induction of environmental anoxia, whereas in vivo, the heart becomes progressively hypoxic until true anoxia is reached when the fish's large O2 stores are depleted; this could take a considerable period of time given the hagfish's low haemoglobin P50 value [∼7–12 mmHg PO2 (Wells et al., 1986; Forster et al., 1992)], large blood volume and low metabolic rate (Cox et al., 2010,, 2011). This fact made the in vivo exposures especially important to testing our hypothesis that metabolic activity is maintained in the hagfish heart during prolonged periods of anoxia.

List of symbols and abbreviations
     
  • Hct

    haematocrit

  •  
  • P50

    O2 tension at which hemoglobin is 50% saturated

  •  
  • PCO2

    partial pressure of dissolved CO2

  •  
  • PO2

    partial pressure of dissolved O2

  •  
  • pHi

    intracellular pH

  •  
  • SO2

    haemoglobin oxygen saturation

  •  
  • TO2

    whole blood total O2 content

Hagfish

Hagfish [Eptatretus stoutii (Lockington 1878)] were captured at ∼100 m depth in Barkley Sound, British Columbia, Canada and transported to the Department of Fisheries and Oceans-University of British Columbia (UBC) Centre for Aquaculture and Environmental Research, West Vancouver, British Columbia. They were housed in 1100 litre tanks with aerated, flow-through seawater (∼30‰ at 10±1°C) and fed squid once a week. Experimental animals were fasted for 1 week prior to use. All procedures were approved by the UBC Animal Care Committee (A08-0312) and conducted in accordance with their guidelines.

Calorimeter

The calorimeter used in our study was built in the lab of Lord Rothschild (Clarke, 1957). The instrument consists of a glass chamber embedded into either end of an aluminium cylinder (140 mm long ×70 mm diameter). One chamber is the experimental chamber the other is the reference chamber. When sealed with a rubber stopper, the volume of each glass chamber was 3.4 ml. Two stainless steel needles (23 gauge) were inserted through each stopper and connected to polyethylene (PE) 50 tubing so that each chamber had inflow and outflow lines. This allowed for physiological saline to be circulated through each chamber and for the cannulated heart to be perfused. A third needle was inserted through the rubber stopper of the experimental chamber for the insertion of a thermocouple so that temperature could be monitored. The PE tubing serving the calorimetry chambers was connected to a single peristaltic pump (Gilson Minipuls 3) that circulated saline from an Erlenmeyer flask through the system. Flow rate was set at 9.4 ml h−1 and flow was identical between chambers. An air-stone in the Erlenmeyer flask was used to aerate the saline or to remove oxygen using compressed nitrogen. An oxygen probe (Oakton DO 6+) placed in the Erlenmeyer flask was used to monitor dissolved oxygen concentrations throughout the experiment. The calorimeter was placed within a styrofoam cooler that was held within an insulated ice chest (Coleman 6-Day Xtreme, Golden, CO, USA). Inside the cooler, the calorimeter was placed on top of a layer of aluminium blocks that were used to draw heat from the experimental and reference chambers. The entire system was held within a temperature-controlled (8±1°C) environmental chamber for the duration of the study.

Heat calibration, completed using the procedure described by Regan et al. (2013), showed that 1 mV was equal to 0.67 mW. The calibration also revealed that the calorimeter achieved thermal equilibrium after approximately 2 h. This time lag was factored out of the heat traces using an ‘instantaneous correction’ (Bartholomew et al., 1981) for flow-through respirometry corrected for calorimetric heat data. The equation used was:
(1)

where Qeq is the time-corrected data point (in mW), QeT is the uncorrected data point at time T, QeT−1 is the uncorrected data point at time T−1, dT is the time interval between successive data points (in min) and Tlag is the calorimeter's time lag (in min).

Heart dissection for calorimetry

The animals were killed using a sharp blow to the head followed by decapitation. The systemic heart was quickly excised and then washed in ice-cold hagfish saline (in mmol l−1: 450 NaCl, 9 KCl, 10 CaCl2, 50 MgCl2 and 15 HEPES-Na buffer, pH 7.8). The bulbus arteriosus was removed so that saline could flow easily from the ventricle. A 3 cm cannula made from PE50 tubing was inserted through the atrium into the ventricle via the atrial ventricular valve. Surgical silk secured the cannula in place and the cannula was attached to the steel needle inserted through the rubber stopper of the experimental chamber. This allowed the heart to be perfused when sealed inside the chamber. Observation of the heart, prior to insertion into the system, revealed that the ventricle would inflate and then contract, releasing saline in a continuous pulsatile manner at a rate of approximately 6 bpm.

Measuring metabolic heat production

Once the heart was mounted inside the experimental chamber, aerated saline at 8°C was pumped through the system. This saline did not contain glucose because we wanted to accurately monitor the use of glycogen stores by the heart. To determine if we could measure the heat of metabolism with the calorimeter, we ran a number of preliminary experiments using conditions that would influence the metabolic rate of the heart. These included alternating periods of anoxia and normoxia, as well as the application of the poison iodoacetate to inhibit anaerobic glycolysis during anoxia exposure. In preliminary experiments we observed that if iodoacetate was not applied, the excised hagfish heart continued to contract for more than 24 h of experimentation, including prolonged periods of anoxia (Fig. S1).

To quantify the change in metabolic heat production with anoxia exposure, the hearts were perfused with aerated saline until a stable baseline was achieved (4.1±0.42 h). At this point, anoxia was induced by bubbling nitrogen into the circulating saline for the following 16 h and then iodoacetate was applied to the anoxic saline. An example of the data from such a protocol is Fig. 1A. The baseline voltage value varied slightly from run to run. It was therefore measured for each run once the output voltage stabilized after the application of iodoacetate. This baseline value was then subtracted from all measurements in the experiment. The dry mass of each heart was also obtained.

Converting metabolic heat production to ATP equivalents

To convert joules to approximate ATP equivalents in the normoxic hearts, we assumed the hearts were oxidizing glucose and yielded 36 moles of ATP per mole of substrate. Next, we calculated the theoretical quantity of waste heat in the entire aerobic ATP generation pathway by summing the standard free energy change for aerobic ATP production (7.3 kcal mol−1 ATP×36 moles ATP=+263 kcal) with the standard free energy change of glucose oxidation to CO2 and H2O (−686 kcal mol−1 glucose). This yielded a total of 423 kcal of heat produced per mole of glucose oxidized or 11.8 kcal of waste heat per mole of ATP produced via aerobic metabolism (423 kcal ×36 mol−1 ATP=11.8 kcal mol−1 ATP) equating to 49.2 kJ mol−1 ATP produced. We then used this number to convert the kJ detected by our calorimetry experiments to ATP equivalents produced by the normoxic hearts.

To convert joules to approximate ATP equivalents in the anoxic hearts, we assumed an entire reliance on anaerobic glycolysis and a range of glycolytic ATP yields, from 2 to 3 moles ATP per mole of substrate. Using the method described above and assuming a glucose-to-lactate free energy change of −47 kcal mol−1 glucose, we calculated 67.8 kJ mol−1 ATP (2 moles ATP per mole substrate) and 35.2 kJ mol−1 ATP (3 moles ATP per mole substrate) produced via anaerobic metabolism, which we used to convert the kJ detected by our calorimetry experiments to ATP equivalents produced by the anoxic hearts.

Heart preparation for tissue metabolite analysis

We used the same exposure protocols as described for the calorimetry studies to monitor changes in tissue pH and metabolite concentrations in the excised hearts at multiple time points during anoxia exposure. The only difference was that the glass test tubes containing the cannulated hearts were held in a beaker of saline instead of the calorimeter. The hearts were sampled after 6 h of normoxia, and after 1.5 h and 8 h of anoxia: 6 h of normoxia was chosen because this was the time required for the excised hearts to reach an initial thermal equilibrium in normoxia; 1.5 h of anoxia was chosen because the maximum decrease in heat production was seen at 2.38±0.83 h of anoxia and we wanted to identify the basis for suppression of aerobic metabolism; and 8 h of anoxia was selected as the third sampling time because the rate of metabolic heat production was still increasing at this point. It is important to note that all hearts were observed to be beating when sampled. To preserve the tissues, hearts were quickly removed from the system, rinsed in ice-cold saline, blotted dry and then quickly frozen on a ∼2 kg piece of aluminium stored at −80°C. These samples were then kept at −80°C until analysis.

Metabolite analyses

We measured cardiac concentrations of lactate, glycogen and glucose by first sonicating ∼40 mg of frozen ground cardiac tissue in 10 volumes of a metabolic inhibitor solution consisting of 150 mmol l−1 KF and 6 mmol l−1 nitrilotriacetic acid (Portner, 1990). We then measured the pHi of this homogenate using an Accumet microcombination pH electrode (Fisher Scientific, 13-620-96). Next, we prepared the remaining homogenate for metabolite measurements by acidifying it in 30% perchloric acid, neutralizing it in 3 mol l−1 potassium carbonate, then centrifuging it at 10,000 g for 10 min (4°C). We extracted the resulting supernatant and stored it at −80°C for future analyses. Once thawed, the supernatant was centrifuged at 10,000 g for 10 min (4°C) before measuring lactate, glycogen and glucose levels according to Bergmeyer et al. (1983).

In vivo anoxia exposure

To examine the influence of anoxia exposure on the blood composition and heart metabolites of hagfish, live animals were exposed to either 36 h of normoxia or anoxia at 8°C, as previously described (Cox et al., 2011). Following treatment or control exposure, samples of mixed venous blood were taken from the large caudal blood sinus, as previously described (Cox et al., 2011) and placed on ice. The fish were then killed and their hearts were quickly frozen. The blood was analysed immediately after sampling for haemoglobin (Hb) concentration, haematocrit (Hct), partial pressure of dissolved CO2 (PCO2), Hb O2 saturation (SO2) and whole blood pH. Plasma was also frozen and later analysed for glucose and lactate. The hearts were stored at −80°C until they were used to determine tissue metabolites and pHi.

Measurement of plasma composition

Hct was measured in triplicate using micro Hct tubes (10 μl) after centrifuging at 17,000 g for 3 min. Hb concentration was measured in triplicate with a Shimadzu UV-1800 spectrophotometer (Kyoto, Japan) using the cyanomethaemoglobin method and was calculated based on absorption measurements at 540 nm and using an extinction coefficient of 11. Whole-blood pH measurements were performed using a Radiometer BMS 3 Mk2 system, maintained at 5°C with a Radiometer acid–base analyzer PHM71 (Copenhagen, Denmark). Whole-blood total O2 content (TO2) was measured according to Tucker (1967). SO2 was calculated from TO2 after subtracting physically dissolved O2 according to Boutilier et al. (1984). Total CO2 in the blood plasma was measured using a Corning 965 CO2 analyzer (Corning, NY, USA), after centrifuging whole blood for 5 min at 5000 g. Plasma HCO3 concentration and blood PCO2 were calculated from the plasma pH and TCO2 using the Henderson–Hasselbalch equation and the physicochemical parameters reported by Boutilier et al. (1984). The concentrations of plasma lactate and glucose were measured using a YSI 2300 Stat Plus analyzer (Yellow Springs, OH, USA).

Statistical analysis

The differences in blood parameters, plasma components, metabolites in the hearts and pHi between the control and anoxia-exposed animals were tested using an unpaired, two tailed t-test (P<0.05). The differences in tissue metabolites, pHi, heat loss and ATP turnover rates from the calorimetry studies were tested using a single-factor ANOVA followed by a Bonferroni post hoc test (P<0.05).

Influence of anoxia exposure on the metabolic rate of the excised heart

Fig. 1A is a representative trace from one of the calorimetry experiments. The measured voltage is directly proportional to the heat generated by the heart and detected by the thermopiles surrounding the experimental chamber. This figure demonstrates that perfusion of the heart with anoxic saline causes an initial decrease in metabolic heat production, but that this recovers with time. This trace also demonstrates a decrease in metabolic heat production caused by the application of iodoacetate during anoxia exposure, supporting the idea that the heat being detected by the calorimeter was indeed metabolic. The inset in Fig. 1A shows the voltage before and after the heart was perfused with anoxic saline, and demonstrates that the heart had reached a steady heat output before it was made anoxic. The time required for the excised hearts to reach an aerobic steady state prior to the anoxia exposure was 4.1±0.42 h. The red line on Fig. 1A is the voltage output from a control experiment where anoxic saline was perfused through the chamber in the absence of a heart. The rapid deoxygenation of the chamber has an observable and consistent effect on the calorimeter's thermal equilibrium, probably because of a slight decrease in the temperature of the circulating saline caused by the compressed N2 bubbling in the Erlenmeyer flask. All voltage traces were corrected to compensate for this experimental artefact. The black line on Fig. 1B is the average metabolic heat produced by the hearts during anoxia exposure calculated directly from the voltage output and corrected for the initial decrease in output voltage caused by the application of anoxic saline to the system. The grey line in Fig. 1B is that generated when the average metabolic heat production data are de-smeared to compensate for the time lag of the calorimeter. The average rates of ATP and metabolic heat production, at discrete time points in Fig. 1C, result from this correction.

The metabolic heat production of the excised hearts at the initial aerobic steady state condition was 0.14±0.03 mW g−1. This equalled an estimated ATP turnover rate of 10.3±2.4 µmol g−1 h−1 (Fig. 1C). The de-smeared data shown in Fig. 1B illustrate that the application of anoxic saline causes metabolic heat production to rapidly drop to zero. After approximately 1 h, metabolic processes began to recover, so that by 1.5 h of anoxia, metabolic heat production was 0.168±0.05 mW g−1, which, depending on ATP yield per mole substrate, was equal to an estimated ATP turnover rate of 8.9±2.5 µmol g−1 h−1 (2 moles ATP per mole substrate) to 17.2±4.9 µmol g−1 h−1 (3 moles ATP per mole substrate) (Fig. 1C). Metabolic heat production continued to recover during the anoxia exposure so that by 8 h it was 0.18±0.05 mW g−1, equal to an estimated ATP turnover rate of 9.9±2.8 µmol g−1 h−1 (2 moles ATP per mole substrate) to 19.1±5.3 µmol g−1 h−1 (3 moles ATP per mole substrate) (Fig. 1C). At 16 h of anoxia exposure, just prior to the application of iodoacetate, heat production was 0.19±0.06 mW g−1, which is equal to an estimated ATP turnover rate of 10.3±3.04 µmol g−1 h−1 (2 moles ATP per mole substrate) to 19.9±5.9 µmol g−1 h−1 (3 moles ATP per mole substrate) (Fig. 1C). It is important to note that these heat-derived ATP-turnover calculations are estimations based on standard free energy changes and assumptions provided in the Materials and methods. Cellular concentrations of reactants and products will affect these estimations of ATP turnover.

Influence of anoxia on pHi and metabolic substrates

As the saline used to perfuse the excised hearts did not contain glucose, it is not surprising that tissue glycogen decreased during the anoxia exposure. This also confirms that there was metabolic activity. Between 1.5 h and 8 h of anoxia, glycogen stores decreased by 78% (Fig. 2A). These levels are lower than that measured after 6 h of perfusion with aerated saline (i.e. normoxia on Fig. 2A). These differences could not be tested, as we were only able to measure glycogen in a single heart that had been perfused with aerated saline. However, the glycogen level measured in this heart was similar to that measured in the hearts of the control animals in the in vivo study (Fig. 2A). The glycogen levels in the hearts sampled after 8 h of anoxia in the in vitro study were 87% lower than those in the hearts sampled prior to anoxia exposure, whereas the levels of glycogen in the hearts of the fish exposed to anoxia for 36 h in the in vivo study were 97% lower than those in the control fish (Fig. 2A).

Lactate concentrations in the excised hearts were 3- and 5-fold higher at 1.5 h and 8 h anoxia than they were at 0 h anoxia (Fig. 2B). However, this level of lactate accumulation is small compared with that of the in vivo hearts and that expected from the measured rates of glycogen depletion (Fig. 2A), and is probably a function of lactate loss to the perfusate. In the hearts from the in vivo study, the concentration of lactate was 0 in the normoxic group but 17.9±1.3 μmol g−1 in the hearts of the anoxic group (Fig. 2B).

The excised hearts from the in vitro study had pHi values of 7.03, 6.95 and 6.96 at 0 h, 1.5 h and 8 h of anoxia, respectively (Fig. 2C). While statistically significant, an absolute change in pH of 0.08 pH units over 8 h is relatively small for a completely anoxic tissue and suggests the heart was buffering and/or extruding the protons resulting from anaerobic glycolysis. Similarly, the hearts of the anoxia-exposed hagfish from the in vivo study were 0.07 pH units lower compared with those of the normoxia-exposed hagfish (Fig. 2C).

Influence of anoxia on behaviour and blood composition of hagfish

Normoxic hagfish remained in a coiled position within the exposure chamber and periodically changed position. Under anoxic conditions, the hagfish were uncoiled and displayed no body movement, consistent with what has been seen previously (Cox et al., 2010,, 2011).

Thirty-six hours of anoxia significantly altered blood composition (Table 1). For example, haemoglobin content was 1.9-times higher in anoxic hagfish than in normoxic hagfish (Table 1), probably largely due to the 1.7-times higher haematocrit of anoxic hagfish (Table 1). Haemoglobin oxygen saturation was 0% in the anoxic hagfish and 68.6% in the normoxic hagfish, suggesting that 36 h of anoxia was enough to exhaust the blood's oxygen stores (Table 1). PCO2 was 7-times higher in the anoxic fish and plasma lactate increased from ∼0 mmol l−1 in the normoxic hagfish to 28.0 mmol l−1 in the anoxic hagfish. Also, blood pH was reduced from 7.81 in the normoxic fish to 6.97 in the anoxic fish. Finally, the 3.3-times higher plasma glucose concentrations of anoxic fish suggest they were mobilizing glucose stores (Table 1).

The results of this study illustrate the significant ability of the Pacific hagfish heart to maintain metabolic function during prolonged periods of anoxia both in vitro and in vivo. Excised hearts initially responded to anoxia with a significant reduction in metabolic heat production, but recovered after 1.5 h to levels that translated to ATP production rates similar to those of normoxic control hearts. These metabolic rates confirm the hypothesis that the low metabolic demand of the hagfish heart allows routine levels of cardiac metabolism to be supported anaerobically.

Influence of anoxia exposure on the metabolic rate of excised hearts

The general description of cardiac function in the hagfish is ‘low and slow’. This is supported by the comparatively low activity of enzymes associated with both aerobic and anaerobic metabolism (Hansen and Sidell, 1983), ventricular blood pressures that are 5- to 10-times lower in hagfish than in any other fish (Cox et al., 2010; Forster et al., 1988; Johnsson and Axelsson, 1996) and a resting heart rate of E. stoutii in normoxia of approximately 8 bpm (Cox et al., 2010). In trout, heart rate is 31 bpm at 15°C (Altimiras and Larsen, 2000) and 45 bpm in bluefin tuna under routine conditions (Clark et al., 2013). This low level of cardiac work in the hagfish – the lowest of any fish species studied to date – led Forster (1991) to predict that it could be supported solely by the tissue's glycolytic potential. Such a strategy, fuelled by glycogen stores, would enable the maintenance of cardiac function during anoxia. The results of our study support this idea.

Integration of the heart glycogen measurements into the calorimetry results provides some interesting insight into fuel utilization by the excised hearts during anoxia exposure. First, it appears that the rate of glycogen utilization decreased during anoxia exposure. This is suggested by comparing the rate of glycogen utilization between sampling points. If we assume that the average amount of glycogen in the excised hearts at time 0 was similar to that measured in the normoxic hearts from the in vivo experiment (31.1±4.35 μmol g−1), the rate of glycolysis between time 0 and 1.5 h of anoxia was approximately 7 µmol glucosyl units g−1 h−1 whereas that between 1.5 h and 8 h of anoxia was 2.5 µmol glucosyl units g−1 h−1. These rates of glycogen utilization translate into ATP production rates of 14.0–21.0 (2 and 3 moles of ATP, respectively) and 4.9–7.4 µmol g−1 h−1, respectively (Fig. 3). Interestingly then, these calculated rates also suggest that as the period of anoxia continues, the rate of anaerobic metabolism estimated by calorimetry cannot be supported solely by the glycogen stored in the heart. The metabolic heat measurements at times 1.5 h, 8 h and 16 h of anoxia translate into ATP production rates of 8.9±2.5 (to 17.2±4.9), 9.9±2.8 (to 19.1±5.3) and 10.3±3.0 (to 19.9±5.9) µmol g−1 h−1, respectively. In addition, based on the assumption that the glycogen in the heart at time 0 was 31.1±4.35 μmol g−1, the amount measured in the hearts at 8 h indicates an 85% decrease in glycogen content. If we then assume that the glycogen stores were completely depleted over the next 8 h (8–16 h anoxia), the rate of ATP production during this time would have been approximately 1.2 µmol g−1 h−1. However, the sustained levels of metabolic heat over the 8–16 h time frame suggest that no such metabolic depression occurs (Fig. 3). We therefore think that another fuel, in addition to glycogen, may have been supporting metabolic activity during this time frame. Cox et al. (2010,, 2011) have previously suggested that the hagfish heart may use an alternative fuel source during anoxia to supplement glycogen stores, which must be the case given that the metabolic rate is maintained during anoxia in the current study (Fig. 3). One alternative fuel may be glycerol, which could be sourced from the breakdown of triglycerides stored in the heart and enter glycolysis via the enzyme glycerol kinase. We are currently examining this possibility.

The red-eared slider turtle Trachemys scripta and the crucian carp Carassius carassius are highly tolerant of anoxia but employ different survival strategies (Galli et al., 2013; Hicks and Farrell, 2000; Lau and Richards, 2011; Stecyk et al., 2008; van Waversveld et al., 1989). In the turtle, metabolic rate and cardiac output decrease rapidly with anoxia exposure and then remain low for the duration of the anoxia exposure (Galli et al., 2013; Hicks and Farrell, 2000; Lau and Richards, 2011; Stecyk et al., 2008). In the crucian carp, anoxia exposure causes a rapid decrease in cardiac output (metabolic rate has not been measured) that returns to pre-anoxia levels with the upregulation of anaerobic glycolysis, allowing the animal to remain active for at least 5 days (Stecyk et al., 2004). The anoxic response of the hagfish heart in this study is therefore similar to the anoxia response of the crucian carp, where metabolic rate is initially depressed and then recovers as a result of the upregulation of anaerobic ATP production pathways. Our results suggest that this metabolic recovery is quicker in the hagfish than in the crucian carp.

Cox et al. (2010) demonstrate that exposure of E. stoutii to anoxia causes heart rate to decrease by 55%, but that stroke volume nearly doubles. The net result is that cardiac output only decreases by ∼33% with anoxia exposure. This work also demonstrated that power output from the heart decreased by only 25%. As power output can be used as a surrogate for energy requirement, this demonstrates that there are still relatively high energy demands of the heart during anoxia exposure.

There are a number of factors that must be taken into account when considering the results of the calorimetry study. First, it was not possible to monitor the rate of contraction, but we did see the hearts contracting throughout the metabolite study under identical conditions and previous work has found that the excised heart of Pacific hagfish continues to beat during 24 h of anoxia exposure (Wilson, 2014). We did not set pre- or after-load on the excised heart because of the technical limitations associated with the calorimetry system. However, in vivo, the hagfish heart works at extremely low pressures with ventral aortic pressure being ∼1 kPa during anoxia exposure (Cox et al., 2010). Therefore, it would not take much resistance within the tubing to generate a physiological afterload for the cardiac contractions. The fact that we did not control the afterload on the heart may not be that significant an issue. The energy requirements of the hearts in our excised experiments, however, may not have been as high as that found in vivo. This may explain why Cox et al. (2010) observed a mass-specific cardiac power output (the product of cardiac output and ventral aortic blood pressure, expressed as mW g−1) that was approximately 50% higher than the metabolic rates (expressed in mW g−1) measured here in both normoxia and 12 h anoxia. Interestingly, while the excised hearts in our study showed a significant reduction in metabolic rate over the first hour of anoxia, the in vivo hearts measured in the Cox et al. (2010) study showed a significant increase in cardiac power output. This increased power output was primarily the result of an increase in ventral aortic blood pressure. Since the excised hearts in our system were held at a constantly low pressure, this may in part explain the different findings of the two studies.

The in vivo heart becomes progressively rather than abruptly anoxic. For example, Forster et al. (2001) estimated that the large blood volume of the hagfish has an O2 storage capacity of 137.9 ml O2 kg−1 for E. stoutii, a number twice that of any other fish. This represents a significant reserve of on-board oxygen that has been estimated to support the hagfish's routine metabolic rate [assuming an O2 consumption rate of 20.6±1.0 mg O2 kg−1 h−1 (Cox et al., 2010)] in complete anoxia for approximately 25 min; for a hagfish with a reduced aerobic metabolism, this would be even longer, with estimations between 3 and 6 h for the final molecules of O2 to be used by the animal (Cox et al., 2010). Certainly by 36 h of anoxia, as indicated by the Hb O2 saturation levels measured in the in vivo anoxia-exposed fish in our study, the last of the blood O2 stores had been used up. Thus, the initial response of the excised heart to anoxia is indeed an initial response to anoxia; there are no behavioural responses that could influence cardiac function as has been suggested for previous in vivo studies (Cox et al., 2010; Forster, 1990; Forster et al., 1992). Furthermore, apart from the relatively small myoglobin-bound O2 reserves in the cardiac cells themselves (Davison et al., 1990), there are no large stores of O2 that could obfuscate the onset of cardiac anoxia. The initial metabolic response of the hagfish heart to anoxia is therefore an immediate depression in metabolic rate – probably as a result of rapid depression of aerobic ATP production – followed by a re-establishment of metabolic rate to approximately routine levels, which is probably due to upregulation of anaerobic glycolysis .

Anaerobic metabolites of the heart and blood

Cardiac glycogen depletion in both the anoxia-exposed excised and in vivo hearts indicates that anaerobic glycolysis is taking place. The rate of glycogen depletion appears to be faster in the excised hearts than the in vivo hearts. This is likely attributable to two aspects of the calorimetry set-up: first, as mentioned above, the onset of anoxia in the excised hearts was probably quicker than in the in vivo hearts; second, the cardiac glycogen stores in the in vivo hearts were probably augmented by the mobilization of glycogen stores from elsewhere in the animal, as seen by the elevated glucose levels of anoxic hagfish.

Tissue lactate levels did not increase significantly in the anoxia-exposed excised hearts, but they did in the in vivo hearts. This could be because the lactate produced by the myocytes of the excised heart was lost into the perfusing saline and thus diluted by its relatively large volume (approximately 25-fold that of the blood volume in vivo). The lactate therefore never accumulated in the excised heart itself, but rather continually moved down its concentration gradient and into the saline over the course of the experiment. In the in vivo exposures, the accumulation of lactate in the plasma from all of the glycolytically active tissues of the body might decrease the rate at which it could move from the heart into the blood. It would therefore accumulate in the heart, and this could be why our results show an accumulation of lactate in the in vivo hearts but not in the excised hearts.

In addition to this accumulation of lactate, the in vivo hearts showed significantly lower glycogen concentrations compared with the normoxic control hearts. These concentrations were so low by the 36 h anoxic time point as to be essentially depleted, a result consistent with Cox et al. (2011), who found that cardiac glycogen stores of in vivo hagfish hearts were depleted after 24 h of anoxia. Like our fish, their fish survived until the 36 h time point, suggesting that the glycogen required by the hearts was sourced from storage tissues such as the liver for at least 12 h and/or was augmented by another fuel besides glycogen.

pH of the heart and blood

Reductions in blood and intracellular pH are concomitant with anaerobic glycolysis (Robergs et al., 2004). Consistent with this, the blood pH of in vivo anoxia-exposed fish was 6.97, 0.84 pH units lower than the blood of the normoxia-exposed fish. Despite this reduction, the pHi of the anoxic hearts was only 0.07 pH units lower than the normoxic hearts, and the pHi of the excised hearts sampled at 1.5 h and 8 h anoxia was almost identical to that seen in the hearts from the in vivo experiment after 36 h of anoxia exposure. This suggests an ability of the hagfish cardiac myocytes to buffer and/or actively regulate intracellular H+ concentrations during a metabolic acidosis, and is consistent with regulation of pHi in the hagfish heart in response to a respiratory acidosis (Baker et al., 2015). The contractile element in cardiac myocytes is extremely sensitive to pH change (Gillis et al., 2000), so this defence of pHi is essential to the maintenance of cardiac function during anoxia exposure.

There may have been more than just glycolytically sourced H+ contributing to the blood acidosis measured in the in vivo anoxia-exposed fish. Although we did not measure it, it is possible that the chambers that held the fish for the duration of the anoxia exposures became hypercarbic (i.e. elevated environmental CO2 levels). Plasma PCO2 measurements of approximately 13 mmHg for the anoxia-exposed fish (compared with approximately 2 mmHg for the normoxia-exposed fish at the same 36 h time point) support this possibility. We believe the source of this CO2 was the fish within the chambers. More precisely, as the flow of water through the chamber was only stopped when PO2 measurements reached 0 mmHg, the source of this CO2 must have been the fish's O2 reserves. Interestingly, by combining an estimation for O2 storage capacity of E. stoutii (137.9 ml O2 kg−1) (Forster et al., 2001) with the total mass of tissue in the chambers, we calculate a PCO2 for the 2.5 litre chambers of 15.1 mmHg, which is very similar to the PCO2 of 13.4 mmHg that we measured in the blood of the anoxic fish. So the hagfish's large O2 stores probably caused the chambers' hypercarbia. But how did the hypercarbia influence the hagfish's response to anoxia? It is possible that an accumulation of H+ in the blood stemming from the hydration of CO2 slowed the rate of H+ diffusion from the tissues into the blood. This would slow the production of lactate in the tissues, and subsequently, the production rates of glycolytically sourced ATP (e.g. Cochran and Burnett, 1996). In any case, the preferred E. stoutii habitats of sandy, muddy burrows (Martini, 1998) and their penchant for feeding within the rotting remains of large, dead vertebrates (Martini, 1998) suggest that when these animals encounter anoxia in the natural environment, they very likely encounter hypercarbia too. Anoxia together with hypercarbia may then be the most ecologically relevant anoxic exposure for these animals.

Summary and perspective

The hagfish heart appears to be an incredibly robust organ capable of maintaining metabolic activity during prolonged periods of anoxia. These characteristics, coupled with a large blood volume to store metabolic wastes, enable cardiac function in the animals to be maintained during extended periods of anoxia with minimal detrimental effects to the animal.

This study represents the first time direct calorimetry has been performed on an excised, functional organ. This technique allowed us to measure the direct and immediate response of the heart to complete anoxia, and highlighted two results in the process that we might not have predicted from in vivo experiments alone. First, there was a significant and immediate depression in cardiac metabolic rate with the induction of anoxia, followed by a re-establishment of metabolic rate to near-routine levels when expressed as ATP production rates. This finding, illustrating that metabolic activity is sustained in the heart during anoxia exposure, is supported by substrate use and the accumulation of metabolic wastes. We suspect that the initial metabolic depression is the result of an immediate arrest of aerobic respiration, while the re-establishment of metabolic rate is the result of the activation of anaerobic metabolism supported by glycogen. And second, calorimetry revealed a maintenance of metabolic heat long after the cardiac glycogen levels were likely to have been depleted, suggesting that the hagfish heart utilizes an alternative fuel during prolonged periods of anoxia. This is clearly an area worthy of further study.

We thank the staff at InSEAS at UBC for technical expertise and Janice Oaks at Fisheries and Oceans Canada for assistance with animal husbandry. The authors would also like to thank Dr Eric Clelland (Bamfield Marine Sciences Center), Alyssa Weinrauch, Alex Clifford and Dr Greg Goss (University of Alberta) for assistance in obtaining the hagfish. The authors would like to thank Dr John Gosline (UBC) for the loan of the calorimeter and for his insight while we were getting it to work with an isolated heart.

Author contributions

T.E.G., M.D.R. and A.P.F. designed the study. T.E.G., M.D.R., G.K.C. and T.S.H. completed the experiments. M.D.R. performed metabolite analyses. T.E.G., M.D.R. and T.S.H. analysed the data and all authors wrote the paper.

Funding

This work was supported by an Natural Sciences and Engineering Research Council of Canada (NSERC) PGS and Zoology Graduate Fellowship to M.D.R. and NSERC Discovery grants to T.E.G., C.J.B., J.G.R. and A.P.F. A.P.F. is supported by the Canada Research Chairs program and C.J.B. by an NSERC Accelerator Supplement.

Altimiras
,
J.
and
Larsen
,
E.
(
2000
).
Non-invasive recording of heart rate and ventilation rate in rainbow trout during rest and swimming. Fish go wireless!
J. Fish Biol.
57
,
197
-
209
.
Baker
,
D. W.
,
Sardella
,
B.
,
Rummer
,
J. L.
,
Sackville
,
M.
and
Brauner
,
C. J.
(
2015
).
Hagfish: champions of CO2 tolerance question the origins of vertebrate gill function
.
Sci. Rep.
5
,
11182
.
Bartholomew
,
G. A.
,
Vleck
,
D.
and
Vleck
,
C. M.
(
1981
).
Instantaneous measurements of oxygen consumption during pre-flight warm-up and post-flight cooling in Sphingid and Saturniid moths
.
J. Exp. Biol.
90
,
17
-
32
.
Bergmeyer
,
H. U.
,
Bergmeyer
,
J. R.
and
Grassl
,
M.
(
1983
).
Methods of Enzymatic Analysis
.
Weinheim, FL
:
Verlag Chemie
.
Boutilier
,
R. G.
,
Heming
,
T. A.
and
Iwama
,
G. K.
(
1984
).
Appendix: physicochemical parameters for use in fish respiratory physiology
.
Fish Physiol.
10
,
403
-
430
.
Clark
,
T. D.
,
Farwell
,
C. J.
,
Rodriguez
,
L. E.
,
Brandt
,
W. T.
and
Block
,
B. A.
(
2013
).
Heart rate responses to temperature in free-swimming Pacific bluefin tuna (Thunnus orientalis)
.
J. Exp. Biol.
216
,
3208
-
3214
.
Clarke
,
E. W.
(
1957
).
Anaerobic heat production of bull spermatozoa
.
Proc. R. Soc. B Biological Sciences
147
,
316
-
331
.
Cochran
,
R. E.
and
Burnett
,
L. E.
(
1996
).
Respiratory responses of the salt marsh animals, Fundulus heteroclitus, Leiostomus xanthurus, and Palaemonetes pugio to environmental hypoxia and hypercapnia and to the organophosphate pesticide, azinphosmethyl
.
J. Exp. Mar. Biol. Ecol.
195
,
125
-
144
.
Cox
,
G. K.
,
Sandblom
,
E.
and
Farrell
,
A. P.
(
2010
).
Cardiac responses to anoxia in the Pacific hagfish, Eptatretus stoutii
.
J. Exp. Biol.
213
,
3692
-
3698
.
Cox
,
G. K.
,
Sandblom
,
E.
,
Richards
,
J. G.
and
Farrell
,
A. P.
(
2011
).
Anoxic survival of the Pacific hagfish (Eptatretus stoutii)
.
J. Comp. Physiol. B
181
,
361
-
371
.
Davison
,
W.
,
Baldwin
,
J.
,
Davie
,
P. S.
,
Forster
,
M. E.
and
Satchell
,
G. H.
(
1990
).
Exhausting exercise in the Hagfish, Eptatretus Cirrhatus: the anaerobic potential and the appearance of lactic acid in the blood
.
Comp. Biochem. Physiol. A Physiol.
95
,
585
-
589
.
Farrell
,
A. P.
and
Stecyk
,
J. A.
(
2007
).
The heart as a working model to explore themes and strategies for anoxic survival in ectothermic vertebrates
.
Comp. Biochem. Physiol. A Mol. Integr. Physiol.
147
,
300
-
312
.
Forster
,
M. E.
(
1990
).
Confirmation of the low metabolic rate of hagfish
.
Comp. Biochem. Physiol. A Physiol.
96
,
113
-
116
.
Forster
,
M. E.
(
1991
).
Myocardial oxygen-consumption and lactate release by the hypoxic hagfish heart
.
J. Exp. Biol.
156
,
583
-
590
.
Forster
,
M. E.
,
Davie
,
P. S.
,
Davison
,
W.
,
Satchell
,
G. H.
and
Wells
,
R. M. G.
(
1988
).
Blood pressures and heart rates in swimming hagfish
.
Comp. Biochem. Physiol. A Physiol.
89
,
247
-
250
.
Forster
,
M. E.
,
Davison
,
W.
,
Axelsson
,
M.
and
Farrell
,
A. P.
(
1992
).
Cardiovascular responses to hypoxia in the hagfish, Eptatretus cirrhatus
.
Respir. Physiol.
88
,
373
-
386
.
Forster
,
M. E.
,
Russell
,
M. J.
,
Hambleton
,
D. C.
and
Olson
,
K. R.
(
2001
).
Blood and extracellular fluid volume in whole body and tissues of the Pacific hagfish, Eptatretus stouti
.
Physiol. Biochem. Zool.
74
,
750
-
756
.
Galli
,
G. L. J.
,
Lau
,
G. Y.
and
Richards
,
J. G.
(
2013
).
Beating oxygen: chronic anoxia exposure reduces mitochondrial F1FO-ATPase activity in turtle (Trachemys scripta) heart
.
J. Exp. Biol.
216
,
3283
-
3293
.
Gillis
,
T. E.
,
Marshall
,
C. R.
,
Xue
,
X. H.
,
Borgford
,
T. J.
and
Tibbits
,
G. F.
(
2000
).
Ca(2+) binding to cardiac troponin C: effects of temperature and pH on mammalian and salmonid isoforms
.
Am. J. Physiol.
279
,
R1707
-
R1715
.
Hansen
,
C. A.
and
Sidell
,
B. D.
(
1983
).
Atlantic hagfish cardiac muscle: metabolic basis of tolerance to anoxia
.
Am. J. Physiol.
244
,
R356
-
R362
.
Hicks
,
J. M.
and
Farrell
,
A. P.
(
2000
).
The cardiovascular responses of the red-eared slider (Trachemys scripta) acclimated to either 22 or 5 degrees C. I. Effects of anoxic exposure on in vivo cardiac performance
.
J. Exp. Biol.
203
,
3765
-
3774
.
Johnsson
,
M.
and
Axelsson
,
M.
(
1996
).
Control of the systemic heart and the portal heart of Myxine glutinosa
.
J. Exp. Biol.
199
,
1429
-
1434
.
Lau
,
G. Y.
and
Richards
,
J. G.
(
2011
).
AMP-activated protein kinase plays a role in initiating metabolic rate suppression in goldfish hepatocytes
.
J. Comp. Physiol. B
181
,
927
-
939
.
Lesser
,
M. P.
,
Martini
,
F. H.
and
Heiser
,
J. B.
(
1997
).
Ecology of the hagfish, Myxine glutinosa L. in the Gulf of Maine I. Metabolic rates and energetics
.
J. Exp. Mar. Biol. Ecol.
208
,
215
-
225
.
Martini
,
F. H.
(
1998
).
The ecology of hagfishes
. In
The Biology of Hagfishes
,
Chap. 5 (ed. J. M. Jorgensen, J. P. Lomholt, R. E. Weber and H. Malte), pp. 57-77. London: Chapman and Hall
.
Portner
,
H.-O.
(
1990
).
Determination of intracellular buffer values after metabolic inhibition by fluoride and nitrilotriacetic acid
.
Respir. Physiol.
81
,
275
-
288
.
Regan
,
M. D.
,
Gosline
,
J. M.
and
Richards
,
J. G.
(
2013
).
A simple and affordable calorespirometer for assessing the metabolic rates of fishes
.
J. Exp. Biol.
216
,
4507
-
4513
.
Robergs
,
R. A.
,
Ghiasvand
,
F.
and
Parker
,
D.
(
2004
).
Biochemistry of exercise-induced metabolic acidosis
.
Am. J. Physiol.
287
,
R502
-
R516
.
Stecyk
,
J. A. W.
,
Stenslokken
,
K.-O.
,
Farrell
,
A. P.
and
Nilsson
,
G. E.
(
2004
).
Maintained cardiac pumping in anoxic crucian carp
.
Science
306
,
77
.
Stecyk
,
J. A. W.
,
Galli
,
G. L.
,
Shiels
,
H. A.
and
Farrell
,
A. P.
(
2008
).
Cardiac survival in anoxia-tolerant vertebrates: an electrophysiological perspective
.
Comp. Biochem. Physiol. C
148
,
339
-
354
.
Tucker
,
V. A.
(
1967
).
Method for oxygen content and dissociation curves on microliter blood samples
.
J. Appl. Physiol.
23
,
410
-
414
.
van Waversveld
,
J.
,
Addink
,
A. D. F.
and
Van Den Thillart
,
G.
(
1989
).
Simultaneous direct and indirect caloimetry on normoxic and anoxic goldfish
.
J. Exp. Biol.
142
,
325
-
335
.
Wells
,
R. M. G.
,
Forster
,
M. E.
,
Davison
,
W.
,
Taylor
,
H. H.
,
Davie
,
P. S.
and
Satchell
,
G. H.
(
1986
).
Blood oxygen transport in the free-swimming hagfish, Eptatretus cirrhatus
.
J. Exp. Biol.
123
,
43
-
53
.
Wilson
,
C. M.
(
2014
).
Cardiac control in the Pacific hagfish (Eptatretus stoutii)
.
PhD thesis
, pp.
170
.
The University of British Columbia: University of British Columbia
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information