- Research article
- Open Access
Effects of vagotomy on cardiovascular and heart rate variability alterations following chronic normobaric hypoxia in adult rabbits
Biological Researchvolume 51, Article number: 57 (2018)
chronic hypoxia increases basal ventilation and pulmonary vascular resistance, with variable changes in arterial blood pressure and heart rate, but it’s impact on heart rate variability and autonomic regulation have been less well examined. We studied changes in arterial blood pressure, heart rate and heart rate variability (HRV) in rabbits subjected to chronic normobaric hypoxia (CNH; PB ~ 719 mmHg; FIO2 ~ 9.2%) for 14 days and assess the effect of autonomic control by acute bilateral vagal denervation.
exposure to CNH stalled animal weight gain and increased the hematocrit, without affecting heart rate or arterial blood pressure. Nevertheless, Poincaré plots of the electrocardiographic R–R intervals showed a reduced distribution parallel to the line of identity, which interpreted as reduced long-term HRV. In the frequency domain, CNH reduced the very-low- (< 0.2 Hz) and high-frequency components (> 0.8 Hz) of the R–R spectrograms and produced a prominent component in the low-frequency component (0.2–0.5 Hz) of the power spectrum. In control and CNH exposed rabbits, bilateral vagotomy had no apparent effect on the short- and long-term HRV in the Poincaré plots. However, bilateral vagotomy differentially affected higher-frequency components (> 0.8 Hz); reducing it in control animals without modifying it in CNH-exposed rabbits.
These results suggest that CNH exposure shifts the autonomic balance of heart rate towards a sympathetic predominance without modifying resting heart rate or arterial blood pressure.
Hypoxic challenges induce physiological homeostatic responses that reduce the mismatch between oxygen supply and tissue demand. Peripheral and central chemoreceptors are responsible of generating systemic reflex ventilatory and cardiovascular adjustments to cope with a hypoxic challenge [1,2,3], although hypoxia is mainly sensed by the peripheral chemoreceptors (carotid and aortic bodies). In response to acute hypoxia, ventilation increases in frequency and/or volume, which in turn increased PaO2, reduced PaCO2 and increased pH, thus increasing oxygen gradient and tissue supply [3, 4]. Concomitantly, an increase in systemic arterial blood pressure and heart rate augments the oxygen delivery rate to the tissues, aiding to restore the normal oxygen supply. If the hypoxic challenge is maintained, as occurs in high altitude exposure, there is an increase in resting ventilation  and heart rate , and the development of pulmonary hypertension [7, 8]. However, systemic arterial blood pressure has been reported to remain unchanged in chronic hypobaric hypoxia [7, 9] or increased in chronic normobaric hypoxia [8, 10, 11]. Thus, cardiovascular response during long-term hypoxia appear to be modified by barometric pressure that may modify the physiologic responses to hypoxia [12,13,14]. It is generally accepted that the hypoxia-induced changes in cardiovascular parameters may occur in tandem with an increase in the sympathetic drive . Heart rate variability (HRV) indexes in the time and frequency domains have been widely used in humans and animals to estimate the autonomic balance to the heart rate without any invasive intervention . Spectral analysis of HRV has two major oscillatory components defined as: (1) a low-frequency (LF) band, related to the combined sympathetic and parasympathetic influences, and (2) a high-frequency (HF) band, related to vagal influences and respiratory sinus arrhythmia [16,17,18]. Spectral analysis has shown that chronic sustained hypobaric hypoxia modifies the LF and HF bands of HRV in humans  and experimental animals [20,21,22], even in the absence of systemic hypertension. Moreover, cats and rats exposed to chronic intermittent normobaric hypoxia, a condition that mimics the obstructive sleep apnea syndrome, showed early shifts in LF and HF components of the HRV that preceded the development of systemic hypertension [17, 18, 23]. Therefore, HRV analysis estimates autonomic changes following exposure to hypoxia. Despite the studies showing HRV alterations in response to hypobaric hypoxia less is known about the effect of chronic sustained normobaric hypoxia on HRV. Moreover, different physiological responses to the same low PO2 level [12, 14] or in response to exercise [24, 25] have been described in hypobaric versus normobaric hypoxic conditions, suggesting that hypoxia and hypobaria could evoke different physiological responses. Therefore, we studied the effects of chronic hypoxia (O2 inspiratory fraction (FIO2) ~ 9.2%) on autonomic control of heart rate in rabbits maintained in normobaric conditions (PB ~ 719 mmHg) for 2 weeks and assessed the balance in autonomic regulation of HRV with the elimination of vagal parasympathetic control after acute bilateral vagotomy.
Experiments were performed in male New Zealand White rabbits, 11 control and 8 exposed to chronic normobaric hypoxia (CNH) for 14 days. After recording the cardiovascular variables in baseline normoxic conditions both vagus were cut, in 6 control and 6 CHN rabbits to assess contribution of the parasympathetic tone on physiological parameters. Experiments were conducted in accordance to the guidelines of the National Fund for Scientific and Technological Research (FONDECYT, Chile) and the Guide for the Care and Use of Laboratory Animals (National Research Council of the National Academies, USA). Bio-Ethics Committees from the Facultad de Ciencias of the Universidad de Chile, and Facultad de Ciencias Biológicas of the Pontificia Universidad Católica de Chile approved the experimental protocols.
Chronic normobaric hypoxia exposure
In each experiment, two rabbits were simultaneously exposed to normobaric hypoxia as previously described  (Fig. 1a), with a total of eight rabbits exposed to CNH. Briefly, individually caged (D × H × W; 55 cm × 42 cm × 38 cm) rabbits were placed in a 300 L (60 cm × 50 cm × 100 cm) Lucite® chamber. The oxygen content (FIO2) was monitored continuously by an oxygen sensor (AX300, Teledyne Analytical Instruments, USA) whose output was connected to a programmable automatic controller (Zelio SR2 B121BD, Schneider Electric, France) that activated solenoid valves (2026BV172, Jefferson Solenoid Valves, USA), admitting N2 or compressed air to the chamber if FIO2 values were over ~ 9.8% or below ~ 8.7%, respectively; a relief output valve opened simultaneously with the input valves and remained opened for 40 s after the closure of the input valves. Additionally, mechanic relief valves opened whenever pressure inside the chamber exceeded 17 mmHg above the ambient pressure. Thus, the mean FIO2 was maintained near 9.2%. The atmosphere of the chamber was mixed by 4 fans. CO2 produced by the ventilation of the rabbits was trapped using CaCO3 (250 g) and the urinary ammonia with H3BO3 (60 g). The chamber was opened every other day for approximately 5 min to clean the cages, to replenish the food and water and to replace CaCO3 and H3BO3 when necessary. The pressure inside the chamber was measured continuously with a gauge transducer (Statham P20). Throughout the conditioning period (14 d), chamber pressure and FIO2 signals were recorded by an analog-to-digital acquisition system sampling at 1 Hz (DI-158U, DATAQ Instruments Inc., USA), and the chamber’s temperature was recorded at 5-min intervals by a data logger (EL-USB-2, Lascar Electronics Inc., USA). Reference barometric pressure values for Santiago (717.2 ± 0.3 mmHg; altitude: 567 m) and relative humidity (63.3 ± 2.5%) were obtained daily from the web database (http://188.8.131.52/climatologia/) of the Dirección Meteorológica de Chile of the Dirección General de Aeronáutica Civil.
The CNH chamber design and control system, FIO2, pressure and temperature dynamics over a 24 h period are display in Fig. 1. At the beginning of each hypoxic exposure, flushing the chamber for 4 min with N2 decreased the FIO2 from 21.11 ± 0.04% (8 experiments) to 9.87 ± 0.30%, with a time constant of 2.15 ± 0.05 min (n = 8), reaching a stable value of 9.14 ± 0.08% after one additional min. When the FIO2 level fell below ~ 8.7% (Fig. 1b), the system automatically injected compressed air in the chamber for 20 s, thus maintaining a mean FIO2 value of 9.18 ± 0.03% (14 periods of 24 h). The admission of air (or N2 when FIO2 increased over ~ 9.8%) induced a small rapid increase in the pressure inside the chamber (Fig. 1b) with a mean maximal value of 14.7 ± 0.1 mmHg (n = 980 admissions in 13 periods of 24 h), which returned to the baseline level within 10 min (Fig. 1b). The mean interval between consecutive air admissions to the chamber was 16.3 ± 1.4 min (10 periods of 48 h), with a maximal mean interval of 43.8 ± 2.2 min (10 periods of 48 h). Because of the admission of compressed gases into the chamber, the mean pressure increased in 1.3 ± 0.1 mmHg (13 periods of 24 h) within the chamber, resulting in a mean pressure of 718.5 ± 0.4 mm Hg (13 periods of 24 h) during the hypoxic challenge. On the other hand, the mean temperature within the chamber remained constant (23.0 ± 0.2 °C; 14 periods of 24 h; Fig. 1c), presenting periodic variations in a single experiment but with no constant pattern throughout the different experiments (Fig. 1c).
Rabbits were anesthetized initially with ketamine/xylazine (75/7.5 mg/kg, i.m.) and thereafter maintained with sodium pentobarbitone (6–9 mg/kg every 30 min, i.v.). The rabbits exposed to CNH (n = 8) were recorded within 5 h after their removal from the hypoxic chamber between 8:30 am and 01:00 pm, while control animals (n = 11) were recorded at 8:30 am. The rabbits were placed in supine position on a heated pad. We cannulated their trachea and connected the cannula to a heated tachograph which output was connected to a differential pressure transducer. The right saphenous vein and the left femoral artery were cannulated for drug application (i.e. sodium pentobarbitone) and arterial blood pressure measurement, respectively. The arterial cannula was connected to a low volume pressure transducer. The electrocardiogram (ECG) was recorded using the II lead. Three stainless steel 22G-needle electrodes were localized in the insertion of the right (G1) and left (GND) front legs, and in the left (G2) rear leg. The signals from the pressure transducer (pulse pressure) and the ECG were sampled at 2 kHz and recorded by our data acquisition system (PowerLab8, ADInstruments, Australia). Systolic, diastolic, pulse and mean blood pressure and heart rate were analyzed offline using the blood pressure module of data analysis software (LabChart Pro, ADInstruments, Australia). After recording the cardiovascular variables in baseline normoxic conditions, bilateral transection of both vagus nerves rostral to the nodose ganglia were performed in a subset of control (n = 6) and CNH (n = 6) rabbits chosen randomly to analyze the contribution of the parasympathetic tone on the measured parameters.
Heart rate variability
The ECG were recorded for 10 min in each animal and the R–R intervals series were exported from LabChart Pro as text file and processed using Kubios Software (v1.1 SP1) from the Biomedical Signal Analysis Group, University of Kuopio, Finland . Time domain analysis of the HRV was performed using Poincaré plots. Briefly, the standard deviation of the (R–R)n intervals respect to the next (R–R)n+1 interval was computed. Short-term variability was defined by the points above the contingency line (SD1) and the long-term variability correspond to the points along the contingency line (SD2). Power spectral analysis was used to calculate HRV in the frequency domain. Briefly, the R–R time series was detrended using a smoothness prior algorithm, with α = 800 and an interpolation rate of 8 Hz. The spectral bandwidth was 0–2.1 Hz, with a spectral resolution of 0.1 Hz, using a window width of 420 points with no overlap. The spectral analysis was performed over the whole spectrum without assuming specific ranges within the bandwidth. Spectral analysis of HRV has two major oscillatory components defined as a low-frequency (LF) band, thought to be related to the combine sympathetic and parasympathetic influences, and a high-frequency (HF) band related to vagal influences and the respiratory sinus arrhythmia [16,17,18]. Thus, the LF-to-HF ratio (LF/HF) of HRV is used as an indirect index of cardiac sympathovagal balance.
The frequency bands of R–R intervals spectra were defined according to the power spectra obtained from control rabbits: (i) very-low frequency (VLF), DC-0.2 Hz, (ii) low frequency (LF), 0.2–0.5 Hz and (iii) high frequency (HF): 0.5–2.1 Hz. Calculations considered the relative power of the LF band and the HF band normalized to the total spectral power (standardized units, S.U.); the LF/HF was calculated from these values.
Blood glucose levels
Venous blood was withdrawn using a Teflon catheter from the ear marginal vein, and blood glucose level was measured using commercial strips and a glucose meter (Accu-Chek® Active, Roche, Germany).
All data is presented as mean ± SEM. Statistical analysis of data was performed with GraphPad Prism (version 6.07 for Windows, GraphPad Software, San Diego California, USA). Non-normally distributed variables: body weight, hematocrit, SD1, SD2, LF, HF, and LF/HF ratio. Normally distributed variables: blood glucose, heart rate, cardiac cycle duration, cardiac contraction time, systolic pressure, diastolic pressure, HRV total power. Mean arterial pressure as well as HRV power bands were normally distributed before but not after vagotomy. Comparison between two groups was performed using parametric (Student´s t-tests) or nonparametric (Mann–Whitney or Wilcoxon test) tests, according to the data structure. Comparison between more than two groups was performed using nonparametric analysis of variance (Kruskal–Wallis) test. Statistical decision level was set at p < 0.05.
CNH effects on body weight, blood glucose and hematocrit
The mean body weight of control animals (1.87 ± 0.08 kg; n = 11) was like the initial mean weight of animals exposed to CNH (1.79 ± 0.10 kg; n = 8). After 14 days, the body weight of control rabbits increased significantly to 2.27 ± 0.10 kg (p < 0.05). In contrast, the average body weight of CNH exposed rabbits was not significantly modified after 2-wks of CNH exposure (1.79 ± 0.10 kg vs. 1.78 ± 0.11 kg). Moreover, at the end of the CNH exposure rabbits had lower body weights than predicted values based on unexposed rabbits (2.27 ± 0.10 kg, p < 0.01).
CNH affected hematocrit but did not alter blood glucose levels. The mean hematocrit increased from 40.3 ± 0.7% at the beginning of the CNH period to 55.8 ± 0.9% after 14 days of CNH (p < 0.001). The mean blood glucose levels, measured after anesthetic induction and prior to physiological recordings, were similar in control and CNH rabbits (7.99 ± 0.51 mmol/L, n = 5 and 8.55 ± 0.47 mmol/L, n = 5, respectively; p > 0.05).
Effects of CNH on resting cardiovascular variables in basal conditions
Table 1 contains the mean values for the cardiovascular variables for control and CNH rabbits measured during a 10-min recording of consecutive heart beats after the induction of anesthesia. Rabbits exposed to 14 days of CNH had no statistically significant differences (p > 0.05) in arterial pressure, heart rate, cardiac cycle duration and contraction time, compared to control rabbits. Specifically, heart rate and mean arterial as well as systolic and diastolic arterial pressures were slightly but not significantly increased in CNH rabbits compared to the values in control rabbits (Table 1). Similarly, arterial pulse pressure was slightly but not significantly reduced in CNH rabbits with respect to control animals (Table 1).
Cardiovascular parameters after bilateral vagotomy
Most cardiovascular variables were affected by bilateral vagotomy (Table 2) in control rabbits. After vagotomy, mean heart rate increased significantly from 252.10 ± 24.58 bpm to 281.00 ± 24.35 bpm (p < 0.05; Fig. 2a). Similarly, after vagotomy, mean systolic, diastolic, and pulse arterial pressure increased (p < 0.05), and mean arterial blood pressure also increased (p < 0.05). Conversely, bilateral vagotomy did not modify significantly these cardiovascular variables in the CNH rabbits (Table 2).
Changes in heart rate after bilateral vagotomy depended on the resting heart rate. In control rabbits, heart rate was almost unaffected at the extremes of the pre-vagotomy range (low, 194 or high, 364 bpm) but was increased maximally near midrange (Fig. 2b, c). In contrast, changes in heart rate in CNH rabbits were linearly related to the heart rate prior to vagotomy (Fig. 2b, c).
Heart rate variability (HRV): time-domain analysis
Analysis in intact conditions
Poincaré plots of the R–R interval series were constructed and the short-term (SD1; minor axis) and long-term (SD2; major axis) HRV variabilities were calculated. Chronic normobaric exposure did not modify the short-term, but reduced the long-term variability (Fig. 3a, b). Control rabbits had a SD1 mean value of 1.14 ± 0.18 ms and rabbits exposed to CNH had a mean SD1 value of 0.83 ± 0.13 ms (p > 0.05, Fig. 3c). On the other hand, CNH rabbits showed a 1.8-fold reduction in SD2 compared to control rabbits (3.54 ± 0.54 ms versus 6.42 ± 1.1 ms; p < 0.05, Fig. 3c).
Analysis after acute bilateral vagotomy
The effects of bilateral vagotomy on the time-domain analysis of HRV are illustrated in Fig. 4. Vagotomy in control (Fig. 4a) and CNH rabbit (Fig. 4b) had small effects on short-term (SD1), and long-term (SD2) HRV. In control rabbits, the mean value of SD1 was not modified by vagotomy, with an initial average of 1.32 ± 0.26 ms before and 0.99 ± 0.09 ms after denervation (p > 0.05, Fig. 4c). Control rabbits had mean SD2 values of 7.59 ± 1.80 ms before and 12.58 ± 3.12 ms after vagotomy (p > 0.05, Fig. 4c). Similarly, CNH rabbits had mean SD1 values of 0.72 ± 0.32 ms before and 0.98 ± 0.34 ms after bilateral vagotomy (p > 0.05, Fig. 4c). The mean values of SD2 in CNH rabbits were 3.57 ± 1.60 ms before and 9.75 ± 7.43 ms (p > 0.05, Fig. 4c). Thus, acute bilateral vagotomy did not modify HRV in the time-domain for both control and CNH-exposed rabbits.
Heart rate variability (HRV): frequency domain analysis
Analysis in intact conditions
Individual spectrograms of the variation in RR Intervals of each of the control and CNH rabbits are displayed in Fig. 5. We found that the spectra of control rabbits consistently had maximal power peaks in the very-low (< 0.2 Hz) and high- (0.8–1.9 Hz) frequency bands (Fig. 5a). Spectrograms from rabbits exposed to CNH showed a marked flattening of their power spectrum and a reduction of frequency components below 0.2 Hz and above 0.8 Hz (Fig. 5b), compared to spectra obtained from control rabbits (Fig. 5a). However, HRV total power was similar in control and CNH rabbits (220.8 ± 45.2 ms2 vs. 162.7 ± 35.0 ms2, control vs. CNH, respectively; p > 0.05). Total power was differentially disturbed in control and CNH-conditioned rabbits. Compared to control rabbits, the mean spectrogram of CNH rabbits had decreases (p < 0.01) in the power of the very low-frequency band (< 0.2 Hz, Fig. 6) and in the high frequency band (> 0.8 Hz, Fig. 6), and a concomitant increase (p < 0.01) in the power of the low frequency range (0.2–0.5 Hz, Fig. 6). Accordingly, the LF/HF ratio was greater in rabbits exposed to CNH compared to control rabbits (0.77 ± 0.26, 0.21 ± 0.07, respectively; p < 0.05). Indeed, rabbits exposed to CNH had a ~ 3.5-fold increase in the LF/HF ratio compared to control rabbits (Fig. 6).
Analysis after acute bilateral vagotomy
In control rabbits, the total power of the spectra remained comparable before and after vagotomy (319.2 ± 72.7 ms2 vs. 272.0 ± 96.8 ms2, control vs. CNH, respectively; p > 0.05). However, the total power was redistributed. Compared to their pre-vagotomy spectrograms (Fig. 7a), bilateral vagotomy: (i) decreased power in the high-frequency band, especially > 0.7 Hz; (ii) increased power slightly in the low frequency band, 0.2–0.5 Hz (Fig. 7b), and (iii) increased power in the very-low-frequency band, < 0.2 Hz (p < 0.01, Fig. 7c). Acute vagotomy did not affect the mean low-frequency peak (0.056 ± 0.002 Hz before and 0.055 ± 0.004 Hz after vagotomy; p > 0.05, Fig. 7c). The mean high-frequency peak decreased from 1.33 ± 0.17 Hz in the intact condition to 0.94 ± 0.22 Hz after denervation (p > 0.05). Thus, the LF band and the LF/HF ratio increased from 0.05 ± 0.01 and 0.11 ± 0.03 to before to 0.07 ± 0.02 and 0.25 ± 0.06 respectively, after bilateral vagotomy (p < 0.05, Fig. 7c).
In CNH rabbits, comparing their spectra pre and post vagotomy (Fig. 8), the total power of the spectra increased due to increases in power at very-low- and low- frequency bands (< 0.6 Hz). Average total power of the spectra increased after bilateral vagotomy from 134.7 ± 35.3 ms2 to 742.3 ± 210.9 ms2 (p < 0.05). The mean power spectrum increased in peak amplitude of the very-low-frequency band (0.07 ± 0.005 Hz vs. 0.8 ± 0.21 Hz, pre vs. post vagotomy, respectively). Two additional peaks were evident in the low frequency band at 0.3 and 0.4 Hz (Fig. 8b, c). Nevertheless, the power content of the high frequency band (0.6–2.0 Hz) remained small and unaffected by vagotomy (0.39 ± 0.16 ms2 vs. 0.28 ± 0.05 ms2, pre vs. post vagotomy, respectively; p > 0.05; Fig. 8c). Moreover, the LF band and the LF/HF ratio increased from 0.25 ± 0.08 and 0.87 ± 0.37 before vagotomy to 0.54 ± 0.12 and 7.18 ± 3 after vagotomy (p < 0.05; Fig. 8c).
The main findings of this study are that exposure of rabbits to CNH for 2 weeks: (i) prevented weight gain, (ii) decreased long-term variability of the HRV, (iii) shifted the spectral distribution of the HRV from the HF band towards LF band suggesting a reduction in parasympathetic modulation on HRV despite the absence of changes in arterial blood pressure and heart rate. Together, these results suggest that CNH induced autonomic deregulation characterized by sympathetic/parasympathetic imbalance.
We found that body weight did not increase significantly during the 14 days of CNH exposure. However, growth curves for control rabbits indicate that the animals should have increased their weight around 20% after 14 days [27, 28]. A slight decrease in weight has been observed in mice , rats  and piglets  under a similar CNH paradigm, with a fast recovery of weight and weight gain after 7 and 14 days, respectively, after returning to normoxia . Thus, the CNH exposure appears not to impair permanently the metabolism of adult rabbits as it does to adult rats exposed to postnatal hypobaric hypoxia . Stagnation of body weight and pulmonary hypertension has been associated with acute and chronic hypoxia [7, 33]. Pulmonary arterial pressure increases with hypoxia , but only hypobaric hypoxia increases lung fluid and protein flow  producing edema. Thus, animals exposed to CNH should not develop pulmonary edema that could confound the observed weight stagnation during the CNH exposure.
Mean hematocrit was similar to arterial hematocrit values reported for rabbits [34,35,36,37,38] and increased to approximately 56% after 14 days of CNH. Similar hematocrit increases have been reported in adult New Zealand White rabbits after 7–10 days of hypoxia [36, 38]. On the other hand, blood glucose levels were similar between control and CNH rabbits and like those previously reported for healthy rabbits behaving or under light sedation . Thus, rabbits exposed to CNH appear not to be subjected to a major stress in the hypoxic chamber. Moreover, because hyperglycemia reduces hypertension and modifies autonomic control in mice exposed to intermittent hypoxia , rabbits subjected to CNH have glycemic levels that should not interfere with responses to hypoxia.
Resting arterial blood pressure and heart rate were not significantly affected by the CNH exposure with respect to the control controls and were similar to those previously reported in rabbits [41,42,43,44,45,46]. Exposure of rats to 14 days of CNH does not affect the heart rate , and adult rats previously exposed to chronic hypobaric hypoxia (FIO2 = 0.12) for 3 weeks present the same heart rate as control animals in normoxia . Thus, heart rate appears not to be significantly modified by either normobaric or hypobaric hypoxic challenges. Systolic, pulse and mean arterial pressure increase in adult rats exposed to hypobaric hypoxia (FIO2 = 0.12) during 10 days after birth, without modification of diastolic pressure or heart rate , although more moderate hypoxic levels (FIO2 = 0.14) appear to have no effects on these cardiovascular parameters . Similarly, CNH increases both mean arterial blood and pulmonary pressure in rats . Acute hypoxia increases arterial blood pressure and reduces rabbit heart rate linearly with blood oxygen saturation , but these changes appear not to modify the cardiovascular function of rabbits exposed to CNH.
In control animals, short-term (SD1) and long-term (SD2) HRV values were lower to those previously reported for awaken New Zealand White rabbits  although our values are consistent with previous reports that long-term variability was also greater than short-term variability . This difference may arise from the cardiovascular and ventilatory effects of the anesthetic used in our experiments [48,49,50]. Rabbits exposed to CNH presented no significant changes in short-term variability but a significant reduction in long-term variability. A reduction in long-term variability without modification of short-term variability has been observed after heart isolation in rabbits . Thus, this reduction of long-term variability observed in CNH rabbits may indicate a reduction in cardiac autonomic control. Accordingly, we used HRV spectrogram analysis as a surrogate of autonomic control to the heart, as previously described [42, 44, 51,52,53,54]. In the present study, we found that HRV spectrograms of control rabbits showed the largest peaks in frequencies below 0.2 Hz and above 0.8 Hz, with mean maximal around 0.06 Hz and 0.97 Hz, respectively. Similar frequency band distribution has been reported in anesthetized rabbits . In our experiments, the changes in the higher frequency (> 0.8 Hz) components displayed by CNH rabbits appear not to be related to altered resting ventilation, since both control and CNH animals have similar ventilatory frequencies . Exposure to CNH results in flattened HRV spectra, with a reduction of peak amplitude below 0.2 Hz, an almost complete absence of peaks above 0.8 Hz and a concomitant increase of the 0.2–0.8 Hz band power. The mean HRV spectrum of CNH-rabbits shows a large and significant reduction in the power of the highest and lowest peak frequencies of the studied range, with a significant increase in the power in the 0.2–0.5 Hz range, with no changes in the total power content of the spectrum with respect to control animals. It has been shown that parasympathetic blockade reduces the power content of the higher frequency spectral components (up to 0.6 Hz), while sympathetic blockade reduces the power of the lower part of the spectrum (up to 0.44 Hz) in unrestrained awaken rabbits . Therefore, modification of the mean spectra observed in CNH rabbits suggests an important change in the balance of sympathetic-parasympathetic activity, giving rise to a new set-point in autonomic heart control.
Limitations of the study
All the physiological recordings performed in our study were made under anesthesia. It has been previously shown that mostly all anesthetic agents normally depress cardiorespiratory function. Then, it is plausible that our results may underestimate the cardiovascular consequences of CNH exposure. However, it is worth noting that for both control and CNH-exposed animals, special care was taken when using anesthesia to have both experimental groups under the same anesthesia depth (Stage III plane 2). Therefore, both groups were comparable and even in the presence of anesthesia we observed an effect of CNH exposure on cardiovascular regulation. Nevertheless, we believe that future studies should address the consequences of CNH exposure on freely moving awake rabbits. Another limitation is that we did not directly measure sympathetic or parasympathetic discharges to the heart. We used HRV analysis as an indirect method to assess cardiac autonomic function. However, we performed bilateral vagotomy to assess the contribution of vagal modulation on spectral HRV bands. New and future studies should completely unveil the effects of CNH on cardiac autonomic drive by using selective pharmacological blockage of cardiac sympathetic/parasympathetic activity.
Exposure to CNH produced no significant changes in cardiovascular variables but modified HRV. The HRV analysis showed that exposure to CNH is associated with increased long-term variability and a decrease in higher frequency components of the power spectra, suggesting a cardiac autonomic imbalance. Additionally, acute vagotomy further increase the HRV shift towards lower frequencies and larger power content strengthening the notion of autonomic imbalance.
chronic normobaric hypoxia
- FIO2 :
oxygen inspiratory fraction
high-frequency HRV band
heart rate variability
low-frequency HRV band
- PaO2 :
arterial oxygen partial pressure
- PaCO2 :
arterial carbon dioxide partial pressure
power spectral density
very-low frequency HRV band
Guyenet PG, Stornetta RL, Bayliss DA. Central respiratory chemoreception. J Comp Neurol. 2010;518:3883–906.
Iturriaga R, Alcayaga J. Neurotransmission in the carotid body: transmitters and modulators between glomus cells and petrosal ganglion nerve terminals. Brain Res Rev. 2004;47:46–53.
Teppema LJ, Dahan A. (The ventilatory response to hypoxia in mammals: mechanisms, measurement, and analysis. Physiol Rev. 2010;90:675–754.
Gonzalez C, Almaraz L, Obeso A, Rigual R. Carotid body chemoreceptors: from natural stimuli to sensory discharges. Physiol Rev. 1994;74:829–98.
Powell FL, Milsom WK, Mitchell GS. Time domains of the hypoxic ventilatory response. Respir Physiol. 1998;112:123–34.
Klausen K. Cardiac output in man in rest and work during and after acclimatization to 3,800 m. J Appl Physiol. 1966;21:609–16.
Rabinovitch M, Gamble W, Nadas AS, Miettinen OS, Reid L. Rat pulmonary circulation after chronic hypoxia: hemodynamic and structural features. Am J Physiol Heart Circ Physiol. 1979;236:H818–27.
Schwenke DO, Tokudome T, Shirai M, Hosoda H, Horio T, Kishimoto I, et al. Exogenous ghrelin attenuates the progression of chronic hypoxia-induced pulmonary hypertension in conscious rats. Endocrinology. 2008;149:237–44.
Hsieh YH, Jacono FJ, Siegel RE, Dick TE. Respiratory modulation of sympathetic activity is attenuated in adult rats conditioned with chronic hypobaric hypoxia. Respir Physiol Neurobiol. 2015;206:53–60.
Flor KC, Silva EF, Menezes MF, Pedrino GR, Colombari E, Zoccal DB. Short-term sustained hypoxia elevates basal and hypoxia-induced ventilation but not the carotid body chemoreceptor activity in rats. Front Physiol. 2018;9:134.
Moraes DJA, Bonagamba LGA, Costa KM, Costa-Silva JH, Zoccal DB, Machado BH. Short-term sustained hypoxia induces changes in the coupling of sympathetic and respiratory activities in rats. J Physiol. 2014;592:2013–33.
Richard NA, Koehle MS. Differences in cardio-ventilatory responses to hypobaric and normobaric hypoxia: a review. Aviat Space Environ Med. 2012;83:677–84.
Savourey G, Launay JC, Besnard Y, Guinet-Lebreton A, Alonso A, Sauvet F, et al. Normo or hypobaric hypoxia tests: propositions for the determination of the individual susceptibility to altitude illnesses. Eur J Appl Physiol. 2007;100:193–205.
Savourey G, Launay JC, Besnard Y, Guinet A, Travers S. Normo- and hypobaric hypoxia: are there any physiological differences? Eur J Appl Physiol. 2003;89:122–6.
Wolfel EE, Selland MA, Mazzeo RS, Reeves JT. Systemic hypertension at 4,300 m is related to sympathoadrenal activity. J Appl Physiol. 1994;76:1643–50.
Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. Circulation. 1996;93:1043–65.
Rey S, Del Rio R, Alcayaga J, Iturriaga R. Chronic intermittent hypoxia enhances cat chemosensory and ventilatory responses to hypoxia. J Physiol. 2004;560:577–86.
Rey S, Tarvainen MP, Karjalainen PA, Iturriaga R. Dynamic time-varying analysis of heart rate and blood pressure variability in cats exposed to short-term chronic intermittent hypoxia. Am J Physiol Regul Integr Comp Physiol. 2008;295:R28–37.
Cornolo J, Mollard P, Brugniaux JV, Robach P, Richalet JP. Autonomic control of the cardiovascular system during acclimatization to high altitude: effects of sildenafil. J Appl Physiol. 2004;97:935–40.
Fauchier L, Melin A, Eder V, Antier D, Bonnet P. Heart rate variability in rats with chronic hypoxic pulmonary hypertension. Ann Cardiol Angeiol. 2006;55:249–54.
Melin A, Fauchier L, Dubuis E, Obert P, Bonnet P. Heart rate variability in rats acclimatized to high altitude. High Alt Med Biol. 2003;4:375–87.
Meyer M, Rahmel A, Marconi C, Grassi B, Cerretelli P, Skinner JE. Stability of heartbeat interval distributions in chronic high altitude hypoxia. Integr Physiol Behav Sci. 1998;33:344–62.
Del Rio R, Moya EA, Iturriaga R. Carotid body and cardiorespiratory alterations in intermittent hypoxia: the oxidative link. Eur Respir J. 2010;36:136–42.
Boos CJ, O’Hara JP, Mellor A, Hodkinson PD, Tsakirides C, Reeve N, et al. A four-way comparison of cardiac function with normobaric normoxia, normobaric hypoxia, hypobaric hypoxia and genuine high altitude. PLoS ONE. 2016;11(4):e0152868.
Netzer NC, Rausch L, Eliasson AH, Gatterer H, Friess M, Burtscher M, et al. SpO2 and heart rate during a real hike at altitude are significantly different than at its simulation in normobaric hypoxia. Front Physiol. 2017;8:81.
Icekson G, Dominguez CV, Dedios VP, Arroyo J, Alcayaga J. Petrosal ganglion responses to acetylcholine and ATP are enhanced by chronic normobaric hypoxia in the rabbit. Respir Physiol Neurobiol. 2013;189:624–31.
Blasco A, Piles M, Varona L. A Bayesian analysis of the effect of selection for growth rate on growth curves in rabbits. Genet Sel Evol. 2003;35:21–41.
Larzul C, De Rochambeau H. Comparison of ten rabbit lines of terminal bucks for growth, feed efficiency and carcass traits. Anim Res. 2004;53:535–45.
Mosqueira M, Willmann G, Zeiger U, Khurana TS. Expression profiling reveals novel hypoxic biomarkers in peripheral blood of adult mice exposed to chronic hypoxia. PLoS ONE. 2012;7:e37497.
Neckář J, Szárszoi O, Herget J, Ošťádal B, Kolář F. Cardioprotective effect of chronic hypoxia is blunted by concomitant hypercapnia. Physiol Res. 2003;52:171–5.
Ananthakrishnan M, Barr FE, Summar ML, Smith HA, Kaplowitz M, Cunningham G, et al. l-Citrulline ameliorates chronic hypoxia-induced pulmonary hypertension in newborn piglets. Am J Physiol Lung Cell Mol Physiol. 2009;297:L506–10.
Ross B, McIntosh M, Rodaros D, Hébert TE, Rohlicek CV. Systemic arterial pressure at maturity in rats following chronic hypoxia in early life. Am J Hypertens. 2010;23:1228–33.
Levine BD, Kubo K, Kobayashi T, Fukushima M, Shibamoto T, Ueda G. Role of barometric pressure in pulmonary fluid balance and oxygen transport. J Appl Physiol. 1988;64:419–28.
Jones RT. Normal values for some biochemical constituents in rabbits. Lab Anim. 1975;9:143–7.
Harcourt-Brown FM, Baker SJ. Parathyroid hormone, haematological and biochemical parameters in relation to dental disease and husbandry in rabbits. J Small Anim Pract. 2001;42:130–6.
Jelkmann W, Beckman B, Fisher JW. Enhanced effects of hypoxia on erythropoiesis in rabbits following beta-2 adrenergic activation with albuterol. J Pharmacol Exp Ther. 1979;211:99–103.
Ouaknine-Orlando B, Samama CM, Riou B, Bonnin P, Guillosson JJ, Beaumont JL, et al. Role of the hematocrit in a rabbit model of arterial thrombosis and bleeding. Anesthesiology. 1999;90:1454–61.
Prestwich KN, Buss DD, Posner P. A new method for raising neonatal rabbits in a hypoxic environment. J Appl Physiol. 1984;57:1913–6.
Harcourt-Brown FM, Harcourt-Brown S. Clinical value of blood glucose measurement in pet rabbits. Vet Rec. 2012;170:674.
Chalacheva P, Thum J, Yokoe T, O’Donnell CP, Khoo MCK. Development of autonomic dysfunction with intermittent hypoxia in a lean murine model. Respir Physiol Neurobiol. 2013;188:143–51.
Aydin MD, Bayram E, Atalay C, Aydin N, Erdogan AR, Gundogdu C, et al. The role of neuron numbers of the petrosal ganglion in the determination of blood pressure: an experimental study. Minim Invasive Neurosurg. 2006;49:358–60.
Janoušek O, Ronzhina M, Scheer P, Nováková M, Provazník I, Kolářová J. HRV in isolated rabbit hearts and in vivo rabbit hearts. Comput Cardiol. 2010;37:923–6.
Korner PI, Edwards AWT. The immediate effects of acute hypoxia on the heart rate, arterial pressure, cardiac output and ventilation of the unanӕsthetized rabbit. Exp Physiol. 1960;45:113–22.
Manzo A, Ootaki Y, Ootaki C, Kamohara K, Fukamachi K. Comparative study of heart rate variability between healthy human subjects and healthy dogs, rabbits and calves. Lab Anim. 2009;43:41–5.
Mokrý J, Remeňová T, Javorka K. Changes in respiratory rate, blood pressure and heart rate variability in rabbits during orthostasis. Acta Vet Brno. 2006;75:3–12.
Sun S-Y, Wang W, Zucker IH, Schultz HD. Enhanced peripheral chemoreflex function in conscious rabbits with pacing-induced heart failure. J Appl Physiol. 1999;86:1264–72.
Rohlicek CV, Matsuoka T, Saiki C. Cardiovascular response to acute hypoxemia in adult rats hypoxemic neonatally. Cardiovasc Res. 2002;53:263–70.
Komatsu T, Singh PK, Kimura T, Nishiwaki K, Bando K, Shimada Y. Differential effects of ketamine and midazolam on heart rate variability. Can J Anaesth. 1995;42:1003–9.
Latson TW, McCarroll SM, Mirhej MA, Hyndman VA, Whitten CW, Lipton JM. Effects of three anesthetic induction techniques on heart rate variability. J Clin Anesth. 1992;4:265–76.
Zeller A, Arras M, Jurd R, Rudolph U. Identification of a molecular target mediating the general anesthetic actions of pentobarbital. Mol Pharmacol. 2007;71:852–9.
Moguilevski V, Oliver J, McGrath BP. Sympathetic regulation in rabbits with heart failure: experience using power spectral analysis of heart rate variability. Clin Exp Pharmacol Physiol. 1995;22:475–7.
Mokra D, Tonhajzerova I, Petraskova M, Calkovska A. Effects of dexamethasone on cardiovascular functions in acute phase in meconium-injured rabbits. Pediatr Int. 2009;51:132–7.
Nadareishvili KSh, Meskhishvili II, Kakhiani DD, Ormotsadze GL, Nazarishvili GT, Gvasalia MG, et al. Heart rate variability in Chinchilla rabbits. Bull Exp Biol Med. 2002;134:568–70.
Ronzhina M, Janousek O, Scheer P, Novakova M, Provaznik I, Kolarova J. Determination of the frequency bands for heart rate variability: studies on the intact and isolated rabbit hearts. Comput Cardiol. 2010;37:927–30.
Iwao T, Yonemochi H, Nakagawa M, Takahashi N, Saikawa T, Ito M. Effect of constant and intermittent vagal stimulation on the heart rate and heart rate variability in rabbits. Jpn J Physiol. 2000;50:33–9.
Alcayaga J, Del Rio R, Moya EA, Freire M, Iturriaga R. Rabbit ventilatory responses to peripheral chemoexcitators: effects of chronic hypoxia. Adv Exp Med Biol. 2012;758:307–13.
JA, RDR and RI conceived and designed the experiments. JA, RDR, EM and MF performed the experiments. JA and RDR analyzed the data. JA and RI contributed reagents/materials/analysis tools. JA wrote the draft paper, all authors contributed to the final version. All authors read and approved the final manuscript.
We thank Mr. Jorge Arroyo for the care and maintenance of the animals subjected to chronic normobaric hypoxia.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated and analyzed during this study are included in this published article.
Consent for publication
Ethics approval and consent to participate
All experiments were conducted in accordance to the guidelines of the National Fund for Scientific and Technological Research (FONDECYT, Chile) and the Guide for the Care and Use of Laboratory Animals (National Research Council of the National Academies, USA). Bio-Ethics Committees from the Facultad de Ciencias of the Universidad de Chile, and Facultad de Ciencias Biológicas of the Pontificia Universidad Católica de Chile approved the experimental protocols.
This work was supported by Grants 1090157 to J.A. and 1180172 to R.D.R. from the National Fund for Scientific and Technological Development (FONDECYT) of Chile.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
- Heart rate variability
- Chronic normobaric hypoxia
- Cardiovascular regulation
- Sympathetic control
- Vagal activity
- Poincaré plot
- Spectral components