|
|
||||||||
Toronto Rehabilitation Institute Sleep Research Laboratory, Department of Medicine, Mount Sinai Hospital, Department of Medicine, Toronto General Hospital/University Health Network, and Centre for Sleep Medicine and Circadian Biology, University of Toronto, Toronto, Ontario, Canada
Submitted 6 January 2005 ; accepted in final form 27 June 2005
| ABSTRACT |
|---|
|
|
|---|
sleep apnea; heart rate variability; atrioventricular node
The VR in AF is not completely random. For example, the VR displays a circadian variation, being higher during the day and lower at night (13). These nonrandom alterations in rate have been attributed to modulation of the electrophysiological properties of the AV node, mediated by the autonomic nervous system (13, 14).
During AF, the VR is determined mainly by the electrophysiological properties of the atrioventricular (AV) node, chiefly the AV refractory period and the degree of concealed AV conduction. The AV refractory period determines the shortest possible time between successive atrial impulses that can be successfully transmitted to the ventricles (36). Concealed AV conduction is the phenomenon whereby atrial impulses incompletely penetrate the AV junction and do not reach the ventricles but affect the transmission of subsequent impulses (16, 42). In turn, changes in AV node refractory period and concealed conduction are considered reflections of alterations in autonomic input to the AV node (13, 14, 24, 45).
Our laboratory has previously shown, in subjects with systolic left ventricular heart failure and sinus rhythm, that Cheyne-Stokes respiration with central sleep apnea (CSR-CSA) entrained HRV at the very low frequency (VLF, <0.05 Hz) of periodic breathing (17, 18, 37). Because sinoatrial rate and AV nodal conduction are modulated by the autonomic nervous system in a similar fashion (14, 24, 45), this finding suggested that CSR-CSA might also be capable of organizing the seemingly random ventricular response to AF by entraining corresponding cyclical changes in AV nodal conduction. We, therefore, hypothesized that this breathing pattern would influence VR in AF at VLF by causing cyclical changes in the electrophysiological properties of the AV node. If so, this would be the first demonstration of respiratory modulation of HRV in the frequency domain in AF. To test this hypothesis, we studied patients with systolic left ventricular heart failure who also had AF and CSR-CSA during sleep. We studied patients with systolic heart failure because they have a high prevalence of both CSR-CSA and AF.
| METHODS |
|---|
|
|
|---|
|
Data acquisition. All data were acquired by a polysomnographic computer program [Sandman Elite; Nellcor Puritan Bennett (46), Ottawa, ON] and analyzed offline by a customized computer program (LabVIEW; National Instruments, Austin, TX). All signals were sampled at 200 Hz, except for the electrocardiograph signal, which was sampled at 1,000 Hz.
Frequency spectral analysis. Spectral analysis is a form of frequency domain analysis that makes use of the fact that component signals can be described by sinusoidal functions of differing amplitude and frequency. The component signals are extracted from the signal being analyzed and used to generate a plot of power (variability) distributed as a function of frequency, referred to as the power spectrum. In the present study, the heart rate [expressed as its inverse, R-wave to R-wave (R-R) interval] and ventilatory signals were analyzed with frequency spectral analysis by using Fourier transforms. This allowed study of the signals and of their interactions at separate frequencies. Segments 13 min 40 s in length were used to provide acceptable frequency resolution at the periodic breathing frequency. Analyses were confined to periods during which there was a stable pattern of breathing (either regular breathing or CSR-CSA) during stage 2 non-rapid eye movement sleep to maximize stationarity of the R-R interval data.
Once data were acquired, they were analyzed over two stages. The first stage (preprocessing) involved extracting the R-R intervals. The R-R intervals and all other signals were then passed through a low-pass digital filter (5th order, cutoff frequency = 0.8 Hz) and resampled at 10 Hz. In addition to allowing reduction of the data set, the low-pass filter suppressed the spectral power at the heartbeat frequency. The second stage of processing involved subdividing the 13-min 40-s data set, composed of 8,192 points, into overlapping segments of 2,048 points each, with half-overlapping of each segment. A Blackman-Harris window was applied to minimize spectral leakage. The cross-power spectrum and autopower spectrum were then ensemble averaged over the seven segments to reduce spectral variance. Power spectra were obtained over the frequency range of 0.0050.8 Hz.
With instantaneous lung volume (ILV) derived from the VT signal taken as the input variable, and R-R interval and BP taken as the output variables, we calculated spectral power, coherence, and phase angles at the respiratory (i.e., high) frequency during periods of regular breathing and at the VLF of periodic breathing during CSR-CSA. Coherence is a measure (from 0 to 1) of the linear correlation between input and output signals. A "significant" coherence implies a strong relationship between two signals at the frequency of interest. The value of 0.5 to indicate significant coherence was first proposed by de Boer et al. (7, 8) and has subsequently been adopted by most authors who use spectral analysis techniques. However, the value is somewhat arbitrary and does not constitute a true test of statistical significance. Accordingly, we used a method first described by Manjarrez et al. (19), wherein a coherence value is considered significant if it exceeds the 95% confidence interval of the coherence function across the entire spectrum from 0.005 to 0.8 Hz. As an internal control, we also examined spectral power and coherence of VR and respiration at VLF (measured at the frequency of the previously observed periodic breathing for each patient) during regular breathing.
Electrophysiological evaluation. The ventricular response to AF is largely determined by the electrophysiological properties of the AV node: specifically, the AV nodal refractory period and the degree of concealed AV conduction. Both of these parameters can be estimated by using information derived from the R-R interval time series (2, 3, 6, 13). To ensure the robustness of our findings, we employed both the original technique of R-R interval analysis described by Billette et al. (3), as well as a separate technique devised by Chishaki et al. (6) and modified by Hayano et al. (13), which employs Lorenz plots to account for the rate dependency of AV nodal conduction. First, the R-R interval series for each patient was divided into sequences occurring during hyperpneic and apneic phases. Wide complex beats were identified visually and excluded from the analysis by two investigators who were blinded to the respiratory signal. Mean R-R intervals and estimates of AV conduction were then compared between apnea and hyperpnea.
The AV nodal refractory period is easily estimated by the minimum R-R interval (3), whereas the scatter (SD) of R-R intervals estimates the degree of concealed conduction (13). Accordingly, minimum R-R interval and R-R interval scatter (corresponding to AV nodal refractory period and degree of concealed AV conduction, respectively) were compared between hyperpnea and apnea. The minimum R-R interval was taken as the fifth percentile of the R-R interval distribution. The fifth percentile value rather than the absolute minimum R-R interval is considered a more reliable estimate of the AV nodal refractory period because it represents more than a single event and has been validated against direct measurements of minimum AV nodal refractory period obtained by rapid atrial pacing (3).
Because both AV nodal refractory period and concealed conduction vary with heart rate, the rate-dependent refractory period and degree of concealed conduction were estimated according to the Lorenz plot technique described by Hayano and coworkers (13). Separate Lorenz plots were generated for R-R interval sequences (2, 6) occurring in the hyperpneic and apneic phases, wherein each R-R interval was plotted on the vertical axis against the immediately preceding R-R interval on the horizontal axis. For each Lorenz plot, the R-R intervals occurring after a particular length of preceding R-R interval have a distribution with a distinct lower limit, corresponding to the rate-dependent AV nodal refractory period. The set of lower limits associated with each length of preceding R-R interval defines a straight line (the lower envelope), above which longer R-R intervals are scattered (Fig. 1).
|
Statistical analysis. Results are reported as means (SD). In a given subject, respiration (ILV) was considered to exert a significant influence on VR (R-R intervals) and BP when the coherence between the two signals exceeded the 95% confidence interval (19) across the frequency range from 0.005 to 0.8 Hz. The prevalences of significant coherences between ILV and HRV at HF during regular breathing and at VLF during both CSR-CSA and regular breathing were compared by Fishers exact tests, and average values of coherence and peak power spectral densities at HF and VLF during each breathing condition were compared within subjects using repeated-measures analysis of variance with post hoc Tukey's test. Mean R-R intervals and estimates of AV conduction were compared between apnea and hyperpnea using paired t-tests. The statistical software used was Sigmastat 2.03 (SPSS). Results were considered significant at P < 0.05.
| RESULTS |
|---|
|
|
|---|
Frequency spectral analysis. Figure 2 shows tracings from a representative subject demonstrating no significant respiratory arrhythmia during regular breathing. The power spectrum shows only a "white noise" pattern characteristic of AF (14, 44). Figure 3 illustrates periodic respiratory-related oscillations in VR R-R intervals at the VLF of CSR-CSA in the same subject. The R-R intervals increase (i.e., VR falls) during apnea, and the R-R intervals decrease (i.e., VR rises) during hyperpnea. Evidence of respiratory arrhythmia at HF was absent in all but one subject, and no subject had evidence of respiratory-related HRV at VLF during regular breathing (Table 2). In contrast, during CSR-CSA, significant respiratory-related HRV at VLF was observed in 12 subjects (92%). The group-average coherence between respiration and R-R interval at VLF during CSR-CSA was significantly greater than that between respiration and R-R interval at either VLF or HF during regular breathing (P < 0.001 for both). Peak power spectral densities of both respiration (input power) and R-R intervals (output power) were approximately sevenfold higher at VLF during CSR-CSA [11.1 l2/Hz (SD 9.5) and 764 ms2/Hz (SD 674), respectively] than at HF during regular breathing [1.71 l2/Hz (SD 2.62) and 103 ms2/Hz (SD 95), respectively, P < 0.05 for both]. The phase angle between VLF oscillations in respiration and HRV was 33° (SD 27), indicating that peak VR followed peak respiration by 7.7 s (SD 3.8).
|
|
|
|
-blockers [68.1 (SD 17) vs. 62.1 beats/min (SD 16), P = 0.006] or not [67.9 (SD 19) vs. 61.0 beats/min (SD 20), P < 0.001], or digoxin [65.2 (SD 16) vs. 60.0 beats/min (SD 16), P < 0.001] or not [57.4 (SD 7) vs. 50.3 beats/min (SD 6.5)]. Both the minimum R-R interval (non-rate dependent) and intercept of the lower envelope at 1.0 s (rate dependent) estimates of AV nodal refractory period were significantly shorter during hyperpnea than apnea [684 (SD 126) vs. 735 ms (SD 147), P < 0.001, and 595 (SD 113) vs. 668 ms (SD 114), P = 0.018, respectively]. Similarly, both estimates of concealed conduction, R-R interval scatter (non-rate dependent), and scattering index (rate dependent) were lower during hyperpnea than during apnea [178 (SD 56) vs. 246 ms (SD 72), P = 0.001, and 431 (SD 120) vs. 494 ms (SD 158), P = 0.029, respectively].
|
|
| DISCUSSION |
|---|
|
|
|---|
Frequency spectral analysis. In the frequency domain, the respiratory variations of CSR-CSA differ from those of regular breathing in two ways: they are of higher magnitude and of lower frequency. It is possible that respiratory variation must exceed a certain level to influence AV nodal conduction. Indeed, the transformation from regular breathing to CSR-CSA resulted in an approximately sevenfold increase in the input power of respiration (ILV), in association with a sevenfold increase in output power of R-R intervals.
The lower frequency of CSR-CSA compared with regular breathing may also be important. Recently, some groups have applied spectral analysis to AF and reported that the ventricular response to AF is not completely random, particularly at lower frequencies. While the power spectrum of AF resembles random "white noise" at HF, it is indistinguishable from the power spectrum associated with sinus rhythm at lower frequencies, with the same 1/frequency fractal component (14). Given the similarity of the power spectrum of HRV in AF to that of sinus rhythm at lower frequencies, VR in AF may also be subject to the same low-frequency influences. Our findings reinforce this concept. VR in AF is chaotic at HF, but it is, nonetheless, strongly influenced by respiration at lower frequencies, responding in much the same way as sinus rhythm (18, 37). In this regard, the AV node in AF may act as a physiological low-pass filter for VR, because it responds only to those influences below a certain frequency threshold.
AV node refractoriness and concealed conduction. AV node refractoriness determines the lower limit of the R-R intervals (which are manifestations of successive impulses successfully conducted by the AV node) during AF and is the most important determinant of mean VR (2, 3). However, the characteristic "irregular irregularity" of VR in AF cannot be explained by AV node refractoriness. Because an atrial impulse is always immediately available at the end of the AV node refractory period, a relatively regular ventricular response might be expected, with each R-R interval being close to the length of the refractory period. The chaotic rhythm characteristic of AF arises not from AV node refractoriness but mainly from concealed AV conduction. This is the phenomenon whereby atrial impulses penetrating incompletely into the AV junction do not reach the ventricles but affect the transmission of subsequent impulses (16, 42). Concealed conduction, by creating a high level of variation in the number and frequency of penetrating impulses, results in the classic chaotic rhythm of AF. The greater the degree of concealed conduction, the greater the scatter of the R-R intervals.
We obtained estimates of AV nodal refractory period and concealed conduction during hyperpnea and apnea using both rate-independent and rate-dependent techniques. Using either technique, the findings were consistent: hyperpnea results in shorter R-R intervals, a shorter AV node refractory period, and a lesser degree of concealed conduction than apnea.
These cyclical variations in AV conduction most likely result from changes in the autonomic neural modulation of the electrophysiological properties of the atria and AV node (13, 14, 24, 45), since it has been demonstrated under a variety of conditions that AV refractory period is shortened by sympathetic activity and lengthened by parasympathetic activity (20, 24, 36, 45). Similarly, parasympathetic stimulation of the atria has been shown to increase the degree of concealed conduction (23). AV conduction has been shown to be altered by hypoxia but only at levels of hypoxia much more extreme than we observed in our subjects (31, 47). Although AV conduction has not been reported to be influenced by direct effects of hypercapnia or cardiac stretch, it is remotely possible that these factors may influence AV conduction.
Oscillations in ventricular response were present, despite almost all of the subjects being on medications that can influence conduction through the AV node and slow VR in AF (26). Mean VR was higher during hyperpnea than during apnea, regardless of whether or not patients were taking
- blockers or digoxin. The lack of effect of these AV nodal blocking agents on our findings may be due to our inability to detect small changes due to the small number of subjects in our study, which precluded assessment of overlapping subgroups. Nor does the failure of our findings to be influenced by medications shed light on the physiological mechanisms, since
-blockers act on the sympathetic nervous system, whereas digoxin acts largely through the parasympathetic nervous system.
Oscillations in autonomic outflow during CSR-CSA might be the direct result of oscillations in central respiratory drive (12, 41). Indeed, sympathetic nervous system activity has been shown to be elevated and to oscillate during CSR-CSA (22, 40). There is evidence for direct connections between respiratory and cardiovascular sympathetic neurons in the brain stem. Activation of the respiratory neurons can coactivate these sympathetic neurons (12) in animal preparations. In intact humans, there is evidence both for and against respiration directly causing sympathoexcitation. Isocapnic hyperventilation has been shown by one group to diminish baroreflex suppression of muscle sympathetic nerve activity (MSNA) (39), but another did not detect any increase in MSNA with increasing levels of respiratory motor output (33). However, we have previously reported, in a heart failure patient with AF and CSR-CSA, that MSNA rose during hyperpnea and fell during apnea in concert with rises and falls in VR and BP, respectively, (5), and that these changes in VR and BP are not eliminated by correction of hypoxia (17).
While changes in both sympathetic and parasympathetic tone can influence AV conduction, our findings are more compatible with sympathetic influences. Our subjects had heart failure, which is associated with greatly diminished parasympathetic modulation of heart rate (11). Indeed, the loss of the normal parasympathetic modulation of heart rate was the earliest autonomic abnormality described in both experimental and human heart failure (1, 4, 9, 30). Considering that respiratory modulation of both sympathetic and parasympathetic outflow has been shown to be dependent on the level of preexisting tone (10), the observed changes in AV properties that correspond with oscillations in BP are more likely to be the result of sympathetic than of parasympathetic modulation. However, the larger VT values accompanying CSR-CSA might conceivably engage parasympathetic reflexes (35) not active during regular breathing in these patients.
Because hyperpnea-induced increases in BP following apnea have been shown to be due to increases in sympathetic outflow (15), our observation of synchronous BP oscillations in three patients lends further support for a sympathetic mechanism being responsible for our findings.
Another way in which we might distinguish between the sympathetic and parasympathetic nervous systems is to examine the latency of the responses. The oscillations in mean VR and BP lagged changes in respiration by
8 s. Thus there appears to be a relatively long time constant between generation of central respiratory drive, ventilation, and the subsequent VR and vasomotor responses. This 8-s delay is again more in keeping with the behavior of the sympathetic than the parasympathetic nervous system, which has a shorter time constant (27). On the other hand, the observed differences in concealed conduction have been linked more closely to the effects of parasympathetic (25, 42) than sympathetic (43) stimulation, implying that both arms of the autonomic nervous system might be influenced by CSR-CSA.
Alterations in blood gases might be a necessary intermediate step by which CSR-CSA exerts its cyclical effects on the autonomic nervous system. The sympathoexcitatory effects of CO2 are likely to be particularly relevant in patients with CHF and CSR-CSA, in whom there is evidence of an enhanced muscle sympathetic nerve response to CO2 (21). Dips in oxygen saturation also accompany CSR-CSA, and, because of the circulatory delay, such dips are maximal in the middle of the hyperpnea and should cause stimulation of sympathetic activity at that time. Thus, although the hypoxia we observed was relatively mild [nadir SaO2 90% (SD 4)], we cannot exclude its role in causing the observed oscillations in VR, particularly since the peripheral chemoreceptors are sensitized in CHF (32, 34). Nevertheless, our laboratory has previously shown that oscillations in heart rate and BP occur during periodic breathing even in the absence of hypoxia (17).
Limitations. Our study is limited by the lack of direct measurements of sympathetic and parasympathetic activity. Instead, the activity of the autonomic nervous system has been surmised by analyzing its putative effect on AV nodal electrophysiology. Moreover, we cannot exclude the role of intermediate steps between the generation of central respiratory drive and the observed electrophysiological changes, such as alterations in blood gases, lung stretch reflexes, or direct effects on cardiac stretch.
It is remotely possible that CSR-CSA influences the electrophysiology of the AV node by a completely nonautonomic mechanism, such as a direct effect of local hypoxia or cardiac stretch. However, no such direct effects of these factors on AV nodal conduction have been reported, with the exception of much more extreme levels of hypoxia (31, 47).
In conclusion, in subjects with heart failure and AF, VLF oscillations in VR occur in time with periodic oscillations in ventilation during CSR-CSA, such that mean VR is higher during hyperpnea than during apnea. In addition, we have demonstrated that the electrophysiological mechanism mediating these oscillations in VR is cyclical alterations in AV nodal refractory period and concealed conduction that occur in time with oscillations in ventilation. The precise means by which CSR-CSA influences these properties of the AV node remain to be determined, but changes in sympathetic or parasympathetic activity are the most likely candidates. To determine whether these effects are mediated centrally, or peripherally through alterations in blood gases, or engagement of lung stretch reflexes will require further investigation.
| GRANTS |
|---|
|
|
|---|
| 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 |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. Bollmann, D. Husser, L. Mainardi, F. Lombardi, P. Langley, A. Murray, J. J. Rieta, J. Millet, S. B. Olsson, M. Stridh, et al. Analysis of surface electrocardiograms in atrial fibrillation: techniques, research, and clinical applications. Europace, November 1, 2006; 8(11): 911 - 926. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |