Dataset Information

Genome-Wide Association Study of Suicide Death and Polygenic Prediction of Clinical Antecedents.

ABSTRACT:

Objective

Death by suicide is a highly preventable yet growing worldwide health crisis. To date, there has been a lack of adequately powered genomic studies of suicide, with no sizable suicide death cohorts available for analysis. To address this limitation, the authors conducted the first comprehensive genomic analysis of suicide death using previously unpublished genotype data from a large population-ascertained cohort.

Methods

The analysis sample comprised 3,413 population-ascertained case subjects of European ancestry and 14,810 ancestrally matched control subjects. Analytical methods included principal component analysis for ancestral matching and adjusting for population stratification, linear mixed model genome-wide association testing (conditional on genetic-relatedness matrix), gene and gene set-enrichment testing, and polygenic score analyses, as well as single-nucleotide polymorphism (SNP) heritability and genetic correlation estimation using linkage disequilibrium score regression.

Results

Genome-wide association analysis identified two genome-wide significant loci (involving six SNPs: rs34399104, rs35518298, rs34053895, rs66828456, rs35502061, and rs35256367). Gene-based analyses implicated 22 genes on chromosomes 13, 15, 16, 17, and 19 (q<0.05). Suicide death heritability was estimated at an h²_SNP value of 0.25 (SE=0.04) and a value of 0.16 (SE=0.02) when converted to a liability scale. Notably, suicide polygenic scores were significantly predictive across training and test sets. Polygenic scores for several other psychiatric disorders and psychological traits were also predictive, particularly scores for behavioral disinhibition and major depressive disorder.

Conclusions

Multiple genome-wide significant loci and genes were identified and polygenic score prediction of suicide death case-control status was demonstrated, adjusting for ancestry, in independent training and test sets. Additionally, the suicide death sample was found to have increased genetic risk for behavioral disinhibition, major depressive disorder, depressive symptoms, autism spectrum disorder, psychosis, and alcohol use disorder compared with the control sample.

SUBMITTER: Docherty AR

PROVIDER: S-EPMC7872505 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Multiple sources of evidence suggest that genetic factors influence variation in clinical features of schizophrenia. The authors present the first genome-wide association study (GWAS) of dimensional symptom scores among individuals with schizophrenia.Based on the Lifetime Dimensions of Psychosis Scale ratings of 2,454 case subjects of European ancestry from the Molecular Genetics of Schizophrenia (MGS) sample, three symptom factors (positive, negative/disorganized, and mood) were identified with exploratory factor analysis. Quantitative scores for each factor from a confirmatory factor analysis were analyzed for association with 696,491 single-nucleotide polymorphisms (SNPs) using linear regression, with correction for age, sex, clinical site, and ancestry. Polygenic score analysis was carried out to determine whether case and comparison subjects in 16 Psychiatric GWAS Consortium (PGC) schizophrenia samples (excluding MGS samples) differed in scores computed by weighting their genotypes by MGS association test results for each symptom factor.No genome-wide significant associations were observed between SNPs and factor scores. Most of the SNPs producing the strongest evidence for association were in or near genes involved in neurodevelopment, neuroprotection, or neurotransmission, including genes playing a role in Mendelian CNS diseases, but no statistically significant effect was observed for any defined gene pathway. Finally, polygenic scores based on MGS GWAS results for the negative/disorganized factor were significantly different between case and comparison subjects in the PGC data set; for MGS subjects, negative/disorganized factor scores were correlated with polygenic scores generated using case-control GWAS results from the other PGC samples.The polygenic signal that has been observed in cross-sample analyses of schizophrenia GWAS data sets could be in part related to genetic effects on negative and disorganized symptoms (i.e., core features of chronic schizophrenia).

Project description:Suicidal behavior (SB) has a complex etiology of genes, environment or both. One of the genetic components in SB could be copy number variations (CNVs), since CNVs are implicated in a range of neurodevelopmental disorders. However, a recently published genome-wide and case-control study failed to observe a significant role of CNVs in SB (see E-MTAB-3519). Here we complement those initial observations by conducting a brief CNV-association study, for the first time in a family-based trio-sample with severe suicide attempt (SA) outcome in offspring (n=660 trios; the \GISS\ sample, see http://www.ncbi.nlm.nih.gov/pubmed/23422793 and refs therein). For association testing, we here used the FBAT-CNV methodology (http://www.ncbi.nlm.nih.gov/pubmed/18228561), which allows for CNV association testing directly on the raw intensity values (Illumina \log R ratio\), evaluating any type of CNVs (e.g. de novo or inherited) without reliance on CNV calling algorithms and robust to control selection biases. Here we have deposited these raw logR-values for 88,450 on-chip CNV-loci of the HumanOmni1-Quad_v1 chip, arranged in a standard pedigree format used as input for FBAT-CNV analysis in SVS software (goldenhelix.com): column 1 is the family ID, column 2 is the Subject type (1 = offspring, 2 = mother, 3 = father), column 3 is the father ID (? for missing), column 4 is the mother ID (? for missing), column 5 is the sex (1 = female, 0 = male), column 6 is the affection status (1 = suicide attempt, 0 = not affected) and columns 7 and onward are the logR-values for each CNV-loci as is listed in the header row. We observed experiment-wide significant association (P ≤ 5.6 x 10-7) of two proximal CNVs at chromosome 14 (Illumina markers cnvi0108946 and cnvi0118308, with NCBI 36.3 start positions 22794779 and 22582820). However, these CNV-associations mapped to T-cell receptor regions, and therefore most likely reflect inter-individual variation in somatic rearrangements occurring in white blood cells (as our source of DNA was blood), rather than association with SA. In conclusion, our results did not suggest any major role of CNVs in SA etiology, in line with the results of the recent case-control study (see E-MTAB-3519).