A Social Gradient in Fatal Opioids and Cocaine Related Overdoses?

Abstract

Background

To determine the existence of a social gradient in fatal overdose cases related to non-prescribed opioids and cocaine use, recorded in Luxembourg between 1994 and 2011.

Methods

Overdose cases were individually matched with four controls in a nested case-control study design, according to sex, year of birth, drug administration route and duration of drug use. The study sample, composed of 272 cases and 1,056 controls, was stratified according to a Social Inequality Accumulation Score (SIAS), based on educational attainment, employment, income, financial situation of subjects and the professional status of their father or legal guardian. Least squares linear regression analysis on overdose mortality rates and ridit scores were applied to determine the Relative Index of Inequality (RII) of the study sample.

Results

A negative linear relationship between the overdose mortality rate and the relative socioeconomic position was observed. We found a difference in mortality of 29.22 overdose deaths per 100 drug users in the lowest socioeconomic group compared to the most advantaged group. In terms of the Relative Inequality Index, the overdose mortality rate of opioid and cocaine users with lowest socioeconomic profiles was 9.88 times as high as that of their peers from the highest socioeconomic group (95% CI 6.49–13.26).

Conclusions

Our findings suggest the existence of a marked social gradient in opioids and cocaine related overdose fatalities. Harm reduction services should integrate socially supportive offers, not only because of their general aim of social (re)integration but crucially in order to meet their most important objective, that is to reduce drug-related mortality.

About the Authors:

Alain Origer

* E-mail: [email protected]

Affiliations Drug Coordination Office, Ministry of Health, Luxembourg, Luxembourg, Research Unit INSIDE, Institute Health & Behaviour, University of Luxembourg, Walferdange, Luxembourg

Etienne Le Bihan

Affiliation: Research Unit INSIDE, Institute Health & Behaviour, University of Luxembourg, Walferdange, Luxembourg

Michèle Baumann

Affiliation: Research Unit INSIDE, Institute Health & Behaviour, University of Luxembourg, Walferdange, Luxembourg

Introduction

Latest estimates from the United Nations Office on Drugs and Crime [1] and the World Health Organisation [2] refer to 210,000 and 245,000 drug-related deaths per year worldwide and approximately half of these cases are attributed to acute overdoses. Fatal overdoses caused by the use of illicit drugs in Europe account for an estimated 4% of premature death in young adults aged 15–39 years [3].

Social and economic conditions shape risk behaviors and the health of drug users as they impact on drug use patterns and affect the availability of resources and access to social welfare systems [4]. Nonetheless, fatal overdoses are avoidable to a large extend [5].

The existence of inequalities in health and mortality between groups with different socioeconomic profiles has been extensively documented since the 1980s [6–9].

Research focusing on social and economic determinants of drug-related mortality appears to be less prevalent, although previous studies have reported associations between fatal overdoses and individual parameters such as poor educational attainment [10–12], unemployment [13,14], low income [14,15], poverty status [15–17], homelessness [18,19] and poor psychosocial functioning [20].

Davoli et al. [21] investigated risk factors for drug overdose mortality, including 81 cases and 324 controls, matched on sex and year of birth and found no association with the educational status of victims. Galea et al. (2003) [22] applied a multi-level case-control study including 725 accidental overdose deaths and 453 accidental deaths due to other causes in various neighborhoods in New York City. Their findings showed increased odds of dying from overdose in neighborhoods with unequal income distribution as second level measure. In a further ecological study, including 59 residential community districts of New York City, Marzuk et al. [16] showed that the mortality rates due to overdoses involving cocaine and opiates are associated with poverty status. By comparing mortality rates according to causes and educational attainment in a cross-sectional study design, Borrell et al. [23] found highest relative risks of fatal drug overdose among victims with lowest educational profiles in the City of Barcelona. The study, however, included overdose deaths caused by both, illicit and prescription drugs.

The objective of our study is to investigating the existence of a social gradient in fatal drug overdose cases related to non-prescribed opioids’ and cocaine use by addressing the question whether the odds of a fatal overdose incident increase proportionally with the accumulation of social and economic disadvantages. Our findings should contribute to improving knowledge on the pathways that may lead to overdose incidents and to reduce their occurrence and consequences by implementing evidence-based public health, social and harm reduction policies.

Methods

In the present study, a fatal overdose case is defined as an intentional or unintentional death for which an acute adverse reaction after the recent use of products containing non- prescription opioids and/or cocaine was reported as the primary cause of death according to toxicological and forensic evidence. Concomitant use of other substances is not an exclusion criterion, if the use of opioids and/or cocaine is reported as the primary cause of death by forensic authorities.

Data sources

The following data sources have been linked and cross-examined in order to complete overdose victims’ profiles, drug use patterns and life histories:

Drug misuse surveillance data.

The national drug monitoring system (RELIS) is operated by the Luxembourg focal point of the European Monitoring Centre on Drugs and Drug Addiction (EMCDDA) and indexes drug-related contacts with both, drug demand and drug supply reduction institutions in a single and integrated database. The RELIS network includes national psychiatric departments of general hospitals, specialised in- and outpatient drug care services, including opioid substitution treatment offers and harm reduction services, prisons and special drug law enforcement agencies.

Indexed drug users are digitally anonymised by means of an encryption algorithm approved by the National Commission on Data Protection. The derived attribution code allows respondents’ tracing within single and between multiple data sources while guarantying anonymity of subjects, as requested by national data protection authorities. The RELIS data protocol includes routine items on socio-demographics, educational, socioeconomic and health status, drug use patterns and histories, treatment records and contacts with the penal system.

Forensic evidence.

In case of a suspicious death case, the public prosecutor’s office orders a toxicological investigation and an autopsy of the victim. The forensic department of the National Laboratory of Health (LNH) reports drug-related death cases also to the Ministry of Health for public health surveillance and statistical purposes. Forensic reports contain data on substances detected in the victims as well as an assessment of the association of detected substances and the occurrence of death. Autopsy reports also contain contextual information and elements of victims’ life history. Toxicological reports have been anonymised and made available to the research team by the LNH, following the authorisation of the Public Prosecutors Office.

National law enforcement records of fatal overdose cases.

National judicial police authorities exhaustively list overdose deaths occurring on the national territory after forensic confirmation. Toxicological evidence as well as related police reports are compiled for each case.

Linkage of data sources was possible as RELIS codes have been calculated and attributed to all overdose victims, previously de-identified. Since persons included in the national RELIS database or in the Police overdose record are anonymised by the same coding routine, matching cases could be reliably detected.

Matching procedures in the case-control setting

A nested case-control study design was applied. Cases are defined as victims of a fatal overdose having occurred in Luxembourg between 1994 and 2011. Matched controls refer to persons indexed by the epidemiological surveillance system on drug users in contact with national services (RELIS). Cases were matched with four controls for sex, year of birth, drug administration route and duration of drug use. Previous research has shown that matching more than four controls with each case does not significantly add statistical power to the analysis [24,25]. Controls were never matched to more than one case and deceased users were excluded from the RELIS database in order to avoid that they were matched with cases.

Matching variables are defined as follows:

* Sex (Male, female).

* Year of birth (In case no perfect match was available, a difference of one year more or less between the case and controls was accepted. Beyond age-specific risk differences, this variable also allowed to take account of potential changes in educational, economic and employment standards and contexts during the observation period).

* Route of administration (Injecting drug use, non-injecting drug use).

* Duration of illicit drug use (Matching according to duration of drug use is an essential condition to avoid the selection bias of selective survival. Risk exposure parameters of cases and controls were set to be comparable and in case no perfect match could be found, controls with the closest duration of drug use were chosen).

We used the Relative Index of Inequality (RII) to determine socioeconomic inequalities within the study sample, taking into account both, the size of sub-groups and the relative position of individuals within. In the present study, the RII measures the relative risk of a fatal overdose for the least advantaged group compared to the most advantaged group. The higher the RII, the more prevalent are socioeconomic inequalities within the study population.

The RII is a regression based index and only applies to linear changes in incidence (mortality) rates according to socioeconomic position parameters. To calculate the RII, we stratified our study sample according to a Social Inequality Accumulation Score (SIAS), based onfactors and values ^presented in Table 1. Scores attributed to each value have been weighted to ensure that the total weight of each factor is equal.

[Figure omitted. See PDF.]

Table 1. Factors and values with associated weighted Social Inequality Accumulation Scores (SIASs).

https://doi.org/10.1371/journal.pone.0125568.t001

SIASs have been attributed to cases and to controls to determine their respective socioeconomic position within the sample. The lowest possible SIAS equals to 0, whereas the most beneficial socioeconomic position is quantified by a SIAS of 5. Knowing that the study population is not comparable to the national general population in terms of socioeconomic parameters [12], mortality (by overdose) could not be assessed in absolute terms but based on the relative socioeconomic position that each individual occupies within the distribution of socio-economic parameters of the study sample, composed of cocaine and/or opioids’ users only.

In order to determine the RII, the sample has been stratified into 4 groups according to SIASs situated within the following intervals: ≤1; >1 and ≤ 2; >2 and ≤ 3; >3. We calculated respective overdose mortality rates as well as relative frequency (RF) and cumulative frequency (CF) of survivors, both required to determine a parameter called ridit (r) [26] for each of the 4 SIAS groups the study population (cases and controls) is composed of, according to the following formula:Weighted least squares linear regression analysis on overdose mortality rates and ridit scores were applied and the RII was calculated according to the formula provided by Kunst and Mackenbach [27]. Sex and age were not applied as adjustment variables to determine the RII since they were used as matching variables in the case-control selection process.

In order to account for missing values, we performed multiple imputations using the fully conditional specification approach [28] on the entire data set and generated 10 imputation data sets. Results have been combined across the imputed sets of data according to the method recommended by Rubin [29] and the overall RII estimate is the average of the individual RII estimates. Statistical analysis was conducted using SPSS, version 21.0.

Results

Our sample was composed of 1,328 persons including 272 fatal overdose cases and 1,056 controls. Six cases had to be excluded from the sample as no matching controls were found. For 3, 4 and 15 cases, respectively 1, 2 or 3 controls only could be matched. Eighty-one percent of controls matched cases for all matching variables, 4% were closest matches (+/- 1 year) for date of birth, 13% for duration of drug use career and 2% for both matching variables. The distribution of missing values for predictor variables is as follows: professional status of father or legal guardian (26.5%), debt status (16.8%), legal income, salary (13.7%), educational attainment (10.7%), and professional inactivity/unemployment (6.9%).

The overall characteristics of national fatal overdose victims and the case-control study sample were further described in previous publications [12,30].

Descriptive results (Table 2) reveal that FOD victims concentrate in lower SAIS classes. The mean SIAS of fatal overdose victims is 1.22, compared to 1.85 for survivors (t = 25.089; p = .000).

[Figure omitted. See PDF.]

Table 2. Mean SIAS and distribution within SIAS classes according to gender and fatal overdose (FOD) status.

https://doi.org/10.1371/journal.pone.0125568.t002

Table 3 presents overdose mortality rates according to SIAS classes as well as associated ridit scores.

[Figure omitted. See PDF.]

Table 3. Overdose mortality rates and ridit scores associated to SIAS classes.

https://doi.org/10.1371/journal.pone.0125568.t003

A negative linear relationship is observed between the overdose mortality rate and the relative socioeconomic position, determined by ridit scores, associated to SIAS (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Overdose Mortality Rate (OMR) according to relative socioeconomic position determined by ridit scores associated to SIAS.

Fig 1 shows that the overdose mortality rates decrease with increasing relative socioeconomic positions of opioids’ and/or cocaine users.

https://doi.org/10.1371/journal.pone.0125568.g001

We found an absolute difference in mortality of 29.22 overdose deaths per 100 drug users between the lowest SIAS class compared to the most advantaged class (Table 4). In terms of Relative Inequality Index, the overdose mortality rate of opioids’ and cocaine users with lowest socioeconomic profiles is 9.88 times as high [95% CI 6.49–13.26] as that of their peers from the highest socioeconomic group.

[Figure omitted. See PDF.]

Table 4. Results of least squares linear regression analysis on imputed data, imputed RIIs and final RII.

https://doi.org/10.1371/journal.pone.0125568.t004

In order to lay out the evolution of the RII and overdose mortality rates over time, a separate regression analysis was performed for two consecutive 9 years’ periods. Results showed a marked decrease in the RII from 17. 27 [95% CI 3.22–31.32] for 1994 to 2003 to 5.79 [95% CI 3.17–8.40] in reference to the period 2003 to 2011. Between the same time periods, a decrease in national overdose mortality from 6.37 to 4.22 cases/100.000 inhabitants aged 15 to 64 has been observed.

Discussion

Key results

This study investigated whether the odds of dying from opioid and/or cocaine related overdose increase gradually down the socioeconomic positions users of these drugs occupy. Results suggest the existence of a marked social gradient in opioids and cocaine related overdose fatalities, extending beyond the mere binary difference in mortality between the most and the least advantaged on the socioeconomic scale. In other words, adverse socioeconomic conditions affect the overdose mortality risk according to a ‘dose-response’ scheme.

In terms of prevention and harm reduction strategies, this is a relevant finding as it suggests that any measure aiming at reducing social disparities, even if isolated or targeted, may have a positive and dynamic impact that counter the negative cumulative effect in terms of survival. This strongly supports the idea that drug-related harm reduction strategies exclusively focusing on reducing direct harm of drug use fail to address harm associated with the context of drug use such as homelessness, violence and poverty [31]. Also, harm reduction services should integrate socially supportive offers, not only to reach their general goal of social reintegration but crucially, because they can contribute to save lives, generate years of life without major impairments and have a sustainable positive impact on drug-related mortality, and thus meet their first and overall objective: reducing drug-related morbidity and mortality.

Added value and generalizability

Studies on the association between single socioeconomic indicators and fatal overdoses in reference to the general population [14,32] and both, fatal and non-fatal overdoses in cross-sectional settings [33] have provided valuable outcomes in recent years. However, we could not identify previous case-control studies addressing inequalities in fatal overdose cases based on multiple socioeconomic parameters, combined into a cumulative measurement tool.

Moreover, the present study included probably for the first time, drug administration route and duration of drug use as matching variables in order to take due account of the risk exposition of both, cases and controls. These methodological aspects, as well as the large number of cases and controls and the cross-analysis of longitudinal data from different sources and settings, should add validity to our findings.

Limitations and further research

The scarcity of studies on socioeconomic determinants of fatal drug overdoses may partly be due to the difficulty to gather data on social profiles of drug users. Parameters such as education, income, financial wealth, professional activity and socioeconomic situation of the family of origin are often not included in routine drug registration or monitoring systems.

Also, in a retrospective study design, poor flexibility is given as to the selection of study variables. Additional data on income sources of subjects and the socioeconomic status of their parents could have contributed to a more refined stratification of subjects according to their respective socioeconomic position.

Furthermore, research on drug users mostly rely on episodic data, covering the period from first or established drug use to the moment of data collection. In a life-course perspective, it would be beneficial to have access to early life and life trajectory data in order to allow a more in-depth analysis of cumulative effects of socioeconomic determinants. Therefore, a prospective cohort study design would have been most appropriate to address the present research topic, bearing in mind, however, the resources’ and time consuming nature of such studies.

We found that the RII has markedly decreased during the observation period along with a parallel reduction of the average overdose mortality rate. Future research should, beyond investigating the existence of socioeconomic inequalities in drug using populations, further address the association between the amplitude and distribution of these inequalities and mortality risks.

This said, social and economic parameters should not be seen as independent determinants and even more so should the concept of social inequalities be widened when applied to drug-related mortality. Whereas the social status seems to hold a major role here; family, social and societal environments, new communication means, mobility, migration and acculturation contexts as well as the emergence of new psychoactive substances, new consume patterns and related risks are at play in the attempt to explain differences in terms of morbidity and mortality associated to drug use. The challenge for future research will thus lie in the capacity to take into account the changing context and environments that next generations will experience and the way they will, or will not, succeed to sharing resources and distributing wealth in a ‘healthy’ way.

Acknowledgments

Thanks are due to the Department of Toxicology of the National Laboratory of Health for providing forensic data and the staff of the EMCDDA Focal Point Luxembourg (Public Research Centre for Health) for its relentless efforts to ensure and maintain the quality of data provided by the national RELIS surveillance system and to all involved agencies.

Author Contributions

Conceived and designed the experiments: AO ELB MB. Performed the experiments: AO ELB MB. Analyzed the data: AO ELB MB. Contributed reagents/materials/analysis tools: AO ELB MB. Wrote the paper: AO ELB MB. Conceived and carried out the study, and had the main responsibility for the study and writing the manuscript: AO. Realized the statistical analysis and participated in data analysis: ELB. Participated in data analysis and in writing the manuscript: MB.

Citation: Origer A, Le Bihan E, Baumann M (2015) A Social Gradient in Fatal Opioids and Cocaine Related Overdoses? PLoS ONE 10(5): e0125568. https://doi.org/10.1371/journal.pone.0125568

References

1. UNODC (2013) World Drug Report 2013 (United Nations Publication), Vienna UNODC.

2. WHO (2009) Global Health Risks: Mortality and Burden of Disease Attributable to Selected Major Risks. Geneva: WHO.

3. EMCDDA (2011) Mortality related to drug use in Europe: public health implications. Lisbon: EMCDDA.

4. Galea S, Vlahov D (2002) Social determinants and the health of drug users: socioeconomic status, homelessness, and incarceration. Public Health Rep 117 Suppl 1: S135–145. pmid:12435837

5. Sporer KA (2003) Strategies for preventing heroin overdose. BMJ 326: 442–444. pmid:12595388

6. Marmot MG, McDowall ME (1986) Mortality decline and widening social inequalities. Lancet 2: 274–276. pmid:2874291

7. Wilkinson RG (1989) Class mortality differentials, income distribution and trends in poverty 1921–1981. J Soc Policy 18: 307–335. pmid:10296719

8. Illsley R (1990) Health inequities in Europe. Comparative review of sources, methodology and knowledge. Soc Sci Med 31: 229–236. pmid:2218603

9. Marmot MG, Kogevinas M, Elston MA (1987) Social/economic status and disease. Annu Rev Public Health 8: 111–135. pmid:3555518

10. Arias LC, Borrell C (1998) Mortality inequalities according to education in the City of Barcelona. Med Clin (Barc) 110: 161–166. pmid:9547718

11. Farrell M, Neeleman J, Griffiths P, Strang J (1996) Suicide and overdose among opiate addicts. Addiction 91: 321–323. pmid:8867195

12. Origer A, Le Bihan E, Baumann M (2014) Social and economic inequalities in fatal opioid and cocaine related overdoses in Luxembourg: A case-control study. Int J Drug Policy.

13. Harlow KC (1990) Patterns of rates of mortality from narcotics and cocaine overdose in Texas, 1976–87. Public Health Rep 105: 455–462. pmid:2120721

14. Pasarin M, Borrell C, Plasencia A (1999) Two patterns of social inequalities in mortality in Barcelona, Spain? Gac Sanit 13: 431–440. pmid:10619869

15. Torralba L, Brugal MT, Villalbi JR, Tortosa MT, Toribio A, Valverde JL (1996) Mortality due to acute adverse drug reactions: opiates and cocaine in Barcelona, 1989–93. Addiction 91: 419–426. pmid:8867204

16. Marzuk PM, Tardiff K, Leon AC, Hirsch CS, Stajic M, Portera L (1997) Poverty and fatal accidental drug overdoses of cocaine and opiates in New York City: an ecological study. Am J Drug Alcohol Abuse 23: 221–228. pmid:9143635

17. Jones R, Gruer L, Gilchrist G, Seymour A, Black M, Oliver J (2002) Recent contact with health and social services by drug misusers in Glasgow who died of a fatal overdose in 1999. Addiction 97: 1517–1522. pmid:12472635

18. Langendam MW, van Brussel GH, Coutinho RA, van Ameijden EJ (2001) The impact of harm-reduction-based methadone treatment on mortality among heroin users. Am J Public Health 91: 774–780. pmid:11344886

19. O'Driscoll PT, McGough J, Hagan H, Thiede H, Critchlow C, Alexander ER (2001) Predictors of accidental fatal drug overdose among a cohort of injection drug users. Am J Public Health 91: 984–987. pmid:11392946

20. Darke S, Zador D (1996) Fatal heroin 'overdose': a review. Addiction 91: 1765–1772. pmid:8997759

21. Davoli M, Perucci CA, Forastiere F, Doyle P, Rapiti E, Zaccarelli M (1993) Risk factors for overdose mortality: a case-control study within a cohort of intravenous drug users. Int J Epidemiol 22: 273–277. pmid:8505184

22. Galea S, Ahern J, Vlahov D, Coffin PO, Fuller C, Leon AC (2003) Income distribution and risk of fatal drug overdose in New York City neighborhoods. Drug Alcohol Depend 70: 139–148. pmid:12732407

23. Borrell C, Rodriguez M, Ferrando J, Brugal MT, Pasarin MI, Martínez V et al. (2002) Role of individual and contextual effects in injury mortality: new evidence from small area analysis. Inj Prev 8: 297–302. pmid:12460966

24. Ury HK (1975) Efficiency of case-control studies with multiple controls per case: continuous or dichotomous data. Biometrics 31: 643–649. pmid:1100136

25. Cai T, Zheng Y (2012) Evaluating prognostic accuracy of biomarkers in nested case-control studies. Biostatistics 13: 89–100. pmid:21856652

26. Bross ID (1958) "How to Use Ridit Analysis". Biometrics 14 (1):18–38

27. Mackenbach JP, Kunst AE (1997) Measuring the magnitude of socio-economic inequalities in health: an overview of available measures illustrated with two examples from Europe. Soc Sci Med 44: 757–771. pmid:9080560

28. van Buuren S (2007) Multiple imputation of discrete and continuous data by fully conditional specification. Stat Methods Med Res 16: 219–242. pmid:17621469

29. Rubin DB (1987) Multiple Imputation for Nonresponse in Surveys. New York: John Wiley & Sons. pmid:24122822

30. Origer A, Lopes da Costa S, Baumann M (2013) Opiate- and Cocaine-Related Fatal Overdoses in Luxembourg from 1985 to 2011: A Study on Gender Differences. Eur Addict Res 20: 87–93. pmid:24192492

31. Pauly B (2008) Harm reduction through a social justice lens. Int J Drug Policy 19: 4–10. pmid:18226520

32. Arias LC, Borrell C (1998) [Mortality inequalities according to education in the city of Barcelona]. Med Clin (Barc) 110: 161–166. pmid:9547718

33. Bergenstrom A, Quan VM, Van Nam L, McClausland K, Thuoc NP, Celentano D et al. (2008) A cross-sectional study on prevalence of non-fatal drug overdose and associated risk characteristics among out-of-treatment injecting drug users in North Vietnam. Subst Use Misuse 43: 73–84. pmid:18189206