3.1 Clinical characteristics of 619 patients with AECOPD.
A total of 619 patients with AECOPD were included into the retrospective study in accordance with the inclusion and exclusion criteria, and the baseline characteristics of these patients were summarized in Table 1. 55 patients died during hospitalization, reflecting an in-hospital mortality rate of 8.89%. The included patients consisted of 516 (83.4%) male patients and 103 (16.6%) female patients, with a median age of 70 years old. The median COPD history of the patients was 5.28 years, and 308 (49.8%) patients had no smoking history. The median height, weight, and body mass index (BMI) of the patients were 163 cm, 56 kg, and 21.4 kg/m2, respectively. 11 (1.8%), 117 (18.9%), 330 (53.3%), and 161 (26.0%) patients with AECOPD were classified as GOLD stages Ⅰ, Ⅱ, Ⅲ, and Ⅳ, respectively. The median FEV1% predicted, FVC% predicted, and FEV1/FVC were 38.5%, 67.5%, and 58.5%, respectively. 92 (14.9%) patients were evaluated as respiratory failure. 72 (11.6%), 180 (29.1%), 199 (32.1%), and 168 (27.1%) patients with AECOPD had 0, 1, 2, and ≥ 3 comorbidities, respectively. 67 (10.8%), 250 (40.4%), 59 (9.5%), 187 (30.2%), and 42 (6.8%) patients were accompanied by asthma, pneumonia, PHD, hypertension, and diabetes mellitus, respectively. The median PLR, NLR, MLR, BLR, and ELR were 185.8, 4.712, 0.465, 0.022, and 0.087, respectively.
55 (8.89%) and 564 (91.1%) patients were divided into dead and survival groups, respectively. The age, PaCO2, WBC, neutrophil, PLR, NLR, MLR, D-dimer, UA, Scr, BUN, and LDH were significantly higher in the dead group compared to the survival group (all P < 0.05). In contrast, the FEV1% predicted, FVC% predicted, FEV1/FVC, PH, PaO2, lymphocyte, eosinophil, basophil, ELR, and albumin were significantly lower in the dead group compared with the survival group (all P < 0.05). In addition, the smoking history, GOLD stages, number of comorbidities, asthma, pneumonia, and PHD were related to in-hospital mortalities (all P < 0.05). The respiratory failure rate of the patients was marginally elevated in the dead group compared to the survival group (P = 0.055, Table 1).
3.2 The levels of systemic inflammatory factors in COPD patients during stable and exacerbation periods and healthy volunteers.
For the survival patients in hospital (n = 564), FEV1% predicted (P = 2.69e-88), FVC% predicted (P = 1.57e-81), and FEV1/FVC (P = 4.73e-18) were dramatically elevated in the stable period compared to the exacerbation period. In the stable period, PaO2 was significantly increased (P = 4.64e-38) while PaCO2 was significantly decreased (P = 5.02e-19), as compared to the exacerbation period (Table S2). Among the survival patients, 410 COPD patients in stable condition performed blood routine examinations. Blood cell counts, Hb, PLR, NLR, MLR, BLR, and ELR in COPD patients between stable and exacerbation conditions were compared in Table 2. PLR (P = 0.015), NLR (P = 5.11e-13), MLR (P = 0.007), and BLR (P = 7.37e-16) were significantly decreased in stable condition compared to exacerbation condition, while ELR was significantly increased in stable condition compared with exacerbation condition (P = 0.018). In addition, demographic characteristics, blood cell counts, Hb, PLR, NLR, MLR, BLR, and ELR in patients with stable COPD and healthy volunteers were also summarized in Table 2. PLR (P = 3.69e-38), NLR (P = 1.00e-69), MLR (P = 2.76e-66), BLR (P = 4.48e-25) and ELR (P = 1.18e-14) were dramatically decreased in the healthy volunteers compared to the patients with stable COPD.
3.3 Relationships between systemic inflammatory factors and airflow limitation in patients with AECOPD.
Kruskal-Wallis test was used evaluate the differences in PLR, NLR, MLR, BLR, and ELR among AECOPD patients with stages Ⅰ, Ⅱ, Ⅲ, and Ⅳ. As a result, there were significant differences in PLR (P = 0.003), NLR (P = 3 × 10-6), MLR (P = 1.4 × 10-5), and ELR (P = 0.025) among AECOPD patients with stages Ⅰ ~ Ⅳ. The median (25th–75th centile) PLR values were 136.3 (125.4 – 171.1), 171.6 (117.4 – 239.0), 184.3 (134.9 – 275.8), and 215.5 (138.4 – 323.1) in AECOPD patients with stages Ⅰ, Ⅱ, Ⅲ, and Ⅳ, respectively. The patients classified as GOLD stage Ⅳ had significantly higher PLR than the patients classified as GOLD stage Ⅱ (Bonferroni’s P = 0.009, Fig 1A). The median (25th–75th centile) NLR values were 3.405 (1.319 – 8.685), 3.631 (2.141 – 6.326), 4.635 (2.827 – 8.692), and 5.870 (3.465 – 12.457) in AECOPD patients with GOLD stages Ⅰ, Ⅱ, Ⅲ, and Ⅳ, respectively. NLR were increased in patients with more severe airflow limitation (Fig 1B). The corresponding median (25th–75th centile) MLR values were 0.281 (0.232 – 1.023), 0.362 (0.262 – 0.563), 0.474 (0.322 – 0.820), and 0.560 (0.390 – 0.858), respectively. The patients classified as GOLD stage Ⅱ had significantly lower MLR compared to the patients classified as GOLD stage Ⅲ (Bonferroni’s P = 6.5 × 10-4) as well as the patients classified as GOLD stage Ⅳ (Bonferroni’s P = 7 × 10-6, Fig 1C). For ELR, the median (25th–75th centile) values were 0.043 (0.034 – 0.083), 0.105 (0.041 – 0.212), 0.090 (0.036 – 0.181), and 0.075 (0.021 – 0.155). The patients with stage Ⅳ had significantly lower ELR than the patients with stage Ⅱ (Bonferroni’s P = 0.037, Fig 1D).
In addition, FEV1% predicted negatively correlated with PLR (r = -0.166, P = 3.2 × 10-5), NLR (r = -0.223, P = 1.99e-8), MLR (r = -0.197, P = 7.69e-7), and BLR (r = -0.088, P = 0.029), whereas ELR presented a positive correlation with FEV1% predicted (r = 0.087, P = 0.030) (Fig S1). With respect to FEV1/FVC, neither BLR (r = -0.073, P = 0.068) nor ELR (r = -0.077, P = 0.055) had significant correlations with FEV1/FVC. PLR (r = -0.129, P = 0.001), NLR (r = -0.117, P = 0.003), and MLR (r = -0.178, P = 8 × 10-6) all presented negative correlations with FEV1/FVC (Fig 2).
3.4 Influences of systemic inflammatory factors on the hospital length of stay and C-reactive protein in patients with AECOPD.
Spearman correlation coefficients were used to evaluate the correlations of systemic inflammatory factors (PLR, NLR, MLR, BLR, and ELR) with hospital length of stay (LOS) as well as CRP due to abnormal distributions of these clinical data. The median (25th–75th centile) hospital LOS was 8 (6 – 11) days. As shown in Fig 3, PLR (r = 0.152, P = 1.5 × 10-4), NLR (r = 0.279, P = 1.52e-12), and MLR (r = 0.262, P = 3.30e-11) were all positively correlated with hospital LOS, while ELR was negatively correlated with hospital LOS (r = -0.117, P = 0.004). BLR had no correlation with hospital LOS (r = 0.035, P = 0.390). The median (25th–75th centile) CRP was 10.3 (2.2 – 45.7) mg/L. As shown in Fig S2, CRP positively correlated with PLR (r = 0.280, P = 2.12e-12), NLR (r = 0.464, P = 1.06e-33), and MLR (r = 0.456, P = 1.88e-32), but negatively correlated with ELR (r = -0.166, P = 4.1 × 10-5). Moreover, BLR had no correlation with CRP (r = 0.039, P = 0.332).
3.5 Receiver operating characteristic (ROC) curves for in-hospital mortality in patients with AECOPD.
ROC curves were plotted to assess the sensitivity, specificity, and accuracy of the systemic inflammatory factors in patients with AECOPD, and the area under the curve (AUC) was used to indicate the predictive ability for in-hospital mortality. The AUCs (95% CI) for PLR, NLR, MLR, BLR, and ELR were 0.650 (0.569 – 0.731), 0.715 (0.646 – 0.785), 0.721 (0.651 – 0.791), 0.540 (0.448 – 0.632), and 0.644 (0.568 – 0.720), respectively (Fig 4). The results indicated that both NLR and MLR had higher predictive value for in-hospital mortality compared to PLR as well as ELR. In addition, the results suggested that MLR had the highest predictive value for in-hospital mortality. The AUC for BLR was close to 0.5, implying that BLR almost had no predictive ability for in-hospital mortality. In addition, the AUCs (95% CI) for GOLD stages, FEV1% predicted, FVC% predicted, FEV1/FVC, and CRP were 0.677 (0.607 – 0.746), 0.729 (0.662 – 0.795), 0.714 (0.641 – 0.787), 0.679 (0.602 – 0.737), and 0.696 (0.638 – 0.753), respectively, suggesting that FEV1% predicted had the highest predictive ability for in-hospital mortality. The AUCs for NLR, MLR, FEV1% predicted, and FVC% predicted all exceeded 0.7, indicating that these clinical characteristics had good predictive ability for in-hospital mortality.
3.6 Association of clinical characteristics with in-hospital mortality in patients with AECOPD.
The age, FEV1% predicted, FVC% predicted, FEV1/FVC, PH, PaO2, PaCO2, WBC, neutrophil, lymphocyte, eosinophil, basophil, PLR, NLR, MLR, ELR, CRP, D-dimer, UA, Scr, BUN, albumin, and LDH as continuous variables had differences between survival and dead patients with AECOPD, and had potential to predict in-hospital mortality (Table 1). Association of the clinical parameters with in-hospital mortality of patients with AECOPD was analyzed using univariate binary logistic regression. As shown in Table 3, these clinical parameters, except for basophil and LDH, were significantly associated with in-hospital mortality in patients with AECOPD. Among these clinical parameters, the age, PaCO2, WBC, neutrophil, PLR, NLR, MLR, CRP, D-dimer, UA, Scr, and BUN were all risk factors for in-hospital mortality in patients with AECOPD (all OR values > 1, Table 3).
In view of the correlations between clinical characteristics (Table S3), multivariate binary logistic regression analysis included the age, smoking history, FEV1% predicted, number of comorbidities, asthma, pneumonia, PHD, PH, PaCO2, PaO2, CRP, D-dimer, UA, Scr, BUN, albumin, WBC, NLR, MLR, and ELR. As shown in Table 4, smoking history, FEV1% predicted, pneumonia, PHD, UA, albumin, and MLR were significant independent predictors for in-hospital mortality in AECOPD patients. Through multivariate analysis, increased ELR was marginally associated with decreased in-hospital mortality (P = 0.070). Among these predictors, smoking history, pneumonia, PHD, UA, and MLR were all risk factors for in-hospital mortality.
3.7 Establishment and validation of a predictive model.
The predictors, including smoking history, FEV1% predicted, pneumonia, PHD, UA, albumin, MLR, and ELR were used to establish a nomogram to predict in-hospital mortality in AECOPD patients (Fig 5). Each predictor reflected a designated score presented on the “Points” line (top line) of the nomogram. The total score of a patient with AECOPD was a summation of the scores for all the predictors. The probability of in-hospital mortality in a patient with AECOPD was estimated on the basis of the total score. The index of concordance (C-index) for the nomogram was 0.850 (95% CI: 0.799 – 0.901), indicating that the nomogram had high predictive ability. Moreover, the nomogram without MLR or ELR was also plotted with a C-index of 0.826 (95% CI: 0.772 – 0.879, Fig S3). The C-index and predictive ability of the nomogram were validated by ROC curves (Fig 6A). The AUC for the nomogram with MLR and ELR was significantly higher than the AUC for nomogram without MLR or ELR (P = 0.015), indicating that the introduction of MLR and ELR significantly improved the predictive ability of the nomogram. Decision curve analysis (DCA) was performed to assess the clinical applicability of the nomogram as we previously described [29]. The nomogram with MLR and ELR added high net benefits within a wide range of threshold probabilities (Fig 6B). Calibration curves were plotted through bootstrap sampling 2000 times. The calibration curves for the nomogram with MLR and ELR and the nomogram without MLR or ELR were close to the reference lines, indicating that the predicted values by both the nomograms were consistent with the actual observed values (Fig 6C and 6D).
A clinical impact curve (CIC) was used to further evaluate the predictive value and clinical applicability of the nomogram. As shown in Fig 7A, the nomogram with MLR and ELR had superior net benefits within a wide range of threshold probabilities. DAC and CIC analysis indicated that the nomogram with MLR and ELR possessed significant predictive value. Moreover, the nomogram without MLR or ELR also had good predictive value (Fig 7B). The net reclassification improvement (NRI) and integrated discriminatory index (IDI) were used to compare the predictive value of the two nomograms. Compared to the nomogram without MLR or ELR, the nomogram with MLR and ELR had a NRI of 0.323 (95% CI: 0.048 – 0.599; P = 0.022) and IDI of 0.048 (95% CI: 0.012 – 0.084; P = 0.009), indicating improved discriminatory performance.