Buprenorphine is an effective treatment for opioid use disorder, contributing to significant reductions in all-cause and opioid-associated mortality amid the current U.S. opioid epidemic (
1). Among patients with opioid use disorder who are taking buprenorphine, benzodiazepine use is highly prevalent. Some studies have estimated that as many as 30% of patients with opioid use disorder who receive opioid maintenance treatment also receive benzodiazepine prescriptions (
2), with more than one-third of these patients endorsing past-month problematic use of benzodiazepines (
3). Although benzodiazepines are frequently prescribed for treatment of comorbid mood and anxiety disorders, which are common in this patient population, recent research has shed light on respiratory depression, overdose risk, and addictive potential associated with benzodiazepine use in patients taking opioids chronically, such as buprenorphine (
4 –
7). It is unclear whether the risks associated with benzodiazepine use outweigh the treatment benefits of buprenorphine.
The relationship between benzodiazepine use and buprenorphine treatment outcomes is poorly characterized. Existing observational studies that examined mortality risk associated with benzodiazepine prescriptions among patients with opioid use disorder have produced contradictory results (
8 –
10). While some findings suggest that benzodiazepines may enhance retention in buprenorphine maintenance treatment (
4,
11), benzodiazepines have also been associated with increases in drug-related poisonings (
11 –
14), all-cause mortality (
11), nonoverdose deaths (
8), decreased retention in treatment (
12,
15 –
17), and accidental injury-related emergency department visits (
5). To our knowledge, no studies have specifically examined potential additive or interactive effects between buprenorphine and benzodiazepines.
Another limitation of the current research base for benzodiazepine-related morbidity and mortality is that the different pharmacologic properties of sedative/hypnotics, such as half-life and potency, have seldom been explored (
18). This is an important clinical and scientific gap, given that shorter-acting benzodiazepines, such as alprazolam, are thought to have higher addictive potential than longer-acting medications, such as clonazepam (
6,
19). In fact, previous studies of populations without opioid use disorder have shown significant intra-benzodiazepine variation in health outcomes; for example, certain benzodiazepines appear to be associated with greater injury risk in the elderly than others (
18). Additionally, little research has investigated drug-related poisoning associated with selective benzodiazepine receptor modulators (Z-drugs: zolpidem, zaleplon, and eszopiclone). These medications have increasingly been found to have a spectrum of adverse effects similar to that for benzodiazepines, with similar dose-response effects on all-cause mortality in the general population (
20,
21). Because studies to date have relied on relatively small numbers of benzodiazepine users, they have precluded examination of how drug subclasses (e.g., short-acting compared with long-acting benzodiazepines, Z-drugs) or dosing regimens affect opioid use disorder outcomes.
In light of the research gap on opioid use disorder outcomes among patients taking both buprenorphine and benzodiazepines or Z-drugs, we used the IBM MarketScan Research Databases (2006–2016; IBM, Armonk, N.Y.) in the present study to quantify odds of nonfatal drug-related poisoning (including overdoses) associated with benzodiazepine or Z-drug use in patients with opioid use disorder; to determine whether benzodiazepine or Z-drug use improves, nullifies, or reverses the protective effect of buprenorphine in patients with opioid use disorder; and to evaluate whether different sedative/hypnotic subtypes (e.g., long-acting compared with short-acting agents and Z-drugs) correspond to different poisoning risks.
Methods
Data Set and Study Subjects
We used pharmaceutical claims data on 304,676 individuals, ages 12–64, in the IBM MarketScan Commercial and Multi-State Medicaid Databases who received buprenorphine treatment for opioid use disorder (
Figure 1). The data, spanning from January 1, 2006, to December 31, 2016, represent claims for insured active employees, early (non-Medicare) retirees, and dependents insured by employer-sponsored plans, as well as Medicaid insurance holders, spanning all 50 states (
22). Details on data extraction are summarized in the Methods section in the
online supplement. Because all data were de-identified, our analysis was exempt from human subjects review by the institutional review board at Washington University, St. Louis.
Individuals ages 12 to 64 with insurance claims indicating an opioid use disorder diagnosis, at least one buprenorphine prescription, and at least one nonfatal drug-related poisoning were included for analysis. Buprenorphine prescriptions were defined by associated diagnostic codes selected using established methods (
22) (for the list of codes, see Table S1 in the
online supplement). We included persons who had periods of buprenorphine-free treatment as long as they received the medication at another time point. Buprenorphine use was characterized in terms of strength, quantity, and days of supply in order to calculate a daily milligram dose. This was further stratified into daily buprenorphine doses >12 mg and ≤12 mg, given previous analyses suggesting differences in treatment retention associated with this dose (
23). Because each individual serves as his or her own control in case-crossover designs, all individuals in the sample must experience the outcome at least once. Therefore, we excluded all individuals who did not experience a drug-related poisoning.
Study Design
We used a case-crossover study design that exploited within-subject variation in exposure and outcome, which addresses unobserved time-invariant confounding by allowing each individual to serve as his or her own control. Units of observation were person-days, denoting days during which patients were enrolled in a health insurance plan. Case periods were days when a patient experienced nonfatal drug-related poisonings. Control periods were nearby days without poisoning events (
24,
25). We characterized each person-day of observation by the presence or absence of benzodiazepine or Z-drug treatment and the presence or absence of buprenorphine treatment. Each study subject could have multiple poisoning events as long as these events fell within the interval of ≤365 days before and after the index poisoning (
Figure 2). Individuals with fewer observation days on either side of the index event were included with missing days treated as censored. We evaluated study subjects during a time period up to 1 year before and the year after their first drug-related poisoning (index date), thus limiting subjects to a maximum of 2-year periods of observation to reduce heterogeneity in observation time per person. We included periods both before and after the index poisoning in order to define the poisoning as a repeatable event, allowing for the inclusion of time as a covariate (
25,
26). Thus, the unit of analysis was days, stratified within persons, and the comparison of interest was drug exposure status concurrent with drug-related poisoning (case period) compared with exposure status during referent periods when poisoning events did not occur (control period).
Given their high prevalence in clinical practice and overlapping indications with benzodiazepines, selective serotonin reuptake inhibitors (SSRIs) (sertraline, fluoxetine, escitalopram, and citalopram) were included in our conditional logistic models as an active comparator analysis. This is a validated approach employed by previous investigators (
27) to control for unobserved confounding. In other words, the SSRI term assesses whether underlying conditions for which patients were receiving benzodiazepines or Z-drugs contributed to drug-related poisonings.
Ascertainment of Outcomes and Exposures
The primary outcome was nonfatal drug-related poisonings (including overdoses) (
Figure 2). Controls were adjacent person-days when a poisoning event did not occur. Poisonings were defined using ICD codes and Healthcare Common Procedure Coding System codes for naloxone reversal. In accordance with guidelines compiled by the Centers for Disease Control and Prevention consensus recommendations for poisoning surveillance (
28), the codes we employed encompassed not only opioids but also alcohol, benzodiazepines, psychotropic medication, and overdose with other substances (see Table S1 in the
online supplement). This was intended to capture conditions commonly associated with chronic drug misuse, as well as to avoid misclassification bias and underestimation of overdose risks. We searched all emergency department visits, outpatient ambulatory visits, and hospitalizations for relevant diagnostic codes (
28).
The primary exposure was benzodiazepine or Z-drug prescriptions, ascertained through pharmacy files. A person was assumed to be exposed if he or she was covered by a prescription for a drug on a given day. We characterized initial benzodiazepine, Z-drug, and buprenorphine prescriptions in terms of strength, quantity dispensed, and days of supply. To standardize the daily dose, we calculated the strength of each benzodiazepine or Z-drug in terms of total diazepam-equivalent milligrams (
29) using known pharmacologic conversion factors (see Table S2 in the
online supplement) (
30 –
32). We calculated a daily diazepam-equivalent dose by multiplying the number supplied by strength (in diazepam-equivalent milligrams) and dividing by days of supply. Benzodiazepine and Z-drug dosage was stratified into high-dose (diazepam-equivalent-mg dose >30 mg) and low-dose (≤30 mg) based on established thresholds (
33). Benzodiazepine exposure was categorized by the duration of action, namely, short-acting (half-life ≤24 hours) or long-acting (half-life >24 hours), guided by established classifications (
6). Specific classifications of benzodiazepines and Z-drugs, based on previously published definitions, are presented in Table S2 in the
online supplement (
34,
35).
We also collected data on age, sex, relationship of patient to the primary beneficiary, and insurance status. Because cases were self-matched in a case-crossover design, our models did not warrant adjustment for time-invariant measured confounders (e.g., sex and race).
Statistical Analysis
Analyses were conducted on a day-level data set in which each medication coverage date was identified. All analyses were conducted using SAS, version 9.4 (SAS Institute, Cary N.C.). Before evaluating variation in medication exposure between case and control days, we calculated descriptive statistics for the primary analytical sample (i.e., patients who had nonfatal drug-related poisonings). To compare the characteristics of opioid use disorder treatment between patients with and without poisoning events, we obtained a random sample of persons who never experienced poisonings and compared descriptive statistics between them and those who had poisoning events (see the Methods section and Tables S3 and S4 in the online supplement).
To evaluate the effect of medication exposure, we estimated conditional logistic regression models stratified by subject and modeled the risk of poisoning as a function of drug exposure by days with or without treatment. We first built a crude model containing benzodiazepine or Z-drug prescription status as the only predictor variable. We then built additional models that included predictor variables for both benzodiazepine or Z-drug and buprenorphine prescription status. Simultaneous inclusion of benzodiazepines or Z-drugs and buprenorphine permitted us to model additive or interactive effects of benzodiazepines or Z-drugs and buprenorphine in association with drug-related poisonings. SSRIs were included in our models as an active comparator analysis. Subgroup analyses were conducted to assess the effect of buprenorphine treatment days, compared with days without treatment, on drug-related poisoning among patients who received benzodiazepine or Z-drug prescriptions and those who did not. Controls for both calendar time and time from index poisoning were included using cubic spline methods, described in detail in the Methods section in the online supplement. Several sensitivity analyses were conducted to evaluate the possibility of persistent user bias and to assess the robustness of our findings to alternative time samplings (see Figure S1, the Methods section, and Table S5 in the online supplement).
Results
Study Sample and Treatment Characteristics
Of 304,676 individuals ages 12–64 with an opioid use disorder diagnosis, prescription drug coverage, and at least one opioid use disorder claim during the period of enrollment, we excluded individuals without drug-related poisonings, individuals who never received medication for opioid use disorder, and individuals without days of naltrexone and methadone treatment and days of observation outside a maximum of a 1-year period before and after the index poisoning (
Figure 1). Our final analytic sample comprised 23,036 patients (mean age, 30 years [SD=12.15]; 51% male; mean observation time, 299 days [SD=107.88]), spanning 14,213,075 person-days of insurance coverage (
Table 1). A total of 2,210,927 person-days (15.6%) entailed claims for buprenorphine (mean daily dose, 15.4 mg [SD=7.31]). A total of 1,968,944 person-days (13.9%) entailed claims for benzodiazepines or Z-drugs, of which 474,181 person-days entailed concurrent buprenorphine treatment. We calculated the mean daily dose of any benzodiazepine or Z-drug to be 23.4 diazepam-milligram equivalents and the mean daily dose for short-acting benzodiazepines, long-acting benzodiazepines, and Z-drugs to be 25.3, 31.3, and 4.9 diazepam-milligram equivalents, respectively (
Table 2).
Benzodiazepines or Z-Drugs and Drug-Related Poisonings
Buprenorphine treatment days were associated with 37% lower odds of drug-related poisoning (95% CI=0.60, 0.66) compared with nontreatment days, whereas odds of poisoning increased 81% on days when patients were treated with benzodiazepines or Z-drugs (95% CI=1.73, 1.91; model 1) (
Table 3). When we evaluated benzodiazepines and Z-drugs separately, Z-drug treatment days were associated with increased odds of poisoning events (odds ratio=1.29, 95% CI=1.19, 1.39), but this was notably lower than the odds associated with benzodiazepine treatment days (odds ratio=1.88, 95% CI=1.78, 1.98; model 2). We found no association between SSRI treatment days and drug-related poisonings (odds ratio=0.95, 95% CI=0.90, 1.00; model 3). We observed no difference in the magnitude of protective effect against poisoning conferred by buprenorphine treatment days when conducting stratified analyses of patients who used benzodiazepines or Z-drugs (odds ratio=0.64, 95% CI=0.60, 0.69, model 4) and those who never used benzodiazepines or Z-drugs during the study’s observation period (odds ratio=0.64, 95% CI=0.59, 0.69; model 5).
Subgroup and Sensitivity Analyses
We stratified benzodiazepines by mechanism and observed similarly elevated odds of drug-related poisoning for short-acting benzodiazepine treatment days (odds ratio=1.86, 95% CI=1.75, 1.97; model 6) and long-acting benzodiazepine treatment days (odds ratio=1.68, 95% CI=1.54, 1.83; model 6). When we stratified benzodiazepines by low and high doses, we found the increases in the odds of a poisoning event to be 78% (95% CI=1.67, 1.88; model 6) and 122% (95% CI=2.03, 2.43; model 7), respectively. To assess whether this dose-response relationship remained significant across all sedative/hypnotic drugs, we grouped benzodiazepines and Z-drugs together and stratified them into low-dose and high-dose strata, observing a similar pattern of increased poisoning odds associated with dose (low-dose: odds ratio=1.86, 95% CI=1.77, 1.79; high-dose: odds ratio=2.53, 95% CI=2.35, 2.73; model 8). In contrast to benzodiazepines, we found similar protective effects of buprenorphine irrespective of the daily dose threshold (low-dose: odds ratio=0.62, 95% CI=0.57, 0.67; high-dose: odds ratio=0.63, 95% CI=0.59, 0.67; model 9).
Interaction Between Benzodiazepines or Z-Drugs and Buprenorphine
The odds of drug-related poisoning associated with benzodiazepine or Z-drug use, stratified by benzodiazepine or Z-drug dose and buprenorphine course, are presented in
Table 3 (model 10). Stratification on these variables allowed for observation of possible interaction effects between benzodiazepine or Z-drug dose and buprenorphine use. For high-dose benzodiazepines or Z-drugs, poisoning odds associated with buprenorphine course were higher than would have been expected based on the assumption of additive risks (Wald χ
2=4.02, p=0.045). Odds ratios associated with low-dose and high-dose benzodiazepine or Z-drug course with buprenorphine use were 1.11 (95% CI=1.00, 1.23) and 1.64 (95% CI=1.39, 1.93), respectively. Notably, these effects were lower than those for low-dose and high-dose benzodiazepines or Z-drugs in the absence of buprenorphine (low-dose: odds ratio=1.69, 95% CI=1.60, 1.79; high-dose: odds ratio=2.23, 95% CI=2.04, 2.45, respectively). In other words, individuals with opioid use disorder taking both buprenorphine and benzodiazepines or Z-drugs had a lower net risk of poisoning than those taking benzodiazepines or Z-drugs without buprenorphine.
Discussion
In this study, we used the IBM MarketScan databases to analyze nonfatal drug-related poisoning (including overdose events) associated with specific benzodiazepines or Z-drugs and dosing regimens in patients with opioid use disorder. Even though buprenorphine treatment days conferred a nearly 40% reduction in poisonings, benzodiazepine or Z-drug treatment days corresponded to a near-doubling in poisoning risk. While individuals taking both buprenorphine and benzodiazepines or Z-drugs were at elevated risk of poisoning, they still had a lower net risk than those taking benzodiazepines or Z-drugs without buprenorphine. We found no association between SSRIs and drug-related poisoning, making it less likely that underlying conditions for which patients were receiving benzodiazepines or Z-drugs contributed to overdoses. Additionally, we found no dose dependence for the association between buprenorphine and poisoning risk.
Our findings may have significant implications for clinical practice. First, our results show a clear dose-dependent pattern of worsened overdose-related outcomes associated with benzodiazepine or Z-drug use, indicating that among patients with opioid use disorder for whom benzodiazepine or Z-drug cessation is risky, lower doses and shorter treatment duration of sedative/hypnotics may reduce risk. Second, we found slightly lower risk with long-acting benzodiazepines compared with short-acting benzodiazepines and substantially lower risk associated with Z-drugs compared with either long- or short-acting benzodiazepines, which may be related to lower mean standardized dosages observed for Z-drugs. Overall, these results suggest that switching benzodiazepine users from short-acting to long-acting agents or to Z-drugs may hold promise in lowering overdose risk.
Importantly, buprenorphine’s beneficial effect was observed among both benzodiazepine or Z-drug users and nonusers. This is an important finding because the safety of benzodiazepine or Z-drug use among buprenorphine users has been a topic of contentious debate in the past, leading to restrictions on buprenorphine treatment among benzodiazepine users that were later reversed by the U.S. Food and Drug Administration (
4,
36). Our results demonstrate that even though benzodiazepines and Z-drugs may increase drug-related poisonings, buprenorphine’s protective effect is not eliminated by benzodiazepine or Z-drug treatment. Overall, our findings suggest that for patients taking both benzodiazepines or Z-drugs and buprenorphine, dose reduction of benzodiazepines or Z-drugs—while maintaining buprenorphine treatment—may have the advantage of decreasing overdoses associated with both benzodiazepines and the benzodiazepine-buprenorphine interaction. Because previous studies have found benzodiazepines to be associated with improved treatment retention in patients receiving buprenorphine (
4,
11), dose reduction of benzodiazepines or Z-drugs may be preferable to abrupt cessation.
There are several limitations of this study. First, despite use of an active comparator and case-crossover design, we cannot exclude the possibility of residual confounding by indication. For example, benzodiazepines may constitute a proxy for underlying anxiety symptoms that lead to sedative/hypnotic prescriptions, substance use, and overdose, as opposed to directly contributing to poisoning. Unmeasured exposures, such as illicit substances and nonprescribed benzodiazepines, have commonly been noted in the opioid user population (
37) and warrant further investigation. Second, secular time trends in exposure and outcome may introduce confounding into case-crossover designs (
24,
25). We cannot rule out the possibility of unmeasured time-varying factors associated with drug-related poisoning and benzodiazepine or Z-drug exposure. Mitigating this, we made efforts to control for temporal variation and reduce heterogeneity in observation time per person using calendar time and time from event as a covariate and restricting study subjects to a maximum of 2-year periods of observation bracketing the index event. In addition to an active comparator analysis, we used bidirectional sampling (before and after index poisoning), which has been found to reduce overlap bias resulting from control period selection as a function of event times (
25).
In addition, our study is limited by its focus on nonfatal drug-related poisonings as opposed to poisoning deaths. While our reliance on nonfatal poisoning allows for the creation of a more robust clinical risk profile that can be used to prevent fatal cases, further research evaluating the effect of benzodiazepine or Z-drug use on overdose deaths is warranted in light of our findings. Additionally, although the MarketScan data comprise a large, nationally representative sample with strong longitudinal follow-up at the patient level, the database is limited by measurement error, insofar as medication coverage does not always reflect the actual dose consumed. The generalizability of the MarketScan data is also limited to insured patients with observed drug-related poisonings. Our findings may not be generalizable to lower-risk patients who did not experience poisoning events within the study observation period.
A key strength of our study was its use of a case-crossover approach to harness within-person variation (i.e., individuals serving as their own controls) and to estimate the degree to which poisonings were reduced on days when patients were taking benzodiazepines or Z-drugs, as opposed to nontreatment days. Because each patient acted as his or her own control, we reduced the selection and sampling bias that would result from recruitment of different case and control subjects in conventional observational study designs. This also allowed us to examine the relationship between common exposures and less common acute outcomes, such as overdoses (
18,
38,
39). Furthermore, this strategy allowed us to investigate subtle changes in exposure (i.e., variation by benzodiazepine or Z-drug dose and subtype) that may result in poisoning (
18). We also used SSRIs as an active comparator to control for unobserved confounders and to serve as a proxy for underlying conditions leading to benzodiazepine prescriptions, given their overlapping indications with benzodiazepines and no known association with drug-related poisoning or mortality in adults. The lack of association between SSRI prescriptions and poisoning is notable, suggesting that increased odds of poisoning following receipt of a benzodiazepine prescription is more attributable to benzodiazepines than to underlying conditions leading to a prescription for a benzodiazepine or an SSRI. Finally, our findings are pharmacologically informed, because we used the high statistical power of the MarketScan database to identify subgroups of patients receiving buprenorphine who may have been at elevated risk of treatment failure, such as those taking short-acting or higher-potency benzodiazepines.
Conclusions
Patients with opioid use disorder who are treated with buprenorphine have double the odds of a drug-related poisoning when benzodiazepines or Z-drugs are coprescribed, with short-acting and high-dose benzodiazepines conferring additional risk. Buprenorphine remains protective against overdose in this high-risk population. Our findings contribute to emerging evidence that patients with opioid use disorder who receive buprenorphine may benefit from cessation or dose reduction of benzodiazepines or Z-drugs while remaining on buprenorphine.