Data Source
The NHIRD, provided by the National Health Insurance Administration (NHIA) of Taiwan, is a nationwide claims-based database of the National Health Insurance (NHI) program. The NHI program, launched in 1995, is a compulsory insurance program that provides reimbursement for most medical services and over 30 000 prescription drugs. The data used in the study were collected between 2000 and 2015 and were maintained by the Health and Welfare Data Science Center (HWDC), Ministry of Health and Welfare, Executive Yuan, Taiwan. The NHIRD files include inpatient, outpatient and drug prescription claims and use the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) and the Anatomical Therapeutic Chemical (ATC) system to define the patients had a specific disease diagnosis or drug prescription. To verify the accuracy of the diagnoses and the rationale for treatments, the NHIA also routinely samples and reviews a proportion of NHI claims. In addition, hospitals and clinics are penalized if they provide any unnecessary medical treatment to patients. Besides, each patient has a unique encrypted identifier which can be linked to National Death Registry under the regulation of HWDC. This study was approved by the Joint Institutional Review Board of Taipei Medical University (approval no. N201808075).
Study Cohort
The initial cohort consisted of new-onset COPD patients with diabetes between 2003 and 2014. COPD patients were defined if the patients had at least three disease diagnostic claims within a year of follow-up and at least a three-year washout period was adopted to ensure the patients were newly diagnosed with COPD. We then excluded the patients had unknown sex, were not a citizen in Taiwan, or were aged less than 40-year old (1); had no COPD prescription claim within a year after the first date of COPD diagnosis (2); had a disease history of asthma, malignant tumor, chronic kidney disease, and renal dialysis (3); had a diagnosis of type 1 diabetes or had no antidiabetic drug prescription claim or had metformin monotherapy and had insulin therapy before the first COPD diagnosis (4). The later exclusion was made to increase the homogeneity of the study cohort.
Case and Control Patient Selection
There is no general agreement on definition of AECOPD. Mostly, the definition of AECOPD was based on increasing symptoms and/or increased health care utilization. Thus, in this claim-based study, we used the following approach to identify the patients with AECOPD modified according to previous study [12, 13] if they 1) had a hospital admission or an emergency visit due to COPD and also requiring oral or injection corticosteroid (CS), or 2) added on oral or injection CS therapy in a new visit. To increase the comparability, matched controls were selected based on incidence density sampling which involved in matching each AECOPD case to a sample of those potential controls who were at risk at the time of case occurrence, resulting obtaining unbiased estimates of relative risk. Before matching, we additionally excluded the patients had monotherapy and then included subjects with double or triple combination therapy of OADs which regimen were validated by clinical trials and meta-analyses [14]. Finally, each case was matched to 4 randomly selected controls according to the propensity score estimation by sex, age, the year of COPD diagnosis, the initial year of DM status, previous and coexisting disease conditions, Charlson comorbidity index (CCI), level complexity of COPD and the COPD medication use three month prior to the date of AECOPD. The initial year of DM status was defined based on the first claim year of the patients initially received 2nd line OADs continuously for at least three month. Since AECOPD did not occur in control patients, we randomly assigned of proxy event dates to control patients which corresponded to the index date of their matched cases. When applying this method, we created a basis of comparison for the OADs exposure between case and control patients.
Exposure to oral antihyperglycemic drugs (OADs)
We examined all OAD prescription records within three month before the index date of AECOPD of cases and pseudo-AECOPD date of controls, respectively. We investigated the type of OADs, including metformin (MET), sulfonyurea (SU), a-glucosidase inhibitors (AGI), thiazolidinediones (TZD) and dipeptidyl peptidase-4 inhibitor (DPP-4i). The aim of our study is to answer what is the best drug as add-on OADs to monotherapy for progressive T2DM in the context of considering effect on COPD outcomes. Then, we further categorized T2DM-COPD patients by using a double or triple combination of OADs.
Potential Confounding Variables
Previous or coexisting medical conditions were recorded if patients were diagnosed with chronic artery disease (CAD), hypertension (HTN), congestive heart failure (CHF), pneumonia, chronic liver disease (CLD), dementia/Parkinson and osteoporosis. Additionally, CCI which is represented the severity of comorbid conditions of patients was also considered a major risk and the CCI in this current study has been modified since all patients had diagnosed with both diabetes and COPD but did not have a history of malignant neoplasm. To adjust the severity of COPD itself, we categorized patients into low, moderate and high complexity according to the previous study [15] and further grouped the patients into low and moderate/high complexity group due to a small sample size of high complexity. Besides, we also considered the history of COPD medication use of AECOPD cases and non-AECOPD controls, respectively [15], including short acting beta agonists (SABA), short-acting muscarinic antagonists (SAMA), long-acting beta agonists (LABAs) , long-acting muscarinic antagonists (LAMAs) and inhaled corticosteroids (ICS).
Statistical Analysis
The baseline differences between case and control patients were measured by standardized mean difference (SMD). Conditional logistic regression was used to estimate the odd ratios (OR), adjusted odds ratios (aOR), and 95% confidence intervals (CI) for the association of AECOPD risk and OADs use. The statistical analyses were performed using SAS/STAT, Version 9.4, (SAS Institute, Cary, NC, USA) and STATA 13 (Stata Corp, College Station, TX, USA). A P value < 0.05 and SMD >0.1 was set as the level of statistical significance.