This clinical study was conducted in compliance with the ethical principles outlined in the Declaration of Helsinki12 and received approval from the Ethics Committee of Bishan Hospital, Chongqing Medical University, prior to study initiation (Approval number: KYLL2022005). The Ethics Committee waived the requirement for written informed consent due to the retrospective nature of the study and the use of anonymized clinical data. The reagents and the equipment used are listed in the Table of Materials.
1. Study design and participants
This retrospective cohort study included 150 patients diagnosed with COPD at Bishan Hospital of Chongqing Medical University between January 2022 and December 2024. Data were extracted from the hospital’s electronic medical record system. Patients were divided into two groups according to the treatment they actually received: the control group (formoterol monotherapy, n = 75) and the observation group (formoterol combined with budesonide, n = 75). The study flow is illustrated in Figure 1. Data acquisition and analysis were performed by researchers blinded to the clinical treatment assignments.
From the initial pool of eligible patients, 75 patients were randomly selected from each group to achieve balanced sample sizes for statistical comparison. This selection was performed after data extraction, not prospectively. All selected patients had complete follow-up records at all time points (10 days, 3 months, 6 months, 12 months) as documented in routine clinical notes and laboratory databases; therefore, no loss to follow-up occurred within the analyzed cohort.
- Inclusion criteria
Patients were included if they met all of the following criteria: (1) aged 50–80 years; (2) confirmed COPD diagnosis based on symptoms, risk exposure history, physical examination, and post-bronchodilator FEV₁/FVC <70%13; (3) prescribed treatment with either LABA alone (formoterol) or LABA/ICS combination (formoterol/budesonide). (4) COPD severity was classified according to the Global Initiative for Chronic Obstructive Lung Disease (GOLD) 2023 criteria, based on post-bronchodilator forced expiratory volume in one second (FEV₁) as a percentage of predicted value (% predicted). All enrolled patients had moderate to very severe COPD (GOLD stages 2-4).
- Exclusion criteria
Exclusion criteria were: (1) prior diagnosis of type 1 or type 2 diabetes mellitus; (2) fasting plasma glucose ≥7.0 mmol/L or HbA1c ≥6.5% at baseline14; (3) other conditions affecting glucose metabolism (Cushing’s syndrome, thyroid dysfunction, pancreatic disease); (4) concomitant non-COPD respiratory or systemic diseases (bronchial asthma, bronchiectasis, interstitial lung disease, lung cancer, myocardial infarction, osteoporosis); (5) recent use of systemic glucocorticoids, or long-term use of diuretics, immunosuppressants, or biologic agents; (6) end-stage COPD, acute exacerbation at enrollment, severe hepatic or renal impairment, or inability to cooperate with follow-up.
2. Data source and treatment groups
Patients were divided into groups based on their documented treatment regimens retrieved from electronic medical records. The control group received formoterol monotherapy (4.5 µg twice daily via dry powder inhaler). The observation group received a fixed-dose combination of formoterol (4.5 µg) and budesonide (160 µg) per inhalation, twice daily. The total daily dose was 320 µg budesonide and 9 µg formoterol. Treatment adherence was determined by reviewing telephone follow-up records (which were documented as part of routine clinical care every 3–6 months), and patients who changed their medication regimen during the follow-up period were excluded from the final analysis (none met this exclusion criterion).
- Outcome measures
Baseline data collected included age, sex, body mass index (BMI), smoking history, disease duration, FEV₁, FEV₁/FVC, CAT score, FPG, HbA1c, ALT, AST, SCr, arterial blood gases (PaO₂, PaCO₂, pH), and BMD.
Primary outcomes were obtained as follows: FEV₁ and FEV₁/FVC measured at admission (baseline), 10 days, 3 months, and 12 months were retrieved from the pulmonary function test database. FPG and HbA1c measured at baseline and 12 months were obtained from the laboratory information system. FPG was measured by the hexokinase method using commercial kits and an automatic biochemical analyzer.
Secondary outcomes were extracted from routine clinical records: (1) CAT score15at baseline, 10 days, 3 months, and 12 months; (2) Serum creatinine (SCr)16 at baseline, 6 months, and 12 months. Blood samples (5 mL fasting venous blood) were centrifuged at 1500 × g for 10 min; SCr was measured at 546 nm using an automatic biochemical analyzer; (3) ALT and AST17 at baseline, 6 months, and 12 months. Serum was separated by centrifugation at 1500 × g for 10 min; ALT and AST were measured at 340 nm on the same analyzer. The ALT/AST ratio was calculated; (4) Arterial blood gases (PaO2, PaCO2, pH)18at baseline and 12 months. Radial artery blood (1.5 mL) was collected and analyzed immediately using a portable blood gas analyzer; (5) Bone mineral density (BMD)19at baseline and 12 months, measured by dual-energy X-ray absorptiometry (DXA); (6) Adverse reactions (palpitation, headache, skeletal muscle spasm, tremor)20 recorded throughout the treatment period. Incidence (%) was calculated as (number of patients with adverse events / total patients in group) × 100. All measurements were extracted from routine clinical records; no scheduled study-specific visits were added.
3. Statistical analysis
Statistical analyses were performed using a standard statistical software package. All tests were two-tailed, and a P-value <0.05 was considered statistically significant.
Normality of continuous data was assessed using the Shapiro–Wilk test. Normally distributed data were expressed as mean ± standard deviation (SD); non-normally distributed data as median with interquartile range (P25–P75).
Between-group comparisons at baseline and each follow-up time point were performed using independent-samples t‑tests for normally distributed variables (FEV₁, FEV₁/FVC, FPG, HbA1c, CAT scores) and Mann–Whitney U tests for non-normally distributed variables (ALT, ALT/AST ratio, SCr, arterial blood gas parameters, BMD). Within-group changes over time were analyzed using paired t‑tests (for normally distributed variables) or Wilcoxon signed-rank tests (for non-normally distributed variables).
To account for repeated measurements (FEV₁, FEV₁/FVC, CAT scores assessed at four time points; ALT/AST ratio, SCr at three time points), a linear mixed-effects model with unstructured covariance was additionally applied, with treatment group, time, and group‑by‑time interaction as fixed effects and patient as a random intercept. This approach addresses the correlation of repeated observations within the same patient.
All 150 patients had complete records for the scheduled follow-up time points (10 days, 3 months, 6 months, and 12 months) as extracted from routine clinical notes and laboratory databases. Patients with incomplete follow-up records were excluded prior to analysis; no patient met this exclusion criterion. Thus, no missing data were present for any of the reported outcomes.
Effect sizes for between-group comparisons are reported as Cohen’s d for t‑tests and as rank-biserial correlation for Mann–Whitney U tests.