Introduction

Acute respiratory failure (ARF) is one of the most common complications in critically ill patients admitted to the intensive care unit (ICU) and usually leads to high mortality1. Even with mechanical ventilation and extracorporeal membrane oxygenation (ECMO), approximately 36% of patients died consequently during hospitalization2. Although there are some scoring systems used to predict the clinical outcomes of patients in ICU, no one is specifically available for patients with AFR. Several studies have demonstrated that hypoproteinemia, malnutrition, predecease functional status, and Acute Physiology and Chronic Health Evaluation II (APACHE II) score3,4 are related to the need for mechanical ventilation in ARF, but most of these studies were conducted in patients with acute exacerbation of chronic obstructive pulmonary disease (AECOPD) who made up only a small proportion of patients with ARF. Delayed recognition of clinical deterioration in AFR is common, which leads to poor clinical outcomes and increased in-hospital mortality5. Therefore, it is necessary to find useful indicators to predict the clinical outcomes of ARF.

Red blood cell distribution width (RDW) is a parameter of red blood cell volume heterogeneity obtained by a standard complete blood count6. It is commonly used to identify different types of anemia and reactive bone marrow states and it is also a simple, widely used, and inexpensive parameter. In recent years, some clinical studies have shown that there is a correlation between RDW and some acute diseases, such as brain infarction, sepsis, infective endocarditis and diabetic ketoacidosis7,8,9,10,11. In addition, several studies have reported the predictive effect of RDW on lung disease12,13,14,15,16, while the predictive effect of RDW on acute respiratory failure has not been reported. Therefore, our objective is to study whether RDW is useful to predict clinical outcomes of patients with ARF.

Results

Characteristics of patients

There were 2999 patients with ARF included in the study, with 1654 (55.15%) patients died within 3 years. We divided all the patients into three groups, according to the tertiles of RDW: the low group (RDW < 14.1%), the middle group (14.1% ≤ RDW < 15.7%) and the high group (RDW ≥ 15.7%). Patients with high RDW were more likely to have high APS III score, OASIS score and more comorbidities. They also had a high 1-year and 3-year mortality. The characteristics and hematologic laboratory data of the study participants were displayed in Table 1.

Table 1 Characteristics of the enrolled patients according to RDW.

Association of RDW levels and clinical outcomes

The univariate Cox regression analysis (Table S2) indicated that age, gender, APS III score, OASIS score, qSOFA score, SIRS score, SID30, mechanical ventilation, congestive heart failure, chronic renal failure, liver disease, and metastatic cancer were significantly associated with 3‑year all-cause mortality. The above variables also were associated with 1-year mortality (Table S3).

Figure 1 presents the association between baseline RDW levels and risk of 1-year and 3-year mortality (log HR). Taking RDW on admission as a continuous variable, there was a significant positive association with the clinical outcomes in different multivariable models. We further evaluated this finding by the multivariable Cox regression analysis, which was shown in Table 2. As a continuous variable, an SD increase in RDW levels was associated with a 10% higher risk of 3-year mortality in the adjusted model (HR 1.10, 95% CI 1.07, 1.12, P < 0.0001), and there was a similar trend for 1-year mortality (HR 1.10, 95% CI 1.07, 1.12, P < 0.0001). When RDW was assessed as tertiles, there was a significant high risk of 3-year mortality (HR 1.73, 95% CI 1.52, 1.98, P < 0.0001) and 1-year mortality (HR 1.74, 95% CI 1.51, 2.00, P < 0.0001) in patients in high group compared with patients in low group.

Figure 1
figure 1

Association between RDW and clinical outcomes for patients with ARF in different multivariable models. (A,C,E) association between RDW and 1-year mortality for patients with ARF. (B,D,F) association between RDW and 3-year mortality for patients with ARF. Adjusted I for age, gender and ethnicity; Adjusted II for age, gender, ethnicity, liver disease, metastatic cancer, congestive heart failure, renal failure, APS III score, SID30, OASIS, qSOFA, SIRS and mechanical ventilation. HR indicates hazard risk, ARF acute respiratory failure, RDW red cell distribution width, SID30 Elixhauser Comorbidity Index, APS III acute physiology score III, OASIS Oxford acute severity of illness score, qSOFA quick sequential organ failure assessment score, SIRS Systemic inflammatory response syndrome.

Table 2 Association of RDW with 1-year mortality and 3-year mortality.

Subgroup analyses were performed according to age, sex, ethnicity, pulmonary circulatory disease, chronic pulmonary disease, renal failure, liver disease, metastatic cancer, congestive heart failure and mechanical ventilation for the primary and secondary outcomes. The association between RDW on admission and the clinical outcomes was similar for all strata (Tables 3, S4).

Table 3 Subgroup analysis of the associations between RDW and 3-year all-cause mortality by multivariable Cox regression.

ROC curve analysis

ROC curves of 3-year mortality and 1-year mortality generated using the indicated variables (RDW, APS III, OASIS, qSOFA and SIRS) are plotted in Fig. 2A,B. The AUC of 3-year all-cause mortality for RDW was 0.651 (95% CI 0.631, 0.670), which was significantly higher than qSOFA and SIRS scores (P < 0.0001, Fig. 2B). And for 1-year mortality, the AUC was 0.652 (95% CI 0.632, 0.671), which is also higher than qSOFA and SIRS scores. According to the CO2 levels in arterial blood, respiratory failure is divided into hypoxic respiratory failure and hypercarbic respiratory failure. We further analyzed the long-term predictive value of RDW for different types of respiratory failure, and the results were consistent (Fig. 2C–F). The AUC and P value of RDW and different scoring systems for clinical outcomes were shown in Table 4.

Figure 2
figure 2

ROC analyses of predictors of red cell distribution width (RDW) for 3-year and 1-year mortality in critical ill patients, which were compared with APS III, OASIS, qSOFA and SIRS. (A,B) ROC for all patients; (C,D) ROC for patients with hypoxic respiratory failure; (E,F) ROC for patients with hypercarbic respiratory failure. APS III acute physiology score III, OASIS Oxford acute severity of illness score, qSOFA quick sequential organ failure assessment score, SIRS Systemic inflammatory response syndrome, ROC receiver operating curve.

Table 4 ROC of RDW and severity scales for different types of respiratory failure.

Survival status of patients with different admission RDW level

The K–M survival curves showed that patients in the high RDW group had significantly shorter survival time and higher mortality than patients in the middle and low RDW group (log-rank test: P < 0.001), which was presented in Fig. 3A. And for different types of respiratory failure, we obtain the same results (Fig. 3B,C).

Figure 3
figure 3

Kaplan–Meier (K-M) survival curves for 3-year all-cause mortalities. (A) 3-year all-cause mortalities for all patients; (B) 3-year all-cause mortalities for patients with hypoxic respiratory failure; (C) 3-year all-cause mortalities for patients with hypercarbic respiratory failure.

Discussion

Our results demonstrated that there is a significant positive correlation between RDW on admission and long-term mortality of patients with ARF. High RDW was associated with high 3-year and 1-year mortality, and this correlation was also demonstrated in all subgroup analyses. We further analyzed the relationship between the number of reticulocyte and the RDW grouping, and the results indicated that the count of reticulocyte increased gradually according to the three groups (low, medium and high) of RDW patients, with a significant difference (Table S5). The results suggest a fact that hypoxemia lead to the increase of red blood cells, thus affecting red blood cell size.

RDW represents the size variation of circulating red blood cells. Any physiologic process that affects the morphology of red blood cells and causes the early release of young cells into circulation can lead to the increase of RDW. First, ARF can lead to severe hypoxemia which induces EPO release that can reach a level of several 100-fold17. After the onset of ARF, hypoxia can lead to the production of hypoxia-inducible factors (HIFs) which may induced cell-type specific gene expression changes, thus, increasing the production of EPO in the renal and liver under the action of hypoxia-induced transcription factor 2 (HIF-2)17. EPO increases regardless of hypoxemia from a variety of causes, including COPD, pulmonary hypertension, interstitial lung disease, or heart failure18,19,20. The EPO not only increases the rate of formation of RBCs, but also increases the volume of RBCs leading to an increase in RDW18,21. Second, Previous studies have shown that RDW is associated with the outcomes of several inflammatory diseases, including acute pancreatitis, sepsis, and chronic heart failure22,23,24,25. Abnormal RDW elevation is closely related to proinflammatory factors24,26. Proinflammatory cytokines can decrease the half-life of red blood cells and change the morphology of red blood cells27. Moreover, inflammation can retard the maturation of red blood cells, causing an upregulation of reticulocytosis and the release of a large number of reticulocytes into the peripheral circulation, leading to the increase of RDW28,29. Third, in patients with mechanical ventilation and acute lung injury, oxidative stress presents to reactive active oxygen free radicals which can curtail the lifespan of red blood cells, thus promoting the release of young cells into the circulation30. Renal failures, hyperglycemia and vitamin D3 deficiency have also been suggested as factors linking RDW to poor clinical outcomes6.

Combined with the results of reticulocyte and previous studies, we proposed the following hypothesis: after the onset of ARF, hypoxemia can promote the synthesis of EPO in the kidney and liver, thereby increasing the rate of formation of RBCs and affecting the morphology of red blood cells. Inflammation and oxidative stress prevent the maturation of red blood cells, resulting in the release of large numbers of reticulocytes into the peripheral circulation and leading to the increase of RDW.

To our knowledge, this is the first study to demonstrate high RDW as an independent predictor for the long-term outcomes of patients with ARF. Previous studies have shown that high RDW is associated with a lower PaO2/FiO2 rate, a lower ventilation-free day and a greater need for ventilation in acute respiratory failure. Tiffany et al.28 found that in a cohort of 637 patients with ARF, participants with high RDW at admission had a 32% reduction in total length of 30-day period without ventilator compared with the patients with RDW in normal range (RR: 0.68; 95% CI 0.55–0.83, P < 0.001). Tom et al.14 showed that after adjusting for anemia, age, and disease severity (PRISM III), high RDW was an independent risk factor for mechanical ventilation of children (OR 2.6; 95% CI 1.4–4.9; P = 0.004). Danny et al.31 constructed survival curves after adjusting for Charlson comorbidity index and reported that 60-day readmission for AECOPD was higher in patients with high RDW than patients with normal RDW (P = 0.0038). Moreover, results from a retrospective cohort study of MIMIC-III data including 404 eligible ARDS patients suggest that RDW ≥ 14.5% was an independent predictor of 30-day (OR 1.91, 95% CI 1.08, 3.39) and 90-day mortality (OR 2.56, 95% CI 1.50, 4.37)13. These findings are similar to our results and high RDW is associated with poor prognosis in patients with respiratory failure.

In addition, we compared the value of RDW and other severity scale scores to predict long-term mortality. The results showed that, the AUC of RDW has a significant higher predictive value than the qSOFA and SIRS scores. In our study, only APS III had predictive accuracy comparable to RDW on long-term clinical outcomes for ARF. However, the clinical application of APS III is very limited, which requires more patients’ data to calculate, including vital signs, blood gas analysis, laboratory test results, etc. On the contrary, RDW is a routine parameter of blood routine. It can be obtained easily without any additional cost, which may make it more convenient for clinical use to early detect the patients with poor long-term prognosis, and early to take different treatment and nursing measures.

The advantage of the present study is that we analyzed the long-term prognosis of ARF which is difficult to be observed in other studies, because we had a complete follow-up database and a large sample size. However, our study has some limitations. First, as a retrospective cohort study, it is impossible to adjust for all confounders. We have adjusted for known confounders as much as possible, but there are still some unmeasured variables that may affect our results. Due to several variables were not recorded in MIMIC-III, we lacked some indicators about inflammatory response (such as C-reactive protein and interleukin 6) which may affect our results. Therefore, prospective studies are needed to confirm these results. Second, we only measured RDW levels of patients on admission; we did not investigate the trends, which could reveal more information. Third, this is a single-center study, and the results of the study should also be interpreted with caution when implicating in other populations and areas. Fourth, considering the effects of blood transfusion and anemia on RDW, patients with anemia and received blood transfusion were excluded according to the ICD9 code of MIMIC database. Therefore, our conclusions could not interpret these patients.

Conclusion

In conclusion, our data suggested an association between the RDW and survival time of 3-year follow-up, particularly a high RDW on admission was associated with an increased risk of long-term mortality in patients with ARF. RDW may provide an alternative indicator to predict the prognosis and disease progression and more it is easy to get. There are several theoretical mechanisms to explain the increase of RDW, and further studies are needed to validate the relationship between EPO, pro-inflammatory factors, reticulocytes, and RDW in the serum of patients with ARF.

Patients and methods

Data sources

This is a retrospective cohort study with the data of patients obtained from the Medical Information Mart for Intensive Care III (MIMIC-III). We gained access to the database by taking an online course at the National Institutes of Health and passing the 'protection Human Research Participants’ exam (no. 6182750). MIMIC III is a single-center and freely accessible database which contains 53,423 adult patients (aged over 16) during 2001 to 2012 in the ICU of the Beth Israel Deaconess Medical Center in Boston32. The establishment of the MIMIC III database was approved by the institutional review board of Beth Israel deacons Medical Center and Massachusetts Institute of Technology. Informed consent was not required because all protected health information has been de-identified.

Patients

All the patients who were over the age of 18 and hospitalized in ICU for more than 48 h were included in the study. Patients who did not have an RDW data within 24 h after admission were excluded. Patients with anemia or received blood transfusion were also excluded. If a patient was admitted repeatedly during the study period, we used only the record of his first hospital admission (Fig. 4).

Figure 4
figure 4

Flow chart of the current study. ARF acute respiratory failure, RDW red cell distribution width, ICU intensive care unit.

Variables and outcome measures

The following variables were extracted or calculated: age, sex, comorbidities, ethnicity, admission type, 1-year and 3-year mortality, length of ICU and hospital stay, RDW (within 24hour of ICU admission), mechanical ventilation, Acute Physiology Score III (APS III), the Oxford Acute Severity of Illness Score (OASIS), the quick Sequential Organ Failure Assessment score (qSOFA), the Systemic Inflammatory Response Syndrome score (SIRS) and the Elixhauser Comorbidity Index (SID30). The APS III, OASIS, qSOFA and SIRS scores were estimated for all patients within 24 h of ICU admission. The comorbidities were identified and extracted by ICD-9 code in MIMIC-III database. Detailed ICD-9 codes and disease names of each comorbidities are presented in S1. The codes of Structured Query Language which was used to extract the data were obtained from https://github.com/MIT-LCP/mimic-website32. The primary outcome was the 3-year mortality following hospital admission, and the secondary outcome was 1-year mortality.

Statistical analysis

The whole participants were divided into three groups according to the tertiles of RDW. The categorical variables were presented as percentages and the continuous variables were expressed as the mean (SD) or IQR. The baseline characteristics and the scoring systems were compared by chi-square test when they were categorical variables and by Kruskal–Wallis test when they were continuous variables for the three groups. Multivariable Cox regression model and generalized linear models with a logit link were used to test the independent effects of RDW on 3-year and 1-year all-cause mortality with crude and full model. The adjusting variables included age, gender, ethnicity, liver disease, metastatic cancer, congestive heart failure, renal failure, APS III score, ECI (SID30), OASIS, qSOFA, SIRS and mechanical ventilation. We selected those confounders based on their associations with the outcomes of interest or a change in effect estimate of more than 10%.

Additionally, a subgroup analysis was performed to determine whether there were differences of each subgroup in RDW prediction of clinical outcomes. The effects of RDW (as a continuous variable) on the risks of ARF were estimated using COX regression models among subgroups classified by adjusting for age, gender, ethnicity, liver disease, metastatic cancer, congestive heart failure, renal failure, APS III score, SID30, OASIS, qSOFA, SIRS and mechanical ventilation, if not stratified. The Kaplan–Meier (K–M) method and log-rank tests were used to compare the differences in survival rate between each group of patients with RDW at admission. We performed the Receiver operating characteristic (ROC) curves to assess the predictive value of RDW for 3-year and 1-year mortality of patients.

All the data were processed and analyzed by EmpowerStats software (www.empowerstats.com version R.3.4.3) and statistical software package R. A two-tailed P value < 0.05 was considered statistically significant.