Journal of Applied Physiology  AJP: Regulatory, Integrative and Comparative Physiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 102: 286-293, 2007. First published August 17, 2006; doi:10.1152/japplphysiol.00536.2006
8750-7587/07 $8.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
102/1/286    most recent
00536.2006v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zwemer, C. F.
Right arrow Articles by D'Alecy, L. G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Zwemer, C. F.
Right arrow Articles by D'Alecy, L. G.

Strain differences in response to acute hypoxia: CD-1 versus C57BL/6J mice

Charles F. Zwemer,1,3,* Michael Y. Song,2,* Katari A. Carello,3 and Louis G. D'Alecy3,4,5

1Department of Biology, Dickinson College, Carlisle, Pennsylvania; 2College of Literature, Science, and The Arts, University of Michigan, Ann Arbor; 3Department of Molecular and Integrative Physiology and 4Department of Surgery (Vascular), University of Michigan Medical School, Ann Arbor; and 5William Beaumont Hospital, Department of Surgery, Royal Oak, Michigan

Submitted 12 May 2006 ; accepted in final form 9 August 2006


    ABSTRACT
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Some mammals respond to hypoxia by lowering metabolic demand for oxygen and others by maximizing efficiency of oxygen usage: the former strategy is generally held to be the more effective. We describe within the same species one outbred strain (CD-1) that lowers demand and another inbred strain (C57BL/6J) that maximizes oxygen efficiency to markedly extend hypoxic tolerance. Unanesthetized adult male mice (Mus musculus, CD-1 and C57BL/6J) between 20 and 35 g were used. Sham-conditioned (SC) C57BL/6J mice survived severe hypoxia (4.5% O2, balance N2) roughly twice as long as SC CD-1 mice (median 211 and 93.5 s, respectively; P < 0.0001). Following acute hypoxic conditioning (HC), C57BL/6J mice survived subsequent hypoxia 10 times longer than HC CD-1 mice (median 2,198 and 238 s respectively; P < 0.0001). Therefore, C57BL/6J mice are both naturally more tolerant to hypoxia and show a greater increase in hypoxic tolerance in response to hypoxic conditioning. Indirect calorimetry indicates that CD-1 mice lower mass-specific oxygen consumption (V'O2 in ml O2·kg–1·min–1) and carbon dioxide production (V'CO2 in ml CO2·kg–1·min–1) in response to HC (P = 0.002 and P < 0.0001, respectively), but C57BL/6J mice maintain V'O2 and V'CO2 after HC. Respiratory exchange ratio and fluorometric assay of plasma ketones suggest that C57BL/6J mice rapidly switch to ketone metabolism, a more efficient substrate, while CD-1 mice reduce overall metabolic activity. We conclude that under severe hypoxia in mice, switching fuel, possibly to ketones, while maintaining V'O2, may confer a greater survival advantage than simply lowering demand.

hypoxic survival time; hypoxic conditioning; indirect calorimetry; respiratory exchange ratio; D-beta-hydroxybutyrate


WHEN FACED WITH HYPOXIA, mammals tend to survive by either reducing metabolic demand to levels supported by oxygen delivery or by maximizing delivery and metabolic efficiencies to sustain utilization (14, 21). Hochachka et al. (21) and Gautier (14) hold that reducing demand most commonly confers the greatest tolerance. Furthermore, Frappell et al. (13) suggest that the degree of drop in metabolic rate is inversely related to body mass such that the smaller an adult mammal species, the more significant the drop in metabolic rate on exposure to hypoxic hypoxia. Our laboratory has studied hypoxia tolerance in Swiss CD-1 mice (5, 24, 27, 28, 3034, 37), and we were the first to identify significant protection from hypoxia conferred by brief, intermittent hypoxia, now referred to as hypoxic conditioning (HC) (37). Our subsequent work linked much of this HC effect to activation of centrally active, {delta}-specific, opioid pathways (30–33). Decreases in core temperature and the seeking of cooler microenvironments in a thermocline in response to HC suggested that part of the conditioning effect may also be linked to induced decreases in metabolic demand (33). While tempting to ascribe a reduction in metabolic demand as a tolerance mechanism, linkage of hypometabolism to the mechanism of protection in our model has remained unresolved. Indeed, if the hypometabolic state is associated with reduced overall function, the overall utility for survival is limited. By contrast, a fuel shift with enhanced efficiency would be expected to maximize survival value of the acute changes.

In this study, we tested, using indirect calorimetry, the hypotheses that HC induces a rapid decrease in metabolic rate and that this hypometabolism could thus play a major role in improving hypoxic survival time. Swiss CD-1 outbred mice (CD-1) were chosen to further complement our previous work with the strain, while C57BL/6J inbred mice (C57) were chosen because of their popularity of use in studies of acute hypoxic tolerance (2, 7, 8, 20, 22, 36, 41) and induced hypometabolism (3). Additionally, C57 mice are the most common inbred strain used for deriving knockout models because of their robust tolerance of developmental manipulation and well-characterized genome (4, 6). The unexpected finding of marked strain differences in metabolic response to hypoxia opens the path to explore the genomic basis for the profound phenotypic difference within the same species.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Adult male mice (CD-1, outbred and C57, inbred) of comparable age were obtained from Charles River Laboratories, Wilmington, MA, and subjected to either sham conditioning (SC) or hypoxic conditioning (HC) and then tested to determine hypoxic survival time (HST) or plasma D-beta-hydroxybutyrate (BHB) concentrations. Additional cohorts of CD-1 and C57 were evaluated pre-, peri-, and post-hypoxic/sham conditioning for mass-specific oxygen consumption (V'O2), mass-specific carbon dioxide production (V'CO2), and respiratory exchange ratio (RER) but not tested for HST. All experiments were carried out during daylight hours (7:00 AM to 5:00 PM EST) at room temperature (~21–23°C). The University of Michigan's University Committee for Animal Care and Use approved all protocols (approval no. 8693 and no. 07124-V).

Hypoxic conditioning.   Two or three mice were placed in either of two 500-ml airtight flow-through chambers and subjected to HC or SC before HST testing or blood sampling. Both chambers were initially flushed with room air (20.93% O2) at 1.3 l/min. HC mice received four sequential exposures of 0.75, 1.5, 2, and 2.5 min of 4.6% O2-balance N2 separated by 5 min of room air including 5 min after the final hypoxic conditioning (see Fig. 1). SC mice received normoxic room air for the same duration (21.75 min) as the HC mice. A similar HC or SC protocol was followed for the cohort used in the indirect calorimetry protocols described below; however, one mouse at a time was conditioned (HC or SC) and subsequently evaluated with indirect calorimetry.


Figure 1
View larger version (10K):
[in this window]
[in a new window]

 
Fig. 1. Time course of the indirect calorimetry experiment is shown. The hypoxic conditioning (HC) period is expanded to highlight the normoxic and hypoxic periods. The end of the HC period is set as zero min. Both the preconditioning and postconditioning periods were performed under normoxia. Five min post-HC is indicated (arrow) to show where hypoxic survival time (HST) testing would have started in a separate cohort of mice.

 
Hypoxic survival time.   Groups of five (2 SC and 3 HC; or 3 SC and 2 HC) mice were transferred to five individual airtight flow-through chambers (3.7 cm ID x 10.5 cm in length) arranged in parallel for simultaneous HST testing. Hypoxic hypoxia was induced by flushing the five-chamber system for 20 s at 1.71 l/min with premixed 8.6% O2-balance N2, then by a constant flow of 4.6% O2-balance N2 at 1.71 l/min. HST, as used in this and all previous studies (11, 25, 26, 28, 3033), is the time from onset of 4.6% O2 to the cessation of spontaneous ventilation. Within 60–90 s of the onset of hypoxia, the animal falls into a state that is reported in humans as euphoric (16), loses consciousness, undergoes seizures, and then ceases spontaneous ventilation. At the cessation of spontaneous ventilation, the heart is still beating, and induction of artificial ventilation could resuscitate the animal. However, since we did not ventilate the mice, cessation of spontaneous ventilation in the unconscious mouse leads irrevocably to cardiac death and hence is a precise measurement of the time to death due to hypoxia. Hypoxic hypoxia is not asphyxia, and death by acute hypoxic exposure is humane and approved as an end point for use in this model by the university Institutional Animal Care and Use Committee that fully embraces all the articles of the recently published Principles for Establishment of Humane End Points (10).

Indirect calorimetry.   V'O2 (ml O2·kg–1·min–1), V'CO2 (ml CO2·kg–1·min–1), and RER (V'CO2/V'O2; dimensionless) were calculated using an open-flow, indirect calorimeter designed to closely match conditions used in the HC, SC, and HST phases described above. Chamber inflow rate (VI) was measured using a standard temperature and pressure (STP) dry (STPD)-corrected mass flowmeter [National Institute of Standards and Technology-traceable, GFM-17 (air, from 0 to 500 standard cubic centimeters per minute), Aalborg Instruments and Controls, Orangeburg, NY]. Fractional inlet (FIO2) and fractional outlet oxygen concentrations (FEO2) were measured separately using two high-temperature zirconia electrochemical cells (model S-3A/II oxygen analyzer, AEI Technologies, Pittsburgh, PA). Fractional outlet carbon dioxide concentration (FECO2) was measured using an infrared analyzer (CD-3A carbon dioxide analyzer, AEI Technologies) arranged in series before the FEO2 measurement. The majority flow-through system was plumbed using -in. stainless steel or copper tubing to minimize diffusional equilibration with room air while the balance was plumbed with plastic (Tygon) tubing. For normoxic exposures, compressed room air was filtered and dried before entering the mouse test chamber. For hypoxic exposures, calibrated premixed gases (Cryogenic Gasses, Detroit, MI) were filtered and dried before entering the test chamber. Test chamber outlet gases were dried and filtered before entering gas analyzers. All measured and calculated variables were recorded continuously using a commercially available analog-to-digital data logger and software (MP150 system, Biopac Instruments, Santa Barbara, CA) and laptop computer (Apple G4 Powerbook, Cupertino, CA). The processed data were analyzed in 30-s average blocks bracketing standardized time points at preconditioning, 2, 4, 5, 6, 8, 10, 20, and 30 min. All calculated calorimetry data were determined using Haldane's transformation of Fick's consumption/production equations (1, 12, 18, 19). Please see APPENDIX for a complete derivation of these formulas.

BHB assay.   Plasma concentration of BHB was determined for control and HC mice of both strains by using a fluorometry assay developed in our laboratory (29). Blood was collected by decapitation, and cervical effluent (mixed arterial and venous blood) was centrifuged at 3,000 rpm and 25°C to separate the plasma. Fifty microliters of plasma was deproteinized by using 200 µl of chilled 1 M HClO4, vortexed, placed on ice for 10 min, and finally centrifuged at 3,000 rpm and 4°C. Two-hundred microliters of supernatant was neutralized to pH 5.8–6.5 using 163 µl (determined empirically) of 0.5 M K2CO3. One-hundred microliters of neutralized sample was added to a cuvette containing a solution containing hydrazine (H-0883, Sigma Chemical, St. Louis, MO). Using the change in fluorescence and a standard curve, plasma BHB concentration was calculated (29).

Statistical analyses.   Sample sizes were sufficient to give a significance of P ≤ 0.05 and power of 0.95 or greater. Survival curves were created and analyzed by using the Kaplan and Meier product limit method and the log-rank test. Indirect calorimetry variables were analyzed by using repeated measures, two-factor (conditioning treatment and strain) ANOVA. Specifically, interaction between conditioning treatments and interaction between strains was evaluated. Plasma BHB data were analyzed using two-factor ANOVA, and interaction between strains and sampling times was evaluated. One-factor ANOVA was used to evaluate changes over sampling time within strain. ANOVA of indirect calorimetry and plasma BHB data also included Bonferroni posttests to test for differences at selected time points between treatments and strains. Where appropriate, unpaired t-tests were used to ascertain individual differences either between strains or between treatments within strains. Linear regression of HC indirect calorimetry data were analyzed by using least-squares linear regression analysis to determine if the slope was not different from zero. Statistical software consisted of Microsoft Excel and Graphpad Prism 4.0.


    RESULTS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
Strain differences in SC and HC hypoxic survival time.   Consistent with well-established strain data, CD-1 mice weighed ~15% more than comparably aged C57 mice (P < 0.01). Survival curve analysis showed that SC CD-1 mice had a median HST of 94 s, and this was significantly lower (P < 0.0001, log-rank test) than SC C57 mice with a median survival time of 211 s (Fig. 2A). HC CD-1 mice had a median HST of 238 s, and this was significantly lower (P < 0.0001) than HC C57 mice with a median survival time of 2,198 s (Fig. 2B). Survival curve analysis indicates HC prolongs survival in both CD-1 (P = 0.007) and C57 (P < 0.001) over SC control groups of the same strain (compare Fig. 2, A and B).


Figure 2
View larger version (14K):
[in this window]
[in a new window]

 
Fig. 2. Survival analysis following either sham conditioning (SC; A) or HC (B) for C57BL/6J (C57; bullet) and CD-1 ({circ}) mice (CD-1: SC n = 8, HC n = 5; C57: SC n = 11, HC n = 12). Two HC C57 mice that survived more than 1.5 h were deemed survivors and not included in statistical analysis of significance. Log-rank analysis of survival curves showed statistically significant differences (P < 0.0001) between CD-1 and C57 in A and B.

 
Strain differences during preconditioning.   The preconditioning V'O2, V'CO2, and RER were analyzed by two-factor ANOVA, and the preconditioning values were found to be not different between treatments (SC or HC); therefore the preconditioning V'O2, V'CO2, and RER were combined between treatments to obtain an average strain value (n = 16). While the RER was not different between strains, the V'O2 and V'CO2 were statistically higher in the CD-1 strain than in the C57 strain (Fig. 3).


Figure 3
View larger version (13K):
[in this window]
[in a new window]

 
Fig. 3. The preconditioning (baseline) values are pooled within strain for the groups that would receive SC or HC, producing a total of 16 mice of each strain. The pooled preconditioning mass-specific oxygen consumption and carbon dioxide production (mass-specific Formula 12O2 and Formula 12CO2) are lower in C57 mice than in the CD-1 mice; however, the respiratory exchange ratio (RER) is not statistically different.

 

Figure 4
View larger version (18K):
[in this window]
[in a new window]

 
Fig. 4. Indirect calorimetry measurement of mass-specific Formula 12O2 (A) and mass-specific Formula 12CO2 (B) in CD-1 ({circ}) and C57 mice (bullet) are shown. Values are means ± SE; n = 8 for each group. The solid lines represent the HC groups, and the dotted lines represent the SC groups. Five minutes postconditioning is indicated by a vertical dashed line that is the start time for HST testing in separate cohorts of mice with survival curves in Fig. 2. Two-factor (conditioning treatment and strain) ANOVA tested for interaction between treatments (HC vs. SC, exact P values at right of symbol key) and interaction between strains (C57 vs. CD-1, exact P values at left of symbol key).

 

Figure 5
View larger version (16K):
[in this window]
[in a new window]

 
Fig. 5. A: RER in CD-1 ({circ}, n = 8) and C57 mice (bullet, n = 8). The solid lines represent the HC groups, and the dotted lines the SC groups. Five min postconditioning is indicated by a vertical dashed line that is the start time for HST testing in separate cohorts of mice with survival curves in Fig. 2. B: plasma D-beta-hydroxybutyrate concentration ([BHB]) for C57 and CD-1 mice at control and 5 and 30 min postconditioning. Values are means ± SE with sample size from 6–10 mice in each group. In CD-1 mice, [BHB] was unchanged from control at 5 min but increased significantly at 30 min. In C57 mice, [BHB] was unchanged from control at 5 min but decreased significantly at 30 min. *Statistically significant difference (P = 0.006) between CD-1 and C57 mice [BHB] at 30 min.

 

Figure 6
View larger version (17K):
[in this window]
[in a new window]

 
Fig. 6. Differences in indirect calorimetry data at 5 min postconditioning are shown for C57 and CD-1 mice. The 5-min postconditioning values were analyzed separately in that this is the time when HST testing (Fig. 1) would have been initiated in a separate cohort of mice. Mass-specific Formula 12O2 and mass-specific Formula 12CO2 decrease with HC in the CD-1 mice but are not different with HC in the C57 mice. RER is unchanged with HC in CD-1 mice and significantly reduced in the C57 mice.

 
Strain differences during HC.   During HC (conditioning-timing illustrated in Fig. 1), V'O2 and V'CO2 appeared to decrease during each of the four sequential HC exposures before returning toward baseline with the restoration of normoxia after each exposure (data not shown). ANOVA of the final 30 s of the second, third, and fourth HC periods indicated no significant strain differences over subsequent conditioning steps for V'O2 and RER but produced a P < 0.0001 difference between strains for V'CO2, suggesting CO2 production was maintained higher in the C57 mice than the CD-1 mice (see Fig. 7). The longest HC exposure gave the most stable response, which was further analyzed for average values during the final 30 s of the fourth (final) HC period. During this final conditioning period, the CD-1 group had a significantly lower V'O2 and V'CO2 (P = 0.05 and P = 0.0065, by unpaired t-test) than the C57 mice. The V'O2 and V'CO2 of the CD-1 mice was 15.3 ± 3.0 ml O2·kg–1·min–1 and 34.9 ± 3.3 ml CO2·kg–1·min–1, respectively, and for the C57 was 25.5 ± 3.7 ml O2·kg–1·min–1 and 48.7 ± 2.8 ml CO2·kg–1·min–1, respectively.


Figure 7
View larger version (14K):
[in this window]
[in a new window]

 
Fig. 7. HC indirect calorimetry data for mass-specific Formula 12O2 (A), mass-specific Formula 12CO2 (B), and RER (C) are shown for CD-1 ({circ}) and C57 mice (bullet). The data were 30-s averages obtained at the end of each 5-min normoxia-recovery periods between hypoxic exposure during the conditioning period. Data are presented as means ± SE with n = 8 per group. A 95% confidence interval is shown for the CD-1 mice. Linear regression analysis showed a decrease in both mass-specific Formula 12O2 and mass-specific Formula 12CO2 in CD-1 mice over the conditioning period. In contrast, C57 mice mass-specific Formula 12O2 remained unchanged and mass-specific Formula 12CO2 trended downward but was not statistically significant. RER in CD-1 mice remained stable but decreased significantly in C57 mice. These results indicate the CD-1 mice decrease metabolic demand for oxygen, while C57 mice maintain metabolic demand and switch to an alternative fuel, perhaps ketones.

 
Strain differences following HC.   Analysis of the indirect calorimetry data between strains for the entire 30-min period following HC revealed that V'O2, V'CO2, and RER responded significantly differently (showed an interaction over time) in the CD-1 mice compared with the C57 mice, but with SC no such difference was detected (exact strain P values are shown at left of symbol key in Figs. 4 and 5A). For the CD-1 mice HC altered the response of V'O2 and V'CO2 but not RER. For the C57 mice, HC altered the response of V'CO2 and RER but not V'O2 (exact treatment P values are shown at right of symbol key in Figs. 4 and 5A).

Plasma BHB levels were evaluated by using two-factor (strain and time) ANOVA. Significant interaction between the strains was detected (P = 0.0015). On the basis of this significant difference, we then tested differences between (unpaired t-tests) and within strains (one-way ANOVA) at all time points. No significant differences were detected between strains at control or 5 min postconditioning, but at 30 min postconditioning, CD-1 BHB was significantly higher (unpaired t-test, P = 0.006). Within strains, BHB in CD-1 at 30 min was more than double that of control and 5-min samples (1-factor ANOVA, overall P = 0.013 and Bonferroni's posttest, P < 0.05 for each), while C57 BHB at 30 min was roughly one-half of the C57 BHB at the 5-min point (1-factor ANOVA, overall P = 0.011 and Bonferroni's posttest, P < 0.05).

Treatment effect at start of HST.   In Figs. 4 and 5, the 5-min postconditioning time point is indicated by a vertical dashed line indicating that it is the start time for HST testing. HST testing was done in separate cohorts of mice but was not done in the indirect calorimetry groups of mice. At this critical 5-min postconditioning time point, indirect calorimetry data (Fig. 6) revealed that HC in CD-1 mice produced significantly lower V'O2s and V'CO2s than SC. HC and SC values were not different in C57 mice. However, C57 mice significantly decreased RER 5 min post-HC, but CD-1 mice did not reduce their RER by 5 min post-HC.

Treatment effect during HC.   V'O2, V'CO2, and RER perihypoxic values were analyzed by using standard least-squares linear multiple regression analysis (Fig. 7). Data for each variable were taken at control, over the last 30 s of each 5-min normoxic recovery period during HC, and 5 min following cessation of HC (Fig. 1). Linear regression of perihypoxic C57 mice V'O2 (Fig. 7A) revealed a slope no different from zero, suggesting no measurable conditioning effect. The same treatment of C57 mice V'CO2 data (Fig. 7B) suggested a moderately negative slope (–0.32 ± 0.35 ml CO2·kg–1·min–2, r2 = 0.22) that was also not significantly different from zero. Over the same time period, however (Fig. 7, A and B), CD-1 mice demonstrated a significant conditioning-dependent decrease in V'O2 and V'CO2 (V'O2 slope: –1.29 ml ± 0.21 ml O2·kg–1·min–2, r2 = 0.93; V'CO2 slope: –1.33 ± 0.26 ml CO2·kg–1·min–2, r2 = 0.90). RER (Fig. 7C) decreased significantly with HC in C57 mice at a rate of –0.005 ± 0.0005 dimensionless units/min (r2 = 0.97) and was significantly different from a zero slope (P = 0.0028). CD-1 mice RER values tended to decrease at roughly one-half the rate of C57 (–0.0025 ± 0.0019 dimensionless units/min, r2 = 0.37) and did not statistically separate from a slope of zero.


    DISCUSSION
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The main findings of this study are 1) the inbred C57 strain was both naturally and hypoxia-induced more tolerant to hypoxia than the outbred CD-1 strain, and 2) since both strains displayed different metabolic changes when subjected to hypoxic conditioning, the metabolic strategy of the inbred strain, i.e., a switching toward another metabolic fuel (possibly ketones), may be more efficient than the one adopted by the outbred strain, i.e., a lowering of the metabolic demand to face the severe hypoxic challenge.

The hypothesis that CD-1 mice would decrease metabolic demand, as measured by indirect calorimetry, following HC was supported. As seen in Fig. 2, HC conferred roughly a 2.5-fold increase in median HST in CD-1 mice, and that was accompanied by a 40% and 57% decrease in V'O2 and V'CO2, respectively, at 5 min following HC, the time at which mice would be subjected to HST testing (Fig. 6). This observation is also supported by the deceleration (negative slope) of V'O2 and V'CO2 of CD-1 mice over the hypoxic conditioning period (Fig. 7, A and B). Extrapolating these data, a rough halving of aerobic demand as reflected by V'O2 could reasonably account for much of this improved HST in CD-1 mice. Using the same model system and mouse strain, Mayfield and D'Alecy (30) showed that protection afforded by HC is mediated through {delta}-selective opioid activity. The current study confirms this protection and increased HT induced by HC in CD-1 mice and extends the basic observation to a second strain of mice (C57). Additional recent work in this area suggests that hibernation and many survival-related responses in mammals may exploit a similar set of mechanisms via activation of Gi-protein receptors ultimately opening sarcolemmal and mitochondrial ATP-sensitive K channels (40). The resultant hyperpolarization of the cell membrane (in excitable tissues) putatively serves to reduce probability of depolarization and thus spares overall cell energy demand. It seems likely, then, that the Swiss CD-1 strain, when hypoxically conditioned, initiates a series of metabolic changes aimed at reducing total aerobic demand to a level supportable by delivery. This torporlike state has survival utility but only in a protected setting where minimal activity does not increase susceptibility to environmental or predatory threats.

By contrast, naive (SC) C57 mice had a median HST only 11% lower than CD-1 mice that had already increased their HST by hypoxic conditioning (Fig. 2, A and B). Furthermore, HC C57 mice survived almost 10 times longer than HC-CD-1 mice (with 2 C57 mice survivors lasting more than 1.5 h at FIO2 = 4.5%). In stark contrast to the response in CD-1 mice, C57 mice dramatically increased HC survival time without decreases in V'O2 and V'CO2 measured during conditioning (Fig. 7, A and B) or at 5 min following HC (Fig. 6) when HST testing would have started. We conclude from these data that C57 mice do not respond to HC by reducing metabolic demand; rather they must be relying on increases in delivery and or increased efficiencies in carbon fuel utilization. It is as if the C57 strain responds as though it were experiencing significantly less relative stress despite having been exposed to exactly the same FIO2.

An alternative approach to assessing HT would be to expose the different strains of mice to differing FIO2s and record at what FIO2 spontaneous ventilation ceased. This would presumably yield a different FIO2 for the CD-1 and the C57 mouse strains, reflecting their difference in overall hypoxic tolerance. This is a very different set of experiments exploring a graded hypoxic stress rather than the fixed hypoxic challenge we used and perhaps could be evaluated in future protocols.

Work on hypoxic ventilatory responses in inbred mouse strains by Tankersley et al. (42) shows that C57 mice markedly increase ventilation by increasing tidal volume (VT) and breathing frequency in response to hypoxia (FIO2 = 10% O2) and hypercapnic hypoxia (FIO2 = 10% and FICO2 = 8%), and these were among the largest responses of eight inbred strains studied. Comparative morphology of lung structure suggests C57 mice have relatively large mass-specific lung volume compared with other inbred strains (23). Using a similar hypoxia depth and time course of acute exposure per Tankersley et al. (42), Campen et al. (8) have recently shown that C57 mice decrease low-frequency (LF) power heart rate variability, while high-frequency power (HF) heart rate variability remains unchanged. Traditional interpretation of these power data suggests a diminishing sympathetic LF response along with an unchanging parasympathetic HF response and may explain the relative drop in mean arterial blood pressure and heart rate with hypoxia that these authors described (7, 8). Using an HST model similar in design to our previously established HC-HST model, Zhang et al. (45) have recently shown that HC resulted in significantly improved blood gas profiles when C57 mice were exposed to FIO2 = 7%. These authors also showed that HC significantly reduced pulmonary vascular permeability, suggesting a protection of pulmonary diffusive capacitance (45). Taken together, these strain-specific differences may work to improve hypoxic ventilatory response in C57 mice, but the exact role here remains highly speculative as measuring these differences was beyond the scope of our study. Furthermore, while an increase in oxygen delivery may have played a part in the modest enhancement in hypoxia tolerance of SC C57 mice, it is unlikely that this increase in delivery could account for a 10-fold increase in hypoxia tolerance with HC in C57 mice.

Prass et al. (35) also identified strain differences in mice; however, they focused on differences between C57 and SV129 mice. Using a 5-day exposure to hyperbaric 100% oxygen (HBO), they demonstrated that HBO induced tolerance to permanent but not transient focal ischemia in SV129 mice; however, HBO did not induce tolerance against permanent or transient focal cerebral ischemia in C57BL/6 mice. While this is relevant in that it does represent a strain difference, it is not fully comparable in that they did not, in an acute setting, use an intermittent hypoxic hypoxia to induce an increase in hypoxic survival time.

In addition to potential increases in oxygen delivery, changes in metabolic fuel selection provoked by either HC or simply in response to acute hypoxia may also have contributed to greater HST seen in C57 mice. The RER of C57 mice was depressed well below 0.6 (see proposed explanation below) in the first 2 min following HC (Fig. 5A). This was driven by a large decrease in V'CO2 and no change in V'O2 (Fig. 4, A and B). By 5 min post-HC, this suppression was gone, and V'CO2 in C57 mice was no different from the control value (Fig. 6). CD-1 mice did not demonstrate the same magnitude of change immediately following HC. Instead, CD-1 mice had decreased RER significantly by 10 min after the conditioning, in contrast to C57 mice that were increasing RER to mid-0.80 levels. We infer that C57 mice may be more agile in rate of change of fuel mobilization and selection than CD-1 mice and that ketones may be the fuel of choice in hypoxia and following HC.

We sampled whole blood for D-beta-hydroxybutyrate (BHB), arguably the most relevant of ketones, at control, 5 min, and 30 min post-HC (Fig. 5B). BHB in C57 mice tended to increase at 5 min following HC but by 30 min had decreased (by ~ 46%) significantly from 5-min values. In opposition to this response, BHB in CD-1 mice remained unchanged at 5 min but by 30 min following HC had increased by more than 1.5 times that of control. Previous studies from our laboratory (24–28) have drawn a consistent correlation between increased plasma BHB and increased hypoxic tolerance; however, regardless of the consistency and strength of the correlation, this does not establish causality.

Together, these plasma ketone and RER data suggest that C57 mice may be producing and consuming BHB rapidly during and immediately following HC. Placing this in context, Simonson and DeFronzo (39) noted that when BHB is produced from palmitate (one of the most common mammalian fatty acid chains), the calculated RER is zero since oxygen is used but no carbon dioxide is produced. If a systemic switch to palmitate-based BHB production was forced by HC, it is possible that this was happening in C57 mice at 2 min post-HC and earlier, and thus the RER was depressed well below 0.6 (Fig. 5A). By 5 min post-HC, consumption of BHB was nearly matched by production. By 10 min, BHB consumption rate slightly surpassed production rate such that between 10 and 30 min, the RER of C57 mice closely approximated the theoretical BHB RER of 0.89 (39). By 30 min, consumption surpassed production, resulting in the decrease in BHB in the plasma of C57 mice; alternatively, a production decrease could also account for the same. Along a similar line of reasoning, reducing V'O2 (and no change in RER), CD-1 mice were likely not producing BHB as early as C57 mice. By 10 min post-HC, however, they may have begun BHB production, as suggested by sub-0.8 RERs. Thus, by 30 min, net production of BHB may have exceeded consumption (RER still <0.8), and thus its concentration in the blood was elevated from control. Simple plasma concentration measurements cannot resolve this, and radiolabeled tracers would be required. It is noteworthy that the reduction in V'CO2 even during the hypoxic phases of HC in the C57 mice was statistically less compared with the changes CD-1 mice. This is most likely interpreted as a better maintenance of carbon dioxide production by oxidative metabolism in the C57 mice.

Alternatively, it is possible that BHB levels increased in CD-1 mice because of HC stress-induced systemic lipolysis. Since beta-oxidation of free fatty acids can decrease measured aerobic metabolism (oxygen consumption) this could have resulted in acetyl-CoA accumulation and thus increased total acetoacetate and BHB by mass action. C57 mice, by virtue of their relatively constant metabolic rate, may have utilized ketones at roughly the same rate ketones were made. Since we did not measure plasma free fatty acid concentrations and because BHB was massively increased by 30 min post-HC, and not earlier, it is difficult to accept this alternative as a mechanism in an animal with such a high mass-specific resting metabolic rate.

Ketones as primary fuel sources are highly efficient, and only recently has their unique role in stress and survival physiology been appreciated. In an elegant study utilizing an isolated/perfused working rat heart model, Sato et al. (38) showed that addition of BHB to the perfusion buffer increased pumping-work-to-VO2 ratio by ~30%. The authors attributed much of the ketone effect to an increase in the free energy ({Delta}G' ) of cytosolic ATP hydrolysis, the result of a BHB-dependent increase in redox energy of electron transport chain couples between sites I and II (38, 43, 44). BHB use as a fuel increased work output by roughly one-third while holding VO2 unchanged or slightly decreased (34). Veech et al. (43) suggests that beyond the increased efficiency of work, the greater {Delta}G' of ATP hydrolysis at the Na+-K+ ATPase pump widens the potential energy of extra/intracellular ionic gradients and thus hyperpolarizes highly excitable cell membranes and reduces their unwanted discharge. Whether BHB-mediated increases in {Delta}G' of ATP hydrolysis is at work in our model remains speculative. Additionally, other possible mechanisms have recently been reviewed and include the role of activation of several different signaling pathways and modulation of potential second messengers that could contribute to the preconditioning response (15).

Our analysis of the data also suggests that the rate change of RER [|d(RER)/dt| in dimensionless units/min2] during HC or immediately following HC (first 5 min) may be mechanistically linked to two strategies we observed between the strains, with a greater relative |d(RER)/dt| indicating greater metabolic agility and perhaps enhanced hypoxic tolerance. It holds, then, that measured |d(RER)/dt| on hypoxic challenge could indicate which strategy of hypoxic tolerance is to be used by the animal, and it may be a relatively sensitive metric predictive of qualitative HST.

We conclude that CD-1 mice improve HST following HC primarily by lowering metabolic demand as measured by indirect calorimetry. Quite differently, C57 mice did not decrease metabolic demand during HST testing to match oxygen delivery but instead appear to switch metabolic fuel, perhaps to BHB, and thereby gain putative efficiencies not possible with other substrates. This fuel change, combined with a putative increase in oxygen delivery as suggested by Zhang et al. (45), may explain the dramatic improvement in HST seen in HC C57 mice (41) compared with CD-1 mice. This profound phenotypic strain difference within the same species presents an excellent opportunity of exploring the genotypic basis for hypoxic tolerance. In addition, the demonstrated capacity of C57 mice for agile switching of fuel choice suggests enhanced ketone metabolism should be evaluated as a potential therapeutic target for management of acute low oxygen states.


    APPENDIX
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
The following are the abbreviations used in derivations of the V'O2, V'CO2, and RER equations found below:

Glossary

VO2
volumetric rate (STP) of oxygen consumption (ml O2/min)

V'O2
mass-specific volumetric rate (STP) of oxygen consumption (ml O2·kg–1·min–1)

VCO2
volumetric rate (STP) of carbon dioxide production (ml CO2/min)

V'CO2
mass-specific volumetric rate (STP) of carbon dioxide production (ml CO2·kg–1·min–1)

VI
volumetric rate (STP) of chamber inlet airflow (ml air/min)

VE
volumetric rate (STP) of chamber outlet airflow (ml air/min)

FIO2
chamber inlet fraction of oxygen (decimal %)

FEO2
chamber outlet fraction of oxygen (decimal %)

FICO2
chamber inlet fraction of carbon dioxide (decimal %)

FECO2
chamber outlet fraction of carbon dioxide (decimal %)

FIN2
chamber inlet fraction of nitrogen (decimal %)

FEN2
chamber outlet fraction of nitrogen (decimal %)

FO2
fraction of oxygen (decimal %)

FCO2
fraction of carbon dioxide (decimal %)

FN2
fraction of nitrogen (decimal %)

RER
respiratory exchange ratio (dimensionless)

Haldane's transformation of Fick's equation for oxygen consumption proceeds as follows, where

Formula 1(1)

Formula 2(2)
and

Formula 3(3)
Thus if

Formula 4(4)
and

Formula 5(5)
and

Formula 6(6)
then by substitution

Formula 7(7)
So,

Formula 8(8)
or in mass-specific terms and using distributive property for Formula 8I:

Formula 9(9)
For mass-specific carbon dioxide production,

Formula 10(10)
and, making the same corrections [Haldane transformation (17–19)] as used above for oxygen and dividing through by mass, then:

Formula 11(11)
Respiratory exchange ratio was calculated as the following:

Formula 12(12)


    GRANTS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
C. F. Zwemer was supported in part by a sabbatical supplement grant from the Dickinson College Dean's Office. M. Y. Song was supported by the University of Michigan Undergraduate Research Opportunities Program.


    ACKNOWLEDGMENTS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 
We thank Steven E. Whitesall and Mary C. Lloyd for expert and excellent advice and assistance throughout this project. We also thank Prof. Henry D. Prange for help with revisions of this manuscript.


    FOOTNOTES
 

Address for reprint requests and other correspondence: L. G. D'Alecy, 7744 Medical Sciences Bldg. II, Dept. of Molecular and Integrative Physiology, Univ. of Michigan Medical School, Ann Arbor, MI 48109-0622 (e-mail: lgdalecy{at}umich.edu)

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

* C. F. Zwemer and M. Y. Song contributed equally to this study. Back


    REFERENCES
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 GRANTS
 ACKNOWLEDGMENTS
 REFERENCES
 

  1. Bassett DR Jr, Howley ET, Thompson DL, King GA, Strath SJ, McLaughlin JE, Parr BB. Validity of inspiratory and expiratory methods of measuring gas exchange with a computerized system. J Appl Physiol 91: 218–224, 2001.[Abstract/Free Full Text]
  2. Beckmann N. High resolution magnetic resonance angiography non-invasively reveals mouse strain differences in the cerebrovascular anatomy in vivo. Magn Reson Med 44: 252–258, 2000.[CrossRef][ISI][Medline]
  3. Blackstone E, Morrison M, Roth MB. H2S induces a suspended animation-like state in mice. Science 308: 518, 2005.[Abstract/Free Full Text]
  4. Blake JA, Richardson JE, Bult CJ, Kadin JA, Eppig JT. The Mouse Genome Database (MGD): the model organism database for the laboratory mouse. Nucleic Acids Res 30: 113–115, 2002.[Abstract/Free Full Text]
  5. Bofetiado DM, Mayfield KP, D'Alecy LG. Alkaloid delta agonist BW373U86 increases hypoxic tolerance. Anesth Analg 82: 1237–1241, 1996.[Abstract]
  6. Bult CJ, Blake JA, Richardson JE, Kadin JA, Eppig JT, Baldarelli RM, Barsanti K, Baya M, Beal JS, Boddy WJ, Bradt DW, Burkart DL, Butler NE, Campbell J, Corey R, Corbani LE, Cousins S, Dene H, Drabkin HJ, Frazer K, Garippa DM, Glass LH, Goldsmith CW, Grant PL, King BL, Lennon-Pierce M, Lewis J, Lu I, Lutz CM, Maltais LJ, McKenzie LM, Miers D, Modrusan D, Ni L, Ormsby JE, Qi D, Ramachandran S, Reddy TB, Reed DJ, Sinclair R, Shaw DR, Smith CL, Szauter P, Taylor B, Vanden Borre P, Walker M, Washburn L, Witham I, Winslow J, Zhu Y. The Mouse Genome Database (MGD): integrating biology with the genome. Nucleic Acids Res 32: D476–D481, 2004.[Abstract/Free Full Text]
  7. Campen MJ, Shimoda LA, O'Donnell CP. Acute and chronic cardiovascular effects of intermittent hypoxia in C57BL/6J mice. J Appl Physiol 99: 2028–2035, 2005.[Abstract/Free Full Text]
  8. Campen MJ, Tagaito Y, Jenkins TP, Balbir A, O'Donnell CP. Heart rate variability responses to hypoxic and hypercapnic exposures in different mouse strains. J Appl Physiol 99: 807–813, 2005.[Abstract/Free Full Text]
  9. Dalton J. Experimental enquiry into the proportion of the several gases or elastic fluids, constituting the atmosphere. Memoirs Lit Philosoph Soc Manchester 1: 244–258, 1805.
  10. Demers G, Griffin G, De Vroey G, Haywood JR, Zurlo J, Bedard M. Animal research. Harmonization of animal care and use guidance. Science 312: 700–701, 2006.[Abstract/Free Full Text]
  11. Eiger SM, Kirsch JR, D'Alecy LG. Hypoxic tolerance enhanced by beta-hydroxybutyrate-glucagon in the mouse. Stroke 11: 513–517, 1980.[Abstract/Free Full Text]
  12. Fick A. On the measurement of blood quantity in the ventricles of the heart. In: Proceedings of the Physiological Medical Society of Wurzburg. Wurzburg, Germany, 1870.
  13. Frappell P, Lanthier C, Baudinette RV, Mortola JP. Metabolism and ventilation in acute hypoxia: a comparative analysis in small mammalian species. Am J Physiol Regul Integr Comp Physiol 262: R1040–R1046, 1992.[Abstract/Free Full Text]
  14. Gautier H. Interactions among metabolic rate, hypoxia, and control of breathing. J Appl Physiol 81: 521–527, 1996.[Abstract/Free Full Text]
  15. Gidday JM. Cerebral preconditioning and ischaemic tolerance. Nat Rev Neurosci 7: 437–448, 2006.[CrossRef][ISI][Medline]
  16. Gradwell RHD. Aviation Medicine. Oxford, UK: Butterworth Heinemann, 1999.
  17. Haldane JS. Respiration. New Haven, CT: Yale Univ. Press, 1922.
  18. Haldane JS, Graham JI. Methods of Air Analysis. London, UK: Griffin, 1912.
  19. Haldane JS, Priestley JG. Respiration. New Haven, CT: Yale Univ. Press, 1935.
  20. Han F, Subramanian S, Dick TE, Dreshaj IA, Strohl KP. Ventilatory behavior after hypoxia in C57BL/6J and A/J mice. J Appl Physiol 91: 1962–1970, 2001.[Abstract/Free Full Text]
  21. Hochachka PW, Buck LT, Doll CJ, Land SC. Unifying theory of hypoxia tolerance: molecular/metabolic defense and rescue mechanisms for surviving oxygen lack. Proc Natl Acad Sci USA 93: 9493–9498, 1996.[Abstract/Free Full Text]
  22. Iwase M, Izumizaki M, Kanamaru M, Homma I. Effects of hyperthermia on ventilation and metabolism during hypoxia in conscious mice. Jpn J Physiol 54: 53–59, 2004.[CrossRef][ISI][Medline]
  23. Kida K, Fujino Y, Thurlbeck WM. A comparison of lung structure in male DBA and C57 black mice and their F1 offspring. Am Rev Respir Dis 139: 1238–1243, 1989.[ISI][Medline]
  24. Kirsch JR, D'Alecy LG. Effect of altered availability of energy-yielding substrates upon survival from hypoxia in mice. Stroke 10: 288–291, 1979.[Free Full Text]
  25. Kirsch JR, D'Alecy LG. Glucagon stimulates ketone utilization by rat brain slices. Stroke 15: 324–328, 1984.[Abstract/Free Full Text]
  26. Kirsch JR, D'Alecy LG. Hypoxia induced preferential ketone utilization by rat brain slices. Stroke 15: 319–323, 1984.[Abstract/Free Full Text]
  27. Kirsch JR, D'Alecy LG. Role of tissue lactate and substrate availability in 1,3-butanediol-enhanced hypoxic survival in the mouse. Stroke 14: 971–976, 1983.[Abstract/Free Full Text]
  28. Kirsch JR, D'Alecy LG, Mongroo PB. Butanediol induced ketosis increases tolerance to hypoxia in the mouse. Stroke 11: 506–513, 1980.[Abstract/Free Full Text]
  29. Lloyd MC, D'Alecy LG. Assay for the determination of D-beta-hydroxybutyrate using the TKO 100 Mini-Fluorometer. In: Technical Bulletin 139. San Francisco, CA: Hoefer Scientific Instruments, 1992.
  30. Mayfield KP, D'Alecy LG. Delta-1 opioid agonist acutely increases hypoxic tolerance. J Pharmacol Exp Ther 268: 683–688, 1994.[Abstract/Free Full Text]
  31. Mayfield KP, D'Alecy LG. Delta-1 opioid receptor dependence of acute hypoxic adaptation. J Pharmacol Exp Ther 268: 74–77, 1994.[Abstract/Free Full Text]
  32. Mayfield KP, D'Alecy LG. Role of endogenous opioid peptides in the acute adaptation to hypoxia. Brain Res 582: 226–231, 1992.[CrossRef][ISI][Medline]
  33. Mayfield KP, Hong EJ, Carney KM, D'Alecy LG. Potential adaptations to acute hypoxia: Hct, stress proteins, and set point for temperature regulation. Am J Physiol Regul Integr Comp Physiol 266: R1615–R1622, 1994.[Abstract/Free Full Text]
  34. Moursi MM, Luyckx BA, D'Alecy LG. The role of ethanol in diluents of drugs that protect mice from hypoxia. Stroke 14: 791–796, 1983.[Abstract/Free Full Text]
  35. Prass K, Wiegand F, Schumann P, Ahrens M, Kapinya K, Harms C, Liao W, Trendelenburg G, Gertz K, Moskowitz MA, Knapp F, Victorov IV, Megow D, Dirnagl U. Hyperbaric oxygenation induced tolerance against focal cerebral ischemia in mice is strain dependent. Brain Res 871: 146–150, 2000.[CrossRef][ISI][Medline]
  36. Price ER, Han F, Dick TE, Strohl KP. 7-Nitroindazole and posthypoxic ventilatory behavior in the A/J and C57BL/6J mouse strains. J Appl Physiol 95: 1097–1104, 2003.[Abstract/Free Full Text]
  37. Rising CL, D'Alecy LG. Hypoxia-induced increases in hypoxic tolerance augmented by beta-hydroxybutyrate in mice. Stroke 20: 1219–1225, 1989.[Abstract/Free Full Text]
  38. Sato K, Kashiwaya Y, Keon CA, Tsuchiya N, King MT, Radda GK, Chance B, Clarke K, Veech RL. Insulin, ketone bodies, and mitochondrial energy transduction. FASEB J 9: 651–658, 1995.[Abstract]
  39. Simonson DC, DeFronzo RA. Indirect calorimetry: methodological and interpretative problems. Am J Physiol Endocrinol Metab 258: E399–E412, 1990.[Abstract/Free Full Text]
  40. Sonneborn JS, Gottsch H, Cubin E, Oeltgen P, Thomas P. Alternative strategy for stress tolerance: opioids. J Gerontol A Biol Sci Med Sci 59: 433–440, 2004.
  41. Subramanian S, Han F, Erokwu BO, Dick TE, Strohl KP. Do genetic factors influence the Dejours phenomenon? Adv Exp Med Biol 499: 209–214, 2001.[ISI][Medline]
  42. Tankersley CG, Fitzgerald RS, Kleeberger SR. Differential control of ventilation among inbred strains of mice. Am J Physiol Regul Integr Comp Physiol 267: R1371–R1377, 1994.[Abstract/Free Full Text]
  43. Veech RL. The therapeutic implications of ketone bodies: the effects of ketone bodies in pathological conditions: ketosis, ketogenic diet, redox states, insulin resistance, and mitochondrial metabolism. Prostaglandins Leukot Essent Fatty Acids 70: 309–319, 2004.[CrossRef][ISI][Medline]
  44. Veech RL, Chance B, Kashiwaya Y, Lardy HA, Cahill GF Jr. Ketone bodies, potential therapeutic uses. IUBMB Life 51: 241–247, 2001.[ISI][Medline]
  45. Zhang SX, Miller JJ, Gozal D, Wang Y. Whole-body hypoxic preconditioning protects mice against acute hypoxia by improving lung function. J Appl Physiol 96: 392–397, 2004.[Abstract/Free Full Text]




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
102/1/286    most recent
00536.2006v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar