Male-specific association of the 2p25 region with suicide attempt in bipolar disorder

We previously conducted a genome-wide association study (GWAS) of attempted suicide within bipolar disorder, which implicated common variation in the 2p25 region primarily in males. The top association signal from our GWAS occurred in an intergenic region of 2p25 (p = 5.07 × 10⁻⁸) and was supported by two independent studies. In the current study, to better characterize the association of the 2p25 region with attempted suicide, we sequenced the entire 350kb 2p25 region in 476 bipolar suicide attempters and 473 bipolar non-attempters using targeted next-generation sequencing. This fine-mapping project identified 4,681 variants in the 2p25 region. We performed both gene-level and individual-variant tests on our sequencing results and identified 375 variants which were nominally significant (p < 0.05) and three common variants that were significantly associated with attempted suicide in males (corrected p = 0.035, odds ratio (OR) = 2.13). These three variants are in strong linkage disequilibrium with the top variant from our GWAS. Our top five variants are also predicted expression quantitative trait loci (eQTL) for three genes in the 2p25 region based on publicly available brain expression databases. Our sequencing and eQTL data implicate these three genes – SH3YL1, ACP1, and FAM150B – and three additional pathways – androgen receptor, Wnt signaling, and glutamatergic/GABAergic signaling – in the association of the 2p25 region with suicide. The current study provides additional support for an association of the 2p25 region with attempted suicide in males and identifies several candidate genes and pathways that warrant further investigation to understand their role in suicidal behavior.