Introduction

Although the impairment of neural regulation in diabetes has been well documented and defined during the past four decades, the precise conditions in which such alterations occur and develop are still unclear. In particular, the extent and relative proportion to which the abnormalities observed are the result of morphological neuropathy or are functional in nature remain unclear.

Defining the threshold between functional and organic alteration could be of great practical importance. While many clinical trials failed to demonstrate improvements in diabetic autonomic neuropathy based on various rationales that assumed neural damage [1, 2], interventions based on the idea of functional abnormalities (e.g. physical training) proved to be effective [3, 4]. On a shorter timescale, it was shown that even a simple procedure such as slow breathing, an intervention known to reduce sympathetic drive and increase baroreflex sensitivity (BRS) [5, 6], increased BRS in diabetic patients with type 1 diabetes and subnormal baseline values [7].

Several possibilities may explain why cardiovascular modulation can be functionally impaired in patients with type 1 diabetes. These involve insulin treatment, hyperglycaemia, abnormal glycaemic control, previous hypoglycaemia and s-cortisol [812]. Additionally, hypoxia has recently attracted the attention of diabetologists as it seems to be implicated in the origin of important complications, such as diabetic nephropathy and diabetic retinopathy [1315]. During acute exposure to high altitude, hypoxia induces direct vasodilation and reflex sympathetic activation, chemoreflex activation and a reduction in BRS [16].

In this study, we tested the hypothesis that BRS, an established prognostic index in cardiovascular diseases [1719], is restrained by hypoxia in patients with type 1 diabetes. We reasoned that if this were the case, short-term inhalation of oxygen should restore BRS, and diabetic participants would be more responsive to the same intervention than healthy control participants. Conversely, in the presence of a neural lesion, and also if hypoxia were not responsible for the reduced BRS, diabetic individuals would not be responsive to hypoxia or else would respond to an extent similar to but not greater than the healthy participants. Several studies have shown that oxygen administration reduces heart rate and increases indices of heart rate variability related to parasympathetic activity [2023]. Finally, oxygen administration increases ventilation by a dose-dependent stimulation of the parasympathetic lung afferents [24, 25]. Two studies reported an increase in BRS in healthy participants [21] or patients with chronic obstructive pulmonary disease (COPD) [22], but none measured ventilatory changes or BRS responses to hyperoxia in diabetes. Inhalation of oxygen may be a simple intervention to test the hypothesis of a resting hypoxic condition in diabetes, and test the efficiency of the autonomic function with a simple clinical test.

Methods

Participants

We studied 96 participants with type 1 diabetes and 40 age-matched healthy controls. The participants were recruited through the register of The Social Insurance Institution that comprises all patients entitled to special reimbursement of insulin or antidiabetic medication in Finland. Selection criteria were diabetes (E10 in ICD-10; www.who.int/classifications/icd/en/) diagnosed before the age of 35 years, and age 18–35 years at the time of inclusion.

Type 1 diabetes was defined as C-peptide deficiency (<0.03 nmol/l) and initiation of permanent insulin treatment within 1 year after the diagnosis of diabetes. None of the patients showed clinical signs of cardiovascular disease. However, six patients were given laser treatment because of diabetic retinopathy.

Nineteen patients received antihypertensive medication (14 with an ACE inhibitor, one with a combination of an ACE inhibitor, a calcium channel blocker and a diuretic, one with a combination of an ACE inhibitor and a calcium channel blocker, and three with angiotensin-2 receptor blockers). The healthy control participants were recruited by e-mail advertisements among university students and staff. Only individuals with a normal fasting glucose level and without first-degree relatives with diabetes mellitus were included. Before participation, all participants gave their written informed consent. The study protocol was approved by the Ethics Committee of Helsinki University Hospital, and the study was carried out in accordance with the principles of the Declaration of Helsinki as revised in 2000.

The participants underwent a clinical examination, resting ECG, laboratory testing, overnight urine collections, and standard autonomic function evaluation by four cardiovascular tests: the expiration/inspiration ratio of RR interval during slow deep breathing, the maximum/minimum 30:15 ratio of RR interval during a Valsalva manoeuvre and active standing, and the systolic blood pressure (SBP) response to standing. Cardiovascular autonomic neuropathy was defined as the presence of two or more abnormal tests [26].

Each participant completed a detailed questionnaire on lifestyle, smoking habits and family history.

Protocol

All participants were investigated in a quiet room, at a temperature between 19°C and 23°C, between 08:00 and 14:00 hours. The participants received instructions to refrain from alcohol for 36 h, and caffeinated beverages and cigarettes for 12 h prior to the examination. A light meal was permitted 2 h before testing. If a participant reported or measured symptoms or values of hypoglycaemia in the previous 24 h, the test was postponed. The ECG was recorded using a bipolar precordial lead. Continuous blood pressure was monitored with a Finapres 2300 digital plethysmograph (Ohmeda, Louisville, CO, USA) from the middle finger of the right arm held at heart level. Two respiratory signals were obtained by inductive plethysmography (Z-rip; Pro-Tech, Mukilteo, WA, USA), from belts positioned around the chest and the abdomen. Pulse oximetry and expired carbon dioxide partial pressure (Cosmo; Novametrix, Wallingford, CT, USA) were also obtained.

The signals were simultaneously recorded in the supine position during 5 min spontaneous breathing, during 2 min controlled breathing at a frequency similar to normal spontaneous breathing (15 breaths/min), and during 2 min slow deep breathing at the rate of six cycles/min. Subsequently, the participants repeated the entire protocol while breathing 5 l/min oxygen. Signal recordings started after the first 5 min of oxygen administration to allow stabilisation of oxygen saturation and ventilation. The sequence of breathing rate (spontaneous, 15/min and 6/min) was randomised within each session (normoxia and hyperoxia). All signals were simultaneously acquired on a personal computer with an analogue-to-digital converter with a 12-bit resolution at a sampling rate of 200 Hz (WinAcq data acquisition system; Absolute Aliens Ltd, Turku, Finland).

Assessment of BRS

From the original data, the time series of RR interval (from each of two consecutive R waves of the ECG) and SBP were obtained. Previous studies have shown a poor correlation between different indices of BRS, while, on the other hand, no method has shown clear superior performance over the other [27]. Accordingly, we computed a set of six different tests and used their average [7, 28].

BRS was determined from spontaneous fluctuations in the RR interval and SBP during the spontaneous, 15/min and 6/min recordings using the positive and negative sequence methods [29], the alpha coefficient in the low- and high-frequency bands and its average [30], and the transfer function technique [31]. In the sequence methods, BRS was estimated by identifying spontaneously occurring sequences of three or more consecutive heartbeats in which both the SBP and the subsequent RR intervals changed in the same direction. The minimum criteria for change were 1 mmHg for SBP and 5 ms for the RR intervals. For identified positive and negative sequences with a correlation coefficient between the RR intervals and the SBP exceeding 0.85, the regression slopes (the slope of the regression line between SBP and RR intervals) were calculated, and the average was taken as a measure of BRS positive and negative slopes, respectively. The other four BRS methods were calculated by autoregressive uni- and bivariate spectral analysis [5]. The alpha coefficient was calculated as the square root of the ratio of the powers of RR intervals and SBP in the low frequency range (0.04–0.15 Hz) and in the respiratory (0.15–0.40 Hz) high-frequency range when coherence was greater than 0.5 and the phase difference between the SBP and RR intervals was negative. In the transfer function method, BRS was calculated as the average value of SBP–RR cross-spectrum divided by the SBP spectrum in the low-frequency range (0.04–0.15 Hz), when coherence exceeded 0.5.

Analysis of respiration

The signals from the inductive plethysmographic belt signals were analysed by an automated program with interactive supervision, to identify for each breath the positive and negative respiratory peaks, together with the respiratory period. The sum of the signals obtained by the two belts was taken as a relative index of tidal volume. The same program automatically identified the end-expiratory (end-tidal) value in the carbon dioxide signal. Using the inductive belt data, a semi-quantitative intra-subject analysis of ventilation could be obtained, by comparing the relative changes in tidal volume and minute ventilation induced by oxygen inhalation or different breathing patterns. Although the device used for the present study does not allow us to obtain tidal volume and minute ventilation in absolute values (ml and l/min, respectively), we took advantage of the strong linear relationship between tidal volume and the inductive belt signals [32], enabling us to obtain ventilation in relative units. The limitation of the semi-quantitative analysis is compensated by the lack of interference with the spontaneous breathing, typical of mouthpieces [33]. We therefore set the minute ventilation obtained during spontaneous breathing in normoxia (our baseline) as 100% in each participant, and calculated the minute ventilation or tidal volume in percentage changes from that value for each recording.

Laboratory tests

Venous blood samples were obtained after a light breakfast and were analysed for HbA1c, lipids and serum creatinine. HbA1c concentrations were determined by an immunoturbidimetric immunoassay (Medix Biochemica, Kauniainen, Finland). Serum lipids (cholesterol, triacylglycerols and HDL-cholesterol) and creatinine were measured by enzymatic methods. Urinary AER was measured from three consecutive timed urine collections, one 24 h and two overnight collections. Normal AER was defined as values persistently <20 μg/min or <30 mg/24 h, microalbuminuria as AER ≥20 <200 μg/min or ≥30 <300 mg/24 h, and macroalbuminuria as AER ≥200 μg/min or ≥300 mg/24 h in at least two out of three urine collections [34].

Statistical analyses

Data are presented as mean ± 95% CI, unless differently stated. A global index of heart rate variability was assessed using the standard deviation of the RR interval (SDNN), as this variable has a more normal distribution as compared with other indices of variability (e.g. variance). Differences between the two groups and between conditions were expressed as mean ± 95% CI. Statistical differences between participants with type 1 diabetes and healthy control participants and interventions (6/min vs 15/min controlled breathing and oxygen) were tested using a linear normal model. As outcome, we modelled the different continuous variables (BRS and blood pressures), and included conditions (normoxia/hyperoxia), breathing patterns (spontaneous, controlled at 15 breaths/min and 6 breaths/min) and participant group (healthy control/diabetic) as categorical covariates. For each continuous variable, we assessed the interactions between conditions, groups and breathing patterns [35]. The Sheffé test was used to test for significances between different breathing rates. All tests were repeated after adjusting for the effect of age. Linear regression analysis was used to test for association between variables. Statistical significance was defined as a p value ≤0.05.

Results

The clinical characteristics of the participants studied are shown in Table 1. The diabetic participants studied had reasonable glycaemic control, one participant had macroalbuminuria, six had microalbuminuria, 19 had one abnormal autonomic test (orthostatic hypotension was abnormal in two participants), and deep breathing was also reduced in five controls. Average values of mean RR interval, systolic and diastolic blood pressure, and oxygen saturation are shown in Table 2 with the oxygen-induced differences and their estimates. The effects of oxygen on BRS are shown in Fig. 1.

Table 1 Clinical characteristics and laboratory measurements of participants with type 1 diabetes and healthy control participants (mean ± SEM)
Table 2 Effects of oxygen and respiratory patterns on cardiorespiratory variables (mean and 95% CI)
Fig. 1
figure 1

Effect of oxygen and slow breathing on BRS. Data are expressed as means ± SEM. White circles, spontaneous breathing; black circles, slow breathing. *p < 0.05, **p < 0.01, ***p < 0.001. Under spontaneous breathing in ambient conditions (normoxia), BRS is depressed in individuals with type 1 diabetes. However, in hyperoxia, type 1 diabetic participants increase BRS more than controls, and the difference disappears. Slow breathing increases BRS in normoxia in both groups to an extent similar to hyperoxia, and eliminates the difference between diabetic and control participants. Hyperoxia did not further increase the BRS in control participants, indicating that the two effects are probably related

Resting data and effect of respiratory manoeuvres

As expected, the patients with diabetes had lower (p < 0.05) BRS, and a trend to a reduced RR interval and heart rate variability, compared with the control participants. In addition, diabetic participants had higher resting systolic and diastolic blood pressures (p < 0.05). All other variables, including resting oxygen saturation, were similar to controls (Table 2).

During slow breathing, BRS increased more in diabetic than in control participants, and as a result the statistical difference in BRS present at rest disappeared (Fig. 1). Controlled breathing at 15 breaths/min induced a significant increase in heart rate (a reduction in RR interval), a reduction in heart rate variability and also a significant reduction in BRS (from 14.8 ± 0.8 to 13.5 ± 0.9 ms/mmHg in diabetic participants, and from 18.5 ± 1.4 to 15.3 ± 1.2 in controls; p < 0.05 for both). Controlled breathing, and particularly slow breathing, increased oxygen saturation (p < 0.001) in normoxia, but while controlled breathing at 15 breaths/min nearly doubled resting ventilation in both groups (p < 0.001), slow breathing induced only a modest (+10%), non-significant increase in ventilation (Table 2).

Effect of oxygen inhalation

Inhalation of 5 l/min oxygen increased oxygen saturation (p < 0.001) to nearly 99% on average (Table 2) during spontaneous breathing in both groups. The effect was attenuated but still significant during controlled breathing at 15 breaths/min (p < 0.01 in diabetic and p < 0.001 in control participants) and during slow breathing. The attenuation may be explained by the increased level of oxygen saturation due to the controlled or slow breathing per se in normoxia (Table 2).

Oxygen inhalation reduced heart rate (increase in RR interval), and increased heart rate variability (SDNN). Oxygen inhalation increased systolic and diastolic blood pressure and BRS during spontaneous breathing (Table 2), with a similar trend in diabetic and control participants. However, the extent of the change was greater and significant only in the diabetic individuals (Table 2), with a significant (p < 0.05 or better) interaction effect (groups × interventions × breathing) for systolic and diastolic blood pressure. For BRS, there was a significant interaction effect between interventions and groups during spontaneous breathing (Table 2). As a consequence, the significant differences (reported with their confidence limits) observed at baseline between diabetic and control participants increased during oxygen administration (p < 0.001 for both systolic and diastolic blood pressures), or decreased for BRS (p = NS; Table 2 and Fig. 1) during spontaneous breathing. The differences with their estimates can be seen in Table 2. It clearly appears that the behaviour of blood pressure and BRS differed between the diabetic and control participants, as shown by a significant interaction. The BRS values obtained during spontaneous breathing approximately predicted the BRS increase induced by oxygen administration (r = −0.216, p = 0.012; BRS increase from baseline to oxygen = 4.77–0.175 × resting BRS), and showed a reverse effect for higher BRS values; that is, for values >28 ms/mmHg (calculated with the present methodology), oxygen administration would have reduced BRS. This lack of increase with oxygen in fact occurred during slow breathing, when the BRS reached the highest values already in normoxia (Fig. 1). During controlled breathing at 15 breaths/min, the increase in BRS was still present but did not reach statistical significance.

Oxygen administration increased tidal volume during spontaneous breathing in diabetic but not in control participants. However, there was no significant interaction. There was also a trend towards a relative increase in ventilation that was more evident in the diabetic participants. This trend was also confirmed by a reduction in end-tidal carbon dioxide. During controlled breathing at 15 breaths/min, oxygen administration induced only minor changes in ventilatory variables (Table 2). Adjusting for age did not modify the results.

The BRS response to oxygen was not influenced by glycaemic control or albuminuria. All these results remained unchanged even after removing the 19 participants treated with drugs affecting angiotensin II.

Discussion

These results show that the autonomic imbalance and particularly the reduction in BRS of patients with type 1 diabetes can be partially reversed by oxygen, suggesting a possible role of tissue hypoxia.

This observation does not exclude several other factors, such as insulin treatment, hyperglycaemia, antecedent hypoglycaemia, cortisol and poor glycaemic control, from being also relevant. The presence of other possible determinants could explain why the correction remained partial after oxygen administration. However, hypoxia appears to be an important factor that restricts the BRS and modifies the autonomic function of these patients. This finding suggests that the autonomic dysfunction might be part of a more general modification of autonomic reflexes, possibly involving an altered control of ventilation. Previous studies in diabetes showed abnormalities not only in the cardiac vagal regulation, but also in the control of blood pressure [36, 37], microcirculation [38] and respiration [39, 40]. Respiratory and cardiovascular control are tightly intertwined [41]. Accordingly, any modification of respiratory control will influence cardiovascular control, and vice versa, if the autonomic abnormalities are functional. Conversely, in neural damage, this interaction is minimal or non-existent, as has been shown in patients with severe autonomic neuropathies [42]. Our findings of an important improvement in BRS with oxygen in type 1 diabetes supports our previous findings [7] of a predominantly functional disorder in patients without complications.

Effect of oxygen inhalation

Oxygen administration reduces heart rate and increases indices of heart rate variability related to parasympathetic activity [20, 23]. Two studies [21, 22] have also reported an increase in BRS in healthy participants and patients with COPD. These findings were interpreted as if oxygen inhalation reduces the chemoreflex activity, even in normoxia. This in turn would increase BRS, as these reflex systems have reciprocal effects. Activation of chemoreflexes normally induces sympathetic activation, an increase in ventilation and a depression of BRS (whereas baroreflex activation is associated with parasympathetic stimulation [41] and opposite effects on ventilation). However, there was no demonstration of any such phenomena in the previous studies on oxygen administration, as ventilation (or chemoreflex) was not evaluated [21, 22]. In our study, we did not measure ventilation or the chemoreflex in absolute values, as the use of a mouthpiece or a face mask would have altered the spontaneous ventilation [33]. However, hyperoxia showed a definite trend towards increased ventilation during spontaneous breathing. This well-documented increase occurs after a transient reduction [24, 25, 43], and has been ascribed to a vagal reflex [43, 44], in connection with an increased generation of oxygen free radicals. There was also some evidence of a more pronounced effect in diabetic participants, probably as a consequence of a greater vagal stimulation induced by oxygen. Unfortunately, the lack of a precise definition of resting ventilation precluded more precise assessment, although the carbon dioxide data confirmed that the increase in ventilation was, at least, as large as in controls. All these findings contrast with the concept of predominant neural damage (that would invariably blunt responses to parasympathetic stimuli) and again suggests an important functional abnormality. Higher pre-existing hypoxia and consequent endothelium-related vascular changes are a more logical explanation for the response in patients with type 1 diabetes.

Another interesting and novel finding of the present study is the marked oxygen-induced increase in blood pressure seen in diabetic individuals. Oxygen administration increases blood pressure, probably by a direct local effect. This might have stimulated the BRS and parasympathetic activity. This response is reciprocal to the well-known effect of hypoxia (direct vasodilation, and a reflex increase in sympathetic activity and reduction in BRS). If this were the case, the stronger effect of oxygen seen in type 1 diabetes could again be suggestive of resting tissue hypoxia. An alternative possibility could be an effect of reactive oxygen species on the dysfunctional endothelium typical of diabetes. However, previous findings [45] (and the notion also that hypoxia increases free radicals but causes arterial vasodilation) do not completely support this possibility, which remains open for future research.

Sun et al. [23] were the first and to our knowledge the only group that has tested the effect of hyperoxia in diabetes. They found an increase in heart rate variability in diabetic patients after 4 weeks of hyperbaric hyperoxia. To the extent that their study and our present study can be compared, the results are in complete agreement, but we now provide the first evidence of an improvement in BRS in type 1 diabetes in response to oxygen inhalation. We therefore suggest that the common observation of a low BRS in diabetes may be explained in part by resting tissue hypoxia. The results did not change after removing the participants treated with drugs affecting angiotensin II, confirming no relevant confounding effects of these drugs on BRS [46].

Interaction of oxygen inhalation and respiratory pattern

As expected from our previous studies, slow breathing increased BRS in uncomplicated type 1 diabetes [7] and improved oxygen saturation even in normoxia, similar to what has been shown in heart failure [47] or hypoxia-dependent diseases [48]. It is possible that the effect of slow breathing in increasing parasympathetic activity and BRS could be mediated by an increase in arterial oxygen pressure, as evidenced by the increased oxygen saturation. Accordingly, the increase in BRS seen during slow breathing was not further augmented by oxygen (Fig. 1). This effect could also have been predicted by the relationship between BRS during spontaneous breathing in normoxia and its changes with oxygen. This suggests that after a maximal value is reached, BRS cannot be further increased (a ‘saturation’ effect). It is also likely that the high BRS observed during slow breathing prevented a further increase in blood pressure as a result of oxygen in control participants, whereas it could not fully block the increase in diabetic patients, possibly as a consequence of tissue hypoxia and lower resting values.

The response to oxygen was markedly reduced when breathing was paced at faster respiratory rate (15/min). This is likely to be due to the very large increase in ventilation induced by paced breathing, and the associated increase in sympathetic activity typically occurring with controlled breathing at faster rates. As such, the increase in ventilation due to voluntary control could have counteracted/blunted the effects oxygen administration.

Implications and conclusions

This study has theoretical and practical implications. Diabetic patients are thought to have autonomic abnormalities due to neural damage. In a previous study [7], the low BRS of type 1 diabetic participants was partially restored in the short term by a simple functional manoeuvre (slow breathing), suggesting a functional component. We now provide evidence suggesting that tissue hypoxia may be responsible, in part, for the autonomic dysfunction observed in type 1 diabetes. Hypoxia per se is a potent stimulator of the sympathetic nervous system and of ventilation, which should in turn correct the hypoxia. However, our results indicate that hypoxia persists despite the high sympathetic activity present at rest, and are also suggestive of a disorder in respiratory control.

The importance of uncorrected hypoxia is now being recognised as a key factor in the origin and development of diabetic complications [13, 14], and insufficient compensatory responses to long-term hypoxia have been found in experimental diabetes [49]. Accordingly, subclinical hypoxia could be detected by an improvement in autonomic function after a simple clinical manoeuvre like oxygen administration. Our findings confirm that hyperoxia is not simply a suppression of the chemoreflex stimulus, but provides an additional parasympathetic stimulus, capable of increasing ventilation. The presence of hypoxia, the possible coexistence of abnormal respiratory control and the functional component of these abnormalities highlight the possibility of correction using simple strategies like physical training. Physical training has been shown to improve both respiratory and cardiovascular control in patients with heart failure [50], and is also known to improve cardiovascular and autonomic nervous function in type 2 diabetes [3, 4].

In conclusion, resting tissue hypoxia appears as one of several possible causes of functional autonomic abnormalities in patients with type 1 diabetes. These findings may have relevant implications for the understanding of the origin of diabetic complications, as hypoxia may in fact represent a possible link between the autonomic and the other major diabetic complications. Our findings support interventions such as physical activity, and thus may be of practical help in the prevention and management of diabetic complications in patients with type 1 diabetes.