Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

A genome-wide association study of suicide attempts in the million veterans program identifies evidence of pan-ancestry and ancestry-specific risk loci

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

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Manhattan plots summarizing the GWAS results for the pan-ancestry meta-analysis and ancestry-specific GWAS.

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

  1. World Health Organization. Suicide in the world: Global estimates. 2019. https://www.who.int/publications/i/item/suicide-in-the-world.

  2. 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.

  3. 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.

  4. Voracek M, Loibl LM. Genetics of suicide: A systematic review of twin studies. Wien Klin Wochenschr. 2007;119:463–75.

    CAS  PubMed  Google Scholar 

  5. 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.

    PubMed  PubMed Central  Google Scholar 

  6. 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.

    CAS  Google Scholar 

  7. 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.

    PubMed  PubMed Central  Google Scholar 

  8. 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.

    PubMed  Google Scholar 

  9. 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.

    PubMed  Google Scholar 

  10. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. 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.

    PubMed  PubMed Central  Google Scholar 

  12. Center for Statistical Genetics. Minimac4. https://genome.sph.umich.edu/wiki/Minimac4.

  13. 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.

    PubMed  PubMed Central  Google Scholar 

  14. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Viechtbauer W. Conducting meta-analyses in R with the metafor package. Journal of Statistical Software. 2010;36:1–48.

    Google Scholar 

  16. Higgins JP, Thompson SG. Quantifying heterogeneity in a meta-analysis. Stat Med. 2002;21:1539–58.

    PubMed  Google Scholar 

  17. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. 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.

    CAS  PubMed  Google Scholar 

  20. Euesden J, Lewis CM, O’Reilly PF. PRSice: Polygenic Risk Score Software. Bioinformatics. 2015;31:1466–8.

    CAS  PubMed  Google Scholar 

  21. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Schizophrenia Psychiatric Genome-Wide Association Study Consortium. Genome-wide association study identifies five new schizophrenia loci. Nat Genet. 2011;43:969–76.

    Google Scholar 

  24. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. 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.

    PubMed  PubMed Central  Google Scholar 

  26. UK Biobank. http://www.nealelab.is/uk-biobank.

  27. 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.

    CAS  PubMed  Google Scholar 

  28. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  29. 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.

  30. Suicide Prevention Resource Center. 2019. https://www.sprc.org/scope/racial-ethnic-disparities.

  31. 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.

  32. 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.

    CAS  PubMed  Google Scholar 

  33. Comings DE, Blum K. Reward deficiency syndrome: genetic aspects of behavioral disorders. Prog Brain Res. 2000;126:325–41.

    CAS  PubMed  Google Scholar 

  34. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  36. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  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.

    CAS  Google Scholar 

  39. 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.

    CAS  PubMed  Google Scholar 

  40. Human Protein Atlas. http://www.proteinatlas.org/.

  41. 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.

  42. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. 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.

    CAS  Google Scholar 

  44. 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.

    Google Scholar 

  45. 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.

    CAS  PubMed  Google Scholar 

  46. 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.

    CAS  PubMed  Google Scholar 

  47. Fazel S, Runeson B. Suicide. N Engl J Med. 2020;382:266–74.

    PubMed  PubMed Central  Google Scholar 

  48. 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.

    PubMed  Google Scholar 

  49. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. 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.

    Google Scholar 

  51. 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.

    CAS  PubMed  Google Scholar 

  52. 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.

    CAS  PubMed  Google Scholar 

  53. 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.

    CAS  Google Scholar 

  54. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  55. 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.

    Google Scholar 

  56. Liao X, Li Y. Genetic associations between voltage-gated calcium channels and autism spectrum disorder: A systematic review. Molecular Brain. 2020;13:96.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. 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.

    CAS  Google Scholar 

  58. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. 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.

    CAS  PubMed  Google Scholar 

  60. 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.

    PubMed  Google Scholar 

  61. 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.

    PubMed  Google Scholar 

  62. 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.

    CAS  PubMed  Google Scholar 

  63. 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.

    Google Scholar 

  64. 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.

    CAS  PubMed  Google Scholar 

  65. 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.

    CAS  PubMed  Google Scholar 

  66. 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.

    Google Scholar 

  67. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  68. Walker WH 2nd, Walton JC, DeVries AC, Nelson RJ. Circadian rhythm disruption and mental health. Trans Psychiatry. 2020;10:28.

    Google Scholar 

  69. 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.

    CAS  PubMed  Google Scholar 

  70. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

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

Authors

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

Correspondence to Nathan A. Kimbrel.

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41380-022-01472-3

This article is cited by

Search

Quick links