| Literature DB >> 34794949 |
James Wilton1, Younathan Abdia1,2, Mei Chong1, Mohammad Ehsanul Karim2,3, Stanley Wong1, Aaron MacInnes4,5, Rob Balshaw6, Bin Zhao1, Tara Gomes7,8,9, Amanda Yu1, Maria Alvarez1, Richard C Dart10,11, Mel Krajden1,12, Jane A Buxton13,2, Naveed Z Janjua1,2, Roy Purssell1,14.
Abstract
OBJECTIVE: To assess the association between long term prescription opioid treatment medically dispensed for non-cancer pain and the initiation of injection drug use (IDU) among individuals without a history of substance use.Entities:
Mesh:
Substances:
Year: 2021 PMID: 34794949 PMCID: PMC8600402 DOI: 10.1136/bmj-2021-066965
Source DB: PubMed Journal: BMJ ISSN: 0959-8138
Fig 1Study flowchart. IDEAs=Integrated Data and Evaluative Analytics data platform; PO=prescription opioid. *Opioid naive or acute use cohort contains the pool of potential opioid naive or acute use individuals available for matching to the chronic use category; see methods for more detailed definitions of opioid use episodes
Characteristics of study participants by prescription opioid use category, matched cohort (n=59 804). Data are number (%) of participants
| Chronic use | Episodic use | Acute use | Opioid naive (n=14 951) | |
|---|---|---|---|---|
| Matched variables (sex, birth year, year of index date, approximate age (±1 year), ethnic origin) | Approximately the same across categories* | |||
| Health authority (British Columbia) | ||||
| Vancouver Coastal | 2931 (19.6) | 3230 (21.6) | 4087 (27.3) | 4843 (32.4) |
| Fraser | 4904 (32.8) | 5220 (34.9) | 5228 (35.0) | 4998 (33.4) |
| Vancouver Island | 2666 (17.8) | 2549 (17.1) | 2366 (15.8) | 2095 (14.0) |
| Interior | 3216 (21.5) | 2564 (17.2) | 2216 (14.8) | 2130 (14.3) |
| Northern | 1234 (8.3) | 1388 (9.3) | 1054 (7.1) | 885 (5.9) |
| Urbanicity | ||||
| Urban | 12 644 (84.6) | 12 696 (84.9) | 12 999 (86.9) | 13 201 (88.3) |
| Rural | 2254 (15.1) | 2186 (14.6) | 1904 (12.7) | 1694 (11.3) |
| Missing | 53 (0.4) | 69 (0.5) | 48 (0.3) | 56 (0.4) |
| Material deprivation† | ||||
| 1 (least) | 2405 (16.1) | 2873 (19.2) | 3444 (23.0) | 3678 (24.6) |
| 5 (most) | 3295 (22.0) | 2969 (19.9) | 2439 (16.3) | 2394 (16.0) |
| Missing | 97 (0.7) | 87 (0.6) | 67 (0.5) | 72 (0.5) |
| Social deprivation† | ||||
| 1 (least) | 2222 (14.9) | 2496 (16.7) | 2648 (17.7) | 2671 (17.9) |
| 5 (most) | 4044 (27.1) | 3622 (24.2) | 3481 (23.3) | 3658 (24.5) |
| Missing | 97 (0.7) | 87 (0.6) | 67 (0.5) | 72 (0.5) |
| Chronic pain | 8452 (56.5) | 6905 (46.2) | 5232 (35.0) | 3197 (21.4) |
| Major mental illness | 2937 (19.6) | 2338 (15.6) | 1837 (12.3) | 1362 (9.1) |
| Alcohol use problem | 996 (6.7) | 686 (4.6) | 409 (2.7) | 255 (1.7) |
| Elixhauser comorbidity index† | ||||
| 0 | 7171 (48.0) | 8793 (58.8) | 10 206 (68.3) | 11 082 (74.1) |
|
| 2234 (14.9) | 1341 (9.0) | 838 (5.6) | 649 (4.3) |
Index date for measurement of variables was initiation date of long term use (for people in episodic and chronic use categories) or index date of the matched chronic use individual (for people in opioid naive and acute use categories). Individuals were matched on sex, birth year, and ethnic origin.
Data for sex, birth year, year of index date, approximate age (±1 year), and ethnic origin were similar across all categories (male (50.4%), female (49.6%); <1950s (17.4%), 1950-60s (53.0%), and >1970s (29.6%); 2000-04 (29.3%), 2005-09 (34.1%), and 2010-15 (36.6%); age <25 (7%), 25-34 (17%), 35-44 (23%), 45-54 (28%), and 55-65 (26%) years; and East Asian (2.2%), South Asian (5.6%), and other residents of British Columbia (92.2%), respectively).
All levels not shown for ease of presentation.
Number of participants, follow-up time, and rate of injection drug use (IDU) initiation by prescription opioid use category
| Category | No of participants | Follow-up (person years) | No of incident IDU initiations | Rate of IDU initiation | 5 year cumulative probability of IDU initiation (%) | |
|---|---|---|---|---|---|---|
| Total | Median (interquartile range) | |||||
| Full matched cohort* | 59 804 | 350 131 | 5.8 (3.0-9.4) | 1149 | 3.3 | 1.6 |
| Opioid naive | 14 951 | 87 939 | 5.8 (3.1-9.4) | 73 | 0.8 | 0.4 |
| Acute | 14 951 | 77 878 | 4.9 (2.4-8.1) | 98 | 1.3 | 0.7 |
| Long term (all) | 29 902 | 194 314 | 6.3 (3.4-10.0) | 978 | 5.0 | 2.6 |
| Episodic | 14 951 | 93 780 | 6.5 (3.5-10.0) | 248 | 2.6 | 1.3 |
| Chronic | 14 951 | 90 534 | 6.2 (3.2-9.8) | 730 | 8.1 | 4.0 |
| Average daily dose (chronic use only; morphine equivalents)† | ||||||
| <50 | 9887 | 58 499 | 5.9 (3.0-9.7) | 291 | 5.0 | 2.4 |
| 50-89 | 2574 | 15 259 | 6.0 (3.3-9.0) | 139 | 9.1 | 4.7 |
| 90-199 | 1695 | 11 187 | 7.1 (4.0-10.0) | 175 | 15.6 | 7.9 |
| ≥200 | 795 | 5589 | 7.9 (4.6-10.0) | 125 | 22.4 | 11.0 |
| Age (chronic use only; years)† | ||||||
| <25 | 1083 | 6953 | 6.6 (3.8-10.0) | 129 | 18.6 | 9.4 |
| 25-34 | 2502 | 15 703 | 6.4 (3.7-10.0) | 209 | 13.3 | 6.6 |
| 35-44 | 3382 | 23 540 | 7.8 (4.4-10.0) | 193 | 8.2 | 4.0 |
| 45-65 | 7984 | 44 338 | 5.4 (2.7-8.7) | 199 | 4.5 | 2.3 |
See methods for more details on opioid use categories.
Average daily dose calculated as the cumulative prescribed morphine equivalents during episode divided by the number of episode days covered by drug supply. Age measured at initiation of chronic opioid use.
Fig 2Cumulative incidence of injection drug use initiation by prescription opioid use category, age at initiation of chronic episode (chronic use category only), and average daily dose of chronic episode (chronic use category only). PO=prescription opioid; MEQ=morphine equivalents; IDU=injection drug use. See methods for more details on opioid use categories. Average daily dose calculated as the cumulative morphine equivalents dispensed during episode divided by number of episode days covered by drug supply
Association between prescription opioid use category and initiation of injection drug use in Cox models
| Prescription opioid use category* | Full matched cohort (primary analysis; n=59 804) | Matched cohort of individuals with history of chronic pain (sensitivity analysis; n=33 564) | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| No of individuals | 5 year cumulative | Bivariable (unweighted; HR (95% CI)) | Bivariable | Multivariable | No of individuals | 5 year cumulative | Bivariable (unweighted; HR (95% CI)) | Bivariable | Multivariable | ||
| Opioid naive | 14 951 | 0.4 | Reference | Reference | Reference | 8391 | 0.3 | Reference | Reference | Reference | |
| Acute | 14 951 | 0.7 | 1.5 (1.1 to 2.1) | 1.5 (1.0 to 2.0) | 1.4 (1.0 to 2.0) | 8391 | 0.5 | 1.8 (1.1 to 2.9) | 1.6 (1.0 to 2.7) | 1.6 (1.0 to 2.6) | |
| Episodic | 14 951 | 1.3 | 3.2 (2.4 to 4.1) | 2.8 (2.1 to 3.7) | 2.7 (2.1 to 3.6) | 8391 | 1.0 | 4.2 (2.7 to 6.3) | 3.3 (2.2 to 5.1) | 3.3 (2.1 to 5.0) | |
| Chronic | 14 951 | 4.0 | 9.7 (7.6 to 12.4) | 8.4 (6.5 to 11.0) | 8.4 (6.4 to 10.9) | 8391 | 3.4 | 12.6 (8.5 to 18.6) | 9.9 (6.7 to 14.8) | 9.7 (6.5 to 14.5) | |
HR=hazard ratio; aHR=adjusted hazard ratio; IPTW=inverse probability treatment weighting; IDU=injection drug use.
The multinomial logistic regression model used to calculate the generalised propensity scores for IPTW included the following variables: sex, age, ethnic origin, health authority, material deprivation, chronic pain (full cohort only), major mental illness, and alcohol use problems. The multivariable IPTW models were adjusted for the same variables as the multinomial logistic regression model used to generate propensity scores. All models used a robust sandwich covariate matrix to account for intracluster correlation within each matched group of four.
See methods on more details on opioid use categories.
Characteristics associated with initiation of injection drug use in the multivariable IPTW Cox model, full matched cohort (n=59 804)
| Adjusted hazard ratio (95% CI) | |
|---|---|
| Prescription opioid use category ( | |
| Acute | 1.4 (1.0 to 2.0) |
| Episodic | 2.7 (2.1 to 3.6) |
| Chronic | 8.4 (6.4 to 10.9) |
| Male sex ( | 1.6 (1.4 to 1.8) |
| Age ( | |
| <25 | 4.8 (3.7 to 6.2) |
| 25-34 | 3.3 (2.6 to 4.2) |
| 35-44 | 2.0 (1.6 to 2.6) |
| 45-54 | 1.1 (0.9 to 1.5) |
| Ethnic origin ( | |
| South Asian | 0.5 (0.3 to 0.7) |
| East Asian | 0.4 (0.2 to 0.9) |
| Material deprivation ( | |
| 2 | 1.3 (1.0 to 1.6) |
| 3 | 1.4 (1.1 to 1.8) |
| 4 | 1.4 (1.1 to 1.8) |
| 5 (most) | 1.6 (1.3 to 2.0) |
| Missing† | 2.1 (1.2 to 3.7) |
| Chronic pain | 0.8 (0.7 to 0.9) |
| Alcohol use problem | 2.7 (2.2 to 3.3) |
| Major mental illness | 1.6 (1.3 to 1.8) |
IPTW=inverse probability of treatment weighting.
Model also adjusted for health authority (no statistically significant associations). The multinomial logistic regression model used to calculate the generalised propensity scores for IPTW included all the variables included in the multivariable Cox model. The Cox model used a robust sandwich covariate matrix to account for intracluster correlation within each matched group.
See methods on more details on opioid use categories.
Missing material deprivation was retained in the model because it could be a proxy for homelessness.