Skip to main content

Comparing the effects of decreasing prescription opioid shipments and the release of an abuse deterrent OxyContin formulation on opioid overdose fatalities in WV: an interrupted time series study

Abstract

Introduction

The 2010 release of an abuse deterrent formulation (ADF) of OxyContin, a brand name prescription opioid, has been cited as a major driver for the reduction in prescription drug misuse and the associated increasing illicit opioid use and overdose rates. However, studies of this topic often do not account for changes in supplies of other prescription opioids that were widely prescribed before and after the ADF OxyContin release, including generic oxycodone formulations and hydrocodone. We therefore sought to compare the impact of the ADF OxyContin release to that of decreasing prescription opioid supplies in West Virginia (WV).

Methods

Opioid tablet shipment and overdose data were extracted from The Washington Post ARCOS (2006–2014) and the WV Forensic Drug Database (2005–2020), respectively. Locally estimated scatterplot smoothing (LOESS) was used to estimate the point when shipments of prescription opioids to WV began decreasing, measured via dosage units and morphine milligram equivalents (MMEs). Interrupted time series analysis (ITSA) was used to compare the impact LOESS-identified prescription supply changes and the ADF OxyContin release had on prescription (oxycodone and hydrocodone) and illicit (heroin, fentanyl, and fentanyl analogues) opioid overdose deaths in WV. Model fit was compared using Akaike Information Criteria (AIC).

Results

The majority of opioid tablets shipped to WV from 2006 to 2014 were generic oxycodone or hydrocodone, not OxyContin. After accounting for a 6-month lag from ITSA models using the LOESS-identified change in prescription opioid shipments measured via dosage units (2011 Q3) resulted in the lowest AIC for both prescription (AIC = -188.6) and illicit opioid-involved overdoses (AIC = -189.4), indicating this intervention start date resulted in the preferred model. The second lowest AIC was for models using the ADF OxyContin release as an intervention start date.

Discussion

We found that illicit opioid overdoses in WV began increasing closer to when prescription opioid shipments to the state began decreasing, not when the ADF OxyContin release occurred. Similarly, the majority of opioid tablets shipped to the state for 2006–2014 were generic oxycodone or hydrocodone. This may indicate that diminishing prescription supplies had a larger impact on opioid overdose patterns than the ADF OxyContin release in WV.

Background

Opioid-involved overdose deaths are a major public health problem in the United States (US), with more than 644,000 fatal overdoses occurring from 1999 to 2021 [1]. The opioid epidemic has been characterized by multiple waves of overdoses associated with different drug classes and routes of administration [2, 3]. The first wave was associated with overdoses due to prescription opioid medications, such as oxycodone and its brand name extended release (ER) formulation OxyContin [4]. Beginning in the late 1990’s, these and other prescription opioids were prescribed and dispensed at increasing rates throughout the US. As a result, fatal overdoses of prescription opioids increased in tandem [2, 5].

Actions taken to decrease rates of prescription medication diversion, misuse, and overdose included efforts to reduce opioid prescribing rates [6] and prescription reformulations aiming to restrict injecting or snorting tablets, including the August 2010 release of an “abuse-deterrent formulation” (ADF) of OxyContin [7]. These targeted measures largely succeeded in decreasing rates of opioid overdoses involving prescribed opioids. However, they also may have had the unintended consequence of diverting those suffering from untreated opioid use disorder (OUD) from prescription opioids to cheaper and more available dangerous illicit alternatives, such as heroin and fentanyl [2, 8].

The ADF OxyContin release in August 2010 was widely promoted as the solution to diversion and misuse of prescription opioids. Immediately following the ADF release, sales of the non-ADF OxyContin brand ceased and OxyContin was solely prescribed in the ADF formulation [9]. Previously published research suggests that this formulation change was one reason that individuals transitioned from prescription to illicit opioid use [10, 11] and that it was the primary reason for an uptick in subsequent heroin use [12]. This hypothesis is supported by decreasing rates of OxyContin misuse and overdose after the ADF release [11]. However, a recent analysis suggests that most individuals misusing OxyContin simply switched to generic ER oxycodone following the ADF release and that falling rates of generic oxycodone prescriptions were more predictive (in comparison to the ADF OxyContin release) of subsequent increases in illicit opioid overdose at the state level [13]. Moreover, some previous studies on this topic have not assessed supplies of commonly prescribed opioids such as hydrocodone in their analyses [14], an opioid which was prescribed and misused at rates comparable to OxyContin before its ADF release [15,16,17].

Given the importance of supply side drivers of overdose [2], omitting widely prescribed opioids from analyses limits current understanding of the transition between prescription and illicit opioid overdose in the US. Thus, using fatal opioid overdose data from West Virginia (WV), we examined the impact of decreasing opioid shipments to WV on the transition from prescription to illicit opioid overdoses in the state. However, no fixed date for the start of declining shipments was available and current time series methods require a fixed starting point. Using a data-driven approach, we used locally estimated scatterplot smoothing (LOESS) regression to identify an approximate start date for the decline in oxycodone and hydrocodone tablet shipments, measured both via dosage units (tablets) and morphine milligram equivalents (MMEs), to WV given naturally occurring quarterly variation in data. These points were used to inform an interrupted time series analysis (ITSA) of fatal prescription and illicit opioid overdoses in the state. ITSA assesses the impact of public health events by quantifying trends before and after an intervention with a known start date [18]. We also assessed the impact of the August 2010 ADF OxyContin release for comparison to our LOESS-informed analysis. In addition to elucidating the role of prescription opioid supply changes in the transition from prescription to illicit opioid overdoses in WV, this study provides a plausible framework for using ITSA to assess the impact of public health events with no intervention start date.

Methods

Data sources

Data on opioid prescription shipments to WV for 2006–2014 were obtained from a subset of the Drug Enforcement Agency’s (DEA) Automation of Reports and Consolidated Orders System (ARCOS) data publicly available through The Washington Post [19]. ARCOS tracks the flow of all schedule I/II and select schedule III/IV substances through their manufacture and subsequent distribution to points of dispersion (i.e., retail pharmacies, hospitals, practitioners, etc.). The Washington Post ARCOS subset contains data on individual oxycodone and hydrocodone tablet shipments or dosage units, including the addresses of each shipment manufacturer/distributor and recipient, date of shipment, number of dosages (i.e., tablets) in each shipment, strength of each dose in milligrams, and morphine milligram equivalents (MME) conversion factors for each shipment (1 for hydrocodone, 1.5 for oxycodone). While this dataset only contains information on oxycodone and hydrocodone tablet shipments, The Washington Post reports that the prescription opioids excluded were shipped and diverted for misuse in much smaller quantities throughout the period reported [20]; this is congruent with another published research finding that the majority of dispensed opioids are oxycodone or hydrocodone at the state level [21]. Quarterly dosage units were calculated using an ARCOS dataset column corresponding to the number of tablets in each shipment, while quarterly MMEs were calculated using the formula MME = Quantity × Strength × Conversion Factor (Supplemental Table 1) [21,22,23].

WV opioid-involved overdose death data for 2005–2020 were obtained from a forensic drug database (FDD) maintained at West Virginia University through an agreement with the WV Office of the Chief Medical Examiner (OCME). The WV OCME uses a state-level, centralized death investigation system which includes comprehensive drug screenings and toxicology testing for suspected drug deaths [24]. Counts of opioid-involved overdose deaths from January 2005 to December 2020 were aggregated to the quarterly level for drug-related deaths involving prescription (oxycodone or hydrocodone) or illicit opioids (heroin or a synthetic opioid other than methadone, including fentanyl, fentanyl analogs, 4-anpp, and u-47,700). While overdoses involving hydrocodone or oxycodone often occur in individuals who obtained them illicitly (i.e., through diversion), these are labeled as “prescription” here as at the time these tablets are manufactured licitly for prescription purposes in the US. Similarly, although some licitly-manufactured fentanyl is diverted, the majority of fentanyl involved in synthetic opioid overdoses is illicitly manufactured outside of the US, particularly after 2013 [25,26,27]. Data from WV has shown that from 2015 to 2017, only 1.7% of decedents in whom fentanyl was involved in the deaths had a prescription for the drug, compared to almost 24% of the fentanyl related deaths from 2005 to 2014 [28]. Proportions of deaths involving prescription or illicit opioids were calculated by dividing the quarterly aggregate of either category by the total number opioid-involved deaths in each quarter.

Statistical analysis

All statistical analyses were performed in RStudio version 4.2.2 [29]. ITSA was used to assess the impact of decreasing opioid prescription shipments to WV on the proportion of opioid overdose deaths associated with prescription and illicit opioids [18, 30]. ITSA is a robust statistical approach in which the impact of an intervention is measured using segmented linear regression. Three potential interventions were investigated. First, two interventions denoted peak tablet shipments to WV, measured via both dosage units and MMEs. As shipments varied by quarter with no defined point when the decline began, we identified a plausible decline start point using LOESS of quarterly ARCOS data. LOESS fits weighted least squared regression to data in several independent variable intervals and requires no global function, providing clear graphical representations of non-linear relationships. As a result, LOESS is often used to identify inflection points (i.e., changes in slope) in non-linear data [31, 32], including in time series data of opioid prescriptions [33] and overdoses [34]. Next, an intervention for the introduction of the ADF OxyContin (released in August 2010 or 2010 Q3) was informed by previous literature; non-ADF OxyContin prescriptions ceased the same month that the ADF was released [9]. As there is likely a temporal lag between changes in prescription opioid shipments and related variations in opioid overdose rates, we lagged each intervention by transition periods of three and six months (i.e., one and two quarters, respectively) and tested each lag/transition period in separate models. This approach has been used in previous time series studies of prescription opioid supply changes during this timeframe [35,36,37,38].

An ITSA intervention may be modeled using the equation:

$${y_t} = {\beta _0} + {\beta _1}t + {\beta _2}P + {\beta _3}D + \varepsilon$$

where yt is an outcome of interest (e.g. quarterly proportion of opioid overdose deaths associated with prescription or illicit opioids), \({\beta }_{0}\) is the model intercept, t is time, \(P\) is a variable representing time since the intervention (zero before the intervention, slope of one afterwards), and \(D\) is a dummy variable representing the immediate effect of the intervention [39, 40]. \({\beta }_{1}\), \({\beta }_{2}\), and \({\beta }_{3}\) represent the pre-intervention slope, the sustained post-intervention effect (i.e., a slope change impact known as a “ramp” variable), and the immediate post-intervention effect of an intervention (i.e., a “step-change” impact), respectively. Finally, error is denoted by \(?\) and may include autoregressive integrated moving average (ARIMA) terms when data violate the linear regression assumption of data independence (i.e., the data is serially correlated). ARIMA models include lagged values of a time series’ dependent variable and/or its error terms and is recommended for use in ITSA when data are not independent [18, 41]. ARIMA terms were fit to opioid overdose data via inspecting autocorrelation and partial autocorrelation plots. ITSA and ARIMA parameters for each intervention were included or excluded based on the need for control of serial correlation, preservation of model parsimony, and minimization of Akaike Information Criteria (AIC) [30, 42]. Since lagged intervention points were included to account for slow changes in prescription drug supplies, and the subsequent impact on fatal overdose would be expected to be gradual rather than immediate on a statewide scale, step change (i.e., immediate impact) variables were not included [41]; measuring slope changes only is an approach similar to published research on this topic [7, 14]. Final ITSA models were assessed for proper fit via inspection of each model’s ACF and PACF plots, as well as inspecting the significance of each model’s Ljung-Box statistic (with a non-significant value considered a properly fitting model) [43]. The AIC of each final ITSA model was used for model comparison. To compare models using AIC, we abided by the convention that when comparing two models, the model with lower AIC is better fit and that a difference of two or more AIC units is meaningful [44].

Results

From 2005 to 2020, a total of 9419 opioid-involved overdose deaths were identified in the WV FDD. The proportion of deaths involving illicit opioids (synthetic opioids including fentanyl or heroin) was 0.48, while the proportion involving prescription opioids (oxycodone or hydrocodone) was 0.37; the majority of remaining fatal overdoses involved methadone (data not shown). Graphical representation of the quarterly proportion of opioid-involved overdose deaths involving prescription or illicit opioids is presented in Fig. 1. During the first two quarters of the study period, illicit opioid overdoses occurred at a rate comparable to prescription opioids; this is potentially related to a multi-state fentanyl overdose outbreak that occurred from 2005 to 2007 [45]. After this, prescription opioid overdoses occurred at a rate greater than those associated with illicit opioids until approximately 2015; at this time (2015 and later), overdose patterns changed substantially, with the majority of overdoses being associated with illicit opioids.

Fig. 1
figure 1

The quarterly proportion of opioid overdoses in WV associated with prescription and illicit opioids. a Prescription opioid overdoses were defined as those associated with oxycodone or hydrocodone, while illicit overdoses were defined as those involving heroin and synthetic opioids other than methadone, including fentanyl, fentanyl analogues, 4-anpp, and u-47,700. Data from the West Virginia Forensic Drug Database, which compiles data from the West Virginia Office of the Chief Medical Examiner

LOESS regression of 2006–2014 ARCOS data indicated that maximum quarterly shipments of oxycodone and hydrocodone tablets occurred in 2011 Q3 when measured via dosage units (i.e., number of tablets) and 2012 Q4 when measured via MMEs. Graphical representation of quarterly dosage units and MMEs are presented in Fig. 2 along with estimated peak total prescription shipments. From their peaks to the end of the available ARCOS data (2014 Q4), dosage unit and MME shipments decreased 15.2% and 9.7%, respectively. Measured via dosage units, hydrocodone was shipped to WV in highest quantities, followed by non-ER oxycodone and brand name oxycodone tablets including OxyContin. Measured via MMEs, hydrocodone was shipped to WV in highest quantities until approximately 2012, when it was surpassed by non-ER oxycodone. OxyContin was shipped in third-highest quantities (measured via MMEs) after approximately 2007. The opioid products included in each category are available in Supplemental Tables 26.

Fig. 2
figure 2

Quarterly opioid tablet shipments to West Virginia, both total and by individual opioid product and estimated change point in total opioid shipments indicated by red dot. a Data are presented measured via dosage units and morphine milligram equivalents (MMEs) and are smoothed using locally estimated scatterplot smoothing (LOESS) regression to allow for the visualization of overall trends. Peak total quarterly dosage units and MMEs were identified via LOESS and are denoted using a red dot. Data are form the Drug Enforcement Agency’s Automation of Reports and Consolidated Orders System (ARCOS) database and were obtained from The Washington Post

Graphical representation of each ITSA model with a 6-month lag and a corresponding counterfactual (i.e., no intervention) scenario is presented in Fig. 3; the blue dotted line indicates the intervention date of interest while the grey shaded area denotes a transition period of two quarters post-intervention. For all ITSA models, serial correlation was adequately controlled with an AR [1] ARIMA term (Tables 1 and 2). ITSA of illicit opioids was best modeled using only a ramp function, indicating there was no pre-intervention trend present (Table 2). For both prescription and illicit opioids, dosage units-informed ITSA models had lowest AIC, indicating best model fit; dosage-unit informed models with a 3-month lag had AIC of -189.2 and − 188.7 for prescription and illicit opioids, respectively, while models with a 6-month lag had AIC of -188.6 and − 189.4 for prescription and illicit opioids, respectively. For both 3- and 6-month lagged ITSA models of either prescription or illicit opioids, MME-informed models showed second lowest AIC and ADF OxyContin release third lowest AIC. For 3-month lagged models, dosage units-informed models were more than two AIC units lower than ADF release and peak MME-informed models, while MME- and ADF OxyContin-informed models were not meaningfully different (less than two AIC unit difference). For six-month lagged ITSA models, AIC difference between the three models was greater than two units for both prescription (peak dosage unit-informed model 3.1 and 10.3 units lower than ADF release and peak MME-informed models, respectively; Table 1) and illicit opioids (peak dosage units-informed model 3.8 and 8.7 units lower than ADF release and peak MME-informed models, respectively; Table 2), suggesting meaningful differences between each model’s performance. Intervention points informed using LOESS-identified changes in hydrocodone (2009 Q1 for both dosage units and MMEs; Supplementary Tables 7 and 8) and oxycodone (2014 Q1 for dosage units, 2014 Q2 for MMEs; Supplementary Tables 9 and 10) did not result in meaningfully better-fitting models (compared to the total opioid shipment informed models shown in Tables 1 and 2) for either prescription or illicit opioid overdose deaths.

Table 1 Interrupted time series results of opioid-involved overdoses in West Virginia involving prescription opioids.a
Table 2 Interrupted time series results of opioid-involved overdoses in West Virginia involving illicit opioids.a

Discussion

Using ITSA, this study compared the impact of decreasing opioid tablet shipments, measured via dosage units (tablets) and MMEs, to the release of ADF OxyContin, which is generally acknowledged as the factor initiating a transition from prescription to illicit opioid use in the US [11, 12, 46]. Our findings suggest that in WV, overdose patterns began changing closer to the time when prescription opioid shipments (measured via dosage units) began decreasing in 2011 Q3. We also accounted for three- and six-month lag/transition periods (i.e., one and two quarter lag/transition periods, respectively) between shipment to observable effect on fatal overdose rates, indicating the impact observed via ITSA was a full 1.25 to 1.5 years after the ADF OxyContin release.

Fig. 3
figure 3

Graphical representation of interrupted time series analysis (ITSA)a of the proportion quarterly opioid-involved overdose deaths involving prescription and illicit opioidsb

a Upper and lower sub-figures represent ITSA of prescription and illicit opioid-involved overdose rates, respectively, while each column represents a unique ITSA intervention. Red lines represent estimated intervention impacts while dotted red lines represent estimated counterfactual (i.e., no intervention) trends. Blue dotted lines represent intervention start dates while grey shaded areas represent a six-month (two-quarter) transition period after which the intervention impact is theorized to have begun; a three-month (one-quarter) lag/transition period, which was modeled and tested separately, is not shown

b Prescription opioid overdoses were defined as those associated with oxycodone or hydrocodone, while illicit opioid overdoses were defined as those involving heroin and synthetic opioids other than methadone, including fentanyl, fentanyl analogues, 4-anpp, and u-47700. Data are from the West Virginia Forensic Drug Database, which compiles data from the West Virginia Office of the Chief Medical Examiner

Changes in prescription opioid supplies have a measurable impact on use of both prescription and illicit opioids [2, 47,48,49]. Therefore, it is possible that decreasing prescription opioid supplies in WV during our study period contributed to increasing rates of overdoses involving illicit opioids. From 2006 to 2011, the WV opioid prescription dispensing rate was the highest of any US state [50]. However, WV opioid prescription rates, measured either via dosage units or MMEs, decreased more quickly than the US average during our study period [6, 50]. These drastic supply changes may explain our finding that the ITSA study based on overall opioid shipments to WV (measured via dosage units) was the preferred model for both prescription and illicit overdose rates in our study. Our data indicated that the majority of prescription opioid tablets shipped to the state were not OxyContin (Fig. 2); this is congruent with national data showing that the majority of oxycodone prescribed post-ADF OxyContin release was generic and therefore prone to misuse [13]. Thus, those who used OxyContin in WV likely had many other opioid options to begin using following the release of its ADF formulation. Moreover, those with OUD often prefer immediate release prescription opioid formulations as opposed to ER formulations [51]. This may be why the ADF OxyContin release did not fit our data well as a dosage unit-based intervention.

Soon after its release, the ADF OxyContin formulation was cited as a major contributor to subsequent increases in illicit opioid use and overdose [46]. However, analyses since the 2019 ARCOS data release by The Washington Post support our conclusion that changing supplies of other prescription opioids had a greater influence. For instance, Zhang and Guth assessed ARCOS, substance use data from The National Survey on Drug Use and Health (NSDUH), and opioid mortality data and found that the majority of OxyContin users in their sample transitioned to generic oxycodone after the ADF release [13]. The authors also noted that heroin mortality was highest in states with previously high generic oxycodone use and that illicit opioid overdose rates began increasing nearly two years after the ADF release. Our findings support and expand on those from Zhang and Guth by using ITSA, quarterly as opposed to annual data, and by including data on shipments of hydrocodone, a prescription opioid that was prescribed and misused at rates similar to OxyContin and oxycodone before the ADF OxyContin release [15,16,17].

While we used ARCOS data to assess prescription opioid supplies throughout our study period, analyses of other data sources support our conclusions. For example, NSDUH data indicate those using misusing OxyContin before its ADF release had 58% lower odds of heroin initiation than those using other prescription opioids [52]. Similarly, in a large sample of individuals screened for substance misuse, oxymorphone and buprenorphine use rates increased after the ADF OxyContin release while heroin use rates did not change significantly [10]. Another study using linked health insurance and National Death Index data did not find an overall effect on fatal and non-fatal overdoses; the study did find a small decrease in OxyContin overdose rates [9]. Finally, a 2020 US Food and Drug Administration (FDA) joint meeting of the Drug Safety and Risk Management Advisory Committee and the Anesthetic and Analgesic Drug Products Advisory Committee used multiple data sources to assess the national post-market impact of the ADF OxyContin release. Among the meeting’s findings were that while non-oral OxyContin use decreased after the ADF release, oral use increased significantly [36] and misuse of hydrocodone and other schedule II opioids increased significantly relative to OxyContin [37].

It is worthwhile to briefly discuss why decreasing prescription opioid supplies may have driven a mass transition to illicit opioid use. Previously published qualitative studies indicate that prescription opioid users who began using heroin during our study period did so because of the latter’s decreased cost relative to former, as well as ease of access to heroin [53]. While trends in national illicit opioid supplies are difficult to accurately quantify relative to prescription opioids, Customs and Border Patrol seizure data suggest that illicit heroin shipments to the US increased by more than 50% between 2012 and 2015; fentanyl seizures increased dramatically from 2015 to 2017 [54]. It has been suggested that these two trends (decreasing prescription opioids supplies and increasing illicit opioids supplies) are tied through supply-side economic principles; as restrictive or prohibitive constraints are placed on an addictive substance, availability of their illicit counterparts will increase to meet the pre-existing demand [8]. Of particular concern is the unpredictable variation in these illicit opioids’ potency and mixing with non-opioid substances, such as xylazine [55] and methamphetamine [24], which do not respond to opioid overdose reversal medications such as naloxone.

In addition to elucidating the transition from prescription to illicit opioid overdose in WV, this study expands on ITSA literature seeking to identify intervention dates using data that is related to, yet separate from, a time series of interest. Notably, Gilmour et al. used previously published survey data to identify a plausible ITSA start date of the Australian heroin shortage [42]. Similarly, Lopez Bernal et al. used a widely accepted definition for the beginning of an economic recession (the point at which gross domestic product growth rate is negative compared to previous quarter) to assess the impact of the late 2000’s financial crisis on suicide rates in Spain [56]. While these studies provided innovative approaches, the methods used to identify interventions were statistically descriptive. Using the inferential LOESS regression approach in our study, future studies might more accurately determine an intervention start date, when one is not easily defined, for ITSA study use.

This study has several strengths. For example, to the author’s knowledge, we are the first to demonstrate the feasibility of LOESS regression to determine a plausible intervention start date for use in an ITSA study of substance use data. As previous substance use research has identified the difficulty in studying interventions with no known start date [42], this is an interesting addition to current ITSA literature. Additionally, our study assessed medical examiner data from WV, which has a highly specific drug death investigation system relative to other states [57]. Our study also has several limitations. First, although ITSA is a powerful statistical study design useful in many situations in which a public health intervention has no control group, it remains an ecologic study design that cannot infer causality. Despite this, we believe the methodology in this study provides a robust approach towards strengthening evidence for a specific intervention’s impact since it uses ITSA to quantify the impact of several intervention points. Second, while changes in opioid tablet shipments occurred gradually over time, ITSA can only incorporate interventions at a single point in time. However, we tested each intervention at pre-specified, literature informed lag/transition periods (i.e., three and six-months) to account for the slow decline in opioid tablet distribution. Third, given urban/rural differences in prescription opioid misuse rates, our results may not be generalizable outside of WV, a largely rural state. Fourth, the medical examiner’s data used in this study relies on toxicology reports that cannot differentiate between formulations of the same drug. We therefore cannot assess which formulation of oxycodone (OxyContin or generic) or hydrocodone contributed to overdose rates throughout our study period. There is also the possibility for some missing or inaccurate data entries into the Forensic Drug Database, although this is believed to be minimal. Finally, ARCOS is limited to hydrocodone and oxycodone; while data suggest these were the primary opioids contributing to substance use disorder during the early years of the opioid epidemic in WV, other prescription opioids played a role.

Conclusion

These results suggest that the transition from prescription to illicit opioid overdose in WV may have been affected to a greater extent by decreasing rates of prescription opioid shipments, as compared to the release of ADF OxyContin previously reported in studies of national data. The large quantity of hydrocodone shipped to WV (relative to OxyContin) might also have reduced the apparent impact on overdose deaths from the ADF OxyContin release. Future research should explore the impacts of supply-side changes in the availability of other prescription opioids in WV on not only drug-related deaths, but also on other outcomes, including substance use, substance use treatment, and related co-morbidities such as acute hepatitis C infection.

Data availability

The ARCOS and West Virginia forensic drug datasets used in this study are available from the authors upon reasonable request.

Abbreviations

ACF:

Autocorrelation function

ADF:

Abuse-deterrent formulation

AIC:

Akaike Information Criteria

ARCOS:

Automation of Reports and Consolidated Orders System

ARIMA:

Autoregressive integrated moving average

DEA:

Drug Enforcement Agency’s

ER:

Extended release

FDD:

Forensic drug database

LOESS:

Locally estimated scatterplot smoothing

ITSA:

Interrupted time series analysis

MME:

Morphine milligram equivalent

NSDUH:

National Survey on Drug Abuse and Health

OCME:

Office of the Chief Medical Examiner

OUD:

Opioid use disorder

PACF:

Partial autocorrelation function

US:

United States

WV:

West Virginia

References

  1. Centers for Disease Control and Prevention. Wide-ranging Online Data for Epidemiologic Research (WONDER). 2022 [cited 2022 Oct 19]. Available from: https://wonder.cdc.gov

  2. Ciccarone D. The triple wave epidemic: supply and demand drivers of the US opioid Overdose crisis. Vol. 71, International Journal of Drug Policy. Elsevier B.V.; 2019. p. 183–8.

  3. Jenkins RA. The fourth wave of the US opioid epidemic and its implications for the rural US: a federal perspective. Prev Med (Baltim). 2021;152:106541.

    Article  Google Scholar 

  4. Van Zee A. The Promotion and Marketing of OxyContin : commercial triumph, Public Health Tragedy. Health Policy and Ethics. 2009;99(2):221–7.

    Google Scholar 

  5. Ciccarone D. The rise of illicit fentanyls, stimulants and the fourth wave of the opioid Overdose crisis. Vol. 34, Current Opinion in Psychiatry. Lippincott Williams and Wilkins; 2021. p. 344–50.

  6. Schieber LZ, Guy GP, Seth P, Young R, Mattson CL, Mikosz CA, et al. Trends and patterns of Geographic Variation in Opioid Prescribing practices by State, United States, 2006–2017. JAMA Netw Open. 2019;2(3):e190665.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Sessler NE, Downing JM, Kale H, Chilcoat HD, Baumgartner TF, Coplan PM. Reductions in reported deaths following the introduction of extended-release oxycodone (OxyContin) with an abuse-deterrent formulation. Pharmacoepidemiol Drug Saf. 2014;23(12):1238–46. Available from: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pubmed/24916486

  8. Beletsky L, Davis CS. Today’s fentanyl crisis: prohibition’s Iron Law, revisited. International Journal of Drug Policy. 2017;46:156–9.

    Article  PubMed  Google Scholar 

  9. Beachler DC, Hall K, Garg R, Banerjee G, Li L, Boulanger L, et al. An evaluation of the Effect of the OxyContin®Reformulation on Unintentional Fatal and non-fatal Overdose. Clinical Journal of Pain. 2022;

  10. Cassidy TA, DasMahapatra P, Black RA, Wieman MS, Butler SF. Changes in prevalence of prescription opioid abuse after introduction of an abuse-deterrent opioid formulation. Pain Medicine. 2014;15(3):440–51.

    Article  PubMed  Google Scholar 

  11. Cicero TJ, Ellis MS. Abuse-deterrent formulations and the prescription opioid abuse epidemic in the United States. JAMA Psychiatry. 2015;72(5):424.

    Article  PubMed  Google Scholar 

  12. Evans WN, Lieber EMJ, Power P. How the Reformulation of OxyContin Ignited the Heroin Epidemic. Rev Econ Stat. 2019;101(1):1–15. Available from: https://direct.mit.edu/rest/article/101/1/1/58660/How-the-Reformulation-of-OxyContin-Ignited-the

  13. Zhang S, Guth D. The OxyContin Reformulation Revisited: New Evidence From Improved Definitions of Markets and Substitutes. 2021 [cited 2023 Jan 8];1–49. Available from: http://arxiv.org/abs/2101.01128

  14. Coplan PM, Kale H, Sandstrom L, Landau C, Chilcoat HD. Changes in oxycodone and heroin exposures in the National Poison Data System after introduction of extended-release oxycodone with abuse‐deterrent characteristics. Pharmacoepidemiol Drug Saf. 2013;22(12):1274–82. Available from: https://0-onlinelibrary-wiley-com.brum.beds.ac.uk/doi/https://0-doi-org.brum.beds.ac.uk/10.1002/pds.3522

  15. Cicero TJ, Inciardi JA, Muñoz A. Trends in Abuse of OxyContin® and Other Opioid Analgesics in the United States: 2002–2004. J Pain. 2005;6(10):662–72. Available from: https://0-linkinghub-elsevier-com.brum.beds.ac.uk/retrieve/pii/S1526590005006565

  16. Cicero TJ, Surratt H, Inciardi JA, Munoz A. Relationship between therapeutic use and abuse of opioid analgesics in rural, suburban, and urban locations in the United States. Pharmacoepidemiol Drug Saf. 2007;16(8):827–40.

    Article  CAS  PubMed  Google Scholar 

  17. Kenan K, Mack K, Paulozzi L. Trends in prescriptions for oxycodone and other commonly used opioids in the United States, 2000–2010. Open Med. 2012;6(2):e41-7. Available from: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pubmed/23696768

  18. Lopez Bernal J, Cummins S, Gasparrini A. Interrupted time series regression for the evaluation of public health interventions: a tutorial. Int J Epidemiol. 2016;46(1):dyw098. Available from: https://0-academic-oup-com.brum.beds.ac.uk/ije/article-lookup/doi/https://0-doi-org.brum.beds.ac.uk/10.1093/ije/dyw098

  19. The Washington Post. The Washington Post. 2020 [cited 2022 Oct 19]. Drilling into the DEA’s pain pill database. Available from: https://www.washingtonpost.com/graphics/2019/investigations/dea-pain-pill-database/?utm_term=.927f44f26372&itid=lk_inline_manual_2

  20. Rich S, Diez MS, Vongkiatkajorn K. The Washington Post. 2019 [cited 2023 Feb 21]. How to download and use the DEA pain pills database. Available from: https://www.washingtonpost.com/national/2019/07/18/how-download-use-dea-pain-pills-database/

  21. Winstanley EL, Zhang Y, Mashni R, Schnee S, Penm J, Boone J, et al. Mandatory review of a prescription drug monitoring program and impact on opioid and benzodiazepine dispensing. Drug Alcohol Depend. 2018;188:169–74. Available from: https://0-linkinghub-elsevier-com.brum.beds.ac.uk/retrieve/pii/S0376871618302369

  22. Centers for Disease Control and Prevention. Calculating Total Daily Dose of Opioids for Safer Dosage. [cited 2022 Oct 19]. Available from: https://www.cdc.gov/drugoverdose/pdf/calculating_total_daily_dose-a.pdf

  23. Sedney CL, Khodaverdi M, Pollini R, Dekeseredy P, Wood N, Haggerty T. Assessing the impact of a restrictive opioid prescribing law in West Virginia. Subst Abuse Treat Prev Policy. 2021;16(1):1–12.

    Article  Google Scholar 

  24. Dai Z, Abate MA, Groth CP, Rucker T, Kraner JC, Mock AR, et al. Fentanyl and other opioid involvement in methamphetamine-related deaths. Am J Drug Alcohol Abuse. 2022;48(2):226–34.

    Article  CAS  PubMed  Google Scholar 

  25. Barry CL. Fentanyl and the Evolving Opioid Epidemic: What Strategies Should Policy Makers Consider? Psychiatric Services. 2018;69(1):100–3. Available from: https://0-doi-org.brum.beds.ac.uk/10.1176/appi.ps.201700235

  26. Felbab-Brown V. Fentanyl and geopolitics Controlling opioid supply from China. 2020 Jan. Available from: www.brookings.edu/about-us/annual-report/.

  27. O’Donnell J, Gladden RM, Kariisa M, Mattson CL. Using death scene and toxicology evidence to define involvement of heroin, pharmaceutical morphine, illicitly manufactured fentanyl and pharmaceutical fentanyl in opioid overdose deaths, 38 states and the District of Columbia, January 2018–December 2019. Addiction. 2022;117(5):1483–90. Available from: https://0-onlinelibrary-wiley-com.brum.beds.ac.uk/doi/https://0-doi-org.brum.beds.ac.uk/10.1111/add.15768

  28. Dai Z, Abate MA, Smith GS, Kraner JC, Mock AR. Fentanyl and fentanyl-analog involvement in drug-related deaths. Drug Alcohol Depend. 2019;196:1–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. RStudio Team. RStudio: Integrated Development for R. Boston, MA: RStudio, Inc.; 2022 [cited 2022 Mar 29]. Available from: http://www.rstudio.com/

  30. Schaffer AL, Dobbins TA, Pearson SA. Interrupted time series analysis using autoregressive integrated moving average (ARIMA) models: a guide for evaluating large-scale health interventions. BMC Med Res Methodol. 2021;21(1):58. Available from: https://0-bmcmedresmethodol-biomedcentral-com.brum.beds.ac.uk/articles/https://0-doi-org.brum.beds.ac.uk/10.1186/s12874-021-01235-8

  31. Jacoby WG. Electoral inquiry section Loess: a nonparametric, graphical tool for depicting relationships between variables. Electoral Studies. 2000;19. Available from: www.elsevier.com/locate/electstud

  32. Ryan SE, Porth LS. A tutorial on the piecewise regression approach applied to bedload transport data. 2007. Available from: https://www.fs.usda.gov/treesearch/pubs/27004

  33. Ahmedani BK, Peterson EL, Wells KE, Lanfear DE, Williams LK. Policies and events affecting prescription opioid use for non-cancer pain among an insured patient population. Pain Physician. 2014;17(3):205–16. Available from: http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pubmed/24850102

  34. Stevens JP, Wall MJ, Novack L, Marshall J, Hsu DJ, Howell MD. The critical care crisis of opioid overdoses in the United States. Ann Am Thorac Soc. 2017;14(12):1803–9.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Larochelle MR, Zhang F, Ross-Degnan D, Wharam JF. Rates of opioid dispensing and Overdose after introduction of abuse-deterrent extended-release oxycodone and withdrawal of Propoxyphene. JAMA Intern Med. 2015;175(6):978.

    Article  PubMed  Google Scholar 

  36. Mallama C. FDA review of PMR 3051-1 & 3051-3: Treatment Center Data to Assess the Impact of OxyContin Reformulation on Non-Oral and Overall Abuse of OxyContin. 2020 [cited 2023 Feb 19]. Available from: https://web.archive.org/web/20210306201006/https://www.fda.gov/media/141973/download

  37. Secora A. FDA Review of PMR 3051-2 & 3051-4: Poison Control Center Study and Opioid Overdose Study Using Administrative Claims Dat. 2020 [cited 2023 Feb 19]. Available from: https://web.archive.org/web/20210306201006/https://www.fda.gov/media/141973/download

  38. Severtson SG, Ellis MS, Kurtz SP, Rosenblum A, Cicero TJ, Parrino MW, et al. Sustained reduction of diversion and abuse after introduction of an abuse deterrent formulation of extended release oxycodone. Drug Alcohol Depend. 2016;168:219–29.

    Article  CAS  PubMed  Google Scholar 

  39. Nyugen M. Interrupted Time Series. In: A Guide on Data Analysis. 2022 [cited 2023 Jan 3]. Available from: https://bookdown.org/mike/data_analysis/interrupted-time-series.html

  40. Penfold RB, Zhang F. Use of interrupted time series analysis in evaluating health care quality improvements. Acad Pediatr. 2013;13(6 SUPPL.).

  41. Lopez Bernal J, Soumerai S, Gasparrini A. A methodological framework for model selection in interrupted time series studies. J Clin Epidemiol. 2018;103:82–91. Available from: https://0-linkinghub-elsevier-com.brum.beds.ac.uk/retrieve/pii/S0895435617314117

  42. Gilmour S, Degenhardt L, Hall W, Day C. Using intervention time series analyses to assess the effects of imperfectly identifiable natural events: a general method and example. BMC Med Res Methodol. 2006;6.

  43. LJUNG GM, BOX GEP. On a measure of lack of fit in time series models. Biometrika. 1978;65(2):297–303. Available from: https://0-academic-oup-com.brum.beds.ac.uk/biomet/article-lookup/doi/https://0-doi-org.brum.beds.ac.uk/10.1093/biomet/65.2.297

  44. Burnhamn KP, Anderson DR. Model selection and Multimodel Inference: a practical information-theoretic Approach. 2nd ed. New York, NY: Springer-Verlag New York, Inc.; 2002.

    Google Scholar 

  45. DEA Diversion Control Division. Fentanyl. 2016 [cited 2022 Sep 26]. Available from: https://www.dea.gov/sites/default/files/drug_chem_info/fentanyl.pdf

  46. Cicero TJ, Ellis MS, Surratt HL. Effect of abuse-deterrent formulation of OxyContin. New England Journal of Medicine. 2012;367(2):187–9.

    Article  CAS  PubMed  Google Scholar 

  47. Fischer B, Pang M, Jones W. The opioid mortality epidemic in North America: do we understand the supply side dynamics of this unprecedented crisis? Vol. 15, Substance Abuse: Treatment, Prevention, and Policy. BioMed Central Ltd.; 2020.

    Google Scholar 

  48. Greenfield VA, Paoli L. Research as due diligence: what can supply-side interventions accomplish and at what cost? International Journal of Drug Policy. 2017;41:162–3.

    Article  PubMed  Google Scholar 

  49. Werle N, Zedillo E. We can’t go cold Turkey: why suppressing drug markets endangers society. Journal of Law, Medicine and Ethics. 2018;46(2):325–42.

    Article  Google Scholar 

  50. Centers for Disease Control and Prevention NC for IP and C. U.S. Opioid Dispensing Rate Maps. 2021.

  51. Cicero TJ, Ellis MS, Kasper ZA. Relative preferences in the abuse of immediate-release versus extended-release opioids in a sample of treatment-seeking opioid abusers. Pharmacoepidemiol Drug Saf. 2017;26(1):56–62.

    Article  CAS  PubMed  Google Scholar 

  52. Wolff C, Dowd WN, Ali MM, McClellan C, Meinhofer A, Glos L, et al. The impact of the abuse-deterrent reformulation of extended-release OxyContin on prescription pain reliever misuse and heroin initiation. Addictive Behaviors. 2020;105.

  53. Mars SG, Bourgois P, Karandinos G, Montero F, Ciccarone D. ‘Every “never” i ever said came true’: transitions from opioid pills to heroin injecting. International Journal of Drug Policy. 2014;25(2):257–66.

    Article  PubMed  Google Scholar 

  54. Finklea K. Illicit Drug Flows and Seizures in the United States: What Do We [Not] Know? [Internet]. 2019. Available from: www.crs.gov.

  55. Sibbesen J, Abate MA, Dai Z, Smith GS, Lundstrom E, Kraner JC, et al. Characteristics of xylazine-related deaths in West Virginia—Xylazine‐related deaths. Am J Addict. 2023;32(3):309–13.

    Article  PubMed  Google Scholar 

  56. Lopez Bernal JA, Gasparrini A, Artundo CM, McKee M. The effect of the late 2000s financial crisis on suicides in Spain: an interrupted time-series analysis. Eur J Public Health. 2013;23(5):732–6.

    Article  PubMed  Google Scholar 

  57. Warner M, Hedegaard H. Identifying opioid Overdose deaths using vital Statistics Data. Am J Public Health. 2018;108(12):1587–9.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Funding

BH and GSS effort on this project was supported by NIH NIGMS award U54GM104942.

Author information

Authors and Affiliations

Authors

Contributions

EWL, CPG, BH, ELW, and GSS contributed to the conceptualization of this study. ZD aided in data extracted for West Virginia opioid-involved overdoses. ELW and MA provided expert opinion regarding opioid use and overdose. EWL performed all statistical analyses and drafted the manuscript for review. All authors approved this manuscript and the results presented within it before submission.

Corresponding author

Correspondence to Eric W. Lundstrom.

Ethics declarations

Ethics approval

This study was submitted to the Institutional Review Boards of West Virginia University, who determined it did not meet the definition of human subject research (IRB protocol number: 2211674154).

Patient consent for publication

Not applicable.

Competing interests

None declared.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

13011_2023_587_MOESM1_ESM.docx

Supplementary Material 1: Supplemental Table 1: Washington Post ARCOS data variables used in the calculation of dosage units and MMEs. 

Supplemental Table 2: Hydrocodone products shown in Fig. 2, measured via total dosage units and morphine milligram equivalents (MMEs) shipped to West Virginia for 2006–2014. 

Supplemental Table 3. Generic extended-release oxycodone products shown in Fig. 2, measured via total dosage units and morphine milligram equivalents (MMEs) shipped to West Virginia for 2006–2014. 

Supplemental Table 4: Generic non-extended-release oxycodone products shown in Fig. 2, measured via total dosage units and morphine milligram equivalents (MMEs) shipped to West Virginia for 2006–2014. 

Supplemental Table 5: OxyContin products shown in Fig. 2, measured via total dosage units and morphine milligram equivalents (MMEs) shipped to West Virginia for 2006–2014. Bolded rows denoted abuse-deterrent formulations. 

Supplemental Table 6: Brand name oxycodone products (other than OxyContin) shown in Fig. 2, measured via total dosage units and morphine milligram equivalents (MMEs) shipped to West Virginia for 2006–2014. 

Supplemental Table 7: Interrupted time series results of opioid-involved overdoses in West Virginia involving prescription opioids; informed via LOESS regression of hydrocodone tablet shipments to WV. 

Supplemental Table 8: Interrupted time series results of opioid-involved overdoses in West Virginia involving illicit opioids; informed via LOESS regression of hydrocodone tablet shipments to WV. 

Supplemental Table 9: Interrupted time series results of opioid-involved overdoses in West Virginia involving prescription opioids; informed via LOESS regression of oxycodone (including OxyContin) tablet shipments to WV. 

Supplemental Table 10: Interrupted time series results of opioid-involved overdoses in West Virginia involving illicit opioids; informed via LOESS regression of oxycodone (including OxyContin) tablet shipments to WV.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lundstrom, E.W., Dai, Z., Groth, C.P. et al. Comparing the effects of decreasing prescription opioid shipments and the release of an abuse deterrent OxyContin formulation on opioid overdose fatalities in WV: an interrupted time series study. Subst Abuse Treat Prev Policy 19, 4 (2024). https://0-doi-org.brum.beds.ac.uk/10.1186/s13011-023-00587-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13011-023-00587-2

Keywords