Skip to main content
Download PDF

Association between circadian variation of heart rate and mortality among critically ill patients: a retrospective cohort study

BMC Anesthesiology volume 22, Article number: 45 (2022) Cite this article

Abstract

Background

Heart rate (HR) related parameters, such as HR variability, HR turbulence, resting HR, and nighttime mean HR have been recognized as independent predictors of mortality. However, the influence of circadian changes in HR on mortality remains unclear in intensive care units (ICU). The study is designed to evaluate the relationship between the circadian variation in HR and mortality risk among critically ill patients.

Methods

The present study included 4,760 patients extracted from the Multiparameter Intelligent Monitoring in Intensive Care II database. The nighttime mean HR/daytime mean HR ratio was adopted as the circadian variation in HR. According to the median value of the circadian variation in HR, participants were divided into two groups: group A (≤ 1) and group B (> 1). The outcomes included ICU, hospital, 30-day, and 1-year mortalities. The prognostic value of HR circadian variation was investigated by multivariable logistic regression models and Cox proportional hazards models.

Results

Patients in group B (n = 2,471) had higher mortality than those in group A (n = 2,289). Multivariable models revealed that the higher circadian variation in HR was associated with ICU mortality (odds ratio [OR], 1.393; 95% confidence interval [CI], 1.112–1.745; P = 0.004), hospital mortality (OR, 1.393; 95% CI, 1.112–1.745; P = 0.004), 30-day mortality (hazard ratio, 1.260; 95% CI, 1.064–1.491; P = 0.007), and 1-year mortality (hazard ratio, 1.207; 95% CI, 1.057–1.378; P = 0.005), especially in patients with higher SOFA scores.

Conclusions

The circadian variation in HR might aid in the early identification of critically ill patients at high risk of associated with ICU, hospital, 30-day, and 1-year mortalities.

Peer Review reports

Background

Heart rate (HR) is regulated by the autonomic nervous system, and therefore, the autonomic nervous system function could be reflected by HR variations, such as heart rate variability (HRV), heart rate turbulence (HRT), resting heart rate (RHR), and nighttime mean heart rate (NHR). Among them, HRV is measured by the variation in the beat-to-beat interval between consecutive heart beats to evaluate the autonomic nervous mediation of the cardiovascular system [1, 2]. Several studies have reported that the decrease in HRV was associated with impaired sympathetic modulation, which led to early in-hospital deterioration of sepsis and multiple organ dysfunction syndrome (MODS), and a higher mortality in the intensive care unit (ICU) [3,4,5,6,7,8,9,10]. However, Monfredi et al. have put forward a novel controversial discussion about the exponential decay relationship between HR and HRV, that is, HRV is primarily dependent on HR, with HRV increasing when HR slows, and vice versa [11]. Therefore, HRV cannot be used in any simple way to assess the autonomic nerve activity to the heart. HRT describes the physiological short-term oscillation of beat-to-beat intervals after spontaneous ventricular premature beats [12], which has been introduced as an autonomic predictor for cardiac mortality following myocardial infarction or heart failure [12, 13]. Besides, there have been previous reports in literature that increased resting HR can be associated with increased short-term and long-term mortality from critical illness [14,15,16,17,18,19]. Although the above parameters have been verified to have a higher predictive value for mortality, they either have their own limitations or ignore the circadian rhythm of the HR.

Most of physiological parameters, such as HR, blood pressure, spontaneous motor activity, cortisol, sleep cycles, the immune system, the autonomic nervous system, and body temperature follow 24-h stable circadian rhythms in healthy individuals [20]. However, these circadian rhythms may be interrupted or impaired by the loss of external zeitgebers when critically ill patients are often exposed to different degrees of artificial light, noise, continuous therapeutic intervention including ventilation, parenteral nutrition, and medications in the ICU [20,21,22,23]. Therefore, we took the circadian rhythms of HR into account to study the relationship between circadian variation of HR and mortality among critically ill patients admitted to the ICU in the Multiparameter Intelligent Monitoring in Intensive Care II (MIMIC-II) database.

Methods

Study design

This was an observational study conducted using individual medical information from the MIMIC-II database, which is maintained by the Laboratory for Computational Physiology at the Massachusetts Institute of Technology (MIT; Cambridge, MA) and contains de-identified data on patients hospitalized at an ICU at Beth Israel Deaconess Medical Center (Boston, MA) from 2001 to 2008 [24]. All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional research committee (Cardif Metropolitan University ethics committee—17/4/02R) and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. Since the study did not impact clinical care and all protected health information was de-identified, the requirement for individual patient consent was waived by the Institutional Review Board of Beth Israel Deaconess Medical Center. All authors have attended the “protecting human subjects training” and obtained the certificate which permits access to the data.

Date extraction

Patient data were exacted from MIMIC-II (version 2.6) using SAS version 9.4 (SAS Institute, Cary, NC). The extracted data, including patient identifiers, demographics, clinical diagnoses, medications effecting HR, scoring systems, length of stay in ICU, and complete HR records within the first 24 h after ICU admission, were acquired from 2001 to 2008 at Beth Israel Deaconess Medical Center. Age, sex, and ethnicity were included in the demographics and clinical diseases were diagnosed according to the International Classification of Diseases-9 diagnoses (ICD, 9th Editor), which includes respiratory failure, renal failure, liver cirrhosis, shock, diabetes uncomplicated, diabetes complicated, acquired immune deficiency syndrome (AIDS), lymphoma, metastatic cancer, coagulopathy, rheumatoid arthritis, infection, poisoning, hypoferric anemia, and a series of cardiovascular diseases. Severity-of-illness scores systems, including the Simplified Acute Physiology Score-I (SAPS-I) score and Sequential Organ Failure Assessment (SOFA) score, were calculated using physiological measurements and clinical information according to published recommendations and accepted formulae [25]. The medications included catecholamine, β-blockers, sedatives, and opioid analgesics.

Patient selection

The data of a total of 4,760 adult patients were extracted from the MIMIC-II V2.6 database. These patients met the following characteristics: (1) single ICU admission; (2) the length of ICU stay ≥ 24 h; (3) no history of cardiovascular diseases, including cardiac arrhythmias, cardiac arrest, hypertension, heart failure, and coronary heart disease, or no use of cardiac pacemaker due to any cause; (4) complete HR and covariate records; (5) no outliers in HR records. Outliers were defined as values exceeding the mean ± 3 times the SD [26].

Circadian rhythm of HR

HR was measured using a bedside monitor (Component Monitoring System IntelliVue MP-70; Philips Healthcare, Andover, MA) during the first 24 h in the ICU. In the present study, we focused on differences between the mean nighttime and daytime HR. Thus, the circadian rhythm of HR is calculated as the mean nighttime HR value divided by the mean daytime HR value. Patients were divided into two groups according to the median of the circadian rhythm of HR: the lower circadian rhythm of HR group (group A ≤ 1) and the higher circadian rhythm of HR group (group B > 1). The mean daytime HR value was computed as the average of all HR values from 7 am to 11 pm, and the mean nighttime HR was the average of all HR values from 11 pm to 7 am on the next day.

Endpoints

The primary outcome measures in our study were ICU mortality and hospital mortality. The secondary outcome measures were 30-day mortality and 1-year mortality.

Statistical analysis

Continuous variables are presented as the mean with SD or median with interquartile ranges. The one-way analysis of variance (ANOVA), or Kruskal–Wallis test was used as appropriate. Categorical variables are presented as frequency with percentage and compared using the Chi-Square test. To assess the association of circadian variation of HR with the study outcomes, we used univariate and multivariable logistic regression models for ICU and hospital mortality, cox proportional hazards model for the 30-day and 1-year mortality. All variables in the Cox proportional hazards model must meet the proportional hazard assumption. The multivariable models were adjusted for the following factors: (1) unadjusted in Model 1; (2) age, sex, and ethnicity in Model 2; (3) Model 2 plus clinical diagnoses (including respiratory failure, renal failure, liver cirrhosis, shock, diabetes uncomplicated, diabetes complicated, AIDS, lymphoma, metastatic cancer, coagulopathy, rheumatoid arthritis, infection, poisoning, and hypoferric anemia) in Model 3; (4) Model 3 plus medication taking on day 1 (including sedatives, catecholamines, β-blockers, and opioid analgesics) in Model 4; (5) Model 4 plus the first SOFA and SAPS-I score after ICU admission in Model 5; (6) Model 5 plus length of stay in ICU in Model 6; (7) Model 6 plus the first 24-h average HR in Model 7. Besides, based on the median of the first SOFA score, all patients were divided into the low SOFA score (0–5) group and the high SOFA score (6–24) group. Logistic regression and Cox regression were used to evaluate the relationship between the circadian rhythm of HR and outcomes in the two subgroups. P values of < 0.05 were considered to be statistically significant. All statistical analyses were performed with the SPSS software package version 22.0 (IBM Corp., Armonk, NY).

Results

Baseline characteristics

The records of 4760 patients were included in this study. Most of the study population came from surgical ICU (40.3%), followed by medical ICU (37.5%) (Supplementary Table 1). There were no significant differences in age, sex, ethnicity, and clinical diseases between the lower circadian rhythm of HR group (group A ≤ 1) and the higher circadian rhythm of HR group (group B > 1). However, patients in group B had higher SOFA scores [5.00 (2.00, 8.00) vs. 4.00 (2.00, 7.00), P = 0.001], SAPS-I scores [12.00 (9.00, 16.00) vs. 12.00 (8.00, 16.00), P = 0.006], and length of stay in ICU [2.52 (1.69, 4.96) vs. 2.31 (1.54, 4.26), P = 0.001) than those in group A. Moreover, patients in group B had a higher exposure rate for catecholamines (10.8% vs. 8.7%, P = 0.014) and opioid analgesics (2.5% vs. 1.3%, P = 0.002) than those in group A (Table 1). Besides, the univariate analysis of factors associated ICU, hospital, 30-day and 1-year mortality was presented in the Supplementary Table 2 and Supplementary Table 3. Furthermore, compared to ICU survival group, the dead patients during ICU stay were also more serious, with higher SOFA score [10.00 (7.00, 14.00) vs. 4.00 (2.00, 7.00), P < 0.001] and SAPS-I score [19.00 (15.00, 23.00) vs. 12.00 (8.00, 15.00), P < 0.001] (Table 2).

Table 1 Baseline Characteristics of the Study Subjects
Full size table
Table 2 Baseline Characteristics of the factors associated with ICU mortality and ICU survival
Full size table

Association of the circadian variation of HR with mortality

Patients in group B had higher ICU mortality (9.6% vs. 6.6%, P < 0.001), hospital mortality (12.3% vs. 8.7%; P < 0.001), 30-day mortality (13.9% vs. 10.3%; P < 0.001) and 1-year mortality (21.2% vs. 17.3%, P = 0.001) than that in group A patients (Table 1 and Fig. 1). These trends were also verified by the logistic regression analysis and the Cox regression model. Taking group A as the reference, after adjusting for all covariates, including age, sex, ethnicity, ICD-9 diagnoses, medication usage on day 1 of admission to the ICU, the first SOFA score, the first SAPS-I score, length of stay in the ICU, and the first 24-h average HR, the logistic regression results showed that group B, that is, the higher circadian variation of HR group was significantly associated with ICU mortality (odds ratio [OR], 1.393; 95% confidence interval [CI], 1.112–1.745; P = 0.004) and hospital mortality (OR, 1.393; 95% CI, 1.112–1.745; P = 0.004). Similarly, after adjusting for the same variables, the Cox proportional hazards model also showed that the higher circadian variation of HR group was significantly associated with 30-day mortality (hazard ratio, 1.260; 95% CI, 1.064–1.491; P = 0.007) and 1-year mortality (hazard ratio, 1.207; 95% CI, 1.057–1.378; P = 0.005) (Table 3).

Fig. 1
figure 1

The difference in ICU and hospital mortality in patients with different circadian variation of HR. Patients in group B had higher ICU (P < 0.001) and hospital mortality (P < 0.001) than patients in group A

Full size image
Table 3 Association of circadian variation of heart rate with Mortality in Overall Patients
Full size table

Interaction analysis

Considering the SOFA score to be a strong predictor of mortality in ICU patients, the interaction effect between the circadian variation of HR and the SOFA score was also analyzed. The results revealed no interaction between them (P > 0.05), which verified the absence of synergistic or antagonistic effect between the two in the prediction of mortality and showed that they were independent of each other in predicting mortality.

Sensitivity analysis

To diminish the influence of severity-of-illness on mortality from clinical practice, all individuals were also divided into two subgroups based on the median SOFA score of all the included patients: lower SOFA score group (0 ≤ SOFA score ≤ 5) and higher SOFA score group (6 ≤ SOFA score ≤ 24). Interestingly, in the higher SOFA score subgroup, taking group A as a reference, group B showed a remarkable association with short-term and long-term mortality in the multivariable regression model after adjusting a series of covariates. The specific statistical results were as follows: ICU mortality: OR, 1.444, 95% CI, 1.084–1.924, P = 0.012; hospital mortality: OR, 1.442, 95% CI, 1.110–1.873, P = 0.006; 30-day mortality: hazard ratio, 1.331, 95% CI, 1.091–1.624, P = 0.005; 1-year mortality: hazard ratio, 1.263, 95% CI, 1.069–1.493, P = 0.006. However, in the lower SOFA score subgroup, there was no statistical significance between the circadian variation of HR and mortality (Table 4). Kaplan–Meier survival curves also revealed that the short-term as well as long-term probability of survival was higher in group A than that in group B in the higher SOFA score subgroup (Fig. 2).

Table 4 Association of circadian variation of heart rate with Mortality in different SOFA score subgroups
Full size table
Fig. 2
figure 2

Kaplan–Meier survival analysis plot for mortality with different circadian variation of HR. Patients in group B had lower rates of 30-day (c) (P < 0.001) and 1-year (d) (P < 0.001) survival than patients in group A in the higher SOFA score subgroups. However, in the lower SOFA score subgroup, there was no statistical significance between the circadian variation of HR and 30-day (a) and 1-year (b) mortality

Full size image

Discussion

In this study, we investigated the relationship between circadian variation of HR and mortality in ICU patients. Our results show that patients with a higher circadian variation of HR had higher ICU mortality, hospital mortality, 30-day mortality and 1-year mortality than that in patients with a lower circadian variation of HR after adjusting for demographic parameters, clinical comorbidities, medications affecting heart rate, severity scoring systems including SAPS-I and SOFA score, length of stay in the ICU, and the first 24-h average HR.

A circadian pattern was observed in the variation of many biological and psychological processes, most of which comply a rhythmically change with a period length of approximately 24 h. Previous studies have reported the existence of circadian variations in cognitive functions [27, 28] as well as the underlying brain processes [29] in healthy individuals. Recently, there’s growing interest in the research of disturbed circadian rhythms in patients situated in special environments settings without clearly delineated day or night, such as ICUs or long-term care homes [23, 30, 31]. Patients from ICUs commonly suffer from impaired sleep quality, which may result from the unnatural light rhythm and light intensity, further leading to impaired immune responses and ability of recovery from illness [23, 32], ultimately causing different prognosis.

Paul and Lemmer showed that the 24-h profile parameters, including cortisol, blood pressure, HR, body temperature, and spontaneous motor activity, were all disturbed in 24 critically ill sedated patients compared with the normal rhythmic 24-h patterns in healthy controls [21]. However, clinical complications that occur when circadian rhythms are disturbed have been generally unappreciated, especially in ICU patients. Therefore, our team has been working on the study of the influence of circadian rhythms of physiological signals on the prognosis of ICU patients. Previous research has found that reduced circadian fluctuations of mean arterial pressure were associated with hospital mortality and ICU mortality among critically ill patients [33]. Besides, we also found that rise in the nocturnal mean arterial pressure within the first 24 h of admission might serve as an independent risk factor for short-term and long-term mortality; this may help in early risk stratification and personalized mean arterial pressure management in the ICU [34].

Previous studies have reported multiple HR variations. Among them, HRV, a temporal beat-to-beat variation in successive RR intervals on an electrocardiographic (ECG) recording, was considered as a reflection of the regulation of HR by the autonomic nervous system [35]. In recent years, the prognostic value of HRV was evaluated outside a specific cardiological context. The reduction in HRV was associated with a higher mortality in critically ill patients as reported in several studies [36]. More specifically, descending HRV was found as a common feature in septic patients. Especially, impaired sympathetic modulation, represented by low frequency power in spectrum analysis of HRV, could be the early marker of MODS for ICU patients with sepsis [8]. In addition, the low-/high-frequency power ratio in spectrum analysis of HRV, considered as a representation of sympathovagal imbalance, may predict early clinical deterioration [9], providing an early indication for the identification of impending septic shock among emergency department patients [10]. In the meantime, upon ICU admission, a short-term reduction of HRV was also reported as a predictive factor for the mortality in septic patients [37].

HRT is another parameter that describes the physiological short-term oscillation of beat-to-beat intervals after spontaneous ventricular premature beats [38]. HRT also could independently predict cardiac mortality in patients with diabetes mellitus or simple ischemic heart disease [12, 39]. Notably, a meta-analysis reveals that elevated resting HR, an indicator of myocardial oxygen consumption, may increase the mortality of acute coronary syndrome patients in the coronary intervention era [40]. Besides, a cohort study on the Japanese general population reports that nighttime mean HR is associated with increased risk of all-cause mortality in epidemiological settings [18]. Furthermore, by using pervasive sensing technology and artificial intelligence in ICU, the autonomous and granular monitoring system presented that the delirious patients in ICU had higher average heart rate [41]. These indices can be easily calculated from a 24-h Holter ECG. An increasing amount of evidence also supports the association of these parameters with adverse outcomes [42].

In this study, the circadian variation in HR was proposed as a powerful parameter for assessing HR change and emphasizes the circadian rhythm, which should be able to better predict the occurrence of adverse events combined with previous parameters in critically ill patients. In normal conditions, the circadian decline rate of HR is within 10%-20%, called the nocturnal dipping of HR, and is mostly determined by the endogenous neuroendocrine rhythm [18, 43]. Critically ill patients experience circumstance alterations in light–dark cycles, auditory disruption, iatrogenic treatment, and psychological reactions when they are admitted to the ICU. These external changes disturb the endogenous neuroendocrine rhythm and sleep–wake cycle, which may contribute to the nocturnal non-dipping or anti-dipping of HR. The nighttime higher HR or reduced dipping of sleep HR can represent the status of sympathetic overdrive. On the one hand, higher nighttime HR caused by sympathetic overdrive could be associated with insulin resistance or direct atherosclerotic lesions via hemodynamic disturbances [44, 45]. On the other hand, sleep HR was a better marker of mortality than awake or clinical HR [46]. Therefore, we proposed that circadian variations of HR (calculated as the ratio of NHR to daytime mean HR) reflected impaired autonomic nerve regulation. Abnormal circadian variation of HR might be associated with adverse prognosis.

The SOFA score is an acknowledged parameter that is robustly correlated with patient outcomes in the ICU setting [25]. Therefore, we also conducted a sensitivity analysis in which the patients were divided into two subgroups according to the median of SOFA scores. Surprisingly, the association between circadian variation of HR and mortality persisted in the higher SOFA score subgroup but not in the lower after adjusting for a series of variables. The regression analysis verified the association of circadian variation of HR with mortality again. Therefore, a higher circadian fluctuation of HR might serve as an independent risk factor for both short-term mortality and long-term mortality in critically ill patients, especially in the higher SOFA score individuals.

The main strength of our study is the high number of general ICU patients included. To our knowledge, this is the first report that circadian variation of HR is associated with short-term mortality and long-term mortality in ICU patients. The results reveal the importance of circadian rhythm in ICU patient management, which reminds the clinicians to attach great importance to the changes of the circadian rhythm parameters in ICU patients and reduce the external disturbance at nighttime; besides, considering that cardiovascular disease could have a significant effect on HR or HR variations, and the mechanism between them has been not well understood, we excluded these patients to reduce interference of confounding factors. It is important to carry out risk stratification in the early stage to reduce the mortality rate and improve the prognosis of ICU patients.

Several limitations of this study warrant discussion. Firstly, limited by the observational design, causal conclusions cannot be reached. However, a wide range of clinically important confounders were included to reach valid results. Secondly, based on data from a single academic medical center, the applicability of our findings was limited in other sites. Thirdly, confined by the availability of data, SAPS I score, instead of SAPS II, was adopted in the present study in assessing the disease severity along with SOFA score. Fourthly, the MIMIC II database is relatively old. We will continue to mine the updated MIMIC database in the subsequent study to further verify our research results. Finally, as limitations of a retrospective study, undetected potential confounders may also exist. Although multivariable models were adjusted for as many confounding factors as possible to diminish the possible influences, residual confounding factors may still exist and need to be investigated in the future.

Conclusions

The present study revealed the short-term and long-term prognostic value of the higher circadian variation of HR among critically ill patients, especially in individuals with higher SOFA scores. Therefore, the circadian variation of HR, as a cost-effective and readily available parameter, may serve as a potentially indicator of prognosis in the management of ICU patients. For the clinicians, further management, such as reducing the number of operations, reducing ward noise, or lowering light intensity level at night, may be considered to improve the circadian HR rhythm of ICU patients to improve their prognosis. Certainly, our findings are also needed to be confirmed by future studies, especially large prospective studies with a longer follow-up.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

AIDS:

Acquired immune deficiency syndrome

CI:

Confidence interval

ECG:

Electrocardiographic

HR:

Heart rate

HRT:

Heart rate turbulence

HRV:

Heart rate variability

ICD:

International Classification of Diseases-9 diagnoses

ICU:

Intensive care units

MIMIC-II:

Multiparameter Intelligent Monitoring in Intensive Care II

MIT:

Massachusetts Institute of Technology

MODS:

Multiple organ dysfunction syndrome

NHR:

Nighttime mean heart rate

OR:

Odds ratio

RHR:

Resting heart rate

SAPS-I:

Simplified Acute Physiology Score-I

SD:

Standard deviation

SOFA:

Sequential Organ Failure Assessment

References

  1. Camm AJ, Malik M, Bigger JT, Breithardt G, Cerutti S, Cohen RJ, et al. Heart rate variability: standards of measurement, physiological interpretation and clinical use. Task force of the European society of cardiology and the North American society of pacing and electrophysiology. Circulation. 1996;93(5):1043–65.

  2. Kleiger RE, Stein PK, Bigger JT. Heart rate variability: measurement and clinical utility. Ann Noninvasive Electrocardiol. 2005;10(1):88–101.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Goldstein B, Fiser DH, Kelly MM, Mickelsen D, Ruttimann U, Pollack MM. Decomplexification in critical illness and injury: relationship between heart rate variability, severity of illness, and outcome. Crit Care Med. 1998;26(2):352–7.

    Article  CAS  PubMed  Google Scholar 

  4. Winchell RJ, Hoyt DB. Spectral analysis of heart rate variability in the ICU: a measure of autonomic function. J Surg Res. 1996;63(1):11–6.

    Article  CAS  PubMed  Google Scholar 

  5. Yien HW, Hseu SS, Lee LC, Kuo TB, Lee TY, Chan SH. Spectral analysis of systemic arterial pressure and heart rate signals as a prognostic tool for the prediction of patient outcome in the intensive care unit. Crit Care Med. 1997;25(2):258–66.

    Article  CAS  PubMed  Google Scholar 

  6. Godin PJ, Fleisher LA, Eidsath A, Vandivier RW, Preas HL, Banks SM, et al. Experimental human endotoxemia increases cardiac regularity: results from a prospective, randomized, crossover trial. Crit Care Med. 1996;24(7):1117–24.

    Article  CAS  PubMed  Google Scholar 

  7. Annane D, Trabold F, Sharshar T, Jarrin I, Blanc AS, Raphael JC, et al. Inappropriate sympathetic activation at onset of septic shock: a spectral analysis approach. Am J Respir Crit Care Med. 1999;160(2):458–65.

    Article  CAS  PubMed  Google Scholar 

  8. Pontet J, Contreras P, Curbelo A, Medina J, Noveri S, Bentancourt S, et al. Heart rate variability as early marker of multiple organ dysfunction syndrome in septic patients. J Crit Care. 2003;18(3):156–63.

    Article  PubMed  Google Scholar 

  9. Barnaby D, Ferrick K, Kaplan DT, Shah S, Bijur P, Gallagher EJ. Heart rate variability in emergency department patients with sepsis. Acad Emerg Med. 2002;9(7):661–70.

    Article  PubMed  Google Scholar 

  10. Chen WL, Kuo CD. Characteristics of heart rate variability can predict impending septic shock in emergency department patients with sepsis. Acad Emerg Med. 2007;14(5):392–7.

    Article  PubMed  Google Scholar 

  11. Monfredi O, Lyashkov AE, Johnsen AB, Inada S, Schneider H, Wang R, et al. Biophysical characterization of the underappreciated and important relationship between heart rate variability and heart rate. Hypertension. 2014;64(6):1334–43.

    Article  CAS  PubMed  Google Scholar 

  12. Miwa Y, Miyakoshi M, Hoshida K, Yanagisawa R, Abe A, Tsukada T, et al. Heart rate turbulence can predict cardiac mortality following myocardial infarction in patients with diabetes mellitus. J Cardiovasc Electrophysiol. 2011;22(10):1135–40.

    Article  PubMed  Google Scholar 

  13. Stein PK, Barzilay JI. Relationship of abnormal heart rate turbulence and elevated CRP to cardiac mortality in low, intermediate, and high-risk older adults. J Cardiovasc Electrophysiol. 2011;22(2):122–7.

    PubMed  Google Scholar 

  14. Grander W, Mullauer K, Koller B, Tilg H, Dunser M. Heart rate before ICU discharge: a simple and readily available predictor of short- and long-term mortality from critical illness. Clin Res Cardiol. 2013;102(8):599–606.

    Article  PubMed  Google Scholar 

  15. Hoke RS, Muller-Werdan U, Lautenschlager C, Werdan K, Ebelt H. Heart rate as an independent risk factor in patients with multiple organ dysfunction: a prospective, observational study. Clin Res Cardiol. 2012;101(2):139–47.

    Article  PubMed  Google Scholar 

  16. Shashikumar SP, Stanley MD, Sadiq I, Li Q, Holder A, Clifford GD, et al. Early sepsis detection in critical care patients using multiscale blood pressure and heart rate dynamics. J Electrocardiol. 2017;50(6):739–43.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Sandfort V, Johnson AEW, Kunz LM, Vargas JD, Rosing DR. Prolonged elevated heart rate and 90-day survival in acutely Ill patients: data from the MIMIC-III database. J Intensive Care Med. 2018;34(8):622-9.

  18. A H, Inoue R, Ohkubo T, Kikuya M, Metoki H, Asayama K, et al. Predictive value of ambulatory heart rate in the Japanese general population: the Ohasama study. J Hypertens. 2008;26(8):1571–6.

    Article  CAS  Google Scholar 

  19. Hasanpour Mir M, Yousefshahi F, Abdollahi M, Ahmadi A, Nadjafi A, Mojtahedzadeh M. The predictive value of resting heart rate following osmotherapy in brain injury: back to basics. DARU. 2012;20(1):102.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Gazendam JAC, Van Dongen HPA, Grant DA, Freedman NS, Zwaveling JH, Schwab RJ. Altered circadian rhythmicity in patients in the ICU. Chest. 2013;144(2):483–9.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Paul T, Lemmer B. Disturbance of circadian rhythms in analgosedated intensive care unit patients with and without craniocerebral injury. Chronobiol Int. 2007;24(1):45–61.

    Article  CAS  PubMed  Google Scholar 

  22. Korompeli A, Muurlink O, Kavrochorianou N, Katsoulas T, Fildissis G, Baltopoulos G. Circadian disruption of ICU patients: a review of pathways, expression, and interventions. J Crit Care. 2017;38:269–77.

    Article  PubMed  Google Scholar 

  23. Telias I, Wilcox ME. Sleep and circadian rhythm in critical illness. Crit Care. 2019;23(1):82.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Saeed M, Villarroel M, Reisner AT, Clifford G, Lehman LW, Moody G, et al. Multiparameter Intelligent Monitoring in Intensive Care II: a public-access intensive care unit database. Crit Care Med. 2011;39(5):952–60.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Singer M, Deutschman CS, Seymour CW, Shankar-Hari M, Annane D, Bauer M, et al. The third international consensus definitions for sepsis and septic shock (Sepsis-3). JAMA. 2016;315(8):801–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Zheng CF, Liu WY, Zeng FF, Zheng MH, Shi HY, Zhou Y, et al. Prognostic value of platelet-to-lymphocyte ratios among critically ill patients with acute kidney injury. Crit Care. 2017;21(1):238.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Santhi N, Lazar A, McCabe P, Lo J, Groeger JA, Dijk DJ. Sex differences in the circadian regulation of sleep and waking cognition in humans. Proc Natl Acad Sci U S A. 2016;113(19):E2730-9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Schmidt C, Collette F, Cajochen C, Peigneux P. A time to think: circadian rhythms in human cognition. Cogn Neuropsychol. 2007;24(7):755–89.

    Article  PubMed  Google Scholar 

  29. Muto V, Jaspar M, Meyer C, Kussé C, Chellappa SL, Degueldre C, et al. Local modulation of human brain responses by circadian rhythmicity and sleep debt. Science. 2016;353(6300):687–90.

    Article  CAS  PubMed  Google Scholar 

  30. Patel M, Chipman J, Carlin BW, Shade D. Sleep in the intensive care unit setting. Crit Care Nurs Q. 2008;31(4):309–18 quiz 319-20.

    Article  PubMed  Google Scholar 

  31. van Someren EJ, Hagebeuk EE, Lijzenga C, Scheltens P, de Rooij SE, Jonker C, et al. Circadian rest-activity rhythm disturbances in alzheimer’s disease. Biol Psychiatry. 1996;40(4):259–70.

    Article  PubMed  Google Scholar 

  32. Bryant PA, Trinder J, Curtis N. Sick and tired: does sleep have a vital role in the immune system? Nat Rev Immunol. 2004;4(6):457–67.

    Article  CAS  PubMed  Google Scholar 

  33. Gao Y, Wang Q, Li J, Zhang J, Li R, Sun L, et al. Impact of mean arterial pressure fluctuation on mortality in critically Ill patients. Crit Care Med. 2018;46(12):e1167–74.

    Article  PubMed  Google Scholar 

  34. Li J, Li R, Gao Y, Zhang J, Zhao Y, Zhang X, et al. Nocturnal mean arterial pressure rising is associated with mortality in the intensive care unit: a retrospective cohort study. J Am Heart Assoc. 2019;8(19):e012388.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Nicolini P, Ciulla MM, De Asmundis C, Magrini F, Brugada P. The prognostic value of heart rate variability in the elderly, changing the perspective: from sympathovagal balance to chaos theory. Pacing Clin Electrophysiol. 2012;35(5):622–38.

    Article  PubMed  Google Scholar 

  36. de Castilho FM, Ribeiro ALP, Nobre V, Barros G.de Sousa MR. Heart rate variability as predictor of mortality in sepsis: a systematic review. PLoS One. 2018;13(9):e0203487.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. de Castilho FM, Ribeiro ALP, da Silva JLP, Nobre V.de Sousa MR. Heart rate variability as predictor of mortality in sepsis: a prospective cohort study. PLoS One. 2017;12(6):e0180060.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Bauer A, Malik M, Schmidt G, Barthel P, Bonnemeier H, Cygankiewicz I, et al. Heart rate turbulence: standards of measurement, physiological interpretation, and clinical use: international society for holter and noninvasive electrophysiology consensus. J Am Coll Cardiol. 2008;52(17):1353–65.

    Article  PubMed  Google Scholar 

  39. Marynissen T, Flore V, Heidbuchel H, Nuyens D, Ector J, Willems R. Heart rate turbulence predicts ICD-resistant mortality in ischaemic heart disease. Europace. 2014;16(7):1069–77.

    Article  PubMed  Google Scholar 

  40. Xu T, Zhan Y, Xiong J, Lu N, He Z, Su X, et al. The relationship between heart rate and mortality of patients with acute coronary syndromes in the coronary intervention era: meta-analysis. Medicine (Baltimore). 2016;95(46):e5371.

    Article  Google Scholar 

  41. Davoudi A, Malhotra KR, Shickel B. Intelligent ICU for autonomous patient monitoring using pervasive sensing and deep learning. Sci Rep. 2019;9(1):8020.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  42. Goldberger JJ, Cain ME, Hohnloser SH, Kadish AH, Knight BP, Lauer MS, et al. American heart association/American college of cardiology foundation/heart rhythm society scientific statement on noninvasive risk stratification techniques for identifying patients at risk for sudden cardiac death: a scientific statement from the American heart association council on clinical cardiology committee on electrocardiography and arrhythmias and council on epidemiology and prevention. Circulation. 2008;118(14):1497–518.

    Article  PubMed  Google Scholar 

  43. Shaw E, Tofler GH. Circadian rhythm and cardiovascular disease. Curr Atheroscler Rep. 2009;11(4):289–95.

    Article  PubMed  Google Scholar 

  44. Palatini P, Julius S. Heart rate and the cardiovascular risk. J Hypertens. 1997;15(1):3–17.

    Article  CAS  PubMed  Google Scholar 

  45. Palatini P, Julius S. Association of tachycardia with morbidity and mortality: pathophysiological considerations. J Hum Hypertens. 1997;11(Suppl 1):S19-27.

    PubMed  Google Scholar 

  46. Ben-Dov IZ, Kark JD, Ben-Ishay D, Mekler J, Ben-Arie L, Bursztyn M. Blunted heart rate dip during sleep and all-cause mortality. Arch Intern Med. 2007;167(19):2116–21.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

We appreciate the researchers at the MIT Laboratory for Computational Physiology and collaborating research groups to keep MIMIC-II databases available.

Funding

This study was funded by the National Natural Science Foundation of China (81670049).

Author information

Authors and Affiliations

  1. Department of Critical Care Medicine, the Second Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China

    Jingjing Zhang, Linyun Du, Jiamei Li, Ruohan Li, Xuting Jin, Jiajia Ren, Ya Gao & Xiaochuang Wang

Authors
  1. Jingjing Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Linyun Du
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Jiamei Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ruohan Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Xuting Jin
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Jiajia Ren
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Ya Gao
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Xiaochuang Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

XCW, JJZ, and LYD take responsibility for the integrity of the data, interpretation, and analysis. XCW, JJZ, LYD, JML, RHL, XTJ, JJR and YG contributed substantially to the study design, data interpretation, and writing of the manuscript. XCW, JJZ, and LYD performed the statistical analysis and data synthesis. All authors approved the final version of the manuscript.

Corresponding author

Correspondence to Xiaochuang Wang.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional research committee (Cardif Metropolitan University ethics committee—17/4/02R) and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. Because the project did not impact clinical care and all protected health information was deidentified, informed consent was waived for the present study by the Institutional Review Board from the Second Affiliated Hospital of Xi’an Jiaotong University. The study was designed and conducted in accordance with relevant guidelines and regulations (Declaration of Helsinki), and the study protocol was approved by Institutional Review Board from the Second Affiliated Hospital of Xi’an Jiaotong University.

Consent for publication

Consent for publication is not applicable.

Competing interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table 1.

The ICU admission types of the present study subjects.

Additional file 2: Table 2.

The univariate analysis of factors associate ICU, hospital mortality in the wholepopulations.

Additional file 3: Table 3.

The univariate analysis of factors associate 30-day and1-year in the whole populations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, J., Du, L., Li, J. et al. Association between circadian variation of heart rate and mortality among critically ill patients: a retrospective cohort study. BMC Anesthesiol 22, 45 (2022). https://doi.org/10.1186/s12871-022-01586-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12871-022-01586-9

Keywords

BMC Anesthesiology

ISSN: 1471-2253

Contact us