Skip to main content
  • Research article
  • Open access
  • Published:

High-frequency power of heart rate variability can predict the outcome of thoracic surgical patients with acute respiratory distress syndrome on admission to the intensive care unit: a prospective, single-centric, case-controlled study

Abstract

Background

The morbidity and mortality of acute respiratory distress syndrome (ARDS) remains high, and the strategic focus of ARDS research has shifted toward identifying patients at high risk of mortality early in the course of illness. This study intended to identify the heart rate variability (HRV) measure that can predict the outcome of patients with ARDS on admission to the surgical intensive care unit (SICU).

Methods

Patients who had lung or esophageal cancer surgery were included either in the ARDS group (n = 21) if they developed ARDS after surgery or in the control group (n = 11) if they did not. The ARDS patients were further stratified into survivors and non-survivors subgroups according to their outcomes. HRV measures of the patients were used for statistical analysis.

Results

The mean RR interval (mRRI), high-frequency power (HFP) and product of low-/high-frequency power ratio tidal volume and tidal volume (LHR*VT) were significantly lower (p < 0.05), while the normalized HFP to VT ratio (nHFP/VT) was significantly higher in the ARDS patients (p = 0.011). The total power (TP), low-frequency power (LFP), HFP and HFP/VT were all significantly higher in the non-survived ARDS patients, whereas Richmond Agitation-Sedation Scale (RASS) was significantly lower in the non-survived ARDS patients. After adjustment for RASS, age and gender, firth logistic regression analysis identified the HFP, TP as the significant independent predictors of mortality for ARDS patients.

Conclusions

The vagal modulation of thoracic surgical patients with ARDS was enhanced as compared to that of non-ARDS patients, and the non-survived ARDS patients had higher vagal activity than those of survived ARDS patients. The vagal modulation-related parameters such as TP and HFP were independent predictors of mortality in patients with ARDS on admission to the SICU, and the HFP was found to be the best predictor of mortality for those ARDS patients. Increased vagal modulation might be an indicator for poor prognosis in critically ill patients following thoracic surgery.

Peer Review reports

Background

Acute respiratory distress syndrome (ARDS) is a life threatening lung condition that affects both medical and surgical patients [1, 2]. Postoperative ARDS is a recognized complication after lung resection or esophagectomy for cancer [3,4,5]. The incidence of ARDS after these procedures was 1–3% [3, 6], and the overall hospital mortality rates associated with ARDS were over 40% [5,6,7,8].

Many scoring systems have been developed and used in the intensive care unit (ICU) to monitor the treatment and predict the outcomes of patients with ARDS. For instance, Navarrete-Navarro et al. [9] showed that the ICU mortality of trauma patients with ARDS was related to the Acute Physiology and Chronic Health Evaluation (APACHE) III score and PaO2/FIO2 on the 3rd day. Luecke et al. [10] showed that the APACHE II was the only clinical predictor for mortality on admission to the ICU. Bauman et al. [11] showed that the lung injury score was reliable for predicting mortality in critically ill surgical patients.

The autonomic nervous system is important for the maintenance of homeostasis. A wide range of pathophysiologic conditions can influence the balance in the autonomic nervous system. Heart rate variability (HRV) refers to the variation of interbeat intervals due to respiration and other physiological inputs. The assessment of HRV is based on the analysis of consecutive normal R-R intervals, and may provide clinicians with quantitative information about cardiac autonomic nervous modulation. HRV can be quantified by time domain and frequency domain parameters. Thus, HRV analysis has become a useful non-invasive technique for evaluation of cardiac autonomic nervous modulation in many kinds of illnesses [12].

Previous studies have reported that HRV measures are decreased in survivors of cardiac arrest [13]. Some studies found that the standard deviation of normal-to-normal QRS intervals (SDNN), one of six commonly used time-domain HRV measures, is a significant parameter for long-term prognosis in patients with acute myocardial infarction [14, 15]. Tibby et al. [16] showed that the loss of HRV occurs with increasing organ failure, and that this effect is better demonstrated by the power-law model of HRV than by other measures of HRV. Annane et al. [17] reported that the low-frequency (LF) component of the oscillations in heart rate and diastolic blood pressure variability are dramatically reduced in patients with septic shock, and that the decrease in the LF component may identify septic patients with a high level of sympathetic activation. Pontet et al. [18] demonstrated that reduction of HRV in patients on ICU admission may be useful in identifying septic patients at risk for the development of multiple organ dysfunction syndromes. In patients with severe sepsis, depressed sympathetic modulation has been found to be indicative of poor outcome [19]. In patients in the emergency department who have been successfully resuscitated after out-of-hospital cardiac arrest, the initial HRV measures resemble those of severe sepsis and are capable of predicting 24-h mortality [20, 21]. In the emergency department, the characteristics of HRV can predict impending septic shock in patients with sepsis [22]. These studies suggest that HRV measures might be useful for the monitoring of patients with various kinds of diseases, including ARDS, and for predicting their prognosis.

Despite significant advances in the management of patients with ARDS in the past fifty years, the morbidity and mortality of ARDS remains high, and the strategic focus of ARDS research has shifted toward identifying patients with or at high risk of ARDS early in the course of their illnesses [23]. Thus, the aim of this study was to evaluate the possibility of using HRV measures to predict the outcome in patients with ARDS on admission to the surgical intensive unit (SICU).

Methods

Research Design

This was a prospective, single-centric, case-controlled study. The study protocol has been approved by the Institutional Review Boards of National Taiwan University Hospital (NTUH200808065R) and Taipei Veterans General Hospital (VGHIRB97–01-02A), and written informed consent was obtained from the patients or the next of kin of the patients before their enrollment in the study.

Study Participants

This study was conducted in the SICU of National Taiwan University Hospital. All the recruited patients were over 18 years old. They were partitioned into an ARDS group and a control group. The ARDS group consisted of patients who had received thoracic surgery because of lung or esophageal cancer and were transferred to the SICU for intensive care because of ARDS later. The control group consisted of patients who had received thoracic surgery because of lung or esophageal cancer and were transferred to the SICU for post-operative care. The ARDS was diagnosed according to the Berlin Definition [2]. A draft definition proposed 3 mutually exclusive categories of ARDS based on the degree of hypoxemia: mild (200 mmHg < PaO2/FIO2 ≤ 300 mmHg), moderate (100 mmHg < PaO2/FIO2 ≤ 200 mmHg), and severe (PaO2/FIO2 ≤ 100 mmHg) with 4 ancillary variables for severe ARDS: radiographic severity, respiratory system compliance (≤ 40 mL/cm H2O), positive end-expiratory pressure (≥ 10 cm H2O), and corrected expired volume per minute (≥ 10 L/min). Patients with severe coronary artery disease, persistent arrhythmia, cardiac pacing, diabetes mellitus, cerebral vascular accident, or major diseases of kidney or autoimmune system were excluded from the study.

Physiological Measurements

All patients required intubation and mechanical ventilation. The acute lung injury score (ALIS) [24], alveolar-arterial oxygen difference (AaDO2), APACHE II [25], demographic data, vital signs, medications, ventilator parameters including respiratory rate (RR), inspiratory pressure (Pinsp), inspiration time (Tinsp), tidal volume (VT), positive end-expiratory pressure (PEEP), minute ventilation (MV), fraction of inspired oxygen (FIO2), dynamic compliance (Cdyn), and arterial blood gases data were recorded within 4 h of admission to the SICU. Electrocardiographic (ECG) signals were recorded for 12 min in supine position using the MP35 multichannel recorder (BIOPAC Systems, Goleta, CA, USA). The primary outcome measure was survival at discharge from the SICU. Patients who survived the ARDS were classified as the survivors, while those patients who did not survive the ARDS were classified as non-survivors. Clinical information obtained in the SICU and the outcomes of the patients were blinded to the investigator who performed the HRV analysis and clinical data collection.

HRV Analysis

Power spectral density (PSD) analysis by using parametric autoregressive modeling or non-parametric Fourier transform provides the basic information about the power (variance) distributes of experimental data as a function of frequency. The Fourier transform decomposes a signal into its constituent frequencies. The Fourier transform of a signal is a complex-valued function of frequency, whose absolute value represents the amount of that particular frequency present in the signal. The power of a particular frequency band is calculated by integrating the PSD within that frequency band. Thus, the Fourier transform is the frequency domain representation of the original signal [26].

The method used for the HRV analysis adhered to the standards laid down by the Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology [12]. All ectopic beats were removed and the missing data were replaced by interpolated beats derived from the nearest valid data. Any patient with more than 5% of ectopic beats was excluded from the study. The last 512 stationary RR intervals (RRI) were used for HRV analysis.

The time domain HRV measures included mRRI, standard deviation of RRI (SDRR), coefficient of variation of RRI (CVRR), and root mean squared successive difference of RRI (RMSSD). The power spectrum of RRI was obtained by means of fast Fourier transformation (Mathcad 15, Mathsoft Inc.). The areas under the spectral peaks within the range of 0.01–0.04 Hz, 0.04–0.15 Hz, 0.15–0.5 Hz, and 0.01–0.5 Hz were defined as the very low-frequency power (VLFP), LFP, high frequency power (HFP), and TP, respectively. The normalized VLFP (nVLFP = VLFP/TP × 100) was used as the index of renin-angiotensin-aldosterone modulation, thermal regulation and vagal withdrawal; the normalized LFP (nLFP = LFP/TP × 100) as the index of combined sympathetic and vagal activities; the normalized HFP (nHFP = HFP/TP × 100) and HFP as the index of cardiac vagal activity; and the low-/high-frequency power ratio (LHR = LFP/HFP) as the index of sympathovagal balance. Since the HFP and nHFP can be influenced by carbon dioxide, respiratory rate and tidal volume (VT) [27], the HFP/VT and nHFP/VT were used as the VT-corrected indices of vagal modulation, and the LHR*VT was used as the VT-corrected index of sympathovagal balance [20].

Statistical Analysis

Continuous variables were compared between the two groups of patients using Mann–Whitney U test. Fisher’s exact test was used for comparisons of categorical data. Pearson product moment analysis was performed to find the correlations among HRV measures and clinical variables of the patients. The Firth logistic regression and multiple linear regression models were employed to identify the most important independent predictors of mortality and ICU stay day, respectively. Firth’s penalized likelihood approach was used to reduce the bias of the parameter estimates due to small sample size. Significant HRV measures in the univariate analyses of differences between the two subgroups of ARDS patients were included in the multivariate analyses. In this study, the 10 events per variable (EPV − 10) rule [28, 29] was used in logistic regression model to avoid the overfitting problem, which is a generally accepted criterion in multivariable analysis. With a sample size of 32 patients (11 patients in the control group and 21 patients in the ARDS group), the inclusion of only three variables in the multivariate analysis to ensure adequate statistical power is reasonable. The binary receiver operating characteristic (ROC) curves for statistically independent variables associated with mortality were drawn. All statistical analyses were performed using R software (version i386 3.3.2, https://www.r-project.org/) and the contributed R package logistic for Firth’s penalized-likelihood logistic regression. A p < 0.05 was considered statistically significant.

Results

During the 2 years study period, 22 patients developed ARDS among the 817 patients that received lung or esophageal surgery. Thus the incidence rate of ARDS after thoracic surgery for lung or esophageal surgery was 2.69%. Among those 22 ARDS patients, 15 were included in the ARDS group while 7 were not because of arrhythmia. Six patients who were not operated during this admission but had thoracic surgery within past 1 year and developed ARDS due to rapidly progressing pneumonia were also admitted and included in the ARDS group. Thus, there were 11 patients in the non-ARDS group and 21 patients in the ARDS group.

Fentanyl (5~ 10 μg·hr-1) was administered to the patients during the post-operative period for post-surgical pain killing. The Fentanyl was given to the non-ARDS patients for 24–48 h, while it was given to the ARDS patients for more than 48 h. Midazolam (1~ 2 mg/h) was used to maintain Richmond Agitation-Sedation Scale (RASS) 0 ~ − 2 in both groups, and dopamine and norepinephrine (0.05 ~ 0.11 μg·hr− 1·min− 1) were used in 8 patients in the ARDS group to keep mean blood pressure > 60 mmHg to maintain proper cerebral perfusion pressure. The ECG recordings were taken within 4 h of transfer to the SICU in all patients.

In the control group, 6 patients had lung cancer and 5 patients had esophageal cancer. In the ARDS group, 9 patients had lung cancer, 5 patients had esophageal cancer, 1 patient had esophageal rupture, and 6 patients had pneumonia-induced ARDS after the thoracic surgery. Six patients out of 11 patients in the control group were classified as American Society of Anesthesiologists (ASA) I, and 5 patients as ASA II. For the 21 patients in the ARDS group, 1 patient was classified as ASA I, 5 patients as ASA II, 4 patients as ASA III. Six patients in the ARDS group had no ASA classification because they were not operated during this admission, but had thoracic surgery within past 1 year, and were admitted due to ARDS caused by rapidly progressing pneumonia (Table 1).

Table 1 Comparison between patients with and without ARDS

Ten patients in the control group received video-assisted thoracic surgery (VATS) for lobectomy or esophageal resection, and 1 patient received VATS for laminectomy because of lung cancer metastasis. In the ARDS group, 12 patients received VATS for lobectomy or esophageal resection, 3 patients received open lung surgery, and 6 patients received thoracic surgery within the past 1 year who were admitted to the SICU because of ARDS caused by pneumonia (Table 1).

In the control group, 2 patients received bronchus blocker, and 9 patients received double-lumen intubation. In the ARDS group, 15 patients received thoracic surgery during this admission; among them, 5 patients received bronchus blocker and 10 patients received double-lumen intubation. All the 26 patients who received thoracic surgery during this admission were operated under Desflurane and Sevoflurane anesthesia (minimum alveolar concentration: 0.8–1.1 Vol%). Five out of the 11 patients in the control group received blood transfusion; the transfusion rate in the control group was 45.5%. In the ARDS group, 6 out of the 15 patients who received thoracic surgery during this admission had blood transfusion. Therefore, the transfusion rate was 40% in the ARDS group. The duration of surgery in the control group was 255 (245–450) minutes (median and interquartile range), whereas the duration of surgery in the ARDS group was 250 (165–375) minutes (Table 1).

Tables 1 shows that the ARDS patients had higher RR, required higher FIO2 and PEEP during ventilation. The VT, Cdyn, PaO2 and PaO2/FIO2, intraoperative intake/output (I/O) of the ARDS patients was significantly lower (p < 0.01), whereas the arterial partial pressure of carbon dioxide (PaCO2), bicarbonate (HCO3−), ASA classification, and AaDO2 of the ARDS patients were significantly higher than those of non-ARDS patients (p < 0.05). The mRRI, HFP and LHR*VT were significantly lower (p < 0.05), while the nHFP/VT was significantly higher in the ARDS patients (p = 0.011). Table 2 shows that there were no significant differences in most ventilator settings, disease severity scores including ALIS, AaDO2, APACHE II, intraoperative anesthesia control, blood transfusion rate, fluid balance, and blood gas data between the survivors and non-survivors of ARDS. However, the Pinsp, PaCO2, HCO3−, TP, LFP, HFP, and HFP/VT of the non-survivors were higher than those of the survivors (p < 0.05), and RASS of the non-survivors was significantly lower than that of the survivors (p = 0.040).

Table 2 Comparison between survivors and non-survivors in patients with ARDS

Pearson product moment analysis among the significantly different HRV measures between the survived and non-survived ARDS patients showed that the TP, LFP, HFP, and HFP/VT also correlated significantly and positively with one another (Table 3), but these HRV measures did not correlate with RASS.

Table 3 Correlation coefficient in the Pearson product moment correlation analyses of HRV parameters and RASS in patients with ARDS

Firth logistic regression analysis was used to assess the risk of mortality (dependent variable) in ARDS patients. The ASA classification, intraoperative anesthesia, fluid balance, drugs, etc., were not significantly different between the survivors and non- survivors of ARDS; only the RASS of the non-survivors was significantly lower than that of the survivors (Table 2). In univariate analysis, the TP, HFP and RASS were found to be the significant predictors of mortality for ARDS patients in the SICU (p < 0.05) (Table 4). Since the TP and HFP were highly correlated with each other while the RASS did not correlate with either TP or HFP (Table 3), two separate models in multivariate analysis including HFP and TP adjusted for RASS, age and gender were used. After adjustment for RASS, age and gender, the HFP (p = 0.025) and TP (p = 0.024) were found to be the significant independent predictors of mortality in ARDS patients in the SICU (Table 4).

Table 4 Firth logistic regression analyses for predicting the mortality of patients with ARDS

As depicted in Fig. 1, the AUCs and 95% confidence intervals (CI) of TP and HFP in predicting mortality in ARDS patients on admission to the SICU were 0.850 (95% CI = 0.685–1.000; p = 0.021) and 0.888 (95% CI = 0.741–1.000; p = 0.010), respectively. Among all variables assessed in this study, HFP was the best predictor in predicting mortality in ARDS patients on admission to the SICU.

Fig. 1
figure 1

The areas under the ROC curves (AUC) and 95% confidence intervals (CI) of TP and HFP in predicting the mortality of ARDS patients on admission to the SICU were 0.850 (95% CI = 0.685–1.000; p = 0.021), and 0.888 (95% CI = 0.741–1.000; p = 0.010), respectively

Linear regression analysis was used to assess the association between ICU stay days (dependent variable) and HRV measures in ARDS patients. It was found that the TP (p = 0.681), LFP (p = 0.740), HFP (p = 0.686) and RASS (p = 0.108) were not significantly associated with the days of ICU stay of ARDS patients in the SICU (Table 5).

Table 5 Linear regression analyses for predicting the ICU stay days in patients with ARDS

Discussions

Spectral HRV analysis allows us to differentiate between branches of the autonomic nervous system [12] and to assess the autonomic modulation of critically ill patients [16,17,18,19,20,21,22]. This study compared the clinical characteristics and HRV measures between non-ARDS and ARDS patients, and between ARDS patients who survived and did not survive. It was found that the HFP was significantly lower (p = 0.043) while the nHFP/VT was significantly higher (p = 0.011) in the ARDS group compared with those in the non-ARDS group, and that the non-survived ARDS patients had higher TP, LFP, HFP, and HFP/VT than the survived ARDS patients had. After adjustment for RASS, age and gender, firth logistic regression analysis identified HFP and TP as the independent predictors of mortality in ARDS patients on admission to the SICU. The HFP was found to be the best predictor of mortality in ARDS patients.

Many studies have proposed that clinical laboratory data and HRV measures can be used to predict the prognosis of certain diseases. For instance, Luhr et al. [30] showed that age and acute physiologic score were associated with mortality in the ARDS patients. Phillips et al. [31] had shown that the lung injury score of non-survivors was significantly higher than that of survivors in patients with ARDS. Swaroopa et al. [32] had shown that the APACHE II score was significantly higher in non-survivors as compared to survivors in patients with ARDS. In this study it was found that the APACHE II, ALIS, PaO2/FIO2 and AaDO2 were not significantly different between the survived and non-survived ARDS patients; instead, the HFP and TP were found to be the independent predictors of mortality in ARDS patients on admission to the SICU (p < 0.05). This result suggested that the non-invasive HRV indices such as HFP and TP were more sensitive than the conventional scores or indices such as APACHE II, ALIS, PaO2/FIO2 and AaDO2 in predicting the outcome of surgical patients with ARDS. Because of its non-invasiveness and close relation with the autonomic control of heart rate, HRV is often used in the monitoring and prognosis of critically ill patients. For instance, Longin et al. [33] studied the differences in HRV variables during extracorporeal membrane oxygenation (ECMO) therapy between survivors (n = 6) and non-survivors (n = 7) and found that the survivors showed a tendentious, but not significant, decrease in HRV variables (ECMO-start vs. ECMO-stop), while the non-survivors had elevated HRV values. Winchell et al. [34] showed that both low TP and LHR were associated with increased mortality, while sympathetic predominance favored survival in the ICU setting. Chen et al. [35] found that a worse oxygenation status is associated with increased cardiac vagal and decreased cardiac sympathetic activities in patients with chronic obstructive pulmonary disease (COPD). For patients with sepsis in the emergency department, Chen et al. [36] also found that the SDNN, TP, VLFP, LFP, and LHR were all significantly lower while the nHFP was significantly higher in the non-survivors of sepsis. In another study, Chen et al. [20] further showed that the LFP, nLFP and LHR were significantly lower while the RMSSD, HFP, HFP/VT, nHFP, and nHFP/VT were significantly higher in the non-survivors of patients with out-of-hospital cardiac arrest. These results suggested that HRV measures can be used in the monitoring and prognosis of critically ill patients, and that increased cardiac vagal and decreased cardiac sympathetic activities might be indicative of poor prognosis. In accordance with these studies, the present study also found that the ARDS patients had higher vagal modulation and smaller sympathetic modulation than the non-ARDS patients, and that the non-survived ARDS patients had higher vagal modulation than the survived ARDS patients. Further, Firth logistic regression analysis identified HFP and TP as the independent predictors of mortality in ARDS patients on admission to the SICU (Table 4). The HFP, an indicator of vagal modulation, appeared to be the best predictor of mortality in ARDS patients (p = 0.02). Thus the HFP might be useful in the monitoring and prediction of outcome in various critically ill patients, including ARDS, and be used in the prediction of mortality on their admission to the SICU.

The mechanism underlying the increase in vagal modulation in critically ill patients has not been well understood yet. It has been shown that direct electrical stimulation of the peripheral vagus nerve in vivo during lethal endotoxaemia in rats could inhibit tumor necrosis factor (TNF) synthesis in liver, attenuated peak serum TNF amounts, and prevented the development of shock [37]. Thus, the increased vagal modulation in ARDS patients might be one of the pathophysiological responses of the patient to counteract the cytokines synthesis, attenuate the amount of cytokines in the serum, and prevent the cytokines-induced systemic inflammation in ARDS. Since the non-survived ARDS patients had increased vagal modulation relative to the survived ARDS patients, the use of anticholinergic and sympathomimetic agents to correct their autonomic dysfunction might be considered in patients with ARDS after thoracic surgery. Further studies are warranted to investigate the underlying mechanism of increased vagal modulation in critically ill patients, including ARDS.

The HRV parameters in ARDS patients can be influenced by many interventions performed in the ICU, such as sedation, catecholamine therapy, volume management, ventilator settings, etc. Previous study has indicated that the increase in VT can increase the HFP component of HRV [26]. To minimize the effect of VT on the indicator of vagal modulation, a new indicator nHFP/VT was used to represent the VT-corrected indicator of vagal modulation. It was found that the ARDS patients still had higher vagal modulation than the non-ARDS patients. Some studies had reported that midazolam and inotropes might change the autonomic nervous activities of the patients [38, 39]. Though the midazolam and norepinephrine were used more frequently in the ARDS patients than the non-ARDS patients, the use of midazolam and norepinephrine were not significantly different between survivors and non-survivors of ARDS patients. Thus, the significantly increase in vagal modulation in the non-survivors of ARDS should not be caused by the use of midazolam and norepinephrine. In ARDS patients, PEEP has been widely used to improve arterial oxygenation; however, the increase in PEEP can reduce the baroreflex sensitivity, and can alter the autonomic nervous function in these patients [40]. The present study found that ARDS patients required higher PEEP during mechanical ventilation than the non-ARDS patients. Since the level of PEEP was not significantly different between the survivors and non-survivors of ARDS patients, the significantly increased vagal modulation in the non-survivors of ARDS patients should not be caused by the use of higher PEEP in those patients.

The major limitation of this study was the small number of ARDS patients. Since the incidence of ARDS is not very high in the post-operative patients with lung or esophageal cancer, it was not easy to collect a large number of cases for statistical analysis in a study period of 2 years.

Conclusions

The vagal modulation of thoracic surgical patients with ARDS was enhanced compared with non-ARDS patients. Non-survived ARDS patients had higher vagal activity than survived ARDS patients. The vagal modulation-related parameter such as HFP and TP might be used as predictors to identify ARDS patients with high risk of mortality on admission to the SICU, especially the HFP. Increased vagal modulation might be an indicator for poor prognosis in critically ill patients following thoracic surgery, and anticholinergic and sympathomimetic agents might be tried to correct the autonomic dysfunction of ARDS patients.

Abbreviations

AaDO2 :

Alveolar-arterial oxygen difference

ALIS:

Acute lung injury score

APACHE:

Acute Physiology and Chronic Health Evaluation

ARDS:

Acute respiratory distress syndrome

ASA:

American Society of Anesthesiologists

Cdyn:

Dynamic compliance

COPD:

Chronic obstructive pulmonary disease

CVRR :

Coefficient of variation of RRI

ECG:

Electrocardiographic

ECMO:

Extracorporeal membrane oxygenation

EPV − 10:

10 events per variable

FIO2 :

Fraction of inspired oxygen

HCO3 − :

Bicarbonate

HFP:

High-frequency power

HRV:

Heart rate variability

I/O:

Intake/output

ICU:

Intensive care unit

LF:

Low-frequency

LFP:

Low-frequency power

LHR:

Low-/high-frequency power ratio

mRRI:

Mean RR interval

MV:

Minute ventilation

nHFP:

Normalized HFP

nLFP:

Normalized LFP

nVLFP:

Normalized VLFP

PaCO2 :

Arterial partial pressure of carbon dioxide

PaO2 :

Arterial partial pressure of oxygen

PEEP:

Positive end-expiratory pressure

Pinsp :

Inspiratory pressure

PSD:

Power spectral density

RASS:

Richmond Agitation-Sedation Scale

RMSSD:

Root mean squared successive difference of RRI

ROC:

Receiver operating characteristic

RR:

Respiratory rate

RRI:

RR intervals

SDNN:

Standard deviation of normal-to-normal QRS intervals

SDRR :

Standard deviation of RRI

SICU:

Surgical intensive care unit

Tinsp :

Inspiration time

TNF:

Tumor necrosis factor

TP:

Total power

VATS:

Video-assisted thoracic surgery

VLFP:

Very low-frequency power

VT :

Tidal volume

References

  1. Ashbaugh DG, Bigelow DB, Petty TL. et. al. Acute respiratory distress in adults. Lancet. 1967;2:319–23.

    Article  CAS  PubMed  Google Scholar 

  2. Definition Task Force ARDS, Ranieri VM, Rubenfeld GD, et al. Acute respiratory distress syndrome: the Berlin Definition. JAMA. 2012;307:2526–33.

    Google Scholar 

  3. D'Journo XB, Michelet P, Avaro JP, et al. Respiratory complications after oesophagectomy for cancer. Rev Mal Respir. 2008;25:683–94.

    Article  PubMed  Google Scholar 

  4. Kim JB, Lee SW, Park SI, et al. Risk factor analysis for postoperative acute respiratory distress syndrome and early mortality after pneumonectomy: the predictive value of preoperative lung perfusion distribution. J Thorac Cardiovasc Surg. 2010;140:26–31.

    Article  PubMed  Google Scholar 

  5. Kometani T, Okamoto T, Yoshida S, et al. Acute respiratory distress syndrome after pulmonary resection. Gen Thorac Cardiovasc Surg. 2013;61:504–12.

    Article  PubMed  Google Scholar 

  6. Dulu A, Pastores SM, Park B, et al. Prevalence and mortality of acute lung injury and ARDS after lung resection. Chest. 2006;130:73–8.

    Article  PubMed  Google Scholar 

  7. Ware LB, Matthay MA. The acute respiratory distress syndrome. N Engl J Med. 2000;342:1334–49.

    Article  CAS  PubMed  Google Scholar 

  8. Tandon S, Batchelor A, Bullock R, et al. Peri-operative risk factors for acute lung injury after elective oesophagectomy. Br J Anaesth. 2001;86:633–8.

    Article  CAS  PubMed  Google Scholar 

  9. Navarrete-Navarro P, Ruiz-Bailén M, Rivera-Fernández R, et al. Acute respiratory distress syndrome in trauma patients: ICU mortality and prediction factors. Intens Care Med. 2000;26:1624–9.

    Article  CAS  Google Scholar 

  10. Luecke T, Muench E, Roth H, Friess U, et al. Predictors of mortality in ARDS patients referred to a tertiary care centre: a pilot study. Eur J Anaesthesiol. 2006;23:403–10.

    Article  PubMed  Google Scholar 

  11. Bauman ZM, Gassner MY, Coughlin MA, et al. Lung injury prediction score is useful in predicting acute respiratory distress syndrome and mortality in surgical critical care patients. Crit Care Res Pract. 2015;2015:157408. https://doi.org/10.1155/2015/157408.

    PubMed  PubMed Central  Google Scholar 

  12. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. Heart rate variability: standards of measurement, physiological interpretation and clinical use. Eur Heart J. 1996;17(3):354-81. Circulation. 1996;93:1043–65.

    Article  Google Scholar 

  13. Huikuri HV, Linnaluoto MK, Seppanen T, et al. Circadian rhythm of heart rate variability in survivors of cardiac arrest. Am J Cardiol. 1992;70:610–5.

    Article  CAS  PubMed  Google Scholar 

  14. Liu PY, Tsai WC, Lin LJ, et al. Time domain heart rate variability as a predictor of long-term prognosis after acute myocardial infarction. J Formos Med Assoc. 2003;102:474–9.

    PubMed  Google Scholar 

  15. Erdogan A, Coch M, Bilgin M, et al. Prognostic value of heart rate variability after acute myocardial infarction in the era of immediate reperfusion. Herzschrittmacherther Elektrophysiol. 2008;19:161–8.

    Article  PubMed  Google Scholar 

  16. Tibby SM, Frndova H, Durward A, Cox PN. Novel method to quantify loss of heart rate variability in pediatric multiple organ failure. Crit Care Med. 2003;31:2059–67.

    Article  PubMed  Google Scholar 

  17. Annane D, Trabold F, Sharshar T, et al. Inappropriate sympathetic activation at onset of septic shock: a spectral analysis approach. Am J Respir Crit Care Med. 1999;160:458–65.

    Article  CAS  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  19. Yien HW, Hseu SS, Lee LC, et al. 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:258–66.

    Article  CAS  PubMed  Google Scholar 

  20. Chen WL, Tsai TH, Huang CC, et al. Heart rate variability predicts short-term outcome for successfully resuscitated patients with out-of-hospital cardiac arrest. Resuscitation. 2009;80:1114–8.

    Article  PubMed  Google Scholar 

  21. Chen WL, Shen YS, Huang CC, et al. Postresuscitation autonomic nervous modulation after cardiac arrest resembles that of severe sepsis. Am J Emerg Med. 2012;30:143–50.

    Article  PubMed  Google Scholar 

  22. 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:392–7.

    Article  PubMed  Google Scholar 

  23. Yadav H, Thompson BT, Gajic O. Fifty years of research in ARDS: Is acute respiratory distress syndrome a preventable disease? Am J Respir Crit Care Med. 2017;195:725–36.

    Article  CAS  PubMed  Google Scholar 

  24. Murray JF, Matthay MA, Luce JM, et al. An expanded definition of the adult respiratory distress syndrome. Am Rev Respir Dis. 1988;138:720–3.

    Article  CAS  PubMed  Google Scholar 

  25. Knaus WA, Draper EA, Wagner DP, Zimmerman JE. APACHE II: a severity of disease classification system. Crit Care Med. 1985;13:818–29.

    Article  CAS  PubMed  Google Scholar 

  26. McClellan JH, Schafer RW, Yoder MA. Signal Processing First. New Jersey: Prentice Hall; 2003.

    Google Scholar 

  27. Pöyhönen M, Syväoja S, Hartikainen J, et al. The effect of carbon dioxide, respiratory rate and tidal volume on human heart rate variability. Acta Anaesthesiol Scand. 2004;48:93–101.

    Article  PubMed  Google Scholar 

  28. Peduzzi P, Concato J, Kemper E, et al. A simulation study of the number of events per variable in logistic regression analysis. J Clinl Epidem. 1996;49:1373–9.

    Article  CAS  Google Scholar 

  29. Hosmer DW, Lemeshow S. Applied Logistic Regression. 2nd ed. Hoboken, NJ: Wiley Online Library; 2000.

    Book  Google Scholar 

  30. Luhr OR, Karlsson M, Thorsteinsson A, Rylander C, Frostell CG. The impact of respiratory variables on mortality in non-ARDS and ARDS patients requiring mechanical ventilation. Intens Care Med. 2000;26:508–17.

    Article  CAS  Google Scholar 

  31. Phillips CR, Chesnutt MS, Smith SM. Extravascular lung water in sepsis-associated acute respiratory distress syndrome: Indexing with predicted body weight improves correlation with severity of illness and survival. Crit Care Med. 2008;36:69–73.

    Article  PubMed  Google Scholar 

  32. Swaroopa D, Bhaskar K, Mahathi T, et al. Association of serum interleukin-6, interleukin-8, and Acute Physiology and Chronic Health Evaluation II score with clinical outcome in patients with acute respiratory distress syndrome. Indian J Crit Care Med. 2016;20:518–25.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Longin E, Schaible T, Demirakca S, et al. Heart rate variability during extracorporeal membrane oxygenation and recovery in severe neonatal disease. Early Hum Dev. 2006;82:135–42.

    Article  PubMed  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  35. Chen WL, Chen GY, Kuo CD. Hypoxemia and autonomic nervous dysfunction in patients with chronic obstructive pulmonary disease. Respir Med. 2006;100:1547–53.

    Article  PubMed  Google Scholar 

  36. Chen WL, Chen JH, Huang CC, et al. Heart rate variability measures as predictors of in-hospital mortality in ED patients with sepsis. Am J Emerg Med. 2008;26:395–401.

    Article  PubMed  Google Scholar 

  37. Borovikova LV, Ivanova S, Zhang M, et al. Vagus nerve stimulation attenuates the systemic inflammatory response to endotoxin. Nature. 2000;405:458–62.

    Article  CAS  PubMed  Google Scholar 

  38. Win NN, Fukayama H, Kohase H, et al. The different effects of intravenous propofol and midazolam sedation on hemodynamic and heart rate variability. Anesth Analg. 2005;101:97–102.

    Article  CAS  PubMed  Google Scholar 

  39. Kasaoka S, Nakahara T, Kawamura Y, et al. Real-time monitoring of heart rate variability in critically ill patients. J Crit Care. 2010;25:313–6.

    Article  PubMed  Google Scholar 

  40. Van de Louw A, Medigue C, Papelier Y, et al. Positive end-expiratory pressure may alter breathing cardiovascular variability and baroreflex gain in mechanically ventilated patients. Respir Res. 2010;11:38.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This study was supported by a grant VGHUST99-P1–09 from the Veterans General Hospitals and University System of Taiwan Joint Research Program, Taiwan, and a grant 96VN-010 from the Taipei Veterans General Hospital-National Taiwan University Hospital Joint Research Program, Taipei, Taiwan. The funding bodies did not participate in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

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

Author information

Authors and Affiliations

Authors

Contributions

ICC recruited the patients, collected the data, performed preliminary data analysis and drafted the manuscript. CTK performed detailed statistical analysis and prepared the fig. JK, CHL, JZT and WJK participated in the discussion and improved the manuscript. CDK made substantial contributions to the original idea and design, analyses and interpretation of data as well as revising the manuscript. CDK is the corresponding author and is responsible for the finalization of the manuscript. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Cheng-Deng Kuo.

Ethics declarations

Ethics approval and consent to participate

The study protocol has been approved by the Institutional Review Boards of National Taiwan University Hospital (NTUH200808065R) and Taipei Veterans General Hospital (VGHIRB97–01-02A), and written informed consent was obtained from the next of kin of the patients before their enrollment in the study.

Consent for publication

All included patients or their family members signed the informed consent form to report individual patient data. All authors have confirmed the manuscript and approved the publication of the manuscript.

Competing interests

The authors declare that they have no competing interest.

Publisher’s Note

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

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, IC., Kor, CT., Lin, CH. et al. High-frequency power of heart rate variability can predict the outcome of thoracic surgical patients with acute respiratory distress syndrome on admission to the intensive care unit: a prospective, single-centric, case-controlled study. BMC Anesthesiol 18, 34 (2018). https://doi.org/10.1186/s12871-018-0497-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12871-018-0497-5

Keywords