Abstract
Objective
Acute hypobaric hypoxia is associated with autonomic changes that bring a global reduction of linear heart rate variability (HRV). Although changes in nonlinear HRV can be associated with physiologic stress and are relevant predictors of fatal arrhythmias in ischemic heart disease, to what extent these components vary in sudden hypobaric hypoxia is not known.
Methods
Twelve military pilots were supplemented with increasing concentrations of oxygen during decompression to 8230 m in a hypobaric chamber. Linear and nonlinear HRV was evaluated at 8230 m altitude before, during, and after oxygen flow deprivation. Linear HRV was assessed through traditional time-domain and frequency-domain analysis. Nonlinear HRV was quantified through the short-term fractal correlation exponent alpha (αs) and the Sample Entropy index (SampEn).
Results
Hypoxia was related to a decrease in linear HRV indexes at all frequency levels. A non-significant decrease in αs (basal, 1.39 ± 0.07; hypoxia, 1.11 ± 0.13; recovery, 1.41 ± 0.05; P = .054) and a significant increase in SampEn (basal, 1.07 ± 0.11; hypoxia, 1.45 ± 0.12; recovery, 1.43 ± 0.09; P = .018) were detected.
Conclusions
The observed pattern of diminished linear HRV and increased nonlinear HRV is similar to that seen in subjects undergoing heavy exercise or in patients with ischemic heart disease at high risk for ventricular fibrillation.
Introduction
Autonomic responses are crucial for adaptation to hypoxia. Acute hypoxia triggers several autonomic mechanisms, mainly in the respiratory system (eg, the hyperventilation started by the peripheral chemoreceptor reflex) and in the cardiovascular system (eg, sympathetic activation with heart rate and cardiac output increase to improve pulmonary and tissue perfusion). 1 Previous studies have evaluated heart autonomic response using heart rate variability (HRV) analysis in subjects exposed to acute hypobaric hypoxia in real 2 as well as simulated settings,3,4 usually introducing the changes in inspired O2 concentration gradually and at altitudes below 6400 m (21 000 ft).
Heart rate variability analysis is based on the fact that the autonomic control of the sinus rhythm adds variability to the heart rate signal at different frequency levels. Among these, in short-term studies (less than 10 minutes), the high-frequency (HF) component (0.15–0.4 Hz) is related to respiratory sinus arrhythmia mediated by parasympathetic activity, and the low-frequency (LF) component (0.04–0.15 Hz) is related to the baroreflex control including sympathetic and parasympathetic influences. In this context, the very-low-frequency (VLF) component (<0.04 Hz) is a yet undefined measure that may reflect changes in the trend of data. 5
In addition to the cyclic fluctuations previously described, there are complex fluctuations of the sinus rhythm that can be analyzed with nonlinear dynamics methods. It is believed that nonlinear fluctuations are determined by interactions of electrophysiologic, hemodynamic, and humoral variables, as well as by autonomic and central nervous regulation.5,6 Changes of nonlinear HRV are associated with several cardiovascular stressors. The decreased nonlinear HRV has been related to physical stress, 7 aging, 8 and a variety of pathologic conditions such as stable coronary artery disease. 9 On the other hand, increased nonlinear HRV was also associated with physiologic situations such as high exercise intensity levels 10 and complicated ischemic heart disease where it has been shown to be predictive for ventricular fibrillation and sudden cardiac death. 11
Most of the studies that assessed the relation between HRV and acute hypobaric hypoxia reported a reduced linear HRV pattern consistent with an increased sympathetic tone and decreased parasympathetic tone. 2 –4 Nonlinear indicators of HRV have been sparsely studied in these situations. The results indicate that nonlinear heart rate fluctuations also decrease in acute hypobaric hypoxia.12,13
Altitude-induced changes in cardiac rhythm could explain the significant rate of sudden cardiac death at high altitudes. 14 It has been reported that 30% of all deaths during mountain sports at altitude are attributable to sudden cardiac death, more likely in those who have had prior myocardial infarction, coronary artery disease, or coronary risk factors. 15 In addition, it has been seen that the occurrence of ventricular extrasystoles is proportional to the altitude during acute hypoxia exposure in a healthy elderly man. 16
However, contrary to what would be expected, the HRV pattern of diminished linear and nonlinear HRV observed when subjects are exposed to altitudes lower than 6400 m (21 000 ft) does not resemble that of patients with increased risk of death due to cardiac arrhythmias, who display reduced linear HRV but increased nonlinear HRV. 11 It is possible that this discrepancy is related to an insufficient hypobaric-hypoxic stress; therefore, we hypothesized that exposure to an altitude higher than 6400 m (21 000 ft) may result in an alteration of the HRV pattern more similar to that seen in patients with high risk of arrhythmia-related death.
In the current study, we examine the nonlinear HRV pattern associated with severe acute hypoxia aiming to obtain information that could help to understand the autonomic nervous system mechanisms underlying situations related to increased risk of cardiac arrhythmias. In this regard, the use of training protocols followed by aircrews to learn recognition of their symptoms of hypoxia at simulated altitudes of 8230 m (27 000 ft) in hypobaric chambers 17 allows the possibility to compare the autonomic responses with those already studied at lower altitudes.12,13
Materials and methods
Ethical Approval
The study conformed to the standards set by the Declaration of Helsinki, and it was approved in advance by the institutional committee for ethical review of the Argentine National Aerospace Medicine Institute, Buenos Aires, Argentina. All participants provided voluntary written informed consent before taking part. The data were collected noninvasively during routine training that the subjects had to undertake.
Subjects
We studied 17 healthy male military pilots who underwent a sudden hypoxic exposure as part of their standard training procedure. All subjects had been previously subjected to a full medical and psychological examination.
Hypoxic Protocol
The study was performed in a hyperbaric chamber (The Biggs Boiler Works Co., Akron, OH) adapted for use as a hypobaric chamber (Pittsburgh-Des Moines Steel Co., Des Moines, IA) at installations of the Argentine National Aerospace Medicine Institute. Participating crews entered into the chamber in groups of six together with a physician. Pressure was lowered to simulate an ascent to 8230 m (27 000 ft) at a rate of 640 m/min (2100 ft/min). During the ascent subjects received, through a pressure type mask, a mixture of oxygen and ambient air provided on demand by an oxygen regulator (USAF Type CRU-68/A; Aro Corporation, Bryan, OH).
The ratio of oxygen to air was automatically adjusted to supply increasing oxygen as altitude increased. At 8230 m (27 000 ft), all participants were receiving 100% oxygen. Then, oxygen flow was abruptly interrupted for one participant at a time. To detect psychomotor impairment, each participant had to write on a paper sheet numbers in a descending order starting from 999. As soon as the first symptom of hypoxia appeared, the oxygen flow was restarted, and the pilot resumed the writing. After the subject was fully recovered and reached basal oxygen saturation levels, the next pilot was evaluated. Oxygen saturation was measured continuously during the test by using a Healthdyne 930 pulse oximeter (Healthdyne Technologies Inc., Marietta, GA). After the test was run for all participants, the normal atmospheric pressure was gradually restored.
Heart rate variability was analyzed before (basal), during (hypoxia), and after (recovery) oxygen flow interruption. The beats corresponding with the “basal” stage were taken from the period immediately prior to the interruption of oxygen. The transition beats corresponding with the time between the interruption of the oxygen flow and the stabilization of the heart rate at an increased level (as determined by visual inspection) were not considered for the analysis. The beats corresponding with the “hypoxia” stage were taken from the elapsed time after the ending of the transition period and the restoration of oxygen flow. The transition beats corresponding with the elapsed time between the restoration of the oxygen flow and the stabilization of heart rate at basal levels (as determined by visual inspection) were not considered. Finally, the beats corresponding with the “recovery” stage were taken from the period immediately subsequent to this last transition period (Figure 1).

RR interval recording of one subject along the three epochs considered for analysis, with the transition periods between them.
Within each subject, the number of beats selected for the “basal” and “recovery” epoch was the same as the number of beats selected for the “hypoxia” epoch. We selected the same number of beats for all stages because in short time series of less than 200 beats, the error associated with nonlinear HRV analysis may vary as a function of data length.6,18 As the transition periods between stages were not evaluated, the records were free of pronounced upward or downward trends. 5
Signal Recording and HRV Analysis
Signal recording
Electrocardiogram signal was recorded using a digital Holter device (Holter HCAA 348; Holtech, Servicios Computados S.A., Buenos Aires, Argentina) and stored in a solid-state memory. Ventricular depolarization (R waves) were detected through the device software. The time elapsed between R waves (RR intervals) was then computed. We visually identified and manually tagged premature and lost beats in the original file of RR intervals. These abnormal beats were replaced by RR intervals resulting from linear interpolation. 5 Only those segments with >85% qualified beats were included in the analyses.19,20
Time-domain HRV analysis
Quantitative time series analysis was performed on heart rate by evaluating measures of variation over time. Among these, RRm (mean duration of RR intervals in milliseconds) quantifies the mean heart rate, SDNN (standard deviation of RR intervals in milliseconds) represents a coarse quantification of overall variability, and RMSSD (square root of the mean squared differences of successive normal RR) measures short-term heart rate variations.5,6
Spectral HRV analysis
These measurements provide an evaluation of the power (amplitude) of the contributing frequencies underlying HRV. 6 For the analysis of cyclic HRV components, the wavelet transform was chosen rather than the traditional fast Fourier transform (FFT) because it allowed a better fit of a signal to nonlinear fluctuations. The power of the signal along time was determined at previously defined frequency bands 21 : total area (TA; 0–1.2 Hz), very low frequency (VLF; 0–0.0375 Hz), low frequency (LF; 0.0375–0.15 Hz), and high frequency (HF; 0.15–1.2 Hz). Reported values are expressed as the natural logarithm of TA, HF, LF, and VLF; normalized for the total power; and as the ratio between LF and HF. Details of this procedure can be found elsewhere.21,22
Nonlinear HRV analysis
The short-term fractal correlation exponent αs 8 and the Sample Entropy index (SampEn) 23 were used as nonlinear HRV indexes. The scaling exponent αs, based on the “detrended fluctuation analysis,” quantifies the short-term (<11 beats) fractal correlation properties of the interbeat time data. Values of αs close to 0.5 are associated with uncorrelated RR intervals, whereas values close to 1.5 are associated with strong correlation between RR intervals. Values near 1 are characteristic of fractal-like processes, associated with the dynamic behavior of time series generated by complex systems, such as the autonomic regulation of the sinus rhythm in healthy subjects. 8 SampEn measures the degree of irregularity of the RR interval time series. Regular sequences will result in lower SampEn values, whereas random behavior is associated with larger SampEn values. These methods have been previously described. 6 ,8,23
Statistical Analysis
Values are expressed as mean ± standard error. Normality of distributions was verified by the Kolmogorov-Smirnov test.
Demographic characteristics, time spent in hypoxia, Sa
The HRV indexes obtained during baseline, hypoxia, and recovery periods were compared by means of a repeated measures analysis of variance (ANOVA) followed by a post hoc Bonferroni test to establish intergroup differences.
Additional associations of the time spent in hypoxia with the time spent on 100% O2 before the test and with HRV data were analyzed through a Pearson correlation test.
Results
The final sample included 12 subjects, because five individuals had hypoxia periods with less than 100 RR intervals. The demographic characteristics of included and excluded subjects are shown in Table 1.
Characteristics of subjects
BMI, body mass index.
Values are presented as mean ± standard error. Intergroup differences were assessed by means of Mann-Whitney U test.
During the hypoxia period, the included group showed a mean of 138 ± 13 RR intervals (minimum 100, maximum 253), and the excluded group showed a mean of 70 ± 3 RR intervals (minimum 63, maximum 80). The percentage of premature or missed beats was similar in both groups (less than 2% of the data).
The time that the pilots spent with 100% O2 before the test had no correlation with the time spent in hypoxia. Both included and excluded pilots reported similar hypoxia-related symptoms, including transient visual derangements, flushing in the face or limbs, and psychomotor lack of coordination. Table 2 depicts time, Sa
Stage characteristics in included and excluded subjects
I, initial Sa
Values are presented as mean ± standard error.
Differences between included and excluded subjects were assessed by means of a Mann-Whitney U test. Significant differences were found for the time covered by each stage ( a P = .004; b P = .003; c P = .008) and for the mean heart rate in the hypoxia and recovery stages ( d P = .027; e P = .004).
Heart rate variability analysis is shown in Table 3. Linear analysis revealed a marked increase of heart rate related to hypoxia associated with a global decrease of HRV at all frequency ranges. Once oxygen was restored, mean heart rate and global HRV tended to return to basal levels. This trend was not observed in the power spectral indexes of HRV, where the hypoxia period was different from the basal stage but similar to the recovery stage. No significant differences were observed in normalized units of HRV. Nonlinear analysis revealed a marginally significant effect of hypoxia in decreasing αs and a significant effect of hypoxia in increasing SampEn. Once oxygen was restored, SampEn value remained significantly different than the basal stage. Heart rate variability measurements at basal, hypoxia, and recovery stages had no correlation with the time spent in hypoxia.
Heart rate variability analysis
RRm, mean of RR interval duration; SDNN, standard deviation of RR intervals; RMSSD, square root of the mean squared differences of successive normal RR; TA, total area power; VLF, very-low-frequency power; LF, low-frequency power; HF, high-frequency power; αs, short-term fractal correlation exponent; SampEn, sample entropy.
Values are expressed as mean ± standard error.
Different from each other;
different from basal. ANOVA repeated measures followed by Bonferroni test.
Figure 1 shows a typical HRV recording along the three epochs considered for analysis (basal, hypoxia, and recovery) with the transition periods among them. Note the increase of heart rate (lower RR interval duration) and the decrease in the amplitude of heart rate oscillations in the hypoxia period (RRm = 594 milliseconds, SDNN = 34 milliseconds) in comparison with basal (RRm = 734 milliseconds, SDNN = 107 milliseconds) and recovery (RRm = 893 milliseconds, SDNN = 55 milliseconds) periods.
To account for differences in the mean and amplitude of heart rate fluctuations, Figure 2 shows the recordings corresponding with each epoch with RR intervals rescaled as standard deviation units centered at a mean of zero. When compared with basal period (SampEn = 0.89), a higher irregularity of the oscillations during hypoxia period (SampEn = 1.61) and afterwards (SampEn = 1.58) could be noted.

RR interval recording corresponding with basal, hypoxia, and recovery periods of the same subject of Figure 1. To account for differences in the mean and amplitude of heart rate fluctuations, RR intervals were rescaled as standard deviation units centered around a mean of zero.
Discussion
The main finding of the current study is that acute severe hypoxia is associated with greater nonlinear HRV, as evidenced by the observed higher values of SampEn, a measurement of the irregularity of the RR interval time series. This change was accompanied by a decrease of linear HRV of similar proportions in all frequency levels, as previously described by others. 2 –4,24
The marked fall in oxygen saturation and the appearance of subjective symptoms indicate that the pilots were exposed to considerable hypoxic stress. The current study differs from others in the sudden desaturation and the severity of the hypobaric hypoxia condition.12,13 In addition, comparison of hypoxia effects in each of the subjects and at the same simulated altitude reinforces the evidence that the observed changes were caused by hypoxemia.
Changes observed in the linear components of HRV could be related to the increase of heart rate due to the hypoxia effect, as the amplitude of the oscillations of different frequencies as well as the length of the RR interval diminished. 25 However, changes observed in nonlinear components of HRV may not be directly related to the changes in heart rate, especially when heart rate is high. 25 In fact, physical exercise has been associated with either decreased 7 or increased nonlinear HRV. 10 It was shown that during incremental exercise test, αs increased from rest up to an intensity level of approximately 40% of maximum oxygen consumption, declining thereafter linearly toward values that reflect no correlation between RR intervals.10,26 This pattern can be explained by the varying intensity of sympathetic input to the sinus node during different intensity levels of exercise. 26 Reduced values of αs are usually associated with increased values of SampEn, 25 as seen in the current study. The extreme activation of the sympathetic nervous system underlying severe hypoxia and vigorous exercise may lead to an increase in propensity to ventricular fibrillation. 14 ,27,28 Thus, the common HRV pattern of severe hypoxia described herein and the one of vigorous exercise described by others10,26 might reflect this increased risk.
An alternative explanation for the observed changes of nonlinear HRV dynamics after hypoxia is that they are related to intrinsic cardiac effects resulting from compensatory physiologic mechanisms induced by hypoxia. It is known that pathologic disruption of specific components of the intrinsic cardiac plexus can compromise the autonomic control of the heart, thus causing arrhythmogenesis. 29 It has been proposed that this system would act as a “low-pass filter” to smooth the effect of spurious extracardiac stimuli and to maintain a balance between the electrical and mechanical functions of the heart by means of local feedback circuits. 30 A functional alteration of this “filter” in normal subjects under hypoxia could produce the observed pattern of greater nonlinear irregularity that remained elevated even after oxygen was restored. In any case, an alteration of the complexity of dynamics of heart rate fluctuations as described in the current work can render the heart less adaptable for the changing environmental requirements or may start transient electric alterations favoring the appearance of arrhythmias.29,30 In this respect, it is important to note that our results are similar to the nonlinear HRV pattern described in ischemic heart disease patients with high risk for fatal arrhythmias. 11
Some limitations of the current study must be taken into account. The hypoxia period is a nonspontaneous and nonresting state. A maximal sympathetic activation and vagal withdrawal allowed the cardiovascular system to achieve a maximal constant cardiac output during a defined period of time, after which all variables rapidly decline toward prehypoxic levels. Thus, even if the hypoxia state is unstable, we can assume that the data were stationary in the temporal window considered. 31 In this regard, transition periods between states were not considered for analysis, and data were detrended to avoid possible trends. A further limitation is related to the application of nonlinear HRV indexes to short sets of HRV data. However, it was reported that the application of nonlinear indexes such as SampEn or approximate entropy (ApEn) to 100 data points series yields statistically reliable and reproducible results.6,23 Regarding the use of the short-term scaling exponent αs, it was found that, from a collective point of view, there is statistical agreement in real data between the means of the exponents derived from long and short segments involving the use of 100 data points. 18
In sum, this study shows that severe acute hypoxia is associated with a pattern of diminished linear HRV and increased nonlinear HRV, different from what has been seen in previous studies at lower altitudes. While this pattern is similar to that observed during vigorous exercise or in patients with ischemic heart disease at risk for fatal arrhythmias, future studies are needed to assess whether it is related to an increased risk of sudden death in healthy people.
Footnotes
Acknowledgments
We wish to thank the Uruguayan Air Force for its participation in the study.
Funding
This work was supported by a grant from the Agencia Nacional de Promoción Científica y Technólogia, Argentina (PICT 2007-01045). D.E.V., L.C.N.S., and D.P.C. are Research Career Awardees from the Argentine Research Council (CONICET).
