| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1Laboratoire de Physiologie Respiratoire Expérimentale, Théorique et Appliquée, Université Joseph Fourier, 38700 La Tronche, France; 2Department of Physical Education, McGill University, Montreal, Quebec, Canada H2W 1S4; and 3Laboratoire de Modélisation et Calcul, Université Joseph Fourier, 38041 Grenoble, France
Submitted 10 July 2003 ; accepted in final form 26 January 2004
 |
ABSTRACT |
|---|
 TOP ABSTRACT METHODSRESULTSDISCUSSIONREFERENCES |
|---|
human; heart rate variability; autonomic control of heart rate; individuality of breathing pattern
In addition to this "multidimensional" aspect to parasympathetic cardiac control, the between-individual differences in RSA increase the difficulty of RSA interpretation. In the latest committee report on heart rate (HR) variability (HRV), method, and interpretation (3), this between-individual aspect is mentioned several times leading to the conclusion that "caution needs to be exercised in interpreting RSA, especially for between-subject comparisons." For example, the relationship between RSA amplitude and pharmacologically defined vagal tone when investigated for a number of subjects appears to be less close than that found within a given subject (15, 20). The increase in RSA amplitude with increasing breathing period also varies among individuals (4, 17).
Given the importance of the respiratory modulation of human autonomic rhythms and the prepotency of breathing in generating respiratory frequency rhythms (1), the incidence of the spontaneous breathing rate is worth clarifying because there exists an individuality of the breathing pattern (2, 27) that has been reported in several physiological conditions (12, 26). We investigated the role of spontaneous breathing period in RSA changes induced by breathing period changes. To quantify RSA, we used a breath-by-breath analysis where a sinusoid is fitted to the instantaneous HR for each breath (23). The amplitude and the phase (or delay) of this sinusoid constitute the characteristics of RSA for that breath, and mean values are calculated for each breathing period. We thus quantified RSA in healthy subjects at various breathing rates.
Our hypothesis was that the rate of increase in RSA amplitude with respiratory period might be similar in all subjects, provided a comparable range of breathing rates is examined. We defined this comparable range by starting from the individual spontaneous respiratory frequency of each subject and by choosing breathing rates surrounding this spontaneous breathing rate. Our results show that, even under these conditions, a difference in the rate of increase in RSA amplitude with the breathing period existed between individuals. Similarly, difference between individuals was found in the parameters of the parabola fitted to the delay-breathing period relationship. However, these differences (the rate of increase in RSA on one hand and the curvature of RSA on the other hand) were correlated with the spontaneous breathing period and not with the mean R-R interval.
These results suggest that in all individuals there are changes in RSA characteristics with total respiratory period (TTOT), but there are differences between individuals in the features of these changes. These differences appear to be related to the spontaneous breathing period of the subjects.
In addition, this protocol gives the possibility of performing a transfer function analysis of the RSA control system (25). Indeed, the plot of the normalized RSA [amplitude of RSA/tidal volume (VT)] versus TTOT provides the "gain" of the system at each TTOT for each individual.
|
METHODS |
|---|
|
TOPABSTRACTMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
Twelve healthy volunteers between 18 and 28 yr of age, seven of whom were men, participated in the study (means ± SD: height, 173.5 ± 7.8 cm; weight, 64.7 ± 10.4 kg; age, mean 22.9 ± 3.1 yr). Informed consent was obtained from all subjects. The experimental protocol was examined and approved by the Institutional Ethics Review Board.
Experimental Protocol
Subjects were comfortably seated and wore a face mask on which a flowmeter (Fleish head no. 1) and a differential pressure transducer (163PC01D36, Micro Switch) were mounted. Leaks from around the mask were checked for before the recording was initiated using an infrared CO2 analyzer (Engström Eliza/Eliza MC). End-tidal CO2 fraction (FETCO2) was measured continuously using the same apparatus, and an electrocardiographic trace (ECG) was obtained covering the whole of the recording period. For each subject, six series of 5- to 10-min recordings were performed, with the first corresponding to spontaneous breathing and the following five with randomly sorted sequences: at an imposed frequency fixed at the mean spontaneous respiratory frequency and at the mean spontaneous respiratory frequency +3 and +6 and –3 and –6 breaths/min. For subjects with a low spontaneous respiratory frequency, the –3 and –6 breaths/min recordings were replaced by –2 and –4 breaths/min. The highest breathing rate observed among our subjects was 17 breaths/min, so recordings for that subject were performed for 11, 14, 17, 21, and 23 breaths/min, whereas the lowest breathing rate observed was 7 breaths/min with recordings accordingly performed at 3, 5, 7, 9, and 11 breaths/min.
To impose the breathing rate, an auditory cue was used, which signaled only for the inspiration to begin. Hence, the inspiratory II-to-expiratory time ratio, as well as the VT, was chosen by the subject. However, if FETCO2 departed more than ±0.4% from the control level, subjects were requested to change their VT accordingly. In fact, in most cases, during imposed frequency breathing, subjects are inclined to hyperventilate and thus before the recording was started, they were asked to decrease the VT.
Data Acquisition
The acquisition of the data was carried out on a Macintosh microcomputer equipped with an analog-digital interface card. The sampling rate was 256 Hz. To calculate TTOT and VT, and to study HRV and RSA, a breath-by-breath analysis of all recordings was performed. The ECG signal was processed, and the R-R interval series were extracted and displayed on the computer screen to verify that the signal exhibited no noticeable trend and to show up possible errors. Means and SD of the R-R intervals were calculated for each recording. R-R intervals were interpolated linearly at 0.25-s intervals to obtain equidistant time samples, and spectral analysis was performed using a recording length of at least 1,024 sample data points.
Data Analysis
A Fourier transform procedure was applied to obtain the low-frequency (LF; 0.04–0.15 Hz) and high-frequency (HF; 0.15–0.40 Hz) components. For each recording, a restricted respiratory frequency component identified as the respiratory centered frequency (RCF) component was also calculated using the frequency range corresponding to ±10% of the respiratory rate averaged over the entire recording (23). The power corresponding to the different frequency ranges was expressed as a percentage of the total spectral power minus that corresponding to the very LF (0–0.04 Hz).
A more specific analysis of RSA was performed using a breath-by-breath HRV analysis (23). To quantify the extent of within-respiratory cycle HRV, a sinusoid was fitted to the changes in instantaneous HR within the respiratory cycle. For each breath, the maximum value of the sinusoid was expressed as a percentage of the mean cardiac frequency calculated over that breath. This maximum value is used as a measure of RSA amplitude. The time elapsing between the beginning of the breath and the occurrence of this maximum value is expressed either in terms of the fraction of breath duration (phase) at which it occurs or in seconds (delay = phase x TTOT). Average amplitude, phase, and delay values over several breaths were then calculated for each recording. In addition, for each breath, a normalized RSA amplitude (RSA amplitude divided by the corresponding VT) was calculated.
Statistical Analyses
Values are expressed as means ± SD. Mean values of the various variables between spontaneous and imposed breathing at the same rate were compared using a paired t-test, whereas SDs were compared using a Wilcoxon paired test.
To compare HRs between different recordings for a given subject, mean R-R intervals at different imposed breathing frequencies were compared using ANOVA.
Linear correlation coefficients and regression lines were calculated for each subject from amplitude against breathing period plots using all the values available for this subject (250–300 breaths). To test the hypothesis of a parallel regression line for all subjects, a linear regression line was also calculated using all the values of all subjects. The differences in slope between this common line slope and the individual slopes were calculated, and the sum of the weighted differences was compared using a 
2-test with 11 degrees of freedom as it applies to 12 subjects.
Parabolic fit was calculated for each subject on delay-vs.-breathing period plots. Also, one parabolic shape was adjusted using all available values of all subjects. As in the case of the regression lines, to test the hypothesis of the existence of a common parabola, the sum of the weighted differences between the curvature parameter of the common parabolic shape and those of the individual parabolas were calculated again using a 2-test.
For all tests, significance was set at P < 0.05.
|
RESULTS |
|---|
|
TOPABSTRACTMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
|
Figure 1 shows an example of two recordings on one subject when breathing spontaneously (left) and when breathing at an imposed frequency equal to the spontaneous breathing rate (right). Figure 1 also shows the instantaneous HR, delimited breath by breath, and the corresponding spectral analysis for the whole recording.
|
|
|
RSA amplitude and delay were also compared between the two conditions (Fig. 3, B and C), and there was no significant difference with regard to either mean values or SDs.
RSA Changes for Imposed Breathing Frequencies
Spectral analysis.
Figure 4 shows the RCF and HF components of all the imposed breathing frequency recordings for the 12 subjects. As expected, the HF component falls for respiratory periods longer than 6.78 s (0.15 Hz), whereas the RCF component exhibits a plateau at 
80% of the total power and this plateau is reached for values of TTOT of between 6 and 7 s.
|
2-test. The value of the test was 2 = 698, and thus the hypothesis was rejected (P < 0.001).
|
|
2-test to explore the hypothesis of a parabola common to all subject lead to a rejection of the hypothesis for all the parabola curvature (2 = 93.5, P < 0.01).
|
|
|
|
DISCUSSION |
|---|
|
TOPABSTRACTMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
Spontaneous Versus Imposed Breathing
Before studying RSA corresponding to different imposed breathing periods, we felt that the RSA for spontaneous and imposed breathing at the same rate should be compared to ascertain that the respiratory period was the main explicative factor and thus to justify the study being carried out with imposed breathing rates.
In humans, volitional breathing arises from a corticomotor excitation of the diaphragm, which may act directly on the phrenic motor nucleus via the corticospinal tract "bypassing" brain stem respiratory centers; alternatively, or in addition, this may be achieved indirectly via the respiratory centers and bulbospinal paths. There is no clear evidence for these alternatives, although several studies suggest that the corticomotor excitation of the diaphragm does not transit via the medullary respiratory centers (7). The observation of no difference in RSA between spontaneous and imposed breathing at the same rate does not provide indication in favor of either alternative. It merely suggests that for each individual breathing rate per se has an influence on the characteristics of RSA. HRV and RSA, spontaneous or imposed, at the same breathing rate have been studied by several authors (10, 16, 17, 21, 22). In these studies, no difference was found in either mean HR or HF amplitude between spontaneous and paced, or metronomic, breathing. Our results in agreement with these observations lead to the conclusion that the voluntary control of breathing does not enhance vagal tone or alter the vagal modulation of HR.
Experimental Protocol
We chose to impose breathing frequencies starting off from the spontaneous breathing rate of each subject. This was done because of the great variability in resting breathing rate among subjects (8) and also because this breathing rate appears to be an individual characteristic (2). Thus the recordings were not performed over the same range of breathing frequency for all 12 subjects, and also the width of the frequency range used varied with the individual subject. This can be seen clearly in Fig. 5, where the extremities of a regression line show the extreme values of respiratory period obtained for the subject concerned. These differences may be considered as a bias in the experimentation, but on the other hand it should be pointed out that a frequency of 12 breaths/min represents an increase in breathing rate for some subjects and a decrease for others, and this may influence the regulatory mechanisms brought into play. In this study, for each subject, starting from a breathing rate corresponding to the spontaneous rate two higher and two lower rates were imposed, thus surrounding the spontaneous rate.
R-R Interval: Mean Value and Spectrum for Different Breathing Period
Changing the respiratory period did not change the HR. This observation is in accordance with several authors (1, 4, 17), and, as suggested by Brown et al. (4), different within-subject breathing frequencies and depths distribute vagal firing within the respiratory cycle but do not alter the level of vagal outflow. Kollai and Mizsei (20) found that the mean heart period changed in response to slow breathing and that the nature of the change (increase, decrease, or no change) varied with the individual. Differences in the experimental protocol and in the range of periods explored in each subject may explain this divergence.
The percentages of the power spectrum occurring in the HF and RCF bands of the instantaneous cardiac rate signal are shown in Fig. 4. It can be seen that, as expected, the power in the HF band is greater than power in the RCF band up to the period (6.7 s) corresponding to 0.15 Hz and that beyond this there is a noticeable difference between the two plots: the values of HF fall, whereas a plateau is reached for RCF, which appears to be at 80% of the total power spectrum. It should be noted that if the HF band is not suitable as a means of describing the power corresponding to breathing at low breathing rates, the RCF band has the disadvantage of varying in width with breathing rate.
Nevertheless, as the total surface of the power spectrum represents the variance of the analyzed signal, the percentage of the power in a given band (here, RCF) may be considered to be the variance associated with this frequency range. Thus the fact that the power in the RCF band reaches a plateau suggests that the extent of the HRV dependent on breathing rate does not exceed 80% of the total variability.
RSA Amplitude
All subjects exhibited an increase in RSA amplitude with increasing respiratory period at least in the explored range of periods. This observation is in accordance with the findings of Hirsch and Bishop (17), who estimated RSA amplitude as the difference between the lowest and highest instantaneous HRs in each breath, and also those of other authors who used spectral analysis (10).
Our initial hypothesis, which was that parallel regression lines could be found, representing the RSA amplitude-versus-breathing period relationship for all subjects, was rejected by the statistical test, which indicated interindividual differences in this relationship. Similar results have been reported by Hirsch and Bishop (17), who plotted RSA amplitude versus breathing rate for each subject on a log-log scale. They found a constant RSA for low breathing rates, below 3–7 breaths/min, and then a decreasing relationship, the slope of which was expressed in decibels per decade and which defined the system roll-off. The origin of the decreasing relationship (LF intercept) varied among subjects, as did the roll-off, i.e., the slope of the log RSA-log breathing frequency plot.
In the same way as for the roll-off, the slope of RSA amplitude-breathing period relationship exhibited interindividual differences and this slope was not related to an individual's HR but correlated to their spontaneous breathing period. This slope, which represents the rate of change of RSA amplitude with TTOT, may be considered to be a measure of the responsiveness of the HR modulation mechanism to changes in respiratory period. The significant negative correlation with the spontaneous TTOT suggests that subjects with a low spontaneous breathing period (high rate) will be more responsive; i.e., their RSA amplitude will increase more with increasing period than for subjects with long breathing periods.
This responsiveness is somewhat different from the one defined by Kollai and Miszei (20). They defined an individual RSA responsiveness expressed as the ratio 
HP/RSA, where HP is the change in heart period with increasing respiratory period. Because three types of HP were found [individuals with increasing (A), decreasing (Z), and unchanged (A/Z) HP], the slope of the HP-RSA plot was respectively positive, negative, or close to zero. This ratio was found to be correlated to parasympathetic control (PC), defined as the changes in heart period after complete parasympathetic blockade by the administration of atropine. Although there was a continuum of distribution of subjects along the HP/RSA-against-PC regression line, the A type was mainly distributed on the left-hand side of the graph (i.e., high responsiveness associated with low PC) and the Z type was mainly distributed at the right-hand end of the regression line. Kollai and Miszei concluded that the interindividual differences in RSA have their origin mainly in the differences in PC, although they found that introducing respiratory characteristics (TTOT and VT) improved the degree of RSA-PC correlation.
RSA Delay
The delay is the time elapsing in each breath between the onset of inspiration and the reaching of the maximal value of the sinusoid fitted to the instantaneous HR. We chose the delay (delay = phase x TTOT) rather than the phase because of its time dimension, which might make the interpretation of the results easier. This variable differs somewhat from the phase angle between the P-P interval and respiration reported by Eckberg (9), who measured the phase angle between the onset of inspiration and the heart period shortening. Although both Eckberg's study and the present one are concerned with the time lapse between inspiratory onset and HR acceleration, for Eckberg the phase was defined by the start of the acceleration, whereas in this study the delay was defined by the moment of the maximum acceleration. Nevertheless, in both studies, a change occurred at periods in the range of 8–10 s. Indeed, Eckberg reported that at breathing intervals of 8 and 10 s, P-P shortening began before the onset of inspiration, so that the polarity of the phase changed from positive to negative as the breathing period increased. In our study, the delay-TTOT relationship is parabolic and for 8 of 12 subject the maximum occurred in the 8- to 10-s range (Table 1).
Our hypothesis was that of the existence of a common parabolic shape representing the RSA delay-TTOT relationship for all subjects. This was based on the assumption that the time taken to develop the inspiratory parasympathetic inhibition on the heart and reach the maximum HR acceleration would be similar among subjects. Given our protocol, we expected a parabolic fit only for those subjects with long TTOT. Our results rejected this hypothesis, and, as can be seen in Fig. 7, even for those subjects such as subjects 1 and 2 with short spontaneous TTOT, a parabolic relationship exists between the delay and TTOT, and the delay is reached earlier than at a TTOT value of 8 s. The curvature as well as the other parameters of the parabolas were found to be related to the spontaneous TTOT, indicating the influence of the latter on changes in RSA brought about by changing the breathing period.
One possible interpretation of these results would be that the pattern of the modulatory mechanism is affected, or even determined, by the spontaneous breathing period so that for subjects with a long spontaneous TTOT the modulation will be initiated either later in the breathing period or much more gradually than for those subjects with a shorter spontaneous breathing period.
RSA Transfer Function Analysis
In this study, it can be assumed that in each series of recordings on one subject, the only variable that is changed is the breathing rate. Thus the breathing rate changes may be considered as an input to the HR modulation system, the output of which is RSA. As any change in the breathing rate is associated with changes in VT to compare the gain of the HR modulation system within and between subjects, the amplitude of RSA of each breath has to be normalized, i.e., divided by its VT.
Figure 9 shows that the change in the gain of the system varies among individuals; in some subjects (subjects 7, 9, and 10) there were little changes, in contrast to some others (subjects 1, 3, and 12). Furthermore, over a fixed respiratory period range, there are between-individual differences in the gain of the system. This indicates that these differences are not due to the different amount of the input but rather to the characteristics of the system.
Thus if one subject's line is "higher" than another subject's line, then the former subject has greater gain. If one subject's line intersects with another subject's line, then the former subject has greater gain over a particular range of respiratory period and lesser gain over the remaining respiratory range.
Transfer function analyses have been used in physiological (25) and pathological (13) conditions where the gain of the system was quantified and considered as being a measure of the autonomic tone.
Individuality of RSA?
These results, in addition to several other observations (4, 15, 20) on interindividual differences in RSA changes, suggest that there may exist an individuality of RSA. However, the existence of individuality implies not only differences between individuals but also reproducibility for a given subject. Little data are available on the reproducibility of RSA at a given period. Grossman et al. (14) and Grossman and Kollai (15) showed that behavioral tasks known to influence cardiac vagal tone produce closely corresponding within-subject changes in mean R-R interval and RSA when the respiratory parameters are controlled. The comparison of spontaneous and metronomic breathing as in this study and several other studies suggests that there was no difference in RSA amplitude between these two conditions even in head-up tilt and low body negative pressure situations, as reported by Patwardhan et al. (22). For respiratory periods other than spontaneous, the RSA amplitude corresponding to an increased respiratory period induced by the addition of resistive load was similar to the RSA amplitude for the same imposed respiratory period (5).
In conclusion, in addition to the differences in RSA between individuals, there exist interindividual differences in the RSA control system response to changes in TTOT dependent on 1) the spontaneous breathing period and 2) on the strength of autonomic tone.
Therefore, these results suggest that, in addition to the influence of respiratory characteristics on the gating of sympathetic and vagal motoneurons responsiveness, the individual breathing rate may play a role in the build up of the PC of HR.
|
ACKNOWLEDGMENTS |
|---|
|
FOOTNOTES |
|---|
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.
|
REFERENCES |
|---|
|
TOPABSTRACTMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |
