Abstract
To identify pan-ancestry and ancestry-specific loci associated with attempting suicide among veterans, we conducted a genome-wide association study (GWAS) of suicide attempts within a large, multi-ancestry cohort of U.S. veterans enrolled in the Million Veterans Program (MVP). Cases were defined as veterans with a documented history of suicide attempts in the electronic health record (EHR; N = 14,089) and controls were defined as veterans with no documented history of suicidal thoughts or behaviors in the EHR (N = 395,064). GWAS was performed separately in each ancestry group, controlling for sex, age and genetic substructure. Pan-ancestry risk loci were identified through meta-analysis and included two genome-wide significant loci on chromosomes 20 (p = 3.64 × 10−9) and 1 (p = 3.69 × 10−8). A strong pan-ancestry signal at the Dopamine Receptor D2 locus (p = 1.77 × 10−7) was also identified and subsequently replicated in a large, independent international civilian cohort (p = 7.97 × 10−4). Additionally, ancestry-specific genome-wide significant loci were also detected in African-Americans, European-Americans, Asian-Americans, and Hispanic-Americans. Pathway analyses suggested over-representation of many biological pathways with high clinical significance, including oxytocin signaling, glutamatergic synapse, cortisol synthesis and secretion, dopaminergic synapse, and circadian rhythm. These findings confirm that the genetic architecture underlying suicide attempt risk is complex and includes both pan-ancestry and ancestry-specific risk loci. Moreover, pathway analyses suggested many commonly impacted biological pathways that could inform development of improved therapeutics for suicide prevention.
This is a preview of subscription content, access via your institution
Access options
Subscribe to this journal
Receive 12 print issues and online access
$259.00 per year
only $21.58 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Code availability
Please note that the code used to phenotype suicide attempts and suicidal ideation from VA EHR data in the present study are available through the VA’s Centralized Interactive Phenomics Resource (CIPHER) https://www.research.va.gov/programs/cipher.cfm (VA network access only).
References
World Health Organization. Suicide in the world: Global estimates. 2019. https://www.who.int/publications/i/item/suicide-in-the-world.
Hedegaard H, Curtin S, Warner M NCHS Data Brief No. 398. Feb 2021. In Data Briefs. https://www.cdc.gov/nchs/products/databriefs/db398.htm#:~:text=This%20report%20highlights%20trends%20in,than%20the%20rate%20in%202018.
Department of Veterans Affairs Office of Mental Health and Suicide Prevention. 2020 National Strategy for Preventing Veteran Suicide. https://www.mentalhealth.va.gov/suicide_prevention/docs/Office-of-Mental-Health-and-Suicide-Prevention-National-Strategy-for-Preventing-Veterans-Suicide.pdf.
Voracek M, Loibl LM. Genetics of suicide: A systematic review of twin studies. Wien Klin Wochenschr. 2007;119:463–75.
Mirkovic B, Laurent C, Podlipski MA, Frebourg T, Cohen D, Gerardin P. Genetic association studies of suicidal behavior: a review of the past 10 years, progress, limitations, and future directions. Front Psychiatry. 2016;7:158.
Stein MB, Ware EB, Mitchell C, Chen CY, Borja S, Cai T, et al. Genome-wide association studies of suicide attempts in US soldiers. J Med Genet Part B. 2017;174:786–97.
Kimbrel NA, Garrett ME, Dennis MF, VA Mid-Atlantic Mental Illness Research, Education, and Clinical Center Workgroup, Hauser MA, Ashley-Koch AE, et al. A genome-wide association study of suicide attempts and suicidal ideation in U.S. military veterans. Psychiatry Res. 2018;269:64–9.
Mullins N, Kang J, Campos AI, Coleman JRI, Edwards AC, Galfalvy H, et al. Dissecting the shared genetic architecture of suicide attempt, psychiatric disorders and known risk factors. Biol Psychiatry. 2022;91:313–27.
Gaziano JM, Concato J, Brophy M, Fiore L, Pyarajan S, Breeling J, et al. Million Veteran Program: A mega-biobank to study genetic influences on health and disease. J Clin Epidemiol. 2016;70:214–23.
Hunter-Zinck H, Shi Y, Li M, Gorman BR, Ji SG, Sun N, et al. Genotyping array design and data quality control in the Million Veteran Program. Am J Hum Genet. 2020;106:535–48.
Hoffmire C, Stephens B, Morley S, Thompson C, Kemp J, Bossarte RM. VA Suicide Prevention Applications Network: A National Health Care System-Based Suicide Event Tracking System. Public Health Rep. 2016;131:816–21.
Center for Statistical Genetics. Minimac4. https://genome.sph.umich.edu/wiki/Minimac4.
Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7.
Fang H, Hui Q, Lynch J, Honerlaw J, Assimes TL, Huang J, et al. Harmonizing Genetic Ancestry and Self-identified Race/Ethnicity in Genome-wide Association Studies. Am J Hum Genet. 2019;105:763–72.
Viechtbauer W. Conducting meta-analyses in R with the metafor package. Journal of Statistical Software. 2010;36:1–48.
Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med. 2002;21:1539–58.
Pruim RJ, Welch RP, Sanna S, Teslovich TM, Chines PS, Gliedt TP, et al. LocusZoom: Regional visualization of genome-wide association scan results. Bioinformatics. 2010;26:2336–7.
Machiela MJ, Chanock SJ. LDlink: A web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics. 2015;31:3555–7.
Machiela MJ, Chanock SJ. LDassoc: An online tool for interactively exploring genome-wide association study results and prioritizing variants for functional investigation. Bioinformatics. 2018;34:887–9.
Euesden J, Lewis CM, O’Reilly PF. PRSice: Polygenic Risk Score Software. Bioinformatics. 2015;31:1466–8.
Bulik-Sullivan BK, Loh PR, Finucane HK, Ripke S, Yang J, Schizophrenia Working Group of the Psychiatric Genomics Consortium, et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet. 2015;47:291–5.
Stahl EA, Breen G, Forstner AJ, McQuillin A, Ripke S, Trubetskoy V, et al. Genome-wide association study identifies 30 loci associated with bipolar disorder. Nat Genet. 2019;51:793–803.
Schizophrenia Psychiatric Genome-Wide Association Study Consortium. Genome-wide association study identifies five new schizophrenia loci. Nat Genet. 2011;43:969–76.
Wray NR, Ripke S, Mattheisen M, Trzaskowski M, Byrne EM, Abdellaoui A, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genet. 2018;50:668–81.
Nievergelt CM, Maihofer AX, Klengel T, Atkinson EG, Chen CY, Choi KW, et al. International meta-analysis of PTSD genome-wide association studies identifies sex- and ancestry-specific genetic risk loci. Nat Commun. 2019;10:4558.
UK Biobank. http://www.nealelab.is/uk-biobank.
Zheng J, Erzurumluoglu AM, Elsworth BL, Kemp JP, Howe L, Haycock PC, et al. LD Hub: A centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis. Bioinformatics. 2017;33:272–9.
Wang J, Vasaikar S, Shi Z, Greer M, Zhang B. WebGestalt 2017: A more comprehensive, powerful, flexible and interactive gene set enrichment analysis toolkit. Nucleic Acids Research. 2017;45:W130–W137.
Center for Behavioral Health Statistics and Quality. 2019 National Survey on Drug Use and Health: Detailed Tables. Rockville, MD: Substance Abuse and Mental Health Services Administration; 2020. https://www.samhsa.gov/data/report/2019-nsduh-detailed-tables.
Suicide Prevention Resource Center. 2019. https://www.sprc.org/scope/racial-ethnic-disparities.
Centers for Disease Control and Prevention National Center for Health Statistics. 1999-2019 Wide Ranging Online Data for Epidemiological Research (WONDER), Multiple Cause of Death files [Data file]. 2021. http://wonder.cdc.gov/ucd-icd10.html.
Jenkins D, Caubit X, Dimovski A, Matevska N, Lye CM, Cabuk F, et al. Analysis of TSHZ2 and TSHZ3 genes in congenital pelvi-ureteric junction obstruction. Nephrol Dial Transplant. 2010;25:54–60.
Comings DE, Blum K. Reward deficiency syndrome: genetic aspects of behavioral disorders. Prog Brain Res. 2000;126:325–41.
Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. 2017;8:1826 https://doi.org/10.1038/s41467-017-01261-5.
Ogi T, Walker S, Stiff T, Hobson E, Limsirichaikul S, Carpenter G, et al. Identification of the first ATRIP-deficient patient and novel mutations in ATR define a clinical spectrum for ATR-ATRIP Seckel Syndrome. PLoS Genetics. 2012;8:e1002945.
Hoppa MB, Lana B, Margas W, Dolphin AC, Ryan TA. alpha2delta expression sets presynaptic calcium channel abundance and release probability. Nature. 2012;486:122–5.
Neely GG, Hess A, Costigan M, Keene AC, Goulas S, Langeslag M, et al. A genome-wide Drosophila screen for heat nociception identifies alpha2delta3 as an evolutionarily conserved pain gene. Cell. 2010;143:628–38.
Chung HW, Petersen EN, Cabanos C, Murphy KR, Pavel MA, Hansen AS, et al. A molecular target for an alcohol chain-length cutoff. J Mol Biology. 2019;431:196–209.
Cruchaga C, Karch CM, Jin SC, Benitez BA, Cai Y, Guerreiro R, et al. Rare coding variants in the phospholipase D3 gene confer risk for Alzheimer’s disease. Nature. 2014;505:550–4.
Human Protein Atlas. http://www.proteinatlas.org/.
Gottlieb DJ, O’Connor GT, Wilk JB. Genome-wide association of sleep and circadian phenotypes. BMC Med Genet. 2007;19:S9. https://doi.org/10.1186/1471-2350-8-S1-S9.
Edwards AC, Aliev F, Bierut LJ, Bucholz KK, Edenberg H, Hesselbrock V, et al. Genome-wide association study of comorbid depressive syndrome and alcohol dependence. Psychiatr Genet. 2012;22:31–41.
Lopez-Leon S, Gonzalez-Giraldo Y, Wegman-Ostrosky T, Forero DA. Molecular genetics of substance use disorders: An umbrella review. Neuro Biobehav Rev. 2021;124:358–69.
Gonzalez-Castro TB, Hernandez-Diaz Y, Juarez-Rojop IE, López-Narváez ML, Tovilla-Zárate CA, Genis-Mendoza A, et al. The role of C957T, TaqI and Ser311Cys polymorphisms of the DRD2 gene in schizophrenia: systematic review and meta-analysis. Behav Brain Functions. 2016;12:29.
Pan YQ, Qiao L, Xue XD, Fu JH. Association between ANKK1 (rs1800497) polymorphism of DRD2 gene and attention deficit hyperactivity disorder: A meta-analysis. Neurosci Lett. 2015;590:101–5.
Cade BE, Gottlieb DJ, Lauderdale DS, Bennett DA, Buchman AS, Buxbaum SG, et al. Common variants in DRD2 are associated with sleep duration: The CARe Consortium. Hum Mol Genet. 2016;25:167–79.
Fazel S, Runeson B. Suicide. N Engl J Med. 2020;382:266–74.
Septier M, Stordeur C, Zhang J, Delorme R, Cortese S. Association between suicidal spectrum behaviors and Attention-Deficit/Hyperactivity Disorder: A systematic review and meta-analysis. Neurosci Biobehav Rev. 2019;103:109–18.
Harris LM, Huang X, Linthicum KP, Bryen CP, Ribeiro JD. Sleep disturbances as risk factors for suicidal thoughts and behaviours: A meta-analysis of longitudinal studies. Sci Rep. 2020;10:13888.
Johann M, Putzhammer A, Eichhammer P, Wodarz N. Association of the −141C Del variant of the dopamine D2 receptor (DRD2) with positive family history and suicidality in German alcoholics. Am J Med Genetics, Part B Neuropsychiatric Gen. 2005;132B:46–49.
Suda A, Kawanishi C, Kishida I, Sato R, Yamada T, Nakagawa M, et al. Dopamine D2 receptor gene polymorphisms are associated with suicide attempt in the Japanese population. Neuropsychobiology. 2009;59:130–4.
Genis-Mendoza AD, Lopez-Narvaez ML, Tovilla-Zarate CA, Sarmiento E, Chavez A, Martinez-Magaña JJ, et al. Association between polymorphisms of the DRD2 and ANKK1 genes and suicide attempt: A preliminary case-control study in a mexican population. Neuropsychobiology. 2017;76:193–8.
Smedler E, Bergen SE, Song J, Landen M. Genes, biomarkers, and clinical features associated with the course of bipolar disorder. Eur Neuropsychopharm. 2019;29:1152–60.
Dunn EC, Wiste A, Radmanesh F, Almli LM, Gogarten SM, Sofer T, et al. Genome-wide association study (GWAS) and Genome-wide by environment interaction study (GWEIS) of depressive symptoms in African American and Hispanic/Latina women. Depress Anxiety. 2016;33:265–80.
Zhang T, Zhu L, Ni T, Liu D, Chen G, Yan Z, et al. Voltage-gated calcium channel activity and complex related genes and schizophrenia: A systematic investigation based on Han Chinese population. J Psych Res. 2018;106:99–105.
Liao X, Li Y. Genetic associations between voltage-gated calcium channels and autism spectrum disorder: A systematic review. Molecular Brain. 2020;13:96.
Yin X, Bizon C, Tilson J, Lin Y, Gizer IR, Ehlers CL, et al. Genome-wide meta-analysis identifies a novel susceptibility signal at CACNA2D3 for nicotine dependence. Am J Med Gen Part B Neuropsychiatric Genet. 2017;174:557–67.
Schork AJ, Won H, Appadurai V, Nudel R, Gandal M, Delaneau O, et al. A genome-wide association study of shared risk across psychiatric disorders implicates gene regulation during fetal neurodevelopment. Nat Neurosci. 2019;22:353–61.
Collod-Beroud G, Le Bourdelles S, Ades L, Ala-Kokko L, Booms P, Boxer M, et al. Update of the UMD-FBN1 mutation database and creation of an FBN1 polymorphism database. Human Mutation. 2003;22:199–208.
Djurovic S, Gustafsson O, Mattingsdal M, Athanasiu L, Bjella T, Tesli M, et al. A genome-wide association study of bipolar disorder in Norwegian individuals, followed by replication in Icelandic sample. J Affect Disord. 2010;126:312–6.
Jahangard L, Shayganfard M, Ghiasi F, Salehi I, Haghighi M, Ahmadpanah M, et al. Serum oxytocin concentrations in current and recent suicide survivors are lower than in healthy controls. J Psychiatr Res. 2020;128:75–82.
Jokinen J, Chatzittofis A, Hellstrom C, Nordstrom P, Uvnas-Moberg K, Asberg M. Low CSF oxytocin reflects high intent in suicide attempters. Psychoneuroendocrinology. 2012;37:482–90.
Li CT, Yang KC, Lin WC. Glutamatergic dysfunction and glutamatergic compounds for major psychiatric disorders: Evidence from clinical neuroimaging studies. Front Psych. 2018;9:767.
Gray AL, Hyde TM, Deep-Soboslay A, Kleinman JE, Sodhi MS. Sex differences in glutamate receptor gene expression in major depression and suicide. Mol Psychiatry. 2015;20:1057–68.
Zhao J, Verwer RWH, Gao SF, Qi XR, Lucassen PJ, Kessels HW, et al. Prefrontal alterations in GABAergic and glutamatergic gene expression in relation to depression and suicide. J Psychiatr Res. 2018;102:261–74.
Witt K, Potts J, Hubers A, Grunebaum MF, Murrough JW, Loo C, et al. Ketamine for suicidal ideation in adults with psychiatric disorders: A systematic review and meta-analysis of treatment trials. Aus and New Zealand J Psychiatry. 2020;54:29–45.
McClung CA, Sidiropoulou K, Vitaterna M, Takahashi JS, White FJ, Cooper DC, et al. Regulation of dopaminergic transmission and cocaine reward by the Clock gene. Proc Natl Acad Sci USA. 2005;102:9377–81.
Walker WH 2nd, Walton JC, DeVries AC, Nelson RJ. Circadian rhythm disruption and mental health. Trans Psychiatry. 2020;10:28.
DiBlasi E, Shabalin AA, Monson ET, Keeshin BR, Bakian AV, Kirby AV, et al. Rare protein-coding variants implicate genes involved in risk of suicide death. Am J Med Genet B Neuropsychiatr Genet. 2021;186:508–20.
Le-Niculescu H, Roseberry K, Gill SS, Levey DF, Phalen PL, Mullen J, et al. Precision medicine for mood disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs. Mol Psychiatry. 2021;26:2776–804.
Acknowledgements
This research is based on data from the Million Veteran Program, Office of Research and Development (ORD), Veterans Health Administration (VHA), and was supported by award #I01CX001729 from the Clinical Science Research and Development (CSR&D) Service of VHA ORD. J.C. Beckham was also supported by a Senior Research Career Scientist Award (#lK6BX003777) from CSR&D. Niamh Mullins was supported by a NARSAD Young Investigator Award from the Brain & Behavior Research Foundation. This publication does not represent the views of the Department of Veteran Affairs or the United States Government. We also thank and acknowledge MVP, the MVP Suicide Exemplar Workgroup, and the ISGC for their contributions to this manuscript. A complete listing of contributors from the MVP, MVP Suicide Exemplar Workgroup, and ISGC is provided in the Supplementary Materials. This work was also supported in part by the joint U.S. Department of Veterans Affairs and US Department of Energy MVP CHAMPION program. This manuscript has been co-authored by UT-Battelle, LLC under contract no. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a nonexclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan, last accessed September 16, 2020). The GWAS summary statistics generated from this study will be available via dbGaP. The dbGaP accession assigned to the Million Veteran Program is phs001672.v1.p.
Author information
Authors and Affiliations
Consortia
Contributions
This study was conceived and designed by NAK, AAK, ERH, MAH, and JCB. Statistical analyses were conducted by XJQ, JHL, and MEG, and were supervised by AAK and ERH. NAK, AAK, ERH, MAH, JCB, RKM, DAJ, JAT, HC, ARD, NM, and DMM reviewed and interpreted statistical findings. NAK, AAK, MFD, LPH, JEH, HC, ARD, JK, NM, DMM, PDH, BHM, DWO, ERH, MAH, and JCB were involved in data acquisition and data preparation. NAK, AAK, XJQ, JHL, and MEG wrote the initial draft of the manuscript. All authors reviewed the results and approved the final version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Lists of members and their affiliations appears in the Supplementary Information.
Supplementary information
Rights and permissions
About this article
Cite this article
Kimbrel, N.A., Ashley-Koch, A.E., Qin, X.J. et al. A genome-wide association study of suicide attempts in the million veterans program identifies evidence of pan-ancestry and ancestry-specific risk loci. Mol Psychiatry 27, 2264–2272 (2022). https://doi.org/10.1038/s41380-022-01472-3
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41380-022-01472-3
This article is cited by
-
Genetic factors associated with suicidal behaviors and alcohol use disorders in an American Indian population
Molecular Psychiatry (2024)
-
Association study of the complement component C4 gene and suicide risk in schizophrenia
Schizophrenia (2024)
-
Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder
Molecular Psychiatry (2024)
-
Functional and molecular characterization of suicidality factors using phenotypic and genome-wide data
Molecular Psychiatry (2023)
-
Suicidal behavior across a broad range of psychiatric disorders
Molecular Psychiatry (2023)