Search

Recent large-scale genome-wide association studies (GWAS) have started to identify potential genetic risk loci associated with risk of suicide; however, a large portion of suicide-associated genetic factors affecting gene expression remain elusive. Dysregulated gene expression, not assessed by GWAS, may play a significant role in increasing the risk of suicide death. We performed the first comprehensive genomic association analysis prioritizing brain expression quantitative trait loci (eQTLs) within regulatory regions in suicide deaths from the Utah Suicide Genetic Risk Study (USGRS). 440,324 brain-regulatory eQTLs were obtained by integrating brain eQTLs, histone modification ChIP-seq, ATAC-seq, DNase-seq, and Hi-C results from publicly available data. Subsequent genomic analyses were conducted in whole-genome sequencing (WGS) data from 986 suicide deaths of non-Finnish European (NFE) ancestry and 415 ancestrally matched controls. Additional independent USGRS suicide deaths with genotyping array data (n = 4657) and controls from the Genome Aggregation Database were explored for WGS result replication. One significant eQTL locus, rs926308 (p = 3.24e-06), was identified. The rs926308-T is associated with lower expression of RFPL3S, a gene important for neocortex development and implicated in arousal. Gene-based analyses performed using Sherlock Bayesian statistical integrative analysis also detected 20 genes with expression changes that may contribute to suicide risk. From analyzing publicly available transcriptomic data, ten of these genes have previous evidence of differential expression in suicide death or in psychiatric disorders that may be associated with suicide, including schizophrenia and autism (ZNF501, ZNF502, CNN3, IGF1R, KLHL36, NBL1, PDCD6IP, SNX19, BCAP29, and ARSA). Electronic health records (EHR) data was further merged to evaluate if there were clinically relevant subsets of suicide deaths associated with genetic variants. In summary, our study identified one risk locus and ten genes associated with suicide risk via gene expression, providing new insight into possible genetic and molecular mechanisms leading to suicide., (© 2023. The Author(s).)

Suicide is a significant public health concern with complex etiology. Although the genetic component of suicide is well established, the scope of gene networks and biological mechanisms underlying suicide has yet to be defined. Previously, we reported genome-wide evidence that neurexin 1 (NRXN1), a key synapse organizing molecule, is associated with familial suicide risk. Here we present new evidence for two non-synonymous variants (rs78540316; P469S and rs199784139; H885Y) associated with increased familial risk of suicide death. We tested the impact of these variants on binding interactions with known partners and assessed functionality in a hemi-synapse formation assay. Although the formation of hemi-synapses was not altered with the P469S variant relative to wild-type, both variants increased binding to the postsynaptic binding partner, leucine-rich repeat transmembrane neuronal 2 (LRRTM2) in vitro. Our findings indicate that variants in NRXN1 and related synaptic genes warrant further study as risk factors for suicide death., (© 2021. The Author(s).)

Suicide is the 10th leading cause of death in the United States. Although environment has undeniable impact, evidence suggests that genetic factors play a significant role in completed suicide. We linked a resource of ~ 4500 DNA samples from completed suicides obtained from the Utah Medical Examiner to genealogical records and medical records data available on over eight million individuals. This linking has resulted in the identification of high-risk extended families (7-9 generations) with significant familial risk of completed suicide. Familial aggregation across distant relatives minimizes effects of shared environment, provides more genetically homogeneous risk groups, and magnifies genetic risks through familial repetition. We analyzed Illumina PsychArray genotypes from suicide cases in 43 high-risk families, identifying 30 distinct shared genomic segments with genome-wide evidence (p = 2.02E-07-1.30E-18) of segregation with completed suicide. The 207 genes implicated by the shared regions provide a focused set of genes for further study; 18 have been previously associated with suicide risk. Although PsychArray variants do not represent exhaustive variation within the 207 genes, we investigated these for specific segregation within the high-risk families, and for association of variants with predicted functional impact in ~ 1300 additional Utah suicides unrelated to the discovery families. None of the limited PsychArray variants explained the high-risk family segregation; sequencing of these regions will be needed to discover segregating risk variants, which may be rarer or regulatory. However, additional association tests yielded four significant PsychArray variants (SP110, rs181058279; AGBL2, rs76215382; SUCLA2, rs121908538; APH1B, rs745918508), raising the likelihood that these genes confer risk of completed suicide.