Mortality in the Melbourne injecting drug user cohort study (MIX)
Harm Reduction Journal volume 12, Article number: 55 (2015)
There are few studies of mortality amongst people who inject drugs (PWID) in Australia. In this study, we estimate mortality in a cohort of PWID in Melbourne and examine predictors of mortality including health service use, demographic characteristics, drug use and personal wellbeing.
We linked identifiers from the Melbourne injecting drug use cohort study (MIX; n = 655) to the National Death Index from 2008 to 2012 to estimate standardised mortality ratios (SMRs). Cox regression was used to examine the bivariate relationship between exposures determined at baseline and subsequent mortality. There were 24 (3.6 %) deaths over the study period. The mortality rate in the cohort was 1.0 per 100 PY (95 % CI 0.71–1.57), with an SMR of 17.3 (95 % CI 11.6–25.8). Baseline reports of four or more lifetime incarcerations (HR 3.65, 95 % CI 1.16–11.52), past month ambulance attendance (HR 4.43, 95 % CI 1.76–11.17), past month emergency department presentation (HR 3.44, 95 % CI 1.47–8.03) and past 6-month self-reported heroin overdose (HR 3.14, 95 % CI 1.24–7.96) were associated with increased mortality risk.
Contact with emergency services, particularly for drug overdose, remains a lost opportunity to provide referrals for harm reduction and naloxone training programmes to PWID at greater risk of mortality.
Estimates of mortality among people who inject drugs (PWID) in developed countries (based mostly on treatment samples ) range between 0.5 and 5 deaths per 100 person-years (PY). Stoové et al.  reported a mortality rate of 0.8 per 100 PY among PWID in Melbourne between 1990 and 2006, with most deaths related to drug overdose and other injecting-related causes; unlike other countries, HIV mortality was uncommon .
A range of factors increase mortality risk among PWID, including having a history of non-fatal overdose events . PWID who overdose often come into contact with ambulance and emergency department (ED) services; however, the relationship between the use of these and other health services with mortality has not been previously explored.
We provide updated estimates of mortality among PWID in Australia where HIV prevalence is low and targeted health services for PWID are available . We describe the relationship between mortality and demographic characteristics, patterns of drug use, measures of personal wellbeing and use of health services.
Baseline data from the Melbourne injecting drug user cohort study (MIX) collected in Melbourne between November 2008 and March 2010 was linked to Australia’s National Death Index (NDI). MIX participants were recruited through a combination of respondent-driven sampling, snowball sampling and outreach. Eligible participants were aged 18 and over, reported injecting regularly for at least 6 months and had a valid Medicare card. The baseline MIX questionnaire included questions on socio-demographics, drug use, health service use, and self-reported health and wellbeing. Further details of the recruitment methods and the wider study are described elsewhere . Probabilistic data linkage of participant identifier codes based on letters from first and last names, gender and date of birth was used to link NDI data, with matches validated using full first names, deaths known to the study team and date of last contact. Deaths up to December 2012 were captured across all states and territories. Cause-of-death codes are reported where available.
The study was approved by the Human Research Ethics Committee of the Victorian Department of Health and the Monash University Human Research Ethics Committee, and informed consent was obtained from all participants.
Person-years (PY) were calculated as the time from baseline interview until death or censorship at 31st December 2012. Mortality rates were calculated using the person-time method  and disaggregated by gender. Standardised mortality ratios (SMRs) compared observed age- and sex-specific mortality rates of the cohort with age- and calendar-year-matched Australian population estimates using indirect standardisation. Five-year age groups were assigned. The 95 % confidence intervals (CI) were calculated using the Poisson distribution.
Semi-parametric analysis using Cox regression was used to describe baseline characteristics associated with mortality, including socio-demographic, drug use and health service use data. Validated tools included a variant of the AUDIT-C scale  to assess past month alcohol consumption (categorised as abstinent, 1–7, 8+) and the SF8 quality of life scale summary measurements for physical and mental health  as measures of health-related quality of life (categorised as low (lower quartile scores), average (second and third quartile scores) and high (upper quartile scores) ). The small number of events (deaths) restricted analyses to bivariate relationships between exposures and mortality. Variables were considered significant at p < 0.05. Associations were reported using hazard ratios (HR); 95 % confidence intervals (CIs) for crude mortality rates were calculated using a Poisson distribution. We describe time-varying factors associated with mortality using longitudinal data from the MIX cohort as supplementary analyses to identify missed opportunities for interventions [2, 7]. Analyses were conducted using Stata Version 13.1 (StataCorp LP, TX, USA).
We recruited 688 participants between 2008 and 2010, of whom 23 (3.4 %) were excluded from analysis due to incomplete data on exposures of interest such as past month use of emergency, outpatient and ambulance services. No mortality cases were excluded in this way. Participants with complete data (n = 665) contributed 2277 PYs across a median of three interviews. Table 1 describes baseline characteristics of participants who were young relative to other studies of Australian PWID (median age 29, interquartile range (IQR) 17–40), predominantly male and unemployed. Participants had been injecting for a median of 10 years (IQR 6–13) and over a third were receiving OST. Participants were frequent users of health services and almost one-third had poor self-reported physical health.
Mortality rates, causes of death, and correlates of mortality
There were 24 (3.6 %) deaths in study period; the median age at death was 28 years (range 19–31) and the median time from the last interview to death was 257 days (range 7–1134 days). The estimated mortality rate was 1.0 per 100 PY (95 % CI 0.7–1.6). The estimated SMR was 17.3 (95 % CI 11.6–25.8) overall and there was no evidence of a significant difference between women (23.5, 95 % CI 10.5–52.3) and men (15.9, 95 % CI 10.0–25.3).
Cause of death was available for 23 of the 24: 12 were drug-related, involving opioids and multiple drug use, four were suicides, four were due to physical assault, two were from medical conditions not directly related to injecting drug use, and one resulted from an accident-related injury. Heroin was the most common drug detected in the drug-related deaths (83 %), followed by benzodiazepines (67 %). Most (83 %) drug-related deaths involved multiple substances, including heroin, benzodiazepines, methadone, antidepressants, antipsychotics and alcohol.
Bivariate Cox regression analyses of associations between mortality and exposures are presented in Table 1. Increased mortality risk was associated with frequency of past incarceration (four or more compared to none), reporting past month ambulance attendance, reporting past month ED presentation and reporting a heroin overdose in the past 6 months. The supplementary analyses of time-varying exposures showed the same pattern, but some exposures failed to reach significance (Additional file 1: Table S1).
The excess mortality reported in this study was almost three times higher than previous Australian estimates involving PWID and people on opioid substitution therapy (OST) [2, 11]. The excess mortality was also higher than reported in other countries with a low HIV prevalence . The heightened risk of mortality in this community-based cohort might be due to systematic differences beyond injecting drug use. The cohort was community-recruited rather than treatment-recruited  and was relatively young compared to other cohorts [2, 11]. The excess mortality is comparable to results in contemporaneous cohorts of similar age , and the mortality rate might converge to estimates in the wider literature over a longer observation period potentially reflecting changes in risk behaviour with age and retention in drug treatment observed in other studies [11, 13].
The mortality rate reported here is comparable to that in a recent Australian cohort of PWID recruited from prison settings and followed post-release . Over half the MIX cohort reported a history of incarceration at baseline and we found that frequent incarceration was associated with mortality. Multiple events of incarceration indicate a particularly vulnerable group with an over-representation of low levels of education and socioeconomic status, and unsafe injecting practices . The primary overdose prevention tools potentially available through the criminal justice system are the pre- to post-release continuation of OST and the provision of naloxone at prison release. The impact of such programmes can be enhanced through a continuum of care between prison and community health services that facilitates programme retention.
Ambulance and ED attendance reported in the month prior to the last interview were associated with increased risk of mortality, as was reporting a heroin overdose in the previous 6 months. For three of the five deaths for which a previous ambulance attendance was reported, the attendance related to an overdose (data not shown). Frequent contact with emergency and ambulance services could be used to identify PWID at greater risk of fatal overdoses. Multiple overdose events may indicate reduced tolerance, interactions with other substances and the presence of other comorbidities . Contact with emergency services is an opportunity to provide targeted harm reduction initiatives, such as take-home naloxone.
Our study was limited by a small number of mortality events which affected the precision of estimates and limited adjusted mortality models. In most cases there was only a short observation period available for participants who died which limited our ability to observe the effect of time-varying risk factors, although this strengthens the temporal relevance of the observations in relation to mortality. The quality of the linkage was dependent on the accuracy of participant-provided identifier data. Exposures in cohort studies such as this are assumed to be fixed between interviews or censorship; however, this assumption is not always valid and means the predictors of mortality reported here should be interpreted with caution.
This analysis of mortality in community-recruited Australian PWID found significantly higher excess mortality than previous Australian studies and studies in comparable developed countries. These results highlight the importance of tailored interventions prioritising overdose prevention to young PWID. Our findings point to prison and emergency settings as ideal opportunities to deliver interventions such as harm reduction messages and referrals for naloxone training and prescribing to reduce mortality .
We describe mortality in a cohort of young, community-based Australian PWID
Excess mortality in the cohort was higher than other local estimates and settings where HIV prevalence is low
Recent use of emergency services such as emergency departments and ambulances was associated with mortality
Opportunities remain in prisons and emergency settings to deliver overdose-related interventions
Mathers BM, Degenhardt L, Phillips B, Wiessing L, Hickman M, Strathdee SA, et al. Global epidemiology of injecting drug use and HIV among people who inject drugs: a systematic review. Lancet. 2008;372(9651):1733–45.
Stoové MA, Dietze PM, Aitken CK, Jolley D. Mortality among injecting drug users in Melbourne: a 16-year follow-up of the Victorian injecting cohort study (VICS). Drug Alcohol Depend. 2008;96(3):281–5.
Mathers BM, Degenhardt L. Examining non-AIDS mortality among people who inject drugs. AIDS. 2014;28 Suppl 4:S435–44.
Stoove MA, Dietze PM, Jolley D. Overdose deaths following previous non-fatal heroin overdose: record linkage of ambulance attendance and death registry data. Drug Alcohol Rev. 2009;28(4):347–52.
Nambiar D, Stoove M, Dietze P. A cross-sectional study describing factors associated with utilisation of GP services by a cohort of people who inject drugs. BMC Health Serv Res. 2014;14:308.
Horyniak D, Higgs P, Jenkinson R, Degenhardt L, Stoove M, Kerr T, et al. Establishing the Melbourne injecting drug user cohort study (MIX): rationale, methods, and baseline and twelve-month follow-up results. Harm Reduction J. 2013;10(1):11.
Miller CL, Kerr T, Strathdee SA, Li K, Wood E. Factors associated with premature mortality among young injection drug users in Vancouver. Harm Reduction J. 2007;4:1.
Dietze P, Jenkinson R, Aitken C, Stoove M, Jolley D, Hickman M, et al. The relationship between alcohol use and injecting drug use: impacts on health, crime and wellbeing. Drug Alcohol Depend. 2013;128(1–2):111–5.
Ware JE. How to score and interpret single-item health status measures: a manual for users of the of the SF-8 health survey Lincoln, RI. Boston: QualityMetric, Inc.; Health Assessment Lab; 2001.
Nambiar D, Agius PA, Stoove M, Hickman M, Dietze P. Cessation of injecting drug use: the effects of health service utilisation, drug use and demographic factors. Drug Alcohol Depend. 2015;154:208–13.
Degenhardt L, Larney S, Randall D, Burns L, Hall W. Causes of death in a cohort treated for opioid dependence between 1985 and 2005. Addiction. 2014;109(1):90–9.
Degenhardt L, Bucello C, Mathers B, Briegleb C, Ali H, Hickman M, et al. Mortality among regular or dependent users of heroin and other opioids: a systematic review and meta-analysis of cohort studies. Addiction. 2011;106(1):32–51.
McDonald SA, Hutchinson SJ, Bird SM, Mills PR, Dillon J, Bloor M, et al. A population-based record linkage study of mortality in hepatitis C-diagnosed persons with or without HIV coinfection in Scotland. Stat Methods Med Res. 2009;18(3):271–83.
Degenhardt L, Larney S, Kimber J, Gisev N, Farrell M, Dobbins T, et al. The impact of opioid substitution therapy on mortality post-release from prison: retrospective data linkage study. Addiction. 2014;109(8):1306–17.
Galea S, Vlahov D. Social determinants and the health of drug users: socioeconomic status, homelessness, and incarceration. Public Health Rep. 2002;117 Suppl 1:S135–45.
Warner-Smith M, Darke S, Lynskey M, Hall W. Heroin overdose: causes and consequences. Addiction. 2001;96(8):1113–25.
Strang J, Bird SM, Dietze P, Gerra G, McLellan AT. Take-home emergency naloxone to prevent deaths from heroin overdose. BMJ. 2014;349:g6580.
The authors thank study participants and the research team who conduct data collection and management, and assist with interpretation. We also acknowledge the support received from the Australian Institute of Health and Welfare in conducting the linkage, as well as funding from the National Health and Medical Research Council (NHMRC) and the Colonial Foundation, the Centre for Research Excellence in Injecting Drug Use and the Victorian Operational Infrastructure Support Program.MS and PD receive Fellowship support from the NHMRC.
The authors declare that they have no competing interests.
DN, MS and PD developed the study concept; DN, MS and PA developed the analysis plan; DN conducted the analysis and drafted the manuscript; PA and MH contributed to the analysis; MH, MS and PD contributed to the interpretation of the results and revisions of the final manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Nambiar, D., Agius, P., Stoové, M. et al. Mortality in the Melbourne injecting drug user cohort study (MIX). Harm Reduct J 12, 55 (2015). https://doi.org/10.1186/s12954-015-0089-3