Comparative cardiovascular safety of LABA/LAMA FDC versus LABA/ICS FDC in patients with chronic obstructive pulmonary disease: a population-based cohort study with a target trial emulation framework

Background Use of combinations of long-acting β2 agonists/long-acting muscarinic antagonists (LABA/LAMA) in patients with chronic obstructive pulmonary disease (COPD) is increasing. Nevertheless, existing evidence on cardiovascular risk associated with LABA/LAMA versus another dual combination, LABA/inhaled corticosteroids (ICS), was limited and discrepant. Aim The present cohort study aimed to examine comparative cardiovascular safety of LABA/LAMA and LABA/ICS with a target trial emulation framework, focusing on dual fixed-dose combination (FDC) therapies. Methods We identified patients with COPD who initiated LABA/LAMA FDC or LABA/ICS FDC from a nationwide Taiwanese database during 2017–2020. The outcome of interest was a hospitalized composite cardiovascular events of acute myocardial infarction, unstable angina, heart failure, cardiac dysrhythmia, and ischemic stroke. Cox regression models were used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for composite and individual cardiovascular events after matching up to five LABA/LAMA FDC initiators to one LABA/ICS FDC initiator using propensity scores (PS). Results Among 75,926 PS-matched patients, use of LABA/LAMA FDC did not show a higher cardiovascular risk compared to use of LABA/ICS FDC, with a HR of 0.89 (95% CI, 0.78–1.01) for the composite events, 0.80 (95% CI, 0.61–1.05) for acute myocardial infarction, 1.48 (95% CI, 0.68–3.25) for unstable angina, 1.00 (95% CI, 0.80–1.24) for congestive heart failure, 0.62 (95% CI, 0.37–1.05) for cardiac dysrhythmia, and 0.82 (95% CI, 0.66–1.02) for ischemic stroke. The results did not vary substantially in several pre-specified sensitivity and subgroup analyses. Conclusion Our findings provide important reassurance about comparative cardiovascular safety of LABA/LAMA FDC treatment among patients with COPD. Supplementary Information The online version contains supplementary material available at 10.1186/s12931-023-02545-9.


Introduction
Chronic obstructive pulmonary disease (COPD) poses a paramount clinical burden worldwide.It is the third and the eighth major cause of mortality globally and in Taiwan, respectively [1,2].Long-acting bronchodilators, including long-acting β 2 agonists (LABA) and longacting muscarinic antagonists (LAMA), are the central maintenance treatment in reducing COPD-related symptoms and exacerbations [3].The combination therapy of LABA and LAMA (LABA/LAMA) is further indicated for patients with a high risk of acute exacerbations or with suboptimal response to LABA or LAMA monotherapy [3].On the other hand, because of the pharmacological actions of LABA on β 2 receptors and LAMA on M 3 receptors [4,5], sympathetic activation and potential cardiovascular risk of LABA and LAMA, especially as combination therapy, deserves paying attention.
There may be potential clinical and methodological issues accounting for the aforementioned discrepant findings.For example, the increased cardiovascular risk associated with LABA/LAMA found in the metaanalysis [6] was mainly driven by three efficacy trials in which patients had to experience exacerbation episodes and have apparent respiratory symptoms [9][10][11].However, patients with clinically significant cardiovascular abnormalities were excluded.Moreover, the trial design of abrupt ICS withdrawal at randomization may exacerbate disease control and even cardiovascular outcomes for patients allocated to the LABA/LAMA treatment [12,13].The above cohort studies addressed the cardiovascular risk of LABA/LAMA in the daily practice; however, none of the studies specifically examined the risk among patients with cardiovascular disease or with longer treatment durations [7,8].
With the availability and the increasing use of LABA/ LAMA fixed-dose combination products (FDC, i.e., ≥ 2 active drugs in a single inhaler) in the market [14,15], it becomes important to comprehensively assess cardiovascular safety associated with LABA/LAMA FDC use in the real-world settings.The present populationbased cohort study aimed to examine comparative cardiovascular safety of LABA/LAMA FDC and LABA/ ICS FDC in patients with COPD.LABA/ICS FDC, rather than LABA or LAMA monotherapy, was selected as an active comparison group since both LABA/LAMA FDC and LABA/ICS FDC are dual combination therapies and tend to be comparable in patient characteristics, which may mitigate potential confounding.

Data source
This present study used data from the Taiwan National Health Insurance Research Database (NHIRD), which included de-identified data of approximately 23 million beneficiaries enrolled in a single-payer national health insurance system [16,17].See Additional file 1: eMethods for detailed data source description.

Cohort study with a target trial emulation framework
To our knowledge, no existing trials aimed to examine cardiovascular safety issues of LABA/LAMA FDC versus LABA/ICS FDC among patients with COPD.Therefore, we specified components of a target trial (i.e., a hypothetical trial) and emulated the trial using Taiwan NHIRD.In another word, we conducted a cohort study with a target trial framework in which ascertainment of study population, exposure, outcomes, and baseline covariates were anchored at the cohort entry date (i.e., T 0 ) and appropriate statistical analyses were applied to enhance causal inference estimation using real-world data [18][19][20][21][22]. See Additional file 1: Table S1 for specification of each component and Fig. 1 for the graphic depiction of the design.Corresponding approaches applied in the present study were also mentioned in our previous work [23,24].

Study population and exposure
Our study population comprised patients with COPD who initiated LABA/LAMA FDC or LABA/ICS FDC from the NHIRD between 2017/01/01 and 2020/12/31.COPD was defined as having ≥ 1 outpatient or inpatient relevant diagnosis code in any diagnosis position (see Additional file 1: Table S2 for codes [25]).Use of LABA/ LAMA FDC or LABA/ICS FDC was defined based on outpatient pharmacy dispensing claims (see Additional file 1: Table S3 for codes).The cohort entry date (i.e., T 0 ) was the date of the first dispensing of a LABA/LAMA FDC or LABA/ICS FDC after a diagnosis of COPD.Further exclusion criteria were presented in Additional file 1: eMethods.
The primary analysis applied an on-treatment approach, which followed patients from cohort entry to outcome occurrence, study FDC treatment discontinuation or change, death, or end of data (2020/12/31), whichever came first.Study FDC treatment discontinuation was defined using a grace period of 60 days between the end date of one dispensing and the start date of the next dispensing, if any; and the discontinuation date was 60 days after the end of the last dispensing.Study FDC treatment change was defined as a dispensing of ICS-containing regimens or another LABA/LAMA FDC for LABA/LAMA FDC initiators, and a dispensing of LAMA-containing regimens or another LABA/ICS FDC for LABA/ICS FDC initiators.

Positive control outcome
Trial and real-world data have shown that use of LABA/ LAMA has a lower risk of pneumonia compared to use of LABA/ICS (relative risk of 0.57 [95% CI, 0.42-0.79]and 0.66 [95% CI, 0.50-0.87],respectively) [36,37].Therefore, we chose pneumonia as a positive control outcome and examined if our design could identify a decreased risk of pneumonia associated with LABA/LAMA FDC.This positive control outcome approach is commonly applied in observational studies for assessing internal validity [38,39].The pneumonia outcome was determined using inpatient diagnosis codes recorded in any positions, with the claims-based algorithm exhibiting a positive predictive value of 88% (Additional file 1: Table S4) [40].

Covariates
We assessed > 80 potential baseline confounders including age at cohort entry, sex, COPD duration defined as the duration from the first recorded date of COPD diagnosis (looking back until 2014/01/01) to the cohort entry date, calendar year of the cohort entry date, as well as comorbidities (e.g., cardiovascular disease), other medication use (e.g., cardiovascular medication), influenza or pneumococcal vaccination, and outpatient visits or hospital admissions for any reasons, for COPD, or for cardiovascular disease that may be associated with both the use of study FDC treatment and the risks of study outcomes.These characteristics were determined using diagnosis files or pharmacy dispensing records derived from outpatient or inpatient claims of the NHIRD within 365 days before cohort entry.Additional file 1: Tables S5  and S6 present detailed covariate information.

Statistical analysis
Using all above predefined, claims-based covariates, we estimated the baseline propensity scores (PS), i.e., the probabilities of initiating LABA/LAMA FDC versus LABA/ICS FDC, with a logistic regression model.Because there were many more LABA/LAMA FDC initiators, up to five LABA/LAMA FDC initiators were matched to one LABA/ICS FDC initiator using a nearest-neighbor algorithm without replacement and with a maximum matching caliper of 0.025 on the PS scale (see Additional file 1: eMethods for detailed variable-ratio matching steps) [41].We computed standardized differences for each covariate to evaluate covariate distribution before and after PS matching, with a value < 0.1 in absolute value indicating balance between treatment groups [42].
The incidence rates and 95% CIs for LABA/LAMA FDC and LABA/ICS FDC initiators were estimated according to a Poisson distribution.The cumulative incidence plots were derived by one minus the Kaplan-Meier estimate (i.e., complement of the Kaplan-Meier survival function).Using conventional Cox proportional hazards models, we estimated HRs and 95% CIs of composite cardiovascular events comparing LABA/LAMA FDC with LABA/ICS FDC before and after PS matching.We also examined the risks for the individual components of the composite outcome.All the analyses after PS matching took the variable-ratio matching into consideration.

Sensitivity analyses
For composite cardiovascular events, we conducted four pre-specified sensitivity analyses to examine the robustness of the results comparing LABA/LAMA FDC with LABA/ICS.First, we used an intention-to-treat approach which continued to follow patients regardless of study FDC treatment discontinuation or change during follow-up and censored patients at the first of outcome occurrence, death, or end of data.Second, to account for the influence of competing risk from overall death, we re-draw cumulative incidence plots and applied the Fine-Gray subdistribution hazard model for subdistribution HR estimation, which remains patients who have previously experienced overall death in the risk sets (Additional file 1: eMethods) [43,44].Third, to mitigate potential unmeasured confounding, we conducted highdimensional PS (hd-PS) estimation which identifies and includes additional 100 empirically-identified, claimsbased covariates in the PS model (Additional file 1: eMethods) [45,46].
Finally, besides controlling for claims-based covariates, we attempted to additionally adjust for several important clinical parameters captured from the Taiwan National Health Insurance (NHI) Laboratory Database and the Taiwan COPD pay-for-performance (P4P) Database, both of which were recently established nationwide databases under series of national health policies in Taiwan (see Additional file 1: eMethods for detailed database description [47,48]).Specifically, we measured the following 11 clinical parameters potentially related to the risks of outcomes from outpatient and inpatient health encounter records in the NHI Laboratory Database or the COPD P4P Database within 365 days before cohort entry, including laboratory test results (eosinophil, c-reactive protein, low-density lipoprotein-cholesterol, glycated hemoglobin, and glomerular filtration rate [GFR] or estimated GFR), lung function test results (predicted postdose forced expiratory volume in one second [FEV 1 ], post-dose FEV 1 /forced vital capacity), COPD Assessment Test (CAT) score, systolic blood pressure, and health behavior (body mass index and smoking status) (Additional file 1: Table S7).Because not all the patients had information on these parameters, we conducted multiple imputation to handle missing data issues (Additional file 1: eMethods) [49][50][51].

Subgroup analyses
For composite cardiovascular events, we also conducted three pre-specified subgroup analyses (1) to evaluate potential effect measure modification by patient characteristic (age, sex, COPD duration, history of hospitalized COPD exacerbations, and history of cardiovascular diseases); (2) to assess whether the association was different comparing individual LABA/LAMA FDC with individual LABA/ICS FDC; (3) and to explore potential treatment duration-response relation for LABA/LAMA FDC (1-90, 91-180, 181-365, and > 365 days).We re-estimated the PS and re-matched patients within each subgroup [41,52,53].See Additional file 1: Table S8 for details of each analysis.Besides comparing the subgroup-specific effect estimates, we applied the Wald test to formally test if the results differ materially across patient subgroups, and a p-value < 0.05 was considered statistically significant.
Before PS matching, LABA/LAMA FDC initiators were slightly older (mean age in years: 70 versus 69) and more likely to be male (male %: 85 versus 65) than LABA/ICS FDC initiators.LABA/LAMA FDC initiators also tended to have a history of lung cancer, receive LABA or LAMA monotherapy, and have more frequent outpatient visits for COPD than LABA/ICS FDC initiators.In contrast, LABA/LAMA FDC initiators were less likely to receive ICS than LABA/ICS FDC initiators.The PS model yielded a c-statistic of 0.790.After 5:1 variable-ratio matching, a total of 75,926 patients were included (76% of the study cohort, 80% of LABA/LAMA FDC initiators, and 71% of LABA/ICS FDC initiators).PS matching achieved balance in all baseline claims-based characteristics (Table 1; Additional file 1: Table S9, Fig. S2).

Composite cardiovascular events
Before PS matching and for composite cardiovascular events, the mean follow-up duration was 256 days for LABA/LAMA FDC initiators and 163 days for LABA/ ICS FDC initiators.The crude incidence rates for LABA/ LAMA FDC and LABA/ICS FDC initiators were 24.97 and 30.74 per 1000 person-years, respectively (Table 2), corresponding to a crude HR of 0.92 (95% CI, 0.83-1.02)(Table 3).After PS matching, the incidence rates were 25.22 and 30.90 per 1000 person-years among LABA/ LAMA FDC and LABA/ICS FDC initiators, respectively (Table 2).Use of LABA/LAMA FDC did not show a higher risk compared to LABA/ICS FDC (HR after PS matching, 0.89; 95% CI, 0.78-1.01).The cumulative incidence plots of the composite cardiovascular events were consistent with above findings (Additional file 1: Fig. S3).

Individual cardiovascular events
The majority of cardiovascular events were congestive heart failure, followed by ischemic stroke and acute myocardial infarction (Table 2).Use of LABA/LAMA FDC was not observed with an evident, increased risk of individual outcomes compared to LABA/ICS FDC although there was a numerically elevated risk for unstable angina.

Findings of the positive control outcome
LABA/LAMA FDC initiators had a lower rate of pneumonia compared to LABA/ICS FDC initiators during follow-up, corresponding to an HR after PS matching of 0.65 (95% CI, 0.58-0.74)(Additional file 1: Table S10).

Findings of sensitivity analyses
The intention-to-treat approach yielded mean follow-up days of 651 days for LABA/LAMA FDC initiators and 676 days for LABA/ICS FDC initiators, corresponding to an HR after PS matching of 0.97 (95% CI, 0.90-1.05)(Table 4; Additional file 1: Table S11).After accounting for competing risk from overall death, the cumulative incidence plots tended to be slightly flatter (Additional file 1: Fig. S4) compared to that of the main analysis.However, the Fine-Gray subdistribution HR after PS matching (0.90 [95% CI, 0.79-1.03],see Table 4) did not change materially when compared to what was generated by the traditional Cox model.The hd-PS estimation also yielded a similar finding (HRs after PS matching, 0.88 [95% CI, 0.77-1.00],see Table 4).
A total of 83,025 patients (83% of the eligible patients) had at least one clinical parameter captured from either the NHI Laboratory Database or the COPD P4P Database Information from the NHI Laboratory Database tended to be more comprehensive (approximately 82% of patients having at least one laboratory test result).Information from the COPD P4P Database; however, was less complete (approximately 8% of patients having data of lung function test, respiratory symptoms, blood pressure, or health behavior) (Additional file 1: Table S12).When we applied multiple imputation to incorporate these additional parameters in the PS-matched analysis, the risk estimate did not change apparently from those in the main analysis (HRs of 0.73 [95% CI, 0.31-1.73]when only including data from the NHI Laboratory Database and 0.89 [95% CI, 0.77-1.03]when simultaneously including data from the NHI Laboratory Database and the COPD P4P Database) (Table 4).Distribution of these parameters between FDC treatment groups became balanced as well (Additional file 1: Table S13).

Findings of subgroup analyses
Overall, the risk of composite cardiovascular events associated with LABA/LAMA FDC did not change materially by patient characteristic (Figure 2; Additional file 1: Table S14), individual LABA/LAMA FDC and LABA/ ICS FDC (Additional file 1: Table S15), and treatment duration (Figure 2; Additional file 1: Table S16), although the risk seemed to be lower among some subgroups.
The p-values for test of homogeneity across patient subgroups were > 0.05.

Discussion
This nationwide cohort study with >75,000 eligible patients did not observe an increased cardiovascular risk associated with LABA/LAMA FDC compared to LABA/ ICS FDC.The sensitivity analyses with a longer followup duration, accounting for competing risk from overall death, and addressing potential unmeasured confounding yielded similar findings.The subgroup analyses, stratified by patient characteristic, comparing individual LABA/LAMA FDC to LABA/ICS FDC, and examining different treatment duration, also did not show evident variations in risk estimates across each patient subgroups.Our results provide important reassurance about the comparative safety of dual bronchodilator treatment strategies in COPD in the real-world settings.

Comparison with existing clinical trials
One The discrepant findings generated from these trial data and our study should be elaborated carefully in terms of differences in patient characteristics and methodological concerns across data source.Specifically, given the inclusion and exclusion criteria, patients enrolled in these efficacy trials tended to have more severe COPD disease but had no clinically significant cardiovascular abnormalities.As demonstrated in Additional file 1: Table S17, patients in these trials had and nearly 56-81% of enrollees had baseline ICS use and some may have prior asthma history.In terms of methodological perspectives, the trial design of forced ICS withdrawal at randomization may therefore exacerbate disease control for those who have benefited from ICS treatment but were allocated to the LABA/LAMA group [12,13], potentially leading to subsequent cardiovascular consequences.Moreover, the primary endpoint of the three trials was an annual rate of COPD exacerbations rather than cardiovascular outcomes.This may raise concerns about adjudication and potential misclassification of cardiovascular events in the original trials and in the meta-analysis [54].This may also result in lack of detailed data on time to cardiovascular events for HR estimation.The present Taiwanese cohort study included patients with COPD who initiated LABA/LAMA FDC or LABA/ ICS FDC in the daily practice, who tended to have mild COPD disease but may have major underlying cardiovascular disease compared to enrollees in aforementioned trials.For example, our study patients had mean COPD duration of 2 years (measured in the database) and mean predicted FEV 1 of 63% (mainly based on imputed data); 14% of patients had at least one severe exacerbation episode; the mean CAT score was 13 (mainly based on imputed data); and only 10% of patients ever used ICS at baseline (Additional file 1: Table S17).To mitigate misclassification of cardiovascular outcomes, we ascertained cardiovascular events according to validated algorithms.Our on-treatment and intention-to-treat follow-up approaches yielded similar results, which did not show a higher cardiovascular risk associated with LABA/LAMA FDC use.Our findings highlight that it deserves more caution when directly applying trial data for safety evaluation in the real-world settings.
Using the Taiwan NHIRD with more updated data (2017-2020, n = 75,926), the present cohort study yielded a larger sample size and included more FDC products (three LABA/LAMA FDC [IND/GLY, VIL/UME, OLO/ TIO] and five LABA/ICS FDC [FOR/BEC, FOR/BUD, FOR/FLU, SAL/FLU, VIL/FLU]).Extending findings from prior cohort studies, we comprehensively presented the risk of various cardiovascular outcomes and the risk across patients with different characteristics (such as patients with hospitalized COPD exacerbations or with cardiovascular disease at baseline), with different comparisons of LABA/LAMA FDC versus LABA/ICS FDC, and with longer treatment durations (such as > 365 days).The results provide more informative messages and facilitate physicians for treatment decision in different clinical scenarios.

Strengths and limitation of the present study
Our study has some notable strengths.First, to our knowledge, our study is the largest study which provides real-world evidence on cardiovascular safety of LABA/ LAMA FDC versus LABA/ICS FDC.Second, under a target trial emulation framework, we defined study population, treatment initiation, outcome occurrence, and baseline covariates anchored at the cohort entry date (i.e., T 0 ).This prevented a common shortcoming of "looking forward the future" and accompanying bias in the realworld studies [18][19][20][21][22]. Third, in alignment with the emulation framework, our main analysis extracted diagnosis and medication information from the Taiwan NHIRD derived from a single-payer health insurance system.We used validated claim-based algorithms to determine cardiovascular and pneumonia outcomes.We also replicated the known association between LABA/ICS FDC and pneumonia.To mitigate potential confounding by patient characteristic, we applied the state-of-the-art PS matching analysis accounting for > 80 pre-specified claims-based covariates that may be associated with the exposures and the outcomes.Collectively, these lend support to our study validity.Fourth, we further improved confounding control using the sophisticated hd-PS algorithm with additional empirical claims-based covariates and the multiple imputation approach with additional important clinical parameters.Fifth, the findings of several subgroup analyses enhance the generalizability and the application of our study.We need to recognize some limitations.First, our study yielded a mean follow-up duration of 215 days (based on an on-treatment approach), which tended to be shorter than that in previous clinical trials (52-week follow-up duration based on an intention-to-treat approach) [9][10][11].However, this reflected actual use patterns of LABA/LAMA FDC in the real-world environment and was consistent with previous observational studies [7,8].Our sensitivity analysis based on an intention-to-treat approach (a mean follow-up duration of 660 days, Table 4 and Additional file 1: Table S11) and subgroup analysis with longer treatment duration (such as >365 days, Figure 2 and Additional file 1: Table S16) did not change the results materially.Second, we examine a broad spectrum of cardiovascular events to enhance clinical relevance.However, the current database did not provide information on cardiovascular death for further risk estimation.Third, while we controlled for a large number of potential confounders with several approaches, we could not fully exclude the possibility of residual or unmeasured confounding, which was an inherent limitation of most observational studies.Similarly, we attempted to integrate claims-based covariates and clinical parameters to strengthen confounding adjustment.For example, for adjusting for the influence of COPD severity, we simultaneously captured information on COPD duration, history of hospitalized COPD exacerbations, and lung function test results.However, proportions of patients with clinical parameters from COPD P4P Database tended to be low.This may be because the database was only available since 2017.The availability with corresponding information is expected to increase over time and future studies with accrual data would facilitate replication of our results.

Conclusion
In this nationwide cohort study conducted with a target trial emulation framework, there was no substantial increased risk associated with LABA/LAMA FDC compared to LABA/ICS FDC.The concerns of sympathetic activation and cardiovascular events may not preclude dual bronchodilator treatment for patients with appropriate indication.

Take home message
This population-based cohort study with a target trial emulation framework provides important reassurance about comparative cardiovascular safety of LABA/ LAMA FDC compared with LABA/ICS FDC among patients with COPD.This study was partly supported by the research grant from Yen Tjing Ling Medical Foundation (CI-110-29).The current work is based on Chun-Yu Chen's master thesis.The results were presented, in part and in an abstract, at International Society of Pharmacoepidemiology (ISPE)'s 14th Asian Conference on Pharmacoepidemiology, Tainan, Taiwan, October 21-23, 2022.However, the complete manuscript has not been published or under simultaneous consideration for publication elsewhere.Role of the Funder/Sponsor: The funding source has no role in study design, collection, analysis and interpretation of data, writing of the report, and in the decision to submit the article for publication.

Availability of data and materials
The data used in the current study were obtained from the Applied Health Research Data Integration Service from the National Health Insurance Administration, Taiwan, which are not publicly available given the data protection policy.However, the authors are willing to have further discussion if there are any questions.

Fig. 1
Fig. 1 Graphic depiction of the cohort study design.COPD chronic obstructive pulmonary disease, FDC fixed-dose combinations, ICS inhaled corticosteroids, LABA long-acting β 2 agonists, LAMA long-acting muscarinic antagonists, NHI National Health Insurance

Fig. 2
Fig. 2 Subgroup analyses for risk of composite cardiovascular events comparing LABA/LAMA FDC with LABA/ICS FDC after PS matching, by patient characteritstic and treatment duration a .CI confidence interval, COPD chronic obstructive pulmonary disease, FDC fixed-dose combinations, HR hazards ratio, ICS inhaled corticosteroids, LABA long-acting β 2 agonists, LAMA long-acting muscarinic antagonists, NA not applicable, PS propensity score.a We re-estimated PS and re-matched patients in each patient subgroup.b The HR after PS matching was stratified on the matching ratio.c COPD duration was defined as the duration from the first recorded date of COPD diagnosis (looking back until 2014/01/01) to the cohort entry date.d History of hospitalized COPD exacerbations was defined as having hospitalized COPD exacerbations within 365 days before cohort entry based on any diagnosis positions in the inpatient claims.e History of cardiovascular diseases was defined as having the following cardiovascular diseases within 365 days before cohort entry based on any diagnosis or procedure positions or health services records in the outpatient and inpatient claims, including hypertension, ischemic heart disease or angina, myocardial infarction, coronary revascularization, cardiac dysrhythmia, congestive heart failure, cerebrovascular disease, ischemic stroke, hemorrhagic strike, peripheral vascular disease, diabetes, and hyperlipidemia

Table 1
Selected patient characteristics of the eligible cohort before and after PS matching

Table 1
(continued) ACEI angiotensin converting enzyme inhibitors, ARB angiotensin receptor blockers, COPD chronic obstructive pulmonary disease, FDC fixed-dose combinations, ICS inhaled corticosteroids, LABA long-acting β 2 agonists, LAMA long-acting muscarinic antagonists, NSAID non-steroidal anti-inflammatory drugs, PS propensity score, SD standard deviation a One randomly sampled LABA/LAMA FDC initiator versus one LABA/ICS FDC initiator in each matched subset b COPD duration was defined as the duration from the first recorded date of COPD diagnosis (looking back until 2014/01/01) to the cohort entry date c Cardiovascular disease include hypertension, ischemic heart disease or angina, myocardial infarction, coronary revascularization, cardiac dysrhythmia, congestive heart failure, cerebrovascular disease, ischemic stroke, hemorrhagic stroke, transient ischemic attack, peripheral vascular disease, diabetes mellitus, and hyperlipidemia

Table 2
Number of patients and events, follow-up duration, and incidence rate of composite and individual cardiovascular events among LABA/LAMA FDC initiators and LABA/ICS FDC initiators before and after PS matching CI confidence interval, FDC fixed-dose combinations, ICS inhaled corticosteroids, LABA long-acting β 2 agonists, LAMA long-acting muscarinic antagonists, PS propensity score, SD standard deviation a The unit of incidence rate was per 1000 person-years.The incidence rate after PS matching was weighted by the inverse of the matching ratio

Table 3
Risk of composite and individual cardiovascular events comparing LABA/LAMA FDC versus LABA/ICS FDC before and after PS matching CI confidence interval, FDC fixed-dose combinations, HR hazards ratio, ICS inhaled corticosteroids, LABA long-acting β 2 agonists, LAMA long-acting muscarinic antagonists, PS propensity score a The HR after PS matching was stratified on the matching ratio

Table 4
Sensitivity analyses for risk of composite cardiovascular events comparing LABA/LAMA FDC versus LABA/ICS FDC after PS matching CI confidence interval, FDC fixed-dose combinations, hd high-dimensional, HR hazards ratio, ICS inhaled corticosteroids, LABA long-acting β 2 agonists, LAMA longacting muscarinic antagonists, PS propensity score a The HR after PS matching was stratified on the inverse of the matching ratio