Front. SleepFrontiers in SleepFront. Sleep2813-2890Frontiers Media S.A.10.3389/frsle.2024.1230958SleepOriginal ResearchCardiopulmonary coupling-calculated sleep stability and nocturnal heart rate kinetics as a potential indicator for cardiovascular health: a relationship with blood pressure dippingHilmissonHugi1ThomasRobert Joseph2*MagnusdottirSolveig1*1MyCardio LLC, SleepImage®, Denver, CO, United States2Department of Medicine, Division of Pulmonary, Critical Care and Sleep Medicine, Beth Israel Deaconess Medical Center, Boston, MA, United States
Edited by: Ding Zou, University of Gothenburg, Sweden
Reviewed by: Sven Svedmyr, Sahlgrenska University Hospital, Sweden
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) and the copyright owner(s) 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.
Introduction
High blood pressure (HBP) is an independent, modifiable driver of cardiovascular (CV) morbidity and mortality. Nocturnal hypertension and non-dipping of blood pressure (NdBP) may be early markers of HBP. Similar to patients with NdBP, individuals with non-dipping of heart rate (NdHR) during sleep have an increased risk of CV disease, CV events, and CV-related mortality. The aim of this analysis was to evaluate if cardiopulmonary coupling (CPC) analysis-derived sleep states [stable/unstable non-rapid eye movement (NREM) sleep] and concomitant heart rate (HR) changes can provide information about nocturnal blood pressure (BP).
Method
Plethysmogram (pleth) signals from the HeartBEAT study (NCT01086800) were analyzed for CPC sleep states. Included in the analysis are sleep recordings from participants with acceptable pleth-signal quality at baseline (n = 302) and follow-up (n = 267), all having confirmed CV disease or CV-disease risk factors. The participants had a high prevalence of obstructive sleep apnea (OSA), 98.4% with moderate-OSA [apnea–hypopnea index (AHI) ≥ 15) and 29.6% severe OSA (AHI ≥ 30). A “heart-rate module” was created to evaluate the utility of identifying patients more likely to have BP dipping during sleep. Patients who did not have a decrease of ≥10% in their BP from wake to sleep were defined as NdBP and NdHR if their heart rate during stable-NREM sleep was higher than during unstable-NREM sleep.
Results
The most significant difference in minimum HR (HRmin) was observed when comparing BP dippers [56 ± 4 beats per minute (BPM)] and non-BP dippers (59 ± 4 BPM; p < 0.0001) during diastolic blood pressure in stable-NREM sleep. Higher HRmin were associated with an increased likelihood of being a non-dipper, with the strongest relationship with diastolic BP and stable-NREM sleep. Every increase of 1 BPM during stable-NREM sleep was associated with an ~4.4% increase in the probability of NdBP (p = 0.001). Subjects with NdHR have higher mean BP during sleep and wake periods than HR dippers. When continuous positive airway pressure therapy is efficacious, and a dipping pattern is achieved—physical and mental health is improved.
Conclusion
HR analytics in relation to the sleep period and the CPC spectrogram-estimated sleep states can provide novel and potentially clinically useful information on autonomic health. HR dipping (or not) may be a useful screener of BP dipping or non-dipping to identify individuals who may benefit from a formal assessment of 24-h ambulatory BP. Such a stepped approach may enable a more practical and applicable approach to diagnosing HBP.
Clinical Trial Registration
The Heart Biomarker Evaluation in Apnea Treatment (HeartBEAT) study is registered at clinicaltrials.gov/ct2/show/NCT01086800.
High blood pressure (HBP) is a major independent and modifiable driver of cardiovascular (CV) morbidity and mortality globally (Mills et al., 2020). Therefore, early diagnosis and optimal management of HBP are essential (Williams et al., 2018). Hypertension may be overt, with the physician's estimate and multiple home recordings in agreement, or it may be less obvious. The term masked hypertension (MH) was first introduced in 2002, describing a hypertension phenotype characterized by normal/not-increased office blood pressure (BP) readings and elevated out-of-office BP (Pickering et al., 2002). Identifying the MH phenotype is of interest for multiple reasons: (a) the prevalence is substantial, estimated at ~10%−30% of individuals attending hypertension clinics (O'Brien et al., 2013; Stergiou et al., 2021); (b) MH increases the risk of CV morbidity, with a meta-analysis comparing normotensive individuals and those with MH finding a 2.09 times increase in experiencing CV events for individuals with MH, while individuals with sustained hypertension have a 2.59 times increase in CV events (Pierdomenico and Cuccurullo, 2011); and (c) MH is more difficult to diagnose, as an out-of-office evaluation is required with 24-h ambulatory BP monitoring (24-ABPM). The availability of classic cuff 24-ABPM may be relatively limited and expensive, and the process is cumbersome and uncomfortable for the user, with cuff inflations and the frequency of measurements possibly causing arousals, disturbing sleep, and affecting the nighttime BP (NBP) measurement, which may not represent the true NBP (Pickering et al., 2006). While novel technologies are in development, including wrist pulse wave analyses, these are not yet fully validated. Moreover, such new technologies enable the collection of a large number of BP readings, which are not readily transferable to intermittent cuff inflation. Thus, disease-related validation remains to be done with more continuous forms of BP recording.
Nocturnal hypertension and non-dipping of BP may be early markers of HBP and MH. Nocturnal BP and night-to-day BP ratio have been identified as significant predictors of adverse CV outcomes (O'Brien et al., 1988; Asayama et al., 2023) and better predictors of fatal and non-fatal CV events and organ damage than daytime BP (Ohkubo et al., 2002; Salles et al., 2016; Staplin et al., 2023). Similar to changes in NBP, individuals with non-dipping of their heart rate (HR) during sleep have an increased risk of CV disease, CV events, and CV-related mortality (Eguchi et al., 2009; Kabutoya et al., 2010; Tadic et al., 2018; Nelde et al., 2023). For example, obstructive sleep apnea (OSA) is a common disease capable of disrupting normal BP control (Senaratna et al., 2017; Benjafield et al., 2019). A bidirectional relationship exists between OSA and HBP, with OSA patients having an increased risk of developing HBP, and the prevalence of OSA is higher in patients with HBP (Sawatari et al., 2016). However, other causes of sleep disruption such as insomnia and periodic limb movements are also associated with elevated NBP (Palagini et al., 2013).
Current guidelines highlight the importance of accurate diagnosis of HBP. The gold standard for identifying MH is 24-ABPM (Franklin et al., 2017), but still, the diagnosis is largely based on in-office BP measures, and patients with HBP and MH may therefore remain undiagnosed (Stergiou et al., 2021). Somewhat surprisingly, another readily available signal, HR during sleep [from ambulatory recordings and the millions of laboratory polysomnography (PSG) or home polygraphy done yearly], has not been subjected to rigorous analysis. Arousals from sleep are associated with both HR and BP elevations. It is thus plausible that an analysis of HR during sleep may provide a surrogate of BP or at least identify those who should have a formal 24-ABPM performed. HBP in patients with OSA is often predominantly nocturnal, and non-dipping BP is common (Marin et al., 2012). Therefore, alternative, simple, and less invasive measures that may estimate changes related to HR might be of value to identify non-dippers for further evaluation, including 24-ABPM as a tool that might assist in managing patients with OSA, insomnia, or comorbid insomnia and sleep apnea (Sweetman et al., 2019) and HBP (Tadic et al., 2018; Nelde et al., 2023).
This report targeted nocturnal heart rate (NHR) by analyzing photoplethysmogram (PPG) and oximetry information (SpO2) recorded during a conventional home sleep apnea test (HSAT). The cardiopulmonary coupling (CPC) sleep state and the HR kinetics analysis were computed to estimate alignments of HR with stable- and unstable-NREM sleep (Al Ashry et al., 2021b). It has been noticed that HR profiles across the night can on average drop (“HR dipping”), remain relatively flat (“HR non-dipping”), or even rise (“HR reverse dipping”). HR dipping is maximal during stable-NREM sleep periods [high-frequency coupling (HFC)] but can sometimes rise during such periods, suggesting increased sympathetic drive when normally NREM3/stable NREM2 should be associated with reduced sympathetic drive. We hypothesized that a simple method for calculating HR changes during stable-NREM sleep can be utilized to estimate dipping, non-dipping, or reverse-dipping HR. To test our hypothesis, we utilized the Heart Biomarker Evaluation in Apnea Treatment (HeartBEAT) study, which was designed to evaluate alternative approaches to reduce the risk of heart diseases. This study included classic cuff ABPM (Gottlieb et al., 2014).
2 Materials and methods2.1 Study design
The HeartBEAT study (NCT01086800) was a four-site, randomized, parallel-group trial among patients with high CV risk (Gottlieb et al., 2014). At baseline, patients were screened for OSA with the Berlin questionnaire (Netzer et al., 1999), and if they were at risk for OSA, an HSAT was initiated utilizing a portable sleep monitor, including recording PPG-signal and pulse-oximetry (SpO2) data from a fingertip sensor. Patients with an AHI of ≥15 events per hour of sleep were eligible to participate in the study. Patients with an AHI of >50 and a central index of >5 were excluded from randomization. The primary outcome of the HeartBEAT study was to evaluate changes in 24-h mean arterial blood pressure (MAP). In addition, nocturnal dipping and non-dipping of BP were reported, with non-dipping BP defined as a mean nocturnal BP higher than 90% of the mean daytime BP value. Participants in the HeartBEAT study were randomly assigned to one of three groups: continuous positive airway pressure (CPAP) therapy, nocturnal oxygen therapy, or healthy lifestyle with sleep education. Institutional review board approval was obtained from each participating institution. For this analysis, a data user agreement was obtained from the National Sleep Research Resource (Dean et al., 2016).
2.2 Participants
Patients aged 45–75 with established coronary heart disease or multiple CV-disease risk factors and well-managed hypertension, were recruited from cardiology practices at four participating medical centers. Patients with an AHI in the range of 15–50 were offered the opportunity to participate in the study. This study is based on a CPC analysis of the data derived from the fingertip PPG-sensor (the pleth waveform and oxygen data), as well as information reported about participants' dipping or non-dipping status.
2.3 Follow-up
A detailed description of the methodology and primary results of the trial's outcome have been reported (Gottlieb et al., 2014). In brief, outcomes were measured at baseline and 12 weeks after randomization. The primary outcome was 24-h MAP (measured using the 90207 Ambulatory Blood Pressure Monitor, Spacelabs Healthcare). The mean pressure was calculated at each reading as one-third of systolic pressure plus two-thirds of diastolic pressure. The 24-h mean pressure was calculated as a weighted average of the mean pressure during wakefulness and sleep, with the weights determined by the percentage of reported time spent in each state as recorded in a sleep diary. Nocturnal non-dipping blood pressure was defined as the mean nocturnal BP higher than 90% of the mean daytime value.
2.4 Methods2.4.1 The data set
The HeartBEAT study measured and reported BP dipping and non-dipping. Subjects were stratified into BP dippers and BP non-dippers, where non-dippers did not demonstrate a decrease in BP of ≥10% from wake to sleep.
2.4.2 CPC analysis
The CPC method has been described in detail in several prior publications (Thomas et al., 2005; Al Ashry et al., 2021a). Cardiopulmonary sleep spectrograms were first obtained from a single lead electrocardiogram (ECG). ECG-derived respiration (EDR) is obtained either by using R-S wave amplitudes or variations in the QRS complex area. Ectopic beats are identified and removed, normal sinus–normal sinus (NN) intervals are extracted, and outliers are filtered (Thomas et al., 2005). After extracting the NN interval series on ECG and its associated EDR, the signals are then resampled using cubic splines at 2 Hz. The fast Fourier transform is applied to three overlapping 512-sample sub-windows within the 1,024-coherence window. The 1,024-coherence window is then advanced by 256 samples (2.1 min), and the calculation is repeated until the entire NN interval/EDR series is analyzed. Thus, the cross-spectral power and coherence of these two signals are calculated over a 1,024-sample (8.5-min) window. For each 1,024-sample window, the product of the coherence and cross-spectral power is used to calculate the ratio of coherent cross-power in the low-frequency (0.01–0.1 Hz) band to that in the high-frequency (0.1–0.4 Hz) band. The logarithm of the high-to-low-frequency CPC ratio [log (HFC/LFC)] is then computed to yield a continuously varying measure of CPC sleep stability/instability output metrics. While, originally, the ECG signal was used as input, any signal or signal set that encodes respiration and heart rate variability (HRV) may be used to compute the CPC sleep spectrogram; most conveniently, this signal set can be obtained from the peripheral PPG-signal, which is readily available from current generation oximeters. The current embodiment uses a ring- or fingertip-based oximeter to collect the data coupled with a mobile application and Bluetooth to stream the data for automated analysis. The SleepImage system complies with the Health Insurance Portability and Accountability Act (HIPAA), is cleared by the Food and Drug Administration (K182618), and complies with the EU Medical Device Directive (CE-mark 2862) to automatically generate biomarkers, presented numerically and graphically (Figure 1). The analysis is otherwise essentially identical. The outputs of the CPC algorithm include low-frequency (LFC; 0.01–0.1 Hz), high-frequency (HFC; 0.1–0.4 Hz), and very low-frequency (vLFC; 0.001–0.01 Hz) couplings, and an elevated LFC-broad band (eLFCBB) that is a sleep fragmentation signal biomarker (Thomas et al., 2005). HFC/LFC covary more strongly with an electroencephalographic (EEG) non-cyclic alternating pattern (n-CAP) and CAP, respectively (Thomas et al., 2005), than conventional N3/N2—although most of N3 is HFC, and much of HFC occurs during N2.
Heart rate kinetics during sleep. The arrows indicate heart rate characteristics (green: dipping, black: non-dipping, red: reverse dipping).
Stable-NREM sleep (HFC) is associated with several desirable sleep characteristics, including increased absolute and relative delta power (Thomas et al., 2014), a consolidated NREM sleep <1-Hz slow oscillation, temporally stable breathing, stable arousal thresholds, normal arterial oxygen (O2) and carbon dioxide (CO2) concentrations, and BP dipping (Wood et al., 2020). Unstable NREM (LFC) is characterized by features opposite of stable-NREM (HFC), and ineffective (fragmented) REM sleep takes on LFC coupling signatures, while wake or effective REM sleep shows vLFC pattern (Thomas et al., 2005). HFC covaries better with relative than absolute EEG slow-wave power and is thus less constrained by the “loss” of slow-wave sleep (SWS) with age (Thomas et al., 2014). Specific spectrographic signatures of fragmented sleep (elevated LFC narrow band, eLFCNB) are biomarkers of strong chemoreflex effects on sleep respiration (Thomas et al., 2007), identifying areas of sleep with central apneas and periodic breathing. BP dipping occurs only during periods of HFC (Wood et al., 2020), consistent with the demonstration that non-CAP is the EEG correlate of BP dipping (Iellamo et al., 2004). LFC is associated with hypertension and stroke (Thomas et al., 2009), while HFC is reduced in depression (Yang et al., 2011), heart failure (Yeh et al., 2008), and fibromyalgia syndrome (Thomas et al., 2010). HFC is an independent determinant of the glucose disposition index (Pogach et al., 2012). Pre- and posttreatment effects in sleep apnea are captured via changes in HFC/LFC (Lee et al., 2014). An integrated metric, the Sleep Quality Index (SQI), which is heavily weighted by stable-NREM sleep, is associated with desirable directions of metabolic health and blood pressure (Magnusdottir et al., 2020, 2021, 2022).
2.4.3 HR and CPC analysis
The PPG signal from each polygraphy recording in the HeartBEAT data was processed through the SleepImage System algorithms. A software module (“HR module”) was developed for evaluating the HR data collected during the sleep study, which was then intersected with the sleep-state output from the SleepImage System: (a) A 2-sample-per-second resampled NN series was evaluated by cropping the entire series from sleep onset to sleep conclusion. Then, a linear trendline was fitted to generate a slope coefficient, the associated p-value, and an R2 for model fit. While this method may provide insight that some clinical professionals may value, it was not expected to have much explanatory power, as it attempts to describe HR as a linear trend over the course of the sleep period, lumping all sleep states together. The disproportionate effect of sleep stage on BP dynamics has been explained in the literature (Stein and Pu, 2012). (b) HR rate statistics for each CPC sleep state classification were then calculated using the same 2-sample-per-second resampled NN series.
2.4.4 Primary endpoint
For this analysis, the following variables of interest were extracted from the HeartBEAT data set: (a) MAP non-dipping and dipping, (b) systolic blood pressure (SBP) non-dipping and dipping, and (c) diastolic blood pressure (DBP) non-dipping and dipping, where non-dippers were defined as patients who did not demonstrate a decrease in BP ≥10% from wake to sleep. The primary endpoint was to evaluate (a) the relationship between HR during stable and unstable CPC sleep states and BP and (b) if sleep evaluation combined with this HR module can be utilized to stratify patients to identify patients more likely to have non-BP dipping during sleep for further evaluation.
2.5 Statistical analysis
The primary endpoint, questioning if the HR module can be utilized to identify BP dipping during sleep, was evaluated using a logistic regression, regressing HR parameters from the HR module on the dipping and non-dipping variables extracted from the HeartBEAT data, controlling for gender, age, and race to assess the predictive power on BP dipping.
Basic summary statistics, such as counts, are presented, along with means and standard errors (in parentheses) for gender, race, age, and body mass index (BMI). A one-way analysis of covariance was utilized to compare non-dippers and dippers controlling for age, gender, BMI, and site identifier. Means and standard errors for the HR parameters are presented, along with the p-values for pairwise comparison between the groups. The analysis was performed using Stata version 15.1.
3 Results3.1 Study population
Included in this analysis of the pleth signal from the HeartBEAT data (Gottlieb et al., 2014) are both baseline and follow-up recordings with signal quality defined as average successful peak detection on the overnight recording of no <70% for analysis to evaluate HR during sleep and on data describing dipping and non-dipping BP status. The data set contained sleep recordings with acceptable pleth-signal quality from 302 patients at baseline and 267 at follow-up, or a total of 569 sleep recordings (Supplementary Table S1).
3.2 Primary outcome measures
First, the relationship between the CPC-sleep-state analysis and the software HR module and dipping or non-dipping status of BP was evaluated (Table 1). The most significant difference in minimum HR (HRmin) was observed when comparing HR for dippers (56 ± 4 BPM) and non-dippers (59 ± 4 BPM) during DBP and stable-NREM sleep (p < 0.0001).
Heart rate (HR) software module (“HR parameters”) for dippers and non-dippers, beats per minutes (BPM).
Average HR (BPM) measured during diastolic blood pressure (mm/Hg)
Average HR (BPM) measured during systolic blood pressure (mm/Hg)
Average HR (BPM) measured during mean arterial blood pressure (mm/Hg)
Non-Dipa
Dipa
p-valueb
Non-Dipa
Dipa
p-valueb
Non-Dipa
Dipa
p-valueb
Sleep
Mean HR
66.4 (±5.5)
64.7 (±5.5)
0.063
66.0 (±5.5)
65.1 (±5.5)
0.301
66.6 (±5.5)
65.1 (±5.5)
0.090
Min HR
55.6 (±2.6)
54.31 (±2.6)
0.003
55.6 (±2.6)
54.6 (±2.6)
0.016
55.8 (±2.6)
54.6 (±2.6)
0.003
Wake
Mean HR
69.2 (±5.7)
67.8 (±5.6)
0.126
69.1 (±5.7)
68.1 (±5.6)
0.293
69.6 (±5.7)
68.1 (±5.6)
0.100
Min HR
57.5 (±3.0)
55.9 (±3.0)
0.001
57.7 (±3.0)
56.2 (±3.0)
0.003
57.871 (±3.0)
56.3 (±3.0)
0.001
Stable NREM
Mean HR
67.2 (±7.4)
64.3 (7.3)
0.017
66.1 (±7.4)
64.9 (±7.34)
0.318
66.6 (±7.4)
64.6 (±7.3)
0.150
Min HR
58.6 (±4.1)
56.3 (4.1)
<0.0001
57.6 (±4.2)
56.74 (±4.1)
0.189
58.25 (±4.2)
56.7 (±4.1)
0.022
HR, heart rate; NREM, non-rapid eye movement sleep; BPM, heart beats per minute.
aValues are mean ± standard deviation.
bDifference between groups (p-value).
Table 2 presents the result from different logit models regressing variables of interest (HR metrics; HRmin during sleep and HRmin during stable-NREM sleep), controlling for age and race on the indicator variables for diastolic, systolic, and MAP non-dipping. BMI was explored as a control but did not add to the explanatory/predictive power of any of the models. The results are presented as coefficient estimates from the logit regression with p-values in parentheses. The strongest associations were observed when including HRmin during stable-NREM sleep. Higher HRmin were associated with a higher likelihood of being a non-dipper, with an increase of 1 beat per minute in HRmin during stable NREM being associated with an approximate 4.4% increase in the probability of being a non-dipper (p = 0.001). Being African American and increased age were found to be factors that increased the likelihood of being a non-dipper, while gender did not seem to be a significant predictor. For the sake of brevity, the mean HR in stable-NREM sleep was not included in the results, although it showed a statistically significant effect on diastolic non-dipping, in favor of the minimum, which had a higher coefficient. The variance inflation factor (VIF) was calculated to investigate the presence of collinearity. The average VIF did not exceed 1.08 for any of the models, with the highest calculated VIF of 1.14 for an individual variable. This indicates that multicollinearity is not of concern with the chosen independent variables.
Logistic regression models regressing the HR parameter on diastolic, systolic, and mean arterial blood pressure (MAP) non-dip indicators.
Diastolic non-dipping blood pressurea
Systolic non-dipping blood pressurea
MAP non-dipping blood pressurea
Sleepmin
0.049 (0.010)
0.041 (0.035)
0.049 (0.010)
Stable NREMmin
0.044 (0.001)
0.026 (0.032)
Male
0.092 (0.674)
0.117 (0.591)
0.030 (0.891)
0.153 (0.483)
0.193 (0.378)
Race
African American
0.744 (0.009)
0.720 (0.012)
1.013 (0.001)
0.960 (0.001)
0.921 (0.002)
Other
0.086 (0.815)
−0.002 (0.995)
0.007 (0.983)
0.195 (0.592)
0.117 (0.749)
Age
0.019 (0.146)
0.023 (0.073)
0.035 (0.008)
0.026 (0.049)
0.032 (0.016)
Const.
−4.205 (0.000)
−4.572 (0.001)
−4.400 (0.002)
−3.532 (0.002)
−5.010 (0.000)
NREM, non-rapid eye movement sleep; MAP, mean arterial pressure.
aValues are the logistic regression coefficients with p-values in parenthesis.
Second, the cohort was stratified based on the “HR parameter” into HR dippers and HR non-dippers (Table 3). At baseline, fewer HR non-dippers used calcium-channel blockers (10.7%, p = 0.047) and diuretics (10.9%, p = 0.044) than HR-dippers. Additionally, HR non-dippers presented with a higher mean MAP when awake of 2.2 mmHg (p = 0.046) than HR-dippers. Comparing HR dippers and non-dippers and focusing on BP, HR non-dippers have significantly higher average DBP (DBPave; 85.3 vs. 83.8 mmHg, p = 0.035); DBPsleep (74.6 vs. 72.7 mmHg, p = 0.018); and DBPwake (90.8 vs. 89.0 mmHg, p = 0.024). Additionally, they have a higher average MAP (MAPave; 89.7 vs. 88.1 mmHg, p = 0.050); MAPsleep (79.4 vs. 77.4 mmHg, p = 0.036); and MAPwake (95.2 vs. 93.3 mmHg, p = 0.031), respectively (Table 4).
Baseline characteristics: heart rate dippers vs. non-dippers.
aValues are mean ± standard deviation for continuous variables and percentages for categorical variables.
bDifference between groups (p-value).
Finally, the subgroup that received the CPAP therapy was stratified based on dipping status (HR dipper or HR non-dipper) at baseline and at a 12-week follow-up (Table 5). When comparing participants who were HR non-dippers at baseline and HR dippers at follow-up, significant improvements were observed in depression severity measured using the Patient Health Questionnaire-9, −3.1 (p = 0.005), and the health and quality-of-life indicators that were evaluated using the Short Form (36) Health Survey: (a) vitality, 12.8 (p = 0.036); (b) physical functioning, 20.1 (p = 0.003); and (c) emotional functioning, 16.5 (p = 0.005); social functioning, 17.4 (p = 0.007); and mental health and emotional wellbeing, 13.9 (p = 0.001).
Heart rate dipper vs. non-dippers: comparison of characteristics.
aValues are mean ± standard deviation for continuous variables and percentages for categorical variables.
bP-values based on difference between groups at baseline and follow-up.
4 Discussion
The analysis of HR across the night in relation to CPC sleep state, HR, and BP showed the following statistically significant results: (a) lower HRmin during sleep in participants who demonstrated diastolic, systolic, and MAP dipping when compared to non-dippers; (b) lower HRmin during wake within the sleep period in participants who demonstrated diastolic-, systolic-, and MAP dipping compared to non-dippers; (c) lower HRmin and HRmean during stable-NREM sleep (HFC) in those who demonstrated diastolic dipping compared to non-dippers; (d) the strongest associations were observed when including HRmin during stable-NREM sleep (HFC); (e) when utilizing the HR module to stratify the cohort based on HR dipping and non-dipping, participants with non-dipping of HR had significantly higher MAP and DBP when comparing all readings, readings during sleep, and readings during wake; and (f) in the subgroup of participants who received CPAP-therapy, participants with a HR non-dipping pattern at baseline and HR dipping pattern at follow-up (i.e., treatment responders based on heart rate profiles), significantly improved their subjective mental and physical functioning. The results overall suggest that HR dipping in stable-NREM sleep/HFC is a desirable biological characteristic.
HR kinetics during sleep seems to provide indirect information about BP during sleep, an important cardiovascular health variable. HR is readily available through most systems that track sleep oximetry and could allow for risk stratification; individuals with a non-dipping HR pattern could be directed to selectively undergo conventional ABPM. The finding that HR was specifically influenced by stable-NREM sleep (HFC) as estimated through the pleth spectrogram was not surprising. Autonomic physiology presents an important window into sleep; for example, hemodynamics, HRV, and respiration are markedly dependent on the sleep stage, with vagal dominance, stable breathing, and BP reductions (“dipping”) during SWS/N3 (Javaheri and Redline, 2012). Furthermore, standard reporting of EEG-based stages as a percentage of sleep time is an insensitive metric of sleep fragmentation (Bianchi and Thomas, 2013). Most HSATs do not provide EEG stages, although machine learning applied to ECG and respiratory signals can approximate deep sleep. The ability to evaluate HR dynamics in relation to sleep state without the need for extensive PSG has practical advantages.
A reduction in BP during sleep (BP dipping) is considered a BP-related biomarker of healthy sleep (Routledge et al., 2007; Salles et al., 2016). There is a progressive reduction of BP from wake through SWS, with an increase in REM sleep or transiently in association with arousals. The HR profile follows this basic scheme and is the highest in REM sleep and unstable-NREM sleep or during arousals and lowest in conventional N3 (Javaheri and Redline, 2012). In a prior study by our group, we used PSG with beat-to-beat BP monitoring, ECG-derived CPC analysis, and quantified delta power and the rate of occurrence of the <1-Hz slow oscillation. We found that BP dipping occurred only during periods of stable NREM (HFC), concomitant with slow oscillation/delta power-enriched NREM sleep. HR was lowest in N3, but the small sample size of 11 subjects perhaps explained the lack of HR dipping during HFC; however, the current analysis shows the predicted dip in HR. Mechanisms associated with rising slow wave/delta power and a high grade of electrocortical synchrony are likely the drivers of an integrated response of BP, HR, and stable breathing. Even in conditions associated with fragmented sleep, such as sleep apnea, delta power, and vagal HRV dominance tends to ebb and flow in a correlated manner (Jurysta et al., 2006; Wood et al., 2020).
There is substantial variability in sleep quality in individuals with similar severities of sleep apnea. Such differences may be quantified by subjective symptoms e.g., questionnaire such as the Insomnia Severity Index (Bastien et al., 2001); conventional criteria e.g., N1, N3, total sleep time, sleep efficiency; EEG-based methods such as the odds ratio product, which estimates sleep depth continuously (Younes et al., 2015; Younes, 2023); and ECG/PPG CPC spectrograms and the SQI (Thomas et al., 2014; Hilmisson et al., 2019; Magnusdottir et al., 2020). BP during sleep is another useful measure, while HR analysis could provide a complementary metric for sleep quality.
There are sleep and circadian influences on BP and HR control. Even in forced desynchrony experiments, both metrics are low in the biological night when body temperature is low and melatonin is high. Thus, there are both sleep and circadian components to BP and HR dipping, and a loss of this pattern can occur from either sleep or circadian factors. Any case of sleep fragmentation can flatten or even reverse BP and HR during the biological night. OSA can affect sleep quality and cause non-dipping BP by autonomic dysfunction, transient surges associated with arousals, the upregulation of neurohumoral systems, oxidative stress, and a general lowering of sleep depth. However, when there is non-dipping of either BP or HR during stable (unfragmented) sleep of which HFC is a good biomarker, it likely reflects abnormal autonomic regulation as many other drivers are less active during this state (e.g., breathing and oxygenation are stable). Profiles of HR may be useful when following treatment of sleep apnea, especially with therapies with residual apnea. While CPAP when used can largely normalize breathing, other increasingly used therapies, such as weight loss, hypoglossal nerve stimulation, and oral appliances, are more likely to have residual apnea. Partial CPAP use will also demonstrate residual apnea. HR profiles could be one way to assess the impact of residual disease as successful therapy could be expected to improve HR dipping and even convert a non-/reverse dipper to a dipper as demonstrated in this analysis when looking at the subgroup treated with CPAP.
This analysis has some limitations, including (a) the study population was selected for presence of CV disease or risk and does not readily generalize to the range of medical backgrounds on which sleep apnea occurs; (b) classic sleep staging was not available or possible; (c) only one night of recording was available at any given time point; (d) HR can be modified by numerous factors including stress, anxiety, and pain, the impact of which on the type of analysis we performed is unknown; (e) the implications of HR-pattern analysis for disease prognostics, well established for conventional ambulatory blood pressure, are unknown; (f) the impact of drugs such as beta-blockers or antihypertensives in general on the noted patterns need to be established, and it is plausible that both attenuation or amplification of the HR response during sleep may occur based on cardiovascular functional status; and (g) conditions such as heart failure, advanced renal disease, post-cardiac transplant, or advanced autonomic neuropathy are likely to have relatively unchanging HR across the night and may not allow this analysis. Heart rate analysis would be invalidated by atrial fibrillation and during fixed-rate cardiac pacing for bradyarrhythmias, while other modes of pacing may still allow analysis, but that needs to be directly demonstrated. The parent study design and our current analysis cannot determine if HR is an independent risk factor (beyond nocturnal BP) for cardiovascular outcomes. Additionally, a generalization of its potential usefulness to non-apnea conditions such as insomnia or restless legs and their treatments cannot be determined.
In conclusion, HR analytics in relation to the sleep period and the CPC spectrogram-estimated sleep states can provide novel and potentially clinically useful information on autonomic health. HR dipping (or not) may be a useful screener of BP dipping/non-dipping and identify individuals who may benefit from formal assessment of ambulatory BP and/or evaluate the efficacy of various therapies. Such a stepped approach may enable a more practical, cost-effective, and applicable approach to diagnosing MH.
Data availability statement
Publicly available datasets were analyzed in this study. This analysis was made possible by the National Sleep Research Resource, with access to the HeartBEAT-database—https://www.sleepdata.org/datasets/heartbeat.
Ethics statement
The studies involving humans were approved by Brigham & Women's Hospital and Case Western Reserve University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
HH conceptualized and designed the analysis, analyzed the data and conducted the statistical analysis, interpreted the data output, and reviewed, edited, and approved the overall content of the final manuscript. RT and SM wrote the manuscript, interpreted the data, and reviewed, edited, and approved the overall content of the final manuscript. All authors approved the final version of the manuscript and agreed to be accountable for all aspects of the study in ensuring that questions related to the accuracy or integrity of any part of the study are appropriately investigated and resolved.
Funding
The HeartBEAT study was supported by grants from the National Heart, Lung, and Blood Institute (RC2HL101417, 1R01HL109493, and R21HL108226) and by a grant from the National Center for Research Resources (UL1 RR024989).
Conflict of interest
HH is Director of Research & Development for MyCardio LLC. SleepImage is the brand name of MyCardio LLC, a privately held entity. RT reports patent, license, and royalties from MyCardio LLC; patent, license, and royalties from DeVilbiss-Drive; and general sleep medicine consulting through Guidepoint Global and GLG Councils. SM is Chief Medical Officer for MyCardio LLC and has a partial ownership. SleepImage is the brand name of MyCardio LLC, a privately held entity. MyCardio LLC is a licensee of the Cardiopulmonary coupling (CPC) algorithms, a method to phenotype sleep and sleep apnea, from the Beth Israel Deaconess Medical Center, Boston, MA, USA.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frsle.2024.1230958/full#supplementary-material
ReferencesAl AshryH. S.HilmissonH.NiY.ThomasR. J.Apples Investigators. (2021a). Automated apnea-hypopnea index from oximetry and spectral analysis of cardiopulmonary coupling. Ann. Am. Thorac. Soc.18, 876–883. 10.1513/AnnalsATS.202005-510OC33472017Al AshryH. S.NiY.ThomasR. J. (2021b). Cardiopulmonary sleep spectrograms open a novel window into sleep biology—implications for health and disease. Front. Neurosci.15:755464. 10.3389/fnins.2021.75546434867165AsayamaK.Stolarz-SkrzypekK.YangW. Y.HansenT. W.Brguljan-HitijJ.OdiliA. N.. (2023). What did we learn from the international databases on ambulatory and home blood pressure in relation to cardiovascular outcome?Hypertens. Res.46, 934–949. 10.1038/s41440-023-01191-436737461BastienC. H.VallieresA.MorinC. M. (2001). Validation of the Insomnia Severity Index as an outcome measure for insomnia research. Sleep Med.2, 297–307. 10.1016/S1389-9457(00)00065-411438246BenjafieldA. V.AyasN. T.EastwoodP. R.HeinzerR.IpM. S. M.MorrellM. J.. (2019). Estimation of the global prevalence and burden of obstructive sleep apnoea: a literature-based analysis. Lancet Respir. Med. 7, 687–698. 10.1016/S2213-2600(19)30198-531300334BianchiM. T.ThomasR. J. (2013). Technical advances in the characterization of the complexity of sleep and sleep disorders. Prog. Neuropsychopharmacol. Biol. Psychiatry45, 277–286. 10.1016/j.pnpbp.2012.09.01723174482DeanD. A.2ndGoldbergerA. L.MuellerR.KimM.RueschmanM.MobleyD.. (2016). Scaling up scientific discovery in sleep medicine: the national sleep research resource. Sleep39, 1151–1164. 10.5665/sleep.577427070134EguchiK.HoshideS.IshikawaJ.PickeringT. G.SchwartzJ. E.ShimadaK.. (2009). Nocturnal nondipping of heart rate predicts cardiovascular events in hypertensive patients. J. Hypertens.27, 2265–2670. 10.1097/HJH.0b013e328330a93819834343FranklinS. S.O'BrienE.StaessenJ. A. (2017). Masked hypertension: understanding its complexity. Eur. Heart J.38, 1112–1118. 10.1093/eurheartj/ehw50227836914GottliebD. J.PunjabiN. M.MehraR.PatelS. R.QuanS. F.BabineauD. C.. (2014). CPAP versus oxygen in obstructive sleep apnea. N. Engl. J. Med.370, 2276–2285. 10.1056/NEJMoa130676624918372HilmissonH.LangeN.MagnusdottirS. (2019). Objective sleep quality and metabolic risk in healthy weight children results from the randomized Childhood Adenotonsillectomy Trial (CHAT). Sleep Breath.23, 1197–1208. 10.1007/s11325-019-01802-w30798410IellamoF.PlacidiF.MarcianiM. G.RomigiA.TombiniM.AquilaniS.. (2004). Baroreflex buffering of sympathetic activation during sleep: evidence from autonomic assessment of sleep macroarchitecture and microarchitecture. Hypertension43, 814–819. 10.1161/01.HYP.0000121364.74439.6a14981054JavaheriS.RedlineS. (2012). Sleep, slow-wave sleep, and blood pressure. Curr. Hypertens. Rep.14, 442–448. 10.1007/s11906-012-0289-022846982JurystaF.LanquartJ. P.van de BorneP.MigeotteP. F.DumontM.DegauteJ. P.. (2006). The link between cardiac autonomic activity and sleep delta power is altered in men with sleep apnea-hypopnea syndrome. Am. J. Physiol. Regul. Integr. Comp. Physiol.291, R1165–R1171. 10.1152/ajpregu.00787.200516675631KabutoyaT.HoshideS.IshikawaJ.EguchiK.ShimadaK.KarioK.. (2010). The effect of pulse rate and blood pressure dipping status on the risk of stroke and cardiovascular disease in Japanese hypertensive patients. Am. J. Hypertens.23, 749–755. 10.1038/ajh.2010.4520300073LeeW. H.AhnJ. C.WeJ.RheeC. S.LeeC. H.YunP. Y.. (2014). Cardiopulmonary coupling analysis: changes before and after treatment with a mandibular advancement device. Sleep Breath. 18, 891–896. 10.1007/s11325-014-0961-524699839MagnusdottirS.HilmissonH.ThomasR. J. (2020). Cardiopulmonary coupling-derived sleep quality is associated with improvements in blood pressure in patients with obstructive sleep apnea at high-cardiovascular risk. J. Hypertens.38, 2287–2294. 10.1097/HJH.000000000000255332649638MagnusdottirS.ThomasR. J.HilmissonH. (2021). Can improvements in sleep quality positively affect serum adiponectin-levels in patients with obstructive sleep apnea?Sleep Med.84, 324–333. 10.1016/j.sleep.2021.05.03234225174MagnusdottirS.WitmansM.HilmissonH. (2022). Sleep quality, sleep apnea, and metabolic health in children treated with adenotonsillectomy. Sleep Breath.27, 1527–1540. 10.1007/s11325-022-02747-336434376MarinJ. M.AgustiA.VillarI.FornerM.NietoD.CarrizoS. J.. (2012). Association between treated and untreated obstructive sleep apnea and risk of hypertension. JAMA307, 2169–2176. 10.1001/jama.2012.341822618924MillsK. T.StefanescuA.HeJ. (2020). The global epidemiology of hypertension. Nat. Rev. Nephrol.16, 223–237. 10.1038/s41581-019-0244-232024986NeldeA.KlammerM. G.NolteC. H.StenglH.KramerM.von RennenbergR.. (2023). Data lake-driven analytics identify nocturnal non-dipping of heart rate as predictor of unfavorable stroke outcome at discharge. J Neurol.270, 3810–3820. 10.1007/s00415-023-11718-x37079032NetzerN. C.StoohsR. A.NetzerC. M.ClarkK.StrohlK. P. (1999). Using the Berlin Questionnaire to identify patients at risk for the sleep apnea syndrome. Ann. Intern. Med.131, 485–491. 10.7326/0003-4819-131-7-199910050-0000210507956O'BrienE.ParatiG.StergiouG.AsmarR.BeilinL.BiloG.. (2013). European Society of Hypertension position paper on ambulatory blood pressure monitoring. J. Hypertens.31, 1731–1768. 10.1097/HJH.0b013e328363e96424029863O'BrienE.SheridanJ.O'MalleyK. (1988). Dippers and non-dippers. Lancet2:397. 10.1016/S0140-6736(88)92867-X2899801OhkuboT.HozawaA.YamaguchiJ.KikuyaM.OhmoriK.MichimataM.. (2002). Prognostic significance of the nocturnal decline in blood pressure in individuals with and without high 24-h blood pressure: the Ohasama study. J. Hypertens.20, 2183–2189. 10.1097/00004872-200211000-0001712409956PalaginiL.BrunoR. M.GemignaniA.BaglioniC.GhiadoniL.RiemannD.. (2013). Sleep loss and hypertension: a systematic review. Curr. Pharm. Des.19, 2409–2419. 10.2174/138161281131913000923173590PickeringT. G.DavidsonK.GerinW.SchwartzJ. E. (2002). Masked hypertension. Hypertension40, 795–796. 10.1161/01.HYP.0000038733.08436.9812468559PickeringT. G.ShimboD.HaasD. (2006). Ambulatory blood-pressure monitoring. N. Engl. J. Med.354, 2368–2324. 10.1056/NEJMra06043316738273PierdomenicoS. D.CuccurulloF. (2011). Prognostic value of white-coat and masked hypertension diagnosed by ambulatory monitoring in initially untreated subjects: an updated meta analysis. Am. J. Hypertens.24, 52–58. 10.1038/ajh.2010.20320847724PogachM. S.PunjabiN. M.ThomasN.ThomasR. J. (2012). Electrocardiogram-based sleep spectrogram measures of sleep stability and glucose disposal in sleep disordered breathing. Sleep35, 139–148. 10.5665/sleep.160422215928RoutledgeF. S.McFetridge-DurdleJ. A.DeanC. R.Society Canadian Hypertension. (2007). Night-time blood pressure patterns and target organ damage: a review. Can. J. Cardiol.23, 132–138. 10.1016/S0828-282X(07)70733-X17311119SallesG. F.ReboldiG.FagardR. H.CardosoC. R.PierdomenicoS. D.VerdecchiaP.. (2016). prognostic effect of the nocturnal blood pressure fall in hypertensive patients: the ambulatory blood pressure collaboration in patients with hypertension (ABC-H) meta-analysis. Hypertension67, 693–700. 10.1161/HYPERTENSIONAHA.115.0698126902495SawatariH.ChishakiA.AndoS. I. (2016). The epidemiology of sleep disordered breathing and hypertension in various populations. Curr. Hypertens. Rev.12, 12–17. 10.2174/157340211266616011409330726778203SenaratnaC. V.PerretJ. L.LodgeC. J.LoweA. J.CampbellB. E.MathesonM. C.. (2017). Prevalence of obstructive sleep apnea in the general population: a systematic review. Sleep Med. Rev.34, 70–81. 10.1016/j.smrv.2016.07.00227568340StaplinN.de la SierraA.RuilopeL. M.EmbersonJ. R.VinyolesE.GorostidiM.. (2023). Relationship between clinic and ambulatory blood pressure and mortality: an observational cohort study in 59 124 patients. Lancet401, 2041–2050. 10.1016/S0140-6736(23)00733-X37156250SteinP. K.PuY. (2012). Heart rate variability, sleep and sleep disorders. Sleep Med. Rev.16, 47–66. 10.1016/j.smrv.2011.02.00521658979StergiouG. S.PalatiniP.ParatiG.O'BrienE.JanuszewiczA.LurbeE.. (2021). 2021 European Society of Hypertension practice guidelines for office and out-of-office blood pressure measurement. J. Hypertens.39, 1293–1302. 10.1097/HJH.000000000000284333710173SweetmanA.LackL.BastienC. (2019). Co-morbid insomnia and sleep apnea (COMISA): prevalence, consequences, methodological considerations, and recent randomized controlled trials. Brain Sci.9:371. 10.3390/brainsci912037131842520TadicM.CuspidiC.GrassiG. (2018). Heart rate as a predictor of cardiovascular risk. Eur. J. Clin. Invest. 48. 10.1111/eci.1289229355923ThomasR. J.MietusJ. E.PengC. K.GilmartinG.DalyR. W.GoldbergerA. L.. (2007). Differentiating obstructive from central and complex sleep apnea using an automated electrocardiogram-based method. Sleep30, 1756–1769. 10.1093/sleep/30.12.175618246985ThomasR. J.MietusJ. E.PengC. K.GoldbergerA. L. (2005). An electrocardiogram-based technique to assess cardiopulmonary coupling during sleep. Sleep28, 1151–1161. 10.1093/sleep/28.9.115116268385ThomasR. J.MietusJ. E.PengC. K.GoldbergerA. L.CroffordL. J.ChervinR. D.. (2010). Impaired sleep quality in fibromyalgia: detection and quantification with ECG-based cardiopulmonary coupling spectrograms. Sleep Med. 11, 497–498. 10.1016/j.sleep.2009.09.00320015685ThomasR. J.MietusJ. E.PengC. K.GuoD.GozalD.Montgomery-DownsH.. (2014). Relationship between delta power and the electrocardiogram-derived cardiopulmonary spectrogram: possible implications for assessing the effectiveness of sleep. Sleep Med.15, 125–131. 10.1016/j.sleep.2013.10.00224269134ThomasR. J.WeissM. D.MietusJ. E.PengC. K.GoldbergerA. L.GottliebD. J.. (2009). Prevalent hypertension and stroke in the Sleep Heart Health Study: association with an ECG-derived spectrographic marker of cardiopulmonary coupling. Sleep32, 897–904.19639752WilliamsB.ManciaG.SpieringW.Agabiti RoseiE.AziziM.BurnierM.. (2018). 2018 ESC/ESH Guidelines for the management of arterial hypertension. Eur. Heart J.39, 3021–3104. 10.1093/eurheartj/ehy33930165516WoodC.BianchiM. T.YunC. H.ShinC.ThomasR. J. (2020). Multicomponent analysis of sleep using electrocortical, respiratory, autonomic and hemodynamic signals reveals distinct features of stable and unstable NREM and REM sleep. Front. Physiol.11:592978. 10.3389/fphys.2020.59297833343390YangA. C.YangC. H.HongC. J.TsaiS. J.KuoC. H.PengC. K.. (2011). Sleep state instabilities in major depressive disorder: detection and quantification with electrocardiogram-based cardiopulmonary coupling analysis. Psychophysiology48, 285–291. 10.1111/j.1469-8986.2010.01060.x20624250YehG. Y.MietusJ. E.PengC. K.PhillipsR. S.DavisR. B.WayneP. M.. (2008). Enhancement of sleep stability with Tai Chi exercise in chronic heart failure: preliminary findings using an ECG-based spectrogram method. Sleep Med.9, 527–536. 10.1016/j.sleep.2007.06.00317689142YounesM. (2023). New insights and potential clinical implications of the odds ratio product. Front. Neurol.14:1273623. 10.3389/fneur.2023.127362337885480YounesM.OstrowskiM.SoifermanM.YounesH.YounesM.RaneriJ.. (2015). Odds ratio product of sleep EEG as a continuous measure of sleep state. Sleep38, 641–654. 10.5665/sleep.458825348125