(Circulation. 1999;100:1411-1415.)
© 1999 American Heart Association, Inc.
Clinical Investigation and Reports |
From Service d'Exploration Fonctionnelle CardioRespiratoire, Laboratoire de Physiologie (F.R., V.P., D.D., F.C., J.-C.B.) and Service de Pneumologie et d'Oncologie Thoracique (I.C.-F.), CHU Nord, Faculté de Médecine Jacques Lisfranc, Université Jean Monnet, Saint-Etienne, France; Clinique de Médecine II et Division de Cardiologie, Département de Médecine Interne, Hôpitaux Universitaires, Geneva, Switzerland (J.-M.G.); Ecole Nationale de Statistiques et d'Analyse de l'Information, Rennes, France (P.M.); and Laboratoire de Physiologie, Faculté de Médecine Lyon Sud, Université Lyon I, Oullins, France (J.-R.L.).
| Abstract |
|---|
|
|
|---|
Methods and ResultsIn G1, 39 of 91 patients (42.8%) were
identified as diseased by polysomnography, as were 24 of 52 patients
(46%) in G2. Time-domain HRV variables (SD of NN intervals
[SDNN], mean of the standard deviations of all NN intervals for all
consecutive 5-minute segments of the recording [SDNN index], square
root of the mean of the sum of the squares of differences between
adjacent normal RR intervals [r-MSSD], and SD of the averages of NN
intervals in all 5-minute segments of the recording [SDANN])
were calculated for daytime and nighttime periods, as well as the
differences between daytime and nighttime values (
[D/N]).
Correlations between HRV variables and OSAS status were
analyzed in G1 by use of receiver-operating characteristic
(ROC) curves and logistic regression analysis. By ROC curve
analysis, 7 variables were significantly associated with
OSAS. After adjustment for other variables through multiple
logistic regression analysis,
[D/N]SDNN index and
[D/N] r-MSSD remained significant independent predictors of OSAS,
with ORs of 8.22 (95% CI, 3.16 to 21.4) and 2.86 (95% CI, 1.21 to
6.75), respectively. The classification and regression tree methodology
demonstrated a sensitivity reaching 89.7% (95% CI, 73.7 to 97.7) with
[D/N] SDNN index and a specificity of 98.1% (95% CI, 86.4 to
100) with
[D/N] SDNN using appropriate thresholds. These
thresholds, applied to G2, yielded a sensitivity of 83% using
[D/N] SDNN index and a specificity of 96.5% using
[D/N]
SDNN.
ConclusionsTime-domain HRV analysis may represent an accurate and inexpensive screening tool in clinically suspected OSAS patients and may help focus resources on those at the highest risk.
Key Words: sleep heart rate nervous system, autonomic
| Introduction |
|---|
|
|
|---|
Abnormalities in nocturnal cyclical heart rate (HR) variations have previously been described in sleep-related breathing disorders.3 We hypothesized that a quantification using time-domain analysis of heart rate variability (HRV), which already demonstrated its usefulness in other clinical fields,4 5 could yield more pertinent diagnostic information because of better assessment of the autonomous nervous system activity. Thus, we chose to evaluate the sensitivity and specificity of this seemingly simple tool in patients referred to our clinic for OSAS diagnosis. We correlated HRV variables and diseased status assessed with complete polysomnography to establish pertinent HRV diagnosis thresholds in 1 group of patients; then, the selected HRV thresholds were applied to comparable patients to validate our findings.
| Methods |
|---|
|
|
|---|
Sleep Study and Polysomnography Scoring
Patients of both groups underwent nocturnal polysomnography
recording with standard methods.6 7 The presence
and stages of sleep were monitored by use of 2 pairs of EEG leads
(C4-A1 and
O2-A1), 2 pairs of
electro-oculographic leads, and chin electromyographic leads. Air flow
was measured by an oronasal thermocoupler. Respiratory efforts were
monitored by use of respiratory plethysmography with transducers placed
around the chest and abdomen. Arterial oxygen saturation
was continuously recorded by pulse oximetry (Criticare Systems Inc)
during the entire night.
The polysomnogram was scored manually according to standard
criteria.8 Apnea was defined as the absence of air flow
for >10 seconds in the presence of persistent respiratory efforts.
Hypopnea was defined as the association of a reduction of
50% of the
amplitude of respiratory efforts during
10 seconds and a decrease in
arterial oxyhemoglobin saturation of
4%. The
apnea/hypopnea index (AHI) was defined as the number of episodes of
apnea and hypopnea per hour of sleep; the threshold to accept the
diagnosis of OSAS was chosen as an AHI
10; patients were thus
classified as diseased (OSAS+) or nondiseased (OSAS-). Each episode of
apnea was characterized by measuring apnea duration and mean and
minimal arterial oxyhemoglobin saturation. Total sleep
time, number and duration of REM periods, and number and duration of
arousals were measured. Total duration of arterial
desaturation was quantified as cumulated time with arterial
oxyhemoglobin saturation <90%.
ECG Holter Monitoring HRV Analysis
The Holter tapes were analyzed on a StrataScan 563 (Del
Mar) equipped with an HRV module. All variables were calculated for
the 24-hour, daytime (2:30 to 9:30 PM), and nighttime
(midnight to 7 AM) periods, and the differences between
daytime and nighttime values (
[D/N]) were computed.
HRV was evaluated by use of time-domain analysis. To perform the analysis, each QRS complex was identified, and the length between each QRS complex (RR interval) was calculated. Only normal-to-normal beats were considered for analysis. Several parameters describing the differences between RR intervals were calculated9 : the square root of the mean of the sum of the squares of differences between adjacent normal RR intervals (r-MSSD), SD of NN intervals (SDNN), SD of the averages of NN intervals in all 5-minute segments of the recording (SDANN), and mean of the SD of all NN intervals for all consecutive 5-minute segments of the recording (SDNN index).
Statistical Analysis
Data were analyzed with StatView and JMP (SAS Institute)
software. Differences were considered significant when
P<0.05. Values are expressed as mean and 95% SD
(mean±SD).
In group 1 (the derivation set), statistical analyses were
performed to evaluate the ability of each variable to discriminate
between diseased and nondiseased status. Thus, the dependent
variable was diseased status (OSAS+/-). The independent
variables analyzed were age, sex, body mass index (BMI),
and for HRV variables, their day and night values and the
differences between their night and day values (
[D/N]), as night
mean HR,
[D/N] mean HR, night r-MSSD,
[D/N] r-MSSD, night
SDNN,
[D/N] SDNN, night SDNN index,
[D/N] SDNN index, night
SDANN, and
[D/N] SDANN.
Receiver-operating (ROC) curve analysis was used10 with calculation of the area under the curves (AUC). An AUC value of 0.5 means that the variable distributions are similar in both populations; conversely, an AUC value of 1.0 means that the 2 populations' variable distributions do not overlap at all.
Logistic regression models were used also to analyze the data via a simple logistic regression of the diseased status versus each of the covariates to confirm the previously performed ROC curve analysis. Multiple logistic regression analysis was then used to evaluate the strength of the association of each variable with diseased status after adjustment for the other variables. Also, the ORs were calculated individually for the same most significant independent variables.
Finally, a threshold value was determined for the most significant variables to obtain better separation power between diseased and nondiseased status. Determination of the threshold values was based on the classification and regression tree (CART) methodology,11 which allows continuous variables to be split into 2 groups. All possible values were considered, and the chosen value was the one that maximized the Gini impurity criterion.12 Sensitivity and specificity were also estimated against the diseased and nondiseased status with their simultaneous "exact" 95% CI, ie, the 95% CI for sensitivity and specificity.
In group 2 (the validation set), sensitivity and specificity were calculated for diseased and nondiseased status by use of the 3 most significant independent variables identified by CART methodology in group 1.
| Results |
|---|
|
|
|---|
[D/N] mean HR,
[D/N]
SDNN,
[D/N] SDNN index, and
[D/N] r-MSSD.
|
ROC curves (continuous data) were built for each HRV variable; 7
were able to separate OSAS+ from OSAS- status with statistical
significance (P=0.03 to <0.0001; Table 2
). Once the 7 variables were
classified according to their AUC values, the absolute
[D/N] SDNN
index appeared to be the most powerful separator
(P<0.0001), followed by the night SDNN index,
[D/N]
SDNN, night SDNN,
[D/N] r-MSSD, and night r-MSSD; for each of
these variables, the absolute differences between the day and night
values were more potent predictors than the night values alone.
Finally, the difference between night and day mean HR was the
last significant variable. None of the time-domain values
calculated over the day period reached statistical significance.
|
For each continuous variable, the thresholds giving a chosen
sensitivity or specificity were calculated (Table 3
). For the most powerful variable,
[D/N] SDNN index, the 90% sensitivity threshold was -10 ms and
the 90% specificity threshold reached -43 ms.
|
The highest separation power was obtained by use of the CART
methodology for the most powerful variables to derive a threshold
value with an optimal sensitivity or specificity (Table 4
). We thus recognized
[D/N] SDNN
index as an efficient tool for the detection of OSAS and
[D/N]
SDNN as an efficient parameter to exclude OSAS. For the
[D/N] SDNN index variable, sensitivity reached 89.7% (95%
CI, 73.7 to 97.7) with a threshold value of -11.1; for the
[D/N]
SDNN variable, specificity reached 98.1% (95% CI, 86.4 to 100)
with a threshold value of -34.1.
|
A separate simple logistic regression analysis of diseased
versus nondiseased status was performed for each of the covariates to
verify the ROC curve analysis and to select a reasonable subset
of explanatory variables for further investigation. Maximum
likelihood estimates were obtained for each of the important
explanatory covariates (Table 5
). These
results confirmed the ROC analysis except for the
[D/N]
mean HR.
|
With the use of a multiple logistic regression analysis, only 2
time-domain HRV variables appeared to be independently and
significantly associated with diseased status. After adjustment for the
other variables,
[D/N] SDNN index (P<0.001) and
[D/N] r-MSSD (P<0.01) remained significant predictors
of OSAS (adjusted OR, 8.22 and 2.86; 95% CI, 3.16 to 21.4 and 1.21 to
6.75, respectively).
In group 2, neither BMI (28.7±9.8) nor age (52.7±12) was
significantly different from group 1. Of the 52 patients, 24 (46%)
were diagnosed as OSAS+ with complete polysomnography. The
repartitioning of patients according to associated disease status
(diabetes, hypertension, chronic heart failure, and ischemic
heart disease) was similar to that of group 1; the threshold value of
-11.1 for
[D/N] SDNN index determined a sensitivity of 83%, and
the threshold value of -34.1 for
[D/N] SDNN determined a
specificity of 96.5%.
| Discussion |
|---|
|
|
|---|
Guilleminault et al3 described a strong relationship between the longest apneic events in OSAS and the greatest swings, illustrating what they called cyclical variation HR. These authors already proposed using this feature as a screening tool in OSAS, but this has not been applied. Similarly, Le Heuzey et al13 described greater visual enlargement of the RR trend and a significantly higher difference between minimal and maximal HR during the night in OSAS patients than in control subjects. However, these data were not quantified according to the criteria of the classic time-domain analysis accepted today.9 Meanwhile, both studies suggested using this test as a tool for preliminary identification of OSAS independently of the assistance of a sleep laboratory. At that time, quantification of RR-interval variations were not automated. Quantification of RR-interval variability has recently become a common feature in Holter systems, allowing easy calculation of HRV scores.
Spectral analysis of HRV had already been proposed by Shiomi et
al,14 who used it to evaluate the benefits of the
prosthetic mandibular advancement treatment in OSAS patients;
this treatment decreased the power of the very low frequencies without
other spectral changes. Furthermore, they did not find any significant
change in the time-domain variables. Conversely, we found pertinent
criteria in time-domain variables that allowed OSAS diagnosis by
creating specific new variables with enhanced discriminant
capabilities. These new variables were calculated as the intervals
between the lowest and highest boundary values of 3 classic
variables, benefiting from an increased contrast between day and
night HRV values, particularly
[D/N] SDNN index and
[D/N]
r-MSSD.
The physiological basis of the observed abnormalities lies in the struggle against upper-airway obstruction, first stimulation of the parasympathetic arm of the autonomous nervous system, followed by an abrupt sympathetic activation caused by the resultant hypoxia. Thus, at night, the alternate strong successive parasympathetic and sympathetic drives dramatically enhance RR variability; of course, this enhanced HRV does not correspond to an isolated increase in parasympathetic activity as in healthy subjects. By contrast, the diurnal HRV is known to be reduced in OSAS patients,15 and it has already been demonstrated to be associated with enhanced sympathetic nerve activity.15 16 We hypothesize that the contrast observed between day and night values resulted mainly from the combination of these 2 phenomena.
The strength of our study was deriving and validating new variables to obtain a high sensitivity (89.7%) and high specificity (98.1%) in the diagnosis of OSAS by choosing appropriate thresholds. Clinical use of these criteria should depend both on the population of interest and on the purpose of the investigator, who may want to point toward a higher sensitivity or a higher specificity. A high sensitivity would allow exclusion of the threat of OSAS by use of a single Holter recording in patients at risk, such as obese cardiac patients, whereas a high specificity would help to increase the probability of the disease if performed as a pretest in polysomnography centers.
This method has several limitations. First, some clinical disorders, such as mellitus diabetes, sequelae of myocardial infarction, and chronic heart failure, present with a blunted autonomic reactivity. These disorders, especially heart failure, are quite often associated with OSAS.17 Although such disorders would favor false-negative diagnosis of OSAS by blunting HRV, we were able to classify correctly 8 of 10 OSAS patients in a population of 23 patients having congestive heart failure (n=7), myocardial infarction (n=4), or diabetes mellitus (n=12). Complete polysomnography also carries some inaccuracies18 : transducer position, failure, and displacement; artifacts resulting from movement; and alteration of sleep because of the technical environment. Finally, whereas arterial oxyhemoglobin saturation monitoring appears to be the only reference tool to quantify the deepness of the arterial oxygen desaturation, this variable currently does not reach 100% sensitivity and specificity in OSAS detection because oxyhemoglobin desaturation can reflect other disorders or can be relatively moderate despite severe structural sleep disorders. In this context, HRV abnormalities could be a better representative of sleep-breathing disorders, identified as autonomic microarousals, than a pure apnea hypopnea index. However, our results of HRV analysis are difficult to compare to commercially available diagnosis systems associating snoring, arterial oxyhemoglobin saturation, and HR monitoring because the results of this last algorithm are not convincing and their recordings are limited to night.2
From a technical point of view, to be consistent in the
comparison between our recordings, we selected the same diurnal
and nocturnal periods of HRV analysis. A change in this
methodological choice could eventually lead to different results. The
threshold selected for the reference diagnosis of OSAS was 10 events
per hour, the more commonly used threshold. Importantly, patients who
have
30 events per hour are identified with the same HRV threshold
value with a sensitivity of 96% and a specificity of 99.7%, values
much higher than for patients with
10 events per hour. Such patients
are certainly candidates for specific long-term therapy, such as
nocturnal nasal continuous positive airway pressure.
In summary, time-domain analysis of HR variability, used as the only criterion, could thus represent an efficient tool in OSAS diagnosis with a sensitivity of 90%. The ease of use and interpretation are also of interest because of the high prevalence of the disease in the general population and the need for repeated control.
| Footnotes |
|---|
Received February 10, 1999; revision received June 4, 1999; accepted June 16, 1999.
| References |
|---|
|
|
|---|
2. Stoohs R, Guilleminault C. MESAM 4: an ambulatory device for the detection of patients at risk for obstructive sleep apnea syndrome (OSAS). Sleep. 1992;101:12211227.
3. Guilleminault C, Connoly S, Winkle R, Melvin K, Tilkian A. Cyclical variation of the heart rate in sleep apnea syndrome. Lancet. 1984;21:126131.
4.
Saul JP. Beat-to-beat variations of heart rate reflect
modulation of cardiac autonomic outflow. News Physiol Sci. 1990;5:3237.
5. Stein PK, Bosner MS, Kleiger RE, Conger BM. Heart rate variability: a measure of cardiac autonomic tone. Am J Cardiol. 1994;127:13761381.
6. Rechtschaffen A, Kales A, eds. A Manual of Standardized Terminology, Techniques and Scoring Systems for Sleep Stages of Human Subjects. Los Angeles, Calif: BIS/BRI UCLA; 1968.
7. Guilleminault C. Sleep and breathing. In: Guilleminault C, ed. Sleeping and Waking Disorders: Indications and Techniques. Stoneham, UK: Butterworth Publishers; 1982:155182.
8. Thorpy MJ, for the Diagnostic Classification Steering Committee. The International Classification of Sleep Disorders: Diagnostic and Coding Manuel. Rochester, Minn: American Sleep Disorders Association; 1990;5258.
9.
Task Force of the European Society of
Cardiology and the North American Society of Pacing and
Electrophysiology. Heart rate variability: standards of measurements,
physiological interpretation, and clinical use.
Circulation. 1996;93:10431065.
10.
Hanley JA, McNeil BJ. The meaning and use of the area
under receiver-operating characteristic (ROC) curve.
Radiology. 1982;143:2936.
11. Breiman L, Friedman JH, Olshen RA, Stone CJ. Classification and Regression Trees. Belmont, Calif: Wadsworth International Group; 1994.
12.
Clopper CJ, Pearson ES. The use of confidence or
fiducial limits illustrated in the case of the binomial.
Biometrika. 1934;26:404412.
13. Le Heuzey JY, Romejko P, Fleury B, Scala PJ, Derenne JP, Guize L, Lavergne T, Valty J. Holter monitoring in the diagnosis of sleep apnea syndrome. J Am Coll Cardiol. 1989;13:189A. Abstract.
14. Shiomi T, Guilleminault C, Sasanabe R, Hirota I, Maekawa M, Koboyashi T. Augmented very low frequency component of heart rate variability during obstructive sleep apnea. Sleep. 1996;19:370377.[Medline] [Order article via Infotrieve]
15.
Narkiewicz K, Montano N, Cogliati C, van de Borne PJH,
Dyken ME, Somers VK. Altered cardiovascular variability
in obstructive sleep apnea. Circulation. 1998;98:10711077.
16. Somers VK, Dyken ME, Clary MP, Abboud FM. Sympathetic neural mechanisms in obstructive sleep apnea. J Clin Invest. 1995;96:18971904.
17.
Javaheri S, Parker TJ, Liming JD, Corbett WS, Nishiyama
H, Wexler L, Roselle GA. Sleep apnea in 81 ambulatory male patients
with stable heart failure: types and their prevalences, consequences,
and presentations. Circulation. 1998;97:21542159.
18. Yamashiro Y, Kryger MH. Nocturnal oximetry: is it a screening tool for sleep disorders? Sleep. 1995;18:167171.[Medline] [Order article via Infotrieve]
This article has been cited by other articles:
![]() |
J. Naghshin, K. R. McGaffin, W. G. Witham, M. A. Mathier, L. C. Romano, S. H. Smith, A. M. Janczewski, J. A. Kirk, S. G. Shroff, and C. P. O'Donnell Chronic intermittent hypoxia increases left ventricular contractility in C57BL/6J mice J Appl Physiol, September 1, 2009; 107(3): 787 - 793. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. de Chazal, C. Heneghan, and W. T McNicholas Multimodal detection of sleep apnoea using electrocardiogram and oximetry signals Phil Trans R Soc A, January 28, 2009; 367(1887): 369 - 389. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Szyszko, C. Franceschini, and J. Gonzalez-Zuelgaray Reliability of a Holter-based methodology for evaluation of sleep apnoea syndrome Europace, January 1, 2009; 11(1): 94 - 99. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Schonauer, A. Thomas, S. Morbach, J. Niebauer, U. Schonauer, and H. Thiele Cardiac autonomic diabetic neuropathy Diabetes and Vascular Disease Research, November 1, 2008; 5(4): 336 - 344. [Abstract] [PDF] |
||||
![]() |
V. K. Somers, D. P. White, R. Amin, W. T. Abraham, F. Costa, A. Culebras, S. Daniels, J. S. Floras, C. E. Hunt, L. J. Olson, et al. Sleep Apnea and Cardiovascular Disease: An American Heart Association/American College of Cardiology Foundation Scientific Statement From the American Heart Association Council for High Blood Pressure Research Professional Education Committee, Council on Clinical Cardiology, Stroke Council, and Council on Cardiovascular Nursing In Collaboration With the National Heart, Lung, and Blood Institute National Center on Sleep Disorders Research (National Institutes of Health) Circulation, September 2, 2008; 118(10): 1080 - 1111. [Full Text] [PDF] |
||||
![]() |
V. K. Somers, D. P. White, R. Amin, W. T. Abraham, F. Costa, A. Culebras, S. Daniels, J. S. Floras, C. E. Hunt, L. J. Olson, et al. Sleep Apnea and Cardiovascular Disease: An American Heart Association/American College of Cardiology Foundation Scientific Statement From the American Heart Association Council for High Blood Pressure Research Professional Education Committee, Council on Clinical Cardiology, Stroke Council, and Council on Cardiovascular Nursing In Collaboration With the National Heart, Lung, and Blood Institute National Center on Sleep Disorders Research (National Institutes of Health) J. Am. Coll. Cardiol., August 19, 2008; 52(8): 686 - 717. [Full Text] [PDF] |
||||
![]() |
W. T. McNicholas Diagnosis of Obstructive Sleep Apnea in Adults Proceedings of the ATS, February 15, 2008; 5(2): 154 - 160. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Sumi, K. Chin, K. Takahashi, T. Nakamura, H. Matsumoto, A. Niimi, and M. Mishima Effect of nCPAP therapy on heart rate in patients with obstructive sleep apnoea-hypopnoea QJM, August 1, 2006; 99(8): 545 - 553. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Vazir, M. Dayer, P. C. Hastings, H. F. McIntyre, M. Y. Henein, P. A. Poole-Wilson, M. R. Cowie, M. J. Morrell, and A. K. Simonds Can heart rate variation rule out sleep-disordered breathing in heart failure? Eur. Respir. J., March 1, 2006; 27(3): 571 - 577. [Abstract] [Full Text] [PDF] |
||||
![]() |
M M Shoukri, M H Asyali, and A Donner Sample size requirements for the design of reliability study: review and new results Statistical Methods in Medical Research, August 1, 2004; 13(4): 251 - 271. [Abstract] [PDF] |
||||
![]() |
F. Roche, V. Pichot, E. Sforza, I. Court-Fortune, D. Duverney, F. Costes, M. Garet, and J-C. Barthelemy Predicting sleep apnoea syndrome from heart period: a time-frequency wavelet analysis Eur. Respir. J., December 1, 2003; 22(6): 937 - 942. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. W. Flemons, M. R. Littner, J. A. Rowley, P. Gay, W. M. Anderson, D. W. Hudgel, R. D. McEvoy, and D. I. Loube Home Diagnosis of Sleep Apnea: A Systematic Review of the Literature: An Evidence Review Cosponsored by the American Academy of Sleep Medicine, the American College of Chest Physicians, and the American Thoracic Society Chest, October 1, 2003; 124(4): 1543 - 1579. [Full Text] [PDF] |
||||
![]() |
C. Zamarron, F. Gude, J. Barcala, J. R. Rodriguez, and P. V. Romero Utility of Oxygen Saturation and Heart Rate Spectral Analysis Obtained From Pulse Oximetric Recordings in the Diagnosis of Sleep Apnea Syndrome Chest, May 1, 2003; 123(5): 1567 - 1576. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A. Nelesen, H. Yu, M. G. Ziegler, P. J. Mills, J. L. Clausen, and J. E. Dimsdale Continuous Positive Airway Pressure Normalizes Cardiac Autonomic and Hemodynamic Responses to a Laboratory Stressor in Apneic Patients Chest, April 1, 2001; 119(4): 1092 - 1101. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Crasset, S. Mezzetti, M. Antoine, P. Linkowski, J. P. Degaute, and P. van de Borne Effects of Aging and Cardiac Denervation on Heart Rate Variability During Sleep Circulation, January 2, 2001; 103(1): 84 - 88. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 1999 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |