Whole-genome sequencing analysis of suicide deaths integrating brain-regulatory eQTLs data to identify risk loci and genes

Abstract

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.

Effect of RNA m6A methyltransferase activation by a low molecular weight compound on anxiety- and depression-related behaviours, monoamine neurochemistry and striatal gene expression in the rat

“…It belongs to the group of cadherin-like cell adhesion molecules and knowledge on its function is limited but likely to include the establishment of specific neural connections. Nbl1, a candidate neuroblastoma suppressor gene has also been associated with suicide (Han et al, 2023).…”

Comparative whole-genome resequencing to uncover selection signatures linked to litter size in Hu Sheep and five other breeds

“…According to the reference genome annotation information, there are 39 genes located in these selected chromosome regions. We identified candidate genes associated with immune response ( BCAP29 [ 55 ], ATP8B4 [ 56 ]), reproductive traits ( CC2D1B ) [ 57 ], teat number ( SYNDIG1L ) [ 58 ], and horn phenotype ( RXFP2 ) [ 12 ]. The upregulation of critical genes during the follicular phase reflects the immune system’s role in follicular recruitment, while endocrine changes throughout the estrous cycle regulate immune gene expression, affecting follicular atresia and recruitment [ 59 ].…”

Genome-wide association study of major anxiety disorders in 122,341 European-ancestry cases identifies 58 loci and highlights GABAergic signaling

“…Seven others (ZNF502, ZNF501, STAB1, NT5DC2, GNL3, GLT8D1, NEK4) are located within an extended region on chromosome 3p21. Other nearby genes include prior reports of association with depression (55), schizophrenia (72), bipolar disorder (73), and suicide (74), making this region a "hot spot" for overall psychiatric susceptibility. A prior GWAS of neuroticism (39) reported association with STAB1 suggesting this gene's dual role in anxiety-related phenotypes.…”