Journal of Applied Physiology Watch the video to learn how APS reaches out to developing nations.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 91: 2720-2729, 2001;
8750-7587/01 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
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 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 HighWire
Right arrow Citing Articles via Web of Science (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Fahlman, A.
Right arrow Articles by Kayar, S. R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Fahlman, A.
Right arrow Articles by Kayar, S. R.
Vol. 91, Issue 6, 2720-2729, December 2001

On the likelihood of decompression sickness during H2 biochemical decompression in pigs

Andreas Fahlman1, Peter Tikuisis2, Jeffrey F. Himm1, Paul K. Weathersby1, and Susan R. Kayar1

1 Environmental Physiology Department, Naval Medical Research Center, Silver Spring, Maryland 20910-7500; and 2 Defence and Civil Institute of Environmental Medicine, Toronto, Ontario, Canada M3M 3B9


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

A probabilistic model was used to predict decompression sickness (DCS) outcome in pigs during exposures to hyperbaric H2 to quantify the effects of H2 biochemical decompression, a process in which metabolism of H2 by intestinal microbes facilitates decompression. The data set included 109 exposures to 22-26 atm, ca. 88% H2, 9% He, 2% O2, 1% N2, for 0.5-24 h. Single exponential kinetics described the tissue partial pressures (Ptis) of H2 and He at time t: Ptis = int  (Pamb - Ptis) · tau -1 dt, where Pamb is ambient pressure and tau  is a time constant. The probability of DCS [P(DCS)] was predicted from the risk function: P(DCS) = 1 - e-r, where r = int  (PtisH2 + PtisHe - Thr - Pamb) · Pamb-1 dt, and Thr is a threshold parameter. Inclusion of a parameter (A) to estimate the effect of H2 metabolism on P(DCS): PtisH2 int  (Pamb - A - PtisH2) · tau -1 dt, significantly improved the prediction of P(DCS). Thus lower P(DCS) was predicted by microbial H2 metabolism during H2 biochemical decompression.

probabilistic modeling; Sus scrofa; hydrogen diving; H2 metabolism; Methanobrevibacter smithii


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

MODELING OF DECOMPRESSION SICKNESS (DCS) risk has been impeded by the inability to identify correlated physiological variables. Some studies have tried to find a correlation between DCS risk and variables such as body temperature, body weight, exercise, gender, adiposity, age, serum cholesterol, sensitivity to complement activation, Doppler bubble grades, and patent foramen ovale (4, 12, 16, 23, 28). However, where some studies have found a correlation, others refute those results (5, 8, 16). The only physiological variable that has been undisputedly correlated with DCS risk in rats is body weight (20). Because reliable physiological correlates are lacking, researchers have used a variety of models based solely on the physical history of the compression and decompression sequence to find variables that can predict the probability of DCS (26, 31-34).

The DCS risk assessment used in this study builds on previously published models used in DCS research (26, 31, 33, 34). The goal is to estimate the beneficial effects on DCS risk of the active removal of tissue H2 by injecting H2-metabolizing microbes into the intestines of pigs during simulated H2 dives and to suggest a physiological mechanism for the process called H2 biochemical decompression (19). The metabolism of H2 in the intestine is readily followed by measuring the release of CH4, the metabolic end product of the microbial metabolism (21)
4H<SUB>2</SUB><IT>+</IT>CO<SUB>2</SUB><IT> → </IT>2H<SUB>2</SUB>O<IT>+</IT>CH<SUB>4</SUB> (1)
The model presented here differs from earlier models of DCS (26, 31, 33, 34) in that a parameter for the microbial metabolism of H2 is included. In constructing this model, the microbial metabolism of H2 was considered to have a direct physiological effect by influencing the gas kinetics. The measure of H2 metabolism was based on either the total microbial activity injected into the animals (Inj), or as the CH4 release rate (<A><AC>V</AC><AC>˙</AC></A><SC>ch</SC><SUB>4</SUB>) from individual animals, assuming that there is a direct correlation between the H2 metabolized inside the intestines and the release of CH4 from the intact animal (Eq. 1). This model was used to predict the advantage of biochemical decompression for hyperbaric H2 exposures in general and thus offers a predictive advantage over the descriptive logistic regression model of biochemical decompression that was derived earlier (10, 18).

As a result, this study used a measured physiological variable, namely the rate of removal of H2 by the metabolism of intestinal microbes, together with the physical history of the hyperbaric exposure to predict the DCS incidence after a hyperbaric exposure to H2.


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Data Set and Case Descriptions

The data used in this study were taken from animal experiments previously reviewed and approved and have been reported in detail elsewhere (11, 18). A brief explanation of the experimental procedure will be provided here.

Pigs (Sus scrofa, 17-23 kg) were used for all experiments. The pigs were housed before experiments in an accredited animal care facility and had ad libitum access to water. The pigs were fed once daily with laboratory animal chow (Harlan Teklad, Madison, WI; 2% by body wt). All procedures were approved by an Animal Care and Use Committee. The experiments reported here were conducted according to the principles presented in the Guide for the Care and Use of Laboratory Animals (National Research Council, 1996).

The data set contains 109 well-documented hyperbaric exposures, with 53 severe DCS cases (Table 1). All animals were juvenile male Yorkshire pigs, either castrated (n = 98) or noncastrated (n = 11). There was no difference in DCS incidence between the castrated and noncastrated animals (logistic regression analysis, P > 0.3). Consequently, these two groups were pooled for all subsequent analyses.

                              
View this table:
[in this window]
[in a new window]
 
Table 1.   Experimental protocols

The animals were divided into three groups: untreated control (UC, n = 69, mean body wt 19.6 ± 1.4 kg), surgical control (SC, n = 11, mean body wt 19.6 ± 1.6 kg), and treated (n = 29, mean body wt 19.5 ± 1.3 kg). Treated animals had H2-metabolizing microbes injected into the large intestine (mean 0.92 ± 0.49 mmol CH4/min, range 0.20-2.20 mmol CH4/min), a procedure requiring major surgery. This surgical procedure and the culturing of Methanobrevibacter smithii are similar to prior work in rats (19) and have been described elsewhere (18). In the SC group, animals received intestinal injections of saline bubbled with CO2 to deoxygenate it. The SC group has been described in detail along with the surgical procedure (18). Because DCS incidence was indistinguishable between the SC and UC groups (P > 0.29, Fisher exact test), the two groups were pooled and will be referred to as control animals (C = UC + SC, n = 80).

The experiments were carried out over a period of 30 mo (May 1997-Dec 1999) and include a variety of different pressurization and depressurization sequences (Table 1). The hyperbaric H2 exposures were performed in a dry chamber (WSF Industries, Buffalo, NY) of 5,600-liter volume at 1 atm. Each animal was subjected to only one hyperbaric exposure. The chamber was controlled by a computer that stored the ambient pressure (Pamb, atm), chamber temperature (°C), O2 concentration (%), and elapsed time (min). A gas chromatograph (Hewlett-Packard 5890A, Series II, Wilmington, DE) measured the chamber gases O2, H2, He, N2, and CH4 every 12 min.

Each hyperbaric exposure commenced with a pressurization of the chamber with He to 11 atm, followed by a flush of the chamber with H2 until the H2 concentration was over 60%. The initial pressurization with He was necessary as a safety measure to prevent an explosive mixture of H2 and O2 in the chamber, as described in detail elsewhere (18-20). Pressurization then continued with H2, with addition of O2 as needed to maintain normoxia (0.2-0.4 atm O2). Maximum pressure of the exposures ranged from 22 to 26 atm (absolute pressure; 2.23-2.63 MPa; 700-825 feet of seawater pressure equivalent; Table 1). Final gas composition in the chamber at maximal pressure was roughly 88% H2, 9% He, 2% O2, and 1% N2 for exposures lasting 3 h, with more H2 and less He and N2 present in the chamber over time. Chamber concentrations of CH4 ranged from <1 ppm to 8 ppm after 3 h, and up to 40 ppm after 24 h.

The start of the decompression (time 0) was set as the time that animals were last definitely free of DCS signs (T1). The chamber was decompressed at 0.45-1.8 atm/min to 11 atm while each animal was observed closely. Decompression could not continue past 11 atm because of the need for a normoxic environment and safe handling of H2 and O2 mixtures as described earlier (18-20). On arrival at 11 atm, the animal was made to walk on a treadmill inside the chamber at 5-min intervals for up to 1 h (67-90 min, Tend) or until the animal was declared to have DCS (T2) (18).

Most animals displayed violaceous macular lesions with or without pruritus. On the basis of previous descriptions in the literature, these lesions resembled skin DCS (6, 9). However, skin DCS alone was not considered a DCS case for the purpose of this study. Severe DCS symptoms were of two types: cardiopulmonary or neurological (3, 6, 9). Severe symptoms of DCS included walking difficulties, fore- and/or hindlimb paralysis, falling, seizures, labored breathing, and convulsions. Every diagnosis of DCS was determined by a consensus of at least three observers who discussed their diagnosis as the events occurred. The observers were not blind to the treatments or dive profiles, resulting in the possibility of biased diagnosis. However, the rigorous compression and decompression sequences used throughout this study caused signs of DCS that were in most cases severe enough to leave no ambiguity. In the rare event that the observers did not agree, a detailed description or a viewing of a video recording of the animal was given to a neurologist or diving medical officer to evaluate and pass final judgment.

DCS Risk Assessment Modeling

A model of the probability of DCS [P(DCS)] was defined, using the method of maximum likelihood to search for the best fitting parameters (7, 31). Experience has shown that elevated pressure, longer exposure to elevated pressure, and increasing decompression rate all increase the risk of DCS (14, 34). However, the occurrence is seldom either a certainty or zero for any hyperbaric exposure (34). Therefore, researchers have used probabilistic models to predict the P(DCS) in dives with varying compression and decompression profiles (31, 33, 34). In previous models, the outcome data for each hyperbaric exposure were used to estimate the model parameters. As a result, the parameters are a mathematical composition and have only limited physiological value (2).

The probabilistic models used previously, unlike most logistic regressions, were constructed to be well behaved in the sense that the P(DCS) was predicted to be 0 when no pressure reduction occurred and increased with increasing pressure reduction (34).

The probability of having DCS symptoms at a time t, after a hyperbaric exposure, is defined as
P(DCS)<IT>=</IT>1.0<IT>−</IT>exp<FENCE>−<LIM><OP>∫</OP><LL>0</LL><UL><IT>T</IT></UL></LIM><IT> r </IT>d<IT>t</IT></FENCE> (2)
whereas freedom of symptoms until time t is defined as
P(no DCS)<IT>=</IT>1.0<IT>−P</IT>(DCS)<IT>=</IT>exp<FENCE>− <LIM><OP>∫</OP><LL>0</LL><UL><IT>T</IT></UL></LIM><IT> r </IT>d<IT>t</IT></FENCE> (3)
where r is the instantaneous risk (7), and r >=  0. The r depends on the theory used to describe the mechanism of DCS and can be one of several measures integrated over the dive and postdive period.

The probability of developing DCS for a particular hyperbaric exposure is defined as the integrated r over the exposure period through the end of the postdecompression observation period (Tend). Consequently, if DCS was observed, the following formula was used
P(DCS)<IT>=</IT><FENCE>exp<FENCE>−<LIM><OP>∫</OP><LL>0</LL><UL><IT>T</IT><SUB>1</SUB></UL></LIM><IT> r </IT>d<IT>t</IT></FENCE></FENCE><FENCE>1.0<IT>−</IT>exp<FENCE>−<LIM><OP>∫</OP><LL><IT>T</IT><SUB>1</SUB></LL><UL><IT>T</IT><SUB>2</SUB></UL></LIM><IT> r </IT>d<IT>t</IT></FENCE></FENCE> (4)
This is the product of the probability of DCS not occurring in the interval [0-T1] and the probability of DCS occurring in the interval [T1-T2] (26, 34). T1 is defined as the last time when the diver was definitely free of DCS symptoms, and T2 the time the diver was declared to have definite signs of DCS (26, 34). For an animal with no DCS, Eq. 2 is used to estimate the P(DCS) with the upper limit for integration at the end of the postdecompression observation period (Tend). For an animal displaying DCS, Eq. 4 is used and the integration is performed until the time the animal definitely shows signs of DCS (T2).

The choice of T1 can be made in various ways, as detailed elsewhere (34). Because of subjective aspects inherent in diagnosing DCS, T1 is set to the start of the decompression for this study. This approach is conservative and may result in some loss of temporal information. However, because T1 is difficult to determine in some cases, this approach seems warranted.

Model Variations

Two models were tested, both describing r from a single tissue. One model used the speculation that the maximum r develops rapidly after a decompression step after which it decreases immediately as Pamb > Ptis (Immediate model, Eq. 5, Fig. 1A). The second model (Delayed model, Eq. 6, Fig. 1B) used the hypothesis that r increases more slowly to a maximum value, after which it decreases slowly.


View larger version (18K):
[in this window]
[in a new window]
 
Fig. 1.   Model behavior showing the ambient pressure (Pamb), estimated tissue tension for H2 and He (Ptis), the probability of decompression sickness [DCS; P(DCS)], the instantaneous risk (r1, Eq. 5) for the Immediate (A) model, and the relative supersaturation (r2, Eq. 6) for the Delayed (B) model using a time constant (tau ) and a scale factor (G) for a sample dive to 24 atm. The Pamb, the r1, and the r2 for both models (C) are displayed to show the delayed maximum for the r2 compared with r1. Constant pressure was maintained for 3 h followed by decompression at 0.9 atm/min to 11 atm. The Immediate model (A) used tau  = 0.80 and G = 2.23 whereas the Delayed model (B) used tau  = 0.75 and G = 0.14. The units for r are inverse of time; the scale for r is arbitrary for convenience.

It is widely accepted that the presence of gas bubbles in tissues leads to DCS (14). If one believes that the presence of bubbles triggers the symptoms, it may be most appropriate to use a risk function in which the greatest r occurs immediately after a decompression step (34). Alternatively, if one prefers the theory that DCS develops after the bubbles have grown to a certain size (28) or that they trigger an immune or hematological response (17, 24, 35), r may rise slowly to a maximum and finally decrease slowly. Unrelated to the theory used, when r is zero there should be no occurrence of DCS, and as r increases so should the DCS incidence (34).

Immediate model. The instantaneous risk (r1), is defined as the relative difference between the tissue tension (Ptis, atm) and the absolute Pamb (atm) above a threshold (Thr) (26, 33, 34)
r<SUB>1</SUB>=G<SUB>1</SUB>·(Ptis<IT>−</IT>Pamb<IT>−</IT>Thr)<IT>·</IT>Pamb<SUP>−1</SUP> (5)
where Ptis refers to the sum of the tissue tensions for H2 (PtisH2, atm) and He (PtisHe, atm) and G1 is a scaling factor (min-1) to be determined from the fitting procedure. The inclusion of a threshold parameter has been shown to improve the fit for some human data sets (33), whereas this parameter has not been very successful in describing others (26). The contributions of O2, CO2, and water vapor to Pamb (1) were ignored, as in some other DCS modeling efforts (29). In this model, r1 is constrained to be >= 0, and, accordingly, r1 will be set to 0 at any time Pamb > Ptis (33, 34). The maximum value of r1 is reached at arrival at 11 atm, after which it decreases quickly as Ptis approaches Pamb (Fig. 1A).

Delayed model. The second model has the same structure as the Immediate model, but here the relative supersaturation is integrated over the exposure time (34)
r<SUB>2</SUB>=G<SUB>2</SUB>·<LIM><OP>∫</OP><LL>0</LL><UL>T</UL></LIM> (Ptis<IT>−</IT>Pamb<IT>−</IT>Thr)<IT>·</IT>Pamb<SUP>−1</SUP> d<IT>t</IT> (6)
For this model, the relative supersaturation is integrated from the first occurrence that Ptis > Pamb, followed by integration of both positive and negative values of supersaturation. However, the risk cannot be negative, so r2 is constrained to be >= 0. The value of r2 increases more slowly compared with r1 and reaches its maximum value at the time when Pamb > Ptis (Fig. 1B). As can be seen in Fig. 1C, the r2 maximum value occurs slightly later in the decompression compared with r1. This leads to a delayed effect on the accumulation of P(DCS) for the Delayed model (Fig. 1B) compared with the Immediate model (Fig. 1A).

Tissue Inert Gas Tension

A single exponential gas kinetics model was used to describe the tissue tensions of the inert gases, assuming a single compartment. In other words, the animal was considered to be composed of a single tissue, with a single perfusion rate and tissue gas solubility. A more complex model is often used to describe real tissues, in which the diver is assumed to be composed of several tissues with varying perfusion rates and gas solubilities (15). For this study, the models assumed that the gas uptake and elimination followed symmetrical exponential kinetics (34).

The differential equation used to describe the inert gas tissue tension was as follows
<FR><NU>dPtis</NU><DE>d<IT>t</IT></DE></FR><IT>=</IT><FR><NU>Pblood</NU><DE><IT>&tgr;</IT></DE></FR><IT>−</IT><FR><NU>Ptis</NU><DE><IT>&tgr;</IT></DE></FR> (7)
where Ptis is the tissue tension of the inert gas (atm), Pblood is the arterial blood tension of the inert gas (atm), tau  is the time constant (min) to be determined from the data, and t is time (min). The time constant determines the flux of gases in and out of the tissues. The change in Pblood during compression and decompression is assumed to be instantaneous and therefore equal to the alveolar partial pressure for that inert gas at all times. The alveolar partial pressure in turn is assumed to be equal to the Pamb for that gas.

The effect of a change in Pblood on the Ptis is computationally extensive and has been described in detail elsewhere (26, 33). The inert gas flux during the hyperbaric experiment was calculated by dividing the compression and decompression sequence into pressure-time ramps (33, 34). For the present model, it is assumed that tau  is the same for He and H2. Hence, three parameters need to be determined: tau , G, and Thr.

Effect of Microbial H2 Metabolism on PtisH2

Animals injected with a H2-metabolizing microbe, M. smithii, into the intestines had a significantly lower incidence of DCS compared with control animals (10, 11, 17, 18). It has also been shown that pigs have a native intestinal flora of H2-metabolizing microbes that, if sufficiently active, provides some protection against DCS (11). It is unclear how the reduction of H2 in the cecum and large intestine affects the presence of H2 elsewhere in the body. Figure 2 shows the current working hypothesis on how the conversion of H2 into H2O and CH4 in the intestine (Eq. 1) creates a sink that ultimately reduces the arterial (PbloodH2) and tissue tension (PtisH2) throughout the body. The uptake or removal of H2 in the various regions of the body is dependent on the local perfusion rate and the difference in the arterial and tissue tensions of H2. Therefore, the equation used to describe the H2 tissue tension including H2-metabolism was as follows
<FR><NU>dPtis<SUB>H<SUB>2</SUB></SUB></NU><DE>d<IT>t</IT></DE></FR><IT>=</IT><FR><NU>Pblood<SUB>H<SUB>2</SUB></SUB><IT>−</IT>BUG<IT>∗A</IT></NU><DE><IT>&tgr;</IT></DE></FR><IT>−</IT><FR><NU>Ptis<SUB>H<SUB>2</SUB></SUB></NU><DE><IT>&tgr;</IT></DE></FR> (8)
where BUG is the rate of H2 removal as determined by either the measured CH4 release rate (VCH4 , mmol CH4/min ) or by the total microbial activity injected into the intestines (Inj, mmol CH4/min ). The AVCH4 and AInj [atm · min · (mmol CH4)-1] are the respective parameters to be determined. This equation assumes that the microbial metabolism of H2 is constant throughout the hyperbaric exposure. This simplification appears justified, because our data during 3 h at constant pressure did not show any temporal changes in the metabolism of H2 (18). Consequently, the removal of H2 by microbial metabolism in addition to normal gas kinetics results in up to four parameters that need to be determined for each model: tau , G, Thr, and AVCH4 or AInj.


View larger version (34K):
[in this window]
[in a new window]
 
Fig. 2.   Proposed physiological mechanism for H2 biochemical decompression. For each figure, the number in parentheses represents the net result of the process of gas uptake (compression) or elimination (decompression) by the body at the end of a circulatory pass from the gas exchange surface throughout the body and returning to the gas exchange surface. Units in these figures are hypothetical and are used only to describe the current working hypothesis.   During compression (A, top) in an untreated animal exposed to a hyperbaric environment of 100 units of H2 pressure (PH2), arterial blood is loaded with H2 at the gas exchange surface. Alveolar PH2 is somewhat less than the ambient PH2, until the tissues have reached saturation. The PH2 of alveolar gas is determined by a balance of 2 processes: removal of H2 by capillary blood to the tissues and replenishment by alveolar ventilation. By the time blood leaves the lungs, alveolar PH2 is in complete equilibrium with the blood (Pblood). H2 is delivered to all regions of the body via the arterial blood supply. H2 in the blood diffuses into tissues until the blood and tissues are in equilibrium. Saturation is reached when Pblood, Ptis, and Pamb are all equal.   In the treated animal (A, bottom), the process is slightly different. Microbial metabolism of H2 in the intestines works as a sink, resulting in a mixed venous return that will have a lower PH2 compared with the untreated animal. As venous blood reaches the lung, alveolar PH2 will be somewhat lower than in the untreated animal, assuming that the ventilation is the same in the treated and untreated animal. The result is a slightly lower arterial PH2 compared with the untreated animal. Overall magnitude is dependent on the efficacy of the metabolism of H2. The result is a prolonged time to equilibrium, a somewhat lower PtisH2 at equilibrium, and a chronically subsaturated state.   During decompression (B, top), fluxes of inert gases are reversed. In the case of the untreated animal, H2 is transported from tissues to blood and via the blood to the gas exchange surface. This leads to a continuous decrease in the PtisH2. In a treated animal, there is a dual removal of H2: one from the intestines and one from the gas exchange surface. This enhances the reduction of the PtisH2 in the treated animal. The lower overall PtisH2 of the mixed venous blood in the treated animal (B, bottom) may result in a lower risk of bubble formation compared with the untreated animal, thereby reducing the probability of DCS.

The parameters were determined from the data by fitting the estimated P(DCS) to the actual outcome for each hyperbaric exposure, using Eq. 2 or 4. If the animal did not suffer from DCS by the end of the 1-h postdecompression observation period (Tend), the hyperbaric exposure was considered to be safe throughout the experiment.

Null Models

A special case of the risk function model is a more general and simplified model in which there are no explanatory variables. This model, referred to as Null, or Constant Hazard, considers the risk to be constant (c) at all times. For both the Immediate and Delayed models before decompression begins
r=0 (9A)
and after decompression begins
r=c (9B)
The scaling factor G is set to 1 and not considered a parameter.

Model Analysis

A Marquardt nonlinear parameter estimation routine using maximum likelihood was used to search for best fitting parameters (31). The likelihood ratio test was used to determine significance of parameters compared with the Null models (31) and between nested models (7). In this test, significance is defined by increases in the log-likelihood (LL) values of the models (i.e., significantly smaller negative LL values). A grid search over all floating parameters, involving over 500 unique starting parameter-value sets, was used to increase the likelihood of finding global rather than local LL maxima. Differences were considered significant at the P < 0.05 level.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

The cumulative distribution of the 53 cases of DCS vs. time is shown in Fig. 3. The data show that the occurrence of DCS was tightly clustered within the first half hour after reaching the observation pressure (11 atm, Fig. 3). No cases of DCS were observed during the transition to lower chamber pressure (Fig. 3).


View larger version (15K):
[in this window]
[in a new window]
 
Fig. 3.   Cumulative number of DCS cases with time after observation pressure (11 atm) is reached for the 53 observed cases of DCS out of 109 hyperbaric H2 exposures.

The results for the Immediate (Eqs. 4 and 5) and Delayed (Eqs. 4 and 6) models are summarized in Table 2. Each model was tested with varying numbers and combinations of the model parameters (Imm 1-4 for the Immediate model and Del 1-6 for the Delayed model). Both models showed improvement compared with the Null model (Table 2, Null vs. Imm1 and Null vs. Del1). Therefore, incorporating a specific description of the dive history significantly improved the fit to the data. The LL values cannot be used to compare the Immediate and Delayed models with each other because neither is a subset of the other.

                              
View this table:
[in this window]
[in a new window]
 
Table 2.   Parameter estimates and LL for Imm and Del models for the whole data set

Addition of the threshold parameter Thr to the risk function improved the fit only for the Delayed model (Table 2, Del2, and Del5). The standard errors of the Thr parameter estimates were nearly equal to or greater than the parameter estimates themselves (Table 2). Likelihood ratio analysis of this determined that the 95% confidence limit, as based on change in the LL by 1.95 units, was asymmetric about the Thr parameter (data not shown).

Inclusion of the parameter for H2 metabolic activity injected into the animals (AInj) improved the fit for both the Immediate and Delayed models (Table 2, Imm3, Del3, and Del5). The parameter for CH4 release rate (AVCH4 ) improved the fit only for the Delayed model (Table 2, Del4). Only the Delayed model could successfully be fitted using four parameters, including both Thr and AInj as best fit parameters (Table 2, Del5). There was a trend for an improvement in fit when both AVCH4 and Thr were included in the Delayed model (P < 0.10, Table 2, Del6). Activity injected and VCH4 were significantly correlated with each other (r = 0.60, P < 0.001; Fig. 4); thus no model was tested that included both AInj and AVCH4.


View larger version (17K):
[in this window]
[in a new window]
 
Fig. 4.   Rate of release of CH4 from animals vs. microbial activity injected into animals. Line represents least squares linear regression (y = 41.2 + 0.048x; r = 0.60, P < 0.0001). DCS outcome is indicated for each animal.

The observed vs. predicted DCS incidences for both models were tested by separating the data into groups divided by treatment or by varying compression and decompression sequences (Fig. 5). The Immediate (Imm3) and Delayed (Del4 and Del5) models satisfactorily predicted the number of DCS cases when the data set was divided into two groups (all controls vs. all treated, chi 2 test, P > 0.4-0.8; Fig. 5A), four groups (controls at 3 different decompression rates vs. all treated animals, P > 0.3-0.7; Fig. 5B), or 14 different groups (controls vs. treated using all dive sequences listed in Table 1, P > 0.2-0.7; Fig. 5C). There was no apparent advantage to any one of these models over the other two (Fig. 5).


View larger version (10K):
[in this window]
[in a new window]
 
Fig. 5.   Observed number of cases of DCS vs. number of cases predicted by the Immediate (Imm3) and Delayed (Del4 and Del5) models. A: data divided into control and treated animal groups only. B: data divided into all treated animals and controls at either 0.45, 0.90, or 1.8 atm/min decompression rate. C: data divided into all 14 of the control and treated animal groups at all compression and decompression sequences used. Dotted line in each panel has a slope of unity, representing perfect fit of predictions to observations. All 3 models could satisfactorily predict DCS incidence values that were not different from the observed values, regardless of the manner of assigning groups (P > 0.20).


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

DCS is notoriously difficult to diagnose because its manifestations are so varied and nonspecific (16). Choosing compression and decompression sequences with severe and rapid-onset outcomes as in this study helps reduce the subjectivity of the diagnosis. In this study, most DCS cases were unambiguous seizures or limb paralysis that occurred within the first 20 min after arrival at the lower chamber pressure (Fig. 3). This is similar to other studies on DCS in pigs (9), sheep (2), and rats (20). Previous experience with a pig model has indicated that virtually all manifestations of severe DCS will be evident within 1 h, with observation periods of 4-24 h adding more cases in less than 1% of tested animals (9). In a sheep model of DCS (2), a small number of cases continued to be diagnosed in the second and third hour of observation, but these may have been less severe symptoms. In DCS research in humans, symptoms may emerge as late as 20 h after return to 1 atm (2, 34). The prolonged latency observed in humans but not found in other animals may be an experimental artifact, because only obvious and early signs of DCS can be evaluated in other animals. Furthermore, other animals have only limited ability to report symptoms. Postdecompression, most pigs in this study manifested livid marks on the skin that are typical of skin DCS (6). We cataloged these marks for future reference but did not consider them sufficient to diagnose a case of DCS. This is in keeping with the practice for human divers of not offering recompression treatment for such mild symptoms alone.

Probabilistic models for DCS have been used successfully to model human and sheep data (2, 22, 26, 31, 33, 34). The fitted values were made from the observed outcome and the physical history of the hyperbaric experiment. Physiological interpretation of the parameters from these past studies must be made with care, because no actual physiological measurements were included. In contrast, in this study the proposed mechanism of biochemical decompression has been incorporated into the models. This mechanism is described in detail in Fig. 2. We were able to include in the models a measured physiological variable that is subject to manipulation, namely an estimate of the rate at which intestinal microbes were eliminating a portion of the body burden of H2. The inclusion of such a variable will permit an evaluation of the beneficial effects of H2 biochemical decompression in pigs from any compression and decompression sequence in H2 and will be one of the first demonstrations of a mathematical model using a measured physiological variable to predict DCS.

Various models were tested and not found to have better fitting ability or any additional insights compared with those reported here. Among these were models with more than one tissue, models with a separate tau  for H2 and He, models with a separate tau  during compression and decompression, models with a T1 set by the observers at 11 atm, and models incorporating body weight and chamber temperature.

The successful incorporation of a parameter for H2 metabolism in the models supports our hypothesis (19) that the H2-metabolizing microbes reduced P(DCS) by causing an overall reduction in PbloodH2 (Eq. 8). Two terms for H2 metabolism were tested (Table 2): one relating the removal of H2 in the pigs' tissues to the in vitro measurement of the total activity of H2-metabolizing microbes injected into the intestines of the animals (AInj) and the other relating the removal of H2 to the measured CH4 release rate (AVCH4 ) from the pigs. These two estimates of microbial H2 elimination were significantly correlated with each other (Fig. 4), with VCH4 ~10% of Inj (range 4-30%).

Inclusion of the parameter AInj, along with tau  and G, improved the fit to the data for both the Immediate and Delayed models, compared with models with only two parameters (tau  and G, Table 2, Imm1 and Del1). The parameter AVCH4, on the other hand, improved the fit only for the Delayed model. Using VCH4 had lower predictive power than using Inj in the Delayed model (Table 2, compare LL values for Del3 vs. Del4 and Del5 vs. Del6). This probably indicates a greater accuracy in measuring the metabolic activity of the microbes in vitro before their injection compared with measuring the release of CH4 into the chamber by the animals (18). In the latter case, multiple steps separate the evolution of a CH4 molecule inside the animal from sampling it by the gas chromatograph, with unknown kinetics and sampling errors at each of these steps. Accurately estimating the VCH4 during the hyperbaric experiment is complicated by the large volumes of chamber gas (chamber gas volume at 24 atm is equivalent to 135,000 liters at 1 atm) and the necessarily slow exhaust rate of the chamber (18). Furthermore, the metabolism of H2 within the intestinal ecosystem is more complex than the metabolic pathway shown in Eq. 1. Although methanogenesis is the predominant pathway for microbial H2 metabolism (21, 25), reduction of sulfate and nitrate and formation of acetate may also consume H2 (13, 21, 25). The potential that not all H2 metabolism was due to methanogenesis by M. smithii may mask part of the correlation between VCH4 and DCS outcome (Table 2; Fig. 4).

The total H2 tissue burden eliminated by the intestinal H2-metabolizing microbes is a crucial link between the PtisH2 and the DCS risk. The model (Del5) provides us with a way to assess the fraction of dissolved H2 (Ptis; that is, PtisH2 + PtisHe) removed by means of microbial metabolism. For this calculation, a treated animal compressed to 24 atm for 3 h and decompressed at 0.90 atm/min will be used (Table 1). By using the parameter for AInj, the difference in Ptis for an animal with or without intestinal injections of M. smithii can be estimated. Let us consider an animal with an Inj of 1.50 mmol CH4/min (activity injection in the upper range of those used in the study; Fig. 4). The estimated Ptis from this computation was compared with an animal with an Inj = 0 mmol CH4/min. On the basis of the model computations, it can be shown that the fraction of the total gas burden removed with additional injections of M. smithii is ~5%. In a previous study of H2 biochemical decompression in rats, there was a 50% reduction in DCS incidence associated with an elimination by the intestinal H2-metabolizing microbes of ~5% of the estimated total gas burden (19). The calculation in rats depended on numerous assumptions and simplifications regarding H2 solubility in tissues. The calculation performed above for pigs used fewer assumptions but arrived at a similar numerical result. Calculations based on H2 solubility in pigs have generated results that are on the order of 2-19% reduction of total gas burden for similar decreases in DCS incidence (18).

A graphical example of the possible reduction in the Ptis using Inj and Thr with and without a high dose of H2-metabolizing microbes is shown in Fig. 6. Because risk of DCS is the area under the curve between Ptis and Pamb, where Ptis > Pamb, the beneficial effects of biochemical decompression can be seen (Fig. 6). Even though the two curves for Ptis for control and treated animals are close, the difference in the area under the curve between Ptis and Pamb is appreciable for DCS risk [Fig. 6, compare P(DCS)C vs. P(DCS)T]. This illustrates how elimination of relatively small fractions of dissolved gas may have a surprisingly large impact on the DCS incidence.


View larger version (18K):
[in this window]
[in a new window]
 
Fig. 6.   Ambient pressure (Pamb) and tissue tension of H2 + He for a control animal (PtisC) and a treated animal (PtisT) injected with a total H2-metabolizing activity of 1.50 mmol CH4/min and the estimated probability of DCS for a control [P(DCS)C] and treated animal [P(DCS)T] using the 4-parameter Delayed model (Table 2, Del5). Hyperbaric exposure was to a total pressure of 24 atm for 3 h with a decompression rate of 0.9 atm/min.

The models described here can be used to predict the overall benefit from H2 biochemical decompression for any specified compression and decompression sequence. On the basis of the parameter estimates, it is possible to predict the change in the P(DCS) with varying VCH4 or Inj for any given hyperbaric sequence. For example, the most tested compression and decompression sequence in this study was to a total pressure of 24 atm for 3 h with a decompression rate of 0.9 atm/min (n = 20 control and 18 treated animals, Table 1). The observed DCS incidence was 80% (16 cases) and 39% (7 cases) for control and treated animals, respectively (Table 1). The predicted values for the Delayed (Del5) model with inclusion of Inj and Thr were 73% (14.6 cases) and 41% (7.4 cases) for control and treated animals, respectively (Fig. 5C). The reduction in DCS risk predicted by the Delayed model with inclusion of Inj is thus close to the observed value. The Imm3 and Del4 model predictions were only 1-2 DCS cases different from the Del5 model predictions (Fig. 5C). Naturally, predictions of P(DCS) made for dive exposures or treatment dosages beyond our data set must be made with caution and tested empirically before the full benefit of this model can be evaluated.

In conclusion, a hypothesis of the mechanism supporting biochemical decompression (Fig. 2) has been postulated and incorporated into a physiologically based probabilistic model (Eq. 8). The model can be used to predict the beneficial effects of biochemical decompression from any compression and decompression sequence with pigs in hyperbaric H2. The physiological interpretation of the model gives a foundation for further research including other physiological measurements related to gas transport kinetics and estimates of tissue inert gas content. The model supports our hypothesis that removal of H2 by microbial metabolism reduces the body burden of gas, thereby reducing the DCS risk by as much as 50%.


    ACKNOWLEDGEMENTS

We gratefully thank Diana Temple for editorial assistance and critical reading of this manuscript. We also thank A. Mulligan for graphical help with Fig. 2. We are grateful for the comments by the anonymous referees that we believe helped improve the readability of the paper. Also, we would like to thank C. McClure, H. Hung, and S. Maritz for sharing statistical expertise.


    FOOTNOTES

This work was funded by the Naval Medical Research and Development Command Work Unit no. 61153N MR04101.00D-1103. The opinions and assertions contained herein are the private ones of the authors and are not to be construed as official or reflecting the views of the Navy Department and the naval service at large.

Data were taken from animal experiments previously reviewed and approved by the Institutional Animal Care and Use Committee according to the principles set forth in the Guide for the Care and Use of Laboratory Animals, Institute of Laboratory Animal Resources, National Research Council, National Academy Press, 1996.

Address for reprint requests and other correspondence: S. R. Kayar, Environmental Physiology Dept., Naval Medical Research Center, 503 Robert Grant Ave., Silver Spring, MD 20910-7500.

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.

Received 30 October 2000; accepted in final form 23 August 2001.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

1.   Antonisen, NR, and Fleetham JM. Ventilation: total, alveolar, and dead space. In: Handbook of Physiology. The Respiratory System. Gas Exchange. Bethesda, MD: Am. Physiol. Soc, 1987, sect. 3, vol. IV, chapt. 7, p. 113-129.

2.   Ball, R, Lehner CE, and Parker EC. Predicting risk of decompression sickness in humans from outcomes in sheep. J Appl Physiol 86: 1920-1929, 1999[Abstract/Free Full Text].

3.   Broome, JR, and Dick EJ, Jr. Neurological decompression illness in swine. Aviat Space Environ Med 67: 207-213, 1996[Medline].

4.   Broome, JR, Dutka AJ, and McNamee GA. Exercise conditioning reduces the risk of neurologic decompression illness in swine. Undersea Hyperb Med 22: 73-85, 1995[Medline].

5.   Broome, JR, Pearson RR, and Dutka AJ. Failure to prevent decompression illness in rats by pretreatment with a soluble complement receptor. Undersea Hyperb Med 21: 287-295, 1994[Web of Science][Medline].

6.   Buttolph, TB, Dick EJ, Jr, Toner CB, Broome JR, Williams R, Kang YH, and Wilt NL. Cutaneous lesions in swine after decompression: histopathology and ultrastructure. Undersea Hyperb Med 25: 115-121, 1998[Medline].

7.   Collett, D. Modelling Survival Data in Medical Research. London, UK: Chapman & Hall, 1994.

8.   Cross, S, Jennings K, and Thomson L. Decompression sickness: role of patent foramen ovale is limited. BMJ 309: 743-744, 1994[Free Full Text].

9.   Dromsky, D, Toner CB, Survanshi S, Fahlman A, Parker E, and Weathersby P. The natural history of severe decompression sickness after rapid ascent from air saturation in a porcine model. J Appl Physiol 89: 791-798, 2000[Abstract/Free Full Text].

10.   Fahlman, A, Kayar SR, Becker WJ, Lin WC, and Whitman WB. Decompression sickness risk correlated with activity of H2-metabolizing microbes injected in pigs prior to dives in H2. Undersea Hyperb Med 26: 20, 1999.

11.   Fahlman, A. On the Physiology of Hydrogen Diving and Its Implication for Hydrogen Biochemical Decompression (PhD thesis). Ottawa, Ontario, Canada: Carleton University, 2000.

12.   Germonpré, P, Dendale P, Unger P, and Balestra C. Patent foramen ovale and decompression sickness in sports divers. J Appl Physiol 84: 1622-1626, 1998[Abstract/Free Full Text].

13.   Gibson, GR, Cummings JH, and Macfarlane GT. Competition for hydrogen between sulphate-reducing bacteria and methanogenic bacteria from the human large intestine. J Appl Bacteriol 65: 241-247, 1988[Medline].

14.   Hills, BA. Decompression Sickness. New York: Wiley, 1977.

15.   Homer, LD, Weathersby PK, and Survanshi S. How countercurrent blood flow and uneven perfusion affect the motion of inert gas. J Appl Physiol 69: 162-170, 1990[Abstract/Free Full Text].

16.   Jain, KK. Textbook of Hyperbaric Medicine. Toronto, Ontario, Canada: Hogrefe & Huber, 1990.

17.   Kayar, SR, Aukhert EO, Axley MJ, Homer LD, and Harabin AL. Lower decompression sickness risk in rats by intravenous injection of foreign protein. Undersea Hyperb Med 24: 329-335, 1997[Web of Science][Medline].

18.   Kayar, SR, Fahlman A, Lin WC, and Whitman WB. Increasing activity of H2-metabolizing microbes lowers decompression sickness risk in pigs during H2 dives. J Appl Physiol 91: 2713-2719, 2001[Abstract/Free Full Text].

19.   Kayar, SR, Miller TL, Wolin MJ, Aukhert EO, Axley MJ, and Kiesow LA. Decompression sickness risk in rats by microbial removal of dissolved gas. Am J Physiol Regulatory Integrative Comp Physiol 275: R677-R682, 1998[Abstract/Free Full Text].

20.   Lillo, RS, Parker EC, and Porter WR. Decompression comparison of helium and hydrogen in rats. J Appl Physiol 82: 892-901, 1997[Abstract/Free Full Text].

21.   Miller, TL. Biogenic sources of methane. In: Microbial Production and Consumption of Greenhouse Gases: Methane, Nitrogen Oxides, and Halomethanes, edited by Rogers JE, and Whitman WB.. Washington, DC: Am. Soc. Microbiol., 1991, p. 175-187.

22.   Parker, EC, Survanshi SS, Massell PB, and Weathersby PK. Probabilistic models of the role of oxygen in human decompression sickness. J Appl Physiol 84: 1096-1102, 1998[Abstract/Free Full Text].

23.   Robertson, AG. Decompression sickness risk in women (letter). Undersea Biomed Res 19: 216-217, 1992[Medline].

24.   Stevens, DM, Gartner SL, Pearson RR, Flynn ET, Mink RB, Robinson DH, and Dutka AJ. Complement activation during saturation diving. Undersea Hyperb Med 20: 279-288, 1993[Web of Science][Medline].

25.   Strocchi, AJ, Furne K, Ellis CJ, and Levitt MD. Competition for hydrogen by human faecal bacteria: evidence for the predominance of methane producing bacteria. Gut 32: 1498-1501, 1991[Abstract/Free Full Text].

26.   Thalmann, ED, Parker EC, Survanshi SS, and Weathersby PK. Improved probabilistic decompression model risk predictions using linear-exponential kinetics. Undersea Hyperb Med 24: 255-274, 1997[Web of Science][Medline].

27.   Van Liew, HD, and Hlastala MP. Influence of bubble size and blood perfusion on absorption of gas bubbles in tissues. Respir Physiol 7: 111-121, 1969[Web of Science][Medline].

28.   Vann, RD. Mechanisms and risks of decompression. In: Diving Medicine, edited by Bove AA, and Davis JC.. Philadelphia, PA: Saunders, 1990, p. 29-49.

29.   Weathersby, PK, Hartn BL, Flynn ET, and Walker WF. Role of oxygen in the production of human decompression sickness. J Appl Physiol 63: 2380-2387, 1987[Abstract/Free Full Text].

30.   Weathersby, PK, and Homer LD. Solubility of inert gases in biological fluids and tissues: a review. Undersea Biomed Res 7: 277-296, 1980[Web of Science][Medline].

31.   Weathersby, PK, Homer LD, and Flynn ET. On the likelihood of decompression sickness. J Appl Physiol 57: 815-825, 1984[Abstract/Free Full Text].

32.   Weathersby, PK, Survanshi SS, Hays JR, and MacCallum ME. Statistically based decompression tables III: Comparative Risk Using U. S. Navy, British, and Canadian Standard Schedules. Bethesda, MD: Naval Medical Research Institute, 1986. (Technical Report NMRI 86-50).

33.   Weathersby, PK, Survanshi SS, Homer LD, Hart BL, Nishi RY, Flynn ET, and Bradley ME. Statistically Based Decompression Tables I. Analysis of standard air dives: 1950-1970. Bethesda, MD: Naval Medical Research Institute, 1985, NMRI 85-16.

34.   Weathersby, PK, Survanshi SS, Homer LD, Parker E, and Thalmann ED. Predicting the time of occurrence of decompression sickness. J Appl Physiol 72: 1541-1548, 1992[Abstract/Free Full Text].

35.   Zhang, J, Fife CE, Currie MS, Moon RE, Piantadosi CA, and Vann RD. Venous gas emboli and complement activation after deep repetitive air diving. Undersea Biomed Res 18: 293-302, 1991[Medline].


J APPL PHYSIOL 91(6):2720-2729



This article has been cited by other articles:


Home page
J. Appl. Physiol.Home page
R. T. Mahon, H. M. Dainer, M. G. Gibellato, and S. E. Soutiere
Short oxygen prebreathe periods reduce or prevent severe decompression sickness in a 70-kg swine saturation model
J Appl Physiol, April 1, 2009; 106(4): 1459 - 1463.
[Abstract] [Full Text] [PDF]


Home page
J. Exp. Biol.Home page
A. Fahlman, A. Schmidt, D. R. Jones, B. L. Bostrom, and Y. Handrich
To what extent might N2 limit dive performance in king penguins?
J. Exp. Biol., October 1, 2007; 210(19): 3344 - 3355.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
A. Fahlman, W. C. Lin, W. B. Whitman, and S. R. Kayar
Modulation of decompression sickness risk in pigs with caffeine during H2 biochemical decompression
J Appl Physiol, November 1, 2002; 93(5): 1583 - 1589.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
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 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 HighWire
Right arrow Citing Articles via Web of Science (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Fahlman, A.
Right arrow Articles by Kayar, S. R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Fahlman, A.
Right arrow Articles by Kayar, S. R.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online