SUMMARY

The physiological basis of the aerobic dive limit (ADL), the dive duration associated with the onset of post-dive blood lactate elevation, is hypothesized to be depletion of the muscle oxygen (O2) store. A dual wavelength near-infrared spectrophotometer was developed and used to measure myoglobin (Mb) O2 saturation levels in the locomotory muscle during dives of emperor penguins (Aptenodytes forsteri). Two distinct patterns of muscle O2 depletion were observed. Type A dives had a monotonic decline, and, in dives near the ADL, the muscle O2 store was almost completely depleted. This pattern of Mb desaturation was consistent with lack of muscle blood flow and supports the hypothesis that the onset of post-dive blood lactate accumulation is secondary to muscle O2 depletion during dives. The mean type A Mb desaturation rate allowed for calculation of a mean muscle O2 consumption of 12.4 ml O2 kg–1 muscle min–1, based on a Mb concentration of 6.4 g 100 g–1 muscle. Type B desaturation patterns demonstrated a more gradual decline, often reaching a mid-dive plateau in Mb desaturation. This mid-dive plateau suggests maintenance of some muscle perfusion during these dives. At the end of type B dives, Mb desaturation rate increased and, in dives beyond the ADL, Mb saturation often reached near 0%. Thus, although different physiological strategies may be used during emperor penguin diving, both Mb desaturation patterns support the hypothesis that the onset of post-dive lactate accumulation is secondary to muscle O2 store depletion.

Early studies on forcibly submerged seals revealed a series of physiological responses, referred to as the dive response (Scholander, 1940; Scholander et al., 1942), that support the hypothesis that muscle is the primary source of post-dive blood lactate accumulation. During forced submersions, muscle was isolated from circulation because of severe bradycardia and peripheral vasoconstriction (Bron et al., 1966; Scholander, 1940; Scholander et al., 1942). Muscle metabolism relied on the intrinsic myoglobin (Mb)-O2 store and once that store was consumed, lactate began to accumulate in the muscle (Scholander, 1940; Scholander et al., 1942). During the post-submersion tachycardia, washout of lactate from the muscle elevated blood lactate concentration, reinforcing the concept that muscle was isolated from the circulation and was the primary source of post-dive blood lactate (Scholander, 1940; Scholander et al., 1942).

Recent studies at an isolated dive hole in McMurdo Sound, Antarctica, suggest that freely diving emperor penguins (Aptenodytes forsteri Gray 1844) may have a dive response similar to the classic dive response demonstrated in forced submersion studies. During the initial descent phase of dives, when stroke frequency (muscle workload) is highest, there is an immediate reduction in heart rate, which is consistent with reduced muscle perfusion (Meir et al., 2008; van Dam et al., 2002). Similarly, venous partial pressure of O2 (PO2) values often increase during the descent phase of emperor penguin dives, at a time when significant O2 extraction by working muscle would be expected to lower venous PO2 if MBF were maintained (Ponganis et al., 2007).

Accordingly, we hypothesized that muscle is ischemic during diving and that depletion of the muscle O2 store with subsequent muscle lactate accumulation is the primary mechanism responsible for post-dive blood lactate accumulation in dives beyond the previously measured ADL of 5.6 min (Ponganis et al., 1997b). We predicted that Mb desaturation would follow a monotonic decline until almost complete desaturation by 5.6 min, after which it would plateau near 0% saturation until the end of the dive. To test our hypothesis, we attached a custom-built, near-infrared spectrophotometer (NIRS) instrument to emperor penguins diving at an isolated dive hole. Our goals were to: (1) assess the magnitude and rate of muscle O2 depletion during dives, (2) describe patterns of Mb desaturation, (3) examine Mb desaturation profiles for indications of muscle perfusion or hypoperfusion during dives and (4) determine muscle O2 consumption during dives.

MATERIALS AND METHODS

Instrument design

In order to measure Mb desaturation in the emperor penguin muscle during diving, we developed a small, dual wavelength NIRS instrument. NIRS instruments have frequently been used to measure Mb or Hb saturation (Delpy and Cope, 1997; Guyton et al., 1995; Jobsis, 1977; Jobsis et al., 2001; Sako et al., 2001). The absorption of transmitted NIR radiation by specific chromophores, such as Mb, is measured from the reduction in intensity of reflected radiation due to absorption by Mb. The attenuation in radiation intensity is proportional to the absorption by Mb, which can then be related to the NIR absorption spectra of deoxy-Mb and oxygenated Mb (oxy-Mb) to provide a relative degree of oxygenation. The ratio of two wavelengths, 760 nm (where absorption by deoxy-Mb is significantly higher than absorption by oxy-Mb) and 800 nm (an isosbestic point, where absorption coefficients of both deoxy-Mb and oxy-Mb are equal), is linearly correlated to the relative saturation level of Mb (Guyton et al., 1995; Jobsis, 1977; Mancini et al., 1994).

Probes were constructed using a small circuit board (19×5×3 mm, <1 g) with two surface-mount LEDs (SMT760 & SMT810, Epitex, Inc., Kyoto, Japan) positioned in the middle of the board 3 mm apart (Fig. 1A). The 50 Hz sampling rate, alternately powering the two LEDs, resulted in a 2.8 ms pulse width with 17.2 ms between pulses. With a rise and fall time of less than 100 ns for each LED, it was predicted that the sampling cycle would work well with no overlap of pulses. Two surface-mount photodiodes (PD006-SMT, Epitex, Inc.), for receiving reflectance data, were soldered on each end of the board, 3 mm from each LED (Fig. 1A). Rise and fall times for the photodiodes were 6 ns. Probes were sealed using a two-component epoxy with a spectral transmission of >95% in the NIR range (Epo-Tek 302-3M, Epoxy Technology, Inc., Billerica, MA, USA).

Instrument validation

Prior to using the NIRS instrument in the field, the linearity of the NIRS reflectance readings at different saturation values was verified. In the NIR range, absorption spectra for both Hb and Mb are essentially identical; therefore, Hb was used as a proxy for Mb in the linearity verification experiment. Rat (Rattus norvegicus) whole blood was mixed in a tonometer attached to a Wösthoff gas-mixing pump to obtain 0 and 100% Hb saturation (Tonometer 237, Instrumentation Laboratories, Bornheim, Germany; H. Wösthoff KG, Bochum, Germany). Blood samples for Hb saturation values of 25, 50 and 75% were obtained using a volumetric mixing technique (Scheid and Meyer, 1978). Blood gas analysis (i-STAT blood gas analyzer, Abbott Point of Care, Princeton, NJ, USA) was used to determine PO2 values for samples at each saturation percentage. Accuracy of the mixing technique was confirmed by comparing PO2 values for each sample with PO2 values for each saturation percentage predicted by the rat O2–Hb dissociation curve (Schmidt-Nielsen, 1997). A 2 ml blood sample at each Hb saturation value was transferred anaerobically from the tonometer to a spectrophotometric cuvette. The probe was attached to the side of the cuvette and then covered to prevent light penetration. Readings were immediately taken from the NIRS instrument for each sample.

Field study approach

Non-breeding emperor penguins (8–15 birds per season, N=16 for this experiment, 20–26 kg) were captured on the sea ice of McMurdo Sound, Antarctica, in the austral spring of 2007 and 2008. Penguins were transported to a sea ice camp (Penguin Ranch) set up on McMurdo Sound with two isolated dive holes and a penguin corral as previously described (Kooyman et al., 1992). Penguins dived freely at the isolated dive holes and foraged on sub-ice fish and squid as verified by underwater visual observations and guano deposits. Once experiments were completed, and within 6 weeks of capture, all penguins were released at the sea ice edge. All procedures were approved under a UCSD Animal Subjects Committee protocol and a US Antarctic Treaty permit.

Instrument attachment

The NIRS probe was implanted on the surface of the right pectoral muscle of penguins under general anesthesia (Kooyman and Ponganis, 1994; Kooyman et al., 1992). A small incision (5 cm), located mid-wing level and 10 cm medial from the black–white feather margin, was made in the skin. Muscle was exposed by careful blunt dissection to avoid trauma and bleeding. The probe was sutured to the muscle with three silk sutures previously attached to the back of the probe. Probes were implanted in areas free of major blood vessels. After attachment to the muscle surface and verification of probe function with the recorder, the incision was closed in two layers (subcutaneous and skin) with 2–0 prolene suture. The cable exiting the skin was secured to the feathers with TesaTM tape and LoctiteTM glue. The NIRS recorder was attached to the feathers of the mid-back with a 5 min epoxy glue (Devcon, Danvers, MA, USA), a VelcroTM strip and cable ties as previously described (Stockard et al., 2005). A time–depth recorder (TDR) (MK9, Wildlife Computers, Redmond, WA, USA) (sensitive to 0.5 m, 30 g, 6.5×1.7×1.7 cm, sample rate 1 Hz) was attached above the NIRS instrument midline between the wings as previously described (Stockard et al., 2005). After recovery overnight, instrumented birds were allowed to dive freely for 1 to 2 days. The NIRS instrument, probe and TDR were then removed under anesthesia and all study penguins, observed carefully post-procedure, resumed regular daily diving shortly thereafter.

Zero calibration of probe

During the probe removal procedure, in order to obtain a zero calibration, a portion (∼15×6×4 mm, 200 mg) of the pectoral muscle with the probe still attached to it was excised. The excised muscle sample and probe were immediately placed in a watertight reclosable plastic bag wrapped in black electric tape, which was then kept in a 38°C water bath. An opened corner of the bag allowed exit of the probe cable, which was still connected to the NIRS recorder. Via small Tygon® tubing, 100% nitrogen gas (N2, ultrahigh purity, >99.999%) was bubbled through a saline-filled test tube (38°C) and into the open corner of the bag at a rate of 0.1 l min–1. The NIRS instrument continued to record data until there were no further changes in the reflectance signals in order to ensure a 0% Mb saturation value for each penguin muscle sample. NIRS reflectance ratio data were converted to Mb saturation using the 0% saturation value and an assumed 100% Mb saturation value prior to diving.

NIRS probe evaluation

Effects of blood flow and Hb oxygenation on NIRS probe output from muscle were evaluated on two penguins under general isoflurane-O2 endotracheal anesthesia (Kooyman et al., 1992). The NIRS probe was attached to the pectoral muscle as described above. The birds were also catheterized percutaneously in the wings with a 20 g, 4.5 cm radial artery catheter (RA-04020, Arrow International, Reading, PA, USA) and a 16 g, 4.5 cm brachial vein catheter (BD InsyteTM 381257, Becton Dickinson, Sandy, UT, USA). Arterial blood pressure was transduced by a Hewlett-Packard blood pressure transducer system (78304A/78205D) calibrated by a manometer with the transducer at the level of the penguin's heart. Data were recorded and analyzed in mm Hg, and are expressed so in the text; pressure axes in graphs are labeled in both mm Hg and kPA, using the conversion 0.133 kPA=1 mm Hg. The electrocardiogram (ECG) was recorded with an ECG recorder (RespI/ECG, UFI, Morro Bay, CA, USA) attached to three subcutaneous electrodes placed midline in the back above and below the heart, and laterally in the chest below the heart. Blood pressure and ECG outputs were recorded on a personal computer with Acqknowledge software and a Biopac MP100 interface system (Biopac, Goleta, CA, USA).

Effects of presumed changes in MBF on the NIRS probe output were evaluated with intravenous injections of alpha and beta sympathomimetic agents (phenylephrine and isoproterenol, respectively).

In order to evaluate the effect of Hb saturation on the NIRS probe output, ventilation was maintained with both spontaneous respirations and manually assisted ventilation in order to maintain blood pH in the 7.4 to 7.5 range. The inspired O2 fraction (FIO2) was changed with addition of N2via a second flow meter to the anesthesia circuit and was monitored with an O2 analyzer (Vascular Technology, Inc., Nashua, NH, USA). Arterial and venous blood PO2 and pH were analyzed with an i-STAT blood gas analyzer. Hb saturation was calculated from the pH 7.5 O2–Hb dissociation curve of the emperor penguin (Meir and Ponganis, 2009).

Probes and catheters were removed at completion of the studies, and the birds were released back into the penguin corral after recovery from anesthesia.

Data processing and statistics

Data were processed, graphed and statistically analyzed using Origin (OriginLab Corp., Northampton, MA, USA), SPSS (version 11.5, SPSS, Inc. Chicago, IL, USA) and a custom-developed MATLAB program (The MathWorks, Inc., Natick, MA, USA). NIRS data were processed using several custom MATLAB scripts in which the baseline for each LED was determined and stroke artifact was filtered out (Fig. 2). The 760 nm and 810 nm data were smoothed with a 20-point moving average and the ratio 760 nm/810 nm was calculated. Confirmation that movement artifact from the NIRS signal was an accurate record of strokes was made in several ways. First, NIRS movement artifact data were compared with stroke movements recorded by a two-dimensional accelerometer (UME-D2GT, Little Leonardo Co., Tokyo, Japan) during dives in one penguin, which confirmed that artifact was consistent with stroking. Second, manual movement of the wing while the penguin was under anesthesia was confirmed in the NIRS signal. Third, the pattern and rate of the movement artifact during dives were similar to the stroke frequency profiles observed in a prior study (van Dam et al., 2002). Strokes were then determined from a custom-made peak detection script and visually confirmed. Stroke frequency for each dive was calculated from the total number of strokes for each dive divided by dive duration.

(A) Example of data recorded by the NIRS instrument and time–depth recorder during a dive. (B) A 5 s excerpt from the same dive showing stroke movement artifact and NIRS baseline between strokes. Depth is shown in blue and NIRS reflectance signal is shown in black.

The relationship between NIRS reflectance ratio from the NIRS instrument and Hb saturation was assessed with simple linear regression in the validation study. End-of-dive Mb saturation values and Mb desaturation rate were determined from converted NIRS data. Regression analysis was used to assess the relationship between end-of-dive saturation and dive duration and Mb desaturation rate and stroke frequency. Statistical significance was assumed at P≤0.05. Values are expressed as means ± s.d. unless otherwise noted.

NIRS probe evaluation

In the first experiment, to examine the effects of blood flow on the NIRS signals, the baseline systolic blood pressure was 100 to 120 mm Hg and heart rate was approximately 80 beats min–1. After a 0.6 μg intravenous injection of isoproterenol, systolic blood pressure increased to 160 mm Hg and heart rate increased to 140 beats min–1, but no distinct changes in either the 760 nm or the 810 nm NIRS signals occurred. After return of heart rate and blood pressure to baseline, 100 μg of phenylephrine was injected intravenously. Systolic blood pressure rose dramatically to above 220 mm Hg and heart rate decreased transiently to below 40 beats min–1. Again, no significant changes in NIRS signals were observed.

In the second experiment, to examine the effects of changes in Hb oxygenation on the NIRS signals, FIO2 was reduced from 100 to 50% for 10 min and then to 20% for 4 min while monitoring the ECG record and taking periodic blood gas samples. Heart rate rose slightly during the reduction in FIO2, but was quite variable during the entire procedure (Fig. 4). However, there were no associated changes or trends in either NIRS reflectance signal (Fig. 4). After approximately 4 min at the 20% FIO2, arterial and venous Hb saturation had dropped approximately 40% (from 99 and 97% to 59 and 61%, respectively). During this period, there were minimal, inconsistent changes in the NIRS signals (Fig. 4).

Dive behavior

NIRS data were obtained from three birds for a total of 50 dives >2 min duration (Table 1). Muscle biopsies and subsequent zero calibrations were successfully performed in all three birds. The difficulty of obtaining Mb desaturation profiles on a primary locomotory muscle is demonstrated by the number of birds from which we were unable to recover usable NIRS data. Mb desaturation data were not obtained from 13 birds equipped with the NIRS instrument due to: (1) breakage of the medical cable connecting the probe to the recorder during diving, (2) unsuccessful zero calibration experiments and (3) excessive movement artifact in at least one of the LED's reflectance data, which could not be resolved with filtering.

Dive duration and maximum depth data of the three birds are reported in Table 1. Dive durations from all dives ranged from 2.3 to 11.4 min (Fig. 5). Thirty-one dives were equal to or longer than the 5.6 min ADL. Maximum dive depth for all dives ranged from 7 to 64 m (Fig. 5).

Desaturation rates, stroke frequency and recovery times to 75% resaturation are reported in Table 1. Grand mean stroke frequency for all three birds was 0.66±0.08 Hz. Grand mean Mb desaturation rate during dives was 11.03±2.4% min–1. Mean recovery times ranged from 0.5 min to almost 2 min, with a grand mean of 1.03±0.76 min.

Mb desaturation patterns

Two distinct Mb desaturation patterns were observed in the birds, a monotonic decline (type A) and a mid-dive plateau pattern (type B). Dives were classified into either type A or type B based on visual verification of the Mb desaturation pattern. Six dives could not be classified as either type A or type B. Bird 2 had primarily type A dives, whereas birds 7 and 11 had predominantly type B dives. Dive characteristics, Mb desaturation rates and recovery times of type A and B dives are summarized in Table 2.

Dive characteristics, Mb desaturation rate and recovery time of type A and type B desaturation profiles

End-of-dive Mb saturation was linearly related to dive duration in type A and type B dives (Fig. 6). In type A dives, end-of-dive Mb saturation was close to 0% in dives near the ADL (Fig. 6A). However, in type B dives, end-of-dive Mb saturation was not significantly depleted until dive durations of 8 to 10 min (Fig. 6B). At 5.6 min into type B dives, Mb saturation was close to 50% (Fig. 6B).

Type A dives

In type A dives, Mb saturation generally followed a monotonic decline throughout the dive (Fig. 7A). Eight dives were classified as type A. The desaturation rate was often slightly higher during the initial descent, but the slope then remained fairly constant throughout the rest of the dive. There were only minor changes in slope or desaturation rate during type A dives (Fig. 7A). In a few longer type A dives, the slope of the desaturation rate leveled off at the end (Fig. 7A). Mean desaturation rate in type A dives was 14.4±3.8% min–1. Mean stroke frequency in type A dives was 0.74±0.04 Hz (Table 2). Instantaneous stroke frequency ranged from 0.4 to 2.0 Hz, with higher stroke frequencies associated with descents and ascents during the dive (Fig. 7A). For dives with 5 to 6 min durations (near the ADL), mean end-of-dive Mb saturation was 7.1±5.6% (N=3, range=0.6–10.7%). Stroke frequency was not a significant predictor of Mb desaturation rate (r2=0.26, P=0.2).

Type B dives

Thirty-six dives were classified as having a type B desaturation pattern. During initial descent in type B dives, although stroke frequency was high (1–2 Hz), Mb desaturation rate was sometimes fast, as in type A dives, and sometimes moderate (Fig. 7B, Fig. 8). During the middle phase of the dive, Mb desaturation rate was slow and, sometimes, nearly 0 in the middle of the dive (Fig. 7B, Fig. 8) while stroke frequency was 0.4 to 0.6 Hz and fairly constant (Fig. 7B). Finally, during the ascent phase, the desaturation rate increased but stroke frequency remained low, 0.4–0.6 Hz (Fig. 7B), with higher stroke frequencies only observed in some dives at the end of the ascent. Mean desaturation rate in type B dives was 9.8±2.4% min–1. Mean stroke frequency in type B dives was 0.64±0.09 Hz. For dives with 5 to 6 min durations (near the ADL), mean end-of-dive saturation was 44.8±6.5% (N=6, range= 32.8–50.1%). Stroke frequency was a weak but statistically significant predictor of Mb desaturation rate (r2=0.18, P=0.01).

Type A versus type B dives

The differences between type A and type B Mb desaturation rates, stroke frequencies, dive duration and maximum dive depth were statistically significant (Mann–Whitney U-test, Z=–2.891, P=0.004; Z=–3.228, P=0.001; Z=–2.800, P=0.005; Z=–4.202, P=0.000, respectively; Table 2). There was no characteristic dive profile for either type A or type B dives. Dives with intra-dive ascents to 2–4 m depth followed by descents to near maximum dive depth occurred in both dive types. The difference in recovery to 75% Mb saturation was not statistically different between type A and type B dives (Mann–Whitney U-test, Z=–0.259, P=0.8).

NIRS probe evaluation studies

NIRS signals may include reflectance data from both Hb and Mb saturation because the absorption spectra of Hb and Mb are essentially the same. Thus, the relative contribution of blood flow, Hb oxygenation and Mb oxygenation to NIRS signals may vary (Mancini et al., 1994; Tran et al., 1999; Wilson et al., 1989). Although the NIRS probe in the present study was attached in areas free of visible blood vessels, the potential contributions of Mb and Hb saturation to the NIRS signals were unknown. Two NIRS instrument evaluation experiments were undertaken on anesthetized birds to address this question. During isoflurane anesthesia, MBF is preserved near levels at rest in pigs and horses, and is even elevated in humans (Lundeen et al., 1983; Raisis et al., 2000; Stevens et al., 1971). The first evaluation experiment was designed to examine the effect of different concentrations of Hb on the NIRS signals by changing the blood flow. After injection of 0.6 μg of isoproterenol, a beta-sympathomimetic agent that increases heart rate and cardiac output and causes peripheral arterial vasodilation, no significant changes in the NIRS signals were observed. If Hb saturation were a significant component of the reflectance signals, the NIRS signals would have been expected to change concomitantly with the observed increase in heart rate and systolic blood pressure. Phenylephrine, an alpha-adrenergic agent, vasoconstricts peripheral vessels and raises blood pressure. The lack of any consistent trend in the NIRS reflectance data, despite systolic blood pressure rising to above 220 mm Hg and heart rate falling to 40 beats min–1 immediately after a 100 μg phenylephrine injection, also confirmed that the NIRS data reflected primarily Mb saturation and not changes in the quantity of Hb present due to changes in MBF. Furthermore, because of the transient effect of phenylephrine, the low metabolic rate of muscle at rest and the high O2 content of penguin muscle, one would not expect to see a large decline in Mb saturation during the transient hemodynamic changes after the phenylephrine injection.

The second evaluation experiment demonstrated that changes in Hb oxygenation due to changes in FIO2 had minimal effect on the NIRS reflectance signals. If the NIRS signals were strongly affected by Hb saturation, the 760 nm reflectance signal would decrease as Hb saturation decreased. However, the 760 nm signal did not show a decreased reflectance; rather, the signal varied slightly in both directions despite a 40% drop in arterial and venous Hb saturation (Fig. 4). As a result of these findings, we concluded that the NIRS reflectance signals predominantly represent Mb saturation during dives.

Type A Mb desaturation profiles

In type A Mb desaturation profiles, the initial descent period was often characterized by a sharp decrease in saturation (Fig. 7A). This sharp decrease generally coincided with high stroke frequency and, thus, high muscle workload (van Dam et al., 2002). After the initial descent, there was a consistent monotonic decline of Mb saturation with only minor changes in desaturation rate throughout the rest of the dive. Small changes in desaturation rate during the dive may be related to variations not just in stroke frequency, but also in stroke amplitude and stroke thrust throughout the dive. However, the overall rate of decline in saturation during a dive appeared constant (Fig. 7A). Individual dive Mb desaturation rates did not correlate with stroke frequency during a dive. We suspect that, in addition to stroke frequency, stroke amplitude and thrust will also contribute to muscle metabolic demand during dives.

Several aspects of type A dives are consistent with isolation of muscle from the circulation during dives. First, the overall Mb desaturation rate (14.4% min–1) is consistent with rapid Mb desaturation observed in early forced submersion studies (Scholander et al., 1942). Second, the Mb desaturation rate is considerably higher than median muscle desaturation rates in Weddell seals (5.1% min–1 for short dives and 2.5% min–1 for dives >17 min) (Guyton et al., 1995). In that study, it was proposed that those desaturation rates were consistent with maintenance of MBF during dives (Guyton et al., 1995). Third, in type A profiles, no detectable increases in Mb saturation were observed during dives, in contrast to observations in the Weddell seal study (Guyton et al., 1995). Further, the muscle O2 consumption (discussed below) based on the Mb desaturation rate is consistent with the muscle metabolic rate (12–17 ml O2 kg–1 min–1) calculated based on the O2 capacity of an isolated pectoralis–supracoracoideus muscle complex and the ADL range of 5–7 min (Ponganis et al., 1997a). Finally, Mb saturation reached low values and plateaued near the ADL in longer dives (5–6 min). Therefore, we conclude that type A Mb desaturation patterns are consistent with our original hypothesis that the muscle O2 store is isolated from the circulation during these dives and is significantly or completely depleted at the ADL.

Prior research on emperor penguins has also suggested that muscle is isolated from the circulation during dives. During the initial descent phase of the dive, when stroke frequency and, therefore, muscle workload is at its peak (Sato et al., 2002; van Dam et al., 2002), blood O2 extraction by muscle should be high if there is MBF. However, venous PO2 values often increased rather than decreased during the descent phase of dives (Ponganis et al., 2007; Ponganis et al., 2009). Similarly, penguins often experienced a striking decline in heart rate during early descent (Meir et al., 2008). Again, this is not consistent with maintenance of MBF at that time. Further, the lack of correlation between stroke frequency and heart rate during dives and the rate and pattern of pectoral muscle temperature changes during dives are also consistent with a severe reduction or cessation of MBF (Meir et al., 2008; Ponganis et al., 2003). Finally, although post-dive blood lactate concentrations were elevated in dives 5.6 min or longer (Ponganis et al., 1997b), intra-dive blood samples demonstrated no significant lactate accumulation as far as 10.5 min into a 12.8 min dive (Ponganis et al., 2009). The lack of elevated lactate as late as 10.5 min into a dive suggests that muscle is isolated from the circulation and any accumulated lactate does not wash out into blood until after a dive, when circulation presumably returns to muscle. In the post-dive period, it is also notable that Mb resaturation is extremely rapid, with mean time to return to 75% saturation near 1 min for both type A and type B dives. Rapid resaturation of Mb also occurs in the post-apneic period of sleeping elephant seals (Ponganis et al., 2008). Such a quick return in Mb saturation is consistent with the post-dive tachycardia and rapid replenishment of air sac and blood O2 stores in emperor penguins (Meir et al., 2008; Ponganis et al., 2009; Stockard et al., 2005).

Muscle O2 consumption during dives

If muscle is isolated from the circulation during type A dives and lactate does not accumulate in dives <ADL, then the desaturation rate of the Mb-O2 store can be used to calculate mean muscle O2 consumption during dives. The mean Mb desaturation rate for type A dives was 14.4% min–1. Assuming a Mb concentration of 6.4 g 100 g–1 muscle (Ponganis et al., 1997a) and an O2-binding capacity of 1.34 ml O2 g–1 Mb, mean muscle O2 consumption is 12.4± 3.3 ml O2 kg–1 muscle min–1. It should be noted that this value of 12.4 ml O2 kg–1 muscle min–1 was estimated from dives of less than 6 min duration with a mean stroke frequency near 0.74 Hz. Based on this stroke frequency, the cost per stroke is 1.7 ml O2. Differences in stroke frequency, stroke amplitude and stroke thrust will all contribute to variation in muscle O2 consumption and stroke cost.

Type B Mb desaturation profiles

In type B dives, Mb desaturation rate during the initial descent was either rapid, as in type A dives, or moderate (Fig. 7B, Fig. 8). However, as descent continued, desaturation rates slowed significantly, often leveling off completely until the ascent phase, when desaturation rates increased. We suggest that the mid-dive plateaus in type B Mb desaturation profiles are secondary to maintenance of muscle perfusion during these segments of type B dives. During the plateau periods, which frequently lasted several minutes or more, stroke frequency was constant (0.4–0.6 Hz), indicating that muscle work continued. The lack of change in Mb saturation implies the muscle O2 store was supplemented by the circulation during these periods. During the ascent phase, despite continued stroking at the same frequency, Mb desaturation rate again increased, implying that MBF was reduced during this portion of the dive.

A plateau phase of Mb desaturation, as occurred in type B dives, is similar to that during sleep apnea in elephant seals and during exercise of humans to maximum O2 consumption (Ponganis et al., 2008; Richardson et al., 1995b). In both situations, the decrease in Mb saturation should support blood-to-muscle O2 flux by maintaining or enhancing the blood-to-muscle PO2 gradient. These type B dive profiles were also similar to Mb desaturation profiles in two Weddell seals, which demonstrated a slow decline at the beginning of some dives and then a flat slope for over half of the dive, which suggested blood-to-muscle O2 transfer from the Hb-O2 store (Guyton et al., 1995). However, type B desaturation patterns in emperor penguins are not entirely analogous to results from the prior study on Weddell seals. The increased Mb desaturation rate at the end of type B dives in penguins was not observed in Weddell seals (Guyton et al., 1995).

In addition, the Mb-O2 store in Weddell seals never approached complete depletion. Even in dives longer than the Weddell seal ADL, including a dive of 27 min, Mb was typically still 40–60% saturated (Guyton et al., 1995). In contrast, saturation values in long type B dives often declined to below 10% (Fig. 8). Thus, type B dives are still supportive of the concept that the onset of post-dive lactate accumulation is secondary to muscle O2 depletion.

Muscle blood flow during type B dives

If muscle is perfused during these plateau periods of type B dives, MBF can be estimated based on mean muscle O2 consumption (12.4 ml O2 kg–1 muscle min–1) previously calculated from type A dives. Because muscle continues to work during this portion of the dive, but Mb saturation does not significantly change, the blood O2 store likely supplies the required 12.4 ml O2 kg–1 muscle min–1. Assuming a pectoralis–supracoracoideus muscle mass of 6.25 kg for a 25 kg penguin, 77.5 ml O2 min–1 is required. With a typical arterio–venous (A–V) O2 content difference of 5 mlO2 dl–1 (Ponganis et al., 2007), blood flow to muscle would be approximately 1550 ml min–1 or 250 ml kg–1 muscle min–1. Given a stroke volume of 50 ml (Kooyman et al., 1992), 1500 ml min–1 of blood flow would require 30 beats min–1 in heart rate. This is easily achieved during the transient increases in heart rate frequently observed during the early to mid portions of dives of emperor penguins (Meir et al., 2008). Besides shunting oxygenated blood through peripheral A–V shunts as has been hypothesized during this time period (Ponganis et al., 2009), we propose that emperor penguins also have the option to perfuse muscle instead. Such plasticity in vascular responses could account for the wide variation in previously reported venous PO2 profiles in diving emperor penguins (Ponganis et al., 2009).

O2 consumption during dives: 6.8 ml O2 kg–1 min–1

Mean O2 store depletion rates have now been calculated in the three O2 stores of the emperor penguin. The sum of these mean O2 store depletion rates provides the first calculation of total body O2 consumption during a dive: 6.8 ml O2 kg–1 min–1. It must be emphasized that depletion rates of O2 stores are highly variable among individual dives, and that this value of 6.8 ml O2 kg–1 min–1 is an overall mean value. Furthermore, longer dives will have lower total body O2 consumption. For example, as previously estimated on the basis of complete depletion of all O2 stores (Ponganis et al., 2010), the O2 consumption during a 23.1 min dive was 2.5 ml O2 kg–1 min–1. In addition, even for dives of similar duration, there is a large range of end-of-dive air sac O2 fractions, blood Hb saturations and muscle Mb saturations. Such variability may be related to differences in dive depth and buoyancy, stroke frequency, stroke amplitude, stroke thrust and diving respiratory air volume of deep versus shallow dives.

This mean value of 6.8 ml O2 kg–1 min–1 now allows an estimation of the relative contributions of the individual O2 stores to metabolic rates during dives predominantly under 10 min duration (based on dive durations in studies on air sac PO2 and blood PO2 and the present study) (Meir and Ponganis, 2009; Stockard et al., 2005). The depletion of the primary locomotory muscle O2 store provides 3.1 mlO2 kg–1 min–1 to the total body O2 consumption. The depletion of the non-locomotory muscle O2 store (muscle mass 13% of body mass) (Ponganis et al., 1997a) at a rate of 4 ml O2 kg–1 muscle min–1 (i.e. near a resting rate) would contribute another 0.5 ml O2 kg–1 min–1 to the total body O2 consumption. Thus, the entire muscle mass would contribute 3.6 ml O2 kg–1 min–1 to total body metabolic rate during a dive. This accounts for over half the average total body O2 consumption during diving. The respiratory and blood O2 stores account for 31% (2.1 ml O2 kg–1 min–1) and 16% (1.1 ml O2 kg–1 min–1), respectively, of the total body diving O2 consumption (Meir and Ponganis, 2009; Stockard et al., 2005).

An O2 depletion rate of 6.8 ml O2 kg–1 min–1 supports frequent suggestions that diving O2 consumption is extremely low in emperor penguins (Kooyman and Ponganis, 1994; Nagy et al., 2001). This mean value is approximately one-third of both the field metabolic rate measured at the isolated dive hole with doubly labeled water and the lowest rate measured in penguins swimming in a flume (20 and 20.7 ml O2 kg–1 min–1, respectively) (Kooyman and Ponganis, 1994; Nagy et al., 2001). It is equivalent to resting metabolic rate measured in emperor penguins floating in a flume (6.2–6.7 ml O2 kg–1 min–1) (Kooyman and Ponganis, 1994) or standing in their thermoneutral zone (Dewasmes et al., 1980; LeMaho et al., 1976; Pinshow et al., 1976). This mean O2 consumption is also slightly greater than the resting rate predicted by allometric equations relating metabolic rate to body mass (Aschoff and Pohl, 1970). From this analysis, the actual O2 depletion rate is close to resting values, confirming that diving is not costly in emperor penguins.

CONCLUSIONS

Two different Mb desaturation patterns in the locomotory muscle of diving emperor penguins were revealed in this study. These patterns are consistent with the high level of variability observed in physiological parameters previously measured in freely diving penguins, including air sac PO2 profiles, venous and arterial PO2 profiles, heart rate, stroke frequency and body temperatures.

The type A pattern, a monotonic decline in Mb desaturation, is reminiscent of results from early forced submersion studies and suggests that no MBF occurred during the dives. Complete or near-complete depletion of the primary locomotory muscle O2 store at the previously measured ADL supports the concept that muscle is the primary source of the post-dive blood lactate elevation in dives beyond the ADL. Using the Mb desaturation rate determined in these dives, a mean muscle O2 consumption was calculated at 12.4 ml O2 kg–1 muscle min–1. This low value demonstrates the highly efficient locomotory cost of diving in emperor penguins. We also calculated the first estimate of diving O2 consumption in a freely diving higher vertebrate from O2 store depletion rates measured in this study and previous studies. This value, 6.8 ml O2 kg–1 min–1, is near measured rates of resting O2 consumption in emperor penguins and supports the concept that the bradycardia of diving and efficient hydrodynamics primarily contribute to the slow depletion of O2 stores and to the dive capacity of emperor penguins.

The type B pattern is a more complex Mb desaturation pattern with either a mid-dive plateau or a slow Mb desaturation rate. This suggests maintenance of MBF for at least a segment of type B dives. It was estimated that an MBF of 250 ml kg–1 muscle min–1 would maintain a constant Mb saturation during the mid-dive plateau period. This amount of blood flow would require ∼30 beats min–1 in heart rate, a value consistent with the elevation in heart rate during the first few minutes of dives (Meir et al., 2008). Despite blood O2 supplementation of the muscle O2 store, extensive depletion of the muscle O2 store in dives beyond the ADL still supports the concept that the onset of post-dive lactate accumulation is secondary to muscle O2 depletion.

The two Mb desaturation patterns revealed in this study, and previously observed large variations in venous PO2 profiles, also reinforce the idea that peripheral vascular responses and blood flow patterns are plastic during the early to middle portions of emperor penguin dives (Ponganis et al., 2009). These results also suggest emperor penguins have at least two distinct physiological strategies while diving, to either perfuse muscle to supplement the muscle O2 store, or to not perfuse muscle but rather utilize peripheral A–V shunts to optimize the blood O2 store.

ACKNOWLEDGEMENTS

We are grateful to T. K. Stockard, K. Ponganis, M. Tulis, C. Champagne, K. Shiomi and B. Tribble for support and assistance in the field; J. Powell for assistance with data processing and filtering; P. Jobsis for support and advice on the development of the NIRS recorder; H. Hanish, T. Wright and M. Loughry of UFI for design and construction of the NIRS recorder; D. Alden of the SIO Hydraulics lab for design of the underwater housing unit and the SIO Machine Shop for construction of underwater housing units; and G. L. Kooyman and B. McDonald for support and manuscript review. Deposited in PMC for release after 12 months.

FOOTNOTES

This work was supported by NSF grants OPP 0229638, 0538594 and 0944220. C.L.W. was supported by an NIH Marine Biotechnology Training Program Fellowship and a UC Regents Fellowship. J.U.M. was supported by a Los Angeles Achievement Rewards for College Scientists fellowship and a Philanthropic Educational Organization Scholar Award.

Other journals from The Company of Biologists

Our most recent Editors’ choice is a pair of articles that examine rapid depth perception in hunting archerfish from Caroline P. Reinel and Stefan Schuster (I. The predictive C-starts, II. An analysis of potential cues). Archerfish have an impressive ability to predict where the prey that they topple from perches will land and now it turns out that they are able to calculate the height that a victim will fall in the first 100 ms of the descent.

Animals that live in the air have high levels of CO2 in their bodies, while aquatic species have lower CO2 levels, but what about the internal CO2 levels of animals that evolved in air but returned to the water for part of their lives? Philip Matthews and colleagues have discovered that the CO2 in the haemolymph of aquatic dragonfly nymphs is high, like that of air breathers, even though they live in water.

“So much of science is about things going wrong, failing epically, and I think you have to be able to pick yourself up and figure out what went wrong.”

Erika Eliason is an Assistant Professor at University of California, Santa Barbara, USA, where she studies ecological and evolutionary physiology. She shares how a childhood spent discovering the natural world in the woods where she grew up ignited her love for science and fieldwork, why conferences should be more family-friendly and the importance of being prepared for the unexpected.

Exposure to low temperature often results in death or injury for insects. However, repeated brief warm pulses during periods of low-temperature stress have been shown to confer some protection. Kendra J. Greenlee and colleagues provide an overview of the critical elements that underpin this phenomenon.

Meet the team

Are you going to the 2018 APS Comparative Physiology meeting in October? This meeting is supported by JEB and our Reviews Editor Charlotte Rutledge will be there - stop by and say hello! Advance registration is open until 24 September 2018.