Wenquan Cheng1†
Xingran Cui- 1State Key Laboratory of Bioelectronics, School of Biological Science and Medical Engineering, Southeast University, Nanjing, China
- 2Endocrinology Department, Suzhou Science and Technology Town Hospital, Suzhou, China
- 3Institute of Medical Devices (Suzhou), Southeast University, Suzhou, China
Introduction: Autonomic nervous system (ANS) plays an important role in the exchange of metabolic information between organs and regulation on peripheral metabolism with obvious circadian rhythm in a healthy state. Sleep, a vital brain phenomenon, significantly affects both ANS and metabolic function.
Objectives: This study investigated the relationships among sleep, ANS and metabolic function in type 2 diabetes mellitus (T2DM), to support the evaluation of ANS function through heart rate variability (HRV) metrics, and the determination of the correlated underlying autonomic pathways, and help optimize the early prevention, post-diagnosis and management of T2DM and its complications.
Materials and methods: A total of 64 volunteered inpatients with T2DM took part in this study. 24-h electrocardiogram (ECG), clinical indicators of metabolic function, sleep quality and sleep staging results of T2DM patients were monitored.
Results: The associations between sleep quality, 24-h/awake/sleep/sleep staging HRV and clinical indicators of metabolic function were analyzed. Significant correlations were found between sleep quality and metabolic function (|r| = 0.386 ± 0.062, p < 0.05); HRV derived ANS function showed strengthened correlations with metabolic function during sleep period (|r| = 0.474 ± 0.100, p < 0.05); HRV metrics during sleep stages coupled more tightly with clinical indicators of metabolic function [in unstable sleep: |r| = 0.453 ± 0.095, p < 0.05; in stable sleep: |r| = 0.463 ± 0.100, p < 0.05; in rapid eye movement (REM) sleep: |r| = 0.453 ± 0.082, p < 0.05], and showed significant associations with glycemic control in non-linear analysis [fasting blood glucose within 24 h of admission (admission FBG), |r| = 0.420 ± 0.064, p < 0.05; glycated hemoglobin (HbA1c), |r| = 0.417 ± 0.016, p < 0.05].
Conclusions: HRV metrics during sleep period play more distinct role than during awake period in investigating ANS dysfunction and metabolism in T2DM patients, and sleep rhythm based HRV analysis should perform better in ANS and metabolic function assessment, especially for glycemic control in non-linear analysis among T2DM patients.
1 Introduction
Type 2 diabetes mellitus (T2DM) is a chronic hyperglycemia that causes physiological dysfunction and failure of various organs. Autonomic neuropathy is one of the most common complications of diabetes mellitus, which seriously affects the patients’ life quality and brings about significant morbidity and mortality (Maser et al., 2003; Vinik et al., 2003; Ziegler et al., 2008). The dysfunction of autonomic nervous system (ANS) usually manifests first in vagal nerve (the longest parasympathetic nerve in the body, responsible for nearly three-quarters of parasympathetic activity), and damage to vagal nerve leads to resting tachycardia and an overall decrease in parasympathetic tone (Balcıoğlu and Müderrisoğlu, 2015).
Dysfunction of ANS, caused by T2DM, happens in any part of ANS from early stages of diabetes, which damages the cardiovascular, gastrointestinal, genitourinary and neurovascular systems (Pfeifer et al., 1982; Singh et al., 2000; Goit et al., 2012). The function of healthy ANS has obvious circadian rhythm (Cui et al., 2020), and the dysfunction of ANS can result in the loss of circadian rhythm to varying degrees (Liu et al., 2020). Consequently, diabetes is often followed by sleep disorders (Barone and Menna-Barreto, 2011). Meanwhile, sleep disorders accelerate the development of T2DM by worsening the metabolic control, which forms into a vicious spiral (Barone and Menna-Barreto, 2011). Researchers have found that poor sleep quality (Martyn-Nemeth et al., 2018), insufficient or excessive sleep duration, changes of sleep structure (Koren et al., 2011), decreasing sleep efficiency (Hur et al., 2020), increasing OSA severity (Kent et al., 2014), etc. are associated with poor diabetic control. Also, sleep disorders are strongly related to ANS function of T2DM patients (Jordan et al., 2014). The activation of the sympathetic nervous system (SNS) and the unbalance of ANS in patients with T2DM could be mediated through sleep impairments, including nocturnal breathing disturbances (Punjabi and Polotsky, 2005; López-Cano et al., 2019) and sleep curtailment (Leproult et al., 1997; Spiegel et al., 1999).
Heart rate variability (HRV) is considered to be an effective measure of heart-brain interaction and tension of ANS (Montano et al., 2009), and is widely used as a standard method for assessing autonomic function. Studies have found that T2DM patients’ HRV differs significantly from normal people in time-domain, frequency-domain and non-linear analysis (Barone and Menna-Barreto, 2011; Faust et al., 2012; Martyn-Nemeth et al., 2018; Liu et al., 2020). Therefore, HRV analysis is considered an effective, non-invasive auxiliary detection method for T2DM patients’ autonomic dysfunction. Cardiologists tend to analyze 24-h HRV, whereas internists tend to detect whether HRV is abnormal in different physiological states, known as Ewing test (Gerritsen et al., 2001).
Despite the full utilization in clinical fields, HRV analysis in current researches mainly focused on linear analysis of 24-h or 5–15 min electrocardiogram (ECG) collected in different physiological states, which fails to effectively extract the information in long-term signals and neglects the space-time complexity and fractal properties of heart rate time series. Meanwhile, it was a general problem that HRV metrics were accepted as measures of autonomic function without examination of the underlying physiological patterns (Stein and Pu, 2012) and autonomic pathways. Hence, controversial research results were presented. For example, researchers have found that low-frequency component of HRV was not predominantly influenced by SNS, but vagal nerve system (Reyes del Paso et al., 2013), and reflected multiple sources of variability (Martelli et al., 2014), which also challenged the usefulness of the ratio of low-frequency component to high-frequency component (LF/HF); higher HRV were related to disturbed sleep, such as in sleep disorders as negative feature (Stein and Pu, 2012); imparity of HRV metrics reflecting parasympathetic nervous system (PNS), such as high-frequency power (HF) and root mean square of successive of RR interval differences (RMSSD), occurred during different sleep stages (Glos et al., 2014; Zhang et al., 2020).
Therefore, based on the complex associations among sleep, ANS and T2DM, we monitored 24-h ECG, clinical metabolic indicators, sleep quality as well as sleep staging results of T2DM patients. Our objective was to help optimize the early prevention, post-diagnosis and management of T2DM and its complications. Our hypothesis included: 1) HRV derived ANS function is more associated with T2DM patients’ metabolic function during sleep than during awake or the whole 24-h; 2) ANS function alters during sleep cycles, leading to changes on its interaction with metabolic system, thus we analyzed interaction between HRV derived ANS function with metabolic function in each sleep stage; finally, 3) non-linear HRV analysis during sleep cycles may discover more features of ANS function than conventionally used linear HRV analysis, leading to more extensive association with metabolic function and glycemic control.
2 Materials and methods
2.1 Patients
The study was conducted form July 2019 to January 2021. A total of 64 volunteered inpatients with T2DM from Suzhou Science & Technology Town Hospital took part in this study. Exclusion criteria were: type 1 diabetes mellitus, history of stroke, subacute myocardial infarction, kidney or liver transplant, other systemic disorders, current recreational drug or alcohol abuse and morbid obesity (body mass index >40). Because of the small number of female inpatients in the study (only 4), and considering that women in the same age group may be affected metabolically by perimenopausal syndrome, and decreasing the bias brought by sex (Nunan et al., 2010), data from only 60 male inpatients were included in the follow-up studies. This study was approved by the Institutional Review Board of Suzhou Science & Technology Town Hospital (No. IRB2019045). Prior to the study, all subjects were informed of the experimental protocol and precautions, and signed the written informed consent. Demographic characteristics of the subjects are demonstrated in Table 1. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines for cross-sectional studies were followed for study reporting (Cuschieri, 2019).
TABLE 1. Characteristics of the study cohort. Data are presented as mean ± standard deviation (SD) if the variable is normally distributed, otherwise presented as median (p25, p75).
2.2 Clinical indicators of metabolic function
During hospitalization, subjects underwent clinical examinations to obtain their metabolic function and assess their health status. Upon admission, a fasting blood draw and urinary sample were obtained the next morning for routine glucose, lipid and renal panels. Subjects were also involved in investigation for diabetic complications. Clinical indicators of metabolic function analyzed in this study and diabetic complications statistics are demonstrated in Table 1.
Subsequently, clinical indicator outlier rejection was performed where box plot was utilized. The data were rearranged from largest to smallest, with the difference between the upper quartile U and the lower quartile L defined as IQR. The upper bound was set to U+1.5*IQR and the lower bound was set to L-1.5*IQR, and data exceeding the upper and lower bounds can be considered as outliers. Nevertheless, some data that did not appear to be gross errors, i.e., that were outside the range of normal values but might appear in patients who are in dysfunction states, were retained with the aim of reserving more information underneath the pathological data.
An HbA1c level greater than 9.0% indicates poor control of diabetes, which corresponds to a greater risk of complications and the need for more stringent treatment strategies (Inzucchi et al., 2015; Mosenzon et al., 2016; Ye et al., 2016; Garber et al., 2020). Subjects could be assigned to two groups according to the HbA1c level: one group had poor control of diabetes (HbA1c ≥ 9.0%) and the other did not (HbA1c < 9.0%).
2.3 ECG data recording
ECG recordings were collected by an FDA (U.S. Food and Drug Administration) approved ambulatory electrocardiogram monitor (DynaDx Corporation, Mountainview, CA, United States) with a computer-based data-acquisition system. The ECG recording equipment was a single-lead Holter device that can record ECG for over 24 h. Sampling frequency of ECG monitoring was set to 250 Hz. All subjects were monitored in hospital for 24-h, starting at 10 p.m. on the second day of hospitalization. The 24-h ECG data was analyzed for HRV, which will be reported elsewhere. All ECG recordings were carefully checked with noise level, artifacts, and ectopic beats. Seventeen ECG recordings were discarded due to low data quality, too short recording time, or the fact that the inpatients had atrial fibrillation or server arrhythmia.
2.4 ECG data preprocessing
The band-pass Butterworth filter and zero-phase shift filter with a cut-off frequency of 5 Hz–35 Hz were first used to eliminate the baseline drift, power frequency interference, myoelectric interference, motion artifacts and equipment noise.
At present, plenty of R-peak detection methods have been proposed on ECG signal analysis. However, when applying these detectors on ECG signals collected in long-term recordings especially via wearable single-lead ECG devices, the R-peak detection accuracies were usually unsatisfying. To guarantee the accuracy of R-peak detection, a method for extracting high-quality RR intervals proposed in our previous study (Cui et al., 2020) was applied, which combing five commonly used R-peak detection methods including “Pan-Tompkins” (Pan and Tompkins, 1985; Hamilton and Tompkins, 1986), “jqrs” (Johnson et al., 2015), “mteo” (Solnik et al., 2008), “nqrs” (Li and Yang, 2013), and “sixth power” (Dohare et al., 2014). Outlier rejection for RR interval time series was performed, and box plot was utilized to find and reject outliers.
2.5 HRV analysis and metrics
In this study, we evaluated HRV based on the widely-used linear and non-linear HRV metrics (Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology, 1996; Sassi et al., 2015).
2.5.1 Linear HRV analysis
Linear HRV metrics include time-domain metrics and frequency-domain metrics. Time-domain metrics quantify the amount of variability of the time period between successive heartbeats (RR intervals), and include HR-mean (mean of heart rate), SDNN (standard deviation of RR intervals), RMSSD (root mean square of successive RR intervals differences), SDSD (standard deviation of successive RR interval differences), SDANN (the standard deviation of the average value every 5 min in the RR interval) and pNN50 (the percentage of successive RR intervals that differ by more than 50 ms). Frequency-domain metrics estimate the distribution of absolute or relative power at different frequency bands. Generally, the whole frequency band (whose absolute power is noted as TP) is divided into three frequency bands: high-frequency band (0.15–0.40 Hz, whose absolute power is noted as HF), low-frequency band (0.04–0.15 Hz, whose absolute power is noted as LF) and very-low-frequency band (0.00–0.04 Hz, whose absolute power is noted as VLF). Other frequency-domain metrics include LFP (percentage of LF in TP), HFP (percentage of HF in TP) and LF/HF (ratio of LF-to-HF power). The widely used HRV analysis model for assessing ANS homeostasis includes three core statements (Pagani et al., 1986; Malliani et al., 1991; Montano et al., 2009): 1) the power of the high-frequency component reflects cardiac parasympathetic activity; 2) the power of low-frequency component reflects cardiac sympathetic activity; and 3) the LF-HF-ratio reflects sympathetic parasympathetic balance.
2.5.2 Non-linear HRV analysis
Non-linear analysis utilized in this study includes approximate entropy (ApEn) (Pincus, 1991), sample entropy (SampEn) (Chen et al., 2009), fuzzy entropy (FuzzyEn) (Richman and Moorman, 2000; Azami and Escudero, 2016), detrended fluctuation analysis (DFA, alpha1 and alpha2) (Peng et al., 1994), multi-scale entropy (MSE) (Costa et al., 2002) and multi-scale fuzzy entropy (MFE). Calculation steps are introduced in the supplementary material.
ApEn is a rather conventional measure to quantify irregularity and complexity and reflects the probability of new subsequences, and the more complex time series corresponds to larger ApEn. SampEn is a modification of ApEn, and is more accurate than ApEn in the case of less data. FuzzyEn is improved on the basis of SampEn, which is the entropy of a fuzzy set that contains vagueness and ambiguity uncertainties. When calculating ApEn, the embedding dimension m was set to 2, and the similar tolerance r was set to 0.25 times the SD of the sequence in order to find significant difference more clearly; in SampEn and FuzzyEn, r was set to 0.15 times the SD of the sequence.
MSE analyzes the complexity of time series from multiple time scales. The 25 coarse-grained time series are contracted from the original RR time series, in which we obtained 28 indices, including the SampEn value at each coarse-graining scale and 3 complexity metrics, MSEsum5, MSEsum10 and MSEsum20 (defined as the area of the MSE curve at scales 1–5, 1–10 and 1–20 respectively). The MSE curve usually rises rapidly at lower scales, peaks at scale 5, and then declines slowly, reaching a plateau after scale 20, therefore calculating the area of the MSE curve at scales 1–5 and 1–20 can portray the complexity characteristics of this curve at high and low resolutions. Similarly, MFE were calculated based on FuzzyEn.
DFA is a non-linear fractal analysis tool to discover potential self-similarity in the series and quantify the fractal scale of the time series. The α1 and α2 portray the short- and long-range correlation respectively.
2.6 Sleep quality assessment
Subjective sleep quality was assessed by Pittsburgh Sleep Quality Index (PSQI) and a brief sleep log was used to record sleep duration for the studied night. PSQI includes multiple sleep-related variables over the preceding month, using Likert and open-ended response formats (Spira et al., 2011). The PSQI yields seven component scores: subjective sleep quality, sleep latency, sleep duration, habitual sleep efficiency, sleep disturbances, sleep medication, and daytime dysfunction. Component scores range from 0 to 3 and are summed to obtain a global score, which ranges from 0 to 21. Higher scores suggest greater sleep disturbance (Buysse et al., 1989). During 24-h ECG monitoring, 53 subjects filled in PSQI questionnaire.
Objective sleep quality and sleep staging were assessed by ECG-based cardiopulmonary coupling (CPC) analysis (Thomas et al., 2005), which has been cited in more than 260 publications (Yeh et al., 2008; Schramm et al., 2009; Thomas et al., 2010), and was introduced in the famous international sleep monograph “Principle and Practice of Medicine” (the sixth edition). The ECG recordings during sleep at night were extracted for sleep analysis in this study. CPC analysis is based on mathematical analysis of the coupling between HRV and the respiratory modulation of QRS waveform on a beat-to-beat basis. Major physiological sleep states derived from CPC analysis include stable sleep (indicated by high-frequency coupling, or HFC), unstable sleep (indicated by low frequency coupling, or LFC), and rapid eye movement (REM) or wakeful states (indicated by very low frequency coupling, or VLFC) (Thomas et al., 2005). In this study, 38 CPC sleep reports were available. Sleep quality metrics analyzed in this study are demonstrated in Table 2.
TABLE 2. Sleep quality metrics analyzed in this study. Data are presented as mean ± SD if the variable is normally distributed, otherwise presented as median (p25, p75).
2.7 Statistical analysis
The data were analyzed by SPSS version 27.0 for Windows. The statistical analysis included significance analysis and correlation analysis. The normality of distribution of continuous variables was firstly tested by the Shaprio-Wilk test (significance level α = 0.05). For two normally distributed continuous variables, the Pearson correlation analysis was utilized to estimate the level of correlation, otherwise the Spearman analysis. Partial correlation analysis was also involved in the assessment of correlations between variables. For normally distributed variables, group differences were compared by independent sample t-test (2 groups) and one-way ANOVA (3 or more groups), otherwise Mann-Whitney U test (2 groups) or Kruskal-Wallis test (3 or more groups) were used. A value of p < 0.05 was considered significant. Missing data was completely at random, and pairwise deletion was conducted during the statistical analyses.
There are 25 scales and 3 complexity metrics for MSE and MFE respectively, and in order to demonstrate their overall correlations with clinical indicators in a straightforward way, we define the correlations between MSE and clinical indicators as mean ± SD of the correlations between SampEn at specific scales and clinical indicators. The coefficients of the same significant level (p < 0.05 or p < 0.01) between SampEn at certain scale and clinical indicators were selected to calculate the mean and SD. Similarly, the correlations between MFE and clinical indicators can be calculated. As for complexity metrics, the correlations between MSE complexity and clinical indicators are defined as mean ± SD of the correlations between specific complexity metrics (MSEsum5, MSEsum10 and MSEsum20) of the same significant level and clinical indicators, and similarly the correlations between MFE complexity and clinical indicators could be calculated.
3 Results
3.1 Correlations between sleep quality metrics and clinical indicators
Table 3 demonstrates the correlation coefficients between sleep quality metrics and clinical indicators. PSQI correlated rather weakly with DBP (r = −0.296, p = 0.035), Hb (r = −0.291, p = 0.042) and LDL-C (r = −0.334, p = 0.023), while sleep quality metrics calculated by CPC analysis were more significantly correlated with clinical indicators that were related to metabolic functions of T2DM patients. For example, liver function indicators AST/ALT was associated with TST (r = −0.575, p < 0.001) and UST (r = 0.424, p = 0.009), renal function indicator BUN was associated with TST (r = 0.488, p = 0.005), UST (r = 0.479, p = 0.006), ALUS (r = 0.442, p = 0.011) and USP (r = 0.392, p = 0.027). Moreover, TST and UST were significantly correlated with many clinical indicators, including admission FBG, Hb, ALT, AST/ALT, BUN, TG, HDL-C and LDL-C (e.g., TST & admission FBG, r = −0.383, p = 0.025; TST & TG, r = −0.405, p = 0.019; UST & LDL-C, r = −0.400, p = 0.017).
TABLE 3. Correlation coefficients between sleep quality metrics and clinical metabolic indicators.
3.2 Significance of HRV metrics during sleep
Table 4 demonstrates the mean values of HRV metrics for 40 T2DM patients, which showed significant differences between sleep and awake periods, especially non-linear HRV metrics. During sleep period, subjects had smaller values of HR-mean, SDANN, TP, LF/HF and LFP, while their RMSSD, SDSD, SampEn, FuzzyEn, MSEsum5, MSEsum10, MSEsum20, MFEsum5, MFEsum10 and MFEsum20 were higher than those during awake period.
TABLE 4. HRV metrics for subjects during different recording period.
Table 5 demonstrates correlation coefficients between HRV metrics in 24-h/sleep/awake periods and clinical indicators. During 24-h period, significant correlations could be observed between liver function indicators ALT, AST/ALT and HRV metrics (e.g., ALT & HF, r = 0.688, p < 0.001; AST/ALT & LF, r = −0.707, p < 0.001; AST/ALT & HF, r = −0.721, p < 0.001), and between renal function indicator BUN and HRV metrics (e.g., RMSSD, r = −0.482, p = 0.006; pNN50, r = −0.527, p = 0.002; HF, r = −0.517, p = 0.003); as to HRV metrics, HFP was significantly correlated with many clinical indicators (ALT, r = 0.466, p = 0.004; AST/ALT, r = −0.510, p = 0.001; BUN, r = −0.459, p = 0.007), and since HF component mainly reflects PNS activity (Pagani et al., 1986; Malliani et al., 1991; Montano et al., 2009), it might suggest that PNS function is closely linked to physiological status in T2DM patients. During sleep period, significant correlations could be observed between many HRV metrics and DBP (e.g., SampEn, r = 0.540, p < 0.001; MSE complexity, r = 0.571 ± 0.054, p < 0.001), Hb (e.g., MSE complexity, r = 0.513 ± 0.031, p = 0.001 ± 0.001; DFA-α2, r = −0.524, p < 0.001), ALT (e.g., HF, r = 0.638, p < 0.001; SampEn, r = 0.645, p < 0.001), AST/ALT (e.g., FuzzyEn, r = −0.693, p < 0.001; MFE complexity, r = −0.711 ± 0.021, p < 0.001), and BUN (e.g., pNN50, r = −0.556, p = 0.001; SampEn, r = −0.541, p = 0.002); meanwhile, the significant correlations were observed between many clinical indicators and non-linear HRV metrics (e.g., ALT & SampEn, r = 0.645, p < 0.001; ALT & FuzzyEn, r = 0.645, p < 0.001; AST/ALT & MFE complexity, r = −0.711 ± 0.021, p < 0.001). However, during awake period, the correlations between HRV metrics and clinical indicators decreased observably, and HRV metrics mainly correlated with liver function indicators ALT (e.g., LF, r = 0.632, p < 0.001; HF, r = 0.628, p < 0.001) and AST/ALT (e.g., LF, r = −0.592, p < 0.001; HF, r = −0.606, p < 0.001).
TABLE 5. Correlation coefficients between HRV metrics during different recording period and clinical indicators.
It could be seen that the overall level of correlation coefficients was significantly higher during sleep period than 24-h and awake periods. Correlations between variability-related HRV metrics (RMSSD and pNN50) and clinical indicators, including liver function (ALT, AST/ALT) and renal function (BUN), were more significant in sleep period than in 24-h and in awake period; significant correlations between complexity-related HRV metrics (SampEn, FuzzyEn, MSE complexity, MFE complexity) and clinical indicators were found mainly in sleep period.
3.3 HRV metrics in different sleep stages are more associated with T2DM clinical indicators
To investigate the relationship between HRV metrics (ANS function) in different sleep stages and the clinical indicators (metabolic function) of T2DM patients, the RR interval time series during sleep period were segmented according to sleep stages, including unstable sleep, stable sleep and REM sleep, and series from the same sleep stage were stitched together as a whole. Significant differences were observed among HRV metrics in different sleep stages. For example, the LF/HF reached its maxima in unstable sleep (4.04 ± 3.31), followed by REM sleep (3.36 ± 2.27) and its minima in stable sleep (1.45 ± 0.98). More details are shown in Supplementary Table S2.
The correlation coefficients between HRV metrics in different sleep stages and clinical indicators are demonstrated in Table 6. Linear HRV metrics, especially RMSSD, SDSD, pNN50, LF, HF and HFP were correlated significantly with Hb, ALT, AST/ALT and BUN (e.g., in unstable sleep: SDSD & ALT, r = 0.665, p < 0.001; RMSSD & AST/ALT, r = −0.710, p < 0.001; HF & AST/ALT, r = −0.716, p < 0.001; pNN50 & BUN, r = −0.533, p = 0.002; in stable sleep: HFP & AST/ALT, r = −0.556, p < 0.001; in REM sleep: LF & Hb, r = 0.510, p = 0.001). In non-linear HRV analysis, more metrics were tightly associated with clinical indicators, and the correlations were more significant with those in linear analysis (63 correlations with |r|>0.5, compared with 44 correlations with |r|>0.5 in linear analysis), especially the DFA metric α2 and the complexity metrics MSE and MFE (e.g., in unstable sleep: DFA-α2 & admission FBG, r = −0.584, p < 0.001; MFE & AST/ALT, r = −0.614 ± 0.052, p < 0.001; MFE complexity & AST/ALT, r = −0.680 ± 0.014, p < 0.001; in stable sleep: SampEn & ALT, r = 0.786, p < 0.001; MSE complexity & ALT, r = 0.596 ± 0.076, p < 0.001; MSE & AST/ALT, r = −0.571 ± 0.132, p = 0.002 ± 0.003; FuzzyEn & AST/ALT, r = −0.760, p < 0.001). Compared with the correlation with linear HRV metrics, correlations between admission FBG, DBP, Hb, ALT, AST, AST/ALT, GGT, BUN and non-linear HRV metrics increased considerably (e.g., in unstable sleep: DFA-α2 & admission FBG, r = −0.584, p < 0.001; in stable sleep: SampEn & ALT, r = 0.786, p < 0.001; SampEn & AST, r = 0.630, p < 0.001; FuzzyEn & AST/ALT, r = −0.760, p < 0.001; SampEn & BUN, r = −0.625, p < 0.001; in REM sleep: MSE complexity & DBP, r = 0.560 ± 0.046, p < 0.001; DFA-α2 & Hb, r = −0.549, p < 0.001; MSE complexity & GGT, r = 0.494 ± 0.009, p = 0.002 ± 0.000), and significant association between HbA1c, N%, PLT, HDL-C, LDL-C, UMA and HRV metrics were mainly observed in non-linear analysis (e.g., in unstable sleep: FuzzyEn & HbA1c, r = 0.506, p = 0.001; in stable sleep: MFE & N%, r = −0.474, p = 0.003; in REM sleep: MFE & PLT, r = 0.515, p = 0.001; DFA-α2 & LDL-C, r = −0.495, p = 0.003; MFE & UMA, r = −0.512 ± 0.008, p = 0.003 ± 0.001).
TABLE 6. Correlation coefficients between HRV metrics in three sleep stages and clinical indicators.
The findings hinted that HRV metrics in different sleep stages were more associated with T2DM clinical indicators. During each sleep stage, the correlation coefficients between HRV metrics and clinical indicators demonstrated a higher overall level compared with those during sleep period. Compared with linear HRV metrics, non-linear HRV metrics significantly correlated with more clinical indicators, including LDL-C and UACR.
3.4 Non-linear HRV metrics during sleep are associated with glycemic control
Interestingly, non-linear HRV metrics in different sleep stages were significantly associated with glycemic control, while few weak correlations were found between HRV metrics of 24-h and awake periods and glycemic control indicators. Table 7 demonstrates correlation coefficients between non-linear HRV metrics in three sleep stages and glycemic control indicators including admission FBG, discharge FBG and HbA1c. Most of non-linear HRV metrics were significantly correlated with admission FBG (e.g., in unstable sleep: DFA-α2, r = −0.584, p < 0.001; MSE, r = 0.466, p = 0.006; MSE complexity, r = 0.440, p = 0.009), yet had no or only weak correlation with discharge FBG. Both single- and multiple-scale sample entropies showed significant correlations with HbA1c (in unstable sleep: SampEn, r = 0.438, p = 0.025; MSE, r = 0.411 ± 0.016, p = 0.038 ± 0.008; MSE complexity, r = 0.409 ± 0.017, p = 0.039 ± 0.009; in stable sleep: MSE, r = 0.389, p = 0.050; in REM sleep: SampEn, r = 0.424, p = 0.031; MSE, r = 0.418 ± 0.022, p = 0.035 ± 0.010; MSE complexity, r = 0.442 ± 0.021, p = 0.024 ± 0.007).
TABLE 7. Correlation coefficients between non-linear HRV metrics in three sleep stages and clinical glycemic control indicators including admission FBG, discharge FBG and HbA1c.
Furthermore, the subjects were divided into two groups according to whether the subject had poor control of diabetes (HbA1c ≥ 9.0%, n = 10) or not (HbA1c < 9.0%, n = 16). Table 8 demonstrates between-group differences in non-linear HRV metrics. Significant differences were observed in non-linear metrics in unstable and stable sleep (in unstable sleep: DFA-alpha2, p = 0.002; MSEsum5, p = 0.020; MSEsum10, p = 0.017; MSEsum20, p = 0.037; in stable sleep: SampEn, p = 0.020; DFA-alpha2, p = 0.026; MSEsum5, p = 0.002; MSEsum10, p = 0.004; MSEsum20, p = 0.024).
TABLE 8. Difference in non-linear HRV metrics between the group had poor control of diabetes (HbA1c ≥ 9.0%) or did not (HbA1c < 9.0%). Data are presented as mean ± SD.
3.5 Association among glycemic control, non-linear HRV analysis and sleep
We also applied partial correlation analysis between HbA1c, admission FBG and non-linear HRV metrics in sleep stages as controlling the effect of sleep quality metrics, in order to briefly discover the interaction among glycemic control, non-linear HRV analysis and sleep.
Table 9 demonstrates partial correlations between HbA1c and non-linear HRV metrics in unstable and REM sleep. When nullifying the effect of SST, ALSS and SSP, the correlations between HbA1c and non-linear HRV metrics significantly strengthened (e.g., SampEn in REM sleep: zero order, r = 0.424, p = 0.031; nullified by SST, r = 0.547, p = 0.005; nullified by ALSS, r = 0.549, p = 0.004; nullified by SSP, r = 0.569, p = 0.003). However, the correlations between HbA1c and non-linear HRV metrics decreased and even were lost when corrected for TST (e.g., SampEn in REM sleep: zero order, r = 0.424, p = 0.031; nullified by TST, r = 0.402, p = 0.046; MSEsum10 in REM sleep: zero order, r = 0.428, p = 0.029; nullified by TST, r = 0.400, p = 0.048). The associations were mostly retained when corrected for other sleep metrics (e.g., MSEsum5 in unstable sleep: zero order, r = 0.421, p = 0.032; nullified by AHI, r = 0.428, p = 0.033; nullified by ALUS, r = 0.417, p = 0.038; nullified by USP, r = 0.424, p = 0.035).
TABLE 9. Partial correlation analysis between HbA1c and non-linear HRV metrics in unstable and REM sleep.
Table 10 demonstrates partial correlations between admission FBG and non-linear HRV metrics in 3 sleep stages. When nullifying the effect of SST, the correlations between HbA1c and non-linear HRV metrics at sleep stages significantly strengthened (e.g., SampEn in unstable sleep: zero order, r = 0.400, p = 0.019; nullified by SST, r = 0.464, p = 0.007; FuzzyEn in unstable sleep: zero order, r = 0.378, p = 0.027; nullified by SST, r = 0.466, p = 0.006). When nullifying the effect of TST, the correlations between HbA1c and non-linear HRV metrics decreased and even were lost (e.g., MSEsum5 in unstable sleep: zero order, r = 0.440, p = 0.009; nullified by TST, r = 0.346, p = 0.049; DFA-α2 in REM sleep: zero order, r = −0.524, p = 0.001; nullified by TST, r = −0.465, p = 0.006). The associations between admission FBG and non-linear HRV metrics, especially DFA-α2, were retained when corrected for other sleep metrics (e.g., DFA-α2 in unstable sleep: zero order, r = −0.584, p < 0.001; nullified by AHI, r = −0.585, p < 0.001; nullified by USP, r = −0.582, p < 0.001; nullified by SSP, r = −0.586, p < 0.001; DFA-α2 in REM sleep: zero order, r = −0.524, p = 0.001; nullified by RST, r = −0.537, p = 0.001; nullified by ALUS, r = −0.538, p = 0.001).
TABLE 10. Partial correlation analysis between admission FBG and non-linear HRV metrics in 3 sleep stages.
4 Discussion
This study focused on the relationships among sleep, HRV derived ANS function and metabolic function in T2DM, and the main findings were: 1) metabolic function was significantly correlated with sleep quality; 2) HRV derived ANS function was different between sleep and awake, showing strengthened and distinct distributions of correlations with metabolic function compared to 24-h analysis; 3) HRV derived ANS function and its correlation with clinical indicators of metabolic function altered during sleep cycles; 4) HRV analysis during sleep was strongly associated with clinical indicators of metabolic function, and this also applied to the glycemic control in patients with T2DM in non-linear analysis. These findings demonstrated that sleep rhythm based HRV analysis should be emphasized in ANS and metabolic function assessment and the rationality of non-linear HRV analysis for further application among T2DM patients and investigation of potential interactions among ANS, sleep and metabolic function (Lombardi, 2011; Perkiömäki, 2011).
In sleep quality assessment, PSQI had few and weak correlations with clinical indicators, while TST and UST correlated with more clinical indicators. The results suggest limitations of PSQI from its objectivity and indicate the important role of sleep in regulation of metabolic function. Researchers have found that circadian rhythms and sleep regulate hormones and lipids involved in energy metabolism (VanCauter et al., 1997; Scheer et al., 2009; Nguyen and Wright, 2010; Spiegel et al., 2011; Markwald et al., 2013), and disruption of sleep and circadian rhythms is progressively apparent as an influence factor to damaged physiological function, disease processes and metabolic dysregulation (Spiegel et al., 2005; Turek et al., 2005; Marcheva et al., 2010; Markwald and Wright, 2012).
The enhanced association between ANS and metabolic function during sleep is consistent with previously published results. Researchers have found that most of the biological functions changes during sleep compared to wake, such as heart rate, arterial blood pressure and etc., where ANS plays a vital role (Tobaldini et al., 2017). Results in Table 5, 6 showed that HRV metrics during sleep were associated with clinical indicators of metabolic function more significantly and broadly, compared to 24-h HRV, and there were significant differences between HRV metrics during sleep and awake. Thus, sleep might be a good model to investigate ANS activity and the fluctuation of ANS function caused by intrinsic factors, especially circadian rhythm, without influence factors of daytime activities such as eating behaviour, physical movement, etc. Meanwhile, studies have shown that HRV metrics are significantly influenced by routine life states and ECG recording time (morning, noon, evening), and the RR values are rather stable during night sleeping (Cui et al., 2020), so we analysed the RR interval time series during awake period (i.e., during daytime) as a whole, and at the same time the time of waking and sleeping of subjects was rather uniform, which may decrease bias brought by recording time.
Significant differences observed in many HRV metrics among different sleep stages suggest the dominance of SNS in unstable sleep and the dominance of PNS in stable sleep (Pagani et al., 1986; Malliani et al., 1991; Montano et al., 2009). Researchers have found that ANS tone and modulation are profoundly influenced by sleep related mechanisms (Somers et al., 1993; Trinder et al., 2001; Tobaldini et al., 2014; de Zambotti et al., 2018), while cardiovascular function within each sleep stage tends not to vary across the night (Silvani and Dampney, 2013). For ANS impose regulatory control over the cardiovascular system by SNS and PNS, which leads to HRV, the changes of HRV metrics in different sleep stages and their correlations with clinical indicators may associate more closely with the variability of ANS function during sleep cycles. Researchers have conducted HRV analysis during sleep to confirm different ANS function patterns at sleep stages (Scholz et al., 1997; Brandenberger et al., 2001; Trinder et al., 2001; Stein and Pu, 2012). Studies have shown that, since most of HRV metrics applied in this study were based on the assumption of stationary series (Akselrod et al., 1981; Pagani et al., 1986; Pincus et al., 1993; Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology, 1996; Voss et al., 1996; Cooke et al., 1999; Porta et al., 2000b; Cysarz et al., 2000; Richman and Moorman, 2000; Wessel et al., 2000; Porta et al., 2001), non-stationarities of HRV sequences could have a great influence on the assessment of the balance of ANS by producing a bias toward a higher sympathetic modulation and a lower vagal modulation (Magagnin et al., 2011). We applied Augmented Dickey-Fuller test to all the HRV sequences for stationarity test, and found sequences during 24-h, daytime and sleep were mostly stationary. However, during sleep cycles, the sequences were mostly non-stationary. Hence, the distinct distribution pattern presented in analysis during sleep cycles may demonstrate the shift of sympathovagal balance toward sympathetic predominance since non-stationarities of the HRV sequences will highlight the shift (Magagnin et al., 2011). With enhanced and clarified interaction between ANS function and HRV analysis in each sleep stage, the underlying autonomic pathways (Stein and Pu, 2012) of different HRV metrics may be better distinguished from the remarkably strengthened correlations between clinical indicator of metabolic function and HRV metrics with distinct distribution patterns compared to traditional analysis of sleep quality and HRV, which is crucial to reveal the interaction among sleep, ANS and HRV.
Compared to linear analysis, non-linear HRV metrics in different sleep stages were associated with more T2DM clinical indicators. According to the development of chaotic theory and non-linear dynamics, it is now generally recognized that heart beat intervals (RR intervals) are non-linear and non-stationary time series, caused by complex interactions between physiological systems (Peng et al., 1995; Cerutti et al., 2007). Traditional methods in time- and frequency-domain may not be able to detect subtle but important changes in RR series (Costa et al., 2005). Currently, HRV frequency-domain metrics are extensively used to assess the function of ANS components, leading to controversial results (Stein and Pu, 2012; Reyes del Paso et al., 2013; Martelli et al., 2014). In this study, most comprehensive and strong correlations between HRV derived ANS function and metabolic function with distinct distribution patterns were observed in non-linear analysis, demonstrating the rationality of non-linear analysis in assessing ANS function.
However, a part of non-linear HRV metrics we chose are somehow influenced by some biases. ApEn is a biased statistic, and the bias arises from two separate sources. First, in the calculation of the correlation integral
It is likewise not negligible that the impact of the data length on the ability to predict the health status of ANS (Laborde et al., 2017). Long-term 24-h HRV recordings are acknowledged as the “gold standard” for clinical HRV assessment (Shaffer et al., 2014) and achieve greater predictive power than short-term measurements (Bigger et al., 1989; Fei et al., 1996; Nolan et al., 1998; Kleiger et al., 2005). However, increased HRV may correspond to longer recordings instead of the change of ANS function, and it is inappropriate to compare metrics like SDNN when they are calculated from epochs of different length (Saul et al., 1988; Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology, 1996). Meanwhile, since mechanisms responsible for heart period modulations of a certain frequency, especially LF and HF power components, may not remain unchanged during long-term recording (Furlan et al., 1990), the interpretation of the results of frequency analysis is less well defined, and they may obscure the detailed information about autonomic modulation of RR intervals that is available in shorter recordings (Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology, 1996).
In this study, surrogate tests (Porta et al., 2000a) was applied to all the HRV sequences analysed. The percentages of non-linear dynamics were: 1) 86.05% in 24-h HRV sequences; 2) 90.70% during awake period, and 39.53% during sleep period; 3) 52.63% during stable sleep, 55.26% during unstable sleep, and 68.42% during REM sleep. Therefore, non-linear dynamics in HRV sequences decreased during sleep, possibly due to the suppression of regulation of respiratory rhythm (Mador and Tobin, 1991; González et al., 2000; Francis et al., 2002; Jo et al., 2005), leading to stable patterns. And when HRV sequences were segmented based on sleep stages, the percentage of non-linear dynamics elevated, suggesting better manifestation of the non-linear components in HRV sequences during sleep cycles and demonstrating the effectiveness of non-linear analysis. Additionally, the distribution of non-linear features during different sleep stages, that the maximum of non-linear dynamics appeared during REM sleep, was consistent with non-linear respiratory dynamics (Sako et al., 2001).
However, researchers have found that in the analysis of short-term, stationary and successive HRV sequences, non-linear model-free and linear model-based conditional entropy correlate closely and have similar performance in the assessment of cardiac control (Porta et al., 2017). The conflicting conclusions may be a consequence of differences in stationary and non-linear features of HRV sequences utilized in analysis (Magagnin et al., 2011). Previous study has shown that short-term HRV in healthy young adults at rest is mainly linear (Porta et al., 2007b), while in this study HRV sequences during sleep cycles were mainly non-stationary and had non-linear features, demonstrating that the performance of non-linear analysis is related to the feature of HRV sequences.
Sleep rhythm based HRV analysis presents us significant correlations between ANS and metabolic function of T2DM patients covered by limitations of traditional measures, demonstrating its efficiency on distinguishing potential dynamic characteristics of ANS function and its correlation with metabolic function in different sleep stages. Hence, by analysing the correlation between ANS function patterns and metabolic function of T2DM patients, our findings propose a promising strategy for the evaluation of T2DM patients’ metabolic function and physical condition to improve the early prevention, post-diagnosis and management of T2DM and its complications based on HRV analysis during sleep cycles.
Significant correlations were found among DBP and HRV metrics during sleep (see Table 5), and non-linear HRV metrics during sleep cycles significantly correlated with DBP with distinct distribution patterns (see Table 6), suggesting strong association between cardiac baroreflex control and the fluctuation of HRV metrics, especially non-linear ones, in T2DM patients. According to the SBP (136 ± 14 mmHg) and DBP (85 ± 11 mmHg) of the study cohort shown in Table 1 and “Clinical practice guidelines for the management of hypertension in China” (National Center for Cardiovascular Diseases et al., 2022), plenty of the subjects had hypertension of varying degrees, suggesting impairment of cardiac baroreflex control, associated with impairment of ANS and the presence of other neuropathies (Frattola et al., 1997; Tank et al., 2001; Ruiz et al., 2005; Vinik and Ziegler, 2007), which may be the effect of T2DM on baroreflex control through an impairment of vagal control (de Moura-Tonello et al., 2016), and the significant correlations support the finding that cardiovascular variability indexes are sensitive to autonomic dysfunction in T2DM patients in early stages (de Moura-Tonello et al., 2016). These findings demonstrate the effectiveness of non-linear HRV analysis and promising strategies for the early detection of autonomic dysfunction.
ANS plays an important role in the exchange of metabolic information between organs and regulation on peripheral metabolism (Yamada and Katagiri, 2007). In addition to affecting the activity of ANS, sleep also exerts a significant regulatory effect on glycemic control by influencing the balance and levels of hormones including leptin, growth hormone-releasing peptide, insulin and cortisol (Balbo et al., 2010), which enhances the association between HRV metrics during sleep and glycemic control of T2DM patients. Researchers have found that 24-h HRV metrics in time- and frequency-domain significantly decrease in T2DM patients with good glycemic control (HbA1c < 7.0%) and poor glycemic control (HbA1c > 7.0%), compared to healthy controls (Yu et al., 2020).
In this study, significant associations were only found between MSE complexity during sleep period and HbA1c, and significant group differences between subjects that had poor control of diabetes (HbA1c ≥ 9.0%, n = 10) and did not (HbA1c < 9.0%, n = 16) were observed in many non-linear metrics during unstable and stable sleep. Moreover, non-linear HRV metrics during sleep were significantly associated with FBG level, which is one of the important indicators for glycemic compliance (Association, 2013). The results indicates that non-linear HRV metrics during sleep may better manifest neuronal dysfunction caused by hyperglycemia, which induces oxidative stress and toxic glycosylation products (Pop-Busui, 2010), highlighting the importance of non-linear analysis during sleep.
Interestingly, higher HbA1c level was associated with higher entropy related metrics (see Table 8), the elevation of complexity may indicate worse metabolic function due to poor diabetic control (Costa et al., 2005). Additionally, in each sleep stage, the correlations between non-linear HRV metrics and FBG changed significantly before and after hospitalization, possibly due to the immediate impact of hospitalization on blood glucose control, reflecting the sensitivity of non-linear HRV assessment on glycemic control.
Due to the complex interaction among sleep, diabetic control and ANS function that sleep have a significant effect on both ANS and metabolic function, especially glycemic control of T2DM patients (Koren et al., 2011; Jordan et al., 2014; Kent et al., 2014; Martyn-Nemeth et al., 2018; Hur et al., 2020), partial correlation analysis was required to manifest the association between non-linear HRV metrics in sleep stages and glycemic control indicators clearly. With the nullification of sleep quality metrics, most of the correlations retained, which supported the significant correlations acquired by simple correlation analysis. Whereas the association significantly decreased or even vanished when TST was corrected, indicating the significant role of sleep quality in adjusting ANS and glycemic function of T2DM patients (Jordan et al., 2014; Martyn-Nemeth et al., 2018). Additionally, when stable sleep related metrics, especially SST, were nullified, the correlation between non-linear HRV metrics in sleep stages and glycemic control indicators significantly strengthened, demonstrating the functional significance of stable sleep stage. As well as stable sleep stage (slow wave sleep) is widely recognized as the most restorative of all sleep stages (Van Cauter et al., 2008) and when several important physiological activities occur specifically (Somers et al., 1993), researchers have found that stable sleep associates with brain glucose metabolism (Zoccoli et al., 2002) and suppression of stable sleep without any reduction in TST associates with poor diabetic control and increased risk of T2DM (Tasali et al., 2008).
In our original study, we also used RR interval time series from four time periods (late night, 3:00–4:00; dawn, 6:00–7:00; after lunch, 13:00–14:00; after dinner, 19:00–20:00) to perform HRV analysis and correlation analysis with metabolic indicators (see Supplementary Table S4–S7). Although distinct correlations were found between HRV metrics and metabolic functions during specific time periods, such as correlations between HRV and BUN, DBP or Hb at dawn, between HRV and HDL-C or LDL-C after lunch, and between metabolic indicators and LF/HF or LFP after dinner, the overall results were not satisfactory and systematic enough to be included in the main text. This may be due to the different physiological states of individuals at specific time periods and is influenced by meals and daytime activities.
To the best of our knowledge, few studies have applied time-, frequency-domain and non-linear analysis on HRV sequences in each sleep stage to inquire the association among ANS function, sleep cycle and metabolic function. We compared our results with those of using traditional methods in all-night sleep analysis and short-term linear HRV analysis. Representative results are shown in Table 11 in summary. In comparison with studies of correlations between all-night sleep and metabolic function, conducted among subjects without T2DM (Koren et al., 2011; Martyn-Nemeth et al., 2018; Feng et al., 2021), correlations found in this study were extensively stronger and broadened with a similar or notably smaller sample size. In comparison with studies of correlations between short-term linear HRV analysis and metabolic function, conducted among T2DM patients (Bhati et al., 2019), correlations found in this study were extensively stronger and broadened with a similar sample size, suggesting the potential of HRV analysis during sleep cycles in T2DM patients. However, when compared with a study with larger sample size (Balikai et al., 2022), correlations found in this study were more comprehensive but slightly weaker, which may be due to the limited sample size, exclusion of female subjects, differences in the physiological status of the subjects including duration and severity of diabetes since our subjects were all patients during hospitalization, the bias introduced by non-linear metrics including ApEn and MSE via the contamination of variance and suboptimal filtering procedure (Pincus and Huang, 1992; Nikulin and Brismar, 2004; Porta et al., 2007a; Valencia et al., 2009; Delgado-Bonal and Marshak, 2019) and the selection of long-term HRV sequences since metrics including LF, HF and SDNN would be affected by length of the sequences (Saul et al., 1988; Furlan et al., 1990; Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology, 1996).
TABLE 11. Comparison of related studies on the interactions among sleep, HRV and metabolic function: subjects, methods, results and conclusion.
Therefore, this study provides new insights into research on the interaction among sleep, ANS and metabolic function, for T2DM patients, especially for diabetic control. Researchers have found that the dysfunction of T2DM patients’ ANS manifests first in decreased parasympathetic activity (Pfeifer et al., 1982; Singh et al., 2000; Goit et al., 2012), even before clinical symptoms. The effective use of HRV analysis during sleep cycles to extract dynamic features of PNS and SNS in each sleep stage provides support for early detection of autonomic neuropathy. In addition, the extensively strong correlations between non-linear HRV metrics and glycemic control indicators including FBG and HbA1c demonstrates the validity of HRV analysis in each sleep stage, which can be another promising strategy for diabetic control and early diagnose of T2DM with less computational costs.
For further research, this study proposed to select HRV sequences based on physiological rhythms can locally reveal the dynamic features of HRV and magnify the potential biological origins and kinetic mechanisms of HRV under the influence of many external factors, such as respiration, circadian rhythm, sympathetic and vagal activity, as well as intrinsic factors, such as ion channel fluctuations and other molecular fluctuations, while reducing the computational cost and elevating local resolution of long term signals, providing new insights into long term physiological signal processing.
There are several potential limitations of this study. Mainly middle-aged male inpatients were included in the analysis, who do not represent the whole spectrum of T2DM patients, to eliminate the influence of climacterium on metabolic function. Although we did a further outlier rejection to reduce the influences of specific situations, there are still diversities between people with different physiological factors (age, gender, health status, course of disease, etc.) that were not included in this study, requiring further studies. Despite of the non-linear dynamics in HRV sequences appeared in our study, researchers have found that short-term heart period variability in healthy young adults at rest is mainly linear (Porta et al., 2007b). Since there weren’t thorough restrictions on subjects’ daily behaviour during hospitalization, further studies are required to investigate the production of non-linear components in long-term HRV sequences by applying forcing input to subjects. Lastly, due to the limited sample size, we did not build more accurate statistical models of HRV, ANS, sleep and metabolic function to further investigate their potential biological origins and kinetic mechanisms.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by the Institutional Review Board of Suzhou Science & Technology Town Hospital (No. IRB2019045). The patients/participants provided their written informed consent to participate in this study.
Author contributions
WC, LT, and XC conceived and designed the study. WC and HC performed the experiment and analyzed the data. WC and HC drafted the manuscript. XC and ZM supervised the analysis, reviewed and editing the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported in part by National Key Research and Development Program of China under Grant 2018YFC2001100, in part by the National Natural Science Foundation of China under Grant 82274631, 61807007.
Acknowledgments
The authors thank all the subjects participated in this study, and the members from Suzhou Science & Technology Town Hospital for helping with data collection.
Conflict of interest
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.
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/fphys.2023.1157270/full#supplementary-material
References
Akselrod, S., Gordon, D., Ubel, F. A., Shannon, D. C., Berger, A. C., and Cohen, R. J. (1981). Power spectrum analysis of heart rate fluctuation: A quantitative probe of beat-to-beat cardiovascular control. Science 213(4504), 220–222. doi:10.1126/science.6166045
Association, A. D. (2013). Diagnosis and classification of diabetes mellitus. Diabetes Care 37 (1), S81–S90. doi:10.2337/dc14-S081
Azami, H., and Escudero, J. (2016). Matlab codes for "refined multiscale fuzzy entropy based on standard deviation for biomedical signal analysis. University of Edinburgh, School of Engineering, Institute for Digital Communications.
Balbo, M., Leproult, R., and Van Cauter, E. (2010). Impact of sleep and its disturbances on hypothalamo-pituitary-adrenal Axis activity. Int. J. Endocrinol. 2010, 759234. doi:10.1155/2010/759234
Balcıoğlu, A. S., and Müderrisoğlu, H. (2015). Diabetes and cardiac autonomic neuropathy: Clinical manifestations, cardiovascular consequences, diagnosis and treatment. World J. Diabetes 6 (1), 80–91. doi:10.4239/wjd.v6.i1.80
Balikai, F. A., Javali, S. B., Shindhe, V. M., Deshpande, N., Benni, J. M., Shetty, D. P., et al. (2022). Correlation of serum HDL level with HRV indices using multiple linear regression analysis in patients with type 2 diabetes mellitus. Diabetes Res. Clin. Pract. 190, 109988. doi:10.1016/j.diabres.2022.109988
Barone, M. T. U., and Menna-Barreto, L. (2011). Diabetes and sleep: A complex cause-and-effect relationship. Diabetes Res. Clin. Pract. 91 (2), 129–137. doi:10.1016/j.diabres.2010.07.011
Bhati, P., Alam, R., Moiz, J. A., and Hussain, M. E. (2019). Subclinical inflammation and endothelial dysfunction are linked to cardiac autonomic neuropathy in type 2 diabetes. J. Diabetes & Metabolic Disord. 18 (2), 419–428. doi:10.1007/s40200-019-00435-w
Bigger, J. T., Albrecht, P., Steinman, R. C., Rolnitzky, L. M., Fleiss, J. L., and Cohen, R. J. (1989). Comparison of time- and frequency domain-based measures of cardiac parasympathetic activity in Holter recordings after myocardial infarction. Am. J. Cardiol. 64 (8), 536–538. doi:10.1016/0002-9149(89)90436-0
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 (6), 992–996. doi:10.1016/S1388-2457(01)00507-7
Buysse, D. J., Reynolds, C. F., Monk, T. H., Berman, S. R., and Kupfer, D. J. (1989). The Pittsburgh sleep quality index: A new instrument for psychiatric practice and research. Psychiatry Res. 28 (2), 193–213. doi:10.1016/0165-1781(89)90047-4
Cerutti, S., Esposti, F., Ferrario, M., Sassi, R., and Signorini, M. G. (2007). Long-term invariant parameters obtained from 24-h holter recordings: A comparison between different analysis techniques. Chaos An Interdiscip. J. Nonlinear Sci. 17 (1), 015108. doi:10.1063/1.2437155
Chen, W., Zhuang, J., Yu, W., and Wang, Z. (2009). Measuring complexity using FuzzyEn, ApEn, and SampEn. Med. Eng. Phys. 31 (1), 61–68. doi:10.1016/j.medengphy.2008.04.005
Cooke, W. H., Hoag, J. B., Crossman, A. A., Kuusela, T. A., Tahvanainen, K. U. O., and Eckberg, D. L. (1999). Human responses to upright tilt: A window on central autonomic integration. J. Physiology 517 (2), 617–628. doi:10.1111/j.1469-7793.1999.0617t.x
Costa, M., Goldberger, A. L., and Peng, C. K. (2005). Multiscale entropy analysis of biological signals. Phys. Rev. E 71 (2), 021906. doi:10.1103/PhysRevE.71.021906
Costa, M., Goldberger, A. L., and Peng, C. K. (2002). Multiscale entropy analysis of complex physiologic time series. Phys. Rev. Lett. 89 (6), 068102. doi:10.1103/PhysRevLett.89.068102
Cui, X., Tian, L., Li, Z., Ren, Z., Zha, K., Wei, X., et al. (2020). On the variability of heart rate variability—evidence from prospective study of healthy young college students. Entropy 22 (11), 1302. doi:10.3390/e22111302
Cuschieri, S. (2019). The STROBE guidelines. Saudi J. Anaesth. 13 (1), S31–s34. doi:10.4103/sja.SJA_543_18
Cysarz, D., Bettermann, H., and Leeuwen, P. v. (2000). Entropies of short binary sequences in heart period dynamics. Am. J. Physiology-Heart Circulatory Physiology 278 (6), H2163–H2172. doi:10.1152/ajpheart.2000.278.6.H2163
de Moura-Tonello, S. C. G., Porta, A., Marchi, A., de Almeida Fagundes, A., Francisco, C. d. O., Rehder-Santos, P., et al. (2016). Cardiovascular variability analysis and baroreflex estimation in patients with type 2 diabetes in absence of any manifest neuropathy. PLOS ONE 11 (3), e0148903. doi:10.1371/journal.pone.0148903
de Zambotti, M., Trinder, J., Silvani, A., Colrain, I. M., and Baker, F. C. (2018). Dynamic coupling between the central and autonomic nervous systems during sleep: A review. Neurosci. Biobehav. Rev. 90, 84–103. doi:10.1016/j.neubiorev.2018.03.027
Delgado-Bonal, A., and Marshak, A. (2019). Approximate entropy and sample entropy: A comprehensive tutorial. Entropy 21 (6), 541. doi:10.3390/e21060541
Dohare, A. K., Kumar, V., and Kumar, R. (2014). An efficient new method for the detection of QRS in electrocardiogram. Comput. Electr. Eng. 40 (5), 1717–1730. doi:10.1016/j.compeleceng.2013.11.004
Faust, O., Acharya, U. R., Molinari, F., Chattopadhyay, S., and Tamura, T. (2012). Linear and non-linear analysis of cardiac health in diabetic subjects. Biomed. Signal Process. Control 7 (3), 295–302. doi:10.1016/j.bspc.2011.06.002
Fei, L., Copie, X., Malik, M., and Camm, A. J. (1996). Short- and long-term assessment of heart rate variability for risk stratification after acute myocardial infarction. Am. J. Cardiol. 77 (9), 681–684. doi:10.1016/S0002-9149(97)89199-0
Feng, N., Yang, J., Xu, H., Zhang, C., Wang, F., Wu, X., et al. (2021). The associations between sleep architecture and metabolic parameters in patients with obstructive sleep apnea: A hospital-based cohort study. Front. Neurology 12, 606031. doi:10.3389/fneur.2021.606031
Francis, D. P., Willson, K., Georgiadou, P., Wensel, R., Davies, L. C., Coats, A., et al. (2002). Physiological basis of fractal complexity properties of heart rate variability in man. J. Physiology 542 (2), 619–629. doi:10.1113/jphysiol.2001.013389
Frattola, A., Parati, G., Gamba, P., Paleari, F., Mauri, G., Di Rienzo, M., et al. (1997). Time and frequency domain estimates of spontaneous baroreflex sensitivity provide early detection of autonomic dysfunction in diabetes mellitus. Diabetologia 40 (12), 1470–1475. doi:10.1007/s001250050851
Furlan, R., Guzzetti, S., Crivellaro, W., Dassi, S., Tinelli, M., Baselli, G., et al. (1990). Continuous 24-hour assessment of the neural regulation of systemic arterial pressure and RR variabilities in ambulant subjects. Circulation 81(2), 537–547. doi:10.1161/01.CIR.81.2.537
Garber, A. J., Handelsman, Y., Grunberger, G., Einhorn, D., Abrahamson, M. J., Barzilay, J. I., et al. (2020). Consensus statement by the American association of clinical endocrinologists and American college of endocrinology on the comprehensive type 2 diabetes management algorithm – 2020 executive summary. Endocr. Pract. 26 (1), 107–139. doi:10.4158/CS-2019-0472
Gerritsen, J., Dekker, J. M., TenVoorde, B. J., Kostense, P. J., Heine, R. J., Bouter, L. M., et al. (2001). Impaired autonomic function is associated with increased mortality, especially in subjects with diabetes, hypertension, or a history of cardiovascular disease: The hoorn study. Diabetes Care 24 (10), 1793–1798. doi:10.2337/diacare.24.10.1793
Glos, M., Fietze, I., Blau, A., Baumann, G., and Penzel, T. (2014). Cardiac autonomic modulation and sleepiness: Physiological consequences of sleep deprivation due to 40h of prolonged wakefulness. Physiology Behav. 125, 45–53. doi:10.1016/j.physbeh.2013.11.011
Goit, R. K., Khadka, R., Sharma, S. K., Limbu, N., and Paudel, B. H. (2012). Cardiovascular autonomic function and vibration perception threshold in type 2 diabetes mellitus. J. Diabetes its Complicat. 26 (4), 339–342. doi:10.1016/j.jdiacomp.2012.03.026
González, J. J., Cordero, J. J., Feria, M., and Pereda, E. (2000). Detection and sources of nonlinearity in the variability of cardiac R-R intervals and blood pressure in rats. Am. J. Physiology-Heart Circulatory Physiology 279 (6), H3040–H3046. doi:10.1152/ajpheart.2000.279.6.H3040
Hamilton, P. S., and Tompkins, W. J. (1986). Quantitative investigation of QRS detection rules using the MIT/BIH arrhythmia database. Ieee Trans. Biomed. Eng. 33 (12), 1157–1165. doi:10.1109/tbme.1986.325695
Hur, M. H., Lee, M. K., Seong, K., and Hong, J. H. (2020). Deterioration of sleep quality according to glycemic status. Diabetes & Metabolism J. 44 (5), 679–686. doi:10.4093/dmj.2019.0125
Inzucchi, S. E., Bergenstal, R. M., Buse, J. B., Diamant, M., Ferrannini, E., Nauck, M., et al. (2015). Management of hyperglycemia in type 2 diabetes, 2015: A patient-centered approach: Update to a position statement of the American diabetes association and the European association for the study of diabetes. Diabetes Care 38 (1), 140–149. doi:10.2337/dc14-2441
Jo, J. A., Blasi, A., Valladares, E., Juarez, R., Baydur, A., and Khoo, M. C. K. (2005). Determinants of heart rate variability in obstructive sleep apnea syndrome during wakefulness and sleep. Am. J. Physiology-Heart Circulatory Physiology 288 (3), H1103–H1112. doi:10.1152/ajpheart.01065.2003
Johnson, A. E. W., Behar, J., Andreotti, F., Clifford, G. D., and Oster, J. (2015). Multimodal heart beat detection using signal quality indices. Physiol. Meas. 36 (8), 1665–1677. doi:10.1088/0967-3334/36/8/1665
Jordan, A. S., McSharry, D. G., and Malhotra, A. (2014). Adult obstructive sleep apnoea. Lancet 383 (9918), 736–747. doi:10.1016/s0140-6736(13)60734-5
Kent, B. D., Grote, L., Ryan, S., Pépin, J. L., Bonsignore, M. R., Tkacova, R., et al. (2014). Diabetes mellitus prevalence and control in sleep-disordered breathing: The European sleep apnea cohort (ESADA) study. Chest 146 (4), 982–990. doi:10.1378/chest.13-2403
Kleiger, R. E., Stein, P. K., and Bigger, J. T. (2005). Heart rate variability: Measurement and clinical utility. Ann. Noninvasive Electrocardiol. 10 (1), 88–101. doi:10.1111/j.1542-474X.2005.10101.x
Koren, D., Levitt Katz, L. E., Brar, P. C., Gallagher, P. R., Berkowitz, R. I., and Brooks, L. J. (2011). Sleep architecture and glucose and insulin homeostasis in obese adolescents. Diabetes Care 34 (11), 2442–2447. doi:10.2337/dc11-1093
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
Leproult, R., Copinschi, G., Buxton, O., and Van Cauter, E. (1997). Sleep loss results in an elevation of cortisol levels the next evening. Sleep 20 (10), 865–870. doi:10.1093/sleep/20.10.865
Li, F., and Yang, H. (2013). A real-time QRS complex detection algorithm based on composite filter and curve feature point extraction. Comput. Appl. Softw. 30 (11).
Liu, S.-H., Lo, L.-W., Tsai, T.-Y., Cheng, W.-H., Lin, Y.-J., Chang, S.-L., et al. (2020). Circadian rhythm dynamics on multiscale entropy identifies autonomic dysfunction associated with risk of ventricular arrhythmias and near syncope in chronic kidney disease. J. Cardiol. 76 (6), 542–548. doi:10.1016/j.jjcc.2020.05.017
Lombardi, F., and Stein, P. K. (2011). Origin of heart rate variability and turbulence: An appraisal of autonomic modulation of cardiovascular function. Front. Physiology 2, 95. doi:10.3389/fphys.2011.00095
López-Cano, C., Gutiérrez-Carrasquilla, L., Sánchez, E., González, J., Yeramian, A., Martí, R., et al. (2019). Sympathetic hyperactivity and sleep disorders in individuals with type 2 diabetes. Front. Endocrinol. 10, 752. doi:10.3389/fendo.2019.00752
Mador, M. J., and Tobin, M. J. (1991). Effect of alterations in mental activity on the breathing pattern in healthy subjects. Am. Rev. Respir. Dis. 144 (3), 481–487. doi:10.1164/ajrccm/144.3_Pt_1.481
Magagnin, V., Bassani, T., Bari, V., Turiel, M., Maestri, R., Pinna, G. D., et al. (2011). Non-stationarities significantly distort short-term spectral, symbolic and entropy heart rate variability indices. Physiol. Meas. 32 (11), 1775–1786. doi:10.1088/0967-3334/32/11/S05
Malliani, A., Pagani, M., Lombardi, F., and Cerutti, S. (1991). Cardiovascular neural regulation explored in the frequency domain. Circulation 84(2), 482–492. doi:10.1161/01.CIR.84.2.482
Marcheva, B., Ramsey, K. M., Buhr, E. D., Kobayashi, Y., Su, H., Ko, C. H., et al. (2010). Disruption of the clock components CLOCK and BMAL1 leads to hypoinsulinaemia and diabetes. Nature 466 (7306), 627–631. doi:10.1038/nature09253
Markwald, R. R., Melanson, E. L., Smith, M. R., Higgins, J., Perreault, L., Eckel, R. H., et al. (2013). Impact of insufficient sleep on total daily energy expenditure, food intake, and weight gain. Proc. Natl. Acad. Sci. 110(14), 5695–5700. doi:10.1073/pnas.1216951110
Markwald, R. R., and Wright, K. P. (2012). “Circadian misalignment and sleep disruption in shift work: Implications for fatigue and risk of weight gain and obesity,” in Sleep loss and obesity: Intersecting epidemics. Editors P. Shiromani, T. Horvath, S. Redline, and E. Van Cauter (New York, NY: Springer New York), 101–118.
Martelli, D., Silvani, A., McAllen, R. M., May, C. N., and Ramchandra, R. (2014). The low frequency power of heart rate variability is neither a measure of cardiac sympathetic tone nor of baroreflex sensitivity. Am. J. Physiology-Heart Circulatory Physiology 307 (7), H1005–H1012. doi:10.1152/ajpheart.00361.2014
Martyn-Nemeth, P., Phillips, S. A., Mihailescu, D., Farabi, S. S., Park, C., Lipton, R., et al. (2018). Poor sleep quality is associated with nocturnal glycaemic variability and fear of hypoglycaemia in adults with type 1 diabetes. J. Adv. Nurs. 74 (10), 2373–2380. doi:10.1111/jan.13765
Maser, R. E., Mitchell, B. D., Vinik, A. I., and Freeman, R. (2003). The association between cardiovascular autonomic neuropathy and mortality in individuals with diabetes: A meta-analysis. Diabetes Care 26 (6), 1895–1901. doi:10.2337/diacare.26.6.1895
Montano, N., Porta, A., Cogliati, C., Costantino, G., Tobaldini, E., Casali, K. R., et al. (2009). Heart rate variability explored in the frequency domain: A tool to investigate the link between heart and behavior. Neurosci. Biobehav. Rev. 33 (2), 71–80. doi:10.1016/j.neubiorev.2008.07.006
Mosenzon, O., Pollack, R., and Raz, I. (2016). Treatment of type 2 diabetes: From "guidelines" to "position statements" and back: Recommendations of the Israel national diabetes council. Diabetes Care 39 (2), S146–S153. doi:10.2337/dcS15-3003
National Center for Cardiovascular Diseases, Chinese medical doctor association, hypertension committee of Chinese medical doctor association, Chinese society of Cardiology, Chinese medical association,, and hypertension committee of cross-straits medical exchanges association (2022). clinical practice guidelines for the management of hypertension in China. Chin. J. Cardiol. 50, 1050–1095. doi:10.3760/cma.j.cn112148-20220809-00613
Nguyen, J., and Wright, K. P. (2010). Influence of weeks of circadian misalignment on leptin levels. Nat. Sci. Sleep. 2, 9–18. doi:10.2147/nss.s7624
Nikulin, V. V., and Brismar, T. (2004). Comment on "Multiscale entropy analysis of complex physiologic time series". Phys. Rev. Lett. 92 (8), 089803. author reply 089804. doi:10.1103/PhysRevLett.92.089803
Nolan, J., Batin, P. D., Andrews, R., Lindsay, S. J., Brooksby, P., Mullen, M., et al. (1998). Prospective study of heart rate variability and mortality in chronic heart failure: Results of the United Kingdom heart failure evaluation and assessment of risk trial (UK-heart). Circulation 98(15), 1510–1516. doi:10.1161/01.CIR.98.15.1510
Nunan, D., Sandercock, G. R. H., and Brodie, D. A. (2010). A quantitative systematic review of normal values for short-term heart rate variability in healthy adults. Pacing Clin. Electrophysiol. 33 (11), 1407–1417. doi:10.1111/j.1540-8159.2010.02841.x
Pagani, M., Lombardi, F., Guzzetti, S., Rimoldi, O., Furlan, R., Pizzinelli, P., et al. (1986). Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog. Circulation Res. 59(2), 178–193. doi:10.1161/01.RES.59.2.178
Pan, J., and Tompkins, W. J. (1985). A real-time QRS detection algorithm. Ieee Trans. Biomed. Eng. 32 (3), 230–236. doi:10.1109/tbme.1985.325532
Peng, C. K., Buldyrev, S. V., Havlin, S., Simons, M., Stanley, H. E., and Goldberger, A. L. (1994). Mosaic organization of DNA nucleotides. Phys. Rev. E 49 (2), 1685–1689. doi:10.1103/PhysRevE.49.1685
Peng, C. K., Havlin, S., Stanley, H. E., and Goldberger, A. L. (1995). Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series. Chaos An Interdiscip. J. Nonlinear Sci. 5 (1), 82–87. doi:10.1063/1.166141
Perkiömäki, J. (2011). Heart rate variability and non-linear dynamics in risk stratification. Front. Physiology 2, 81. doi:10.3389/fphys.2011.00081
Pfeifer, M. A., Cook, D., Brodsky, J., Tice, D., Reenan, A., Swedine, S., et al. (1982). Quantitative evaluation of cardiac parasympathetic activity in normal and diabetic man. Diabetes 31 (4), 339–345. doi:10.2337/diab.31.4.339
Pincus, S. M. (1991). Approximate entropy as a measure of system complexity. Proc. Natl. Acad. Sci. 88(6), 2297–2301. doi:10.1073/pnas.88.6.2297
Pincus, S. M., Cummins, T. R., and Haddad, G. G. (1993). Heart rate control in normal and aborted-SIDS infants. Am. J. Physiology-Regulatory, Integr. Comp. Physiology 264 (3), R638–R646. doi:10.1152/ajpregu.1993.264.3.R638
Pincus, S. M., and Huang, W.-M. (1992). Approximate entropy: Statistical properties and applications. Commun. Statistics - Theory Methods 21 (11), 3061–3077. doi:10.1080/03610929208830963
Pop-Busui, R. (2010). Cardiac autonomic neuropathy in diabetes: A clinical perspective. Diabetes Care 33 (2), 434–441. doi:10.2337/dc09-1294
Porta, A., Baselli, G., Guzzetti, S., Pagani, M., Malliani, A., and Cerutti, S. (2000a). Prediction of short cardiovascular variability signals based on conditional distribution. IEEE Trans. Biomed. Eng. 47 (12), 1555–1564. doi:10.1109/10.887936
Porta, A., Gnecchi-Ruscone, T., Tobaldini, E., Guzzetti, S., Furlan, R., and Montano, N. (2007a). Progressive decrease of heart period variability entropy-based complexity during graded head-up tilt. J. Appl. Physiology 103 (4), 1143–1149. doi:10.1152/japplphysiol.00293.2007
Porta, A., Guzzetti, S., Furlan, R., Gnecchi-Ruscone, T., Montano, N., and Malliani, A. (2007b). Complexity and nonlinearity in short-term heart period variability: Comparison of methods based on local nonlinear prediction. IEEE Trans. Biomed. Eng. 54 (1), 94–106. doi:10.1109/TBME.2006.883789
Porta, A., Guzzetti, S., Montano, N., Furlan, R., Pagani, M., Malliani, A., et al. (2001). Entropy, entropy rate, and pattern classification as tools to typify complexity in short heart period variability series. IEEE Trans. Biomed. Eng. 48 (11), 1282–1291. doi:10.1109/10.959324
Porta, A., Guzzetti, S., Montano, N., Pagani, M., Somers, V., Malliani, A., et al. (2000b). Information domain analysis of cardiovascular variability signals: Evaluation of regularity, synchronisation and co-ordination. Med. Biol. Eng. Comput. 38 (2), 180–188. doi:10.1007/bf02344774
Porta, A., Maria, B. D., Bari, V., Marchi, A., and Faes, L. (2017). Are nonlinear model-free conditional entropy approaches for the assessment of cardiac control complexity superior to the linear model-based one? IEEE Trans. Biomed. Eng. 64 (6), 1287–1296. doi:10.1109/TBME.2016.2600160
Punjabi, N. M., and Polotsky, V. Y. (2005). Disorders of glucose metabolism in sleep apnea. J. Appl. Physiology 99 (5), 1998–2007. doi:10.1152/japplphysiol.00695.2005
Reyes del Paso, G. A., Langewitz, W., Mulder, L. J. M., van Roon, A., and Duschek, S. (2013). The utility of low frequency heart rate variability as an index of sympathetic cardiac tone: A review with emphasis on a reanalysis of previous studies. Psychophysiology 50 (5), 477–487. doi:10.1111/psyp.12027
Richman, J. S., and Moorman, J. R. (2000). Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiology-Heart Circulatory Physiology 278 (6), H2039–H2049. doi:10.1152/ajpheart.2000.278.6.H2039
Ruiz, J., Monbaron, D., Parati, G., Perret, S., Haesler, E., Danzeisen, C., et al. (2005). Diabetic neuropathy is a more important determinant of baroreflex sensitivity than carotid elasticity in type 2 diabetes. Hypertension 46(1), 162–167. doi:10.1161/01.HYP.0000169053.14440.7d
Sako, T., Burioka, N., Suyama, H., Nomura, T., Takeshima, T., and Shimizu, E. (2001). Nonlinear behavior of human respiratory movement during different sleep stages. Chronobiology Int. 18 (1), 71–83. doi:10.1081/CBI-100001172
Sassi, R., Cerutti, S., Lombardi, F., Malik, M., Huikuri, H. V., Peng, C.-K., et al. (2015). Advances in heart rate variability signal analysis: Joint position statement by the e-cardiology ESC working group and the European heart rhythm association co-endorsed by the Asia pacific heart rhythm society. EP Eur. 17 (9), 1341–1353. doi:10.1093/europace/euv015
Saul, J. P., Albrecht, P., Berger, R. D., and Cohen, R. J. (1988). Analysis of long term heart rate variability: Methods, 1/f scaling and implications. Comput. Cardiol. 14, 419–422.
Scheer, F. A. J. L., Hilton, M. F., Mantzoros, C. S., and Shea, S. A. (2009). Adverse metabolic and cardiovascular consequences of circadian misalignment. Proc. Natl. Acad. Sci. 106(11), 4453–4458. doi:10.1073/pnas.0808180106
Scholz, U. J., Bianchi, A. M., Cerutti, S., and Kubicki, S. (1997). Vegetative background of sleep: Spectral analysis of the heart rate variability. Physiology Behav. 62 (5), 1037–1043. doi:10.1016/S0031-9384(97)00234-5
Schramm, P., Neville, A. G., Madison, S., Thomas, R. J., and Baker, D. N. (2009). Cardio-pulmonary coupling (CPC) measures correlate to standard sleep variables in a random clinical sample of patients suspected with sleep disordered breathing (SDB). Sleep 32, A189.
Shaffer, F., McCraty, R., and Zerr, C. L. (2014). A healthy heart is not a metronome: An integrative review of the heart's anatomy and heart rate variability. Front. Psychol. 5, 1040. doi:10.3389/fpsyg.2014.01040
Silvani, A., and Dampney, R. A. L. (2013). Central control of cardiovascular function during sleep. Am. J. Physiology-Heart Circulatory Physiology 305 (12), H1683–H1692. doi:10.1152/ajpheart.00554.2013
Singh, J. P., Larson, M. G., O'Donnell, C. J., Wilson, P. F., Tsuji, H., Lloyd-Jones, D. M., et al. (2000). Association of hyperglycemia with reduced heart rate variability (The Framingham Heart Study). Am. J. Cardiol. 86 (3), 309–312. doi:10.1016/s0002-9149(00)00920-6
Solnik, S., DeVita, P., Rider, P., Long, B., and Hortobágyi, T. (2008). Teager-Kaiser Operator improves the accuracy of EMG onset detection independent of signal-to-noise ratio. Acta Bioeng. Biomech. 10 (2), 65–68.
Somers, V. K., Dyken, M. E., Mark, A. L., and Abboud, F. M. (1993). Sympathetic-nerve activity during sleep in normal subjects. N. Engl. J. Med. 328 (5), 303–307. doi:10.1056/nejm199302043280502
Spiegel, K., Knutson, K., Leproult, R., Tasali, E., and Cauter, E. V. (2005). Sleep loss: A novel risk factor for insulin resistance and type 2 diabetes. J. Appl. Physiology 99 (5), 2008–2019. doi:10.1152/japplphysiol.00660.2005
Spiegel, K., Leproult, R., and Van Cauter, E. (1999). Impact of sleep debt on metabolic and endocrine function. Lancet 354 (9188), 1435–1439. doi:10.1016/S0140-6736(99)01376-8
Spiegel, K., Tasali, E., Leproult, R., Scherberg, N., and Van Cauter, E. (2011). Twenty-four-Hour profiles of acylated and total ghrelin: Relationship with glucose levels and impact of time of day and sleep. J. Clin. Endocrinol. Metabolism 96 (2), 486–493. doi:10.1210/jc.2010-1978
Spira, A. P., Beaudreau, S. A., Stone, K. L., Kezirian, E. J., Lui, L.-Y., Redline, S., et al. (2011). Reliability and validity of the Pittsburgh sleep quality index and the epworth sleepiness scale in older men. Journals Gerontology Ser. A 67A (4), 433–439. doi:10.1093/gerona/glr172
Stein, P. K., and Pu, Y. (2012). Heart rate variability, sleep and sleep disorders. Sleep. Med. Rev. 16 (1), 47–66. doi:10.1016/j.smrv.2011.02.005
Tank, J., Jordan, J., Neuke, A., Mölle, A., and Weck, M. (2001). Spontaneous baroreflex sensitivity and heart rate variability are not superior to classic autonomic testing in older patients with type 2 diabetes. Am. J. Med. Sci. 322 (1), 24–30. doi:10.1097/00000441-200107000-00005
Tasali, E., Leproult, R., Ehrmann, D. A., and Van Cauter, E. (2008). Slow-wave sleep and the risk of type 2 diabetes in humans. Proc. Natl. Acad. Sci. 105(3), 1044–1049. doi:10.1073/pnas.0706446105
Task Force of the European Society of Cardiology the North American Society of Pacing Electrophysiology (1996). Heart rate variability. Circulation 93(5), 1043–1065. doi:10.1161/01.CIR.93.5.1043
Thomas, R. J., Mietus, J. E., Peng, C.-K., and Goldberger, A. L. (2005). An electrocardiogram-based technique to assess cardiopulmonary coupling during sleep. Sleep 28 (9), 1151–1161. doi:10.1093/sleep/28.9.1151
Thomas, R. J., Mietus, J. E., Peng, C.-K., Goldberger, A. L., Crofford, L. J., and Chervin, R. D. (2010). Impaired sleep quality in fibromyalgia: Detection and quantification with ECG-based cardiopulmonary coupling spectrograms. Sleep. Med. 11 (5), 497–498. doi:10.1016/j.sleep.2009.09.003
Tobaldini, E., Costantino, G., Solbiati, M., Cogliati, C., Kara, T., Nobili, L., et al. (2017). Sleep, sleep deprivation, autonomic nervous system and cardiovascular diseases. Neurosci. Biobehav. Rev. 74, 321–329. doi:10.1016/j.neubiorev.2016.07.004
Tobaldini, E., Pecis, M., and Montano, N. (2014). Effects of acute and chronic sleep deprivation on cardiovascular regulation. Arch. Ital. De. Biol. 152 (2-3), 103–110. doi:10.12871/000298292014235
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 (4), 253–264. doi:10.1046/j.1365-2869.2001.00263.x
Turek, F. W., Joshu, C., Kohsaka, A., Lin, E., Ivanova, G., McDearmon, E., et al. (2005). Obesity and metabolic syndrome in circadian clock mutant mice. Science 308(5724), 1043–1045. doi:10.1126/science.1108750
Valencia, J. F., Porta, A., Vallverdu, M., Claria, F., Baranowski, R., Orlowska-Baranowska, E., et al. (2009). Refined multiscale entropy: Application to 24-h holter recordings of heart period variability in healthy and aortic stenosis subjects. IEEE Trans. Biomed. Eng. 56 (9), 2202–2213. doi:10.1109/TBME.2009.2021986
Van Cauter, E., Spiegel, K., Tasali, E., and Leproult, R. (2008). Metabolic consequences of sleep and sleep loss. Sleep. Med. 9, S23–S28. doi:10.1016/S1389-9457(08)70013-3
VanCauter, E., Polonsky, K. S., and Scheen, A. J. (1997). Roles of circadian rhythmicity and sleep in human glucose regulation. Endocr. Rev. 18 (5), 716–738. doi:10.1210/er.18.5.716
Vinik, A. I., Maser, R. E., Mitchell, B. D., and Freeman, R. (2003). Diabetic autonomic neuropathy. Diabetes Care 26 (5), 1553–1579. doi:10.2337/diacare.26.5.1553
Vinik, A. I., and Ziegler, D. (2007). Diabetic cardiovascular autonomic neuropathy. Circulation 115(3), 387–397. doi:10.1161/CIRCULATIONAHA.106.634949
Voss, A., Kurths, J., Kleiner, H. J., Witt, A., Wessel, N., Saparin, P., et al. (1996). The application of methods of non-linear dynamics for the improved and predictive recognition of patients threatened by sudden cardiac death. Cardiovasc Res. 31 (3), 419–433. doi:10.1016/0008-6363(96)00008-9
Wessel, N., Ziehmann, C., Kurths, J., Meyerfeldt, U., Schirdewan, A., and Voss, A. (2000). Short-term forecasting of life-threatening cardiac arrhythmias based on symbolic dynamics and finite-time growth rates. Phys. Rev. E 61 (1), 733–739. doi:10.1103/PhysRevE.61.733
Yamada, T., and Katagiri, H. (2007). Avenues of communication between the brain and tissues/organs involved in energy homeostasis. Endocr. J. 54 (4), 497–505. doi:10.1507/endocrj.KR-106
Ye, S., Ruan, P., Yong, J., Shen, H., Liao, Z., and Dong, X. (2016). The impact of the HbA1c level of type 2 diabetics on the structure of haemoglobin. Sci. Rep. 6 (1), 33352. doi:10.1038/srep33352
Yeh, G. Y., Mietus, J. E., Peng, C.-K., Phillips, R. S., Davis, R. B., Wayne, P. M., et al. (2008). Enhancement of sleep stability with Tai Chi exercise in chronic heart failure: Preliminary findings using an ECG-based spectrogram method. Sleep. Med. 9 (5), 527–536. doi:10.1016/j.sleep.2007.06.003
Yu, Y., Hu, L., Xu, Y., Wu, S., Chen, Y., Zou, W., et al. (2020). Impact of blood glucose control on sympathetic and vagus nerve functional status in patients with type 2 diabetes mellitus. Acta Diabetol. 57 (2), 141–150. doi:10.1007/s00592-019-01393-8
Zhang, L., Wu, H., Zhang, X., Wei, X., Hou, F., and Ma, Y. (2020). Sleep heart rate variability assists the automatic prediction of long-term cardiovascular outcomes. Sleep. Med. 67, 217–224. doi:10.1016/j.sleep.2019.11.1259
Ziegler, D., Zentai, C. P., Perz, S., Rathmann, W., Haastert, B., Döring, A., et al. (2008). Prediction of mortality using measures of cardiac autonomic dysfunction in the diabetic and nondiabetic population: The MONICA/KORA augsburg cohort study. Diabetes Care 31 (3), 556–561. doi:10.2337/dc07-1615
Keywords: T2DM, heart rate variability (HRV), sleep, autonomic nervous system (ANS), metabolic function, nonlinear analysis, glycemic control
Citation: Cheng W, Chen H, Tian L, Ma Z and Cui X (2023) Heart rate variability in different sleep stages is associated with metabolic function and glycemic control in type 2 diabetes mellitus. Front. Physiol. 14:1157270. doi: 10.3389/fphys.2023.1157270
Received: 02 February 2023; Accepted: 24 March 2023;
Published: 14 April 2023.
Edited by:
Vitor Engracia Valenti, São Paulo State University, BrazilReviewed by:
Alberto Porta, University of Milan, ItalySzymon Siecinski, Silesian University of Technology, Poland
Copyright © 2023 Cheng, Chen, Tian, Ma and Cui. 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.
*Correspondence: Zhimin Ma, mazhimin01@sina.com; Xingran Cui, cuixr@seu.edu.cn
†These authors have contributed equally to this work and share first authorship