Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program.

Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q's<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (rG = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; rG = 0.78; p = 1.98e-95) and major depressive disorder (MDD; rG = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.

Polygenic risk of social isolation behavior and its influence on psychopathology and personality.

Social isolation has been linked to a range of psychiatric issues, but the behavioral component that drives it is not well understood. Here, a genome-wide associations study (GWAS) was carried out to identify genetic variants that contribute specifically to social isolation behavior (SIB) in up to 449,609 participants from the UK Biobank. 17 loci were identified at genome-wide significance, contributing to a 4% SNP-based heritability estimate. Using the SIB GWAS, polygenic risk scores (PRS) were derived in ALSPAC, an independent, developmental cohort, and used to test for association with self-reported friendship scores, comprising items related to friendship quality and quantity, at age 12 and 18 to determine whether genetic predisposition manifests during childhood development. At age 18, friendship scores were associated with the SIB PRS, demonstrating that the genetic factors can predict related social traits in late adolescence. Linkage disequilibrium (LD) score correlation using the SIB GWAS demonstrated genetic correlations with autism spectrum disorder (ASD), schizophrenia, major depressive disorder (MDD), educational attainment, extraversion, and loneliness. However, no evidence of causality was found using a conservative Mendelian randomization approach between SIB and any of the traits in either direction. Genomic Structural Equation Modeling (SEM) revealed a common factor contributing to SIB, neuroticism, loneliness, MDD, and ASD, weakly correlated with a second common factor that contributes to psychiatric and psychotic traits. Our results show that SIB contributes a small heritable component, which is associated genetically with other social traits such as friendship as well as psychiatric disorders.

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

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.

Genomics-driven screening for causal determinants of suicide attempt.

OBJECTIVE: Each year, around one million people die by suicide. Despite its recognition as a public health concern, large-scale research on causal determinants of suicide attempt risk is scarce. Here, we leverage results from a recent genome-wide association study (GWAS) of suicide attempt to perform a data-driven screening of traits causally associated with suicide attempt. METHODS: We performed a hypothesis-generating phenome-wide screening of causal relationships between suicide attempt risk and 1520 traits, which have been systematically aggregated on the Complex-Traits Genomics Virtual Lab platform. We employed the latent causal variable (LCV) method, which uses results from GWAS to assess whether a causal relationship can explain a genetic correlation between two traits. If a trait causally influences another one, the genetic variants that increase risk for the causal trait will also increase the risk for the outcome inducing a genetic correlation. Nonetheless, a genetic correlation can also be observed when traits share common pathways. The LCV method can assess whether the pattern of genetic effects for two genetically correlated traits support a causal association rather than a shared aetiology. RESULTS: Our approach identified 62 traits that increased risk for suicide attempt. Risk factors identified can be broadly classified into (1) physical health disorders, including oesophagitis, fibromyalgia, hernia and cancer; (2) mental health-related traits, such as depression, substance use disorders and anxiety; and (3) lifestyle traits including being involved in combat or exposure to a war zone, and specific job categories such as being a truck driver or machine operator. CONCLUSIONS: Suicide attempt risk is likely explained by a combination of behavioural phenotypes and risk for both physical and psychiatric disorders. Our results also suggest that substance use behaviours and pain-related conditions are associated with an increased suicide attempt risk, elucidating important causal mechanisms that underpin this significant public health problem.

Investigating rare pathogenic/likely pathogenic exonic variation in bipolar disorder.

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3-2.8, one-sided p = 6.0 × 10-4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ.

RICOPILI: Rapid Imputation for COnsortias PIpeLIne.

SUMMARY: Genome-wide association study (GWAS) analyses, at sufficient sample sizes and power, have successfully revealed biological insights for several complex traits. RICOPILI, an open-sourced Perl-based pipeline was developed to address the challenges of rapidly processing large-scale multi-cohort GWAS studies including quality control (QC), imputation and downstream analyses. The pipeline is computationally efficient with portability to a wide range of high-performance computing environments. RICOPILI was created as the Psychiatric Genomics Consortium pipeline for GWAS and adopted by other users. The pipeline features (i) technical and genomic QC in case-control and trio cohorts, (ii) genome-wide phasing and imputation, (iv) association analysis, (v) meta-analysis, (vi) polygenic risk scoring and (vii) replication analysis. Notably, a major differentiator from other GWAS pipelines, RICOPILI leverages on automated parallelization and cluster job management approaches for rapid production of imputed genome-wide data. A comprehensive meta-analysis of simulated GWAS data has been incorporated demonstrating each step of the pipeline. This includes all the associated visualization plots, to allow ease of data interpretation and manuscript preparation. Simulated GWAS datasets are also packaged with the pipeline for user training tutorials and developer work. AVAILABILITY AND IMPLEMENTATION: RICOPILI has a flexible architecture to allow for ongoing development and incorporation of newer available algorithms and is adaptable to various HPC environments (QSUB, BSUB, SLURM and others). Specific links for genomic resources are either directly provided in this paper or via tutorials and external links. The central location hosting scripts and tutorials is found at this URL: https://sites.google.com/a/broadinstitute.org/RICOPILI/home. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

DNA methylation and inflammation marker profiles associated with a history of depression.

Depression is a common and disabling disorder, representing a major social and economic health issue. Moreover, depression is associated with the progression of diseases with an inflammatory etiology including many inflammatory-related disorders. At the molecular level, the mechanisms by which depression might promote the onset of these diseases and associated immune-dysfunction are not well understood. In this study we assessed genome-wide patterns of DNA methylation in whole blood-derived DNA obtained from individuals with a self-reported history of depression (n = 100) and individuals without a history of depression (n = 100) using the Illumina 450K microarray. Our analysis identified six significant (Sidák corrected P < 0.05) depression-associated differentially methylated regions (DMRs); the top-ranked DMR was located in exon 1 of the LTB4R2 gene (Sidák corrected P = 1.27 × 10-14). Polygenic risk scores (PRS) for depression were generated and known biological markers of inflammation, telomere length (TL) and IL-6, were measured in DNA and serum samples, respectively. Next, we employed a systems-level approach to identify networks of co-methylated loci associated with a history of depression, in addition to depression PRS, TL and IL-6 levels. Our analysis identified one depression-associated co-methylation module (P = 0.04). Interestingly, the depression-associated module was highly enriched for pathways related to immune function and was also associated with TL and IL-6 cytokine levels. In summary, our genome-wide DNA methylation analysis of individuals with and without a self-reported history of depression identified several candidate DMRs of potential relevance to the pathogenesis of depression and its associated immune-dysfunction phenotype.

Genome-wide gene-environment interaction in depression: A systematic evaluation of candidate genes: The childhood trauma working-group of PGC-MDD.

Gene by environment (GxE) interaction studies have investigated the influence of a number of candidate genes and variants for major depressive disorder (MDD) on the association between childhood trauma and MDD. Most of these studies are hypothesis driven and investigate only a limited number of SNPs in relevant pathways using differing methodological approaches. Here (1) we identified 27 genes and 268 SNPs previously associated with MDD or with GxE interaction in MDD and (2) analyzed their impact on GxE in MDD using a common approach in 3944 subjects of European ancestry from the Psychiatric Genomics Consortium who had completed the Childhood Trauma Questionnaire. (3) We subsequently used the genome-wide SNP data for a genome-wide case-control GxE model and GxE case-only analyses testing for an enrichment of associated SNPs. No genome-wide significant hits and no consistency among the signals of the different analytic approaches could be observed. This is the largest study for systematic GxE interaction analysis in MDD in subjects of European ancestry to date. Most of the known candidate genes/variants could not be supported. Thus, their impact on GxE interaction in MDD may be questionable. Our results underscore the need for larger samples, more extensive assessment of environmental exposures, and greater efforts to investigate new methodological approaches in GxE models for MDD.

Genetics of Depression: Progress at Last.

PURPOSE OF REVIEW: We will describe the success of recent genome-wide association studies that identify genetic variants associated with depression and outline the strategies used to reduce heterogeneity and increase sample size. RECENT FINDINGS: The CONVERGE consortium identified two genetic associations by focusing on a sample of Chinese women with recurrent severe depression. Three other loci have been found in Europeans by combining cohorts with clinical diagnosis and measures of depressive symptoms to increase sample size. 23andMe identified 15 loci associated with depression using self-report of clinical diagnosis in a study of over 300,000 individuals. The first genetic associations with depression have been identified, and this number is now expected to increase linearly with sample size, as seen in other polygenic disorders. These loci provide invaluable insights into the biology of depression and exciting opportunities to develop new biomarkers and therapeutic targets.

Genetic relationships between suicide attempts, suicidal ideation and major psychiatric disorders: a genome-wide association and polygenic scoring study.

Epidemiological studies have recognized a genetic diathesis for suicidal behavior, which is independent of other psychiatric disorders. Genome-wide association studies (GWAS) on suicide attempt (SA) and ideation have failed to identify specific genetic variants. Here, we conduct further GWAS and for the first time, use polygenic score analysis in cohorts of patients with mood disorders, to test for common genetic variants for mood disorders and suicide phenotypes. Genome-wide studies for SA were conducted in the RADIANT and GSK-Munich recurrent depression samples and London Bipolar Affective Disorder Case-Control Study (BACCs) then meta-analysis was performed. A GWAS on suicidal ideation during antidepressant treatment had previously been conducted in the Genome Based Therapeutic Drugs for Depression (GENDEP) study. We derived polygenic scores from each sample and tested their ability to predict SA in the mood disorder cohorts or ideation status in the GENDEP study. Polygenic scores for major depressive disorder, bipolar disorder and schizophrenia from the Psychiatric Genomics Consortium were used to investigate pleiotropy between psychiatric disorders and suicide phenotypes. No significant evidence for association was detected at any SNP in GWAS or meta-analysis. Polygenic scores for major depressive disorder significantly predicted suicidal ideation in the GENDEP pharmacogenetics study and also predicted SA in a combined validation dataset. Polygenic scores for SA showed no predictive ability for suicidal ideation. Polygenic score analysis suggests pleiotropy between psychiatric disorders and suicidal ideation whereas the tendency to act on such thoughts may have a partially independent genetic diathesis.

Absence of nonfatal suicidal behavior preceding suicide death reveals differences in clinical risks.

Nonfatal suicidality is the most robust predictor of suicide death. However, only ~10% of those who survive an attempt go on to die by suicide. Moreover, ~50% of suicide deaths occur in the absence of prior known attempts, suggesting risks other than nonfatal suicide attempt need to be identified. We studied data from 4,000 population-ascertained suicide deaths and 26,191 population controls to improve understanding of risks leading to suicide death. This study included 2,253 suicide deaths and 3,375 controls with evidence of nonfatal suicidality (SUI_SI/SB and CTL_SI/SB) from diagnostic codes and natural language processing of electronic health records notes. Characteristics of these groups were compared to 1,669 suicides with no prior nonfatal SI/SB (SUI_None) and 22,816 controls with no lifetime suicidality (CTL_None). The SUI_None and CTL_None groups had fewer diagnoses and were older than SUI_SI/SB and CTL_SI/SB. Mental health diagnoses were far less common in both the SUI_None and CTL_None groups; mental health problems were less associated with suicide death than with presence of SI/SB. Physical health diagnoses were conversely more often associated with risk of suicide death than with presence of SI/SB. Pending replication, results indicate highly significant clinical differences among suicide deaths with versus without prior nonfatal SI/SB.

Analyses of GWAS signal using GRIN identify additional genes contributing to suicidal behavior.

Genome-wide association studies (GWAS) identify genetic variants underlying complex traits but are limited by stringent genome-wide significance thresholds. We present GRIN (Gene set Refinement through Interacting Networks), which increases confidence in the expanded gene set by retaining genes strongly connected by biological networks when GWAS thresholds are relaxed. GRIN was validated on both simulated interrelated gene sets as well as multiple GWAS traits. From multiple GWAS summary statistics of suicide attempt, a complex phenotype, GRIN identified additional genes that replicated across independent cohorts and retained biologically interrelated genes despite a relaxed significance threshold. We present a conceptual model of how these retained genes interact through neurobiological pathways that may influence suicidal behavior, and identify existing drugs associated with these pathways that would not have been identified under traditional GWAS thresholds. We demonstrate GRIN's utility in boosting GWAS results by increasing the number of true positive genes identified from GWAS results.

Clinical, genomic, and neurophysiological correlates of lifetime suicide attempts among individuals with alcohol dependence.

Research has identified clinical, genomic, and neurophysiological markers associated with suicide attempts (SA) among individuals with psychiatric illness. However, there is limited research among those with an alcohol use disorder (AUD), despite their disproportionately higher rates of SA. We examined lifetime SA in 4,068 individuals with DSM-IV alcohol dependence from the Collaborative Study on the Genetics of Alcoholism (23% lifetime suicide attempt; 53% female; 17% Admixed African American ancestries; mean age: 38). We 1) conducted a genome-wide association study (GWAS) of SA and performed downstream analyses to determine whether we could identify specific biological pathways of risk, and 2) explored risk in aggregate across other clinical conditions, polygenic scores (PGS) for comorbid psychiatric problems, and neurocognitive functioning between those with AD who have and have not reported a lifetime suicide attempt. The GWAS and downstream analyses did not produce any significant associations. Participants with an AUD who had attempted suicide had greater rates of trauma exposure, major depressive disorder, post-traumatic stress disorder, and other substance use disorders compared to those who had not attempted suicide. Polygenic scores for suicide attempt, depression, and PTSD were associated with reporting a suicide attempt (ORs = 1.22-1.44). Participants who reported a SA also had decreased right hemispheric frontal-parietal theta and decreased interhemispheric temporal-parietal alpha electroencephalogram resting-state coherences relative to those who did not, but differences were small. Overall, individuals with alcohol dependence who report SA appear to experience a variety of severe comorbidities and elevated polygenic risk for SA. Our results demonstrate the need to further investigate suicide attempts in the presence of substance use disorders.

Distinguishing different psychiatric disorders using DDx-PRS.

Despite great progress on methods for case-control polygenic prediction (e.g. schizophrenia vs. control), there remains an unmet need for a method that genetically distinguishes clinically related disorders (e.g. schizophrenia (SCZ) vs. bipolar disorder (BIP) vs. depression (MDD) vs. control); such a method could have important clinical value, especially at disorder onset when differential diagnosis can be challenging. Here, we introduce a method, Differential Diagnosis-Polygenic Risk Score (DDx-PRS), that jointly estimates posterior probabilities of each possible diagnostic category (e.g. SCZ=50%, BIP=25%, MDD=15%, control=10%) by modeling variance/covariance structure across disorders, leveraging case-control polygenic risk scores (PRS) for each disorder (computed using existing methods) and prior clinical probabilities for each diagnostic category. DDx-PRS uses only summary-level training data and does not use tuning data, facilitating implementation in clinical settings. In simulations, DDx-PRS was well-calibrated (whereas a simpler approach that analyzes each disorder marginally was poorly calibrated), and effective in distinguishing each diagnostic category vs. the rest. We then applied DDx-PRS to Psychiatric Genomics Consortium SCZ/BIP/MDD/control data, including summary-level training data from 3 case-control GWAS ( N =41,917-173,140 cases; total N =1,048,683) and held-out test data from different cohorts with equal numbers of each diagnostic category (total N =11,460). DDx-PRS was well-calibrated and well-powered relative to these training sample sizes, attaining AUCs of 0.66 for SCZ vs. rest, 0.64 for BIP vs. rest, 0.59 for MDD vs. rest, and 0.68 for control vs. rest. DDx-PRS produced comparable results to methods that leverage tuning data, confirming that DDx-PRS is an effective method. True diagnosis probabilities in top deciles of predicted diagnosis probabilities were considerably larger than prior baseline probabilities, particularly in projections to larger training sample sizes, implying considerable potential for clinical utility under certain circumstances. In conclusion, DDx-PRS is an effective method for distinguishing clinically related disorders.

Defining Suicidal Thought and Behavior Phenotypes for Genetic Studies.

Background: Standardized definitions of suicidality phenotypes, including suicidal ideation (SI), attempt (SA), and death (SD) are a critical step towards improving understanding and comparison of results in suicide research. The complexity of suicidality contributes to heterogeneity in phenotype definitions, impeding evaluation of clinical and genetic risk factors across studies and efforts to combine samples within consortia. Here, we present expert and data-supported recommendations for defining suicidality and control phenotypes to facilitate merging current/legacy samples with definition variability and aid future sample creation. Methods: A subgroup of clinician researchers and experts from the Suicide Workgroup of the Psychiatric Genomics Consortium (PGC) reviewed existing PGC definitions for SI, SA, SD, and control groups and generated preliminary consensus guidelines for instrument-derived and international classification of disease (ICD) data. ICD lists were validated in two independent datasets (N = 9,151 and 12,394). Results: Recommendations are provided for evaluated instruments for SA and SI, emphasizing selection of lifetime measures phenotype-specific wording. Recommendations are also provided for defining SI and SD from ICD data. As the SA ICD definition is complex, SA code list recommendations were validated against instrument results with sensitivity (range = 15.4% to 80.6%), specificity (range = 67.6% to 97.4%), and positive predictive values (range = 0.59-0.93) reported. Conclusions: Best-practice guidelines are presented for the use of existing information to define SI/SA/SD in consortia research. These proposed definitions are expected to facilitate more homogeneous data aggregation for genetic and multisite studies. Future research should involve refinement, improved generalizability, and validation in diverse populations.

Polygenic risk of Social-isolation and its influence on social behavior, psychosis, depression and autism spectrum disorder.

Social-isolation has been linked to a range of psychiatric issues, but the behavioral component that drives it is not well understood. Here, a GWAS is carried out to identify genetic variants which contribute to Social-isolation behaviors in up to 449,609 participants from the UK Biobank. 17 loci were identified at genome-wide significance, contributing to a 4% SNP heritability estimate. Using the Social-isolation GWAS, polygenic risk scores (PRS) were derived in ALSPAC, an independent, developmental cohort, and used to test for association with friendship quality. At age 18, friendship scores were associated with the Social-isolation PRS, demonstrating that the genetic factors are able to predict related social traits. LD score regression using the GWAS demonstrated genetic correlation with autism spectrum disorder, schizophrenia, and major depressive disorder. However, no evidence of causality was found using a conservative Mendelian randomization approach other than that of autism spectrum disorder on Social-isolation. Our results show that Social-isolation has a small heritable component which may drive those behaviors which is associated genetically with other social traits such as friendship satisfaction as well as psychiatric disorders.

Polygenic Contributions to Suicidal Thoughts and Behaviors in a Sample Ascertained for Alcohol Use Disorders.

Introduction: Suicidal thoughts and behaviors have partially distinct genetic etiologies. Methods: We used PRS-CS to create polygenic risk scores (PRSs) from GWAS of non-suicidal self-injury, broad-sense self-harm ideation, nonfatal suicide attempt, death by suicide, and depression. Using mixed-effect models, we estimated whether these PRSs were associated with a range of suicidal thoughts and behaviors in the Collaborative Study on the Genetics of Alcoholism (N = 7,526). Results: All PRSs were significantly associated with suicidal ideation and suicide attempt (betas = 0.08-0.44, false discovery rate [FDR] <0.023). All PRSs except non-suicidal self-injury PRS were associated with active suicidal ideation (betas = 0.14-0.22, FDR <0.003). Several associations remained significant in models where all significant PRSs were included as simultaneous predictors, and when all PRSs predicted suicide attempt, the PRS together explained 6.2% of the variance in suicide attempt. Significant associations were also observed between some PRSs and persistent suicidal ideation, non-suicidal self-injury, compounded suicide attempt, and desire to die. Conclusion: Our findings suggest that PRS for depression does not explain the entirety of the variance in suicidal thoughts and behaviors, with PRS specifically for suicidal thoughts and behaviors making additional and sometimes unique contributions.

Identification of Novel, Replicable Genetic Risk Loci for Suicidal Thoughts and Behaviors Among US Military Veterans.

Importance: Suicide is a leading cause of death; however, the molecular genetic basis of suicidal thoughts and behaviors (SITB) remains unknown. Objective: To identify novel, replicable genomic risk loci for SITB. Design, Setting, and Participants: This genome-wide association study included 633 778 US military veterans with and without SITB, as identified through electronic health records. GWAS was performed separately by ancestry, controlling for sex, age, and genetic substructure. Cross-ancestry risk loci were identified through meta-analysis. Study enrollment began in 2011 and is ongoing. Data were analyzed from November 2021 to August 2022. Main Outcome and Measures: SITB. Results: A total of 633 778 US military veterans were included in the analysis (57 152 [9%] female; 121 118 [19.1%] African ancestry, 8285 [1.3%] Asian ancestry, 452 767 [71.4%] European ancestry, and 51 608 [8.1%] Hispanic ancestry), including 121 211 individuals with SITB (19.1%). Meta-analysis identified more than 200 GWS (P < 5 × 10-8) cross-ancestry risk single-nucleotide variants for SITB concentrated in 7 regions on chromosomes 2, 6, 9, 11, 14, 16, and 18. Top single-nucleotide variants were largely intronic in nature; 5 were independently replicated in ISGC, including rs6557168 in ESR1, rs12808482 in DRD2, rs77641763 in EXD3, rs10671545 in DCC, and rs36006172 in TRAF3. Associations for FBXL19 and AC018880.2 were not replicated. Gene-based analyses implicated 24 additional GWS cross-ancestry risk genes, including FURIN, TSNARE1, and the NCAM1-TTC12-ANKK1-DRD2 gene cluster. Cross-ancestry enrichment analyses revealed significant enrichment for expression in brain and pituitary tissue, synapse and ubiquitination processes, amphetamine addiction, parathyroid hormone synthesis, axon guidance, and dopaminergic pathways. Seven other unique European ancestry-specific GWS loci were identified, 2 of which (POM121L2 and METTL15/LINC02758) were replicated. Two additional GWS ancestry-specific loci were identified within the African ancestry (PET112/GATB) and Hispanic ancestry (intergenic locus on chromosome 4) subsets, both of which were replicated. No GWS loci were identified within the Asian ancestry subset; however, significant enrichment was observed for axon guidance, cyclic adenosine monophosphate signaling, focal adhesion, glutamatergic synapse, and oxytocin signaling pathways across all ancestries. Within the European ancestry subset, genetic correlations (r > 0.75) were observed between the SITB phenotype and a suicide attempt-only phenotype, depression, and posttraumatic stress disorder. Additionally, polygenic risk score analyses revealed that the Million Veteran Program polygenic risk score had nominally significant main effects in 2 independent samples of veterans of European and African ancestry. Conclusions and Relevance: The findings of this analysis may advance understanding of the molecular genetic basis of SITB and provide evidence for ESR1, DRD2, TRAF3, and DCC as cross-ancestry candidate risk genes. More work is needed to replicate these findings and to determine if and how these genes might impact clinical care.

Polygenic Scores and Onset of Major Mood or Psychotic Disorders Among Offspring of Affected Parents.

OBJECTIVE: Family history is an established risk factor for mental illness. The authors sought to investigate whether polygenic scores (PGSs) can complement family history to improve identification of risk for major mood and psychotic disorders. METHODS: Eight cohorts were combined to create a sample of 1,884 participants ages 2-36 years, including 1,339 offspring of parents with mood or psychotic disorders, who were prospectively assessed with diagnostic interviews over an average of 5.1 years. PGSs were constructed for depression, bipolar disorder, anxiety, attention deficit hyperactivity disorder (ADHD), schizophrenia, neuroticism, subjective well-being, p factor, and height (as a negative control). Cox regression was used to test associations between PGSs, family history of major mental illness, and onsets of major mood and psychotic disorders. RESULTS: There were 435 onsets of major mood and psychotic disorders across follow-up. PGSs for neuroticism (hazard ratio=1.23, 95% CI=1.12-1.36), schizophrenia (hazard ratio=1.15, 95% CI=1.04-1.26), depression (hazard ratio=1.11, 95% CI=1.01-1.22), ADHD (hazard ratio=1.10, 95% CI=1.00-1.21), subjective well-being (hazard ratio=0.90, 95% CI=0.82-0.99), and p factor (hazard ratio=1.14, 95% CI=1.04-1.26) were associated with onsets. After controlling for family history, neuroticism PGS remained significantly positively associated (hazard ratio=1.19, 95% CI=1.08-1.31) and subjective well-being PGS remained significantly negatively associated (hazard ratio=0.89, 95% CI=0.81-0.98) with onsets. CONCLUSIONS: Neuroticism and subjective well-being PGSs capture risk of major mood and psychotic disorders that is independent of family history, whereas PGSs for psychiatric illness provide limited predictive power when family history is known. Neuroticism and subjective well-being PGSs may complement family history in the early identification of persons at elevated risk.