Gastroesophageal reflux disease predicts 6-month readmission in acute exacerbation of chronic obstructive pulmonary disease: development and validation of a nomogram
Highlight box
Key findings
• This study developed and internally validated a nomogram incorporating five readily available clinical variables (body mass index, acute exacerbation in previous 1 year, gastroesophageal reflux disease (GERD), Global Initiative for Chronic Obstructive Lung Disease stage, and B-type natriuretic peptide) to predict 6‑month readmission risk in patients with acute exacerbation of chronic obstructive pulmonary disease (AECOPD). The model demonstrated acceptable discrimination (area under the curve: 0.708–0.731) and good calibration.
What is known and what is new?
• Existing chronic obstructive pulmonary disease readmission prediction models often require many predictors or lack modifiable risk factors. Currently, there are few studies that incorporate GERD into predictive models.
• Our nomogram requires only five routine clinical indicators and identifies GERD as a novel, independent, and potentially treatable risk factor for readmission (odds ratio =2.553).
What is the implication, and what should change now?
• Clinicians can use this simple nomogram at the bedside to identify AECOPD patients at high risk for 6‑month readmission. For patients with comorbid GERD, a more stringent post‑discharge follow‑up plan and active GERD management should be considered.
Introduction
Chronic obstructive pulmonary disease (COPD) is the most common respiratory disease, characterized by progressive irreversible decline in lung function and extensive lung damage. According to statistics, there are 392 million COPD patients worldwide, causing over 3 million deaths annually, making it the third leading cause of death globally (1,2). In developing countries such as China, due to the high smoking rate among men and air pollution from early industrialization, the prevalence of COPD is higher than that in developed countries (3). Currently, there are approximately 100 million COPD patients in China, ranking as the third leading cause of death after stroke and ischemic heart disease (4). COPD patients are prone to recurrent acute exacerbations, primarily manifesting as dyspnea, cough, and sputum production (5). Approximately 6.2% of COPD patients experience severe acute exacerbations annually, requiring hospitalization (6). Acute exacerbations of COPD (AECOPD) can accelerate lung function decline, worsen quality of life, and increase mortality (7,8). COPD often coexists with multiple comorbidities, including cardiovascular diseases, diabetes mellitus, anemia, osteoporosis, mood disorders, and gastroesophageal reflux disease (GERD) (9). These comorbidities interact with COPD and exhibit a bidirectional causal relationship, significantly increasing the risk of acute exacerbations, hospitalization rates, and mortality in COPD patients. GERD is a common digestive system disorder defined as the reflux of gastric contents (gastric acid, food, bile, etc.) into the esophagus, causing symptoms and/or complications (10). GERD is one of the common comorbidities of COPD but has often been overlooked in the past (9). Microaspiration resulting from reflux of gastric contents can directly damage the airways and lung parenchyma and increase infection risk, thereby increasing the risk of asthma and pneumonia (11,12). However, studies on the impact of GERD on acute exacerbations of COPD remain limited. This study aimed to analyze the risk factors for readmission within 6 months after discharge in patients with AECOPD, to construct a logistic regression prediction model, and to internally validate the model’s predictive performance. We present this article in accordance with the TRIPOD reporting checklist (available at https://jtd.amegroups.com/article/view/10.21037/jtd-2026-0849/rc).
Methods
Study design and patients
This study consecutively enrolled hospitalized patients with AECOPD admitted to the Department of Respiratory Medicine at Anhui No. 2 Provincial People’s Hospital from May 2023 to March 2025. Index hospitalization was defined as the first COPD-related hospitalization during the study period.
Inclusion criteria: (I) primary diagnosis met the diagnostic criteria for COPD according to the Global Initiative for Chronic Obstructive Lung Disease (GOLD) 2023 report (13), with a primary diagnosis of AECOPD [International Classification of Diseases, Tenth Revision (ICD-10) codes J44.001 or J44.101]; (II) received standardized inpatient treatment (length of hospital stay ≥3 days) to ensure adequate clinical observation and completeness of follow-up data; (III) age ≥60 years, as AECOPD admissions in Anhui No. 2 Provincial People’s Hospital occur predominantly in this age group.
Exclusion criteria: (I) concurrent malignant tumors; (II) concurrent chronic kidney disease or chronic liver disease, because these conditions may independently affect B-type natriuretic peptide (BNP) levels and nutritional markers; (III) stroke with significant sequelae; (IV) severe pulmonary hypertension or severe cardiac insufficiency; (V) other neuromuscular diseases affecting swallowing function; (VI) incomplete data; (VII) death during hospitalization or follow-up.
This study adhered to the ethical principles of the Declaration of Helsinki and its subsequent amendments. The study protocol was approved by the Ethics Committee of Anhui No. 2 Provincial People’s Hospital (approval No. 2024154). All patients provided written informed consent.
Data collection
The following data were collected: age, sex, body mass index (BMI), smoking status (never, previous, current), hypertension, diabetes, coronary heart disease (CHD), acute exacerbation in previous 1 year, GOLD stage, inhaled medications for COPD controller [bronchodilators, inhaled corticosteroids (ICS), anticholinergic drugs], GERD, obstructive sleep apnea-hypopnea syndrome (OSAHS), and laboratory testing parameters including PaO2, PaCO2, neutrophil count, lymphocyte count, hemoglobin, red blood cell distribution width (RDW), platelet count, BNP, C-reactive protein (CRP), albumin (Alb), uric acid, triglycerides, high-density lipoprotein (HDL), and low-density lipoprotein (LDL). Pulmonary function test data, including forced expiratory volume in 1 second (FEV1; mL) and forced vital capacity (FVC; mL), were also collected.
Definitions of the variables
GOLD stage: COPD severity was classified into stages 1 to 4 based on the percentage of predicted FEV1 obtained from pulmonary function testing (14). The classification criteria were as follows: GOLD stage 1 (mild): FEV1 ≥80% predicted; GOLD stage 2 (moderate): FEV1 50–79% predicted; GOLD stage 3 (severe): FEV1 30–49% predicted; GOLD stage 4 (very severe): FEV1 <30% predicted. Higher stages indicate more severe airflow limitation.
GERD: the diagnosis of GERD was assessed using the GERD Questionnaire (GerdQ), with a cumulative score of ≥8 used as the cutoff value for diagnosis (15).
Readmission and follow-up
During the follow-up period of 6 months ±15 days, any readmission was recorded, along with its frequency and cause. Readmission due to AECOPD was defined as the endpoint of this study. Hospitalizations for any reason other than COPD were excluded.
Statistical analysis
All statistical analyses were performed using R software (version 4.5.2) and SPSS 27.0. A two-sided P value <0.05 was considered statistically significant. Continuous variables were first tested for normality using the Shapiro-Wilk test. Normally distributed continuous variables were presented as mean ± standard deviation, and comparisons between groups were performed using the independent samples t-test. Non-normally distributed continuous variables were presented as median with interquartile range, and comparisons between groups were performed using the Mann-Whitney U rank-sum test. Categorical variables were presented as frequencies and percentages (n, %), and comparisons between groups were performed using the Chi-squared test.
The patient cohort was randomly assigned to the training group (n=478) and validation group (n=205) in a 7:3 ratio. The least absolute shrinkage and selection operator (LASSO) regression was used to preliminarily screen predictors of readmission risk. Given the large number of candidate variables included, to avoid overfitting, the optimal regularization parameter λ was determined using 10-fold cross-validation. Two criteria were applied to select predictors with non-zero coefficients: λ_min (the λ value with the minimum cross-validation error) and λ_1se (the λ value corresponding to the most parsimonious model within one standard error of the minimum error). Subsequently, the potential predictors identified by LASSO regression were incorporated into multivariate logistic regression analysis, and a predictive model was established using backward stepwise regression. A nomogram was then constructed to predict the risk of readmission within 6 months after discharge in patients with AECOPD.
Model discrimination was evaluated using the receiver operating characteristic (ROC) curve and the area under the curve (AUC). Model calibration was assessed using the Hosmer-Lemeshow goodness-of-fit test and calibration curves. Calibration was quantitatively assessed using calibration intercept, calibration slope, and Brier score. An intercept of 0 and a slope of 1 indicate perfect calibration; a slope <1 suggests overfitting. The Brier score ranges from 0 to 1, with lower values indicating better accuracy. Decision curve analysis (DCA) was employed to evaluate the clinical net benefit of the nomogram model. The lower x-axis (cost: benefit ratio) expresses the threshold probability as the relative weight of false positives (harm) to true positives (benefit): harm : benefit = 1 : c, and threshold probability = 1/(1 + c). Internal validation was performed using the validation set to assess the generalizability of the model.
Results
Patient enrollment
A total of 878 patients with AECOPD were initially enrolled during the study period. After exclusions, 15 patients with malignant tumors, 35 with chronic kidney disease or chronic liver disease, 39 with stroke and significant sequelae, 45 with severe pulmonary hypertension or severe cardiac insufficiency, 19 with other neuromuscular diseases affecting swallowing function, and 42 with incomplete data or death during hospitalization or follow-up were excluded. Ultimately, 683 patients were included in the final analysis, comprising 478 patients in the training group (140 in the readmission group and 338 in the non-readmission group) and 205 patients in the validation group (Figure 1).
Comparison of demographic, clinical, and laboratory characteristics between the training group and validation group
The results showed no statistically significant differences between the two groups in baseline characteristics, including readmission rate (29.29% vs. 30.24%, P=0.87), age (median 77 years in both groups, P=0.69), sex (male proportion 64.23% vs. 60.49%, P=0.40), and BMI (24.40 vs. 24.50 kg/m2, P=0.94). In addition, smoking status, hypertension, diabetes, CHD, history of acute exacerbation in the previous year, GOLD stage, use of inhaled medications, comorbidities (GERD, OSAHS), and all laboratory parameters were similarly distributed between the two groups (all P>0.05). These findings indicate that the training group and validation group were well balanced in terms of demographic, clinical, and laboratory characteristics, satisfying the baseline comparability requirements for subsequent model development and validation (Table 1).
Table 1
| Variable | Total (n=683) | Internal validation group (n=205) | Training group (n=478) | P value |
|---|---|---|---|---|
| Readmissions | 202 (29.58) | 62 (30.24) | 140 (29.29) | 0.87 |
| Age (years) | 77.00 [70.00–82.00] | 78.00 [70.00–83.00] | 77.00 [70.00–82.00] | 0.69 |
| Male | 431 (63.10) | 124 (60.49) | 307 (64.23) | 0.40 |
| BMI (kg/m2) | 24.50 [21.50–26.20] | 24.50 [21.50–26.10] | 24.40 [21.50–26.20] | 0.94 |
| Smoking | ||||
| Never | 241 (35.29) | 74 (36.10) | 167 (34.94) | 0.86 |
| Previous | 383 (56.08) | 112 (54.63) | 271 (56.69) | 0.86 |
| Current | 59 (8.64) | 19 (9.27) | 40 (8.37) | 0.86 |
| Hypertension | 306 (44.80) | 89 (43.41) | 217 (45.40) | 0.69 |
| Diabetes | 172 (25.18) | 52 (25.37) | 120 (25.10) | >0.99 |
| CHD | 213 (31.19) | 59 (28.78) | 154 (32.22) | 0.42 |
| Acute exacerbation in previous 1 year | 191 (27.96) | 58 (28.29) | 133 (27.82) | 0.97 |
| GOLD stage | 0.90 | |||
| Stage 1 | 220 (32.21) | 68 (33.17) | 152 (31.80) | |
| Stage 2 | 253 (37.04) | 75 (36.59) | 178 (37.24) | |
| Stage 3 | 139 (20.35) | 39 (19.02) | 100 (20.92) | |
| Stage 4 | 71 (10.40) | 23 (11.22) | 48 (10.04) | |
| Inhaled medications (COPD controller) | ||||
| Bronchodilator | 193 (28.26) | 60 (29.27) | 133 (27.82) | 0.77 |
| ICS | 188 (27.53) | 57 (27.80) | 131 (27.41) | 0.98 |
| Anticholinergic drugs | 199 (29.14) | 62 (30.24) | 137 (28.66) | 0.74 |
| GERD | 178 (26.06) | 58 (28.29) | 120 (25.10) | 0.43 |
| OSAHS | 172 (25.18) | 51 (24.88) | 121 (25.31) | 0.98 |
| Laboratory testing | ||||
| PaO2 (mmHg) | 88.00 [76.00–91.00] | 88.00 [77.00–91.00] | 88.00 [74.50–91.00] | 0.76 |
| PaCO2 (mmHg) | 41.00 [37.00–45.00] | 42.00 [38.00–46.00] | 41.00 [37.00–45.00] | 0.55 |
| Neutrophil (×109/L) | 5.97 [3.97–8.93] | 6.03 [4.02–8.58] | 5.96 [3.97–8.97] | 0.93 |
| Lymphocyte (×109/L) | 1.01 [0.73–1.54] | 0.99 [0.74–1.44] | 1.03 [0.71–1.55] | 0.74 |
| Hemoglobin (g/L) | 125.82±18.62 | 125.39±19.10 | 126.00±18.43 | 0.69 |
| RDW (%) | 13.99 [13.09–14.72] | 13.95 [13.07–14.63] | 14.00 [13.10–14.74] | 0.45 |
| Platelet (×109/L) | 180.98±31.68 | 180.96±31.92 | 180.99±31.61 | 0.99 |
| BNP (pg/mL) | 167.16 [142.51–190.49] | 170.96 [143.23–192.00] | 165.54 [142.18–189.76] | 0.30 |
| CRP (mg/L) | 15.23±6.06 | 15.53±5.99 | 15.10±6.09 | 0.39 |
| Albumin (g/L) | 36.90 [32.74–40.50] | 36.70 [32.30–40.00] | 36.90 [33.08–40.58] | 0.21 |
| Uric acid (μmol/L) | 295.58±52.94 | 296.85±51.73 | 295.03±53.50 | 0.68 |
| Triglycerides (mmol/L) | 1.80 [1.50–2.19] | 1.80 [1.55–2.16] | 1.80 [1.47–2.19] | 0.71 |
| HDL (mmol/L) | 1.20 [1.11–1.29] | 1.20 [1.10–1.29] | 1.20 [1.11–1.29] | 0.73 |
| LDL (mmol/L) | 2.36 [1.87–2.73] | 2.19 [1.83–2.72] | 2.36 [1.88–2.73] | 0.29 |
Data are presented as median [interquartile range], n (%), or mean ± standard deviation. BMI, body mass index; BNP, B‑type natriuretic peptide; CHD, coronary heart disease; COPD, chronic obstructive pulmonary disease; CRP, C‑reactive protein; GERD, gastroesophageal reflux disease; GOLD, Global Initiative for Chronic Obstructive Lung Disease; HDL, high‑density lipoprotein; ICS, inhaled corticosteroids; LDL, low‑density lipoprotein; OSAHS, obstructive sleep apnea‑hypopnea syndrome; RDW, red blood cell distribution width.
Comparison of demographic, clinical, and laboratory characteristics between the readmission group and non-readmission group in the training group
Compared with the non-readmission group, patients in the readmission group exhibited the following characteristics: significantly lower BMI (22.75 vs. 24.80 kg/m2, P<0.001), significantly higher BNP level (172.44 vs. 162.44, P<0.001), higher proportion of patients with a history of acute exacerbation in the previous year (37.86% vs. 23.67%, P<0.001), and higher prevalence of GERD (39.29% vs. 19.23%, P<0.001). The distribution of GOLD stage differed significantly between the two groups (P=0.01), with a higher proportion of patients with GOLD stages 3–4 in the readmission group (41.43% vs. 26.63%). No significant differences were observed between the two groups in terms of age, sex, smoking status, other comorbidities (hypertension, diabetes, CHD, OSAHS), use of inhaled medications, or other laboratory parameters (all P>0.05) (Table 2).
Table 2
| Variable | Total (n=478) | Non-readmission group (n=338) | Readmission group (n=140) | P value |
|---|---|---|---|---|
| Age (years) | 77.00 [70.00–82.00] | 77.00 [70.00–82.00] | 77.00 [70.00–82.00] | 0.66 |
| Male | 307 (64.23) | 218 (64.50) | 89 (63.57) | 0.93 |
| BMI (kg/m2) | 24.40 [21.50–26.20] | 24.80 [21.80–26.40] | 22.75 [20.70–25.60] | <0.001 |
| Smoking | 0.69 | |||
| Never | 167 (34.94) | 118 (34.91) | 49 (35.00) | |
| Previous | 271 (56.69) | 194 (57.40) | 77 (55.00) | |
| Current | 40 (8.37) | 26 (7.69) | 14 (10.00) | |
| Hypertension | 217 (45.40) | 154 (45.56) | 63 (45.00) | 0.99 |
| Diabetes | 120 (25.10) | 83 (24.56) | 37 (26.43) | 0.75 |
| CHD | 154 (32.22) | 110 (32.54) | 44 (31.43) | 0.89 |
| Acute exacerbation in previous 1 year | 133 (27.82) | 80 (23.67) | 53 (37.86) | <0.001 |
| GOLD stage | 0.01 | |||
| Stage 1 | 152 (31.80) | 119 (35.21) | 33 (23.57) | |
| Stage 2 | 178 (37.24) | 129 (38.17) | 49 (35.00) | |
| Stage 3 | 100 (20.92) | 61 (18.05) | 39 (27.86) | |
| Stage 4 | 48 (10.04) | 29 (8.58) | 19 (13.57) | |
| Inhaled medications (COPD controller) | ||||
| Bronchodilator | 133 (27.82) | 94 (27.81) | 39 (27.86) | >0.99 |
| ICS | 131 (27.41) | 93 (27.51) | 38 (27.14) | >0.99 |
| Anticholinergic drugs | 137 (28.66) | 97 (28.70) | 40 (28.57) | >0.99 |
| GERD | 120 (25.10) | 65 (19.23) | 55 (39.29) | <0.001 |
| OSAHS | 121 (25.31) | 87 (25.74) | 34 (24.29) | 0.82 |
| Laboratory testing | ||||
| PaO2 (mmHg) | 88.00 [74.50–91.00] | 88.00 [72.00–91.00] | 88.00 [76.75–91.00] | 0.83 |
| PaCO2 (mmHg) | 41.00 [37.00–45.00] | 41.00 [37.00–45.00] | 42.00 [38.00–45.25] | >0.99 |
| Neutrophil (×109/L) | 5.96 [3.97–8.97] | 5.97 [3.86–9.24] | 5.95 [4.06–8.22] | 0.72 |
| Lymphocyte (×109/L) | 1.03 [0.71–1.55] | 1.07 [0.73–1.62] | 0.98 [0.67–1.34] | 0.06 |
| Hemoglobin (g/L) | 126.00±18.43 | 125.75±19.28 | 126.59±16.23 | 0.65 |
| RDW (%) | 14.00 [13.10–14.74] | 13.98 [12.97–14.75] | 14.06 [13.45–14.70] | 0.33 |
| Platelet (×109) | 182.00 [161.25–200.75] | 182.00 [159.00–199.00] | 184.50 [167.75–202.00] | 0.22 |
| BNP (pg/mL) | 165.54 [142.18–189.76] | 162.44 [138.50–187.17] | 172.44 [152.30–191.12] | <0.001 |
| CRP (mg/L) | 15.16 [11.13–18.84] | 15.14 [11.25–18.85] | 15.19 [10.85–18.73] | 0.81 |
| Albumin (g/L) | 36.90 [33.08–40.58] | 36.90 [33.38–40.68] | 37.20 [32.95–40.23] | 0.83 |
| Uric acid (μmol/L) | 295.03±53.50 | 295.43±53.67 | 294.07±53.27 | 0.80 |
| Triglycerides (mmol/L) | 1.80 [1.47–2.19] | 1.80 [1.42–2.26] | 1.81 [1.58–2.02] | 0.73 |
| HDL (mmol/L) | 1.20 [1.11–1.29] | 1.19 [1.11–1.29] | 1.21 [1.12–1.30] | 0.19 |
| LDL (mmol/L) | 2.36 [1.88–2.73] | 2.37 [1.91–2.72] | 2.36 [1.85–2.76] | 0.75 |
Data are presented as median [interquartile range], n (%), or mean ± standard deviation. BNP, B-type natriuretic peptide; BMI, body mass index; CHD, coronary heart disease; COPD, chronic obstructive pulmonary disease; CRP, C-reactive protein; GOLD, Global Initiative for Chronic Obstructive Lung Disease; GERD, gastroesophageal reflux disease; HDL, high-density lipoprotein; ICS, inhaled corticosteroids; LDL, low-density lipoprotein; OSAHS, obstructive sleep apnea-hypopnea syndrome; RDW, red blood cell distribution width.
Construction of the prediction model
LASSO regression analysis for variable screening
LASSO regression was used to screen risk factors for readmission in the training group, with 28 clinical and laboratory parameters included as independent variables. The results showed that under the λ_min criterion, GERD (β =0.77), history of acute exacerbation (β =0.35), GOLD stage (β =0.22), HDL (β =0.67), RDW (β =0.07), PaCO2 (β =0.01), BNP (β =0.01), BMI (β =−0.11), and lymphocyte count (β =−0.17) were included in the final model. Under the more conservative λ_1se criterion, only GERD (β =0.50), history of acute exacerbation (β =0.10), GOLD stage (β =0.09), BMI (β =−0.06), and BNP (β =0.01) retained non-zero coefficients. The coefficients of the remaining variables shrank to zero under both λ selection criteria, indicating that they contributed little to the prediction of readmission (Figure 2A,2B).
Model development and visualization
The five variables screened by LASSO regression were incorporated into multivariate logistic regression analysis. The results showed that all five predictors were significantly associated with the risk of readmission in patients with AECOPD (all P<0.05). Among them, GERD [odds ratio (OR) =2.553; P<0.001] and history of acute exacerbation in the previous year (OR =1.629; P=0.03) significantly increased the risk of readmission. The risk also increased with each increment in GOLD stage (OR =1.365; P=0.005) and each unit increase in BNP level (OR =1.008; P=0.006), while BMI (OR =0.875; P<0.001) was the only protective factor (Table 3).
Table 3
| Variable | OR | 95% CI low | 95% CI up | P value |
|---|---|---|---|---|
| BMI | 0.875 | 0.814 | 0.940 | <0.001 |
| Acute exacerbation in previous 1 year | 1.629 | 1.035 | 2.563 | 0.03 |
| GERD | 2.553 | 1.616 | 4.033 | <0.001 |
| GOLD stage | 1.365 | 1.097 | 1.697 | 0.005 |
| BNP | 1.008 | 1.002 | 1.014 | 0.006 |
BMI, body mass index; BNP, B-type natriuretic peptide; CI, confidence interval; GERD, gastroesophageal reflux disease; GOLD, Global Initiative for Chronic Obstructive Lung Disease; OR, odds ratio.
A nomogram predictive model was constructed based on the above five indicators. The total score was calculated by summing the scores corresponding to each factor, and the risk value on the corresponding axis indicated the probability of readmission within 6 months in patients with AECOPD. A higher total score suggested a greater risk of readmission within 6 months after discharge (Figure 3).
Evaluation and validation of the predictive model
The AUC was 0.708 [95% confidence interval (CI): 0.653–0.758] in the training group, indicating good discriminative ability. The Youden index was 0.330, with a sensitivity of 66.4% and a specificity of 66.6%. The AUC in the validation group was 0.731 (95% CI: 0.648–0.784), which was similar to and slightly higher than that in the training group, suggesting that the model did not exhibit overfitting and had good generalizability and stability (Figure 4A,4B).
The Hosmer-Lemeshow goodness-of-fit test yielded a P value of 0.708 in the training group and 0.738 in the validation group (both P>0.05), indicating a good model fit. The calibration curve of the training group showed that the actual curve closely matched the bias-corrected curve, with an intercept close to 0 and a slope close to 1. The Brier score was 0.18, indicating high consistency between the model predictions and observed outcomes, reflecting good calibration. The calibration curve of the validation group also demonstrated a good fit (Figure 5A,5B).
DCA was used to evaluate the clinical utility of the model. The results showed that within a reasonable high-risk threshold range of 0.2–0.6, the nomogram model provided a higher net benefit than the “treat all” or “treat none” strategies.
The cost: benefit ratio displayed on the lower x-axis represents the threshold probability used to calculate the net benefit. For example, applying this model at a threshold probability of 0.2 (cost-benefit ratio of 1:4) implies that physicians consider the cost of “misidentifying a healthy individual as high-risk” to be merely one-fourth of the benefit derived from “correctly identifying a genuine patient”. Consistent trends were observed in both the training group and the validation group (Figure 6A,6B). This finding indicates that the model has good clinical utility and can serve as a reference for risk stratification of readmission within 6 months after discharge and for intervention decision-making in patients with AECOPD.
Discussion
This prospective cohort study enrolled 683 hospitalized patients with AECOPD and identified BMI, Acute exacerbation in the previous year, GERD, GOLD stage, and BNP as independent predictors of 6month readmission. The model demonstrated acceptable discrimination and good calibration in both the training group (AUC =0.708) and the validation group (AUC =0.731). Calibration curves showed good agreement between predicted probabilities and observed outcomes (Hosmer-Lemeshow test P>0.05). Nomogram models can quantify, visualize, and graphically present the results of logistic regression analyses. They provide an intuitive display of variable values and present predicted outcomes as continuous probabilities, and have therefore been widely applied in clinical practice (16). DCA showed that within a threshold probability range of 0.2–0.6, our nomogram provided a higher net benefit than the “treat all” or “treat none” strategies. Clinicians can use the cost: benefit ratio on the lower x-axis to intuitively understand the trade-off between the cost of “misclassifying a healthy individual as high-risk” and the benefit of “correctly identifying a genuine patient”. Nevertheless, the AUC values remain below 0.8, suggesting that clinically meaningful improvements in predictive accuracy could be achieved by incorporating additional biomarkers or clinical parameters in future studies.
Our model showed comparable predictive performance to that of Zhu et al. (17), who developed a 1-year COPD readmission prediction model incorporating white blood cell (WBC) count, disease duration >10 years, the number of acute exacerbations in the past year, and concurrent respiratory failure (AUC: 0.719 in the training set and 0.676 in the validation set), but the generalizability of their model was clearly inferior to ours. Goto et al. (18) reported a 30-day readmission prediction model for COPD patients using the XGBoost machine learning approach; their simplified model with 11 predictors achieved an AUC of 0.746, which is higher than that of our model. However, their model incorporated more predictors and had a shorter prediction window. Our model predicts over a longer time window, which involves more confounding factors and greater prediction difficulty, yet requires only five routine clinical indicators (BMI, history of acute exacerbations, GERD, GOLD stage, BNP), making it simpler and more user-friendly. None of these existing models included GERD as a predictor, whereas our study identifies GERD as an important but easily overlooked predictor in COPD readmission prediction.
A study using the U.S. National Inpatient Sample database, which included 7,159,694 adult hospitalizations, demonstrated that the prevalence of GERD was significantly higher in patients with COPD than in those without COPD (27.8% vs. 14.1%) (19). Reduced lung compliance in COPD patients leads to mediastinal structural deformation and esophageal traction, resulting in lower esophageal sphincter dysfunction (20-22). GERD accelerates lung function decline in COPD patients (with an additional annual decrease in FEV1 of 2–5 mL) and exacerbates airway trapping (23). These studies demonstrate the bidirectional interaction between GERD and COPD. Proton pump inhibitor treatment in COPD patients with comorbid GERD significantly reduces the frequency of acute exacerbations without increasing the risk of pneumonia (24). These findings suggest that controlling GERD symptoms can effectively reduce the risk of COPD exacerbations, and that proton pump inhibitors and gastrointestinal prokinetic agents should be actively used in COPD patients with comorbid GERD.
Our study also found that BMI, Acute exacerbation in previous 1 year, GOLD stage, and BNP were independent predictors of 6-month readmission in patients with AECOPD. Low BMI often reflects malnutrition, which not only impairs respiratory muscle strength and reduces sputum clearance but also leads to overall decreased immune function, making patients more susceptible to infections or other triggers leading to readmission (25,26). Acute exacerbation in previous 1 year is a risk factor for readmission in COPD patients. Previous studies have shown that a history of exacerbations increases the likelihood of subsequent exacerbations and shortens the interval between them (27,28). A large cohort study based on a Chinese population confirmed that patients in GOLD group A1 (one prior exacerbation) had a 1.704-fold higher risk of hospitalization compared with those without a history of exacerbations (7). GOLD stage is an important reference for assessing the severity of COPD; higher GOLD stages are associated with increased rates of readmission and mortality (29). In COPD patients, even in the absence of a history of heart failure, plasma concentrations of BNP and N-terminal proBNP (NT-proBNP) increase proportionally with the severity of right ventricular diastolic dysfunction (30,31). BNP testing in AECOPD patients can increase the detection rate of newly diagnosed or previously unrecognized heart failure by approximately 20% (31).
Limitation
This study has several limitations. First, GERD was assessed using the GerdQ questionnaire rather than gold-standard methods such as pH monitoring or endoscopy, which may introduce subjective bias. Second, this was a single-center study with a limited sample size; multicenter external validation with larger samples is needed in the future. Third, the model achieved only moderate AUC values (0.708–0.731), indicating room for improvement in predictive accuracy. Fourth, the exclusion criteria (e.g., exclusion of patients with chronic kidney disease, chronic liver disease, severe pulmonary hypertension, or severe cardiac insufficiency) may introduce selection bias and affect the generalizability of the findings. Fifth, we did not record patients’ post-discharge medication adherence (particularly for GERD therapy and COPD preventive medications), which may have an impact on readmission outcomes. In addition, our study only included patients aged ≥60 years. Therefore, our findings may not be generalizable to younger COPD populations. Future multicenter studies should validate our model in younger cohorts.
Conclusions
This study selected five easily accessible clinical variables to construct a risk prediction model for 6-month readmission after discharge in patients with AECOPD, and presented the model as a visualized nomogram, facilitating clinical promotion and application. In particular, GERD was identified as an important predictor. For AECOPD patients with comorbid GERD, a more stringent post-discharge follow-up plan should be considered, and proton pump inhibitors and gastrointestinal prokinetic agents should be actively used to control GERD symptoms in order to reduce the risk of readmission.
Acknowledgments
The authors thank the staff of the Department of Respiratory Medicine, Anhui No. 2 Provincial People’s Hospital, for their assistance with data collection.
Footnote
Reporting Checklist: The authors have completed the TRIPOD reporting checklist. Available at https://jtd.amegroups.com/article/view/10.21037/jtd-2026-0849/rc
Data Sharing Statement: Available at https://jtd.amegroups.com/article/view/10.21037/jtd-2026-0849/dss
Peer Review File: Available at https://jtd.amegroups.com/article/view/10.21037/jtd-2026-0849/prf
Funding: None.
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://jtd.amegroups.com/article/view/10.21037/jtd-2026-0849/coif). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. This study adhered to the ethical principles of the Declaration of Helsinki and its subsequent amendments. The study protocol was approved by the Ethics Committee of Anhui No. 2 Provincial People’s Hospital (approval No. 2024154). All patients provided written informed consent.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- de Oca MM, Perez-Padilla R, Celli B, et al. The global burden of COPD: epidemiology and effect of prevention strategies. Lancet Respir Med 2025;13:709-24. [Crossref] [PubMed]
- Prevalence and attributable health burden of chronic respiratory diseases, 1990-2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet Respir Med 2020;8:585-96.
- Thawanaphong S, Nair P. Contemporary Concise Review 2024: Chronic Obstructive Pulmonary Disease. Respirology 2025;30:574-86. [Crossref] [PubMed]
- Wang C, Xu J, Yang L, et al. Prevalence and risk factors of chronic obstructive pulmonary disease in China (the China Pulmonary Health CPH study): a national cross-sectional study. Lancet 2018;391:1706-17. [Crossref] [PubMed]
- Yang T, Cai B, Cao B, et al. Exacerbation in patients with stable COPD in China: analysis of a prospective, 52-week, nationwide, observational cohort study (REAL). Ther Adv Respir Dis 2023;17:17534666231167353. [Crossref] [PubMed]
- Ronconi G, Dondi L, Dondi L, et al. Missed Opportunities in the Therapeutic Patterns of Chronic Obstructive Pulmonary Disease: A Real-World Italian Study on the Impact of Severe Exacerbation on Mortality and Hospital Re-admission. Clin Drug Investig 2026;46:415-27. [Crossref] [PubMed]
- Lin L, Song Q, Cheng W, et al. Impact of exacerbation history on future risk and treatment outcomes in chronic obstructive pulmonary disease patients: A prospective cohort study based on Global Initiative for Chronic Obstructive Lung Disease (GOLD) A and B classifications. J Glob Health 2024;14:04202. [Crossref] [PubMed]
- Daniels K, Lanes S, Tave A, et al. Risk of Death and Cardiovascular Events Following an Exacerbation of COPD: The EXACOS-CV US Study. Int J Chron Obstruct Pulmon Dis 2024;19:225-41. [Crossref] [PubMed]
- Xiang Y, Luo X. Extrapulmonary Comorbidities Associated with Chronic Obstructive Pulmonary Disease: A Review. Int J Chron Obstruct Pulmon Dis 2024;19:567-78. [Crossref] [PubMed]
- Dunbar KB. Gastroesophageal Reflux Disease. Ann Intern Med 2024;177:ITC113-28. [Crossref] [PubMed]
- Wang X, Wang Y, Bu Y, et al. Association of gastroesophageal reflux disease with the incidence of pulmonary disease. Front Cell Dev Biol 2025;13:1552126. [Crossref] [PubMed]
- Bradley JA, Koop A, Lee AS, et al. Unravelling the Links Between Gastroesophageal Reflux and Lung Disease: New Insights. Dig Dis Sci 2026; Epub ahead of print. [Crossref]
- Agustí A, Celli BR, Criner GJ, et al. Global Initiative for Chronic Obstructive Lung Disease 2023 Report: GOLD Executive Summary. Eur Respir J 2023;61:2300239. [Crossref] [PubMed]
- Global Initiative for Chronic Obstructive Lung Disease. Global Strategy for the Diagnosis, Management and Prevention of COPD (2024 Report) [accessed 2024 Feb 29]. Available online: http://www.goldcopd.org
- Jones R, Junghard O, Dent J, et al. Development of the GerdQ, a tool for the diagnosis and management of gastro-oesophageal reflux disease in primary care. Aliment Pharmacol Ther 2009;30:1030-8. [Crossref] [PubMed]
- Zhu S, Liao X, Jiang W, et al. Development and validation of a nomogram model for predicting clinical pregnancy in endometriosis patients undergoing fresh embryo transfer. BMC Pregnancy Childbirth 2023;23:764. [Crossref] [PubMed]
- Zhu J, Lu Z, Song Q, et al. Development and validation of a nomogram model for predicting one-year unplanned readmission in patients with chronic obstructive pulmonary disease. Eur J Med Res 2025;30:698. [Crossref] [PubMed]
- Goto Y, Niwa D, Shibata S, et al. Advanced machine learning models for prediction of readmission and mortality risks in patients with chronic obstructive pulmonary disease using routine clinical data. Fujita Med J 2025;11:121-8. [Crossref] [PubMed]
- Wang X, Wright Z, Wang J, et al. Elucidating the Link: Chronic Obstructive Pulmonary Disease and the Complex Interplay of Gastroesophageal Reflux Disease and Reflux-Related Complications. Medicina (Kaunas) 2023;59:1270. [Crossref] [PubMed]
- Lee AL, Goldstein RS. Gastroesophageal reflux disease in COPD: links and risks. Int J Chron Obstruct Pulmon Dis 2015;10:1935-49. [Crossref] [PubMed]
- Okwara NC, Chan WW. Sorting out the Relationship Between Esophageal and Pulmonary Disease. Gastroenterol Clin North Am 2021;50:919-34. [Crossref] [PubMed]
- Gadel AA, Mostafa M, Younis A, et al. Esophageal motility pattern and gastro-esophageal reflux in chronic obstructive pulmonary disease. Hepatogastroenterology 2012;59:2498-502. [Crossref] [PubMed]
- Baldomero AK, Wendt CH, Petersen A, et al. Impact of gastroesophageal reflux on longitudinal lung function and quantitative computed tomography in the COPDGene cohort. Respir Res 2020;21:203. [Crossref] [PubMed]
- Kang J, Lee R, Lee SW. Effects of gastroesophageal reflux disease treatment with proton pump inhibitors on the risk of acute exacerbation and pneumonia in patients with COPD. Respir Res 2023;24:75. [Crossref] [PubMed]
- Gea J, Orozco-Levi M, Pascual-Guàrdia S, et al. Biological Mechanisms Involved in Muscle Dysfunction in COPD: An Integrative Damage-Regeneration-Remodeling Framework. Cells 2025;14:1731. [Crossref] [PubMed]
- Hu C, Song B, Liu X, et al. The Association Between COPD, Acute Exacerbations of COPD, and Survival in COPD, with Fat-Free Body Mass Index: A Systematic Review and Meta-Analysis. Int J Chron Obstruct Pulmon Dis 2025;20:2589-600. [Crossref] [PubMed]
- Vogelmeier CF, Diesing J, Kossack N, et al. COPD Exacerbation History and Impact on Future Exacerbations - 8-Year Retrospective Observational Database Cohort Study from Germany. Int J Chron Obstruct Pulmon Dis 2021;16:2407-17. [Crossref] [PubMed]
- Alqahtani JS, Njoku CM, Bereznicki B, et al. Risk factors for all-cause hospital readmission following exacerbation of COPD: a systematic review and meta-analysis. Eur Respir Rev 2020;29:190166. [Crossref] [PubMed]
- Cheng W, Zhou A, Zeng Y, et al. Prediction of Hospitalization and Mortality in Patients with Chronic Obstructive Pulmonary Disease with the New Global Initiative for Chronic Obstructive Lung Disease 2023 Group Classification: A Prospective Cohort and a Retrospective Analysis. Int J Chron Obstruct Pulmon Dis 2023;18:2341-52. [Crossref] [PubMed]
- Calzetta L, Orlandi A, Page C, et al. Brain natriuretic peptide: Much more than a biomarker. Int J Cardiol 2016;221:1031-8. [Crossref] [PubMed]
- Hawkins NM, Khosla A, Virani SA, et al. B-type natriuretic peptides in chronic obstructive pulmonary disease: a systematic review. BMC Pulm Med 2017;17:11. [Crossref] [PubMed]


