Drug Safety of Benzodiazepines in Asian Patients With Chronic Obstructive Pulmonary Disease

Purpose: Many comorbidities, including depression, anxiety, and insomnia, occur in patients with chronic obstructive pulmonary disease (COPD). These patients may be prescribed benzodiazepines (BZDs). However, there are some concerns that benzodiazepines increase the risk of drug overdose, hypercapnic respiratory failure, acute exacerbation and increased mortality. The aim of our study was to evaluate the drug safety of BZDs in patients with COPD. Methods: We used the National Health Insurance Research database in Taiwan from 2002 to 2016 to perform a retrospective cohort study. We enrolled patients who were exposed to the first prescription of BZDs, non-BZDs or a combination (mix user) after COPD diagnosis. We performed 1:1:1 propensity score matching in three groups. The outcomes were COPD with acute exacerbation and all-cause mortality. Poisson regression analysis was performed to evaluate the incidence rate ratios for the outcomes in the groups. Results: After propensity score matching, there were 2,856 patients in each group. After adjusting for confounding factors, we found that compared to BZD users, non-BZD and mix users had nonsignificant differences in outpatient management of acute exacerbations, hospitalization management of acute exacerbations, emergency department management of acute exacerbations and all-cause mortality. BZD and mix groups showed significantly increased admission for acute exacerbation of COPD compared with that of the nonuser group, with IRRs of 2.52 (95% CI, 1.52–4.18; p = 0.0004) and 2.63 (95% CI, 1.57–4.40; p = 0.0002), respectively. Conclusion: BZD, non-BZD, and mix users showed increased COPD-related respiratory events compared to nonusers in Asian subjects.


INTRODUCTION
Chronic obstructive pulmonary disease (COPD) is a chronic lung disease that is not fully reversible and is characterized by airflow obstruction. COPD is associated with a substantial health burden worldwide (Soriano and Rodriguez-Roisin, 2011;World Health Organization 2017;Terzano et al., 2010;Decramer et al., 2012). Comorbidities such as insomnia, anxiety, sleep disorder, depression, and psychological disorders are commonly reported in patients with COPD (Light et al., 1985;Kunik et al., 2005;Solano et al., 2006;Cheng et al., 2008;Maurer et al., 2008;Biswas et al., 2017). Benzodiazepines (BZDs) are the most widely prescribed class of sedative-anxiolytic drugs; BZDs cause sedation and muscle relaxation, can lower anxiety levels and are commonly used for the treatment of insomnia, anxiety and sleep disorder. However, previous studies reported that BZDs may have considerable risks or adverse effects and can be fatal in patients with COPD; these adverse effects include hypoxemia, hypercapnia (Block et al., 1984;Beaupre et al., 1988), decreased respiratory muscle strength (Jolly et al., 1996) and respiratory failure (Chen et al., 2015). The risk may increase as the plasma half-life of BZDs increases (Greenblatt et al., 1991;Fisher, 1999).
Benzodiazepine receptor agonists (BZRAs) include traditional BZDs and newer generation non-benzodiazepine receptor drugs that preferentially bind to the ω1-benzodiazepine receptor of the GABAA receptor complex. BZDs nonselectively bind to the receptor (Ebert et al., 2006), whereas non-BZD drugs selectively bind to the benzodiazepine receptor and have a lower affinity for the GABAA receptor; therefore, they lack significant muscle relaxant, anxiolytic, and anticonvulsant activities of traditional BZDs, resulting in fewer pulmonary adverse effects and fewer adverse effects on other systems (Ranlov and Nielsen, 1987;Cohn, 1993;Dämgen and Lüddens, 1999;George, 2000;Roth, 2009).
The Joint American Thoracic Society/European Respiratory Society guidelines do not recommend the use of BZDs for COPD patients, especially when COPD is severe (Celli and MacNee, 2004). Nevertheless, BZRAs are widely used for patients with COPD, with estimates that 44.7%-69% of patients receive these drugs (Halvorsen and Martinussen, 2015;Park et al., 2017). Elderly COPD patients increase the prevalence, exposure duration and dose of BZRAs because of the patient's illness, psychological dependence or other factors (Steinman et al., 2017).
A concern in clinical practice is whether COPD patients using BZDs or non-BZD drugs, or a combination of the two, are potentially at increased risk of respiratory problems. Some studies are available for BZRAs used in patients with COPD, and these studies were limited by a small number of subjects and short-term effects (lasting a few hours) (Murciano et al., 1990;Berry et al., 1992;Murciano et al., 1993).
Therefore, the aim of the study was to evaluate the risk of acute exacerbation in patients with COPD after receiving BZRAs or non-BZRAs and in nonusers by performing a large, populationbased cohort study in an Asian population.

Identification of the COPD Cohort
Using the database, we identified adults aged between 40 and 90 years with newly diagnosed COPD (ICD-9-CM codes 490, 491, 492, and 496) who had received more than one inpatient diagnosis or two or more outpatient diagnoses between January 1, 2002, and December 31, 2015. Patients with COPD should receive one of the following medications for more than one month: longacting β adrenergic agonists (LABAs); long-acting antimuscarinic agonists (LAMAs); or oral methylxanthines combined with short-acting β2-agonists (SABAs) or shortacting muscarinic antagonists (SAMAs) (McKay et al., 1993;Ram et al., 2002;Wang et al., 2018). Patients were excluded if they met any of the following conditions: asthma, lung transplantation, lung cancer or death within one month after diagnosis of COPD (Ekstrom et al., 2013). Moderate COPD exacerbation was defined as a requirement for treatment with antibiotics or systemic corticosteroids or both. Severe exacerbations of COPD was defined as those leading to hospitalization or emergency department visit.

Identification of the Benzodiazepine Receptor Agonists Cohort
The BZRA cohort was identified from COPD patients newly exposed to BZRAs and divided into an oral BZD group and an oral non-BZD BZRA group. Patients were required to continue using the BZRA without change for 30 days. Patients were excluded if they received BZRA injections during the study period or died before the index date. Each patient was followed until an outcome occurred: mortality or the end of the study (December 31, 2016). The BZRA cohort was divided into three groups based on the initial prescription: those taking BZDs (the BZD cohort), those taking non-BZD BZRAs (the non-BZD cohort), and those prescribed both BZDs and non-BZD drugs simultaneously in the initial prescription (the mix cohort).
Frontiers in Pharmacology | www.frontiersin.org December 2020 | Volume 11 | Article 592910 A cohort of patients with COPD who were not taking BZRAs (the nonuser cohort) was also identified.

Exposure and Index Date Assessment
A new prescription for BZRAs was defined as the patient receiving his or her first prescription for BZRAs during the study period following entry into the cohort. We excluded patients who received their first BZRA prescription prior to the diagnosis of COPD. When patients received more than prescription during the study period, only the first was included in the analysis. The index date was the date of the first BZD, non-BZD or Mixed-use prescription after COPD diagnosis. For each patient, we collected any of the following outcomes that occurred during the 30 days following the index date (The index date was the date of the first BZD, non-BZD or Mixed-use prescription after COPD diagnosis and patients were required to continue using the medication without change for 30 days): outpatient visits for respiratory exacerbations, hospital admission for COPD acute exacerbation or for respiratory exacerbations, emergency department attendance for COPD or pneumonia and all-cause mortality. The observation period was short (only 30 days) and only the first prescription was included in the analysis.

Covariates
The baseline characteristics and comorbidities of each patient were collected within one year prior to the index date. These included COPD severity, hypertension, diabetes mellitus, ischemic heart disease, heart failure, atrial fibrillation, arrhythmia, chronic kidney disease, malignancy, depression and anxiety. In addition, the frequency of COPD acute exacerbations was recorded. Severe acute exacerbation was defined as patients requiring admission or visits to the emergency department. Moderate acute exacerbation was defined as a patent prescription for either antibiotics or oral corticosteroids from an outpatient department. Other COPD treatment medications, including LABAs, LAMAs, inhaled corticosteroids, SABAs, SAMAs, systematic beta-2 agonists and methylxanthines, were recorded for analysis.

Propensity Score Matching to Create the Final Groups
Propensity score matching was applied on a 1:1:1 basis to the BZD, non-BZD and mix cohorts to create three final groups (the BZD, non-BZD and mix groups) with balanced baseline characteristics. The matching factors included age group, sex, comorbidities and COPD medication (Lori, 2004). Propensity score matching was also used to match the patients in the three BZRA groups and the patients in the nonuser cohort 1:1 to perform sensitivity analyses.

Outcome Assessment
For each patient, we collected any of the following outcomes that occurred during the 30 days following the index date: outpatient visits for respiratory exacerbations, hospital admission for acute COPD exacerbation or for respiratory exacerbations, emergency department attendance for COPD or pneumonia and all-cause mortality.

Statistical Analysis
Differences in baseline characteristics were evaluated by the chisquare test (for categorical variables). Univariate and multivariate Poisson regression with robust error variance analysis was used to generate crude and adjusted incidence rate ratios (IRRs) with 95% CIs for each type of outcome. Consequently, the IRRs for all the study outcomes and further sensitivity analyses were evaluated by the Poisson regression model with a robust error variance (Zou, 2004) between the BZD, non-BZD and mix groups. Time-toevent analyses were performed using Kaplan-Meier plots and the log-rank test. Various sensitivity analyses were performed. First, we used several different follow-up periods (60, 90, 180, and 365 days) to ensure that the results were consistent with those obtained for the 30-days follow-up. Second, we compared risks between the three BZRA groups and the nonuser group. Third, because BZDs with a longer half-life would be expected to cause more serious side effects than those with shorter half-lives (Greenblatt et al., 1989;Greenblatt et al., 1991), we stratified the BZD group according to half-life to reduce any pharmacokinetic effects. SAS statistical software (Version 9.4; SAS Institute, Inc., Cary, NC) was used to perform all the analyses. A two-sided p value < 0.05 was considered statistically significant.

RESULTS
There are 267,741 patients in COPD cohort. There are 201,624 patients used BZRA and 66,117 (24.69%) patients did not have BZRA. In total, 267,741 patients were included in the COPD cohort (Figure 1), of whom 82,675 were first prescribed BZRA drugs after their COPD diagnosis. After propensity score matching, there were 2,856 patients in each of the three groups. All of the baseline characteristics were balanced after propensity score matching ( Table 1). According to our insurance regulation, benzodiazepines are used for a range of some health issues, including anxiety, sleep disorders, and epilepsy. Most of our study cohort (70%) were used BZRA for sleep disorders. Table 2 summarizes the IRRs. The non-BZD groups experienced significantly fewer outpatient visits for acute exacerbation than the BZD group (IRR 0.84, 95% CI 0.72-0.97, p 0.0207). However, compared with the BZD group, the non-BZD group was not associated with significant decreases in admission for COPD acute exacerbation (p 0.4064), emergency department attendance for COPD acute exacerbation (p 0.5577), or all-cause mortality (p 0.9126). Similarly, compared with the BZD group, the mix group was not associated with significant IRRs for outpatient visits for acute exacerbation (p 0.7957), admission for COPD acute exacerbation (p 0.9660), emergency department visit for COPD acute exacerbation (p 0.9736), or all-cause mortality (p 0.3922).
The clinical outcomes were evaluated further by three sensitivity analyses ( Table 3). The analysis using various follow-up periods showed two inconsistencies with the main Frontiers in Pharmacology | www.frontiersin.org December 2020 | Volume 11 | Article 592910 findings for 30 days. At 90 days, there was no longer a decrease in outpatient visits for respiratory exacerbation in the non-BZD group compared to that in the BZD group, and at 180 days, admission for the acute exacerbation of COPD was inconsistent with the 30-days finding. All the other results in the follow-up sensitivity analysis were consistent with the main findings.

DISCUSSION
This retrospective observational cohort study examined data for 267,741 patients with COPD, of whom 82,675 used BZRAs. The comparison of the outcomes among the three types of BZRAs found that the patients administered the non-BZD drugs Frontiers in Pharmacology | www.frontiersin.org December 2020 | Volume 11 | Article 592910 experienced 0.84-fold fewer outpatient visits for respiratory exacerbation in the first 30 days than did those who received BZDs. This was not influenced by baseline comorbidity or the severity of COPD. The other results showed no statistically significant differences among the use of BZDs, non-BZD BZRA drugs or a combination of these drugs. Previous studies of patients with COPD have shown inconsistent results, with some finding that BZDs did not influence respiratory effects (Beaupre et al., 1988;Murciano et al., 1990), whereas others found that it did (Murciano et al., 1993;Chen et al., 2015). However, none of these studies was designed as a cohort study with a 30days outcome period, and all of the studies focused on differences in lung parameters rather than on COPD-related outcomes.
Therefore, the present study provides new evidence that patients using BZD and non-BZD drugs or a combination were at equal risk of COPD-related exacerbation. This finding allows for flexibility in selecting a BZRA to treat COPD. However, although there were no statistically significant results for most respiratory-related outcomes, the non-BZD group experienced fewer outpatient visits for acute exacerbation. A sensitivity analysis checked whether a longer follow-up period affected the clinical outcome. Although two of the main outcomes changed at particular follow-up periods, most of the outcomes remained consistent with the main outcomes at 30 days. Any side effects of BZRAs would be expected to be observed soon after the drugs were administered, so the 30-day period of   follow-up in this study is acceptable. The results of our second sensitivity analysis that compared the outcomes with those of COPD patients who did not receive BZRA drugs revealed an association with increased respiratory exacerbation, with the BZD and mix groups experiencing greater risk of the occurrence of a clinical outcome than the nonuser group. This was consistent with the results of a cohort study by Vohoris et al., who found that COPD patients using BZDs had a higher association with emergency and outpatient visits for exacerbation than the patients who did not take BZRAs (Vozoris et al., 2014). Other previous studies have also  observed an increased risk of respiratory exacerbation in COPD patients administered BZRA compared to those who were not (Beaupre et al., 1988;Berry et al., 1992;Cohn et al., 1992;Berry et al., 1995;Jolly et al., 1996).
Our study examined population-based, real-world data, adjusted according to COPD severity, and provided evidence related to various follow-up periods and for BZRAs with different half-lives. The non-BZD group showed a risk equal to that of the other groups regarding admission for acute exacerbation or emergency department for respiratory exacerbation but a lower risk of outpatient visits for exacerbation. No previous study has considered the risk of using BZRAs for patients with COPD. Some experimental studies have shown that non-BZD drugs do not affect the lung function parameters of COPD patients in the short term (Murciano et al., 1990;Murciano et al., 1993;Girault et al., 1996). However, one study showed that the non-BZD group showed fewer acute exacerbation events (Chen et al., 2015).
Our study compared the non-BZD and BZD groups based on propensity score matching for age group, sex, COPD severity, comorbidities and COPD medication.
Compared with the BZD group, the non-BZD group underwent fewer outpatient visits for acute exacerbation; however, other clinical COPD-related exacerbation outcomes were similar for the non-BZD, BZD and mix groups, including admission for acute COPD exacerbation, emergency department attendance for COPD exacerbation, and all-cause mortality. In a sensitivity analysis, we examined various follow-up times to check our outcomes and reduce selection bias. Compared to the nonuser group, all three BZRA groups were associated with an increased risk of COPD-related exacerbation. Our results provide treatment information in clinical practice and provide potentially useful data about BZRA risk.

Strengths and Limitations
This study and its findings had some limitations. First, we were unable to obtain direct information about COPD severity from the database. Instead, two surrogate values were used to estimate the severity. There was no information regarding the severity of COPD. We need to evaluate how this correlates with the increased risks. Second, we were unable to establish an indication for using BZRA, which requires a prospective study design. The individual indication of BZRA might have provided information related to underlying diseases of the COPD patients and about particular clinical situations or procedures. Consequently, we performed propensity score matching to balance the severity of COPD between the groups and reduce the bias. Third, the study was based on dispensed prescriptions; therefore, we could not establish the patients' drug adherence. Despite the limitations, our study had some important strengths. The findings were derived from a population-based database. This study evaluated COPD patients divided into BZD, non-BZD and mix cohorts. The observation period is short, therefore the study provides information on short-term effects only, so further studies are needed. Although we used a short period to compare the risk between the groups, we found the same trends in the longer-term sensitivity analysis. We found no differences among the three BZRA groups; however, compared to the nonuser group, the patients taking BZRAs were at increased risk of respiratory events.

CONCLUSIONS
This study has provided clues that the risk of short-term COPD acute exacerbation is similar when using BZDs, non- BZD drugs or a combination of these drugs. In addition, the use of BZDs and non-BZD BZRAs by patients with COPD increased their risk of acute exacerbation compared with that in COPD patients who did not take BZRAs in Asian subjects. Therefore, COPD patients prescribed BZRAs should be observed for potential risk of acute exacerbation in Asian subjects.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

ETHICS STATEMENT
The study was approval by KMUHIRB-EXEMPT-20180043.

AUTHOR CONTRIBUTIONS
Conceived and designed the experiments: C-YC. Performed the analysis: L-YC. Analyzed the data: L-YC and C-YC. Wrote the paper: Y-HL, L-YC, and K-ML. Provided constructive opinions and suggestions and study supervision: C-YC. All authors read and approved the final manuscript.