Baseline Characteristics and Follow-up Outcomes
From January 2022 to October 2023, 248 hospitalized patients with coexisting CCS and COPD were initially assessed. After application of the exclusion criteria, 214 patients were enrolled in the study (Figure 1). During the one-year follow-up, 53 patients experienced the composite endpoint, corresponding to 24.8% of the study cohort. The reported endpoint components comprised 7 all-cause deaths (3.3%), 8 non-fatal myocardial infarctions (3.7%), 17 rehospitalizations for heart failure (7.9%), and 21 rehospitalizations for AECOPD (9.8%). All percentages were calculated using the full cohort of 214 patients as the denominator. As specified in the Methods, patients experiencing more than one qualifying endpoint component were counted only once in the primary composite endpoint analysis, using the date of their first qualifying event.

Figure 1: Flowchart of patient screening, eligibility assessment, and follow-up. A total of 248 hospitalized patients with concomitant chronic coronary syndrome (CCS) and chronic obstructive pulmonary disease (COPD) were initially screened. Of these, 34 patients were excluded because of acute coronary syndrome (n = 18), malignancy (n = 2), severe hepatic or renal dysfunction (n = 11), or substantial missing clinical data (n = 3). The final study cohort comprised 214 patients, of whom 206 (96.3%) completed the one-year follow-up and 8 (3.7%) were lost to follow-up. Patients lost to follow-up were censored at their last documented contact, and their available data were retained in the survival analyses. Please click here to view a larger version of this figure.
The event-free group comprised 161 patients. Accordingly, the study population was classified into an event group (n = 53) and an event-free group (n = 161), based on the occurrence of the composite endpoint. Baseline demographic, clinical, laboratory, and medication characteristics of the two groups are presented in Table 1.
Table 1: Comparison of baseline characteristics between event and event-free groups in patients with CCS and COPD. Baseline demographic characteristics, clinical history, laboratory measurements, and discharge medications were compared between patients who experienced the composite endpoint (event group, n = 53) and those who remained event-free (event-free group, n = 161) during the one-year follow-up. Normally distributed continuous variables are presented as the mean ± standard deviation, non-normally distributed continuous variables as the median (interquartile range), and categorical variables as the number (percentage). Between-group comparisons were performed using Student's t-test, the Mann–Whitney U test, or the chi-square (χ2) test, as appropriate. The reported test statistics are t, Z, and χ2. Statistical significance was defined as P < 0.05. CCS, chronic coronary syndrome; COPD, chronic obstructive pulmonary disease; PCI, percutaneous coronary intervention; GOLD, Global Initiative for Chronic Obstructive Lung Disease; LVEF, left ventricular ejection fraction; TC, total cholesterol; TG, triglycerides; LDL-C, low-density lipoprotein cholesterol; hs-CRP, high-sensitivity C-reactive protein; NT-proBNP, N-terminal pro-B-type natriuretic peptide; ACEI, angiotensin-converting enzyme inhibitor; ARB, angiotensin receptor blocker. Please click here to download this file.
Univariate Cox Regression Analysis
Using the one-year composite endpoint as the outcome variable, univariate Cox proportional hazards regression was performed for each prespecified variable, with P < 0.10 as the criterion for entry into the multivariate model. As shown in Table 2, variables significantly associated with the composite endpoint included older age (P < 0.001), smoking history (P = 0.010), prior PCI (P = 0.033), higher GOLD grade (P < 0.001), lower LVEF (P < 0.001), elevated hs-CRP (P = 0.034), and higher log-transformed NT-proBNP (ln NT-proBNP; P < 0.001). In contrast, sex, hypertension history, diabetes history, COPD duration, TC, TG, LDL-C, β-blocker use, statin use, and antiplatelet agent use did not meet the screening threshold (P ≥ 0.10) and were therefore not carried forward into multivariate analysis.
Table 2: Univariate Cox regression analysis for the composite endpoint. Univariate Cox proportional hazards regression was performed to evaluate the associations between individual baseline variables and the occurrence of the composite endpoint during the one-year follow-up. Regression coefficients (β), standard errors (SE), Wald chi-square (χ2) statistics, P values, hazard ratios (HRs), and 95% confidence intervals (CIs) are presented. The natural logarithm of NT-proBNP was used in the regression analysis. Statistical significance was defined as P < 0.05. COPD, chronic obstructive pulmonary disease; PCI, percutaneous coronary intervention; GOLD, Global Initiative for Chronic Obstructive Lung Disease; LVEF, left ventricular ejection fraction; TC, total cholesterol; TG, triglycerides; LDL-C, low-density lipoprotein cholesterol; hs-CRP, high-sensitivity C-reactive protein; NT-proBNP, N-terminal pro-B-type natriuretic peptide. Please click here to download this file.
Multivariate Cox Proportional Hazards Model
Variables with univariate P < 0.10 were entered into a multivariate Cox proportional hazards model using backward stepwise regression, with consideration of clinical relevance. The final model retained five variables: age, smoking history, prior PCI, GOLD grade, and LVEF (Table 3). Although age was significantly associated with the composite endpoint in univariate analysis (P < 0.001), this association was no longer statistically significant after multivariate adjustment (P = 0.090). Therefore, age was retained in the final Cox model but excluded from the risk score, which was constructed using the four predictors that remained statistically significant after adjustment.
| Variable | β | SE | Wald χ2 | P value | HR | 95% CI lower | 95% CI upper |
| Age | 0.036 | 0.021 | 2.872 | 0.09 | 1.037 | 0.994 | 1.081 |
| Smoking history | 0.674 | 0.298 | 5.100 | 0.024 | 1.962 | 1.093 | 3.521 |
| Prior PCI history | 0.702 | 0.284 | 6.120 | 0.013 | 2.017 | 1.157 | 3.517 |
| GOLD grade | 0.398 | 0.158 | 6.336 | 0.012 | 1.488 | 1.092 | 2.028 |
| LVEF | −0.069 | 0.022 | 9.460 | 0.002 | 0.933 | 0.893 | 0.975 |
Table 3: Multivariate Cox proportional hazards regression analysis for the composite endpoint. Multivariate Cox proportional hazards regression was performed to evaluate the independent associations of selected baseline variables with the occurrence of the composite endpoint during the one-year follow-up. Regression coefficients (β), standard errors (SE), Wald chi-square (χ2) statistics, P values, hazard ratios (HRs), and 95% confidence intervals (CIs) are presented. Hazard ratios for continuous variables represent the relative change in hazard associated with a one-unit increase in the corresponding variable, whereas hazard ratios for smoking history and prior PCI history represent comparisons with patients without the respective history. Statistical significance was defined as P < 0.05. PCI, percutaneous coronary intervention; GOLD, Global Initiative for Chronic Obstructive Lung Disease; LVEF, left ventricular ejection fraction.
Smoking history (HR = 1.962, 95% CI: 1.093–3.521, P = 0.024), prior PCI (HR = 2.017, 95% CI: 1.157–3.517, P = 0.013), GOLD grade (HR = 1.488, 95% CI: 1.092–2.028, P = 0.012), and LVEF (HR = 0.933, 95% CI: 0.893–0.975, P = 0.002) remained independently associated with the composite endpoint. With GOLD grades I–IV coded as an ordinal variable from 1 to 4, each one-grade increase was associated with a 48.8% higher hazard of the composite endpoint, assuming a constant log-hazard increment across successive grades. This association pertains to the composite endpoint and should not be interpreted as specific to cardiovascular outcomes. The proportional hazards assumption was satisfied for all covariates in the multivariate model (Grambsch–Therneau global test, P = 0.179; all individual covariates, P > 0.05).
Model Discrimination
The predictive performance of the multivariate Cox regression model was evaluated using discrimination, calibration, and clinical utility measures. The concordance index (C-index) was 0.736, with a bootstrap-corrected estimate of 0.718. The likelihood ratio test yielded χ2 = 46.033 (P < 0.001), indicating a statistically significant overall model. The one-year area under the receiver operating characteristic curve (AUC) was 0.773 (95% CI: 0.692–0.846), with a sensitivity of 70.9% and specificity of 93.8%. The bootstrap-corrected AUC was 0.753. The Brier score was 0.142, and the bootstrap internal validation calibration slope was 0.901. The calibration plot demonstrated agreement between predicted and observed one-year event probabilities (Figure 2A). Decision curve analysis (DCA) indicated that the prediction model provided a positive net benefit compared with the treat-all and treat-none strategies across threshold probabilities of approximately 0.10–0.25 (Figure 2B). These findings suggest potential clinical utility for risk-based decision-making within this threshold range.

Figure 2: Assessment of the predictive model performance. (A) Calibration plot comparing the predicted and observed probabilities of the composite endpoint at one year. The solid line represents the model calibration curve, whereas the dashed diagonal line represents perfect calibration. (B) Decision curve analysis showing the net benefit of the prediction model across a range of threshold probabilities. The blue solid line represents the prediction model, the red dashed line represents the treat-all strategy, and the black dashed line represents the treat-none strategy. Please click here to view a larger version of this figure.
Subgroup Analysis
A clinical risk score was constructed using the regression coefficients of the four predictors that remained statistically significant in the multivariate Cox regression model:
Risk score = 0.674 × smoking history + 0.702 × prior PCI history + 0.398 × GOLD grade − 0.069 × LVEF.
Smoking history and prior PCI history were coded as 0 for absence and 1 for presence. GOLD grades I–IV were coded as 1–4 and entered as an ordinal numeric variable, whereas LVEF was entered as a continuous variable expressed as a percentage.
The cohort was divided into tertiles according to the risk score distribution: low-risk (n = 72; score range, −4.592 to −3.083), intermediate-risk (n = 71; score range, −3.076 to −2.242), and high-risk (n = 71; score range, −2.219 to −0.604). Kaplan–Meier curves demonstrated one-year composite endpoint-free survival rates of 89.04%, 84.29%, and 52.11% in the low-, intermediate-, and high-risk groups, respectively (Figure 3). Composite endpoint-free survival differed significantly among the three groups (log-rank χ2 = 32.873, P < 0.001), with a significant trend toward poorer outcomes across increasing risk strata (trend test χ2 = 27.268, P < 0.001).

Figure 3: Kaplan–Meier curves of one-year event-free survival stratified by risk group. Kaplan–Meier curves showing the probability of remaining free from the composite endpoint among patients classified into low-risk (n = 72), intermediate-risk (n = 71), and high-risk (n = 71) groups based on tertiles of the clinical risk score. The composite endpoint comprised all-cause death, non-fatal myocardial infarction, rehospitalization for heart failure, or rehospitalization for acute exacerbation of COPD. Event-free survival differed significantly among the three risk groups (log-rank P < 0.001; trend P < 0.001). The numbers of patients at risk at the specified follow-up time points are presented below the survival curves. Please click here to view a larger version of this figure.
Data Availability:
The de-identified participant-level data supporting the findings of this study, including the source data underlying the analyses, tables, and figures, are available in Supplementary File 1, provided with this article.
Supplementary File 1: De-identified patient-level clinical data and variable coding for the study cohort. The supplementary Excel file contains individual-level data for 214 patients with concomitant chronic coronary syndrome (CCS) and chronic obstructive pulmonary disease (COPD). The Data worksheet includes demographic characteristics, clinical history, COPD severity, laboratory measurements, discharge medications, follow-up information, and clinical outcomes. The Code worksheet provides numerical coding definitions for selected categorical variables. Clinical outcomes include all-cause mortality, non-fatal myocardial infarction, heart failure, and acute exacerbation of COPD. The composite endpoint and follow-up time are recorded in the corresponding columns. BMI, body mass index; PCI, percutaneous coronary intervention; GOLD, Global Initiative for Chronic Obstructive Lung Disease; LVEF, left ventricular ejection fraction; TC, total cholesterol; TG, triglycerides; LDL-C, low-density lipoprotein cholesterol; hs-CRP, high-sensitivity C-reactive protein; NT-proBNP, N-terminal pro-B-type natriuretic peptide; ln NT-proBNP, natural logarithm of NT-proBNP; HF, heart failure; AECOPD, acute exacerbation of chronic obstructive pulmonary disease. Please click here to download this file.