To the Editors:
Adults with obstructive sleep apnoea (OSA) display substantial heart rate changes during sleep, associated with obstructive respiratory events 1, 2. These heart rate changes are significantly attenuated or abolished with effective continuous positive airway pressure (CPAP) treatment 3.
Heart rate variability (HRV) refers to the regulation of the sinoatrial node by the sympathetic and parasympathetic branches of the autonomic nervous system. Evaluation of autonomic heart rate control by means of HRV measures has been widely accepted as an important source of information on cardiovascular function 4, 5. HRV is important because it provides a window to observe the heart's ability to respond to normal or abnormal regulatory impulses that affect its rate and rhythm.
Several sophisticated methodologies were developed for measuring HRV, notably those using the heart rate power spectrum and nonlinear dynamics analysis 6. However, based on the literature, no simple methodology is available for comprehensive use in daily medical practice.
However, a Poincaré or Lorenz plot is a simple semi-quantitative method for visualising and calculating HRV 7–10. It has been previously shown that the length and width of a Poincaré plot are a weighted combination of low- and high-frequency powers, so it provides a theoretical link between frequency-domain spectral analysis techniques and time-domain Poincaré plot analysis 7.
Despite its advantages, a Poincaré plot has never been used during sleep studies in order to assess HRV or as a marker of continuous positive airway pressure (CPAP) effectiveness in evaluating autonomic heart changes in OSA patients.
The goal of this preliminary study was to propose a simple qualitative and quantitative methodology for evaluation of autonomic heart rate control using heart rate Poincaré plots during wakefulness and in each sleep stage (non-rapid eye movement (NREM) N1–N2, N3 and rapid eye movement sleep), in patients with OSA before and after CPAP treatment.
Our study included 25 patients (18 males and seven females, aged 33.6±11 yrs, body mass index 30±5.2 kg·m−2) with moderate-severe OSA. Five patients had hypertension (four patients were undergoing treatment with angiotensin converting enzyme (ACE) inhibitor or ACE inhibitor/diuretic combination), two patients were diabetic and one had hypothyroidism. None of the patients had a history of cardiac arrhythmias. Attended overnight polysomnography (Somnostar; Viasys, Hoechberg, Germany) was performed before and after CPAP treatment. During sleep, the apnoea-hypopnoea index, desaturation index and mean oxygen saturation level before and after CPAP treatment were 38.6±19.4 versus 2.8±1.5 events·h−1 (p<0.0001), 37.3±23.9 versus 2.5±1.7 desaturations ≥4% per hour (p<0.0001), and 91±3.6 versus 94.8±2.6% (p<0.0001), respectively. Poincaré plots of sequential R-R intervals (beat to beat variability) were assessed during a 7-h sleep period. We analysed plots both visually and quantitatively using scatter plot length, width and area (calculated from length and width). When estimating the bulk of the plot we excluded obvious isolated outliers, which mainly represent artefacts (fig. 1⇓).
Poincaré plot analysis. ········: maximal heart rate (HR) frequency (minimal RR interval); - - - -: minimal HR frequency (maximal RR interval); ΔRRr: difference on plot diagonal between minimal (RRmin) and maximal (RRmax) RR values; ΔRRt: maximal HR variability as maximal width-difference between of two points at parallel tangential lines determining plot. The square of the plot represents overall HR variability. Isolated outliers (arrows) represent artefacts.
We used the Shapiro–Wilk test to gather evidence about the “non-normality” of the sample. As differences between pairs tended to normality, we used the more sensitive paired t-test to compare the heart rate response before and after CPAP application. One way ANOVA with post hoc Scheffe multiple comparisons was used to evaluate significant differences between Poincaré plot parameters in various sleep stages and wakefulness before CPAP application.
Visual analysis showed a torpedo, fan shaped or complex shaped plots (complex referring to clustering of data points) during sleep stages in OSA patients and mainly comet shaped plots during wakefulness.
ΔRRt (width of the cloud referring to the maximal heart rate variability) of sleep stages NREM N1–2 was found to be different from ΔRRt sleep stage N3 (p = 0.0024) and from wakefulness (p = 0.0046), as well as area plot sleep stage N1–2 versus N3, p = 0.0062 but not from wakefulness (p = 0.077).
Statistically significant differences were observed in the square of the plot (referring to the overall HRV) in sleep stages NREM N1–2 before and after CPAP therapy (241.7 ms2×10−3±22.3 versus 79.8 ms2× 10−3±34.7, respectively) (p = 0.001). ΔRRt showed statistically significant differences in all sleep stages before and after CPAP therapy (fig. 2⇓). The fact that area plots after CPAP in NREM sleep stages 1–2 (data not shown) were found relatively more decreased in comparison to the respective width of the plot (fig. 2⇓), is indirect evidence of a diminished low frequency/high frequency ratio, i.e. a reduction of sympathetic heart activity particularly evident in NREM sleep stages 1–2 (since the area is the product of length and width).
Width of Poincaré dispersed cloud (maximal heart rate variability as maximal width-difference between of two points at parallel tangential lines determining plot) in various sleep stages and wakefulness, before (▪) and after (□) the application of continuous positive airway pressure therapy. NREM: nonrapid eye movement; REM: rapid eye movement. #: p = 0.01; ¶: p = 0.049; +: p = 0.044; §: p = 0.58.
Our study is the first to use the Poincaré plot method in patients with OSA before and after CPAP treatment. Autonomic disorganisation due to repetitive bradytachycardia events, as a result of obstructive events, is a common feature in OSA patients and the resolution of this phenomenon is a marker of effective CPAP treatment. HRV data were used in research studies, although not in daily clinical practice due to the complexity of most of the methods used. Poincaré plot analysis is a simple and possibly superior methodology compared to the classical methods that assess HRV. The square of the plot representing the overall HRV is easy to calculate and does not need special knowledge or training. This method might be applied to clinical practice in OSA patients for purposes such as evaluation of abnormalities in autonomic control or as a practical way to assess the effectiveness of CPAP treatment more specifically as a marker of sympathetic activity reduction. As HRV and the autonomic heart response are an important issue in OSA patients, the influence of CPAP therapy needs to be further assessed.
Statement of interest
None declared.
- © ERS Journals Ltd

















