Suicide is a leading cause of mortality in the United States, with more than 45,000 deaths each year (
1). Over the past decade, the rate of suicide in the U.S. population has risen by more than 20% (
1,
2). Moreover, opioid use disorder (OUD) is an independent risk factor for suicide (
3–
5). Medication treatment has formed the backbone of OUD treatment for several decades (
6). Three medications (methadone, buprenorphine, and naltrexone) are approved by the U.S. Food and Drug Administration (FDA) for the treatment of OUD (
7). Use of medications for opioid use disorder (MOUD) has been shown to increase quality of life, improve health outcomes, and decrease illicit drug use (
8).
Population-based studies have consistently demonstrated that MOUD is associated with reduced all-cause mortality (
9). Two studies have also examined the potential for MOUD to decrease risk of suicide mortality. In an Australian cohort study of over 40,000 individuals with OUD who initiated MOUD with methadone or buprenorphine, Degenhardt et al. (
10) found that time on MOUD was associated with reduced risk of suicide mortality as compared with time off MOUD. Conversely, risk was elevated during MOUD transition periods, including the first 14 days on treatment and the first 14 days off treatment. Similarly, in a Swedish cohort study of over 20,000 individuals with OUD or alcohol use disorder, Molero et al. (
11) observed that methadone, buprenorphine, and naltrexone treatment periods were associated with reduced risk of suicidal behaviors (defined as suicide attempts or deaths by suicide) compared with nontreatment periods, but the finding was significant only for methadone. However, neither study performed adjustment for individual patient characteristics or accounted for contact with the health care system when comparing the various forms of MOUD. A small cohort study of just under 6,000 patients with OUD by Riblet et al. (
12) indicates that some of the mortality benefits for patients undergoing MOUD may be related to a high level of contact with the health care system.
We designed a cohort of Department of Veterans Affairs (VA) patients with OUD who received MOUD between 2003 and 2017. Our goals were to determine whether MOUD is associated with a reduced risk of suicide mortality when controlling for patient characteristics and level of contact with the health care system, as well as to determine whether the findings differed among specific MOUD agents.
Methods
We received approval for this study from the Veterans Institutional Review Board of Northern New England.
Study Population
We identified all VA patients who were diagnosed with OUD at two or more visits on different days within the same year between 2003 and 2017 and received one or more days of MOUD (defined below). Patients entered the study cohort on the first day of MOUD and were followed for up to 5 years, or until the date of their death, or until the end of the predefined study period (December 31, 2017), whichever came first.
Data Sources
We collected demographic and diagnostic information as well as health care utilization and pharmacy data from the VA Corporate Data Warehouse. Because some patients receive portions of their care in the community, we also used linked VA “fee basis” (FB) files and Centers for Medicare and Medicaid Services (CMS) claims files. We identified causes of death by using the VA–Department of Defense Mortality Data Repository. The Mortality Data Repository links mortality data on veterans from the Centers for Disease Control and Prevention (CDC) National Death Index with the VA Corporate Warehouse (
https://www.mirecc.va.gov/suicideprevention/Data/data_index.asp).
Mortality Outcomes
We categorized the causes of death using the CDC National Center for Health Statistics ICD-10 definitions. Because of the high risk of bias in attributing cause of death, we considered three mortality outcomes: suicide mortality, external-cause mortality, and all-cause mortality (
13). External-cause mortality per the CDC definition includes death due to any accident or injury, either intentional or accidental. Common examples include unintentional injuries (e.g., motor vehicle accidents, drowning, and fires), suicide deaths, homicide deaths, and overdoses deaths regardless of intent (self-harm, accident, or undetermined). The outcomes are not mutually exclusive; for example, suicide is nested within external causes.
Covariates
There were three main groups of covariates: variables related to MOUD treatment, patient demographic characteristics, and health care contacts. We derived patient covariates from the VA Corporate Data Warehouse as well as linked CMS and FB files.
MOUD Treatment
We considered MOUD exposure to include the receipt of methadone, buprenorphine, or naltrexone. Methadone for OUD is typically dispensed in liquid formulation from federally licensed opioid treatment programs. Because these doses are not usually recorded as prescriptions, we used service delivery codes as a proxy for dosing. To assess for buprenorphine exposure, we included prescriptions of sublingual buprenorphine and of sublingual combination buprenorphine-naloxone formulations. Because the transdermal, buccal, and intravenous formulations of buprenorphine were only FDA approved for the treatment of pain rather than OUD during our study period (
14), we excluded these formulations from our analysis. To assess for naltrexone exposure, we included oral as well as long-acting injectable formulations. We identified buprenorphine and naltrexone using a combination of prescription fill data (VA Corporate Data Warehouse) and pharmacy claim National Drug Codes (CMS and FB).
We developed several rules to create MOUD coverage periods. First, we assigned patients to a new MOUD coverage period if they received 1 or more days of MOUD after a period of 7 or more days without MOUD. Second, we continued MOUD coverage periods unless patients had more than 7 consecutive days without MOUD. Third, we adapted an approach developed by Pearce et al. (
15) to account for MOUD coverage during inpatient or residential stays. Finally, because previous research suggests that MOUD starting periods (first 14 days on) and stopping periods (first 14 days off) are associated with a higher risk of mortality (
10), we divided time into four MOUD treatment periods: starting, stopping, stable on, and stable off. The stable off period was the reference group for all analyses.
Patient Demographic Characteristics
Demographic characteristics included age, gender, and race/ethnicity. Race/ethnicity was coded as Black Non-Hispanic, White Non-Hispanic, and Hispanic; those with missing race/ethnicity data were assigned to the “unknown or other” race/ethnicity group. To describe physical and mental health comorbidities, we used previously developed categories (
16) representing separate counts of individual physical diagnoses in the Elixhauser comorbidity index and diagnostic categories from DSM-5, excluding the OUD category, as it is a defining characteristic of the cohort.
Health Care Contacts
We identified each day a participant had any health care visits. We provide increased granularity by identifying non-MOUD substance use disorder visits, mental health visits, and all other visits. We also identified any inpatient mental health services use.
Primary Analysis
Patients were entered into a series of proportional hazard models at the time of MOUD initiation and were followed for 5 years or through the end of 2017, whichever came first. We entered data using the counting process format such that each row of data represented a fixed state. New observations represent a new state (e.g., transitioning from on MOUD to a stopping period). We calculated hazard ratios and 95% confidence intervals to compare risk for each of our three outcomes (suicide mortality, external-cause mortality, and all-cause mortality) for MOUD treatment periods (starting, stopping, stable on), using stable off MOUD as the reference. Our initial models were crude. We then adjusted for age, gender, and race. We then added previous-year mental health and physical health diagnoses into our model. We separately categorized previous-year mental and physical health conditions as none, one or two, and three or more. In this way, we included separate variables for number of prior mental health and physical health conditions. Finally, we added health care utilization for the 90 days prior to the current period. We note that because suicide is a rare outcome, the models needed to be parsimonious to avoid overfitting. To achieve this, we categorized race as White or other and age as >50 years or younger. We categorized days of health care contact into <5 days, 5–20 days, and >20 days. We developed a categorical variable for any inpatient stay. The more specific measures of health care utilization (e.g., mental health services contact days) were correlated with total number of contact days and were not included in the models. Thus,
Table 1 provides finer granularity for all covariates and more details for health care utilization than were included in our modeling. Because individual MOUD agents were introduced at different times and suicide rates among veterans have generally been increasing over time, all models included an additional term for three time periods: 2003–2007, 2008–2012, and 2013–2017. Gender and race were fixed over time. Age and year were defined at the start of each MOUD period (row of data). Comorbidities were calculated annually, and we attributed the previous-year values to each MOUD period to ensure that they were preexisting conditions. Health care contacts were based on the 90 days before the current MOUD period.
Sensitivity Analyses
We performed additional analyses to evaluate the robustness of our primary analysis. The first reduced the maximum follow-up time to 3 years instead of 5 years because patients’ use of MOUD tends to decrease over time, and some patients might have improved to the point where they no longer required OUD treatment. The second converted the data into monthly observations where MOUD was classified as either fully on, fully off, or partial. The outcome and covariates were the same as those used in the primary analysis, but here we used logistic regression. The reason for this data structure and statistical model was that it allowed us to calculate E-values to examine the sensitivity of our findings to observational bias. The E-value represents the minimum strength of association that an unmeasured confounder would need to have with both MOUD exposure and mortality to fully explain away effect estimates (
17).
Examining Individual MOUD Agents
Because the three MOUD agents included in the analysis effectively represent very different treatment programs with varying amounts of required health care contact, we performed an additional analysis that evaluated the effect of each compared with periods off MOUD. We employed the identical methods and model from our primary analysis and added the same set of covariates in a stepwise manner. The only model specification change was to code MOUD status as either on or off, as there would not be adequate statistical power to evaluate transition periods. Patients were assigned annually to the MOUD compound used most, as patients rarely received treatment with multiple compounds in a year.
To examine the effect of each individual MOUD agent on mortality outcomes independently, we performed a similar series of proportional hazard models. For this set of analyses, we only compared the MOUD treatment period with the reference period off MOUD. This approach was necessary because each MOUD sample was less than half the total MOUD sample, so we did not have statistical power to examine transition period independently. We then repeated the analysis as in the primary analysis considering each of the three mortality outcomes (suicide, external, and all-cause) and considering potentially confounding patient covariates in a stepwise fashion (demographic characteristics, diagnosis, and mental health treatment) (
17).
We performed all analyses in SAS, version 9.4 (SAS Institute, Cary, N.C.). We considered baseline covariates to be significantly different when p values were less than 0.05 and hazards to be significantly different when the 95% confidence interval did not cross 1.
Discussion
Our study indicates that MOUD is associated with more than a 50% decrease in risk of suicide mortality during stable treatment periods. Given that there was an even larger reduction in risk of external-cause mortality and all-cause mortality, it is unlikely that our findings are due to misattribution of suicide mortality. We found evidence of a consistent risk pattern associated with stages of treatment, namely, that MOUD starting periods were associated with reduced risk of suicide mortality, but to a lesser degree than stable on MOUD periods. At the same time, MOUD stopping periods were associated with an elevated risk across the three mortality outcomes compared with stable off MOUD periods; however, only the analyses of all-cause mortality had adequate power to detect a significantly elevated risk. We found a differential effect among specific MOUD agents. Buprenorphine was associated with a reduction of suicide mortality risk of more than 65%. In stark contrast, naltrexone treatment showed no indication of a reduced risk of suicide mortality. Methadone’s effects on suicide risk were less clear; in crude models, it was associated with a risk reduction that was similar to buprenorphine, and this finding persisted when patient demographic factors were controlled for. However, when patient mental health diagnoses and contacts were considered in the model, the association of methadone with decreased risk of suicide mortality was diminished and no longer statistically significant. This was largely because those who received methadone had far more mental health contacts.
Our results support the previous preliminary findings that found an association between MOUD and decreased risk of suicide mortality or suicide attempts (
10,
11). The magnitude of the effect, an approximate 50% decrease in risk of suicide mortality, was consistent with the only other study in the literature that specifically evaluated suicide mortality (
10). Our findings regarding MOUD risk periods differed from those in the existing literature. We found that MOUD starting periods were associated with reduction of risk of suicide mortality compared with stable off MOUD periods. Degenhardt et al. found that this period had an elevated risk of suicide mortality (
10). While our findings were similar to those of Degenhardt et al. in finding an increased risk of suicide mortality during MOUD stopping periods, the magnitude of that increase was approximately 50% in our sample and 100% in the Degenhardt et al. study. The explanation for the greater risk for transition periods in the Degenhardt et al. study compared with the present study is unclear, but it is possible that some differences in the MOUD care processes between the Australian New South Wales health care system and the U.S. VA health care system could be explanatory (
18). In line with the literature, we found that naltrexone was not protective against suicide mortality. Yet, unlike Molero et al. (
11), we observed that buprenorphine was associated with a large, significant reduction in suicide mortality risk. This finding aligns with emerging evidence that buprenorphine may have rapid antisuicidal properties (
19). These previous studies suggest that buprenorphine’s unique opioid kappa receptor antagonist properties may be responsible for the effect. We also cannot rule out the possibility that buprenorphine and its associated treatment protocols in the United States are simply better tolerated than the other agents.
There are limitations to our study that deserve careful consideration. First, although our sample was large, its demographic characteristics do not reflect the broader population with OUD. Notably, our population was mostly male, predominantly White, and older on average. While we did control for these issues in our statistical models, the lack of diversity in this cohort is notable. That speaks to a broader concern regarding the possibility that some confounder unaccounted for in our adjustments affected our findings. For example, it is possible that patients who were able to remain adherent to MOUD were able to do so because they had access to transportation, greater family support, or a stronger desire to achieve recovery. However, our sensitivity analysis suggests that any unmeasured confounders would have to exert a very large effect to explain our findings. While additional sensitivity analyses in key strata or over shorter time periods would be reassuring, this is challenging even in our large cohort given the relative rarity of suicide mortality and the time-varying nature of our exposure. Thus, while we found associations between MOUD and decreased risk of mortality, including suicide mortality, we cannot conclude that the relationship is causal. Second, we had a limited set of measures of health care utilization. For example, it is possible that when patients were adherent to their MOUD, they were also adherent to treatments for other important conditions, such as depression. However, we did adjust for recent health care utilization. As these adjustments led to very little change in hazards for suicide mortality, it is unlikely that some other treatment with a broad protective effect was provided to patients in our cohort. Lastly, our study was conducted exclusively in patients who elected to at least start MOUD, so we can make no assertions regarding the effects of MOUD in those who do not seek treatment. There may be important differences in that group that make MOUD less associated with reduced suicide mortality.
Adjustment for important risk factors had little effect on our point estimates for MOUD. While this could indicate a large amount of measurement error, this seems unlikely. The more likely explanation for the stability of the estimates is the fact that this is largely a within-person analysis—the changes observed are longitudinal within a fixed cohort, with the comparison between periods and off MOUD.
These results again highlight the importance of providing medication treatment to as many people with opioid use disorder as possible. Doing so may have broad impacts on health outcomes, including suicide mortality. While our work focused on the effects of MOUD for those patients who receive it, rather than exploring efforts to promote use and spread, it is notable that the VA has engaged in several efforts to increase access to MOUD, including academic detailing and models to provide treatment in primary care clinics (
20,
21).
In summary, our results suggest that MOUD is associated with a decreased risk of suicide mortality, and the magnitude of decreased risk was about 50%. Furthermore, our work suggests a clear association of buprenorphine with the findings of decreased risk of suicide mortality and a possible effect for methadone, but no indication that naltrexone shares this effect.