Hostname: page-component-586b7cd67f-t8hqh Total loading time: 0 Render date: 2024-11-25T10:14:23.412Z Has data issue: false hasContentIssue false

The role of depression and use of alcohol and other drugs after partner suicide in the association between suicide bereavement and suicide: cohort study in the Danish population

Published online by Cambridge University Press:  11 March 2024

Alexandra Pitman*
Affiliation:
UCL Division of Psychiatry, 149 Tottenham Court Rd, London W1T 7AD, UK Camden and Islington NHS Foundation Trust, St Pancras Hospital, St Pancras Way, London NW1 0PE, UK
Keltie McDonald
Affiliation:
UCL Division of Psychiatry, 149 Tottenham Court Rd, London W1T 7AD, UK
Yanakan Logeswaran
Affiliation:
UCL Division of Psychiatry, 149 Tottenham Court Rd, London W1T 7AD, UK
Glyn Lewis
Affiliation:
UCL Division of Psychiatry, 149 Tottenham Court Rd, London W1T 7AD, UK Camden and Islington NHS Foundation Trust, St Pancras Hospital, St Pancras Way, London NW1 0PE, UK
Julie Cerel
Affiliation:
Suicide Prevention & Exposure Lab, College of Social Work, University of Kentucky, Lexington, USA
Gemma Lewis
Affiliation:
UCL Division of Psychiatry, 149 Tottenham Court Rd, London W1T 7AD, UK
Annette Erlangsen
Affiliation:
Danish Research Institute for Suicide Prevention – DRISP, Psychiatric Center Copenhagen, Copenhagen, Denmark Copenhagen Research Centre for Mental Health, Mental Health Center Copenhagen, Mental Health Services, Capital Region of Denmark, Copenhagen, Denmark Department of Mental Health, Johns Hopkins School of Public Health, Baltimore, USA Centre for Mental Health Research, The National Centre for Epidemiology and Population Health, ANU College of Health and Medicine, The Australian National University, Canberra, Australia
*
Corresponding author: Alexandra Pitman; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Background

Although suicide bereavement is associated with suicide and self-harm, evidence regarding mechanisms is lacking. We investigated whether depression and substance use (alcohol and/or other drugs) explain the association between partner suicide bereavement and suicide.

Methods

Linkage of nationwide, longitudinal data from Denmark for the period 1980–2016 facilitated a comparison of 22 668 individuals exposed to bereavement by a partner's suicide with 913 402 individuals bereaved by a partner's death due to other causes. Using causal mediation models, we estimated the degree to which depression and substance use (considered separately) mediated the association between suicide bereavement and suicide.

Results

Suicide-bereaved partners were found to have a higher risk of suicide (HRadj = 1.59, 95% CI 1.36–1.86) and of depression (ORadj 1.16, 95% CI 1.09–1.25) when compared to other-bereaved partners, but a lower risk of substance use (ORadj 0.83; 95% CI 0.78–0.88). An increased risk of suicide was found among any bereaved individuals with a depression diagnosis recorded post-bereavement (ORadj 3.92, 95% CI 3.55–4.34). Mediation analysis revealed that depression mediated 2% (1.68%; 95% CI 0.23%–3.14%; p = 0.024) of the association between suicide bereavement and suicide in partners when using bereaved controls.

Conclusions

Depression is a partial mediator of the association between suicide bereavement and suicide. Efforts to prevent and optimize the treatment of depression in suicide-bereaved people could reduce their suicide risk. Our findings might be conservative because we did not include cases of depression diagnosed in primary care. Further work is needed to understand this and other mediators.

Type
Original Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2024. Published by Cambridge University Press

Introduction

Consistent evidence demonstrates an association between suicide bereavement and subsequent risk of suicide (Pitman, Osborn, King, & Erlangsen, Reference Pitman, Osborn, King and Erlangsen2014, Reference Pitman, McDonald, Logeswaran, Lewis, Cerel and Erlangsen2022). Apart from preventing exposure to suicide loss in significant others, intervention hinges upon identification of modifiable risk factors. However, we still lack insights into potential explanatory factors between suicide bereavement and suicide. These include psychiatric disorders (particularly depression), alcohol or other drug use, stigma, grief, loneliness, shared social and economic adversity, suicide suggestion, and genetic factors (inherited or due to homophily) (O'Connor & Nock, Reference O'Connor and Nock2014; Pitman et al., Reference Pitman, Osborn, King and Erlangsen2014). It is important to understand these links to inform the content of post-suicide emotional support (postvention); an established policy adopted widely within suicide prevention strategies (Schlichthorst et al., Reference Schlichthorst, Reifels, Spittal, Clapperton, Scurrah, Kolves and Krysinska2022; World Health Organization, 2018). However, while evidence supports the effectiveness of postvention in reducing depression and anxiety there is no evidence it helps reduce risks of suicide (Andriessen et al., Reference Andriessen, Krysinska, Hill, Reifels, Robinson, Reavley and Pirkis2019; Linde, Treml, Steinig, Nagl, & Kersting, Reference Linde, Treml, Steinig, Nagl and Kersting2017; McDaid, Trowman, Golder, Hawton, & Sowden, Reference McDaid, Trowman, Golder, Hawton and Sowden2008; Szumilas & Kutcher, Reference Szumilas and Kutcher2011).

To date, the mechanisms underlying the elevated risk of suicide after suicide bereavement have not been investigated longitudinally to understand temporality. The two strongest candidate mediators are depression and substance (alcohol and/or other drugs) use, which are both prevalent after a suicide loss (Bolton et al., Reference Bolton, Au, Leslie, Martens, Enns, Roos and Sareen2013; Erlangsen et al., Reference Erlangsen, Runeson, Bolton, Wilcox, Forman, Krogh and Conwell2017; Pitman et al., Reference Pitman, Osborn, King and Erlangsen2014; Spiwak et al., Reference Spiwak, Elias, Sareen, Chartier, Katz and Bolton2020) and each is associated with suicide (Fazel & Runeson, Reference Fazel and Runeson2020). Qualitative studies in Australia, the UK, and the US provide accounts of depressive symptomatology after suicide bereavement (Entilli, Ross, de Leo, Cipolletta, & Kõlves, Reference Entilli, Ross, de Leo, Cipolletta and Kõlves2021; Pitman et al., Reference Pitman, Khrisna Putri, de Souza, Stevenson, King, Osborn and Morant2018) linked to suicidality (Hunt, Young, & Hertlein, Reference Hunt, Young and Hertlein2019), while longitudinal studies of suicide-bereaved partners and children confirm an increased probability of mood disorders (Appel et al., Reference Appel, Johansen, Deltour, Frederiksen, Hjalgrim, Dalton and Bidstrup2013; Erlangsen et al., Reference Erlangsen, Runeson, Bolton, Wilcox, Forman, Krogh and Conwell2017; Kuramoto et al., Reference Kuramoto, Stuart, Runeson, Lichtenstein, Långström and Wilcox2010; Pitman et al., Reference Pitman, Osborn, King and Erlangsen2014; Spiwak et al., Reference Spiwak, Elias, Sareen, Chartier, Katz and Bolton2020; Wilcox et al., Reference Wilcox, Kuramoto, Lichtenstein, Långström, Brent and Runeson2010) and antidepressant treatment (Appel et al., Reference Appel, Johansen, Christensen, Frederiksen, Hjalgrim, Dalton and Bidstrup2016) after suicide loss. Qualitative research in the US and UK also documents excessive use of alcohol and other drugs after being bereaved by suicide (Eng et al., Reference Eng, Drabwell, Stevenson, King, Osborn and Pitman2019; Hunt et al., Reference Hunt, Young and Hertlein2019; Pitman, Stevenson, King, & Osborn, Reference Pitman, Stevenson, King and Osborn2020) and other sudden unnatural causes (Drabwell et al., Reference Drabwell, Eng, Stevenson, King, Osborn and Pitman2020; Pitman et al., Reference Pitman, Stevenson, King and Osborn2020) as a means of coping with overwhelming thoughts and emotions. However, some individuals described reducing their use of alcohol or other drugs to better cope with emotions or to avoid substances perceived as contributing to the death of their relative or friend (Drabwell et al., Reference Drabwell, Eng, Stevenson, King, Osborn and Pitman2020; Eng et al., Reference Eng, Drabwell, Stevenson, King, Osborn and Pitman2019; Pitman et al., Reference Pitman, Stevenson, King and Osborn2020). Longitudinal studies of suicide-bereaved offspring report an increased probability of alcohol or other drug use after suicide bereavement (Brent, Melhem, Donohoe, & Walker, Reference Brent, Melhem, Donohoe and Walker2009; Wilcox et al., Reference Wilcox, Kuramoto, Lichtenstein, Långström, Brent and Runeson2010), possibly transient (Hamdan, Melhem, Porta, Song, & Brent, Reference Hamdan, Melhem, Porta, Song and Brent2013), but longitudinal studies of suicide-bereaved partners find a similar prevalence of alcohol or other drug use disorders to partners bereaved by other causes (Bolton et al., Reference Bolton, Au, Leslie, Martens, Enns, Roos and Sareen2013; Erlangsen et al., Reference Erlangsen, Runeson, Bolton, Wilcox, Forman, Krogh and Conwell2017; Spiwak et al., Reference Spiwak, Elias, Sareen, Chartier, Katz and Bolton2020). Such inconsistencies may be explained by differences in measurement and samples. It is also possible that mechanisms differ by kinship group (Pitman et al., Reference Pitman, McDonald, Logeswaran, Lewis, Cerel and Erlangsen2022).

Given the evidence supporting depression and alcohol and/or other drug use as candidate mediators, large longitudinal studies are required to investigate this, ideally by comparing suicide bereavement to other bereavements to account for the experience of bereavement per se. New methods of counterfactual mediation analysis (MacKinnon, Fairchild, & Fritz, Reference MacKinnon, Fairchild and Fritz2007) are underused in suicide research. These offer a more robust means of assessing the relative magnitude of different pathways and mechanisms by which an exposure may affect an outcome because it makes the causal assumptions explicit, adjusts for confounders of each arm of mediational models, and takes into account interactions with exposure (VanderWeele, Reference VanderWeele2016a). The aim of this study was to estimate the relative proportions of the association between suicide bereavement and suicide mediated by depression and by substance use (alcohol and/or other drugs) when comparing suicide-bereaved partners to other-bereaved partners.

Methods

Study design and participants

We conducted a population-based cohort study, analyzing data from administrative national registers linked using a unique personal identification number assigned to all individuals in Denmark (Erlangsen & Fedyszyn, Reference Erlangsen and Fedyszyn2015). Population data from the Civil Registration System was linked with information on psychiatric and somatic hospital contacts from the Psychiatric Central Research Register (since 1970) (Mors, Perto, & Mortensen, Reference Mors, Perto and Mortensen2011) and the National Patient Register (since 1977) (Schmidt et al., Reference Schmidt, Schmidt, Sandegaard, Ehrenstein, Pedersen and Sørensen2015) but linkage was not available to primary care data.

We included all Danish-born individuals aged 16 years and older who were living in Denmark between 1 January 1980 and 31 December 2016, comparing participants bereaved by the (a) suicide or (b) non-suicide death of a current/former partner (spouses, civil partners, and cohabitees) over that period. Cohort inception in 1980 provided us with over 10 years of previous data on psychiatric confounders as identified in the Psychiatric Central Research Register.

We focused on studying mechanisms in one kinship group (partners) because they offered the largest sample relative to parents, offspring, and siblings (Pitman et al., Reference Pitman, McDonald, Logeswaran, Lewis, Cerel and Erlangsen2022). Exposure to partner bereavement was identified by linking data on individuals in the cohort to deceased current and former spouses/cohabitees, identifying causes of death from the Register of Causes of Death. Partners (opposite and same-sex married and registered partners and opposite sex cohabitees) were identified using an established linkage method based on data recorded in the Civil Registration System (Pedersen, Reference Pedersen2011) (online Supplemental Methods S1). Individuals entered the cohort on the date of the index (first) bereavement and were followed up until the date of their suicide, death from a cause other than suicide, migration out of Denmark, second partner bereavement, or the end of follow-up (31 December 2016).

Measures

Outcome

Our outcome, suicide, was identified from the Register of Causes of Death (Helweg-Larsen, Reference Helweg-Larsen2011) based on International Classification of Diseases (ICD)-8 and ICD-10 codes (online Supplemental Table S1).

Exposure

Relevant ICD-8 and ICD-10 codes were used to identify whether bereavement was due to suicide or other causes. The first partner bereavement over this period was identified as the index exposure, and individuals were censored at any subsequent partner bereavement.

Confounders

We chose eight confounders a priori, based on existing evidence (Pitman et al., Reference Pitman, Osborn, King and Erlangsen2014), which we used in each arm of the mediation model: sex; age; bereavement year; marital status (as distinct from cohabitation status, to capture the effect of divorce); household income level (quartiles); pre-bereavement history of any admission recording self-harm (in psychiatric or physical health settings); pre-bereavement psychiatric disorders recorded on any admission (including alcohol and/or other drug use); and pre-bereavement physical health conditions recorded on medical admission (online Supplemental Table S1). Bereavement year captured period effects and differing inception years for inpatient and outpatient data. All confounders were measured before the index bereavement.

Mediators

Depression was defined as a diagnosis recorded during psychiatric inpatient (since 1980) or outpatient (since 1995) contact. Substance use was defined broadly as an ICD-8/ICD-10 diagnosis of alcohol and/or other drug use disorders recorded during psychiatric inpatient (since 1980) or outpatient (since 1995) contact, including acute intoxication, as well as inpatient medical admission recording specific medical problems linked to alcohol use (Grissa, Rasmussen, Krag, Brunak, & Jensen, Reference Grissa, Rasmussen, Krag, Brunak and Jensen2020). Both mediators were restricted to events recorded after the index bereavement.

Statistical analysis

The demographic and clinical characteristics of the sample were examined using complete case analysis. We used Cox proportional hazards regression to examine the associations between suicide bereavement and time to suicide in unadjusted and adjusted models. We assessed model assumptions including proportionality of hazards (online Supplemental Methods S2). In preliminary analyses we accounted for potential clustering effects, where multiple current/former partners might be bereaved by the same death, using multilevel Cox proportional hazards regression clustered on the individual who died. As the multilevel model provided comparable estimates to the single-level model, we opted for the latter.

Mediation by depression and substance use was assessed separately using the causal inference potential outcomes framework to decompose the total effect (TE) of suicide bereavement on suicide into four components: (a) controlled direct effect (CDE), (b) reference interaction (INTref), (c) mediated interaction (INTmed); and (d) pure indirect effect (PIE), and other derived values, including the total indirect effect (TIE), portion attributable to interaction (PAI), and the portion eliminated (PE) (Richiardi, Bellocco, & Zugna, Reference Richiardi, Bellocco and Zugna2013; Vanderweele, Reference Vanderweele2014; VanderWeele, Reference VanderWeele2016b) (online Supplemental Methods S4; Supplemental Box S1). From this output, we identified: the overall proportion of the association mediated (equivalent to 100 × [TIE/TE]); the proportion of the association attributable to interaction (effectively 100 × [PAI/TE]); and the proportion of the total effect of the exposure on the outcome that would be eliminated in the absence of the mediator (effectively 100 × PE/TE). The last of these is of policy relevance in conveying how much of the effect of the exposure can be prevented by intervening on the mediator.

We fitted the mediation model using the med4way command in Stata (Discacciati, Bellavia, Lee, Mazumdar, & Valeri, Reference Discacciati, Bellavia, Lee, Mazumdar and Valeri2019; StataCorp, 2019). The outcome model (path C in Fig. 1) was fitted using Cox proportional hazards regression, while the mediation model (paths A and B) was fitted using logistic regression because med4way does not accommodate a time-to-event mediator and time-to-event outcome.

Figure 1. Hypothesized model for the mediating effect of (A) depression and (B) substance use on the association between suicide bereavement and suicide.

Legend: Panel A displays the hypothetical model for mediation by depression. Panel B displays the hypothetical model for mediation by substance use. Solid lines represent potential causal pathways. Dashed lines reflect potential confounding pathways. Path A represents the association between suicide bereavement and post-bereavement depression/substance use (modeled using logistic regression). Path B represents the association between post-bereavement depression/substance use (in all those bereaved) and suicide modeled using logistic regression. Path C represents the association between suicide bereavement and suicide (modeled using Cox proportional hazards regression). The same set of confounders was used for all pathways: sex, age, bereavement year, legal marital status, household income level, pre-bereavement history of self-harm, any psychiatric disorders, and any physical disorders.

In our first sensitivity analysis we assessed the influence of competing risks (death from other causes; second bereavement; emigration) in our main association (Supplemental Methods 3). We then conducted an interaction test to investigate whether partner status (ex- versus current) modified this main association. In our second sensitivity analysis we repeated our mediation analyses restricted to records of inpatient admissions to ascertain the influence of having incorporated outpatient contact data within our main analyses. In our third sensitivity analysis we repeated our mediation analyses excluding those with widowed marital status to ascertain the influence of having included those with additional partner bereavements prior to our period of interest.

We repeated our mediation analyses restricted to records of inpatient admissions to ascertain the influence of having incorporated outpatient contact data within our main analyses. We also repeated our mediation analyses excluding those with widowed marital status to ascertain the influence of having included those with additional partner bereavements prior to our period of interest.

Data management was performed using SAS software 9.4 (SAS System for SunOP, 2003) and analyses were performed using Stata 17 software (StataCorp, 2021).

Results

Sample characteristics

A total of 960 272 bereaved partners were identified of whom 24 202 (2.5%) were excluded due to missing data on income level (2.5%) and/or marital status (0.01%) (online Supplemental Table S2). Among the remaining 936 070 (97.5%) partners, 22 668 (2.4%) had been bereaved by suicide and 913 402 (97.6%) by other causes of death (Table 1). The majority of the analytic sample was female (67.3%). Individuals bereaved by a partner's suicide were more likely to be younger, female, to have a higher household income level and a past history of self-harm.

Table 1. Demographic and clinical features of the studied cohort (n = 936 070) according to bereavement status

IQR, interquartile range; PTSD, post-traumatic stress disorder; COPD, chronic obstructive pulmonary disease; SMI, severe mental illness (defined as psychotic disorders, manic episode, bipolar affective disorder, and depression with psychotic symptoms).

a Data are n (%), except age, which is summarized as median (IQR). All values were pre-bereavement unless otherwise specified.

b Household income quartiles represent total income within the household divided by the total number of adults living in the household, then categorized into quartiles based on national annual income averages.

c The widowed category represented people who were bereaved by subsequent partner loss before the population registers had started, yet were bereaved while in a new partnership, so by default were registered as widowed by a former partner's death.

d Data for these variables were from diagnoses (or self-harm) recorded on inpatient admissions (i.e. excluding outpatient data).

e Data for these variables were from diagnoses recorded on inpatient admissions and outpatient contacts.

All p-values were <0.001 apart from depression (<0.01) and PTSD (non-significant) but these values are to be interpreted in the context of a large sample size, and therefore a high probability of statistically significant differences.

Main association

During almost 10 million person-years of follow-up, 3733 individuals died by suicide, of whom 182 (4.9%) had been bereaved by partner's suicide and 3551 (95.1%) bereaved by other causes. Individuals bereaved by a partner's suicide had a higher risk of subsequent suicide compared to those bereaved by other causes (HRadj: 1.59, 95% CI 1.36–1.86) (Table 2). The absolute risk of suicide after partner bereavement was 0.80% in individuals bereaved by a partner's suicide and 0.39% in individuals bereaved by other causes.

Table 2. Hazard ratios and 95% confidence intervals for the association between suicide bereavement and suicide compared with other bereavement (Fig. 1, Path C)

HR, hazard ratio; CI, confidence interval.

a Adjusted for sex, age, bereavement year, legal marital status, and household income level.

b Final model: adjusted for all variables in adjustment 1, plus pre-bereavement history of self-harm, any psychiatric disorders (one binary variable), and any physical disorders (one binary variable).

Mediation pathways

Suicide bereavement was associated with an increased risk of depression (Path A [exposure-mediator]; ORadj 1.16, 95% CI 1.09–1.25). Post-bereavement depression was associated with an increased risk of suicide (Path B [mediator-outcome]: ORadj 3.92, 95% CI 3.55–4.34; Table 3).

Table 3. Odds ratios and 95% confidence intervals for each pathway (modeled separately) in the model for mediation of the association between suicide bereavement and suicide by depression (Fig. 1, Panel A) and substance use (Fig. 1, Panel B)

OR, odds ratio; CI, confidence interval.

a Adjusted for sex, age, bereavement year, marital status, and household income level.

b Final model: adjusted for all variables in adjustment 1, plus pre-bereavement history of self-harm, any psychiatric disorders (one binary variable), and any physical disorders (one binary variable).

c N is solely for the outcome in each path; thus of those who are suicide-bereaved, 938 (4.1%) have depression and 1407 (6.2%) have substance use, and of those who are other-bereaved, 26744 (2.9%) have depression and 27597 (3.0%) have substance use. Of those who have depression, 510 (1.8%) die by suicide, and of those who do not have depression, 3223 (0.4%) die by suicide. Similarly, of those who have substance use, 383 (1.3%) die by suicide, and of those who do not have substance use, 3350 (0.4%) die by suicide.

* p < 0.001.

Suicide bereavement was associated with an increased risk of substance use in the unadjusted model (ORcrude 2.12; 95% CI 2.01–2.25; Table 3) but this attenuated and changed direction when adjusted (ORadj 0.83; 95% CI 0.78–0.88). Post hoc adjustments (online Supplemental Table S3) showed that the attenuation was primarily due to the confounding effect of age, as the suicide-bereaved group was younger than controls, and older people had a lower probability of substance use. Nevertheless, post-bereavement substance use was associated with an increased risk of suicide (ORcrude 3.61; 95% CI 3.25–4.02; ORadj 1.59, 95% CI 1.41–1.79).

Mediation analyses

We estimated that the proportion of the association between suicide bereavement and suicide attributable to mediation by depression (100 × TIE/TE) was 1.68% (95% CI 0.23–3.14%; p = 0.024; Table 4), supporting depression as a partial mediator. We found no evidence to support an interaction between suicide bereavement and depression (100 × PAI/TE = 3.05%; 95% CI −4.47 to 10.59%; p = 0.427; online Supplemental Results S1). The proportion of suicide that would be prevented among the suicide-bereaved if intervening support could prevent depression (100 × PE/TE) was 4.35% (95% CI −3.29 to 12.00%; p = 0.264).

Table 4. Mediation analyses describing the role of depression in the association between suicide bereavement and suicide

RER, relative excess risk; CI, confidence interval; see online Supplemental Box S1 for definitions.

a Adjusted for sex, age, bereavement year, marital status, household income level, pre-bereavement history of self-harm, any psychiatric disorders (one binary variable), and any physical disorders (one binary variable). Note that in a mediation analysis, only adjusted estimates are presented because mediation can only be interpreted in the context of adjusting for all measured confounders (see online Supplemental Methods S4).

b Where no stata label is given, this value was computed from other stata output values.

In a separate model for substance use, we found no evidence to support a mediating role in the association between suicide bereavement and suicide (Table 5).

Table 5. Mediation analyses describing the role of substance use in the association between suicide bereavement and suicide

RER, relative excess risk; CI, confidence interval; see online Supplemental Box S1 for definitions.

a Adjusted for sex, age, bereavement year, marital status, household income level, pre-bereavement history of self-harm, any psychiatric disorders (one binary variable), and any physical disorders (one binary variable). Note that in a mediation analysis, only adjusted estimates are presented because mediation can only be interpreted in the context of adjusting for all measured confounders (see online Supplemental Methods S4).

b Where no stata label is given, this value was computed from other stata output values.

Sensitivity analyses

When accounting for competing risks in our main association (online Supplemental Table S4), similar results were obtained to those from the main analysis (crude sub-distribution hazard ratio [sHR]: 1.99, 95% CI 1.71–2.31; sHRadj: 1.54, 95% CI 1.32–1.80; online Supplemental Table S5). No interaction was found with respect to current/ex-partner status (p = 0.130; online Supplemental Table S6). When we restricted depression to inpatient records only for our main association, suicide bereavement was not associated with an increased risk of depression or substance use (online Supplemental Table S7), therefore we did not conduct mediation analyses using either measure. When we ran our mediation model for depression excluding those with widowed marital status, the proportion mediated rose from 1.68% to 1.74%, but with overlapping confidence intervals.

Discussion

Main findings

Based on complete, nationwide data on all bereaved partners, we found evidence of an elevated risk of depression and of suicide after partner suicide bereavement compared with bereaved controls, but a reduced risk of substance use. Suicide-bereaved partners had almost twice the risk of suicide when compared to partners bereaved by other causes. Bereaved partners with post-bereavement depression had almost a four-fold risk of suicide. Depression mediated at least 2% of the association between suicide bereavement and suicide, but for a number of reasons this could be an underestimate. Substance use was not found to mediate the association. Our definition of depression omitted any diagnoses made in primary care, where the majority of cases of depression in Denmark are treated (Musliner et al., Reference Musliner, Liu, Gasse, Christensen, Wimberley and Munk-Olsen2019), cases diagnosed pre-1995 in outpatient care, and people never treated. Perceived stigma related to suicide might dissuade suicide-bereaved partners from seeking help (Hanschmidt, Lehnig, Riedel-Heller, & Kersting, Reference Hanschmidt, Lehnig, Riedel-Heller and Kersting2016; Yang, Wong, Grivel, & Hasin, Reference Yang, Wong, Grivel and Hasin2017), as could the perception that the system failed the deceased, resulting in under-ascertainment of psychiatric disorder and substance use when compared to controls. Finally, interventions offered during outpatient and inpatient care may also have mitigated suicide risk. For these reasons, it is plausible that depression mediates a greater proportion of the association between suicide bereavement and suicide than our study using secondary care data suggests. Further work is needed to investigate mediators of this association using both primary and secondary care data, as well self-report variables (where available) to capture perceived stigma and reluctance to seek help.

Findings in the context of other studies

This is, to our knowledge, the first formal evaluation of potential mediators of suicide risk after suicide loss. Our findings extend existing evidence regarding the elevated risk of depression (Erlangsen et al., Reference Erlangsen, Runeson, Bolton, Wilcox, Forman, Krogh and Conwell2017) and suicide (Agerbo, Reference Agerbo2003, Reference Agerbo2005; Erlangsen et al., Reference Erlangsen, Runeson, Bolton, Wilcox, Forman, Krogh and Conwell2017) among suicide-bereaved partners compared with bereaved controls. In contrast, our novel finding of a reduced risk of substance use in the suicide-bereaved compared with bereaved controls (primarily due to the confounding effect of age) differs from previous findings of no differences (Bolton et al., Reference Bolton, Au, Leslie, Martens, Enns, Roos and Sareen2013; Erlangsen et al., Reference Erlangsen, Runeson, Bolton, Wilcox, Forman, Krogh and Conwell2017). It is, however, consistent with the reduced risk of liver cirrhosis documented among suicide-bereaved spouses compared with bereaved controls (Erlangsen et al., Reference Erlangsen, Runeson, Bolton, Wilcox, Forman, Krogh and Conwell2017) and the accounts of people bereaved by suicide and other causes reporting a post-bereavement reduction in substance use (Drabwell et al., Reference Drabwell, Eng, Stevenson, King, Osborn and Pitman2020; Eng et al., Reference Eng, Drabwell, Stevenson, King, Osborn and Pitman2019; Pitman et al., Reference Pitman, Stevenson, King and Osborn2020). This may be linked to the enhanced awareness described by suicide-bereaved individuals of their shared vulnerabilities to suicide, engendering a determination to safeguard their mental health after the loss (Pitman et al., Reference Pitman, Nesse, Morant, Azorina, Stevenson, King and Osborn2017). We also noted that Danish partners bereaved by suicide had a significantly higher household income than partners bereaved by non-suicide causes, in keeping with a similar finding for bereaved partners in Canada (Spiwak et al., Reference Spiwak, Elias, Sareen, Chartier, Katz and Bolton2020) but in contrast with the converse finding for bereaved parents in Canada (Bolton et al., Reference Bolton, Au, Leslie, Martens, Enns, Roos and Sareen2013).

Strengths and limitations

Strengths include analyzing longitudinal data on a large population-based cohort of almost a million people, followed for approximately 10 million person-years, with minimal loss to follow-up, low levels of missing data, and avoidance of selection biases or recall biases. The Danish registers have been evaluated as reliable with respect to psychiatric hospital contacts (Tøllefsen et al., Reference Tøllefsen, Helweg-Larsen, Thiblin, Hem, Kastrup, Nyberg and Ekeberg2015) and registration of suicide death (Helweg-Larsen, Reference Helweg-Larsen2011; Tøllefsen et al., Reference Tøllefsen, Helweg-Larsen, Thiblin, Hem, Kastrup, Nyberg and Ekeberg2015) and unnatural deaths (Tøllefsen et al., Reference Tøllefsen, Helweg-Larsen, Thiblin, Hem, Kastrup, Nyberg and Ekeberg2015), with good validity of diagnostic codes and the timing of admission/discharge dates in the Danish National Patient Register (Schmidt et al., Reference Schmidt, Schmidt, Sandegaard, Ehrenstein, Pedersen and Sørensen2015) and of depression diagnoses and admission/outpatient dates in the Psychiatric Central Research Register (Mors et al., Reference Mors, Perto and Mortensen2011). For our study, this was important in separating out mediators and confounders for each pathway. In utilizing the Danish Civil Registration household variable, we were able to identify cohabiting couples as well as those in legal unions (although not opposite-sex cohabiting couples), providing a more realistic representation of Danish family structure and a more comprehensive measure of exposure to partner suicide loss. The limitations of using routine Danish registry data include the under-ascertainment of mental disorder and self-harm (due to lack of secondary care outpatient data pre-1995, primary care data, and data from people who do not seek or receive treatment) and not capturing other influences on suicide risk, such as personality traits, coping style, or social support.

Our application of contemporary causal inference methods overcame the limitations of older methods of mediation by allowing for the adjustment of confounders in each arm of the mediation model, making causal assumptions explicit, allowing us to model interactions, and was particularly suitable for binary outcomes(Rijnhart, Valente, MacKinnon, Twisk, & Heymans, Reference Rijnhart, Valente, MacKinnon, Twisk and Heymans2021). However, mediation assumes the absence of residual confounding between exposure-outcome, exposure-mediator, or mediator-outcome, and we acknowledge the possibility of residual confounding due to unmeasured pre-bereavement depression or self-harm (not captured in settings beyond psychiatric admission). Our mediation models could not assess the role of the other respective putative mediator as an intermediate confounder because med4way can only consider one mediator at a time, for ease of interpretability. Our use of this model also relied on a set of theoretical assumptions, and it is possible that not all were satisfied, introducing the potential for biased estimates of causal effects.

We also acknowledge the potential for collider bias where selecting bereaved individuals for comparison might result in distorted associations between variables (Holmberg & Andersen, Reference Holmberg and Andersen2022). Given the difficulties in aligning follow-up for a cohort comprising individuals bereaved in childhood, early adulthood, mid-life, and late life, and because mechanisms may differ by kinship, we opted only to assess bereaved partners. This diminished statistical power to detect mediation effects. As a putative mediator, our broad substance use definition (including medical complications) relied less on a clinician recording an alcohol and/or other drug use disorder diagnosis e during admission or on treatment-seeking for perceived substance use problems. However, our narrow definition of depression omitted untreated depression, pre-1995 cases of depression in outpatients, cases identified in primary care (as we lacked primary care) or during somatic admissions, or information from prescribing data (available from 1995). Ascertainment bias is also possible if suicide-bereaved individuals are more likely to be in contact with mental health professionals and therefore more likely to be admitted.

Clinical, policy, and research implications

The increased risk of depression after suicide loss and its contribution to suicide risk suggest a benefit of early identification of depression in suicide-bereaved partners (both ex- and current). Suicide risk mitigation might be achieved through appropriate bereavement support (Andriessen et al., Reference Andriessen, Krysinska, Hill, Reifels, Robinson, Reavley and Pirkis2019; Linde et al., Reference Linde, Treml, Steinig, Nagl and Kersting2017; McDaid et al., Reference McDaid, Trowman, Golder, Hawton and Sowden2008; Szumilas & Kutcher, Reference Szumilas and Kutcher2011), improving uptake through collaboration between bereavement services and primary care. Future research should use primary care and prescribing data to capture psychiatric disorders treated in primary and secondary care, investigating depression, substance use, and other putative mediators to advance our understanding of therapeutic targets. This will inform recommendations regarding appropriate support after suicide loss. The inclusion of postvention within a national suicide prevention strategy is not associated with changes in suicide mortality(Schlichthorst et al., Reference Schlichthorst, Reifels, Spittal, Clapperton, Scurrah, Kolves and Krysinska2022). However, a better understanding of mediators has the potential to improve the effectiveness of postvention and reduce suicides.

Conclusions

Suicide risk in suicide-bereaved partners was almost twice that for partners bereaved by other causes. Depression mediated a small proportion (2%) of this association, but this study likely under-estimated the mediating role of depression by omitting depression treated in primary care. Risk of substance use was reduced after partner suicide loss, and substance use did not mediate the association between suicide bereavement and suicide. These findings, from a large, longitudinal representative dataset, identified early treatment of depression as a potential means of reducing the burden of suicide after suicide bereavement.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291724000448

Funding statement

This work was supported by a grant from the American Foundation for Suicide Prevention (SRG-0–111-17). AP, Glyn Lewis and Gemma Lewis are also supported by the National Institute for Health Research (NIHR) University College London Hospital (UCLH) Biomedical Research Centre (BRC). Gemma Lewis is supported by a Sir Henry Dale Fellowship jointly funded by the Wellcome Trust and the Royal Society (grant number 223248/Z/21/Z). The funders were not involved in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Competing interests

None.

Ethical standards

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. The study was approved by the Danish Data Protection Agency (Capital Region of Denmark; P-2020-774) and the UCL (University College London) Research Ethics Committee (approval number 14 075/001). Informed consent was waived.

References

Agerbo, E. (2003). Risk of suicide and spouse's psychiatric illness or suicide: Nested case-control study. British Medical Journal, 327(7422), 10251026. https://doi.org/10.1136/bmj.327.7422.1025CrossRefGoogle ScholarPubMed
Agerbo, E. (2005). Midlife suicide risk, partner's psychiatric illness, spouse and child bereavement by suicide or other modes of death: A gender specific study. Journal of Epidemiology and Community Health, 59(5), 407412. https://doi.org/10.1136/jech.2004.024950CrossRefGoogle ScholarPubMed
Andriessen, K., Krysinska, K., Hill, N. T. M., Reifels, L., Robinson, J., Reavley, N., & Pirkis, J. (2019). Effectiveness of interventions for people bereaved through suicide: A systematic review of controlled studies of grief, psychosocial and suicide-related outcomes. BMC Psychiatry, 19(1), 49. https://doi.org/10.1186/s12888-019-2020-zCrossRefGoogle ScholarPubMed
Appel, C. W., Johansen, C., Christensen, J., Frederiksen, K., Hjalgrim, H., Dalton, S. O., … Bidstrup, P. E. (2016). Risk of use of antidepressants among children and young adults exposed to the death of a parent. Epidemiology (Cambridge, Mass.), 27(4), 578585. https://doi.org/10.1097/EDE.0000000000000481CrossRefGoogle Scholar
Appel, C. W., Johansen, C., Deltour, I., Frederiksen, K., Hjalgrim, H., Dalton, S. O., … Bidstrup, P. E. (2013). Early parental death and risk of hospitalization for affective disorder in adulthood. Epidemiology (Cambridge, Mass.), 24(4), 608615. https://doi.org/10.1097/EDE.0b013e3182915df8CrossRefGoogle ScholarPubMed
Bolton, J. M., Au, W., Leslie, W. D., Martens, P. J., Enns, M. W., Roos, L. L., … Sareen, J. (2013). Parents bereaved by offspring suicide: A population-based longitudinal case-control study. JAMA Psychiatry, 70(2), 158167. https://doi.org/10.1001/jamapsychiatry.2013.275CrossRefGoogle ScholarPubMed
Brent, D., Melhem, N., Donohoe, M. B., & Walker, M. (2009). The incidence and course of depression in bereaved youth 21 months after the loss of a parent to suicide, accident, or sudden natural death. American Journal of Psychiatry, 166(7), 786794. https://doi.org/10.1176/appi.ajp.2009.08081244CrossRefGoogle ScholarPubMed
Discacciati, A., Bellavia, A., Lee, J. J., Mazumdar, M., & Valeri, L. (2019). Med4way: A stata command to investigate mediating and interactive mechanisms using the four-way effect decomposition. International Journal of Epidemiology, 48(1), 1520. https://doi.org/10.1093/ije/dyy236CrossRefGoogle Scholar
Drabwell, L., Eng, J., Stevenson, F., King, M., Osborn, D., & Pitman, A. (2020). Perceptions of the use of alcohol and drugs after sudden bereavement by unnatural causes: Analysis of online qualitative data. International Journal of Environmental Research and Public Health, 17(3), 677. https://doi.org/10.3390/ijerph17030677CrossRefGoogle ScholarPubMed
Eng, J., Drabwell, L., Stevenson, F., King, M., Osborn, D., & Pitman, A. (2019). Use of alcohol and unprescribed drugs after suicide bereavement: Qualitative study. International Journal of Environmental Research and Public Health, 16(21), 4093. https://doi.org/10.3390/ijerph16214093CrossRefGoogle ScholarPubMed
Entilli, L., Ross, V., de Leo, D., Cipolletta, S., & Kõlves, K. (2021). Experiences of parental suicide-bereavement: A longitudinal qualitative analysis over two years. International Journal of Environmental Research and Public Health, 18(2), 564. https://doi.org/10.3390/ijerph18020564CrossRefGoogle ScholarPubMed
Erlangsen, A., & Fedyszyn, I. (2015). Danish nationwide registers for public health and health-related research. Scandinavian Journal of Public Health, 43(4), 333339. https://doi.org/10.1177/1403494815575193CrossRefGoogle ScholarPubMed
Erlangsen, A., Runeson, B., Bolton, J. M., Wilcox, H. C., Forman, J. L., Krogh, J., … Conwell, Y. (2017). Association between spousal suicide and mental, physical, and social health outcomes a longitudinal and nationwide register-based study. JAMA Psychiatry, 74(5), 456464. https://doi.org/10.1001/jamapsychiatry.2017.0226CrossRefGoogle ScholarPubMed
Fazel, S., & Runeson, B. (2020). Suicide. New England Journal of Medicine, 382, 266274.10.1056/NEJMra1902944CrossRefGoogle ScholarPubMed
Grissa, D., Rasmussen, D. N., Krag, A., Brunak, S., & Jensen, L. J. (2020). Alcoholic liver disease: A registry view on comorbidities and disease prediction. PLoS Computational Biology, 16(9), e1008244e1008244. https://doi.org/10.1371/journal.pcbi.1008244CrossRefGoogle ScholarPubMed
Hamdan, S., Melhem, N. M., Porta, G., Song, M. S., & Brent, D. A. (2013). Alcohol and substance abuse in parentally bereaved youth. Journal of Clinical Psychiatry, 74(8), 828833. https://doi.org/10.4088/JCP.13m08391CrossRefGoogle ScholarPubMed
Hanschmidt, F., Lehnig, F., Riedel-Heller, S. G., & Kersting, A. (2016). The stigma of suicide survivorship and related consequences - A systematic review. PLoS ONE, 11(9), E0162688. https://doi.org/10.1371/journal.pone.0162688CrossRefGoogle ScholarPubMed
Helweg-Larsen, K. (2011). The Danish register of causes of death. Scandinavian Journal of Public Health, 39(7), 2629. https://doi.org/10.1177/1403494811399958CrossRefGoogle ScholarPubMed
Holmberg, M. J., & Andersen, L. W. (2022). Collider Bias. JAMA – Journal of the American Medical Association, 327(13), 12821283. https://doi.org/10.1001/jama.2022.1820CrossRefGoogle ScholarPubMed
Hunt, Q. A., Young, T. A., & Hertlein, K. M. (2019). The process of long-term suicide bereavement: Responsibility, familial support, and meaning making. Contemporary Family Therapy, 41(4), 335346. https://doi.org/10.1007/s10591-019-09499-5CrossRefGoogle Scholar
Kuramoto, S. J., Stuart, E. A., Runeson, B., Lichtenstein, P., Långström, N., & Wilcox, H. C. (2010). Maternal or paternal suicide and offspring's psychiatric and suicide-attempt hospitalization risk. Pediatrics, 126(5), e10261032. https://doi.org/10.1542/peds.2010-0974CrossRefGoogle ScholarPubMed
Linde, K., Treml, J., Steinig, J., Nagl, M., & Kersting, A. (2017). Grief interventions for people bereaved by suicide: A systematic review. PLoS ONE, 12(6), e0179496. https://doi.org/10.1371/journal.pone.0179496CrossRefGoogle ScholarPubMed
MacKinnon, D. P., Fairchild, A. J., & Fritz, M. S. (2007). Mediation analysis. Annual Review of Psychology, 58(1), 593614. https://doi.org/10.1146/annurev.psych.58.110405.085542CrossRefGoogle ScholarPubMed
McDaid, C., Trowman, R., Golder, S., Hawton, K., & Sowden, A. (2008). Interventions for people bereaved through suicide: Systematic review. British Journal of Psychiatry, 193(6), 438443. https://doi.org/10.1192/bjp.bp.107.040824CrossRefGoogle ScholarPubMed
Mors, O., Perto, G. P., & Mortensen, P. B. (2011). The Danish psychiatric central research register. Scandinavian Journal of Public Health, 39(7), 5457. https://doi.org/10.1177/1403494810395825CrossRefGoogle ScholarPubMed
Musliner, K. L., Liu, X., Gasse, C., Christensen, K. S., Wimberley, T., & Munk-Olsen, T. (2019). Incidence of medically treated depression in Denmark among individuals 15–44 years old: A comprehensive overview based on population registers. Acta Psychiatrica Scandinavica, 139(6), 548557. https://doi.org/10.1111/acps.13028CrossRefGoogle Scholar
O'Connor, R. C., & Nock, M. K. (2014). The psychology of suicidal behaviour. The Lancet Psychiatry, 1(1), 7385. https://doi.org/10.1016/S2215-0366(14)70222-6CrossRefGoogle ScholarPubMed
Pedersen, C. B. (2011). The Danish civil registration system. Scandinavian Journal of Public Health, 39(7), 2225. https://doi.org/10.1177/1403494810387965CrossRefGoogle ScholarPubMed
Pitman, A., Khrisna Putri, A., de Souza, T., Stevenson, F., King, M., Osborn, D., & Morant, N. (2018). The impact of suicide bereavement on educational and occupational functioning: A qualitative study of 460 bereaved adults. International Journal of Environmental Research and Public Health, 15(4), 643. https://doi.org/10.3390/ijerph15040643CrossRefGoogle ScholarPubMed
Pitman, A., McDonald, K., Logeswaran, Y., Lewis, G., Cerel, J., & Erlangsen, A. (2022). Proportion of suicides in Denmark attributable to bereavement by the suicide of a first-degree relative or partner: Nested case–control study. Acta Psychiatrica Scandinavica, 146(6), 529539. https://doi.org/10.1111/acps.13493CrossRefGoogle ScholarPubMed
Pitman, A., Nesse, H., Morant, N., Azorina, V., Stevenson, F., King, M., & Osborn, D. (2017). Attitudes to suicide following the suicide of a friend or relative: A qualitative study of the views of 429 young bereaved adults in the UK. BMC Psychiatry, 17(1), 400. https://doi.org/10.1186/s12888-017-1560-3CrossRefGoogle ScholarPubMed
Pitman, A., Osborn, D., King, M., & Erlangsen, A. (2014). Effects of suicide bereavement on mental health and suicide risk. The Lancet Psychiatry, 1(1), 8694. https://doi.org/10.1016/S2215-0366(14)70224-XCrossRefGoogle ScholarPubMed
Pitman, A., Stevenson, F., King, M., & Osborn, D. (2020). Self-reported patterns of use of alcohol and drugs after suicide bereavement and other sudden losses: A mixed methods study of 1854 young bereaved adults in the UK. Frontiers in Psychology, 11, 1024. https://doi.org/10.3389/fpsyg.2020.01024CrossRefGoogle Scholar
Richiardi, L., Bellocco, R., & Zugna, D. (2013). Mediation analysis in epidemiology: Methods, interpretation and bias. International Journal of Epidemiology, 42(5), 15111519. https://doi.org/10.1093/ije/dyt127CrossRefGoogle ScholarPubMed
Rijnhart, J. J. M., Valente, M. J., MacKinnon, D. P., Twisk, J. W. R., & Heymans, M. W. (2021). The use of traditional and causal estimators for mediation models with a binary outcome and exposure-mediator interaction. Structural Equation Modeling, 28(3), 345355. https://doi.org/10.1080/10705511.2020.1811709CrossRefGoogle ScholarPubMed
SAS system for SunOP (9.4). (2003). SAS Institute Inc.Google Scholar
Schlichthorst, M., Reifels, L., Spittal, M., Clapperton, A., Scurrah, K., Kolves, K., … Krysinska, K. (2022). Evaluating the effectiveness of components of national suicide prevention strategies. Crisis, 44(4), 318328. https://doi.org/10.1027/0227-5910/a000887CrossRefGoogle ScholarPubMed
Schmidt, M., Schmidt, S. A. J., Sandegaard, J. L., Ehrenstein, V., Pedersen, L., & Sørensen, H. T. (2015). The Danish national patient registry: A review of content, data quality, and research potential. Clinical Epidemiology, 7, 449490. https://doi.org/10.2147/CLEP.S91125CrossRefGoogle ScholarPubMed
Spiwak, R., Elias, B., Sareen, J., Chartier, M., Katz, L. Y., & Bolton, J. M. (2020). Spouses bereaved by suicide: A population-based longitudinal cohort comparison of physician-diagnosed mental disorders and hospitalized suicide attempts. Journal of Psychiatric Research, 130, 347354. https://doi.org/10.1016/j.jpsychires.2020.07.028CrossRefGoogle ScholarPubMed
StataCorp. (2019). Stata statistical software: Release 16. College Station, Texas, USA: StataCorp LLC.Google Scholar
StataCorp. (2021). Stata statistical software: Release 17. College Station, Texas, USA: StataCorp LLC.Google Scholar
Szumilas, M., & Kutcher, S. (2011). Post-suicide intervention programs: A systematic review. Canadian Journal of Public Health, 102(1), 1829. https://doi.org/10.1007/bf03404872CrossRefGoogle ScholarPubMed
Tøllefsen, I. M., Helweg-Larsen, K., Thiblin, I., Hem, E., Kastrup, M. C., Nyberg, U., … Ekeberg, Ø (2015). Are suicide deaths under-reported? Nationwide re-evaluations of 1800 deaths in Scandinavia. BMJ Open, 5(11), e009120. https://doi.org/10.1136/bmjopen-2015-009120CrossRefGoogle ScholarPubMed
Vanderweele, T. J. (2014). A unification of mediation and interaction: A 4-way decomposition. Epidemiology (Cambridge, Mass.), 25(5), 749761. https://doi.org/10.1097/EDE.0000000000000121CrossRefGoogle ScholarPubMed
VanderWeele, T. J. (2016a). Mediation analysis: A practitioner's guide. Annual Review of Public Health, 37, 1732. https://doi.org/10.1146/annurev-publhealth-032315-021402CrossRefGoogle ScholarPubMed
VanderWeele, T. J. (2016b). Commentary: On causes, causal inference, and potential outcomes. International Journal of Epidemiology, 45(6), 18091816. https://doi.org/10.1093/ije/dyw230Google ScholarPubMed
Wilcox, H. C., Kuramoto, S. J., Lichtenstein, P., Långström, N., Brent, D. A., & Runeson, B. (2010). Psychiatric morbidity, violent crime, and suicide among children and adolescents exposed to parental death. Journal of the American Academy of Child and Adolescent Psychiatry, 49(5), 514523. https://doi.org/10.1097/00004583-201005000-00012Google ScholarPubMed
World Health Organization. (2018). National suicide prevention strategies: Progress, examples and indicators. Geneva, Switzerland: World Health Organization.Google Scholar
Yang, L., Wong, L. Y., Grivel, M. M., & Hasin, D. S. (2017). Stigma and substance use disorders: An international phenomenon. Current Opinion in Psychiatry, 30(5), 378388. https://doi.org/10.1097/YCO.0000000000000351CrossRefGoogle ScholarPubMed
Figure 0

Figure 1. Hypothesized model for the mediating effect of (A) depression and (B) substance use on the association between suicide bereavement and suicide.Legend: Panel A displays the hypothetical model for mediation by depression. Panel B displays the hypothetical model for mediation by substance use. Solid lines represent potential causal pathways. Dashed lines reflect potential confounding pathways. Path A represents the association between suicide bereavement and post-bereavement depression/substance use (modeled using logistic regression). Path B represents the association between post-bereavement depression/substance use (in all those bereaved) and suicide modeled using logistic regression. Path C represents the association between suicide bereavement and suicide (modeled using Cox proportional hazards regression). The same set of confounders was used for all pathways: sex, age, bereavement year, legal marital status, household income level, pre-bereavement history of self-harm, any psychiatric disorders, and any physical disorders.

Figure 1

Table 1. Demographic and clinical features of the studied cohort (n = 936 070) according to bereavement status

Figure 2

Table 2. Hazard ratios and 95% confidence intervals for the association between suicide bereavement and suicide compared with other bereavement (Fig. 1, Path C)

Figure 3

Table 3. Odds ratios and 95% confidence intervals for each pathway (modeled separately) in the model for mediation of the association between suicide bereavement and suicide by depression (Fig. 1, Panel A) and substance use (Fig. 1, Panel B)

Figure 4

Table 4. Mediation analyses describing the role of depression in the association between suicide bereavement and suicide

Figure 5

Table 5. Mediation analyses describing the role of substance use in the association between suicide bereavement and suicide

Supplementary material: File

Pitman et al. supplementary material

Pitman et al. supplementary material
Download Pitman et al. supplementary material(File)
File 337.3 KB