Seasonality of suicide: a multi-country multi-community observational study

Abstract

Aims We aimed to investigate the heterogeneity of seasonal suicide patterns among multiple geographically, demographically and socioeconomically diverse populations. Methods Weekly time-series data of suicide counts for 354 communities in 12 countries during 1986–2016 were analysed. Two-stage analysis was performed. In the first stage, a generalised linear model, including cyclic splines, was used to estimate seasonal patterns of suicide for each community. In the second stage, the community-specific seasonal patterns were combined for each country using meta-regression. In addition, the community-specific seasonal patterns were regressed onto community-level socioeconomic, demographic and environmental indicators using meta-regression. Results We observed seasonal patterns in suicide, with the counts peaking in spring and declining to a trough in winter in most of the countries. However, the shape of seasonal patterns varied among countries from bimodal to unimodal seasonality. The amplitude of seasonal patterns (i.e. the peak/trough relative risk) also varied from 1.47 (95% confidence interval [CI]: 1.33–1.62) to 1.05 (95% CI: 1.01–1.1) among 12 countries. The subgroup difference in the seasonal pattern also varied over countries. In some countries, larger amplitude was shown for females and for the elderly population (≥65 years of age) than for males and for younger people, respectively. The subperiod difference also varied; some countries showed increasing seasonality while others showed a decrease or little change. Finally, the amplitude was larger for communities with colder climates, higher proportions of elderly people and lower unemployment rates (p-values < 0.05). Conclusions Despite the common features of a spring peak and a winter trough, seasonal suicide patterns were largely heterogeneous in shape, amplitude, subgroup differences and temporal changes among different populations, as influenced by climate, demographic and socioeconomic conditions. Our findings may help elucidate the underlying mechanisms of seasonal suicide patterns and aid in improving the design of population-specific suicide prevention programmes based on these patterns. ; JY, DY and YC were supported by a Senior Research grant (2019R1A2C1086194) from the National Research Foundation of Korea (NRF), funded by the Ministry of Science, ICT (Information and Communication Technologies) and Government-wide R&D Fund project for Infectious Disease (HG18C0025). YK was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI Grant Number JP19K17104 in Japan. AG and FS were supported by the Medical Research Council UK (Grant ID: MR/R013349/1) and Natural Environment Research Council (Grant ID: NE/R009384/1). NS was supported by the NIEHS-funded HERCULES Center (P30ES019776). ; Peer reviewed