If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Psychosocial and behavioral therapy in conjunction with medication for opioid use disorder: Patterns, predictors, and association with buprenorphine treatment outcomes
Corresponding author at: Institute for Health, Health Care Policy and Aging Research, Rutgers University, 112 Paterson St., New Brunswick, NJ 08901, United States of America.
Institute for Health, Health Care Policy, and Aging Research, Rutgers University, 112 Paterson St., New Brunswick, NJ 08901, United States of AmericaDepartment of Health Behavior, Society and Policy, Rutgers School of Public Health, 683 Hoes Lane West, Piscataway, NJ 08854, United States of America
Department of Psychiatry, New York State Psychiatric Institute, Columbia University Medical Center, 1051 Riverside Dr., New York, NY 10032, United States of America
Institute for Health, Health Care Policy, and Aging Research, Rutgers University, 112 Paterson St., New Brunswick, NJ 08901, United States of AmericaSchool of Social Work, Rutgers University, 120 Albany Street, Tower One - Suite 200, New Brunswick, NJ 08901, United States of America
Department of Psychiatry, New York State Psychiatric Institute, Columbia University Medical Center, 1051 Riverside Dr., New York, NY 10032, United States of AmericaDepartment of Epidemiology, Columbia University Mailman School of Public Health, 722 W. 168th St., New York, NY 10032, United States of America
Most buprenorphine patients receive no psychosocial or behavioral therapy services.
•
Concurrent therapy is associated with indicators of clinical need and complexity.
•
Service patterns suggest therapy intensity aligns with behavioral health problems.
•
Therapy services may improve buprenorphine retention, but future research is needed.
Abstract
Introduction
Current evidence indicates that buprenorphine is a highly effective treatment for opioid use disorder (OUD), though premature medication discontinuation is common. Research on concurrent psychosocial and behavioral therapy services and related outcomes is limited.
The goal of this study was to define patterns of OUD-related psychosocial and behavioral therapy services received in the first 6 months after buprenorphine initiation, identify patients' characteristics associated with service patterns, and examine the course of buprenorphine treatment, including the association of therapy with medication treatment duration.
Methods
We analyzed 2013–2018 MarketScan Multi-State Medicaid claims data. The sample included adults aged 18–64 years at buprenorphine initiation with treatment episodes of at least 7 days (n = 61,976). We used group-based trajectory models to define therapy service patterns and multinomial logistic regression to identify pre-treatment patient characteristics associated with therapy trajectories. Multinomial propensity-score weighted Cox proportional hazards regression estimated time to buprenorphine discontinuation and unweighted Cox proportional hazards models estimated risk of adverse health care events during buprenorphine treatment (all-cause and opioid-related inpatient and emergency department services, overdose treatment).
Results
We identified three trajectories of psychosocial and behavioral therapy services: none (73.8%), low-intensity (17.2%), and high-intensity (9.0%). Compared to those without therapy, low-intensity and high-intensity service patterns were associated with behavioral health diagnoses and medical treatment for opioid overdose in the baseline period prior to buprenorphine initiation. The hazard of buprenorphine discontinuation was significantly lower for low-intensity (HR = 0.55; 95% CI, 0.54–0.57) and high-intensity (HR = 0.71; 95% CI, 0.67–0.74) therapy groups compared to those without therapy services. Yet patients in the high-intensity therapy group had increased risk of opioid-related health care events during buprenorphine treatment, including medical treatment for opioid overdose (HR = 1.29; 95% CI, 1.01–1.64).
Conclusion
Most patients received little or no OUD-related psychosocial and behavioral therapy after initiating buprenorphine treatment. Patients who received therapy had characteristics indicating greater treatment needs as well as more complex treatment courses. Concurrent therapy services may help to address premature buprenorphine discontinuation, particularly for patients with high-risk clinical profiles; however, future prospective research should determine whether therapy is effective for extending buprenorphine retention.
). The strength of the evidence base is reflected in the steady shift away from the term “medication-assisted treatment” (MAT) toward “medications for opioid use disorder” (MOUD). The robust clinical trial data demonstrating the effectiveness of buprenorphine and other medications are also reflected in evidence-based recommendations that emphasize MOUD as central and psychosocial services as complementary in OUD treatment (
). Due to extensive evidence on the effectiveness of pharmacotherapy and high rates of medication discontinuation, retention on MOUD is often selected as a treatment outcome in clinical trials (
Trial data for psychosocial treatment delivered alongside MOUD have some limitations, but largely show no added benefit in treatment retention or outcomes beyond MOUD itself (
). Although select studies of psychosocial interventions focusing specifically on treatment engagement and continuity have demonstrated a positive relationship with buprenorphine retention (
A randomized pilot clinical trial to evaluate the efficacy of community reinforcement and family training for treatment retention (CRAFT-T) for improving outcomes for patients completing opioid detoxification.
), it is unclear if this relationship extends to observational settings and populations that allow for broader outcome measures over a longer timeframe and that may differ from clinical trials in terms of care delivery (i.e., usual treatment without the support services often provided in trials) and patient composition (e.g., more complex medical and behavioral health comorbidities) (
The limitations of the literature may explain, in part, the discrepancy between the overall evidence base and current treatment practices and policies. Expert opinion supports psychosocial and behavioral therapy as particularly important for OUD treatment engagement and retention (
). The perceived importance of therapy and counseling is also represented in third-party payer policies stipulating psychosocial treatment or referral as a prerequisite for MOUD, often including state Medicaid programs (
Substance Abuse and Mental Health Services Administration Medicaid coverage of medication-assisted treatment for alcohol and opioid use disorders and of medication for the reversal of opioid overdose (SMA-18-5093). Substance Abuse and Mental Health Services Administration.
Medicaid and CHIP Payment and Access Commission (MACPAC) Medicaid and the opioid epidemic. Report to congress on Medicaid and CHIP; Medicaid and CHIP Payment and Access Commission (MACPAC).
). Importantly, Medicaid is also a major financer of behavioral health care in the United States, including public services. In addition to covering a disproportionate share of adults with OUD (
Medicaid and CHIP Payment and Access Commission (MACPAC) Medicaid and the opioid epidemic. Report to congress on Medicaid and CHIP; Medicaid and CHIP Payment and Access Commission (MACPAC).
), Medicaid is the primary source of insurance for low-income Americans, who have elevated rates of psychiatric disorders compared to those with higher incomes (
). Identifying treatment patterns in this group is critical to increase understanding of current practices and facilitate improvements in patient outcomes.
One analysis of Medicaid administrative records of adolescents found lower risk of attrition for buprenorphine (with or without psychosocial services) compared to psychosocial treatment alone (
), but this study did not compare combined treatment to buprenorphine alone. A second study of Missouri Medicaid data examined buprenorphine received across settings defined as either “medical” or “psychosocial”, with results suggesting that treatment in medical settings was associated with lower risk of both buprenorphine discontinuation and substance-related hospital visits (
). These two studies are important additions to the literature on variation in patient and treatment characteristics across modalities and settings, but further research could fill gaps on treatment predictors, course, and outcomes for patients receiving concurrent buprenorphine and psychosocial services.
To improve understanding of the role of psychosocial and behavioral therapy in OUD treatment, this study examined services received in the first 6 months after initiating buprenorphine. We focused on this particularly high-risk period for medication discontinuation (
), as it represents a critical window to identify factors associated with treatment retention and could facilitate improvements along the cascade of care (
). The main objectives were to: 1) define longitudinal patterns of therapy services and identify baseline patient characteristics associated with these patterns, 2) describe the course of treatment across patient subgroups, and 3) examine the association between therapy service patterns and buprenorphine treatment retention.
2. Material and methods
2.1 Study design and data
We conducted a retrospective cohort study using the MarketScan Multi-State Medicaid database, with data collected longitudinally from 2013 to 2018 on millions of U.S. enrollees, including enrollment information, demographic characteristics, diagnoses, and comprehensive health services and prescription data.
2.2 Sample
The study sample comprised adults aged 18–64 years at buprenorphine initiation, defined as an index claim following at least a 180-day baseline period with no observed buprenorphine prescriptions. We excluded buprenorphine episodes <7 days in duration, which could represent detoxification rather than maintenance treatment. For patients with multiple treatment episodes, we included the first qualifying episode in analyses. To capture new therapeutic episodes in addition to new medication treatment episodes, we also excluded individuals with any OUD-related psychosocial or behavioral therapy services during the 180- day baseline period prior to initiating buprenorphine. Finally, we excluded older individuals and those otherwise eligible for Medicare as well as individuals with behavioral health and prescription drug coverage not captured in the database (Appendix Fig. A1).
2.3 Measures
2.3.1 OUD treatment outcomes
We identified buprenorphine treatment using National Drug Codes (NDC) from prescription drug claims. The start date for new treatment episodes was the date on the index buprenorphine claim. We measured treatment duration using the number of days' supply listed on buprenorphine claims. Similar to prior studies of buprenorphine treatment episodes (
), we defined the discontinuation outcome as >30 days without buprenorphine supply, using the last date with supply as the end of the treatment episode.
We identified psychosocial and behavioral therapy services using Healthcare Common Procedure Coding System (HCPCS) procedure codes from outpatient claims (
). The study team selected services based on the Healthcare Effectiveness Data and Information Set (HEDIS) measures of initiation and engagement in substance use treatment, excluding services that were outside the scope of psychosocial and behavioral therapy (Appendix Table A1) (
). For example, services included individual, group, and family psychotherapy or counseling services and are referred to hereafter simply as “therapy”. Trajectory analyses were limited to therapy services with an OUD diagnosis to capture treatment episodes specifically focused on opioid addiction (
). However, we did not require OUD to be the primary diagnosis with the aim of capturing comprehensive OUD treatment addressing common psychiatric comorbidities (
), including substance use disorders for which medication treatment is not available. In addition, we conducted a sensitivity analysis, including all services with a mental health or substance use disorder (MHSUD) diagnosis (i.e., not limited to OUD-related therapy). We included patients without documented OUD diagnoses in the baseline period prior to buprenorphine initiation because a substantial proportion of patients receive buprenorphine for opioid maintenance without an OUD diagnosis (
We measured pre-treatment baseline sociodemographic and clinical characteristics in the 180-day period prior to buprenorphine initiation. We identified sociodemographic characteristics using enrollment data on sex (male, female), age (years), race/ethnicity (non-Hispanic white, non-Hispanic Black, Hispanic, other), and insurance plan type (capitation or fee-for-service). Clinical characteristics included baseline medical and behavioral health comorbidities and health service use. We identified comorbidities using primary and secondary diagnosis codes from all service claims.
Mental health indicators included diagnoses of depression, anxiety, bipolar disorder, schizophrenia, and any other mental health disorder. Substance use disorders included indicators for opioids, alcohol, cannabis, tobacco, cocaine, sedatives, amphetamines, and any other drug (Appendix Table A2). Medical indicators included hepatitis C and a weighted composite measure of overall medical burden. We calculated a modified Elixhauser Comorbidity Index score that excluded diagnoses of hepatitis C and behavioral health disorders from the final index, which accounted for 25 medical comorbidities (e.g. HIV, hypertension, obesity, diabetes, cancer) that are associated with hospital-based resource use and mortality (
Baseline health service use included indicators for all-cause inpatient and emergency department (ED) services, opioid overdose, non–opioid drug overdose, and any prescription opioids or psychotropic drugs (i.e., antidepressants, antipsychotics, stimulants, mood stabilizers, benzodiazepines, sedatives).
2.3.3 Treatment characteristics
Because previous research has shown an association with treatment discontinuation (
), we measured the initial buprenorphine dosage, defined as the dosage on the date of the index claim. We measured the proportion of days covered (PDC) by buprenorphine during the treatment period because PDC is a quality measure of the Centers for Medicare and Medicaid Services for other chronic conditions with outcomes that depend on treatment continuity (
Relationship between buprenorphine adherence and relapse, health care utilization and costs in privately and publicly insured patients with opioid use disorder.
To better understand the course of treatment, we also measured adverse health care events observed during the first 6 months of buprenorphine treatment. In addition to all-cause inpatient and ED services, opioid and non–opioid drug overdose, we measured opioid-related events treated in inpatient or ED settings, defined according to methods from the Agency for Healthcare Research and Quality (
). Opioid-related events include inpatient or ED services with diagnosis codes for opioid use disorder, poisoning, and adverse effects.
2.4 Analyses
2.4.1 Therapy patterns and treatment course
We used group-based trajectory modeling (GBTM) to identify patterns of OUD-related therapy services received in the first 6 months after initiating buprenorphine. GBTM is a flexible approach that captures variation in the population to identify longitudinal latent classes with distinctive trajectory paths (
). We estimated trajectories based on the number of days with therapy services in each 30-day period, using a zero-inflated Poisson distribution to account for the large number of individuals with no observed therapy services. The study team based the final model selection on established statistical criteria (Appendix Table A3) (
). Within the sample, separate analyses estimated patterns of MHSUD-related therapy services to identify groups for sensitivity analyses.
We calculated unadjusted sample statistics to show baseline sociodemographic and clinical characteristics across trajectory groups. Adjusted multinomial logistic regression estimated the associations of pre-treatment patient characteristics with trajectory group membership, reported as relative risk ratios (RR).
To examine adverse health care events observed during buprenorphine treatment, we calculated unadjusted rates to determine the number of individuals with events per 1000 person-months of buprenorphine treatment for each trajectory group. Adjusted Cox proportional hazards regressions separately estimated the time to each event. In addition to adjusting for all pre-treatment patient characteristics, models adjusted for the initial buprenorphine dose and the proportion of days covered by buprenorphine from initiation to the time of the event (or the end of the episode for patients with no event).
2.4.2 Buprenorphine treatment duration
To examine the association of therapy patterns with buprenorphine treatment duration, adjusted Cox proportional hazards models estimated time to discontinuation across trajectory groups. We used propensity score weighting to account for differences across groups in observed pre-treatment covariates, including baseline sociodemographic and clinical characteristics, with the goal of improving comparability across groups and reducing bias in the effect estimates of buprenorphine discontinuation across trajectory groups (
). Using established methods for multiple treatments described in detail elsewhere, the study calculated multinomial propensity score weights using generalized boosted modeling (
). This method iteratively fits regression trees to capture nonlinearities in the relationship of pre-treatment covariates to treatment group membership and minimizes imbalance across therapy groups. To assess balance across groups, we calculated the absolute standardized difference in means for each pre-treatment covariate before and after propensity score weighting (Appendix Fig. A2). We used average treatment effect (ATE) weights to generate estimates representing the difference in buprenorphine treatment duration if all patients followed a given therapy trajectory compared to all patients receiving no therapy. We then examined the association between therapy patterns and buprenorphine discontinuation using propensity score weighted Cox models, fully adjusting for baseline characteristics to obtain doubly-robust estimates (
) with additional adjustments for the initial buprenorphine dose and PDC.
We conducted two sensitivity analyses of the association between therapy patterns and buprenorphine treatment duration. In the first, we assessed the robustness of Cox estimates of the association between therapy patterns and buprenorphine treatment duration to unmeasured confounding (Appendix Tables A4 and A5) (
). In the second, we excluded comparisons across therapy trajectory groups, instead including a time-varying indicator for OUD therapy to identify the overall association of OUD therapy with buprenorphine treatment duration in the full sample (Appendix Table A6 and Fig. A3).
3. Results
3.1 Patterns of psychosocial and behavioral therapy services
In trajectory analyses, a three-group specification for therapy service use patterns performed well on all tests of model adequacy (Appendix Table A3). Fig. 1 shows therapy trajectories in the first 6 months after initiating buprenorphine. Most patients were categorized in the group with little to no services (n = 45,720; 73.8%), averaging <1 day with therapy services per month (mean = 0.07). Approximately one-fifth followed a low-intensity trajectory (n = 10,689; 17.2%), averaging about 2 service days per month (mean = 2.06). Nearly one-tenth followed a high-intensity trajectory (n = 5567; 9.0%), averaging about 7 therapy service days per month (mean = 6.94), with higher use in the first month (mean = 8.26) declining over the study period to about 5 services in the sixth and final month (mean = 5.28). Most OUD services were for individual, group, or family therapy/counseling (65.6%) followed by therapy/counseling as part of a behavioral health program (30.3%) (Appendix Table A7).
Fig. 1Patterns of psychosocial and behavioral therapy in the first 180 days after buprenorphine initiation.
Results from sensitivity analyses of all MHSUD-related therapy services were similar to the main findings (Appendix Fig. B1, Appendix Tables B1–B2).
3.2 Patient characteristics associated with receipt of therapy services
The left panel of Table 1 shows unadjusted baseline characteristics across trajectory groups with no therapy, low-intensity, and high-intensity therapy patterns. Across all groups, the sample initiating buprenorphine treatment was mostly female, with average age of 35 years, and overwhelmingly white (90%). However, compared to the group with no therapy, the low- and high-intensity groups had a significantly smaller proportion of females and a higher proportion of whites (Table 1).
Table 1Unadjusted and adjusted baseline characteristics associated with patterns of psychosocial and behavioral therapy services.
Notes: Therapy groups were derived from group-based trajectory models. For unadjusted statistics, bold font indicates significance at p < 0.05 compared to the no therapy group, and * indicates significance at p < 0.05 compared to the low-intensity therapy group. Adjusted estimates were derived from multinomial logistic regression models, fully adjusted for pre-treatment patient characteristics.
The right panel of Table 1 shows fully adjusted multinomial logistic regression estimates of the association between baseline characteristics and trajectory group membership.
Compared to the group with no therapy, low- and high-intensity therapy patterns were associated with capitated Medicaid plans (low-intensity RR = 1.96; 95% CI, 1.87–2.05; high-intensity RR = 2.20; 95% CI, 2.06–2.34), documented OUD (low-intensity RR = 2.36; 95% CI, 2.25–2.47; high-intensity RR = 3.69; 95%CI, 3.44–3.95), cannabis use disorder (low-intensity RR = 1.33; 95% CI, 1.21–1.47; high-intensity RR = 1.69; 95% CI, 1.51–1.89), all-cause ED services (low-intensity RR = 1.13; 95% CI, 1.07–1.19; high-intensity RR = 1.21; 95% CI, 1.13–1.30), and opioid overdose history (low-intensity RR = 1.22; 95% CI, 1.08–1.38; high-intensity RR = 1.72; 95% CI, 1.50–1.96) at baseline prior to treatment initiation. Patients who were less likely to receive therapy included those with “other” race/ethnicity (low-intensity RR = 0.76; 95% CI, 0.67–0.86; high-intensity RR = 0.64; 95% CI, 0.54–0.77), tobacco use disorder (low-intensity RR = 0.91; 95% CI, 0.86–0.96; high-intensity RR = 0.90; 95% CI, 0.84–0.96), and benzodiazepine (low-intensity RR = 0.72; 95% CI, 0.67–0.77; high-intensity RR = 0.76; 95% CI, 0.69–0.83) or sedative (e.g., zolpidem) (low-intensity RR = 0.76; 95% CI, 0.67–0.87; high-intensity RR = 0.71; 95% CI, 0.60–0.85) prescriptions at baseline (Table 1).
Results of sensitivity analyses across MHSUD-related therapy groups were consistent with the main findings (Appendix Table B3), with additional baseline behavioral health diagnoses significantly associated with receipt of therapy in adjusted models (e.g., depression, bipolar disorder, schizophrenia, other mental illness, alcohol use disorder, and amphetamine use disorder).
3.3 Patients' characteristics associated with intensity of therapy services
In unadjusted analyses comparing baseline characteristics of the high-intensity therapy group to the low-intensity group, no significant differences existed in terms of sex or race/ethnicity, but the high-intensity therapy group was slightly younger on average and a greater proportion had capitated Medicaid plans (Table 1). Most clinical characteristics exhibited a gradient, whereby the high-intensity group had significantly higher rates of behavioral health diagnoses, medical comorbidities, all-cause inpatient and ED service use, medically treated overdoses, and prescriptions for antidepressants or antipsychotics.
The far-right panel of Table 1 shows fully adjusted multinomial logistic regression estimates of the association between baseline characteristics and intensity of therapy services. Compared to the low-intensity group, high-intensity therapy was associated with capitated Medicaid plans (RR = 1.12; 95% CI, 1.04–1.21); depression (RR = 1.11; 95% CI, 1.01–1.21); documented OUD (RR = 1.56; 95% CI, 1.45–1.69); and substance use disorder diagnoses for alcohol (RR = 1.55; 95% CI, 1.37–1.75), cannabis (RR = 1.27; 95% CI, 1.12–1.44), cocaine (RR = 1.39; 95% CI, 1.19–1.61), and sedatives (RR = 1.26; 95% CI, 1.03–1.54); hepatitis C (RR = 1.18; 95% CI, 1.05–1.32); all-cause inpatient service use (RR = 1.31; 95% CI, 1.20–1.44); medically treated opioid overdose (RR = 1.41; 95% CI, 1.21–1.64); and prescription antipsychotics (RR = 1.15; 95% CI, 1.03–1.30).
Results from sensitivity analyses comparing MHSUD-related therapy groups were consistent with the main findings (Appendix Table B3), though some adjusted associations with high-intensity therapy were no longer significant (e.g., depression, sedative use disorder, hepatitis C, prescription antipsychotics), and additional associations emerged (e.g., amphetamine use disorder, other drug use disorders).
3.4 Association of therapy services with buprenorphine discontinuation
The left panel of Table 1 shows unadjusted characteristics of buprenorphine treatment duration. The median length of buprenorphine treatment episodes was shortest for the group with no therapy (84 days) but comparable for the low-intensity (145 days) and high-intensity (136 days) groups. Similarly, a lower proportion of patients in the group with no therapy (31.9%) reached the 180-day benchmark for minimum duration of pharmacotherapy than in the low-intensity (43.3%) and high-intensity (42.4%) groups. Overall, about one-third of the sample had 180-day or longer buprenorphine episodes (n = 21,578; 34.8%) and two-thirds of those with minimum duration treatment were in the group with no therapy (n = 14,591; 67.6%).
Propensity score weighting achieved good balance on all covariates across therapy groups (Appendix Fig. A2). Table 2 shows fully adjusted results from Cox proportional hazard models estimating time to buprenorphine discontinuation in the first 180 days after treatment initiation. In unweighted models, both low-intensity (HR = 0.55; 95% CI, 0.53–0.57) and high-intensity (HR = 0.61; 95% CI, 0.58–0.63) therapy were associated with reduced risk of discontinuation. In propensity score weighted models, these associations remained statistically significant, though the hazard of buprenorphine discontinuation for high-intensity therapy was attenuated (HR = 0.71; 95% CI, 0.67–0.74).
Table 2Adjusted hazard of buprenorphine discontinuation in the first 6 months of treatment.
Note: Estimates were derived from fully adjusted Cox proportional hazards regression models estimating time to buprenorphine discontinuation; weighted estimates are based on multinomial propensity scores and represent the average treatment effect if all patients in the sample had followed a given therapy trajectory compared to all patients receiving no therapy.
In sensitivity analyses examining time-varying receipt of OUD-related psychosocial and behavioral services in the full sample (i.e., no comparisons across therapy trajectory groups), therapy was also associated with significantly reduced risk of buprenorphine discontinuation (HR = 0.40; 95% CI, 0.37–0.43; Appendix Table A6; Appendix Fig. A3).
Results from sensitivity analyses of MHSUD-related therapy groups and services were similar to the main findings (Appendix Figs. B2–B3, Appendix Tables B3–B6).
3.5 Additional characteristics of the treatment course across therapy service patterns
In unadjusted analyses (Table 1), a gradient emerged in mean initial buprenorphine dosage (none = 14.5 mg, low-intensity = 13.3 mg; high-intensity = 12.8 mg), though differences were minimal and may not be clinically meaningful. Similarly, the proportion of patients with initial dosage ≥16 mg was highest for the group with no therapy (61.8%), lower for the low-intensity therapy group (52.0%), and lowest for the high-intensity group (48.4%). On average, all groups had a high proportion of days covered (PDC) by buprenorphine throughout the treatment episode, including in the first 180 treatment days (Table 1). Though both the low- and high-intensity groups had significantly lower mean PDC compared to the group with no therapy, these differences were also minimal.
Table 3 shows unadjusted rates and adjusted estimates of adverse health care events observed during the first 6 months after initiating buprenorphine. Both the low- and high-intensity therapy groups had significantly higher unadjusted rates of all-cause and opioid-related ED service use as well as medical treatment for opioid overdose. The high-intensity therapy group also had higher rates of all-cause and opioid-related inpatient service use as well as medical treatment for opioid overdose. In adjusted analyses, only the increased risk of all-cause ED services remained significant for both the low-intensity (HR = 1.05; 95% CI, 1.01–1.09) and high-intensity (HR = 1.16; 95% CI, 1.10–1.22) groups. For the high-intensity group, the increased risk of opioid-related ED services (HR = 1.27; 95% CI, 1.10–1.46) and medical treatment for opioid overdose (HR = 1.29; 95% CI, 1.01–1.64) remained significant. Additionally, the low-intensity group had significantly lower risk of all-cause (HR = 0.92; 95% CI, 0.85–0.99) and opioid-related (HR = 0.82; 95% CI, 0.71–0.94) inpatient services in adjusted analyses.
Table 3Unadjusted and adjusted adverse health care events observed during buprenorphine treatment.
Note: For unadjusted rates, bold font indicates significance at p < 0.05 compared to the no therapy group, and * indicates significance at p < 0.05 compared to the low-intensity therapy group. Adjusted estimates were derived from Cox proportional hazards regression models for each outcome, adjusted for all pre-treatment patient characteristics, initial buprenorphine dosage, and proportion of days covered by buprenorphine from initiation until the event of interest or buprenorphine discontinuation if no event was observed.
In adjusted analyses compared to the low-intensity therapy group, rates of all adverse health care events observed during the first 6 months after buprenorphine initiation were significantly higher for the high-intensity group (Table 3). However, in adjusted analyses, only the increased risk of opioid-related inpatient services (HR = 1.40; 95% CI, 1.18–1.67) and all-cause ED services (HR = 1.11; 95% CI, 1.05–1.17) remained significant (Table 3).
Results of sensitivity analyses across MHSUD-related therapy groups were consistent with the main findings (Appendix Tables B3, B7), with fewer significant associations in adjusted analyses. Compared to the group with no therapy, low-intensity therapy was associated only with higher risk of all-cause ED services, while high-intensity therapy was associated only with higher risk of opioid-related inpatient and ED services. Compared to the low-intensity group, high-intensity therapy was associated only with higher risk of opioid-related inpatient services.
4. Discussion
In this study, we examined the association of concurrent OUD-related psychosocial and behavioral therapy services with buprenorphine treatment retention during a high-risk period for attrition. We found that patients initiating and remaining on buprenorphine for at least 7 days generally followed one of three therapy trajectories in the first six months after initiating medication. The majority of patients had little to no therapy (74%), while roughly 17% and 9% had low-intensity and high-intensity therapy trajectories, respectively. Overall rates of buprenorphine receipt with and without therapy are consistent with prior research of privately insured patients in which about 69% of patients received MOUD alone and 31% received MOUD with psychosocial treatment in a given year from 2013 to 2016 (
). In a study of Missouri Medicaid patients with OUD, about 25% of patients received buprenorphine in a psychosocial setting and 75% in a general clinical setting, where psychosocial services were less prevalent but not completely absent (
In analyses of buprenorphine retention among patients with at least 7 days of treatment, receiving OUD-related psychosocial and behavioral therapy services was associated with lower risk of medication discontinuation in the first 180 days after treatment initiation, raising the possibility that concurrent therapy may help to reduce high attrition rates commonly observed early in the course of treatment (
). Notably, the association of therapy services with reduced risk of buprenorphine discontinuation was consistent across estimates that weighted patients with different service patterns to have similar pre-treatment characteristics and estimates that assessed therapy continuously over time regardless of overall therapy trajectories, strengthening the implications for concurrent psychosocial and behavioral therapy services as one potential approach for addressing premature buprenorphine discontinuation. These findings are consistent with expert opinion that psychosocial services are important for treatment engagement and continuity (
) and with some trial data showing a positive relationship between buprenorphine treatment duration and psychosocial services that focus explicitly on engagement and retention (
A randomized pilot clinical trial to evaluate the efficacy of community reinforcement and family training for treatment retention (CRAFT-T) for improving outcomes for patients completing opioid detoxification.
). However, the observational design of this study limits causal inferences, and future prospective research should determine whether therapy is effective for extending buprenorphine retention.
While fewer individuals in the group with minimal or no therapy services reached the benchmark for minimum duration of pharmacotherapy (
), a substantial proportion had 180 days or more of buprenorphine treatment, indicating that many patients persist in treatment with medication alone. Thus, insurance policies requiring referral or receipt of psychosocial services to receive buprenorphine or other MOUD may create overly restrictive barriers to highly effective medication treatment. Furthermore, many patients have limited or no access to such services due to fragmentation of health care, workforce shortages, underfunding, overburdened providers and programs, and individual barriers to health care (e.g., transportation difficulties, employment/family obligations). Ensuring adequate availability of psychosocial support services is important to facilitate access for patients who could benefit from therapy, particularly with recent expansions in telehealth services that could address current barriers to care (
). To improve alignment of treatment needs and services, more research could help to build understanding of the complex relationships among patient characteristics, care components, and outcomes.
We found that receipt of therapy services was associated with baseline behavioral health diagnoses, including opioid and cannabis use disorders and, for those receiving high-intensity therapy, alcohol, cocaine, and sedative use disorders. In addition, receipt of therapy was associated with medical treatment for opioid overdose at baseline. These findings are consistent with prior research showing higher rates of mental health and substance use problems among patients with combined medication and psychosocial treatment (
), indicating that patients who receive therapy in conjunction with buprenorphine may have more complex behavioral health treatment needs. Furthermore, high-intensity service patterns were associated with comorbid addiction to substances for which no FDA-approved medications are currently available (e.g., cocaine), potentially reflecting the limited treatment options available to patients with multiple substance use problems.
Patients receiving high-intensity therapy also had higher risk of opioid-related health care events during buprenorphine treatment episodes, suggesting that these patients also have more complex treatment courses. Considering the elevated risk profile in this group, these findings could indicate that patients with more complex behavioral health problems require more intensive or longer duration treatment to stabilize. In addition, the increased risk of opioid overdose and opioid-related ED care during buprenorphine treatment may signal a need for emerging hospital-based interventions to include protocols for care coordination and communication with providers of patients currently receiving OUD treatment in addition to OUD treatment induction and referrals. While these analyses adjusted for a range of pre-treatment (baseline) and buprenorphine treatment (initial dosage, PDC) characteristics, they were not propensity-score weighted. Thus, due to potential residual confounding across therapy groups, we cannot draw causal inferences concerning the effects of therapy on risks of these adverse outcomes.
This study and prior research demonstrating the benefits of psychosocial and behavioral therapy in observational settings with more generalizable patient populations may have implications in an era of increasingly personalized medicine, with a growing emphasis on tailoring services to meet individual needs of patients. Overall, patients with greater or more complex treatment needs received more intensive services, which may contribute to improved retention in these groups during a period with high risk of attrition. Future research should continue to examine variation across patient subgroups and therapy services in the relationship of therapy to buprenorphine retention. For example, patients who have comorbid substance use disorders that lack effective medication options could represent a priority group for psychosocial services. This work is critical to align patient needs with evidence-based treatment to improve retention and patient outcomes.
4.1 Limitations
We could not ascertain specific types of therapy because the data do not include codes indicating, for example, cognitive behavioral therapy or contingency management. Some prior studies have compared the effectiveness of various therapeutic approaches (
), but more work is needed in this area. Furthermore, the data do not capture services that are not reimbursed by Medicaid and therefore do not generate a Medicaid claim, such as self-help or peer-support services, brief or limited services delivered as part of medical management, and services funded by block grants or other non-Medicaid sources. However, underestimating service use would have resulted in more conservative effect estimates. In addition, we could not distinguish buprenorphine treatment characteristics across settings because prescription drug data do not include inpatient medication dispensing or the specific outpatient setting (e.g., primary care vs. opioid treatment program). Although we adjusted for a wide range of individual characteristics and propensity score methods achieved good balance across groups, residual confounding could result from unobserved factors. To further address this concern, sensitivity analyses to determine the robustness of results (
) indicated that unmeasured confounding was unlikely to explain away the associations between therapy patterns and buprenorphine retention. Additionally, results may not generalize to patients who discontinue buprenorphine in the first week of treatment, as the sample was restricted to episodes at least 7 days in duration. Finally, findings from Medicaid data may not be generalizable to individuals with no insurance or private coverage and, although the study drew the data from a large multi-state sample, the findings may not be representative of all Medicaid enrollees.
5. Conclusion
This is the first known study characterizing longitudinal trends in psychosocial and behavioral therapy received alongside buprenorphine in the treatment of adults with opioid use disorder. Most patients received little or no therapy services. Among those with concurrent therapy, the intensity of services corresponded to indicators of treatment need, potentially signifying appropriate alignment between patient characteristics and clinical practices. In addition, receipt of psychosocial and behavioral therapy was associated with reduced risk of early buprenorphine discontinuation, suggesting a possible protective role of therapy in improving treatment retention and potentially other outcomes for patients.
Acknowledgements
Financial support for this work was provided by grants from the National Institute on Drug Abuse (NIDA) [grant numbers K01DA049950 and K23DA044342] with additional support from the Foundation for Opioid Response Efforts (FORE) and the National Center for Advancing Translational Sciences [grant number UL1TR003017].
Supplementary material for sensitivity analyses of therapy services with a mental health or substance use disorder (MHSUD) diagnosis (i.e. not restricted to OUD services)
References
Amato L.
Minozzi S.
Davoli M.
Vecchi S.
Psychosocial combined with agonist maintenance treatments versus agonist maintenance treatments alone for treatment of opioid dependence.
A randomized pilot clinical trial to evaluate the efficacy of community reinforcement and family training for treatment retention (CRAFT-T) for improving outcomes for patients completing opioid detoxification.
Relationship between buprenorphine adherence and relapse, health care utilization and costs in privately and publicly insured patients with opioid use disorder.
Substance Abuse and Mental Health Services Administration
Medicaid coverage of medication-assisted treatment for alcohol and opioid use disorders and of medication for the reversal of opioid overdose (SMA-18-5093). Substance Abuse and Mental Health Services Administration.
Guidelines for the psychosocially assisted pharmacological treatment of opioid dependence.
World Health Organization,
2009 (http://books.google.com/books?hl=en&lr=&id=yKV6i3ZK86kC&oi=fnd&pg=PP2&dq=Guidelines+for+the+Psychosocially+Assisted+Pharmacological+Treatment+of+Opioid+Dependen ce&ots=NtatF6p9lt&sig=batQx8UMQBGw- HWDpEa3XFwrj00%5Cnhttp://books.google.com/books?hl=en&lr=&id)