- Open Access
A cohort study evaluating the association between concurrent mental disorders, mortality, morbidity, and continuous treatment retention for patients in opioid agonist treatment (OAT) across Ontario, Canada, using administrative health data
Harm Reduction Journal volume 17, Article number: 51 (2020)
Due to the high prevalence of mental disorders among people with opioid use disorder, the objective of this study was to determine the association between concurrent mental disorders, mortality, morbidity, and continuous treatment retention for patients in opioid agonist treatment in Ontario, Canada.
We conducted a retrospective cohort study of patients enrolled in opioid agonist treatment between January 1, 2011, and December 31, 2015. Patients were stratified into two groups: those diagnosed with concurrent mental disorders and opioid use disorder and those with opioid use disorder only, using data from the Ontario Health Insurance Plan Database, Ontario Drug Benefit Plan Database. The primary outcome studied was all-cause mortality using data from the Registered Persons Database. Emergency department visits from the National Ambulatory Care Database, hospitalizations Discharge Abstract Database, and continuous retention in treatment, defined as 1 year of uninterrupted opioid agonist treatment using data from the Ontario Drug Benefit Plan Database were measured as secondary outcomes. Encrypted patient identifiers were used to link information across databases.
We identified 55,924 individuals enrolled in opioid agonist treatment, and 87% had a concurrent mental disorder diagnosis during this period. We observed that having a mental disorder was associated with an increased likelihood of all-cause mortality (odds ratio (OR) 1.4; 95% confidence interval (CI) 1.2–1.5). For patients diagnosed with mental disorders, the estimated rate of ED visits per year was 2.25 times higher and estimated rate of hospitalization per year was 1.67 times higher than for patients with no mental disorders. However, there was no association between having a diagnosis of a mental disorder and 1-year treatment retention in OAT-adjusted hazard ratio (HR) = 1.0; 95% CI 0.9 to 1.1.
Our findings highlight the consequences of the high prevalence of mental disorders for individuals with opioid use disorder in Ontario, Canada.
The expanding opioid crisis in Canada is a complex issue that is exacerbated by factors such as increased exposure to prescription opioids , the prevalence of untreated mental illness , social isolation , a health system that does not work across silos, and largely unregulated advertising practices by drug companies [4, 5].
Research examining substance use and psychiatric comorbidities reports that approximately 50% of people with opioid use disorder (OUD) receiving treatment have a lifetime psychiatric diagnosis [6,7,8,9,10]. It has been shown that people who use opioids and psychiatric comorbidities have high all-cause mortality . In comparison to southern and urban communities, communities in northern Ontario [11,12,13,14,15], Canada, including First Nations, rural, and remote communities, generally experience high poverty, have access to limited infrastructure and health resources, demonstrate high risk-taking behaviors, and have less control over their environment specific to weather and occupation. Moreover, in Ontario, mental health services are seldom available or even coordinated with OUD services. Such factors pose an increased risk for mental disorders, substance use, and suicide . Given the confluence of emergent factors and the ongoing nature of the relationships between these factors, it is likely that challenges with opioids will continue in Canada.
Opioid agonist treatment (OAT) is currently the intervention with the best evidence for long-term treatment of OUD . However, retention in treatment continues to be a barrier to successful clinical outcomes for individuals with OUD. It is common for OAT patients to cycle through treatment and re-initiation of opioid consumption before they are stabilized in care, which can be dangerous because of changes to opioid tolerance in patients leading to a higher risk of mortality [18,19,20,21,22,23]. Therefore, despite the availability of evidence-based treatment for OUD, opioid-related deaths continue to be a critical issue in Ontario , the most populated province in Canada. For instance, from 2016 and 2018, 9078 Ontarians died of opioid poisoning . There has been a surge in opioid-related hospitalizations and emergency department (ED) visits [26,27,28], increasing from 9.42 per 100,000 population in 2003 to 19.55 per 100,000 population in 2015 .
According to the literature, mental disorders are prevalent among people with OUD [6, 7, 30]. Also, there are problems with the existing literature, including small sample sizes, restricted jurisdictions, and lack of longitudinal studies, which is an issue because OUD and mental disorders are chronic relapsing conditions [31, 32]. Also, previous studies have mainly focused on methadone as a form of OAT. However, buprenorphine prescriptions for the treatment of OUD are increasing, especially in rural areas of Ontario . Lastly, population-level issues such as mortality, ED visits, and hospitalizations [27, 34] and their relation to mental disorders have not been adequately researched.
The goal of this study was to highlight the effects of the chronic nature of opioid use disorder (OUD) and mental disorders on population health and address some of the gaps in the literature related to OUD and mental disorders. We evaluated the associations between concurrent mental disorders and mortality, morbidity, and continuous treatment retention in patients within OAT across Ontario. Our secondary objective was to examine regional variation in health service usage and treatment retention because health care services are materialized differently across Ontario, thus potentially resulting in varying OAT outcomes [27, 35, 36]. We hypothesized that OAT patients with mental disorders had poorer outcomes compared to OAT patients with no mental disorders.
We conducted a retrospective cohort study between January 1, 2011, and December 31, 2015, in Ontario. All health services examined in this study (except prescription drug) are covered by a government single-payer insurance program called the Ontario Health Insurance Plan (OHIP). Therefore, the study included all patients in Ontario who had accessed OAT within the study period. We used the date of the first opioid visit as the index date for all analyses. We designed the study around the index date, and follow-up times were fixed to eliminate bias related to the varying length of time in treatment for each patient. The last date of inclusion was December 31, 2015, and all dependent and independent variables were evaluated for 1 year after enrollment in OAT.
A treatment episode was defined as the time between the last and first OAT billing codes within a period of continuous retention in treatment (no interruptions care of more than 30 days). We used the first episode of OAT to identify patients, meaning that there was no previous history of OAT (including methadone or buprenorphine/naloxone) in the year before the first treatment episode. We chose to only include patients with no history of OAT in the year before the first episode during the study period to eliminate bias associated with cases involving multiple treatment attempts [23, 37, 38].
This study was approved by the Research Ethics Board of Laurentian University in Sudbury. This study is reported following the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines 
The study cohort was created by extracting OAT patients from the Ontario Drug Benefit Plan (ODB) database, which captures data on publicly funded pharmacy billing using drug identification numbers (DIN) (see Additional file 1) and with the (OHIP) database, which captures data on publicly funded physician-based health services using physician billing codes (see Additional file 2). In previous published ICES studies [27, 35, 40], the ODB database was used as the primary source to identify OAT patients. However, in Ontario, residents are only eligible for ODB public drug coverage if they are aged 65 years or older, reside in a long-term care facility, are disabled, are receiving social benefits for income support, or have high prescription drug costs relative to their net household income. Since 2011, new billing codes have helped to identify OAT services  in administrative databases because OHIP coverage is available to all permanent residents of Ontario. Therefore, to avoid excluding a subset of the population and risking selection bias, we used both ODB and OHIP databases to identify the primary patient cohort.
A set of exclusion criteria was applied to define the main cohort. We excluded all patients under 15 years of age (n = 2535 patients). Patients who were identified in the ODB with over 20% of their methadone dose in tablet formulation over 1 year were excluded because methadone in the tablet formulation is not approved for use as an addiction treatment in Canada (n = 5560 patients). Additionally, the following patients were excluded from the study: patients who were not eligible for OHIP (n = 437 patients) and non-Ontario residents (n = 427 patients). Ontario residents may not be eligible for OHIP coverage if they have been absent for more than 30 days during the first 6 months of residence in the province and if they were not in Ontario for at least 153 days in 12 months. Patients identified from ODB (n = 1383 patients), from OHIP (n = 30,124 patients), and patients who were identified in both databases (n = 24,417 patients) were combined to create the primary cohort (n = 55,924 patients). The steps used to create the primary cohort are outlined in Fig. 1.
We obtained anonymized patient-level data from Ontario publicly funded health services by submitting a formal requisition to ICES. ICES, previously known as the Institute for clinical evaluative sciences, is an independent research institute that collects and analyzes health care data for research. Patient information was linked anonymously across databases using encrypted ten-digit health card numbers. The linking protocol is used routinely for health system research in Ontario [42,43,44].
All diagnostic information from physician visits was determined using billing data from OHIP. ED visits were identified using the National Ambulatory Care Reporting System (NACRS). Hospital admissions were identified using the Discharge Abstract Database (DAD). We obtained patients’ location of residence and demographic information, including all-cause mortality from the Ontario Registered Persons Database (RPDB), which contains unique data for each resident who has ever received insured health.
Patients with a diagnosis of one or more mental disorders
Exposure to one or more mental disorders was assumed for all patients who had at least one of the OHIP diagnoses listed in Additional file 4. We excluded any substance-related diagnoses from our mental disorder diagnosis definition, which means that all patients in the mental disorders group had OUD and one or more mental disorders other than substance use disorders. We used OHIP diagnosis codes to identify patients in the mental disorders group. In the database, there is one general substance use diagnosis code. All substance use disorders that fall under code 304 are listed in Additional file 4. All patients in the cohort had that code based on their opioid dependence. Therefore, there was no way to detect those who also had a dependence on other substances such as cocaine, benzodiazepines, or other substances.
Patients were assigned to only one of the following groups: those diagnosed with a mental disorder other than OUD and those not diagnosed with a mental disorder other than OUD. The total time parameter to identify the diagnoses of mental disorders was 2 years for each patient. We identified patients who had a mental disorder diagnosis 1 year before the index date to 1 year after the index date (or the study end date). Mental disorders can be chronic, re-occurring conditions; therefore, the broad time frame was chosen to capture the condition accurately. Such wide time parameters are commonly used in other studies examining mental health comorbidities using administrative databases [45,46,47,48,49].
Study outcomes were defined based on the need to assess the association of all-cause mortality, frequent ED visits, and hospitalizations as a function of concurrent mental disorders and OUD, and not as an exposure leading to an event. All-cause mortality, ED visits, and hospitalizations have been used as indicators of complexity in the OUD population in other studies [7, 27, 50, 51]. Additionally, both frequent ED visits and hospitalizations are metrics used by health system planners and funders in Ontario to understand gaps in services in communities [50,51,52].
We requested the all-cause mortality variable from ICES as a dichotomized variable. At the time of the study, mortality-specific data were not available for the entirety of our study period. We used data from the RPDB database to calculate the number of days to death date from the study index date for each patient in the cohort to create the variable. If the patient had a mortality event anytime between their index date and the end of the study period (December 31, 2016), we assigned a code of 1 (all-cause mortality) or 0 (no all-cause mortality).
We used data from the NACRS database to identify ED visits. We considered a patient as having frequent ED visits if (1) contact with ED was after the index date and (2) a patient had ten or more ED events in a publicly funded Ontario hospital within 1 year. We used the diagnostic code that accompanied each ED event in the NACRS to categorize types of ED visits. Opioid-related, mental health-related, and reasons other than mental health or opioids were included in the analysis. Cutoffs for dichotomization were based on the frequency and distribution of ED visits for the entire cohort. The median number of ED visits was 14 per year, the mode was 6, and the mean was 22 (SD = 30.9). We chose to dichotomize with a value between the mean and mode (code 1 for ten or more ED visits and 0 for less than ten visits per year).
We used the DAD database to identify hospitalization. Hospitalizations were captured in three groups: opioid-related, mental health-related, and for reasons other than mental health or opioids using the primary diagnosis code that accompanied the hospitalization event in the DAD database. Hospitalizations were dichotomized and counted if a hospitalization discharge record appeared after a patient’s index date in a publicly funded Ontario hospital. The cutoff of one hospitalization was decided based on the frequency distribution of the number of hospitalizations for the cohort. The mean number of hospitalization per year was 3.
We conducted a subgroup analysis of 1-year treatment retention using the ODB database (n = 25,800). One-year treatment retention is correlated with a variety of positive health outcomes for patients, including reduced rates of drug use, criminal activity, and an increase in employment . One-year retention in OAT was assessed based on doses dispensed (from the ODB database). If the difference between the last and first days of dispensed medication within a period of continuous retention in treatment (no interruptions in prescribed doses > 30 days) was greater than 365 days, then the patient was considered to be retained for 1 year in OAT. Thirty days was chosen based on the use of this interval in previously published research [18, 53, 54]. The database used for medication dispensing in this study might not capture doses administered in a hospital or provincial correctional settings. However, in Ontario, patients will typically continue to receive methadone or buprenorphine in these settings. Since most hospital admissions or provincial incarcerations are less than 30 days, this approach allowed us to conduct the analysis without misinterpreting such events as treatment interruption.
The covariates available for the study included age; sex; location of residence; income quintile; human immunodeficiency virus infection (HIV); and deep tissue infection including endocarditis (OHIP diagnostic code 429), osteomyelitis (OHIP diagnostic code 730), and septic arthritis (OHIP diagnostic code 711).
Descriptive statistics summarize the characteristics of both the patient groups. We used standardized differences to compare categorical variables and the Wilcoxon rank-sum test to compare continuous variables between exposure groups. Standardized differences are not dependent on sample size; therefore, it was an appropriate measure to summarize characteristics of the large sample of patients in this study . Most researchers suggest that a value of standardized differences (d) ≥ 10% indicates a balance between groups . The group of patients not diagnosed with mental disorders was considered the baseline comparison. Covariates included age, sex, location of residence, income quintile, HIV status, and deep tissue infection.
A logistic regression model was used to test the association between mental disorders and all-cause mortality. We adjusted the model using the covariates and calculated unadjusted and adjusted odds ratios (OR) and 95% confidence intervals (CI). If the OR was equal to one, having a mental disorder did not affect the likelihood of all-cause mortality. If the OR was more than one, having a mental disorder was associated with the likelihood of all-cause mortality. If the OR was less than one, having a mental disorder was associated with a decreased likelihood of all-cause mortality. If the CI spanned one, it was evidence of a lack of association between the exposure and the outcome .
Negative binomial regression models were used to estimate the association between mental disorders, ED visits, and hospitalizations. In a study comparing different models to evaluate health service utilization, the negative binomial regression model was found to be one of the most appropriate because it can address over dispersion that is typically observed in health care utilization data . We used the regression coefficient to calculate OR and 95% CI [58, 59].
A Cox Proportional Hazard Model was used to characterize the risk of treatment discontinuation over time for both methadone and buprenorphine/naloxone for the subgroup of patients identified through the ODB database (n = 28,000). Due to the time frame of the study, we assumed a constant hazard ratio between groups over time. We adjusted the model using the baseline covariates and calculated hazard ratio (HR) and 95% CI. If the HR was equal to one, having a mental disorder did not affect the likelihood of 1-year treatment retention. If the OR was more than one, having a mental disorder was associated with the likelihood of 1-year treatment retention. If the OR was less than one, having a mental disorder was associated with a decreased likelihood of 1-year treatment retention. If the CI spanned one, it was evidence of a lack of association between the exposure and the outcome. All statistical analyses were conducted from the secure server using SAS Version 9.4 . Data was reviewed by ICES to insure privacy standards were met.
A total of 55,924 patients who were enrolled in OAT between 2011 and 2015 were identified. Of those, 48,679 (87.0%) had one or more diagnoses of mental disorder. We followed all patients for 1 year, and there were 2712 (8.2%) deceased patients during follow-up. There were 534 records with missing information on income quintile and 3 records with missing information on the location of residence. Missing income records were re-assigned to the lowest income group, and those with missing information on the location of residence were deleted.
Prevalence of mental disorders
Anxiety disorders, including obsessive compulsive disorder and other anxiety-related disorders (60%) and mood disorders, including bipolar and other depressive-related disorders (20%), were among the most prevalent in the cohort. The results are presented in Fig. 2.
Males represented a lower proportion of individuals in the mental disorders group (63.0%) compared to the group with no mental health disorder (77.0%, (p < 0.001)). The proportion of individuals aged 15–24 in the mental disorders and no mental disorders groups was 17.9% and 20.9%, respectively. Of the individuals with mental disorders, 77.8% lived in southern urban regions, 4.6% in northern rural areas, 9.3% in northern urban areas, and 8.2% in southern rural areas. This compared to individuals who had not been diagnosed with a mental disorder, 71.9% lived in southern urban areas, 21.0% lived in northern rural areas, 10.5% in northern urban areas, and 9.4% in southern rural areas (sd 0.2). We observed significant differences between the two groups with regards to the proportion of patients with HIV (mental disorders group = 0.8% and no mental disorders group = 0.3%, sd 0.1) and deep-tissue infections (mental disorders group = 3.3% and no mental disorders group = 1.3%, p 0.1). Retention in buprenorphine treatment was significantly different between groups (mean days engaged in buprenorphine treatment was 55.0, standard deviation (SD) = 167.8, and 34.8, SD = 121.8). However, there was no statistically significant difference identified in the mean days retained in methadone treatment: 356.8, SD = 441.0 and 344.6, SD = 409.8 for the mental disorders group and no mental disorders group respectively on a subgroup of patients identified from the ODB database (n = 28,500). The patient characteristics are outlined in Table 1. Importantly, there was a significant difference in the proportion of patients with an all-cause mortality event between groups (mental disorders group = 5.1% and no mental disorders group = 3.1%, sd 0.1) and significant differences in the proportion of patients retained for 1 year in OAT (mental disorders group = 37.6%, and no mental disorders group = 35.7%, sd 0.2). The patient characteristics are outlined in Table 1.
Having a diagnosis of one or more mental disorders was associated with an increase in the likelihood of all-cause mortality (adjusted OR (aOR) = 1.4; 95% confidence interval (CI) 1.2 to 1.5) when compared to patients who had not been diagnosed with mental disorders. For patients diagnosed with mental disorders, the estimated rate of ED visits per year was 2.25 times higher than for patients with no mental disorders (factor exponent (0.81) =2.25). For patients who had been diagnosed with mental disorders, the estimated rate of hospitalization per year was 1.67 times higher than for patients with no mental disorders (factor exponent (0.51) = 1.67). The results are presented in Table 2.
We also found that there was no association between having a diagnosis of a mental disorder and 1-year retention in methadone maintenance treatment (adjusted hazard ratio (HR) = 1.0; 95% CI 0.9 to 1.1) in the subgroup of ODB patients. However, having a mental disorder was associated with an increased likelihood of 1-year retention in buprenorphine maintenance treatment (HR = 1.3, 95% CI 1.2–1.3) in the subgroup of ODB patients. The results are presented in Fig. 3.
The highest proportion of patients who died was in southern urban areas (n = 2039, 5.38%) and the lowest in northern rural areas (n = 68, 3.04%, p < 0.001), whereas the highest proportion of patients with frequent ED visits was in northern rural areas (n = 1880, 83.93%) and the lowest in southern urban areas (n = 29,137, 55.81%, p < 0.001) (Fig. 4). Similarly, the highest proportion of patients with hospitalizations was in northern rural areas (n = 1917, 85.58%) and the lowest was in southern urban areas (n = 36,032, 71.72%, p < 0.001). The highest proportion of patients retained for 1 year was in southern rural areas (n = 347, 37.1%), followed by northern urban (n = 456, 36.5%), southern urban (n = 2630, 35.1%), and lastly northern rural (n = 173, 29.1%, p < 0.001). Results are described in Fig. 4.
The findings from this study indicate that patients in OAT with mental disorders have a higher likelihood of mortality, ED visits, and hospitalizations and that patients who live in northern and rural areas of Ontario have a more complex profile of health care utilization than patients who live in southern Ontario. Mortality, ED visits, and hospitalizations are important metrics used to measure health care delivery. Hospitals are acute care facilities that are not designed to treat chronic diseases such as substance use and mental disorders. Therefore, if there are high numbers of patients with OUD going to the hospital, it can be an indication that there are issues with service access in communities.
The results of this study align with the current studies which demonstrate a high prevalence of mental disorders in the OUD population [7, 30, 61]. However, in this study, 87% of the population had a mental disorder, whereas the prevalence in the literature is approximately 50% [7, 30, 61]. The high prevalence identified in our cohort may be because we examined a broad range of disorders and that our time parameter, including study time and follow-up, was much longer than in other studies [6, 7]. The study period was also much more recent and may be more representative of the current issues in Canada relating to patients in OAT. Moreover, the definition of enrollment was different from that of the previous studies because all diagnostic information from physician visits was determined using billing data from OHIP. Similar to other studies, we also found that mental disorders in the OAT population are associated with a more complex profile of health service utilization [7, 62,63,64].
In our study, having one or more mental disorders was associated with an increased likelihood of all-cause mortality when compared to the group that had not been diagnosed with mental disorders. This finding is consistent with a study by Saunders et al. suggesting that individuals with OUD who have psychiatric symptoms have higher rates of overdose, decreased quality of life, and higher rates of continual substance use . It is important to note that the results indicate that individuals with mental health comorbidities are a complex group of patients; however, results should be interpreted with caution. It is possible that either the use of opioids is especially harmful for individuals with mental disorders or that the use of opioids might aggravate a pre-existing condition . For many, OUD is a lifelong illness associated with serious chronic health and social outcomes . Dichotomizing all-cause mortality allowed us to evaluate death associated with chronic mental disorders and OUD and not necessarily as concurrent disorders leading to an event. We believe that future studies should be conducted to evaluate deaths as a function of time in treatment to explore this issue further.
In addition to the high rate of mortality, we observed an association between concurrent mental disorders and high rates of ED visits and hospitalizations. This significant association is not surprising given that mental disorders have been used as a marker for greater clinical complexity by some scholars [6, 68]. The findings of this study do not show causation but associations which suggest that there may be ways to prevent serious adverse events by addressing mental health issues in concert with OUD.
Despite the differences in mortality and morbidity between the exposure groups, having one or more mental disorders was not associated with 1-year retention in methadone but an increased likelihood of retention in buprenorphine/naloxone. Some authors report that mental disorders increase retention [38, 62], and some have reported no effect . It is important to note that buprenorphine/naloxone patients had lower retention rates than methadone patients (mean of approximately 40 days retained on buprenorphine vs. mean of 350 days retained in methadone) in this study and other studies in the literature [69, 70]. We believe that research exploring concepts of severity of dependence, resilience, and motivation would be critical to explain this relationship further.
We identified a trend where patients in the cohort residing in southern urban regions of Ontario had the highest prevalence of all-cause mortality. A cohort study by Gomes et al. on geographic variation in opioid prescribing and opioid-related mortality in Ontario  found that communities with some of the highest rates of opioid-related death were in northern Ontario. One possible factor accounting for higher mortality rates in southern urban regions is the unprecedented increased fentanyl-related deaths in recent years, which started in urban centers and is now making its way into smaller communities . Another possible explanation is that barriers to access OAT are higher in northern Ontario, perhaps resulting in a larger number of deaths among individuals who were never enrolled in OAT (and thus not included in this study). It is also important to note that there may be an increased likelihood of mortality and morbidity closer to 2016 in all regions due to the increased of fentanyl in the drug supply .
Although the results show a lower proportion of individuals who died in northern regions, there was a significantly higher proportion of morbidity measured by ED visits and hospitalizations in rural northern Ontario. Possible factors accounting for the high prevalence of acute care use may be attributed to the well-known limited availability of specialist addiction medicine and psychiatric services in northern rural regions of Ontario [36, 73]. Lastly, with regards to geographical variation, our results demonstrated that the proportion of patients retained for 1 year of OAT treatment was highest for those living in southern rural areas. This finding is counter to the findings published by Eibl et al. , which demonstrated that patients living in northern rural areas were more likely to be retained for 1 year. Notably, our study is more recent than the earlier study by Eibl et al., and currently, there is a broader availability of buprenorphine/naloxone. Further research is needed on region-specific retention for individuals with concurrent mental disorders and OUD, including analysis of changes over time and the impact of changing rates of methadone and buprenorphine/naloxone prescribing.
One of the potential issues inherent to any study of this type is that health administrative data were not collected to do research, which may have led to misclassification of disease prevalence and clinical outcomes. For instance, we were not able to evaluate patients who may have a mental disorder or OUD who have not yet sought out services . We did not include the Ministry of Health-funded community mental health and addiction services and federally funded health services, such as mental health and addiction services provided in First Nation communities, as well as any other mental health and addiction services funded by provincial ministries other than the Ministry of Health. Moreover, since there is a well-known lack of physicians available to diagnose and treat patients in rural areas, we may have underestimated the prevalence of mental disorders and other diagnoses.
With relation to mental diagnoses, further stratification of diagnosis type would have been of interest but was not possible as the numbers became too small to publish for the outcomes of interest. It is essential to highlight the potential underestimation of the prevalence of post-traumatic stress disorder (PTSD) since it is known to be often misdiagnosed as another anxiety, depression or other related disorders [75,76,77,78]. Details such as years of drug use, the amount or type of opioid used, the history of mental health services before 1 year before the first episode of OAT, and the number of times patients were in and out of OAT after their first episode of OAT in this study remain unknown. Additionally, although hepatitis C is highly prevalent in the OAT population and can impact outcomes, data on hepatitis C was not available for this study. In the last few years, there have been several changes concerning hepatitis C treatment availability and coverage in Ontario. Since we used physician billing data and prescription data to ascertain diagnoses in this study, we were not able to define an accurate hepatitis C diagnosis during the time of the study due to all the changes in hepatitis C treatment. Since this study is observational in nature, causality cannot be inferred. Additionally, there is a potential for associations to be significant by chance due to the large sample size.
There are also limitations associated with the way we measured the outcomes for this study. Some authors have stated that creating categories from continuous data can lead to loss of information . More specifically, there are limitations to using all-cause mortality as a dichotomized variable without censoring for time. For instance, for those patients where all-cause mortality occurred during the first year of OAT, the likelihood of 1-year retention is reduced. However, the function of time is a modest bias since we observed an increase in ED visits and hospitalizations in the groups with the highest mortality. The results of this study, which were found to be statistically significant, must be interpreted critically within the context of the population of interest to determine whether the results have a clinical or health system impact.
The outcomes of this study have important implications for those involved in health care planning and policy development because our data suggest that the prevalence of mental disorders in the OAT population is alarmingly high and that mental disorders are associated with severe consequences. Currently, the regulations and model of care for OAT in Ontario promotes access to services but do not incentivize efforts towards coordination with other parts of the health care system, including other addiction or mental health services. Results may be generalizable in regions where OAT programs and health care regulations are similar to those in Ontario. Further study is needed to determine the effectiveness of concurrent delivery of mental health, other substance use, and OAT services.
Availability of data and materials
The datasets used during the current study are not publicly available due to privacy reasons, but aggregated data are included in this published article and its supplementary information files.
Opioid use disorder
Opioid agonist treatment
Formally known as the Institute for Clinical Evaluative Sciences
Ontario drug benefit plan
Statistical Analytics Software
Ontario Health Insurance Plan
Canadian Institute for Health Information
National Ambulatory Care Reporting System
Discharge Abstract Database
Registered Persons Database
Local Health Integration Network
Adjusted odds ratio
International Classification of Disease
Post-traumatic stress disorder
Human immunodeficiency virus
Ontario Drug Policy Research Network
Canadian Institute for Health Information: Pan-Canadian trends in the prescribing of opioids, 2012 to 2016 Ottawa ON 2017 [Available from: https://www.cihi.ca/sites/default/files/document/pan-canadian-trends-opioid-prescribing-2017-en-web.pdf].
Van Ameringen M, Mancini C, Patterson B, Boyle MH. Post-traumatic stress disorder in Canada. CNS Neurosci Ther. 2008;14(3):171–81.
WHO. World Health Organization. Principles of drug dependence treatment Geneva, Switzerland: World Health Organization; 2008 [Available from: http://www.who.int/substance_abuse/publications/principles_drug_dependence_treatment.pdf?ua=1].
Public Health Agency of Canada. Notice of intent to restrict the marketing and advertising of opioids Ottawa Ontario 2018 [Available from: https://www.canada.ca/en/health-canada/services/drugs-health-products/drug-products/announcements/restrict-advertising-opioids.html].
Webster P. Oversight of opioid advertising in Canada remains lax despite new federal guidance. CMAJ. 2018;190(23):E730–E1.
Astals M, Domingo-Salvany A, Buenaventura CC, Tato J, Vazquez JM, Martin-Santos R, et al. Impact of substance dependence and dual diagnosis on the quality of life of heroin users seeking treatment. Subst Use Misuse. 2008;43(5):612–32.
Bogdanowicz KM, Stewart R, Broadbent M, Hatch SL, Hotopf M, Strang J, et al. Double trouble: psychiatric comorbidity and opioid addiction-all-cause and cause-specific mortality. Drug Alcohol Depend. 2015;148:85–92.
MLD G, Dodd DR, Potter JS. Baseline characteristics and treatment outcomes in prescription opioid dependent patients with and without co-occurring psychiatric disorder. Am J Drug Alcohol Abuse. 2015;40:157–62.
Saunders EC, McGovern MP, Lambert-Harris C, Meier A, McLeman B, Xie H. The impact of addiction medications on treatment outcomes for persons with co-occurring PTSD and opioid use disorders. Am J Addict. 2015;24(8):722–31.
Brooner RK, King VL, Kidorf M, Schmidt CW Jr, Bigelow GE. Psychiatric and substance use comorbidity among treatment-seeking opioid abusers. Arch Gen Psychiatry. 1997;54(1):71–80.
Brannen CD, Dyck K, Hardy, C, Mushquash, C. Rural mental health services in Canda: a model for research and practice. In: Kulig JW, A., editor. Health in Rural Canada. Vancouver, BC, Canada: UBC Press; 2012. p. 239-255.
Canadian Institute for Health Information. Canadian Population Health Initiative. How healthy are rural Canadians? An assessment of their health status and health determinants. Canadian Institute for Health Information; 2006.
Ward M. Short report#5: mental health in Northern Ontario. Northern Health Information Partnership; 2005.
Minore B, M. Boone, M.E. Hill, A. Lyubechansky, R. Kuzik, and C. Macdonald. Towards using family health teams to care for chronic mental health clients in Northern Ontario: Barriers, Opportunities and Training Needs. Centre for rural and Northern Health Research, Lakehead University; 2005.
J. H. Rural health and the health of rural communities. 1998.
World Health Organization: Prevention of mental disorders. Summary Report Geneva Switzerland: WHO; [Available from: http://www.who.int/mental_health/evidence/en/prevention_of_mental_disorders_sr.pdf].
Mattick RP, Breen C, Kimber J, Davoli M. Methadone maintenance therapy versus no opioid replacement therapy for opioid dependence. Cochrane Database Syst Rev. 2009;3:CD002209.
Nosyk B, Marsh DC, Sun H, Schechter MT, Anis AH. Trends in methadone maintenance treatment participation, retention, and compliance to dosing guidelines in British Columbia, Canada: 1996-2006. J Subst Abus Treat. 2010;39(1):22–31.
Bell J, Burrell T, Indig D, Gilmour S. Cycling in and out of treatment; participation in methadone treatment in NSW, 1990-2002. Drug Alcohol Depend. 2006;81(1):55–61.
Galai N, Safaeian M, Vlahov D, Bolotin A, Celentano DD, Study A. Longitudinal patterns of drug injection behavior in the ALIVE study cohort,1988-2000: description and determinants. Am J Epidemiol. 2003;158(7):695–704.
Dobler-Mikola A, Hattenschwiler J, Meili D, Beck T, Boni E, Modestin J. Patterns of heroin, cocaine, and alcohol abuse during long-term methadone maintenance treatment. J Subst Abus Treat. 2005;29(4):259–65.
Bovasso G, Cacciola J. The long-term outcomes of drug use by methadone maintenance patients. J Behav Health Serv Res. 2003;30(3):290–303.
Termorshuizen F, Krol A, Prins M, Geskus R, van den Brink W, van Ameijden EJ. Prediction of relapse to frequent heroin use and the role of methadone prescription: an analysis of the Amsterdam Cohort Study among drug users. Drug Alcohol Depend. 2005;79(2):231–40.
Hedegaard H, Warner, M., Miniño, A. M. Drug overdose deaths in the United States, 1999–2016 NCHS Data Brief No. 294 2017 [Available from: https://www.cdc.gov/nchs/products/databriefs/db294.htm].
Public Health Agency of Canada. Statement from the co-chairs of the Special Advisory Committee on the Epidemic of Opioid Overdoses on Updates to Opioid-Related Mortality Data. Ottawa, ON 2018 [Available from: https://www.canada.ca/en/public-health/news/2018/03/statement-from-the-co-chairs-of-thespecial-advisory-committee-on-the-epidemic-of-opioid-o.html].
Gomes TGS, Martins D. Latest trends in opioid-related deaths in Ontario: 1991 to 2015 Toronto: Ontario Drug Policy Research Network; 2017 [Available from: http://odprn.ca/wp-content/uploads/2017/04/ODPRN-Report_Latest-trends-in-opioid-related-deaths.pdf].
Martins DGS, Tadrous M, Paterson M, Bandola D, Singh S, Juurlink D, Mamdani M, Gomes T. Opioid use and related adverse events in Ontario. Ontario Drug Policy Research Network Toronto Ontario: ODPRN; 2016 [Available from: http://odprn.ca/wp-content/uploads/2016/11/ODPRN-Opioid-Use-and-Related-Adverse-Events-Nov-2016.pdf].
Canadian Institute for Health Information, Canadian Centre on Substance Abuse. Hospitalizations and emergency department visits due to opioid poisoning in Canada. Ottawa ON: CIHI; 2016 [Available from: https://secure.cihi.ca/free_products/Opioid%20Poisoning%20Report%20%20EN.pdf].
Kingston Frontenac and Addington Public Health Unit Kingston, Ontario, Canada: KFL&A Public Health Informatics [Available from: https://public.tableau.com/profile/kflaphi#!/vizhome/OntarioOpioidSurveillanceMonitor/ACESEDVisits].
Andrews G, Henderson S, Hall W. Prevalence, comorbidity, disability and service utilisation. Overview of the Australian National Mental Health Survey. Br J Psychiatry. 2001;178:145–53.
Amato L, Minozzi S, Davoli M, Vecchi S. Psychosocial combined with agonist maintenance treatments versus agonist maintenance treatments alone for treatment of opioid dependence. Cochrane Database Syst Rev. 2011;10:CD004147.
Dugosh K, Abraham A, Seymour B, McLoyd K, Chalk M, Festinger D. A systematic review on the use of psychosocial interventions in conjunction with medications for the treatment of opioid addiction. J Addict Med. 2016;10(2):91–101.
Mamakwa S, Kahan M, Kanate D, Kirlew M, Folk D, Cirone S, et al. Evaluation of 6 remote First Nations community-based buprenorphine programs in northwestern Ontario: retrospective study. Can Fam Physician. 2017;63(2):137–45.
Gomes T, Mamdani MM, Dhalla IA, Cornish S, Paterson JM, Juurlink DN. The burden of premature opioid-related mortality. Addiction. 2014;109(9):1482–8.
Gomes T, Juurlink D, Moineddin R, Gozdyra P, Dhalla I, Paterson M, et al. Geographical variation in opioid prescribing and opioid-related mortality in Ontario. Healthcare Quarterly. 2011;14(1):22–4.
Kiepek N, Groom B, Toppozini D, Kakekagumick K, Muileboom J, Kelly L. Evaluation of an inpatient medical withdrawal program in rural Ontario: a 1-year prospective study. Can J Rural Med. 2015;20(3):92–7.
Volkow ND. America’s addiction to opioids: heroin and prescription drug abuse [Available from: https://www.nih.gov/sites/default/files/institutes/olpa/20140514-senate-testimony-volkow.pdf].
Nosyk B, MacNab YC, Sun H, Fischer B, Marsh DC, Schechter MT, et al. Proportional hazards frailty models for recurrent methadone maintenance treatment. Am J Epidemiol. 2009;170(6):783–92.
von Elm E, Altman DG, Egger M, Pocock SJ, Gotzsche PC, Vandenbroucke JP, et al. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. Rev Esp Salud Publica. 2008;82(3):251–9.
Eibl JK, Gomes T, Martins D, Camacho X, Juurlink DN, Mamdani MM, et al. Evaluating the effectiveness of first-time methadone maintenance therapy across northern, rural, and urban regions of Ontario. Canada J Addict Med. 2015;9(6):440–6.
Ministry of Health and Long Term Care. Opioid Agonist Maintenance Program (OAMP) monthly management fee and point of care drug testing health services Branch; 2012 [Available from: http://www.health.gov.on.ca/en/pro/programs/ohip/bulletins/4000/bul4548.pdf].
Mamdani M, Rochon P, Juurlink DN, Anderson GM, Kopp A, Naglie G, et al. Effect of selective cyclooxygenase 2 inhibitors and naproxen on short-term risk of acute myocardial infarction in the elderly. Arch Intern Med. 2003;163(4):481–6.
Juurlink DN. Proton pump inhibitors and clopidogrel: putting the interaction in perspective. Circulation. 2009;120(23):2310–2.
Juurlink DN, Gomes T, Lipscombe LL, Austin PC, Hux JE, Mamdani MM. Adverse cardiovascular events during treatment with pioglitazone and rosiglitazone: population based cohort study. BMJ. 2009;339:b2942.
McConville P, Walker NP. The reliability of case register diagnoses: a birth cohort analysis. Soc Psychiatry Psychiatr Epidemiol. 2000;35(3):121–7.
Frayne SM, Halanych JH, Miller DR, Wang F, Lin H, Pogach L, et al. Disparities in diabetes care: impact of mental illness. Arch Intern Med. 2005;165(22):2631–8.
Frayne SM, Miller DR, Sharkansky EJ, Jackson VW, Wang F, Halanych JH, et al. Using administrative data to identify mental illness: what approach is best? Am J Med Qual. 2010;25(1):42–50.
Mental Health Commission of Canada. Overview of mental health data in Canada: background, needs, and gaps, Calgary, AB 2014 [Available from: www.mentalhealthcommission.ca.
Sajatovic M, Blow FC, Ignacio RV. Psychiatric comorbidity in older adults with bipolar disorder. Int J Geriatr Psychiatry. 2006;21(6):582–7.
Durbin J LE, Rush B. A study of hospital emergency service use, crisis service delivery and police response after mental health system enhancements. 2010.
Genell Andren K, Rosenqvist U. Heavy users of an emergency department--a two year follow-up study. Soc Sci Med. 1987;25(7):825–31.
Althaus F, Paroz S, Hugli O, Ghali WA, Daeppen JB, Peytremann-Bridevaux I, et al. Effectiveness of interventions targeting frequent users of emergency departments: a systematic review. Ann Emerg Med. 2011;58(1):41–52 e42.
Peles E, Linzy S, Kreek M, Adelson M. One-year and cumulative retention as predictors of success in methadone maintenance treatment: a comparison of two clinics in the United States and Israel. J Addict Dis. 2008;27(4):11–25.
Franklyn AM, Eibl JK, Gauthier G, Pellegrini D, Lightfoot NK, Marsh DC. The impact of benzodiazepine use in patients enrolled in opioid agonist therapy in Northern and rural Ontario. Harm Reduct J. 2017;14(1):6.
Yang D DJ. A unified approach to measuring the effect size between two groups using SAS 2012 [Available from: https://support.sas.com/resources/papers/proceedings12/335-2012.pdf].
Austin PC. An introduction to propensity score methods for reducing the effects of confounding in observational studies. Multivar Behav Res. 2011;46(3):399–424.
Szumilas M. Explaining odds ratios. J Can Acad Child Adolesc Psychiatry. 2010;19(3):227–9.
Moineddin R, Meaney C, Agha M, Zagorski B, Glazier RH. Modeling factors influencing the demand for emergency department services in Ontario: a comparison of methods. BMC Emerg Med. 2011;11:13.
Pedan, A.. Analysis of count data using the SAS®System Tewksbury, MA SAS; [Available from: https://support.sas.com/resources/papers/proceedings/proceedings/sugi26/p247-26.pdf.
SAS version 9.4 Cary, North Carolina [Available from: https://www.sas.com/en_us/software/sas9.html].
Astals M, Diaz L, Domingo-Salvany A, Martin-Santos R, Bulbena A, Torrens M. Impact of co-occurring psychiatric disorders on retention in a methadone maintenance program: an 18-month follow-up study. Int J Environ Res Public Health. 2009;6(11):2822–32.
Gelkopf M, Weizman T, Melamed Y, Adelson M, Bleich A. Does psychiatric comorbidity affect drug abuse treatment outcome? A prospective assessment of drug abuse, treatment tenure and infectious diseases in an Israeli methadone maintenance clinic. Isr J Psychiatry Relat Sci. 2006;43(2):126–36.
McLellan ATL, Luborsky L, Woody GE, O’Brien CP, Druley KA. Predicting response to alcohol and drug abuse treatments. Role of psychiatric severity. Arch Gen Psychiatry. 1983;40(6):620–5.
Verthein U, Degkwitz P, Haasen C, Krausz M. Significance of comorbidity for the long-term course of opiate dependence. Eur Addict Res. 2005;11(1):15–21.
Saunders EC, McGovern MP, Lambert-Harris C, Meier A, McLeman B, Xie H. The impact of addiction medications on treatment outcomes for persons with co-occurring PTSD and opioid use disorders. Am J Addict. 2015.
Rounsaville BJ, Kranzler HR, Ball S, Tennen H, Poling J, Triffleman E. Personality disorders in substance abusers: relation to substance use. J Nerv Ment Dis. 1998;186(2):87–95.
Hser YI, Hoffman V, Grella CE, Anglin MD. A 33-year follow-up of narcotics addicts. Arch Gen Psychiatry. 2001;58(5):503–8.
O'Brien CP, Woody GE, McLellan AT. Psychiatric disorders in opioid-dependent patients. J Clin Psychiatry. 1984;45(12 Pt 2):9–13.
Gryczynski J, Mitchell SG, Jaffe JH, Kelly SM, Myers CP, O'Grady KE, et al. Retention in methadone and buprenorphine treatment among African Americans. J Subst Abus Treat. 2013;45(3):287–92.
Hser YI, Saxon AJ, Huang D, Hasson A, Thomas C, Hillhouse M, et al. Treatment retention among patients randomized to buprenorphine/naloxone compared to methadone in a multi-site trial. Addiction. 2014;109(1):79–87.
Agency of Health Protection and Promotion. Opioid related morbidity and mortality Toronto Ontario [Available from: https://www.publichealthontario.ca/en/dataandanalytics/pages/opioid.aspx#/drug].
Opioid-related morbidity and mortality in Ontario [Internet] 2017 [Available from: https://www.publichealthontario.ca/en/dataandanalytics/pages/opioid.aspx#/trends.
Kurdyak P, Stukel TA, Goldbloom D, Kopp A, Zagorski BM, Mulsant BH. Universal coverage without universal access: a study of psychiatrist supply and practice patterns in Ontario. Open Med. 2014;8(3):e87–99.
Davis KA, Sudlow CL, Hotopf M. Can mental health diagnoses in administrative data be used for research? A systematic review of the accuracy of routinely collected diagnoses. BMC Psychiatry. 2016;16:263.
Meltzer EC, Averbuch T, Samet JH, Saitz R, Jabbar K, Lloyd-Travaglini C, et al. Discrepancy in diagnosis and treatment of post-traumatic stress disorder (PTSD): treatment for the wrong reason. J Behav Health Serv Res. 2012;39(2):190–201.
Nutt DS, M, Zohar, J. Post traumatic stress disorder: diagnosis, management and treatment: CRC Press; 2009.
Weinstein D, Staffelbach D, Biaggio M. Attention-deficit hyperactivity disorder and posttraumatic stress disorder: differential diagnosis in childhood sexual abuse. Clin Psychol Rev. 2000;20(3):359–78.
Vickerman KA, Margolin G. Post-traumatic stress in children and adolescents exposed to family violence: II. Treatment Prof Psychol Res Pr. 2007;38(6):620.
Naggara O, Raymond J, Guilbert F, Roy D, Weill A, Altman DG. Analysis by categorizing or dichotomizing continuous variables is inadvisable: an example from the natural history of unruptured aneurysms. AJNR Am J Neuroradiol. 2011;32(3):437–40.
We thank IC/ES Data Analytic Services, more specifically, Ryan Ng, for his assistance with data extraction and database set up. We also thank members of the Patient and Family Advisory Committee for sharing their stories and helping to guide the research project. Lastly, we would like to thank Matt Parotta and Lyne Morvan for editing the final draft of the paper.
We thank our funders (Northern Ontario Academic Medical Association) through the Academic Funding Plan Innovation Fund for their support in completing this project, grant # A-17-05.
Ethics approval and consent to participate
This study was approved by the Research Ethics Board of Laurentian University in Sudbury, Ontario, Canada.
Consent for publication
This study only used administrative health data, there was no primary data collection; therefore, consent for publication from individual patients did not apply.
Dr. David Marsh maintains the following roles: Chief Medical Director at CATC (Canadian Addiction Treatment Center), opioid agonist therapy provider. Dr. Marsh has no ownership stake in the CATC as a stipendiary employee. We do not foresee any conflict of interest as data will be made freely available to the public, and neither the CATC nor the Universities prevent publication and dissemination of knowledge. The authors have no conflicts declared. This does not alter our adherence to the Harm Reduction Journal’s policies on sharing data and materials
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Drug Identification Numbers
Addiction Medicine Fee Codes
Mental Health Service Fee Codes
ICD9 and 10 Codes
Study outcomes related to concurrent mental disorders and all-cause mortality for OAT patients (Regression Output)
About this article
Cite this article
Morin, K.A., Eibl, J.K., Gauthier, G. et al. A cohort study evaluating the association between concurrent mental disorders, mortality, morbidity, and continuous treatment retention for patients in opioid agonist treatment (OAT) across Ontario, Canada, using administrative health data. Harm Reduct J 17, 51 (2020). https://doi.org/10.1186/s12954-020-00396-x
- Opioid use disorder
- Opioid agonist treatment
- Mental disorders