Original Research ARTICLE
Reproducibility of Heart Rate Variability Is Parameter and Sleep Stage Dependent
- 1Preventive Cardiology and Sports Medicine, University Clinic for Cardiology, Bern University Hospital (Inselspital), University of Bern, Bern, Switzerland
- 2Sensory-Motor Systems Lab, Institute of Robotics and Intelligent Systems, ETH Zurich, Zurich, Switzerland
- 3Medical Faculty, University of Zurich, Zurich, Switzerland
- 4Institute of Pharmacology and Toxicology, Zurich Center for Interdisciplinary Sleep Research and Zurich Center for Integrative Human Physiology, University of Zurich, Zurich, Switzerland
Objective: Measurements of heart rate variability (HRV) during sleep have become increasingly popular as sleep could provide an optimal state for HRV assessments. While sleep stages have been reported to affect HRV, the effect of sleep stages on the variance of HRV parameters were hardly investigated. We aimed to assess the variance of HRV parameters during the different sleep stages. Further, we tested the accuracy of an algorithm using HRV to identify a 5-min segment within an episode of slow wave sleep (SWS, deep sleep).
Methods: Polysomnographic (PSG) sleep recordings of 3 nights of 15 healthy young males were analyzed. Sleep was scored according to conventional criteria. HRV parameters of consecutive 5-min segments were analyzed within the different sleep stages. The total variance of HRV parameters was partitioned into between-subjects variance, between-nights variance, and between-segments variance and compared between the different sleep stages. Intra-class correlation coefficients of all HRV parameters were calculated for all sleep stages. To identify an SWS segment based on HRV, Pearson correlation coefficients of consecutive R-R intervals (rRR) of moving 5-min windows (20-s steps). The linear trend was removed from the rRR time series and the first segment with rRR values 0.1 units below the mean rRR for at least 10 min was identified. A 5-min segment was placed in the middle of such an identified segment and the corresponding sleep stage was used to assess the accuracy of the algorithm.
Results: Good reproducibility within and across nights was found for heart rate in all sleep stages and for high frequency (HF) power in SWS. Reproducibility of low frequency (LF) power and of LF/HF was poor in all sleep stages. Of all the 5-min segments selected based on HRV data, 87% were accurately located within SWS.
Conclusions: SWS, a stable state that, in contrast to waking, is unaffected by internal and external factors, is a reproducible state that allows reliable determination of heart rate, and HF power, and can satisfactorily be detected based on R-R intervals, without the need of full PSG. Sleep may not be an optimal condition to assess LF power and LF/HF power ratio.
There is growing interest for assessments of the activity of the autonomic nervous system (ANS) both in research and in applied settings. Few methods exist to measure activity of certain branches of the ANS, such as microneurography to assess muscle sympathetic nerve activity or determination of noradrenalin in the blood to indicate spill-over of the sympathetic nervous system (Grassi and Esler, 1999). HRV is a widely used indirect marker of the cardiac autonomic nervous system activity (CANA) and finds applications in different fields such as psychology (Laborde et al., 2017), sport science (Aubert et al., 2003; Buchheit et al., 2004; Bellenger et al., 2016; Plews et al., 2016) and medicine, namely prenatal diagnostics (Fairchild, 2013), sleep studies (Tobaldini et al., 2013; Dodds et al., 2017), and mortality risk assessment in diseased patients (Huikuri, 1995; Malik et al., 2016; Zhou et al., 2016). In contrast to the other methods, HRV offers a non-invasive and simple method to approximate the CANA.
A fundamental constraint of HRV measurements is the rather poor reproducibility that can be improved when measurements are highly standardized. Standardized measurements are usually performed during 2–5 min in supine position after resting for at least 10 min at the same time of the day (Camm, 1996). However, this time requirement may be problematic in situations where regular (daily) measurements are needed, such as monitoring fatigue and training adaptations in athletes (Plews et al., 2014). Further, resting HRV measurements may also be problematic in young children who cannot stay motionless and relaxed on demand.
Recently, studies have analyzed HRV during sleep to investigate effects of acute stress (Hall et al., 2004; Hynynen et al., 2011), chronic stress (Hynynen et al., 2006), physical activity (Herzig et al., 2017), or diseases (Roumelioti et al., 2010; Amra et al., 2017) on the CANA. Indeed, sleep could provide a highly standardized condition to time-efficiently measure HRV. However, sleep architecture has been shown to affect the CANA in various ways (Somers et al., 1993; Silvani and Dampney, 2013), and also HRV parameters have been shown to vary between different sleep stages and as a function of circadian phase (Trinder et al., 2001; Busek et al., 2005; Vandewalle et al., 2007; Boudreau et al., 2013). Compared to wakefulness, sleep is characterized by a generalized cardiovascular deactivation and a resetting of baroreflex sensitivity (Silvani and Dampney, 2013). Rapid-eye movement (REM) sleep exhibits relatively high muscle sympathetic nerve activity, high lumbar sympathetic nerve activity, and low muscle tone (Silvani and Dampney, 2013). During REM sleep, levels of sympathetic nerve activity above wakefulness values have been recorded when bursts in sympathetic nerve activity induce blood pressure surges and sudden increases in heart rate (Somers et al., 1993). During non-rapid-eye movement (NREM) sleep, low sympathetic activity (Somers et al., 1993; Silvani and Dampney, 2013), reduced cardiac output and lower blood pressure (Somers et al., 1993) have been observed. In sleep stage 2, the occurrence of arousal stimuli accompanied by bursts of sympathetic nerve activity induce a transient increase in blood pressure and heart rate (Trinder et al., 2003). The amplitude and frequency of sympathetic-bursts and blood pressure surges are lower in slow-wave sleep (SWS, N3) than in stage 2 (Somers et al., 1993). SWS is a standardized state with reduced blood pressure variability (Silvani, 2008), constant autonomic activity and a regular breathing frequency, undisturbed by emotional stimuli (Murali et al., 2003). Hence, the absence of sympathetic activity bursts and stationarity of heart rate during SWS offers a highly standardized condition for HRV assessment (Brandenberger et al., 2005). While measurements during SWS are appealing, their validity and usefulness have only sparsely been investigated. In particular, the reliability of HRV measurements has, to our knowledge, never been investigated during specific sleep stages.
Determination of SWS usually requires PSG recordings. However, previous studies have shown close correspondence of HRV parameters and cortical activity (Charloux et al., 1998; Otzenberger et al., 1998; Brandenberger et al., 2001; Dumont et al., 2004). Further, specific HRV characteristics in SWS have been observed, with stationary heart rate and uncorrelated consecutive R-R intervals (Herzig et al., 2016). Recent studies have attempted to classify sleep stages using ECG parameters only (Ebrahimi et al., 2013; Long et al., 2014, 2017; Fonseca et al., 2015; Yoon H. et al., 2017; Yoon H. N. et al., 2017). Moreover, other studies have used a simple method to identify a segment within SWS and analyze HRV (Al Haddad et al., 2009; Herzig et al., 2016, 2017). However, the reliability of these latter methods, relying solely on HRV parameters to select a segment located in a SWS phase, without necessitating PSG recordings, have never been investigated. This method would allow simple HRV measurements in a standardized state that is least influenced by environmental factors.
The aims of this study were to assess the reproducibility of HRV parameters during the different sleep stages by quantifying the variances of HRV parameters between and within subjects, as well as between and within nights of individuals. Furthermore, using the information of PSG sleep staging, we verified the reliability of an algorithm based solely on HRV data to identify a single 5-min SWS segment within each night.
Healthy young males were recruited for a study investigating the effect of vestibular stimulation by a rocking bed on sleep and the sleep EEG (unpublished data). The study was approved by the Institutional Review Board of the Swiss Federal Institute of Technology in Zurich. Written informed consent was provided by all subjects and the study was conducted in accordance with the Declaration of Helsinki.
In the original study, the effect of two different bed rocking protocols on sleep onset was investigated. After a familiarization night, three nights were recorded in all subjects, two with rocking and one with a control condition (randomized cross-over design). The bed rocking protocols were: no rocking, rocking from lights off until sleep onset and rocking from before sleep onset until the end of the first 2 h of sleep. Due to the absence of a significant treatment effect on the mean and variance of any of the analyzed parameters (sleep onset, HRV parameters, heart rate), the data were pooled for the present analyses.
All participants underwent a screening night (polysomnography, PSG) to rule out sleep disorders prior to the study. They were normal sleepers (habitual sleep duration ~8 h) with moderate alcohol and caffeine consumption. Participants had to adhere to regular bed times for 7 days prior to the first experimental night and throughout the entire experimental phase (3 weeks in total). Furthermore, they had to abstain from caffeine and alcohol consumption during the 3 days prior to each experimental night. Compliance with the bedtimes was assessed by actigraphy and sleep logs. Details on the methodology will be published elsewhere.
PSG Recordings/Sleep Staging
PSG data (EEG, EOG, EMG, ECG, and respiration) were continuously recorded throughout the entire 8-h sleep period with a polygraphic amplifier Artisan (Micromed, Mogliano, Veneto, Italy). The signals were sampled at 256 Hz and recorded with the software Rembrandt DataLab (Version 8.0; Embla Systems, Broomfield, CO, USA). The analog signals were filtered with a high pass filter (EEG: −3 dB at 0.16 Hz; EMG: 10 Hz; ECG: 1 Hz) and an anti-aliasing low-pass filter (−3 dB at 67.2 Hz). For further analysis, the EEG signals were re-referenced to the mastoids (A1, A2). The sleep stages were scored visually on a 20-s epoch basis as suggested by Rechtschaffen and Kales (Rechtschaffen and Kales, 1968) according to standard criteria (Iber et al., 2007). For artifact removal, artifacts were identified visually and with a semi-automatic artifact detection (see Lustenberger et al., 2012 for details).
R-peaks were automatically detected from the ECG trace using the ≪nqrsdetect≫ function from the Biosig Matlab (2014a, The Mathworks, Natick, MA) Toolbox (Vidaurre et al., 2011) and inter-beat durations (R-R intervals) were calculated. R-R intervals were analyzed using a Matlab procedure specifically developed for this study. After a stage change, only consecutive 5-min segments consisting of a particular stage were used for analysis (Figure 1). Due to the often reported time delay of varying length between CANA and cortical activity (Jurysta et al., 2003; Long et al., 2015) (with CANA preceding cortical activity), the last 5-min segment of each sleep phase was discarded. Thus, only sleep phases longer than 10 min were included in the analysis. In each 5-min segment, frequency domain and time domain parameters were calculated. Mean heart rate [beats per minute, bpm] was derived from the mean of the R-R intervals using the following formula: heart rate [bpm] = 60,000/(mean R-R intervals [ms]). For spectral analysis, R-R intervals were interpolated (cubic spline interpolation) and resampled at 4 Hz. We applied an advanced smoothness prior approach for detrending of the R-R intervals with a smoothing parameter of λ = 500, which corresponds to a high pass filter with a cut-off frequency of 0.035 Hz (Camm, 1996). We used an artifact correction algorithm that eliminates R-R intervals in case of deviations of 30% or more to adjacent R-R intervals and replaced them using a cubic-spline interpolation. Power spectral density was then calculated using the Welch method (1967) with a Hamming window length of 128 points and an overlap of 50%. Frequency domain parameters were total power (TP, ms2, 0–0.4 Hz), low-frequency power (LF, ms2, 0.04–0.15 Hz), high-frequency power (HF, ms2 0.15–0.4 Hz), and the LF/HF ratio. The following time domain parameters were calculated: the square root of the mean square differences of adjacent R-R intervals (RMSSD, ms) and the standard deviation of all R-R intervals (SDNN, ms). Normalized frequency parameters are not reported in order to avoid redundancy with LF/HF ratio (Massin et al., 1999; Burr, 2007).
Figure 1. Illustration of the selection of the segments for HRV analysis (top). Hypnogram (W, waking; R, REM sleep; 1 to 3, NREM sleep stages N1 to N3). The segments selected for the analysis are indicated by the green boxes (indicated for visual purpose only and are not to scale). Note that the last 5-min segments of a continuous phase of a particular sleep stage were discarded. Time course of the correlation of consecutive R-R intervals (rRR, bottom) calculated with a 5-min windows moved in steps of 20 s. The dotted line represents the cut-off for determining periods of SWS (0.1 unit below the mean of the first 4-h). The identified SWS phases are highlighted in gray. Note that only one 5-min segment in the middle of the first detected SWS phase was used for our analyses.
SWS Segment Identification Using HRV (Segmentlow rRR)
By means of a custom built Matlab script, Pearson's correlation coefficient in the Poincaré plot of the consecutive R-R intervals (rRR, Brennan et al., 2001) were calculated of 5-min windows moved in steps of 20 s over the whole night. The linear trend of rRR time series in the first 4 h was removed and we identified the first period where rRR values were 0.1 units below the mean rRR (which was 0 due to the detrending) of the first 4 h for at least 10 min, the same threshold has been used in a previous study (Herzig et al., 2016). A 5-min segment (Segmentlow rRR) was then put in the middle of the identified period (mid-point of 5-min segment allocated to the midpoint of the identified period) of which frequency and time domain analysis of the R-R intervals were performed. Due to the detrending procedure all segments were placed in the first 4 h of the night.
Verification of the Identified SWS Segment (Segmentlow rRR)
For each night, the Segmentlow rRR (identified by means of the HRV signal) was classified as accurately detected if the complete 5-min segment was located in a phase of SWS (as classified according to the PSG recordings). Further, segments consisting of at least 50% SWS were counted and reported separately.
In order to quantify the variances of heart rate and HRV parameters of the different sleep stages, segments were grouped according to sleep stage as identified by PSG. The study design with three nights within each subject and several segments of different sleep stages within each night allowed the following variance partitioning. Total variance was partitioned into: Between-subjects variance (Varbetween Sub), between-nights variance (Varbetween Nights), and between-segments variance (Varbetween Segments) with total variance being equal to the sum of theses variances (Total Variance = VarBetween Sub + VarBetween Nights + VarBetween Segments). Further, the variance components were expressed as proportions of total variance as well as of the within-subjects variance (Within-Subjects Variance = VarBetween Nights + VarBetween Segments). For this purpose, nested random effect models were applied using the R software package lme4 (Bates et al., 2015) to quantify the variances of interest. Heart rate, HF power, LF power, LF/HF ratio, RMSSD, and SDNN were entered as dependent variables. Standard deviations () were calculated and are reported as a measure of absolute variability.
The between-subjects/within-subjects intra-class correlation coefficient (ICC) was estimated as the proportion of the between-subjects variance to the total variance (ICC = VarBetween Sub/Total Variance) which is equivalent to using a two-way mixed model analysis of variance and the ICC(3,1) formula as defined by Shrout and Fleiss (Shrout and Fleiss, 1979). For the within-subject reliability, the between-nights/within-nights ICC was calculated (ICC = VarBetween Nights/(VarBetween Nights + VarianceBetween Segments)). ICCs were determined for heart rate and all HRV parameters in all sleep stages. Between-nights/within-nights ICC, providing information on whether different nights of the same subject can be distinguished, corresponds to the proportion of the between-nights variance of the within-subject variance. Interpretation of ICCs was based on the classification by Cicchetti (1994).
Statistical analyses were performed using the software R (Version 3.3.1, R Core Team, 2016). Models were tested with regard to satisfaction of underlying statistical assumptions such as normal distribution of residuals and homoscedasticity. Non-normally distributed variables were log transformed and model assumptions were tested again.
Eighteen healthy male participants with mean age 23.7 ± 2.5 (SD) years were included in the study. Data of three subjects had to be excluded. Reasons were insufficient quality of the ECG recording of two subjects and availability of only two nights in another subject. Thus, 15 subjects with 3 recordings each were included in the analysis. Median values of HRV parameters of separate sleep stages are shown in Table 1. While heart rate was highest during REM sleep, comparable heart rates were found in stage 2 and SWS. Similarly, HF power in stage 2 and SWS was lower than during REM sleep, however, differences were small. The lowest values of LF power, RMSSD, and SDNN were observed in SWS, with low LF power also resulting in a low LF/HF ratio.
HRV Segments during Sleep
For the 45 nights included in the analysis, a total of 1792 5-min segments were analyzed (14 in sleep stage 1, 730 in sleep stage 2, 339 SWS and 709 in REM sleep). The low abundance of sleep stage 1 and the generally short duration of this stage resulted in this low number of segments, therefore, only data from sleep stage 2, REM sleep, and SWS were included for further analyses. Mean numbers of analyzed segments per night was 16 for stage 2 (range 3–40), 8 for SWS (range 3–16), and 16 for REM sleep (range 3–24). Of the segments included in the variance analysis, stage 2 segments occurred 265 [IQR 189, 350] min after sleep onset for stage 2, SWS segments 101 [33, 140] min after sleep onset and REM sleep segments after 325 [222, 408] min. Raw data of HF power and LF power over the course of the night in the different sleep stages are presented in the supplementary material (Figure S1).
Effect of Order of the Night
In our sample, there was no evidence of an effect of the order of the night on any of the analyzed parameters (all p > 0.15).
Variance Partitioning of HRV Parameters in the Different Sleep Stages
The between subjects, between nights (within each subject) and between segments (within each night) SDs for each individual sleep stage are reported in Table 2. For the calculated variances, the proportion of the total variance as well as of the within-subject variance are reported in Table 3.
Table 3. Variance decomposition of total variance and within subject variance, as well as intra-class correlation coefficients (indicated in brackets) of the different sleep stages.
For all HRV parameters in all sleep stages and all HRV parameters, with exception of LF/HF ratio and SDNN in stage 2 and in REM sleep, the largest proportion of total variance was explained by the between-subject variance (84% for HF power and RMSSD in SWS). The between segment variances were small in SWS, in particular for heart rate, HF power, and RMSSD, where the between segment variances were between 8 and 12% of the total variance. On the other hand, the between-segment variance was around 40% of the total variance for LF power and LF/HF ratio in all sleep stages and >56% for SDNN in stage 2 and REM.
The between-subjects/within-subjects ICC is equivalent to the percentage of explained between-subject variance (Table 3). The between-nights/within-nights ICC is equivalent to the between-night variance of the within-subject variance (Table 3). In accordance to between subject variances, between-subjects/within-subjects ICCs were good for heart rate, HF power, and RMSSD in SWS. Good to excellent between-subjects/within-subjects ICCs were observed in all sleep stages for heart rate and RMSSD with the highest values in SWS. For LF power and LF/HF ratio, fair ICCs were observed with comparable values in all analyzed sleep stages. Between-nights/within-nights ICCs were generally poor, except for heart rate, where it was good for SWS and fair for stage 2 and REM sleep.
Verification of the Identified SWS Segment (SegmentLow rRR)
Of all analyzed nights, 87% (39 out of 45) of the detected segments were fully located within SWS and where thus correctly identified. Of the remaining 6 segments, 2 were partially located during SWS (>50%), 2 segments were located completely in stage 2 and 2 segments were located across stage 1 and stage 2. An example of the sleep staging of a typical night and the rRR, used to identify the Segmentlow rRR, of the same night are shown in Figure 1.
In the present study aiming at characterizing the reproducibility of HRV variables in the different sleep stages, we found better reproducibility (between-subject/within subject ICC and SDs) for HF power, RMSSD, and heart rate during SWS compared to the other sleep stages. In all sleep stages, variances between subjects generally accounted for the highest percentage of total variance, followed by between-segment variance (within same nights) and between nights within each subject. Between-subject/within subject ICCs were excellent for heart rate in all sleep stages. They were also excellent for HF power and RMSSD and good for SDNN in SWS. Fair ICCs were found for LF and LF/HF ratio in all sleep stages. When a single 5-min segment was located within SWS identified by an algorithm based on HRV data only, 87% were correctly placed within SWS according to PSG.
The variances observed during SWS were comparable to the variances reported by Schroeder et al. (2004) for HRV measurements in wake supine position, using a repeated measure study design similar to ours. The lowest residual variance (within-night variance), reflecting the best reproducibility, was found for heart rate and HF power in SWS, followed by stage 2, and highest residual variance and poorest reproducibility in REM sleep. These results are in accordance with the low cardiovascular variability reported during SWS (Franzini, 2000). Thus, while heart rate, HF power and RMSSD in SWS are highly reproducible under comparable conditions, it remains open whether it is also a stage suited to detect acute or chronic stress. Acute psychological stress was found to be reflected by lower HF power in stress situations both during NREM and REM sleep, and of note, with the difference increasing over the course of the night (Hall et al., 2004). Further research is required to assess how well HRV parameters measured in the different sleep stages may reflect the effects of different stressors.
This is, to the best of our knowledge, the first study to quantify variance components of HRV parameters as well as reproducibility in different sleep stages. Reproducibility of HRV measurements in an awake state between different days has been extensively studied (Schroeder et al., 2004; Pinna et al., 2007; Sacre et al., 2012; Schafer et al., 2015; Silva et al., 2017). During sleep, reproducibility of mean heart rate of a whole night's sleep has been investigated (Waldeck and Lambert, 2003). In these studies, ICC, a measure for relative reproducibility, for HRV measurements between different days ranged from 0.64 to 0.91 when measured in a wake state in supine position. Similar values of reproducibility were found for different HRV parameters in the studies cited above. In our study, ICCs were best for HF power (0.84), RMSSD (0.84), and SDNN (0.73) in SWS, and poorer for the two other sleep stages. For heart rate ICCs were between 0.75 and 0.77 in all sleep stages. For LF power and LF/HF power ratio ICCs were below 0.55 in all sleep stages. The low ICC for LF power and LF/HF ratio during SWS may have been a consequence of the generally very low values of LF power in all subjects, reflecting low sympathetic activity (Somers et al., 1993; Silvani and Dampney, 2013) during deep sleep. Greater heterogeneity of the population regarding the parameter of interest automatically leads to greater ICCs compared to more homogeneous populations. The generally poorer ICCs in stage 2 and REM sleep for all analyzed HRV parameters may be explained by the higher frequency and higher amplitude of sympathetic surges and therefore the lower stationarity of HRV parameters in these two stages.
Based on the results of our study, it seems advisable to assess HF power and/or RMSSD during SWS rather than another sleep stage due to a smaller between segment variance. Given that the variance was greater between segments than between nights (and this was not due to a trend over the course of the night), it seems important to average over several segments within one night. Our findings of a low between-subject/within-subject ICC for LF power and LF/HF ratio in all sleep stages questions the suitability of sleep in general for the assessment of these parameters, which are often used as marker of sympathetic activity.
In previous studies assessing the effects of sleep stages on heart rate and HRV, the highest values of heart rate were observed during REM sleep and lower heart rates during stage 2 and SWS with no significant difference between the latter two stages (Trinder et al., 2001; Boudreau et al., 2013). Similar effects of sleep stages on heart rate were observed in our study. HF power was comparable in SWS, stage 2 and REM sleep which was also the case in previous studies (Busek et al., 2005; Boudreau et al., 2013) while others have observed a reduction in HF power in REM sleep compared to the other sleep stages (Trinder et al., 2001). In our study, LF power was lowest during SWS and highest during REM sleep. This is in accordance with results from other studies and indicates a clear sleep stage dependence of LF power (Trinder et al., 2001; Busek et al., 2005; Cabiddu et al., 2012; Boudreau et al., 2013). Reported differences between sleep stages with regard to absolute LF power (up to 140%) in our and previous studies were much bigger than with regard to absolute HF power (<25%). The large differences in LF power between different sleep stages resulted in a concomitantly large difference in LF/HF ratio, being lowest in SWS, followed by stage 2 and being highest during REM sleep. These findings are also in accordance with previous studies (Scholz et al., 1997; Elsenbruch et al., 1999; Ferri et al., 2000).
Many studies reported a rapid decrease in heart rate concomitant with the wake-sleep transition (Shinar et al., 2006) and a more gradual decrease thereafter throughout the night (Burgess et al., 1999). This has been clearly shown when the different sleep stages were analyzed separately (Cajochen et al., 1994; Trinder et al., 2001). When heart rate is analyzed irrespective of the sleep stages, then often an increase in heart rate toward the morning hours was observed. Reasons for this may be the higher abundance of REM sleep in the later part of sleep, the circadian increase in cortisol concentration in the blood and/or the increase in body temperature toward the morning (Burgess et al., 1997). Further, heart rate has been found to steeply increase at the end of each SWS phase, when an arousal and the concomitant blood pressure surge initiates a transition into REM sleep (Bonnet and Arand, 1997). These autonomic processes precede the central activities according to which the sleep stages are classified. Unexpectedly, in our study, a mean increase in heart rate of 3.2 ± 0.4 (SD) bpm was observed in the first SWS phase. This was followed by a mean decrease in heart rate of (2.9 ± 0.5 bpm) over the rest of the night. The initial increase in heart rate in the first SWS phase has not been reported previously, however, some data seems to be in agreement with our finding (Cajochen et al., 1994; Burgess et al., 1999; Brandenberger et al., 2005). This acceleration of heart rate in the first SWS phase was not accompanied by a concomitant decrease in HF power as would be expected when caused by a decreased vagal activity or by an increased sympathetic activity. We suspect factors other than the CANA to be responsible for this discrepancy.
The accuracy of our algorithm for the detection of a segment in SWS by relying solely on HRV data was satisfactory with most of the segments (87%) being accurately placed within a phase of SWS (as verified by PSG). Different HRV based algorithms to perform sleep staging have been proposed and compared to sleep staging by PSG. Most algorithms have used HRV frequency domain parameters (Shinar et al., 2001; Mendez et al., 2009). Other studies have used a combination of different HRV parameters (time domain, frequency domain, non-linear parameters) with the aim to improve the precision of their classification algorithm (Xiao et al., 2013; Fonseca et al., 2015; Yoon H. et al., 2017; Yoon H. N. et al., 2017). Compared to the traditional sleep staging by PSG, an accuracy of 69% in the classification of sleep stages was achieved in the study by Fonseca et al. (2015). In an attempt to classify sleep into wake-REM-NREM sleep, an accuracy of 79% has been reported using a combination of frequency domain parameters (Mendez et al., 2009). In one early study (Shinar et al., 2001), the use of HRV (LF/HF ratio) to identify SWS resulted in 80% correct identification of SWS. Recently, using HRV parameters besides respiratory signals measured by respiratory inductance plethysmography, an accuracy of 89% was reported for the detection of slow wave sleep (Long et al., 2017). Similarly, using a combination of time domain, frequency domain and non-linear parameters resulted in an Cohens kappa of 0.56 for the detection of SWS (Yoon H. N. et al., 2017). Combining HRV frequency domain parameters with actigraphy resulted in an accuracy of 74% (Muzet et al., 2016), comparable to the methods based on HRV only. In contrast to these studies, our aim was not to perform an actual sleep staging but to verify that a segment identified by HRV analysis was located within SWS This explains the rather high accuracy of our algorithm, based on rRR values, a simple HRV measure. In one of our previous studies in elite athletes (Herzig et al., 2016), the same algorithm resulted in a yield of 75% of segments placed in SWS according to non-conventional PSG sleep staging with a mobile 3-lead EEG system.
A limitation of the present study was the absence of a direct measurement of the CANA. However, the aim of this study was to quantify the reproducibility of HRV parameters without making direct inferences on the CANA. Further, we have only quantified reproducibility under standardized conditions, it will be the task of future studies to assess the sensitivity of HRV parameters in different sleep stages to detect disturbances, such as psychological or physical acute or chronic stress. While cortical activity largely follows autonomic activity (Jurysta et al., 2003) some dissociation may occur. Further, data of the present study was collected for a different study assessing a treatment effect. However, the small between nights variance support the absence of a treatment or night order effect. Last but not least, the present study was conducted in healthy subjects, and findings cannot be extrapolated to diseased populations. Strengths of the present study are the statistical analysis of the nested design separating variances due to different factors. This is the first study to assess the reproducibility of HRV parameters during sleep and to compare reproducibility of these parameters between the different sleep stages. The simple algorithm to identify SWS used in the present study will allow wide usage of this efficient method to measure HRV parameters in SWS.
In conclusion, SWS offers a stationary and reproducible phase for determination of heart rate, HF power and RMSSD, and can satisfactorily be detected based on R-R intervals derived from ECG recordings, without the need of full PSG. LF power and LF/HF ratio had fair to poor reproducibility in all sleep stages. Between segment variance was high compared to between night variance for all HRV parameters except heart rate, indicating that averaging over segments within each night is advisable.
XO, RR, and PA designed the experiment and conducted the data collection. DH, PE, MW, XO, and PA were involved in data analysis. DH and PE composed the manuscript. All authors reviewed and approved the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The reviewer MWR and handling Editor declared their shared affiliation.
We would like to thank all subjects who participated in the study. The study was supported by the ETH Zurich Research Grant ETHIIRA (ETH-18 11-1), the NCCR Transfer Projects of the Swiss National Science Foundation (51NF40-1444639), and a grant from the Swiss National Foundation (32003B_146643).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2017.01100/full#supplementary-material
Figure S1. HF power (Left) and LF power (Right) of 5-min segments over the course of the night of all subjects and all nights. The different colors represent data of different subjects and segments corresponding to the same night are connected by a solid line. The segments of the pooled sleep stages are shown in the top panel, time course of the segments of the different sleep stages are shown in the lower panels. Time 0 corresponds to sleep onset.
Al Haddad, H., Laursen, P. B., Ahmaidi, S., and Buchheit, M. (2009). Nocturnal heart rate variability following supramaximal intermittent exercise. Int. J. Sports Physiol. Perform. 4, 435–447. doi: 10.1123/ijspp.4.4.435
Amra, B., Behjati, M., Penzel, T., Fietze, I., Schoebel, C., and Sarrafzadegan, N. (2017). Nocturnal heart rate variation in diabetic and non-diabetic patients with sleep apnea syndrome. Sleep Med. 29, 57–60. doi: 10.1016/j.sleep.2016.11.003
Bellenger, C. R., Fuller, J. T., Thomson, R. L., Davison, K., Robertson, E. Y., and Buckley, J. D. (2016). Monitoring athletic training status through autonomic heart rate regulation: a systematic review and meta-analysis. Sports Med. 46, 1461–1486. doi: 10.1007/s40279-016-0484-2
Bonnet, M. H., and Arand, D. L. (1997). Heart rate variability: sleep stage, time of night, and arousal influences. Electroencephalogr. Clin. Neurophysiol. 102, 390–396. doi: 10.1016/S0921-884X(96)96070-1
Brandenberger, G., Buchheit, M., Ehrhart, J., Simon, C., and Piquard, F. (2005). Is slow wave sleep an appropriate recording condition for heart rate variability analysis? Auton. Neurosci. 121, 81–86. doi: 10.1016/j.autneu.2005.06.002
Brandenberger, G., Ehrhart, J., Piquard, F., and Simon, C. (2001). Inverse coupling between ultradian oscillations in delta wave activity and heart rate variability during sleep. Clin. Neurophysiol. 112, 992–996. doi: 10.1016/S1388-2457(01)00507-7
Brennan, M., Palaniswami, M., and Kamen, P. (2001). Do existing measures of Poincare plot geometry reflect nonlinear features of heart rate variability? IEEE Trans. Biomed. Eng. 48, 1342–1347. doi: 10.1109/10.959330
Buchheit, M., Simon, C., Piquard, F., Ehrhart, J., and Brandenberger, G. (2004). Effects of increased training load on vagal-related indexes of heart rate variability: a novel sleep approach. Am. J. Physiol. Heart Circ. Physiol. 287, H2813–H2818. doi: 10.1152/ajpheart.00490.2004
Burgess, H. J., Trinder, J., and Kim, Y. (1999). Cardiac autonomic nervous system activity during presleep wakefulness and stage 2 NREM sleep. J. Sleep Res. 8, 113–122. doi: 10.1046/j.1365-2869.1999.00149.x
Cabiddu, R., Cerutti, S., Viardot, G., Werner, S., and Bianchi, A. M. (2012). Modulation of the sympatho-vagal balance during sleep: frequency domain study of heart rate variability and respiration. Front. Physiol. 3:45. doi: 10.3389/fphys.2012.00045
Camm, J. (1996). Heart rate variability: standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Circulation 93, 1043–1065. doi: 10.1161/01.CIR.93.5.1043
Charloux, A., Otzenberger, H., Gronfier, C., Lonsdorfer-Wolf, E., Piquard, F., and Brandenberger, G. (1998). Oscillations in sympatho-vagal balance oppose variations in delta-wave activity and the associated renin release. J. Clin. Endocrinol. Metab. 83, 1523–1528.
Dodds, K. L., Miller, C. B., Kyle, S. D., Marshall, N. S., and Gordon, C. J. (2017). Heart rate variability in insomnia patients: a critical review of the literature. Sleep Med. Rev. 33, 88–100. doi: 10.1016/j.smrv.2016.06.004
Dumont, M., Jurysta, F., Lanquart, J. P., Migeotte, P. F., van de Borne, P., and Linkowski, P. (2004). Interdependency between heart rate variability and sleep EEG: linear/non-linear? Clin. Neurophysiol. 115, 2031–2040. doi: 10.1016/j.clinph.2004.04.007
Ebrahimi, F., Setarehdan, S. K., Ayala-Moyeda, J., and Nazeran, H. (2013). Automatic sleep staging using empirical mode decomposition, discrete wavelet transform, time-domain, and nonlinear dynamics features of heart rate variability signals. Comput. Methods Programs Biomed. 112, 47–57. doi: 10.1016/j.cmpb.2013.06.007
Ferri, R., Parrino, L., Smerieri, A., Terzano, M. G., Elia, M., Musumeci, S. A., et al. (2000). Cyclic alternating pattern and spectral analysis of heart rate variability during normal sleep. J. Sleep Res. 9, 13–18. doi: 10.1046/j.1365-2869.2000.00190.x
Fonseca, P., Long, X., Radha, M., Haakma, R., Aarts, R. M., and Rolink, J. (2015). Sleep stage classification with ECG and respiratory effort. Physiol. Meas. 36, 2027–2040. doi: 10.1088/0967-3334/36/10/2027
Hall, M., Vasko, R., Buysse, D., Ombao, H., Chen, Q., Cashmere, J. D., et al. (2004). Acute stress affects heart rate variability during sleep. Psychosom. Med. 66, 56–62. doi: 10.1097/01.PSY.0000106884.58744.09
Herzig, D., Eser, P., Radtke, T., Wenger, A., Rusterholz, T., Wilhelm, M., et al. (2017). Relation of heart rate and its variability during sleep with age, physical activity, and body composition in young children. Front. Physiol. 8:109. doi: 10.3389/fphys.2017.00109
Herzig, D., Testorelli, M., Olstad, D. S., Erlacher, D., Achermann, P., Eser, P., et al. (2016). Heart rate variability during deep sleep offers a time-efficient alternative to morning supine measurements - a study in world class alpine skiers. Int. J. Sports Physiol. Perform. 12, 648–654. doi: 10.1123/ijspp.2016-0257
Hynynen, E., Konttinen, N., Kinnunen, U., Kyrolainen, H., and Rusko, H. (2011). The incidence of stress symptoms and heart rate variability during sleep and orthostatic test. Eur. J. Appl. Physiol. 111, 733–741. doi: 10.1007/s00421-010-1698-x
Hynynen, E., Uusitalo, A., Konttinen, N., and Rusko, H. (2006). Heart rate variability during night sleep and after awakening in overtrained athletes. Med. Sci. Sports Exerc. 38, 313–317. doi: 10.1249/01.mss.0000184631.27641.b5
Iber, C., Ancoli-Israel, S., Chesson, A., and Quan, S. F. (2007). The AASM Manual for the Scoring of Sleep ans Associated Events: Rules, Terminology and Technical Specifications, 1st Edn. American Academy of Sleep Medicine.
Jurysta, F., van de Borne, P., Migeotte, P. F., Dumont, M., Lanquart, J. P., Degaute, J. P., et al. (2003). A study of the dynamic interactions between sleep EEG and heart rate variability in healthy young men. Clin. Neurophysiol. 114, 2146–2155. doi: 10.1016/S1388-2457(03)00215-3
Laborde, S., Mosley, E., and Thayer, J. F. (2017). Heart rate variability and cardiac vagal tone in psychophysiological research - recommendations for experiment planning, data analysis, and data reporting. Front. Psychol. 8, 213. doi: 10.3389/fpsyg.2017.00213
Long, X., Arends, J. B., Aarts, R. M., Haakma, R., Fonseca, P., and Rolink, J. (2015). Time delay between cardiac and brain activity during sleep transitions. Appl. Phys. Lett. 106, 143702. doi: 10.1063/1.4917221
Long, X., Fonseca, P., Aarts, R. M., Haakma, R., Rolink, J., and Leonhardt, S. (2017). Detection of Nocturnal slow wave sleep based on cardiorespiratory activity in healthy adults. IEEE J. Biomed. Health Inform. 21, 123–133. doi: 10.1109/JBHI.2015.2487446
Long, X., Fonseca, P., Haakma, R., Foussier, J., and Aarts, R. M. (2014). Automatic detection of overnight deep sleep based on heart rate variability: a preliminary study. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2014, 50–53. doi: 10.1109/EMBC.2014.6943526
Lustenberger, C., Maric, A., Durr, R., Achermann, P., and Huber, R. (2012). Triangular relationship between sleep spindle activity, general cognitive ability and the efficiency of declarative learning. PLoS ONE 7:e49561. doi: 10.1371/journal.pone.0049561
Massin, M. M., Derkenne, B., and von Bernuth, G. (1999). Correlations between indices of heart rate variability in healthy children and children with congenital heart disease. Cardiology 91, 109–113. doi: 10.1159/000006889
Mendez, M. O., Matteucci, M., Cerutti, S., Aletti, F., and Bianchi, A. M. (2009). Sleep staging classification based on HRV: time-variant analysis. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2009, 9–12. doi: 10.1109/IEMBS.2009.5332624
Muzet, A., Werner, S., Fuchs, G., Roth, T., Saoud, J. B., Viola, A. U., et al. (2016). Assessing sleep architecture and continuity measures through the analysis of heart rate and wrist movement recordings in healthy subjects: comparison with results based on polysomnography. Sleep Med. 21, 47–56. doi: 10.1016/j.sleep.2016.01.015
Otzenberger, H., Gronfier, C., Simon, C., Charloux, A., Ehrhart, J., Piquard, F., et al. (1998). Dynamic heart rate variability: a tool for exploring sympathovagal balance continuously during sleep in men. Am. J. Physiol. 275(3 Pt 2), H946–H950.
Pinna, G. D., Maestri, R., Torunski, A., Danilowicz-Szymanowicz, L., Szwoch, M., La Rovere, M. T., et al. (2007). Heart rate variability measures: a fresh look at reliability. Clin. Sci. 113, 131–140. doi: 10.1042/CS20070055
Plews, D. J., Laursen, P. B., and Buchheit, M. (2016). Day-to-day Heart Rate Variability (HRV) recordings in world champion rowers: appreciating unique athlete characteristics. Int. J. Sports Physiol. Perform. 12, 697–703. doi: 10.1123/ijspp.2016-0343
Plews, D. J., Laursen, P. B., Le Meur, Y., Hausswirth, C., Kilding, A. E., and Buchheit, M. (2014). Monitoring training with heart rate-variability: how much compliance is needed for valid assessment? Int. J. Sports Physiol. Perform. 9, 783–790. doi: 10.1123/ijspp.2013-0455
Rechtschaffen, A., and Kales, A. (1968). A Manual of Standardized Terminology, Techniques and Scoring System for Sleep Stages and Human Subjects. Los Angeles, CA: Brain Information Service/Brain Research Institute; University of California,
Roumelioti, M. E., Ranpuria, R., Hall, M., Hotchkiss, J. R., Chan, C. T., Unruh, M. L., et al. (2010). Abnormal nocturnal heart rate variability response among chronic kidney disease and dialysis patients during wakefulness and sleep. Nephrol Dial Transplant. 25, 3733–3741. doi: 10.1093/ndt/gfq234
Sacre, J. W., Jellis, C. L., Marwick, T. H., and Coombes, J. S. (2012). Reliability of heart rate variability in patients with type 2 diabetes mellitus. Diabet. Med. 29,e33–e40. doi: 10.1111/j.1464-5491.2011.03557.x
Schafer, D., Nil, M., Herzig, D., Eser, P., Saner, H., and Wilhelm, M. (2015). Good reproducibility of heart rate variability after orthostatic challenge in patients with a history of acute coronary syndrome. J. Electrocardiol. 48, 696–702. doi: 10.1016/j.jelectrocard.2015.04.004
Scholz, U. J., Bianchi, A. M., Cerutti, S., and Kubicki, S. (1997). Vegetative background of sleep: spectral analysis of the heart rate variability. Physiol. Behav. 62, 1037–1043. doi: 10.1016/S0031-9384(97)00234-5
Schroeder, E. B., Whitsel, E. A., Evans, G. W., Prineas, R. J., Chambless, L. E., and Heiss, G. (2004). Repeatability of heart rate variability measures. J. Electrocardiol. 37, 163–172. doi: 10.1016/j.jelectrocard.2004.04.004
Shinar, Z., Akselrod, S., Dagan, Y., and Baharav, A. (2006). Autonomic changes during wake-sleep transition: a heart rate variability based approach. Auton. Neurosci. 130, 17–27. doi: 10.1016/j.autneu.2006.04.006
Silva, C. C., Bertollo, M., Reichert, F. F., Boullosa, D. A., and Nakamura, F. Y. (2017). Reliability of heart rate variability in children: influence of sex and body position during data collection. Pediatr. Exerc. Sci. 29, 228–236. doi: 10.1123/pes.2016-0085
Silvani, A. (2008). Physiological sleep-dependent changes in arterial blood pressure: central autonomic commands and baroreflex control. Clin. Exp. Pharmacol. Physiol. 35, 987–994. doi: 10.1111/j.1440-1681.2008.04985.x
Tobaldini, E., Nobili, L., Strada, S., Casali, K. R., Braghiroli, A., and Montano, N. (2013). Heart rate variability in normal and pathological sleep. Front. Physiol. 4:294. doi: 10.3389/fphys.2013.00294
Trinder, J. Kleiman, J., Carrington, M., Smith, S., Breen, S., Tan, N., et al. (2001). Autonomic activity during human sleep as a function of time and sleep stage. J. Sleep Res. 10, 253–264. doi: 10.1046/j.1365-2869.2001.00263.x
Vandewalle, G., Middleton, B., Rajaratnam, S. M., Stone, B. M., Thorleifsdottir, B., Arendt, J., et al. (2007). Robust circadian rhythm in heart rate and its variability: influence of exogenous melatonin and photoperiod. J. Sleep Res. 16, 148–155. doi: 10.1111/j.1365-2869.2007.00581.x
Vidaurre, C., Sander, T. H., and Schlogl, A. (2011). BioSig: the free and open source software library for biomedical signal processing. Comput. Intell. Neurosci. 2011:935364. doi: 10.1155/2011/935364
Welch, P. D. (1967). The use of fast fourier transform for the estimation of power spectra: a method based on time averaging over short, modified periodograms. IEEE Trans. Audio Electroacoust. 15, 70–73. doi: 10.1109/TAU.1967.1161901
Xiao, M., Yan, H., Song, J., Yang, Y., and Yang, X. (2013). Sleep stages classification based on heart rate variability and random forest. Biomed. Signal Process. Control 8, 624–633. doi: 10.1016/j.bspc.2013.06.001
Yoon, H., Hwang, S. H., Choi, J. W., Lee, Y. J., Jeong, D. U., and Park, K. S. (2017). REM sleep estimation based on autonomic dynamics using R-R intervals. Physiol. Meas. 38, 631–651. doi: 10.1088/1361-6579/aa63c9
Yoon, H. N., Hwang, S. H., Choi, J. W., Lee, Y. J., Jeong, D. U., and Park, K. S. (2017). Slow-wave sleep estimation for healthy subjects and OSA patients using R-R intervals. IEEE J. Biomed. Health Inform. doi: 10.1109/JBHI.2017.2712861. [Epub ahead of print].
Zhou, X., Ma, Z., Zhang, L., Zhou, S., Wang, J., Wang, B., et al. (2016). Heart rate variability in the prediction of survival in patients with cancer: a systematic review and meta-analysis. J. Psychosom. Res. 89, 20–25. doi: 10.1016/j.jpsychores.2016.08.004
Keywords: sleep stages, heart rate, heart rate variability, reproducibility, autonomic nervous system, slow wave sleep
Citation: Herzig D, Eser P, Omlin X, Riener R, Wilhelm M and Achermann P (2018) Reproducibility of Heart Rate Variability Is Parameter and Sleep Stage Dependent. Front. Physiol. 8:1100. doi: 10.3389/fphys.2017.01100
Received: 30 August 2017; Accepted: 13 December 2017;
Published: 10 January 2018.
Edited by:Federico Lombardi, Università degli Studi di Milano, Italy
Reviewed by:Massimo Walter Rivolta, Università degli Studi di Milano, Italy
Heikki Veli Huikuri, University of Oulu, Finland
Antonio Luiz Pinho Ribeiro, Universidade Federal de Minas Gerais, Brazil
Copyright © 2018 Herzig, Eser, Omlin, Riener, Wilhelm and Achermann. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Prisca Eser, firstname.lastname@example.org