s41386-022-01319-z
s41386-022-01319-z
s41386-022-01319-z
com/npp
ARTICLE OPEN
Insomnia and restless leg syndrome (RLS) are associated with increased risk for suicidal behavior (SB), which is often comorbid with
mood or thought disorders; however, it is unclear whether these relationships are causal. We performed a two-sample Mendelian
randomization study using summary-level genetic associations with insomnia symptoms and RLS against the outcomes of risk of
major depressive disorder (MDD), bipolar disorder (BP), schizophrenia (SCZ), and SB. The inverse-variance weighted method was
used in the main analysis. We performed replication and sensitivity analyses to examine the robustness of the results. We identified
outcome cohorts for MDD (n = 170,756 cases/329,443 controls), BP (n = 20,352/31,358), SCZ (n = 69,369/236,642), SB-Cohort-2019
1234567890();,:
(n = 6569/14,996 all with MDD, BP or SCZ; and SB within individual disease categories), and SB-Cohort-2020 (n = 29,782/519,961).
Genetically proxied liability to insomnia symptoms significantly associated with increased risk of MDD (odds ratio (OR) = 1.23, 95%
confidence interval (CI) = 1.2–1.26, P = 1.37 × 10–61), BP (OR = 1.15, 95% CI = 1.07–1.23, P = 5.11 × 10–5), SB-Cohort-2019 (OR = 1.17,
95% CI = 1.07–1.27, P = 2.30 × 10–4), SB-Cohort-2019 in depressed patients (OR = 1.34, 95% CI = 1.16–1.54, P = 5.97 × 10–5), and SB-
Cohort-2020 (OR = 1.24, 95% CI = 1.18–1.3, P = 1.47 × 10–18). Genetically proxied liability to RLS did not significantly influence the
risk of any of the outcomes (all corrected P > 0.05). Results were replicated for insomnia with MDD and SB in Mass General Brigham
Biobank and were consistent in multiple lines of sensitivity analyses. In conclusion, human genetic evidence supports for the first
time a potentially independent and causal effect of insomnia on SB and encourages further clinical investigation of treatment of
insomnia for prevention or treatment of SB.
1
Division of Sleep Medicine, Beth Israel Deaconess Medical Center, Harvard Medical School, Boston, MA, USA. 2Center for Genomic Medicine, Massachusetts General Hospital,
Harvard Medical School, Boston, MA, USA. 3Harvard Medical School, Boston, MA, USA. 4Departments of Psychiatry and Neurology, Massachusetts General Hospital, Harvard
Medical School, Boston, MA, USA. 5Center for Genomic Medicine and Department of Anesthesia, Critical Care and Pain Medicine, Massachusetts General Hospital, Harvard
Medical School, Boston, MA, USA. 6Program in Medical and Population Genetics, Broad Institute, Cambridge, MA, USA. 7Present address: Mesulam Center for Cognitive Neurology
and Alzheimer’s Disease, Northwestern University, Chicago, IL, USA. 8Present address: Department of Neurology, University of California, San Francisco, San Francisco, CA, USA.
9
These authors contributed equally: Malik Nassan, Iyas Daghlas. A full list of members and their affiliations appears in the Supplementary Information.
✉email: [email protected]
genetic variants influencing the risk of BP. Overall, cases needed to meet We then performed multiple sensitivity analyses to assess the
DSM-IV, International Classification of Diseases (ICD)-9, or ICD-10 criteria for robustness of the findings. First, we removed SNPs in LD between the
a lifetime diagnosis of BP (by either structured diagnostic instrument, insomnia symptoms and RLS genetic proxies to reduce any confounding
clinician-administered checklists, or medical record review). Most controls between these exposures (removing SNPs within 10 Mb and with r2 ≥ 0.7).
were examined for the absence of any other lifetime psychiatric disorders. Second, we removed SNPs from the exposure genetic proxies (i.e.,
A total of 30 genome-wide independent variants were associated with insomnia symptoms and RLS) that were in LD (using a more stringent
MDD in the meta-analysis. threshold r2 ≥ 0.01) with an outcome that was significant in the main MR
For SCZ, the most recent meta-GWAS included a total of 69,369 cases IVW analysis (MDD, BP, and SB-Cohort-2020). Third, we performed model-
and 236,642 controls [40]. This was the largest combined cohort yet from based sensitivity analyses that relax various MR assumptions regarding
PGC (included 90 cohorts) and identified 270 independent genome-wide pleiotropy, including: MR-egger regression, weighted median, and MR-
significant genetic loci. Of note, this study included 80% of the sample PRESSO [44, 45]. We also created leave-one-out plots to display the results
from European ancestry and 20% from East Asian ancestry. Cases included from the IVW and Egger regression analyses, to assess for outliers.
diagnoses of SCZ or schizoaffective disorder; details of each of the enrolled We sought to replicate our findings in the Mass General Brigham (MGB)
cohorts are available in the manuscript [40]. A total of 270 genome-wide Biobank (formerly Partners Biobank) for the outcomes of MDD, BP, and SB
independent variants were associated with MDD in the meta-analysis. only (MDD n = 4640/23,849, BP n = 145/28,344, SB n = 1054/27,435) [46].
Two cohorts were used for the outcome of SB (encompassing here both The MGB Biobank is a hospital-based cohort study from the MGB
fatal and non-fatal suicidal attempts). The first cohort is a PGC cohort (SB- healthcare network in Boston, MA with electronic health record (EHR)
Cohort-2019) published in 2019 (cases 6569, controls 14,996) and stratified and genetic data. Recruitment for the Biobank launched in 2010 and
by comorbid psychiatric disorder (MDD, BP, or SCZ) [41]. The subjects were remains ongoing at participating clinics and electronically. Recruitment
obtained from 16 MDD cohorts, 21 BP cohorts, and 9 SCZ cohorts from strategy has been described previously [47]. All recruited patients provided
PGC, where data on suicide attempts had been gathered. Only patients consent written informed upon enrollment. The present study protocol
affected by the three psychiatric disorders were included and all three was approved by the MGB Institutional Review Board (#2018P002276).
psychiatric disorders were defined using structured psychiatric interviews. Effect estimates for MDD, BP, and SB were generated using data for 30,683
All individuals were of European ancestry. Items from structured clinical participants with genetic data and limited to participants of European
interviews offered data on suicidal attempts. Lifetime suicidal attempt was ancestry [48]. Cases of MDD, BP, and SB were determined from EHR using a
characterized across cohorts as an intentional act of self-harm with the validated algorithm based on natural language processing of structured
intent to result in death. Individuals who only endorsed suicidal ideation and unstructured data including coded diagnoses, medications, proce-
were not included as cases. Across the cohorts, there was a sum of 6569 dures, and vital signs [46]. The remaining participants were set as control.
individuals who attempted suicide and 17,232 individuals who had not To determine SNP effects on MDD, BP, and SB, we performed genetic
attempted suicide. No genome-wide significant associations were identi- association analysis in unrelated participants of European ancestry with
fied in the meta-analysis. PLINK logistic regression and an additive genetic model adjusted for age,
The second cohort (SB-Cohort-2020) assessed for SB in 29,782 cases and sex, five principal components, and genotyping array [49].
519,961 controls and is the most recent and largest SB GWAS [pre-print Finally, we ran a reverse MR IVW analysis between MDD and BP (using
released in 2020, and publication was in 2021] [42, 43]. The cohorts included the genome-wide genetic proxies from the same GWAS) as exposures and
samples from European ancestry (the majority of the cases), admixed African insomnia symptoms as outcome (from the same insomnia symptoms
American ancestry (4%), and East Asian ancestry (6%). This GWAS included GWAS but only including the UKB cohort due to lack of public availability
21 cohorts, of which cases were individuals who died by suicide (2 cohorts) of 23andMe data).
or made a non-fatal suicide attempt (19 cohorts) which is defined as a
lifetime act of intentional self-harm with intent to cause one’s own death.
Individuals who only endorsed suicidal ideation were not included as cases. RESULTS
Information on suicidal attempts was obtained via structured clinical Main analyses
interviews for 15 cohorts, self-report questionnaires for 2 cohorts, and ICD
codes or hospital records for 2 cohorts. Cases of death by suicide (2 cohorts) MR analyses showed significant associations of genetically proxied
were obtained from the Medical Examiner’s Office of the Hyogo Prefecture insomnia symptoms with MDD, BP and SB [MDD (OR = 1.23, 95%
and the Division of Legal Medicine, at the Kobe University Graduate School CI = 1.2–1.26, P = 1.37 × 10–61), BP (OR = 1.15, 95% CI = 1.07–1.23,
of Medicine in Japan or the Utah State Office of the Medical Examiner. A P = 5.11 × 10–5), SB-Cohort-2019 (OR = 1.17, 95% CI = 1.07–1.27,
percentage of cases from the Columbia University and iPSYCH cohorts P = 2.30 × 10–4), and SB-Cohort-2020 (OR = 1.24, 95% CI =
included individuals with death by suicide that was established through the 1.18–1.3, P = 1.47 × 10–18)] [Fig. 2]. Analyses for the outcome of
Columbia Classification Algorithm for Suicide Assessment and the Cause of SB in the 2019 cohort stratified by disease status showed that the
Death Register in Denmark, respectively. Two genome-wide significant effect of insomnia symptoms on SB is most robust in the
genetic variants were identified as influencing the risk of SB. depressed population (OR = 1.34, 95% CI = 1.16–1.54, P = 5.97 ×
10–5) [Fig. 2]. The scatterplots for the significant findings are
Mendelian randomization analyses included in Fig. 3 and Supplementary Table 3S. On the other hand,
All analyses were performed in R Version 3.5.3 using the TwoSampleMR genetically proxied liability to RLS was not associated with any of
v0.4.229 package. The inverse-variance weighted (IVW) method was the the study outcomes [Fig. 4 and Supplementary Table 4S].
main MR method used to estimate the effect of genetically proxied liability
to insomnia symptoms or RLS on each of the psychiatric outcomes [44].
The corrected statistical significance threshold is P less than 3.13 × 10–3, Sensitivity analyses
accounting for 16 statistical comparisons across 2 exposures and 8 Results for the effects of genetic liability to insomnia symptoms on
outcomes. the outcomes were unchanged when we removed SNPs in LD
Fig. 2 Forest plot of inverse-variance weighted MR effects of genetically proxied liability to insomnia symptoms on psychiatric
outcomes. Boxes reflect point estimates and surrounding lines reflect 95% confidence intervals. CI confidence interval, OR odds ratio, BP
bipolar disorder, MDD major depressive disorder, RLS restless leg syndrome, SB suicidal behavior, SCZ schizophrenia.
Fig. 3 Mendelian randomization scatterplots for effects of genetic liability to insomnia on psychiatric outcomes. BP bipolar disorder, MDD
major depressive disorder, SB suicidal behavior, SCZ schizophrenia.
between insomnia symptoms and RLS genetic proxies, and when and may reflect heterogeneity in the insomnia discovery GWAS
we removed SNPs in LD between insomnia symptoms and the based on undiagnosed RLS or pleiotropy at this locus [50]. For BP
significant outcomes (MDD, BP, and SB-cohort-2020) [Supplemen- another single outlier SNP (rs9527083 intergenic in chromosome
tary Tables 5S–9S]. 13) was driving the effect in the opposite direction. For SB-Cohort-
Model-based MR sensitivity analyses were performed to assess 2019 no single outlier was driving in the opposite direction. For
the robustness of effects on potential horizontal pleiotropy comparison, we ran leave-one-out analyses for the IVW method
[Table 2]. Results from the weighted median and MR-PRESSO and found no clear outliers [Supplementary Figs. S6–S10].
were overall similar to results from the IVW analysis. However, the We sought to replicate the findings in the MGB Biobank, an
MR-Egger regression for the effect of genetically proxied liability independent clinical biobank, using available data for MDD, BP
to insomnia symptoms on SB had non-significant findings, and a and SB. These analyses showed replication of the findings for the
point estimate in the opposite direction. We performed a leave- effect of genetic liability to insomnia symptoms on MDD (OR =
one-out MR-Egger analysis to assess whether the results were 1.12, 95% CI = 1.03–1.21, P = 0.0096) and SB (OR = 1.19, 95% CI =
driven by an outlier [Supplementary Figs. S1–S5]. This showed that 1.01–1.40, P = 0.03) but not for BP (P = 0.58). The MGB cohort had
a single outlier SNP (rs113851554 in MEIS1) was driving the effect a small BP case sample size (n = 145) and consequently large
in the opposite direction for MDD, SB in MDD, and SB-Cohort- confidence intervals in the analysis, but the effect was in the same
2020. Notably, this MEIS1 SNP is also a strong risk factor for RLS direction (OR = 1.14) [Supplementary Table 10S]. Finally, the
Fig. 4 Forest plot of inverse-variance weighted MR effects of genetically proxied liability to RLS on psychiatric outcomes. Boxes reflect
point estimates and surrounding lines reflect 95% confidence intervals. CI confidence interval, OR odds ratio, BP bipolar disorder, MDD major
depressive disorder, RLS restless leg syndrome, SB suicidal behavior, SCZ schizophrenia.
reverse MR IVW analysis between MDD and BP as exposures and outcomes, indicating that our findings for insomnia were not
insomnia symptoms as outcome showed that genetic liability to driven by RLS. These results were replicated, and consistent across
MDD but not BP is a risk factor for insomnia symptoms several lines of sensitivity analyses. This is the first comprehensive
[Supplementary Table 11S]. study analyzing the causality of insomnia symptoms and RLS for
SB utilizing MR.
We found a robust association between insomnia symptoms
DISCUSSION and MDD (OR = 1.23, P = 1.7 × 10–61), which is consistent with the
In this two-sample MR study we found for the first-time evidence existing literature [34]. In addition, when we sub-classified SB by
for a potentially independent and causal effect of insomnia disorder, the most robust association of insomnia symptoms with
symptoms on SB, and further strengthened the evidence for SB was within the MDD sub-group. However, LD analyses,
insomnia being a potential causal risk factor for MDD and BP. sensitivity MR analyses, and a replication analysis in an
However, genetically proxied RLS (which can be comorbid with independent sample all demonstrated insomnia as an indepen-
insomnia) was not associated with any tested psychiatric dent risk factor for SB independently of MDD. Our results are in
ACKNOWLEDGEMENTS
We thank the GWAS consortia that have made their data publicly available. We thank
the International Suicide Genetics Consortium that shared its summary statistics
with us. This article is licensed under a Creative Commons Attribution 4.0
International License, which permits use, sharing, adaptation,
distribution and reproduction in any medium or format, as long as you give appropriate
AUTHOR CONTRIBUTIONS credit to the original author(s) and the source, provide a link to the Creative Commons
All authors contributed to drafting the work or revising it critically for important licence, and indicate if changes were made. The images or other third party material in
intellectual content, gave final approval of the version to be published; and provided this article are included in the article's Creative Commons licence, unless indicated
an agreement to be accountable for all aspects of the work. MN and ID provided otherwise in a credit line to the material. If material is not included in the article's
substantial contributions to the conception or design of the work, and the Creative Commons licence and your intended use is not permitted by statutory
acquisition, analysis, or interpretation of data for the work. JW, HSD, and RS provided regulation or exceeds the permitted use, you will need to obtain permission directly
substantial contributions to the conception or design of the work, and the from the copyright holder. To view a copy of this licence, visit http://creativecommons.
interpretation of data for the work. International Suicide Genetics Consortium org/licenses/by/4.0/.
provided summary statistics, provided substantial contributions to the interpretation
of data for the work, contributed to revising the work critically for important
intellectual content, and gave final approval of the version to be published. © The Author(s) 2022