Abstract
BACKGROUND: Sleep-disordered breathing (SDB) is critically associated with cardiovascular mortality and morbidity, especially in patients with heart failure (HF). However, the majority of SDB patients remain undiagnosed. In contrast, abnormality in heart rate variability has been reported in patients with SDB. To explore an efficient electrocardiogram (ECG)-based screening tool for SDB, we examined the usefulness of cyclic variation in heart rate score (CVHRS) by Holter ECG in subjects with HF.
METHODS: In this study, 102 subjects with HF were enrolled. We simultaneously performed Holter ECG with overnight portable sleep monitoring, and we measured the respiratory disturbance index (RDI) and CVHRS. We determined the temporal position of the individual dips comprising the CVHRS using time-domain methods. CVHRS was measured as cyclic and autocorrelated dips in smoothed interbeat interval time series.
RESULTS: There were 25 subjects with severe SDB (RDI ≥ 30 events/h) and 77 subjects with none-to-moderate SDB (0 ≤ RDI < 30 events/h). There was a significant positive correlation between CVHRS and RDI (r = 0.60, P < .001). In receiver operating characteristic analysis, CVHRS (cutoff of 30 events/h) identified severe SDB with a sensitivity of 82%, a specificity of 77%, and an area under the curve of 0.83.
CONCLUSIONS: CVHRS determined by Holter ECG is a useful screening index for severe SDB in subjects with HF.
Introduction
Despite recent advancements in its medical management, heart failure (HF) still leads to high morbidity and mortality. The lifetime risk of developing HF is 20% at ≥ 40 y of age, and the absolute mortality rate for HF remains ∼50% within 5 y of diagnosis.1 To improve survival, identification of factors that contribute to increased mortality is required. About 50% of patients with HF have sleep-disordered breathing (SDB), which consists of obstructive sleep apnea (OSA) and Cheyne-Stokes respiration with central sleep apnea. Some large-scale studies have demonstrated that SDB is associated with occurrence of ventricular arrhythmias2,3 and adverse prognosis in subjects with HF.4,5 Studies have demonstrated that treatment of SDB may improve prognosis in subjects with HF.6–9 However, because of the lack of patients' and physicians' awareness of SDB, especially in patients with HF, and limited access to portable sleep monitoring or overnight polysomnography, the majority of SDB patients remain undiagnosed.10 Javaheri et al10 reported that only 2% of subjects with HF received SDB testing. Thus, we are looking for alternative tools that are less labor-intensive and cheaper with better accessibility for patients with HF.
Abnormality in heart rate variability11 has been reported in subjects with HF and/or SDB.3,12,13 SDB is associated with an altered sympathovagal balance determined by nocturnal cyclic alternating of apneas and hyperventilation-bradycardia during apnea, followed by abrupt tachycardia.14 This phenomenon causes cyclic variation in heart rate (CVHR).14 Not only OSA but also Cheyne-Stokes respiration with central sleep apnea demonstrated heart rate oscillations.15 Both types of SDB present cyclic lengthening/shortening in R-R interval during apnea/post-apneic hyperventilation.16 Until now, this ordered variation in heart rate has been applied to the detection of SDB by only a few study groups for the general population who received polysomnography.17–20 Previous studies showed data for only OSA in the general population, excluding subjects with HF. Patients with HF show reduction in heart rate variability power, take β blockers, and have complex SDB, including not only OSA but also Cheyne-Stokes respiration with central sleep apnea. Prevalence of Cheyne-Stokes respiration with central sleep apnea is increased with HF.21 Hence, techniques that could provide a simple screening tool for SDB are needed. The Holter electrocardiogram (ECG) is performed for many patients with HF. Screening for SDB by Holter ECG may be useful for many patients with HF without access to polysomnography. To explore an efficient Holter ECG-based screening tool for SDB in patients with HF, we examined the usefulness of CVHR score (CVHRS) by Holter ECG.
QUICK LOOK
Current knowledge
Sleep apnea and sleep-disordered breathing (SDB) are associated with significant morbidity and mortality in patients with heart failure (HF). SDB in HF is also largely undiagnosed.
What this paper contributes to our knowledge
Cyclic variation in heart rate (CVHR) and respiratory disturbance index were significantly correlated in subjects with severe SDB. Receiver operating characteristic analysis demonstrated a high sensitivity and specificity for heart rate variation and SDB severity. Monitoring the CVHR score may be a useful screening index for SDB in HF.
Methods
Study Subjects and Study Protocol
In this study, we enrolled 110 consecutive subjects with HF who were referred for an overnight portable sleep monitor test with hospitalization regardless of SDB symptoms (such as daytime sleepiness determined by the Epworth sleepiness scale22) at Fukushima Medical University. Inclusion criteria were: (1) presence of symptomatic HF, defined as New York Heart Association class ≥ II1; (2) enforcement of standard pharmacotherapy (including angiotensin-converting enzyme inhibitors, angiotensin II receptor blockers, β blockers, and diuretics); and (3) stable clinical status, which was defined as receiving optimal medical therapy without worsening of HF for at least 2 months before study enrollment. Exclusion criteria were: (1) presence of atrial fibrillation or pacemaker implantation, (2) overnight recordings with > 10% noise or ectopic beats, (3) history of stroke with neurological deficit, (4) acute coronary syndrome and acute decompensated HF, and (5) recent SDB treatment. The final total was 102 subjects.
We performed simultaneous overnight portable sleep monitoring and Holter ECG monitoring. Standard Holter ECG recorders (LS-300, Fukuda Denshi, Tokyo, Japan) were used to acquire data. Polygraphic data and Holter ECG data were matched using polygraphic time synchronized with Holter ECG time before lights out. Two independent physicians, each one blinded to the results of the other, performed heart rate variability analysis and synchronized polygraphic scoring. Subjects were divided into 2 groups based on whether the respiratory disturbance index (RDI) was > 30 events/h (severe SDB, n = 25) or not (none-to-moderate SDB, n = 77) via portable sleep monitor. Written informed consent was obtained from all study subjects. The study protocol was approved by the ethics committee of Fukushima Medical University (number 1050).
Portable Sleep Monitor
All subjects underwent overnight polygraphy with standard techniques.23 An overnight sleep study was performed using the LS-300 cardiopulmonary monitoring system (type 3 polygraph). ECG, thoracoabdominal motion, nasal air flow via an air-flow pressure transducer, and SpO2 were monitored as reported previously.23 Apnea was defined as an absence of air flow for > 10 s. Hypopnea was defined as a > 30% reduction in monitored air flow accompanied by a decrease in SpO2 of > 3%.24 Standard definitions were used for OSA and central sleep apnea based on the presence or absence of rib cage and abdominal excursions with an absence of air flow.24 RDI was defined as the number of apneas and hypopneas/h during time spent in bed.24 The major polygraphic parameters investigated were RDI, central RDI, obstructive RDI, 3% oxygen desaturation index (3% ODI), lowest SpO2, and mean SpO2.23 These data were visually inspected and scored by a single experienced laboratory technician who was blinded to the other results. Our preliminary data from another cohort (n = 16) showed that the RDI obtained with the LS-300 system closely correlated with the apnea-hypopnea index obtained by polysomnography (r = 0.93, P < .01; see the supplementary materials at http://www.rcjournal.com).
ECG Analysis and CVHRS Measurement
ECG signals were obtained via the Holter LS-300 system at a sampling frequency of 125 Hz.25 A Holter automatic arrhythmia analysis system (SCM-8000, Fukuda Denshi) identified all R-wave positions and excluded abnormal beats, such as ventricular ectopic and supraventricular ectopic beats and artifacts. Moreover, the R-R interval of normal beats that exceeded ± 30% of the average R-R interval was removed from the trend signal. Because the R-R interval time series data had irregular intervals after excluding all abnormal beats, a resampling signal process was performed at 2 Hz empirically. For automated detection of cyclic variation in heart rate, we developed a new algorithm. The methods used to create the parameters to detect dip candidates are shown in Figure 1. We created a differentiation waveform, as shown in Figure 1C, from R-R time series data f(i), which was obtained after the resampling signal process, in Figure 1B through a band-pass filter (0.02–0.11 Hz). Three points of intersection (x1, x2, and x3) in Figure 1C and baselines were created as dip starting point y1, lowest dip point y2, and dip ending point y3. We defined y1 and y3 as the maximum points of x1 and x3 in the vicinity of 3 s as y2 as the minimum point of x2 in the vicinity of 3 s (Fig. 1B).
Next, the height of dip h was calculated in each of the dip candidates as: h = 0.8[(y1 − y2) + (y3 − y2)]/2 (Fig. 2B). In another cohort with no cardiovascular disease (n = 36), we postulated that dip h is calculated as: h = [(y1 − y2) + (y3 − y2)]/2. We then defined the width of dip d as the interval points of f(i) = h, as shown in Figure 1B, whereas the true dip d was defined as the inter-apnea interval (time during post-apneic hyperventilation). The coefficient was determined as 0.8 using the least-squares method in the other cohort. This value was adopted as the coefficient since the mean squared error was smallest when the parameter was 0.8 (see the supplementary materials at http://www.rcjournal.com). All h dips were measured from f(i). How significant dips were picked from dip candidates is shown in Figure 2. We determined the window width of each dip candidate as 6 min (± 3 min from the dip candidate) in which at least one cycle of long apnea duration could be included because the longest apnea duration seemed to be 3 min empirically. The mean + 1.5 SD of dip h in a window of 6 min was defined as the threshold for the dips for excluding the dip candidates that might be respiratory changes (noises). To determine the SD coefficient as 1.5, we calculated the CVHRS (events/h), changing the value of the SD coefficient from the database for developing this algorithm. We adopted 1.5 SD because the root mean square of sensitivity and specificity was highest with 1.5 as the SD coefficient using the Youden index method in another cohort (n = 36; see the supplementary materials at http://www.rcjournal.com). We calculated the mean + 1.5 SD of 3 min before and after the vicinity of the dip candidate, and it served as the threshold for the dip candidates; when the height h in a candidate point exceeded the average of + 1.5 SD, this candidate point was recognized as the change from apnea to hyperventilation in the R-R time series data. Thus, we calculated the CVHRS as the mean number of CVHR/h of time in bed, which was defined by the behavior records of subjects. In the other cohort with no cardiovascular disease (n = 36), the CVHRS and RDI obtained with the LS-300 system closely correlated (r = 0.86, P < .01; see the supplementary materials at http://www.rcjournal.com).
Measurement of Laboratory and Echocardiographic Data
Blood samples were obtained the morning after polysomnography while subjects were in a supine position. Plasma brain natriuretic peptide levels were measured using a specific immunoradiometric assay (ShionoRIA BNP kit, Shionogi, Osaka, Japan). The estimated glomerular filtration rate was measured using the Modification of Diet in Renal Disease formula.26 Echocardiography was performed using standard techniques by an experienced echocardiographer at the echo laboratory in our hospital during the day. Echocardiographic parameters investigated included left ventricular ejection fraction and right ventricular fractional area change.27 The left ventricular ejection fraction was calculated using a modification of Simpson's method. The right ventricular fractional area change, defined as (end diastolic area − end systolic area)/end diastolic area × 100, is a measure of right ventricular systolic function.27 All recordings were performed using ultrasound systems (Acuson Sequoia, Siemens, Mountain View, California).
Statistical Analysis
The Kolmogorov-Smirnov test was used for assessment of normal distribution. Normally distributed data are presented as mean ± SD, and non-normally distributed data are presented as the median (interquartile range). Categorical variables are expressed as frequencies and relative frequencies. Characteristics between the 2 groups were compared using the independent Student t test for normally distributed data and the Mann-Whitney U test for non-normally distributed data for continuous variables, whereas the chi-square test was used for categorical variables. Correlations between the RDI and CVHRS were assessed using Pearson correlation analysis. P < .05 was considered significant. To assess the ability of the CVHRS to correctly rank order subjects according to detected severe SDB, we measured the area under the receiver operating characteristic curve. The outcome we detected was severe SDB. These analyses were performed using SPSS 21.0 (IBM, Armonk, New York).
Results
Case Presentation
Figure 3 shows a subject's polygraph and Holter ECG record, including R-R interval, nasal air flow, thoracoabdominal motion, and SpO2. SDB subjects showed cyclic fluctuations in heart rate, consisting of bradycardia during apnea followed by tachycardia upon its cessation, which is called CVHR. CVHR makes dips in R-R interval series. Figure 3A shows recordings from a subject with HF and mild SDB (RDI of 7.2 events/h and CVHRS of 6.3 events/h). Figure 3B shows those from a subject with HF and severe SDB (RDI of 86.9 events/h and CVHRS of 70.9 events/h).
Clinical Characteristics of Study Subjects
Twenty-five subjects had severe SDB (RDI ≥ 30 events/h), and 77 subjects had none-to-moderate SDB (0 ≤ RDI < 30 events/h). Clinical characteristics of subjects with severe SDB and none-to-moderate SDB are compared in Table 1. Standard pharmacotherapy, including angiotensin-converting enzyme inhibitors, angiotensin II receptor blockers, and β blockers, were given to the majority of subjects. There were no significant differences in age, gender, medications, laboratory data, and echocardiographic data between the 2 groups. The RDI, central RDI, and obstructive RDI were significantly higher, whereas the lowest SpO2 and mean SpO2 were significantly lower in severe SDB group compared with the none-to-moderate SDB group.
Correlation Between Polygraphic Data and CVHRS
As shown in Figure 4, there was a significant positive correlation between CVHRS and RDI (r = 0.60, P < .001). As shown in Figure 5, receiver operating characteristic analysis revealed that the CVHRS (cutoff of 30 events/h) identified severe SDB with a sensitivity of 82%, a specificity of 77%, and an area under the curve of 0.83 (95% CI 0.72–0.93).
Discussion
In this study, CVHRS obtained with our automated CVHR detection algorithm showed good performance in identifying severe SDB in subjects with HF. To our knowledge, this is the first study to show that CVHRS is an excellent screening tool to detect severe SDB even in subjects with HF who have a reduction in power of heart rate variability, have taken β blockers, and have complex SDB, including OSA and Cheyne-Stokes respiration with central sleep apnea.
Impact of SDB on Hemodynamic Parameters and Heart Rate Variability in HF
SDB (especially OSA) is characterized by recurrent hypoxia, arousal, and generation of exaggerated negative intrathoracic pressure, which increases sympathetic nervous activity, reduces cardiac parasympathetic activity, and causes repetitive surges in heart rate, blood pressure, and left ventricular preload and afterload, resulting in decreasing stroke volume.28 Yumino et al28 reported the different hemodynamic impact of OSA and Cheyne-Stokes respiration with central sleep apnea on subjects with HF: whereas OSA appears to have a decreasing stroke volume, Cheyne-Stokes respiration with central sleep apnea seems to have little or slightly increased stroke volume. With respect to heart rate, not only OSA but also Cheyne-Stokes respiration with central sleep apnea demonstrated heart rate oscillations.15 Both types of SDB presented cyclic lengthening/shortening in the R-R interval during apnea/post-apneic hyperventilation.16 OSA and Cheyne-Stokes respiration with central sleep apnea are associated with a different autonomic control of heart rate. The respiratory component of the R-R interval variability (respiratory sinus arrhythmia) is reportedly highest during the apneic phase of OSA and is blunted during hyperventilation. Conversely, in central sleep apnea, the respiratory component of R-R interval variability is nearly absent during the apneic phase and recovers with resumption of ventilation.16 Because of these features, episodes of Cheyne-Stokes respiration with central sleep apnea may be detected more consistently with the CVHRS compared with episodes of OSA; this may have led to the tendency for overestimation in subjects with Cheyne-Stokes respiration with central sleep apnea.15,29–31
ECG-Based SDB Screening in Subjects With HF
Cyclic variations in heart rate index accompanying SDB driven by respiratory control dysfunction are known to show more variable and irregular patterns of oscillations.32 Recently, several authors have developed algorithms for automated ECG detection of OSA in the general population using heart rate variability, including time and frequency domain analysis of R-R intervals, and some studies have reported that the combination of parameters of heart rate variability and ECG-derived respiration signals gave the best classification result.17–20 Roche et al33 proposed an increase in the relative power of the very-low-frequency component (0.01–0.05 Hz) of interbeat interval increment as a marker for OSA. They reported an area under the curve of 0.70 for identifying subjects with moderate SDB with 64% sensitivity and 69% specificity. Hayano et al34 used a time-domain and dip-detection method and algorithm (which they called ACAT algorithms) in the general population. The authors reported an area under the curve of 0.91 for identifying subjects with moderate SDB with 83% sensitivity and 88% specificity.
In patients with HF, the power of heart rate variability decreases and is associated with adverse prognosis.35,36 There is convincing evidence that autonomic nervous imbalance caused by SDB contributes to increased mortality in patients with HF.37 Heart rate variability powers in subjects with HF and SDB were reportedly decreased during a 24-h period by increased sympathetic activity and damped circadian oscillation.3,38 Furthermore, patients with HF take β blockers, which affect heart rate variability. Our study differs from previous studies in several ways. (1) We presented a new algorithm named CVHRS, and (2) our study population consisted of subjects with HF, who have a reduction in heart rate variability power, take β blockers, and have complex SDB, including OSA and Cheyne-Stokes respiration with central sleep apnea.
Limitations
Our study has some limitations. First, subjects with atrial fibrillation were excluded because we could not determine the CVHRS in the presence of atrial fibrillation. Second, the influence of medical treatment on these study findings remains uncertain. Moreover, a previous study has shown that β blockers change autonomic nervous balance in subjects with HF.39 Third, the accepted standard for diagnostic testing of SDB is full-channel polysomnography, which provides detailed information on complete differentiation of the types of apnea. Our differentiation of SDB using a portable sleep monitor might be less reliable than that by full-channel polysomnography. Thus, further investigation by full-channel polysomnography may lead to complete clarification of CVHRS in subjects with HF and SDB. Fourth, we cannot exclude the effect of periodic leg movements, which may affect heart rate variability.34,40 Periodic leg movements are characterized by recurrent episodes of repetitive stereotyped limb movements that occur predominantly in the legs.40 Finally, the sample size was relatively small. Large-scale clinical trials may be required to confirm the CVHRS algorithm to determine SDB in subjects with HF.
Conclusions
SDB screening by Holter ECG may be useful for many patients with HF who do not have access to polysomnography. CVHRS determined by Holter ECG is a useful screening index for severe SDB in patients with HF. We do not consider CVHRS to be a final diagnostic tool, but rather a screening tool for SDB. We hope that patients with HF and suspected SDB screened by Holter ECG will then be able to undergo the more precise polysomnography.
Acknowledgments
We thank Ms Yuko Niimura for outstanding technical assistance.
Footnotes
- Correspondence: Akiomi Yoshihisa MD PhD, Department of Cardiology and Hematology, Fukushima Medical University, 1 Hikarigaoka, Fukushima 960-1295, Japan. E-mail: yoshihis{at}fmu.ac.jp.
Supplementary material related to this paper is available at http://www.rcjournal.com.
This study was supported in part by Grant-in-Aid for Scientific Research 25461061 from the Japan Society for the Promotion of Science and grants-in-aid from the Japanese Ministry of Health, Labour, and Welfare (Tokyo, Japan). The Department of Advanced Cardiac Therapeutics is supported by Fukuda Denshi. The authors have disclosed no other conflicts of interest.
- Copyright © 2015 by Daedalus Enterprises