Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder

Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed health problems, prior pharmacological treatments, and polygenic scores (PGS) has potential to inform risk stratification. We examined self-reported SB and ideation using the Columbia Suicide Severity Rating Scale (C-SSRS) among 3,942 SCZ and 5,414 BPI patients receiving care within the Veterans Health Administration (VHA). These cross-sectional data were integrated with electronic health records (EHRs), and compared across lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. PGS were constructed using available genomic data for related traits. Genome-wide association studies were performed to identify and prioritize specific loci. Only 20% of the veterans who reported SB had a corroborating ICD-9/10 EHR code. Among those without prior SB, more than 20% reported new-onset SB at follow-up. SB were associated with a range of additional clinical diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking initiation, suicide attempt, and major depressive disorder were associated with SB. The GWAS for SB yielded no significant loci. Among individuals with a diagnosed mental illness, self-reported SB were strongly associated with clinical variables across several EHR domains. Analyses point to sequelae of substance-related and psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in health records, underscoring the value of regular screening with direct, in-person assessments, especially among high-risk individuals.

Correlates of suicidal behaviors and genetic risk among United States veterans with schizophrenia or bipolar I disorder

Objective

Persons diagnosed with schizophrenia (SCZ) or bipolar I disorder (BPI) are at high risk for self-injurious behavior, suicidal ideation, and suicidal behaviors (SB). Characterizing associations between diagnosed mental and physical health problems, prior pharmacological treatments, and aggregate genetic factors has potential to inform risk stratification and mitigation strategies.

Methods

In this study of 3,942 SCZ and 5,414 BPI patients receiving VA care, self-reported SB and ideation were assessed using the Columbia Suicide Severity Rating Scale (C-SSRS). These cross-sectional data were integrated with electronic health records (EHR), and compared by lifetime diagnoses, treatment histories, follow-up screenings, and mortality data. Polygenic scores (PGS) for traits related to psychiatric disorders, substance use, and cognition were constructed using available genomic data, and exploratory genome-wide association studies were performed to identify and prioritize specific loci.

Results

Only 20% of veterans who self-reported SB had a corroborating ICD-9/10 code in their EHR; and among those who denied prior behaviors, more than 20% reported new-onset SB at follow-up. SB were associated with a range of psychiatric and non-psychiatric diagnoses, and with treatment with specific classes of psychotropic medications (e.g., antidepressants, antipsychotics, etc.). PGS for externalizing behaviors, smoking, suicide attempt, and major depressive disorder were also associated with attempt and ideation.

Conclusions

Among individuals with a diagnosed mental illness, a GWAS for SB did not yield any significant loci. Self-reported SB were strongly associated with clinical variables across several EHR domains. Overall, clinical and polygenic analyses point to sequelae of substance-use related behaviors and other psychiatric comorbidities as strong correlates of prior and subsequent SB. Nonetheless, past SB was frequently not documented in clinical settings, underscoring the value of regular screening based on direct, in-person assessments, especially among high-risk individuals.

Investigation of convergent and divergent genetic influences underlying schizophrenia and alcohol use disorder

BACKGROUND

Alcohol use disorder (AUD) and schizophrenia (SCZ) frequently co-occur, and large-scale genome-wide association studies (GWAS) have identified significant genetic correlations between these disorders.

METHODS

We used the largest published GWAS for AUD (total cases = 77 822) and SCZ (total cases = 46 827) to identify genetic variants that influence both disorders (with either the same or opposite direction of effect) and those that are disorder specific.

RESULTS

We identified 55 independent genome-wide significant single nucleotide polymorphisms with the same direction of effect on AUD and SCZ, 8 with robust effects in opposite directions, and 98 with disorder-specific effects. We also found evidence for 12 genes whose pleiotropic associations with AUD and SCZ are consistent with mediation via gene expression in the prefrontal cortex. The genetic covariance between AUD and SCZ was concentrated in genomic regions functional in brain tissues (p = 0.001).

CONCLUSIONS

Our findings provide further evidence that SCZ shares meaningful genetic overlap with AUD.

Detecting and Adjusting for Hidden Biases due to Phenotype Misclassification in Genome-Wide Association Studies

With the advent of healthcare-based genotyped biobanks, genome-wide association studies (GWAS) leverage larger sample sizes, incorporate patients with diverse ancestries and introduce noisier phenotypic definitions. Yet the extent and impact of phenotypic misclassification on large-scale datasets is not currently well understood due to a lack of statistical methods to estimate relevant parameters from empirical data. Here, we develop a statistical method and scalable software, PheMED, Phenotypic Measurement of Effective Dilution, to quantify phenotypic misclassification across GWAS using only summary statistics. We illustrate how the parameters estimated by PheMED relate to the negative and positive predictive value of the labeled phenotype, compared to ground truth, and how misclassification of the phenotype yields diluted effect-sizes of variant-phenotype associations. Furthermore, we apply our methodology to detect multiple instances of statistically significant dilution in real-world data. We demonstrate how effective dilution biases downstream GWAS replication and heritability analyses despite utilizing current best practices, and provide a dilution-aware meta-analysis approach that outperforms existing methods. Consequently, we anticipate that PheMED will be a valuable tool for researchers to address phenotypic data quality issues both within and across cohorts.

Penetrance and Pleiotropy of Polygenic Risk Scores for Schizophrenia, Bipolar Disorder, and Depression Among Adults in the US Veterans Affairs Health Care System

Importance

Serious mental illnesses, including schizophrenia, bipolar disorder, and depression, are heritable, highly multifactorial disorders and major causes of disability worldwide.

Objective

To benchmark the penetrance of current neuropsychiatric polygenic risk scores (PRSs) in the Veterans Health Administration health care system and to explore associations between PRS and broad categories of human disease via phenome-wide association studies.

Design, Setting, and Participants

Extensive Veterans Health Administration's electronic health records were assessed from October 1999 to January 2021, and an embedded cohort of 9378 individuals with confirmed diagnoses of schizophrenia or bipolar 1 disorder were found. The performance of schizophrenia, bipolar disorder, and major depression PRSs were compared in participants of African or European ancestry in the Million Veteran Program (approximately 400 000 individuals), and associations between PRSs and 1650 disease categories based on ICD-9/10 billing codes were explored. Last, genomic structural equation modeling was applied to derive novel PRSs indexing common and disorder-specific genetic factors. Analysis took place from January 2021 to January 2022.

Main Outcomes and Measures

Diagnoses based on in-person structured clinical interviews were compared with ICD-9/10 billing codes. PRSs were constructed using summary statistics from genome-wide association studies of schizophrenia, bipolar disorder, and major depression.

Results

Of 707 299 enrolled study participants, 459 667 were genotyped at the time of writing; 84 806 were of broadly African ancestry (mean [SD] age, 58 [12.1] years) and 314 909 were of broadly European ancestry (mean [SD] age, 66.4 [13.5] years). Among 9378 individuals with confirmed diagnoses of schizophrenia or bipolar 1 disorder, 8962 (95.6%) were correctly identified using ICD-9/10 codes (2 or more). Among those of European ancestry, PRSs were robustly associated with having received a diagnosis of schizophrenia (odds ratio [OR], 1.81 [95% CI, 1.76-1.87]; P < 10-257) or bipolar disorder (OR, 1.42 [95% CI, 1.39-1.44]; P < 10-295). Corresponding effect sizes in participants of African ancestry were considerably smaller for schizophrenia (OR, 1.35 [95% CI, 1.29-1.42]; P < 10-38) and bipolar disorder (OR, 1.16 [95% CI, 1.11-1.12]; P < 10-10). Neuropsychiatric PRSs were associated with increased risk for a range of psychiatric and physical health problems.

Conclusions and Relevance

Using diagnoses confirmed by in-person structured clinical interviews and current neuropsychiatric PRSs, the validity of an electronic health records-based phenotyping approach in US veterans was demonstrated, highlighting the potential of PRSs for disentangling biological and mediated pleiotropy.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Rare coding variants in ten genes confer substantial risk for schizophrenia

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

Interaction Testing and Polygenic Risk Scoring to Estimate the Association of Common Genetic Variants With Treatment Resistance in Schizophrenia

IMPORTANCE

About 20% to 30% of people with schizophrenia have psychotic symptoms that do not respond adequately to first-line antipsychotic treatment. This clinical presentation, chronic and highly disabling, is known as treatment-resistant schizophrenia (TRS). The causes of treatment resistance and their relationships with causes underlying schizophrenia are largely unknown. Adequately powered genetic studies of TRS are scarce because of the difficulty in collecting data from well-characterized TRS cohorts.

OBJECTIVE

To examine the genetic architecture of TRS through the reassessment of genetic data from schizophrenia studies and its validation in carefully ascertained clinical samples.

DESIGN, SETTING, AND PARTICIPANTS

Two case-control genome-wide association studies (GWASs) of schizophrenia were performed in which the case samples were defined as individuals with TRS (n = 10 501) and individuals with non-TRS (n = 20 325). The differences in effect sizes for allelic associations were then determined between both studies, the reasoning being such differences reflect treatment resistance instead of schizophrenia. Genotype data were retrieved from the CLOZUK and Psychiatric Genomics Consortium (PGC) schizophrenia studies. The output was validated using polygenic risk score (PRS) profiling of 2 independent schizophrenia cohorts with TRS and non-TRS: a prevalence sample with 817 individuals (Cardiff Cognition in Schizophrenia [CardiffCOGS]) and an incidence sample with 563 individuals (Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances [STRATA-G]).

MAIN OUTCOMES AND MEASURES

GWAS of treatment resistance in schizophrenia. The results of the GWAS were compared with complex polygenic traits through a genetic correlation approach and were used for PRS analysis on the independent validation cohorts using the same TRS definition.

RESULTS

The study included a total of 85 490 participants (48 635 [56.9%] male) in its GWAS stage and 1380 participants (859 [62.2%] male) in its PRS validation stage. Treatment resistance in schizophrenia emerged as a polygenic trait with detectable heritability (1% to 4%), and several traits related to intelligence and cognition were found to be genetically correlated with it (genetic correlation, 0.41-0.69). PRS analysis in the CardiffCOGS prevalence sample showed a positive association between TRS and a history of taking clozapine (r2 = 2.03%; P = .001), which was replicated in the STRATA-G incidence sample (r2 = 1.09%; P = .04).

CONCLUSIONS AND RELEVANCE

In this GWAS, common genetic variants were differentially associated with TRS, and these associations may have been obscured through the amalgamation of large GWAS samples in previous studies of broadly defined schizophrenia. Findings of this study suggest the validity of meta-analytic approaches for studies on patient outcomes, including treatment resistance.

Dissecting the Shared Genetic Architecture of Suicide Attempt, Psychiatric Disorders, and Known Risk Factors

BACKGROUND

Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders.

METHODS

We conducted a genome-wide association study (GWAS) of 29,782 suicide attempt (SA) cases and 519,961 controls in the International Suicide Genetics Consortium (ISGC). The GWAS of SA was conditioned on psychiatric disorders using GWAS summary statistics via multitrait-based conditional and joint analysis, to remove genetic effects on SA mediated by psychiatric disorders. We investigated the shared and divergent genetic architectures of SA, psychiatric disorders, and other known risk factors.

RESULTS

Two loci reached genome-wide significance for SA: the major histocompatibility complex and an intergenic locus on chromosome 7, the latter of which remained associated with SA after conditioning on psychiatric disorders and replicated in an independent cohort from the Million Veteran Program. This locus has been implicated in risk-taking behavior, smoking, and insomnia. SA showed strong genetic correlation with psychiatric disorders, particularly major depression, and also with smoking, pain, risk-taking behavior, sleep disturbances, lower educational attainment, reproductive traits, lower socioeconomic status, and poorer general health. After conditioning on psychiatric disorders, the genetic correlations between SA and psychiatric disorders decreased, whereas those with nonpsychiatric traits remained largely unchanged.

CONCLUSIONS

Our results identify a risk locus that contributes more strongly to SA than other phenotypes and suggest a shared underlying biology between SA and known risk factors that is not mediated by psychiatric disorders.

Characterisation of age and polarity at onset in bipolar disorder

BACKGROUND

Studying phenotypic and genetic characteristics of age at onset (AAO) and polarity at onset (PAO) in bipolar disorder can provide new insights into disease pathology and facilitate the development of screening tools.

AIMS

To examine the genetic architecture of AAO and PAO and their association with bipolar disorder disease characteristics.

METHOD

Genome-wide association studies (GWASs) and polygenic score (PGS) analyses of AAO (n = 12 977) and PAO (n = 6773) were conducted in patients with bipolar disorder from 34 cohorts and a replication sample (n = 2237). The association of onset with disease characteristics was investigated in two of these cohorts.

RESULTS

Earlier AAO was associated with a higher probability of psychotic symptoms, suicidality, lower educational attainment, not living together and fewer episodes. Depressive onset correlated with suicidality and manic onset correlated with delusions and manic episodes. Systematic differences in AAO between cohorts and continents of origin were observed. This was also reflected in single-nucleotide variant-based heritability estimates, with higher heritabilities for stricter onset definitions. Increased PGS for autism spectrum disorder (β = -0.34 years, s.e. = 0.08), major depression (β = -0.34 years, s.e. = 0.08), schizophrenia (β = -0.39 years, s.e. = 0.08), and educational attainment (β = -0.31 years, s.e. = 0.08) were associated with an earlier AAO. The AAO GWAS identified one significant locus, but this finding did not replicate. Neither GWAS nor PGS analyses yielded significant associations with PAO.

CONCLUSIONS

AAO and PAO are associated with indicators of bipolar disorder severity. Individuals with an earlier onset show an increased polygenic liability for a broad spectrum of psychiatric traits. Systematic differences in AAO across cohorts, continents and phenotype definitions introduce significant heterogeneity, affecting analyses.

Prognostic value of polygenic risk scores for adults with psychosis

Polygenic risk scores (PRS) summarize genetic liability to a disease at the individual level, and the aim is to use them as biomarkers of disease and poor outcomes in real-world clinical practice. To date, few studies have assessed the prognostic value of PRS relative to standards of care. Schizophrenia (SCZ), the archetypal psychotic illness, is an ideal test case for this because the predictive power of the SCZ PRS exceeds that of most other common diseases. Here, we analyzed clinical and genetic data from two multi-ethnic cohorts totaling 8,541 adults with SCZ and related psychotic disorders, to assess whether the SCZ PRS improves the prediction of poor outcomes relative to clinical features captured in a standard psychiatric interview. For all outcomes investigated, the SCZ PRS did not improve the performance of predictive models, an observation that was generally robust to divergent case ascertainment strategies and the ancestral background of the study participants.

Cooperative Studies Program (CSP) #572: A Study of Serious Mental Illness in Veterans as a Pathway to personalized medicine in Schizophrenia and Bipolar Illness

BACKGROUND

Personalization of psychiatric treatment includes treatment of symptoms, cognition and functional deficits, suicide, and medical co-morbidities. VA Collaborative Study 572 examined a large sample of male and female veterans with schizophrenia (n=3,942) and with bipolar disorder (n=5,414) with phenotyping and genomic analyses. We present the results to date and future directions.

METHODS

All veterans received a structured diagnostic interview and assessments of suicidal ideation and behavior, PTSD, and health. Veterans with schizophrenia were assessed for negative symptoms and lifetime depression. All were assessed with a cognitive and functional capacity assessment. Data for genome wide association studies were collected. Controls came from the VA Million Veteran Program.

RESULTS

Suicidal ideation or behavior was present in 66%. Cognitive and functional deficits were consistent with previous studies. 40% of the veterans with schizophrenia had a lifetime major depressive episode and PTSD was present in over 30%. Polygenic risk score (PRS) analyses indicated that cognitive and functional deficits overlapped with PRS for cognition, education, and intelligence in the general population and PRS for suicidal ideation and behavior correlated with previous PRS for depression and suicidal ideation and behavior, as did the PRS for PTSD.

DISCUSSION

Results to date provide directions for personalization of treatment in SMI, veterans with SMI, and veterans in general. The results of the genomic analyses suggest that cognitive deficits in SMI may be associated with general population features. Upcoming genomic analyses will reexamine the issues above, as well as genomic factors associated with smoking, substance abuse, negative symptoms, and treatment response.

Genome-wide analyses of smoking behaviors in schizophrenia: Findings from the Psychiatric Genomics Consortium

While 17% of US adults use tobacco regularly, smoking rates among persons with schizophrenia are upwards of 60%. Research supports a shared etiological basis for smoking and schizophrenia, including findings from genome-wide association studies (GWAS). However, few studies have directly tested whether the same or distinct genetic variants also influence smoking behavior among schizophrenia cases. Using data from the Psychiatric Genomics Consortium (PGC) study of schizophrenia (35476 cases, 46839 controls), we estimated genetic correlations between these traits and tested whether polygenic risk scores (PRS) constructed from the results of smoking behaviors GWAS were associated with schizophrenia risk or smoking behaviors among schizophrenia cases. Results indicated significant genetic correlations of schizophrenia with smoking initiation (rg = 0.159; P = 5.05 × 10-10), cigarettes-smoked-per-day (rg = 0.094; P = 0.006), and age-of-onset of smoking (rg = 0.10; P = 0.009). Comparing smoking behaviors among schizophrenia cases to the general population, we observe positive genetic correlations for smoking initiation (rg = 0.624, P = 0.002) and cigarettes-smoked-per-day (rg = 0.689, P = 0.120). Similarly, TAG-based PRS for smoking initiation and cigarettes-smoked-per-day were significantly associated with smoking initiation (P = 3.49 × 10-5) and cigarettes-smoked-per-day (P = 0.007) among schizophrenia cases. We performed the first GWAS of smoking behavior among schizophrenia cases and identified a novel association with cigarettes-smoked-per-day upstream of the TMEM106B gene on chromosome 7p21.3 (rs148253479, P = 3.18 × 10-8, n = 3520). Results provide evidence of a partially shared genetic basis for schizophrenia and smoking behaviors. Additionally, genetic risk factors for smoking behaviors were largely shared across schizophrenia and non-schizophrenia populations. Future research should address mechanisms underlying these associations to aid both schizophrenia and smoking treatment and prevention efforts.

Genome-Wide Association Studies of Schizophrenia and Bipolar Disorder in a Diverse Cohort of US Veterans

BACKGROUND

Schizophrenia (SCZ) and bipolar disorder (BIP) are debilitating neuropsychiatric disorders, collectively affecting 2% of the world's population. Recognizing the major impact of these psychiatric disorders on the psychosocial function of more than 200 000 US Veterans, the Department of Veterans Affairs (VA) recently completed genotyping of more than 8000 veterans with SCZ and BIP in the Cooperative Studies Program (CSP) #572.

METHODS

We performed genome-wide association studies (GWAS) in CSP #572 and benchmarked the predictive value of polygenic risk scores (PRS) constructed from published findings. We combined our results with available summary statistics from several recent GWAS, realizing the largest and most diverse studies of these disorders to date.

RESULTS

Our primary GWAS uncovered new associations between CHD7 variants and SCZ, and novel BIP associations with variants in Sortilin Related VPS10 Domain Containing Receptor 3 (SORCS3) and downstream of PCDH11X. Combining our results with published summary statistics for SCZ yielded 39 novel susceptibility loci including CRHR1, and we identified 10 additional findings for BIP (28 326 cases and 90 570 controls). PRS trained on published GWAS were significantly associated with case-control status among European American (P < 10-30) and African American (P < .0005) participants in CSP #572.

CONCLUSIONS

We have demonstrated that published findings for SCZ and BIP are robustly generalizable to a diverse cohort of US veterans. Leveraging available summary statistics from GWAS of global populations, we report 52 new susceptibility loci and improved fine-mapping resolution for dozens of previously reported associations.

Contributions of common genetic variants to risk of schizophrenia among individuals of African and Latino ancestry

Schizophrenia is a common, chronic and debilitating neuropsychiatric syndrome affecting tens of millions of individuals worldwide. While rare genetic variants play a role in the etiology of schizophrenia, most of the currently explained liability is within common variation, suggesting that variation predating the human diaspora out of Africa harbors a large fraction of the common variant attributable heritability. However, common variant association studies in schizophrenia have concentrated mainly on cohorts of European descent. We describe genome-wide association studies of 6152 cases and 3918 controls of admixed African ancestry, and of 1234 cases and 3090 controls of Latino ancestry, representing the largest such study in these populations to date. Combining results from the samples with African ancestry with summary statistics from the Psychiatric Genomics Consortium (PGC) study of schizophrenia yielded seven newly genome-wide significant loci, and we identified an additional eight loci by incorporating the results from samples with Latino ancestry. Leveraging population differences in patterns of linkage disequilibrium, we achieve improved fine-mapping resolution at 22 previously reported and 4 newly significant loci. Polygenic risk score profiling revealed improved prediction based on trans-ancestry meta-analysis results for admixed African (Nagelkerke's R2 = 0.032; liability R2 = 0.017; P < 10-52), Latino (Nagelkerke's R2 = 0.089; liability R2 = 0.021; P < 10-58), and European individuals (Nagelkerke's R2 = 0.089; liability R2 = 0.037; P < 10-113), further highlighting the advantages of incorporating data from diverse human populations.

Support Groups and Individual Mental Health Care via Video Conferencing for Frontline Clinicians During the COVID-19 Pandemic

BACKGROUND

The current coronavirus disease 2019 (COVID-19) pandemic has put an enormous stress on the mental health of frontline health care workers.

OBJECTIVE

Psychiatry departments in medical centers need to develop support systems to help our colleagues cope with this stress.

METHODS

We developed recurring peer support groups via videoconferencing and telephone for physicians, resident physicians, and nursing staff, focusing on issues and emotions related to their frontline clinical work with COVID patients in our medical center which was designated as a COVID-only hospital by the state. These groups are led by attending psychiatrists and psychiatry residents. In addition, we also deployed a system of telehealth individual counseling by attending psychiatrists.

RESULTS

Anxiety was high in the beginning of our weekly groups, dealing with fear of contracting COVID or spreading COVID to family members and the stress of social distancing. Later, the focus was also on the impairment of the traditional clinician-patient relationship by the characteristics of this disease and the associated moral challenges and trauma. Clinicians were helped to cope with these issues through group processes such as ventilation of feelings, peer support, consensual validation, peer-learning, and interventions by group facilitators. People with severe anxiety or desiring confidentiality were helped through individual interventions.

CONCLUSIONS

Our experience suggests that this method of offering telehealth peer support groups and individual counseling is a useful model for other centers to adapt to emotionally support frontline clinical workers in this ongoing worldwide crisis.

Genome-wide association study of cognitive performance in U.S. veterans with schizophrenia or bipolar disorder

Cognitive impairment is a frequent and serious problem in patients with various forms of severe mental illnesses (SMI), including schizophrenia (SZ) and bipolar disorder (BP). Recent research suggests genetic links to several cognitive phenotypes in both SMI and in the general population. Our goal in this study was to identify potential genomic signatures of cognitive functioning in veterans with severe mental illness and compare them to previous findings for cognition across different populations. Veterans Affairs (VA) Cooperative Studies Program (CSP) Study #572 evaluated cognitive and functional capacity measures among SZ and BP patients. In conjunction with the VA Million Veteran Program, 3,959 European American (1,095 SZ, 2,864 BP) and 2,601 African American (1,095 SZ, 2,864 BP) patients were genotyped using a custom Affymetrix Axiom Biobank array. We performed a genome-wide association study of global cognitive functioning, constructed polygenic scores for SZ and cognition in the general population, and examined genetic correlations with 2,626 UK Biobank traits. Although no single locus attained genome-wide significance, observed allelic effects were strongly consistent with previous studies. We observed robust associations between global cognitive functioning and polygenic scores for cognitive performance, intelligence, and SZ risk. We also identified significant genetic correlations with several cognition-related traits in UK Biobank. In a diverse cohort of U.S. veterans with SZ or BP, we demonstrate broad overlap of common genetic effects on cognition in the general population, and find that greater polygenic loading for SZ risk is associated with poorer cognitive performance.

GWAS of Suicide Attempt in Psychiatric Disorders and Association With Major Depression Polygenic Risk Scores

OBJECTIVE

More than 90% of people who attempt suicide have a psychiatric diagnosis; however, twin and family studies suggest that the genetic etiology of suicide attempt is partially distinct from that of the psychiatric disorders themselves. The authors present the largest genome-wide association study (GWAS) on suicide attempt, using cohorts of individuals with major depressive disorder, bipolar disorder, and schizophrenia from the Psychiatric Genomics Consortium.

METHODS

The samples comprised 1,622 suicide attempters and 8,786 nonattempters with major depressive disorder; 3,264 attempters and 5,500 nonattempters with bipolar disorder; and 1,683 attempters and 2,946 nonattempters with schizophrenia. A GWAS on suicide attempt was performed by comparing attempters to nonattempters with each disorder, followed by a meta-analysis across disorders. Polygenic risk scoring was used to investigate the genetic relationship between suicide attempt and the psychiatric disorders.

RESULTS

Three genome-wide significant loci for suicide attempt were found: one associated with suicide attempt in major depressive disorder, one associated with suicide attempt in bipolar disorder, and one in the meta-analysis of suicide attempt in mood disorders. These associations were not replicated in independent mood disorder cohorts from the UK Biobank and iPSYCH. No significant associations were found in the meta-analysis of all three disorders. Polygenic risk scores for major depression were significantly associated with suicide attempt in major depressive disorder (R²=0.25%), bipolar disorder (R²=0.24%), and schizophrenia (R²=0.40%).

CONCLUSIONS

This study provides new information on genetic associations and demonstrates that genetic liability for major depression increases risk for suicide attempt across psychiatric disorders. Further collaborative efforts to increase sample size may help to robustly identify genetic associations and provide biological insights into the etiology of suicide attempt.

Ultra-Rare Genetic Variation in the Epilepsies: A Whole-Exome Sequencing Study of 17,606 Individuals

Sequencing-based studies have identified novel risk genes associated with severe epilepsies and revealed an excess of rare deleterious variation in less-severe forms of epilepsy. To identify the shared and distinct ultra-rare genetic risk factors for different types of epilepsies, we performed a whole-exome sequencing (WES) analysis of 9,170 epilepsy-affected individuals and 8,436 controls of European ancestry. We focused on three phenotypic groups: severe developmental and epileptic encephalopathies (DEEs), genetic generalized epilepsy (GGE), and non-acquired focal epilepsy (NAFE). We observed that compared to controls, individuals with any type of epilepsy carried an excess of ultra-rare, deleterious variants in constrained genes and in genes previously associated with epilepsy; we saw the strongest enrichment in individuals with DEEs and the least strong in individuals with NAFE. Moreover, we found that inhibitory GABA_A receptor genes were enriched for missense variants across all three classes of epilepsy, whereas no enrichment was seen in excitatory receptor genes. The larger gene groups for the GABAergic pathway or cation channels also showed a significant mutational burden in DEEs and GGE. Although no single gene surpassed exome-wide significance among individuals with GGE or NAFE, highly constrained genes and genes encoding ion channels were among the lead associations; such genes included CACNA1G, EEF1A2, and GABRG2 for GGE and LGI1, TRIM3, and GABRG2 for NAFE. Our study, the largest epilepsy WES study to date, confirms a convergence in the genetics of severe and less-severe epilepsies associated with ultra-rare coding variation, and it highlights a ubiquitous role for GABAergic inhibition in epilepsy etiology.

Contribution of Rare Copy Number Variants to Bipolar Disorder Risk Is Limited to Schizoaffective Cases

BACKGROUND

Genetic risk for bipolar disorder (BD) is conferred through many common alleles, while a role for rare copy number variants (CNVs) is less clear. Subtypes of BD including schizoaffective disorder bipolar type (SAB), bipolar I disorder (BD I), and bipolar II disorder (BD II) differ according to the prominence and timing of psychosis, mania, and depression. The genetic factors contributing to the combination of symptoms among these subtypes are poorly understood.

METHODS

Rare large CNVs were analyzed in 6353 BD cases (3833 BD I [2676 with psychosis, 850 without psychosis, and 307 with unknown psychosis history], 1436 BD II, 579 SAB, and 505 BD not otherwise specified) and 8656 controls. CNV burden and a polygenic risk score (PRS) for schizophrenia were used to evaluate the relative contributions of rare and common variants to risk of BD, BD subtypes, and psychosis.

RESULTS

CNV burden did not differ between BD and controls when treated as a single diagnostic entity. However, burden in SAB was increased relative to controls (p = .001), BD I (p = .0003), and BD II (p = .0007). Burden and schizophrenia PRSs were increased in SAB compared with BD I with psychosis (CNV p = .0007, PRS p = .004), and BD I without psychosis (CNV p = .0004, PRS p = 3.9 × 10^-5). Within BD I, psychosis was associated with increased schizophrenia PRSs (p = .005) but not CNV burden.

CONCLUSIONS

CNV burden in BD is limited to SAB. Rare and common genetic variants may contribute differently to risk for psychosis and perhaps other classes of psychiatric symptoms.

Estimation of Genetic Correlation via Linkage Disequilibrium Score Regression and Genomic Restricted Maximum Likelihood

Genetic correlation is a key population parameter that describes the shared genetic architecture of complex traits and diseases. It can be estimated by current state-of-art methods, i.e., linkage disequilibrium score regression (LDSC) and genomic restricted maximum likelihood (GREML). The massively reduced computing burden of LDSC compared to GREML makes it an attractive tool, although the accuracy (i.e., magnitude of standard errors) of LDSC estimates has not been thoroughly studied. In simulation, we show that the accuracy of GREML is generally higher than that of LDSC. When there is genetic heterogeneity between the actual sample and reference data from which LD scores are estimated, the accuracy of LDSC decreases further. In real data analyses estimating the genetic correlation between schizophrenia (SCZ) and body mass index, we show that GREML estimates based on ∼150,000 individuals give a higher accuracy than LDSC estimates based on ∼400,000 individuals (from combined meta-data). A GREML genomic partitioning analysis reveals that the genetic correlation between SCZ and height is significantly negative for regulatory regions, which whole genome or LDSC approach has less power to detect. We conclude that LDSC estimates should be carefully interpreted as there can be uncertainty about homogeneity among combined meta-datasets. We suggest that any interesting findings from massive LDSC analysis for a large number of complex traits should be followed up, where possible, with more detailed analyses with GREML methods, even if sample sizes are lesser.

Identification of shared genetic variants between schizophrenia and lung cancer

Epidemiology studies suggest associations between schizophrenia and cancer. However, the underlying genetic mechanisms are not well understood, and difficult to identify from epidemiological data. We investigated if there is a shared genetic architecture between schizophrenia and cancer, with the aim to identify specific overlapping genetic loci. First, we performed genome-wide enrichment analysis and second, we analyzed specific loci jointly associated with schizophrenia and cancer by the conjunction false discovery rate. We analyzed the largest genome-wide association studies of schizophrenia and lung, breast, prostate, ovary, and colon-rectum cancer including more than 220,000 subjects, and included genetic association with smoking behavior. Polygenic enrichment of associations with lung cancer was observed in schizophrenia, and weak enrichment for the remaining cancer sites. After excluding the major histocompatibility complex region, we identified three independent loci jointly associated with schizophrenia and lung cancer. The strongest association included nicotinic acetylcholine receptors and is an established pleiotropic locus shared between lung cancer and smoking. The two other loci were independent of genetic association with smoking. Functional analysis identified downstream pleiotropic effects on epigenetics and gene-expression in lung and brain tissue. These findings suggest that genetic factors may explain partly the observed epidemiological association of lung cancer and schizophrenia.

Evidence for genetic heterogeneity between clinical subtypes of bipolar disorder

We performed a genome-wide association study of 6447 bipolar disorder (BD) cases and 12 639 controls from the International Cohort Collection for Bipolar Disorder (ICCBD). Meta-analysis was performed with prior results from the Psychiatric Genomics Consortium Bipolar Disorder Working Group for a combined sample of 13 902 cases and 19 279 controls. We identified eight genome-wide significant, associated regions, including a novel associated region on chromosome 10 (rs10884920; P=3.28 × 10-8) that includes the brain-enriched cytoskeleton protein adducin 3 (ADD3), a non-coding RNA, and a neuropeptide-specific aminopeptidase P (XPNPEP1). Our large sample size allowed us to test the heritability and genetic correlation of BD subtypes and investigate their genetic overlap with schizophrenia and major depressive disorder. We found a significant difference in heritability of the two most common forms of BD (BD I SNP-h2=0.35; BD II SNP-h2=0.25; P=0.02). The genetic correlation between BD I and BD II was 0.78, whereas the genetic correlation was 0.97 when BD cohorts containing both types were compared. In addition, we demonstrated a significantly greater load of polygenic risk alleles for schizophrenia and BD in patients with BD I compared with patients with BD II, and a greater load of schizophrenia risk alleles in patients with the bipolar type of schizoaffective disorder compared with patients with either BD I or BD II. These results point to a partial difference in the genetic architecture of BD subtypes as currently defined.

Cross-Disorder Psychiatric Genomics

PURPOSE OF REVIEW

The following review provides some description of the movement in cross-disorder psychiatric genomics toward addressing both comorbidity and polygenicity.

RECENT FINDINGS

We attempt to show how dimensional approaches to the phenotype have led to further addressing the problem of comorbidity of psychiatric diagnoses. And we also attempt to show how a dimensional approach to the genome, with different statistical methods from traditional genome-wide association analyses, has begun to resolve the problem of massive polygenicity.

SUMMARY

Cross-disorder research, of any area in psychiatry, arguably has the most potential to inform clinical diagnosis, early detection and prevention strategies, and pharmacological treatment research. Future research might leverage what we now know to inform developmental studies of risk and resilience.

A simple yet accurate correction for winner's curse can predict signals discovered in much larger genome scans

Motivation: For genetic studies, statistically significant variants explain far less trait variance than ‘sub-threshold’ association signals. To dimension follow-up studies, researchers need to accurately estimate ‘true’ effect sizes at each SNP, e.g. the true mean of odds ratios (ORs)/regression coefficients (RRs) or Z-score noncentralities. Naïve estimates of effect sizes incur winner’s curse biases, which are reduced only by laborious winner’s curse adjustments (WCAs). Given that Z-scores estimates can be theoretically translated on other scales, we propose a simple method to compute WCA for Z-scores, i.e. their true means/noncentralities.

Results:WCA of Z-scores shrinks these towards zero while, on P-value scale, multiple testing adjustment (MTA) shrinks P-values toward one, which corresponds to the zero Z-score value. Thus, WCA on Z-scores scale is a proxy for MTA on P-value scale. Therefore, to estimate Z-score noncentralities for all SNPs in genome scans, we propose FDR Inverse Quantile Transformation (FIQT). It (i) performs the simpler MTA of P-values using FDR and (ii) obtains noncentralities by back-transforming MTA P-values on Z-score scale. When compared to competitors, realistic simulations suggest that FIQT is more (i) accurate and (ii) computationally efficient by orders of magnitude. Practical application of FIQT to Psychiatric Genetic Consortium schizophrenia cohort predicts a non-trivial fraction of sub-threshold signals which become significant in much larger supersamples.

Conclusions: FIQT is a simple, yet accurate, WCA method for Z-scores (and ORs/RRs, via simple transformations).

Availability and Implementation: A 10 lines R function implementation is available at https://github.com/bacanusa/FIQT.

Contact:sabacanu@vcu.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Meta-analysis of Positive and Negative Symptoms Reveals Schizophrenia Modifier Genes

BACKGROUND

Evidence suggests that genetic factors may influence both schizophrenia (Scz) and its clinical presentation. In recent years, genome-wide association studies (GWAS) have demonstrated considerable success in identifying risk loci. Detection of "modifier loci" has the potential to further elucidate underlying disease processes.

METHODS

We performed GWAS of empirically derived positive and negative symptom scales in Irish cases from multiply affected pedigrees and a larger, independent case-control sample, subsequently combining these into a large Irish meta-analysis. In addition to single-SNP associations, we considered gene-based and pathway analyses to better capture convergent genetic effects, and to facilitate biological interpretation of these findings. Replication and testing of aggregate genetic effects was conducted using an independent European-American sample.

RESULTS

Though no single marker met the genome-wide significance threshold, genes and ontologies/pathways were significantly associated with negative and positive symptoms; notably, NKAIN2 and NRG1, respectively. We observed limited overlap in ontologies/pathways associated with different symptom profiles, with immune-related categories over-represented for negative symptoms, and addiction-related categories for positive symptoms. Replication analyses suggested that genes associated with clinical presentation are generalizable to non-Irish samples.

CONCLUSIONS

These findings strongly support the hypothesis that modifier loci contribute to the etiology of distinct Scz symptom profiles. The finding that previously implicated "risk loci" actually influence particular symptom dimensions has the potential to better delineate the roles of these genes in Scz etiology. Furthermore, the over-representation of distinct gene ontologies/pathways across symptom profiles suggests that the clinical heterogeneity of Scz is due in part to complex and diverse genetic factors.

Genome-wide association study reveals greater polygenic loading for schizophrenia in cases with a family history of illness

Genome-wide association studies (GWAS) of schizophrenia have yielded more than 100 common susceptibility variants, and strongly support a substantial polygenic contribution of a large number of small allelic effects. It has been hypothesized that familial schizophrenia is largely a consequence of inherited rather than environmental factors. We investigated the extent to which familiality of schizophrenia is associated with enrichment for common risk variants detectable in a large GWAS. We analyzed single nucleotide polymorphism (SNP) data for cases reporting a family history of psychotic illness (N = 978), cases reporting no such family history (N = 4,503), and unscreened controls (N = 8,285) from the Psychiatric Genomics Consortium (PGC1) study of schizophrenia. We used a multinomial logistic regression approach with model-fitting to detect allelic effects specific to either family history subgroup. We also considered a polygenic model, in which we tested whether family history positive subjects carried more schizophrenia risk alleles than family history negative subjects, on average. Several individual SNPs attained suggestive but not genome-wide significant association with either family history subgroup. Comparison of genome-wide polygenic risk scores based on GWAS summary statistics indicated a significant enrichment for SNP effects among family history positive compared to family history negative cases (Nagelkerke's R(2 ) = 0.0021; P = 0.00331; P-value threshold <0.4). Estimates of variability in disease liability attributable to the aggregate effect of genome-wide SNPs were significantly greater for family history positive compared to family history negative cases (0.32 and 0.22, respectively; P = 0.031). We found suggestive evidence of allelic effects detectable in large GWAS of schizophrenia that might be specific to particular family history subgroups. However, consideration of a polygenic risk score indicated a significant enrichment among family history positive cases for common allelic effects. Familial illness might, therefore, represent a more heritable form of schizophrenia, as suggested by previous epidemiological studies.

JEPEGMIX: gene-level joint analysis of functional SNPs in cosmopolitan cohorts

MOTIVATION

To increase detection power, gene level analysis methods are used to aggregate weak signals. To greatly increase computational efficiency, most methods use as input summary statistics from genome-wide association studies (GWAS). Subsequently, gene statistics are constructed using linkage disequilibrium (LD) patterns from a relevant reference panel. However, all methods, including our own Joint Effect on Phenotype of eQTL/functional single nucleotide polymorphisms (SNPs) associated with a Gene (JEPEG), assume homogeneous panels, e.g. European. However, this renders these tools unsuitable for the analysis of large cosmopolitan cohorts.

RESULTS

We propose a JEPEG extension, JEPEGMIX, which similar to one of our software tools, Direct Imputation of summary STatistics of unmeasured SNPs from MIXed ethnicity cohorts, is capable of estimating accurate LD patterns for cosmopolitan cohorts. JEPEGMIX uses this accurate LD estimates to (i) impute the summary statistics at unmeasured functional variants and (ii) test for the joint effect of all measured and imputed functional variants which are associated with a gene. We illustrate the performance of our tool by analyzing the GWAS meta-analysis summary statistics from the multi-ethnic Psychiatric Genomics Consortium Schizophrenia stage 2 cohort. This practical application supports the immune system being one of the main drivers of the process leading to schizophrenia.

AVAILABILITY AND IMPLEMENTATION

Software, annotation database and examples are available at http://dleelab.github.io/jepegmix/.

CONTACT

donghyung.lee@vcuhealth.org

SUPPLEMENTARY INFORMATION

Supplementary material is available at Bioinformatics online.

DISTMIX: direct imputation of summary statistics for unmeasured SNPs from mixed ethnicity cohorts

Motivation: To increase the signal resolution for large-scale meta-analyses of genome-wide association studies, genotypes at unmeasured single nucleotide polymorphisms (SNPs) are commonly imputed using large multi-ethnic reference panels. However, the ever increasing size and ethnic diversity of both reference panels and cohorts makes genotype imputation computationally challenging for moderately sized computer clusters. Moreover, genotype imputation requires subject-level genetic data, which unlike summary statistics provided by virtually all studies, is not publicly available. While there are much less demanding methods which avoid the genotype imputation step by directly imputing SNP statistics, e.g. Directly Imputing summary STatistics (DIST) proposed by our group, their implicit assumptions make them applicable only to ethnically homogeneous cohorts.

Results: To decrease computational and access requirements for the analysis of cosmopolitan cohorts, we propose DISTMIX, which extends DIST capabilities to the analysis of mixed ethnicity cohorts. The method uses a relevant reference panel to directly impute unmeasured SNP statistics based only on statistics at measured SNPs and estimated/user-specified ethnic proportions. Simulations show that the proposed method adequately controls the Type I error rates. The 1000 Genomes panel imputation of summary statistics from the ethnically diverse Psychiatric Genetic Consortium Schizophrenia Phase 2 suggests that, when compared to genotype imputation methods, DISTMIX offers comparable imputation accuracy for only a fraction of computational resources.

Availability and implementation: DISTMIX software, its reference population data, and usage examples are publicly available at http://code.google.com/p/distmix.

Contact:dlee4@vcu.edu

Supplementary information:Supplementary Data are available at Bioinformatics online.

Genome-wide gene pathway analysis of psychotic illness symptom dimensions based on a new schizophrenia-specific model of the OPCRIT

Empirically derived phenotypic measurements have the potential to enhance gene-finding efforts in schizophrenia. Previous research based on factor analyses of symptoms has typically included schizoaffective cases. Deriving factor loadings from analysis of only narrowly defined schizophrenia cases could yield more sensitive factor scores for gene pathway and gene ontology analyses. Using an Irish family sample, this study 1) factor analyzed clinician-rated Operational Criteria Checklist items in cases with schizophrenia only, 2) scored the full sample based on these factor loadings, and 3) implemented genome-wide association, gene-based, and gene-pathway analysis of these SCZ-based symptom factors (final N=507). Three factors emerged from the analysis of the schizophrenia cases: a manic, a depressive, and a positive symptom factor. In gene-based analyses of these factors, multiple genes had q<0.01. Of particular interest are findings for PTPRG and WBP1L, both of which were previously implicated by the Psychiatric Genomics Consortium study of SCZ; results from this study suggest that variants in these genes might also act as modifiers of SCZ symptoms. Gene pathway analyses of the first factor indicated over-representation of glutamatergic transmission, GABA-A receptor, and cyclic GMP pathways. Results suggest that these pathways may have differential influence on affective symptom presentation in schizophrenia.

Uncovering the hidden risk architecture of the schizophrenias: confirmation in three independent genome-wide association studies

OBJECTIVE

The authors sought to demonstrate that schizophrenia is a heterogeneous group of heritable disorders caused by different genotypic networks that cause distinct clinical syndromes.

METHOD

In a large genome-wide association study of cases with schizophrenia and controls, the authors first identified sets of interacting single-nucleotide polymorphisms (SNPs) that cluster within particular individuals (SNP sets) regardless of clinical status. Second, they examined the risk of schizophrenia for each SNP set and tested replicability in two independent samples. Third, they identified genotypic networks composed of SNP sets sharing SNPs or subjects. Fourth, they identified sets of distinct clinical features that cluster in particular cases (phenotypic sets or clinical syndromes) without regard for their genetic background. Fifth, they tested whether SNP sets were associated with distinct phenotypic sets in a replicable manner across the three studies.

RESULTS

The authors identified 42 SNP sets associated with a 70% or greater risk of schizophrenia, and confirmed 34 (81%) or more with similar high risk of schizophrenia in two independent samples. Seventeen networks of SNP sets did not share any SNP or subject. These disjoint genotypic networks were associated with distinct gene products and clinical syndromes (i.e., the schizophrenias) varying in symptoms and severity. Associations between genotypic networks and clinical syndromes were complex, showing multifinality and equifinality. The interactive networks explained the risk of schizophrenia more than the average effects of all SNPs (24%).

CONCLUSIONS

Schizophrenia is a group of heritable disorders caused by a moderate number of separate genotypic networks associated with several distinct clinical syndromes.

JEPEG: a summary statistics based tool for gene-level joint testing of functional variants

MOTIVATION

Gene expression is influenced by variants commonly known as expression quantitative trait loci (eQTL). On the basis of this fact, researchers proposed to use eQTL/functional information univariately for prioritizing single nucleotide polymorphisms (SNPs) signals from genome-wide association studies (GWAS). However, most genes are influenced by multiple eQTLs which, thus, jointly affect any downstream phenotype. Therefore, when compared with the univariate prioritization approach, a joint modeling of eQTL action on phenotypes has the potential to substantially increase signal detection power. Nonetheless, a joint eQTL analysis is impeded by (i) not measuring all eQTLs in a gene and/or (ii) lack of access to individual genotypes.

RESULTS

We propose joint effect on phenotype of eQTL/functional SNPs associated with a gene (JEPEG), a novel software tool which uses only GWAS summary statistics to (i) impute the summary statistics at unmeasured eQTLs and (ii) test for the joint effect of all measured and imputed eQTLs in a gene. We illustrate the behavior/performance of the developed tool by analysing the GWAS meta-analysis summary statistics from the Psychiatric Genomics Consortium Stage 1 and the Genetic Consortium for Anorexia Nervosa.

CONCLUSIONS

Applied analyses results suggest that JEPEG complements commonly used univariate GWAS tools by: (i) increasing signal detection power via uncovering (a) novel genes or (b) known associated genes in smaller cohorts and (ii) assisting in fine-mapping of challenging regions, e.g. major histocompatibility complex for schizophrenia.

AVAILABILITY AND IMPLEMENTATION

JEPEG, its associated database of eQTL SNPs and usage examples are publicly available at http://code.google.com/p/jepeg/.

CONTACT

dlee4@vcu.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Polygenic dissection of diagnosis and clinical dimensions of bipolar disorder and schizophrenia

Bipolar disorder and schizophrenia are two often severe disorders with high heritabilities. Recent studies have demonstrated a large overlap of genetic risk loci between these disorders but diagnostic and molecular distinctions still remain. Here, we perform a combined genome-wide association study (GWAS) of 19 779 bipolar disorder (BP) and schizophrenia (SCZ) cases versus 19 423 controls, in addition to a direct comparison GWAS of 7129 SCZ cases versus 9252 BP cases. In our case-control analysis, we identify five previously identified regions reaching genome-wide significance (CACNA1C, IFI44L, MHC, TRANK1 and MAD1L1) and a novel locus near PIK3C2A. We create a polygenic risk score that is significantly different between BP and SCZ and show a significant correlation between a BP polygenic risk score and the clinical dimension of mania in SCZ patients. Our results indicate that first, combining diseases with similar genetic risk profiles improves power to detect shared risk loci and second, that future direct comparisons of BP and SCZ are likely to identify loci with significant differential effects. Identifying these loci should aid in the fundamental understanding of how these diseases differ biologically. These findings also indicate that combining clinical symptom dimensions and polygenic signatures could provide additional information that may someday be used clinically.

Molecular validation of the schizophrenia spectrum

BACKGROUND

Early descriptive work and controlled family and adoption studies support the hypothesis that a range of personality and nonschizophrenic psychotic disorders aggregate in families of schizophrenic probands. Can we validate, using molecular polygene scores from genome-wide association studies (GWAS), this schizophrenia spectrum?

METHODS

The predictive value of polygenic findings reported by the Psychiatric GWAS Consortium (PGC) was applied to 4 groups of relatives from the Irish Study of High-Density Schizophrenia Families (ISHDSF; N = 836) differing on their assignment within the schizophrenia spectrum. Genome-wide single nucleotide polymorphism data for affected and unaffected relatives were used to construct per-individual polygene risk scores based on the PGC stage-I results. We compared mean polygene scores in the ISHDSF with mean scores in ethnically matched population controls (N = 929).

RESULTS

The schizophrenia polygene score differed significantly across diagnostic categories and was highest in those with narrow schizophrenia spectrum, lowest in those with no psychiatric illness, and in-between in those classified in the intermediate, broad, and very broad schizophrenia spectrum. Relatives of all of these groups of affected subjects, including those with no diagnosis, had schizophrenia polygene scores significantly higher than the control sample.

CONCLUSIONS

In the relatives of high-density families, the observed pattern of enrichment of molecular indices of schizophrenia risk suggests an underlying, continuous liability distribution and validates, using aggregate common risk alleles, a genetic basis for the schizophrenia spectrum disorders. In addition, as predicted by genetic theory, nonpsychotic members of multiply-affected schizophrenia families are significantly enriched for replicated, polygenic risk variants compared with the general population.

Association study of 83 candidate genes for bipolar disorder in chromosome 6q selected using an evidence-based prioritization algorithm

BACKGROUND

Prior genome-scans of bipolar disorder have revealed chromosome 6q22 as a promising candidate region. However, linkage disequilibrium (LD) mapping studies have yet to identify replicated susceptibility loci.

METHODS

We analyzed 1,422 LD-tagging single nucleotide polymorphisms (SNPs) in 83 genes to test single-marker and locus-wide evidence of association with bipolar disorder in the NIMH Genetics Initiative bipolar pedigrees and the Portuguese Island Collection (PIC) (N = 1,093 in 528 informative pairs). Both studies previously demonstrated significant evidence of linkage to 6q. SNPs were genotyped using an Illumina iSelect genotyping array which employs the Infinium assay. Evidence of single-marker association was assessed using the generalized disequilibrium test (GDT). Empirical estimates of gene-wide significance were obtained by permutation (via 100,000 gene-dropping simulations) of Fisher's combined test of P-values for each locus.

RESULTS

No single variant yielded significant experiment-wide evidence of association, for either the combined sample or in each subsample. Our gene-dropping simulations identified nominally significant gene-wide associations with multiple loci, of which NT5DC1 in the NIMH subsample and CCNC in the PIC were the strongest candidates. However, no one gene consistently exceeded empirical significance criteria in both independent samples or survived Bonferroni correction for the number of genes tested.

CONCLUSIONS

Using a gene-based approach to family-based association, we identified gene-wide associations with several genes, though no single locus was significantly associated with bipolar disorder in both cohorts. This suggests that chromosome 6q may harbor multiple susceptibility loci or that complex patterns of LD in this region may confound approaches based on common SNPs. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.

DIST: direct imputation of summary statistics for unmeasured SNPs

Motivation: Genotype imputation methods are used to enhance the resolution of genome-wide association studies, and thus increase the detection rate for genetic signals. Although most studies report all univariate summary statistics, many of them limit the access to subject-level genotypes. Because such an access is required by all genotype imputation methods, it is helpful to develop methods that impute summary statistics without going through the interim step of imputing genotypes. Even when subject-level genotypes are available, due to the substantial computational cost of the typical genotype imputation, there is a need for faster imputation methods.

Results: Direct Imputation of summary STatistics (DIST) imputes the summary statistics of untyped variants without first imputing their subject-level genotypes. This is achieved by (i) using the conditional expectation formula for multivariate normal variates and (ii) using the correlation structure from a relevant reference population. When compared with genotype imputation methods, DIST (i) requires only a fraction of their computational resources, (ii) has comparable imputation accuracy for independent subjects and (iii) is readily applicable to the imputation of association statistics coming from large pedigree data. Thus, the proposed application is useful for a fast imputation of summary results for (i) studies of unrelated subjects, which (a) do not provide subject-level genotypes or (b) have a large size and (ii) family association studies.

Availability and implementation: Pre-compiled executables built under commonly used operating systems are publicly available at http://code.google.com/p/dist/.

Contact:dlee4@vcu.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Association study of 167 candidate genes for schizophrenia selected by a multi-domain evidence-based prioritization algorithm and neurodevelopmental hypothesis

Integrating evidence from multiple domains is useful in prioritizing disease candidate genes for subsequent testing. We ranked all known human genes (n = 3819) under linkage peaks in the Irish Study of High-Density Schizophrenia Families using three different evidence domains: 1) a meta-analysis of microarray gene expression results using the Stanley Brain collection, 2) a schizophrenia protein-protein interaction network, and 3) a systematic literature search. Each gene was assigned a domain-specific p-value and ranked after evaluating the evidence within each domain. For comparison to this ranking process, a large-scale candidate gene hypothesis was also tested by including genes with Gene Ontology terms related to neurodevelopment. Subsequently, genotypes of 3725 SNPs in 167 genes from a custom Illumina iSelect array were used to evaluate the top ranked vs. hypothesis selected genes. Seventy-three genes were both highly ranked and involved in neurodevelopment (category 1) while 42 and 52 genes were exclusive to neurodevelopment (category 2) or highly ranked (category 3), respectively. The most significant associations were observed in genes PRKG1, PRKCE, and CNTN4 but no individual SNPs were significant after correction for multiple testing. Comparison of the approaches showed an excess of significant tests using the hypothesis-driven neurodevelopment category. Random selection of similar sized genes from two independent genome-wide association studies (GWAS) of schizophrenia showed the excess was unlikely by chance. In a further meta-analysis of three GWAS datasets, four candidate SNPs reached nominal significance. Although gene ranking using integrated sources of prior information did not enrich for significant results in the current experiment, gene selection using an a priori hypothesis (neurodevelopment) was superior to random selection. As such, further development of gene ranking strategies using more carefully selected sources of information is warranted.

The genomic psychiatry cohort: partners in discovery

The Genomic Psychiatry Cohort (GPC) is a longitudinal resource designed to provide the necessary population-based sample for large-scale genomic studies, studies focusing on Research Domain Criteria (RDoC) and/or other alternate phenotype constructs, clinical and interventional studies, nested case-control studies, long-term disease course studies, and genomic variant-to-phenotype studies. We provide and will continue to encourage access to the GPC as an international resource. DNA and other biological samples and diagnostic data are available through the National Institute of Mental Health (NIMH) Repository. After appropriate review and approval by an advisory board, investigators are able to collaborate in, propose, and co-lead studies involving cohort participants.

Genome-wide association study of clinical dimensions of schizophrenia: polygenic effect on disorganized symptoms

OBJECTIVE

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.

METHOD

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.

RESULTS

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.

CONCLUSIONS

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

Experimental validation of candidate schizophrenia gene ZNF804A as target for hsa-miR-137

MicroRNAs (miRNAs) are small non-coding RNAs that mainly function as negative regulators of gene expression (Lai, 2002) and have been shown to be involved in schizophrenia etiology through genetic and expression studies (Burmistrova et al., 2007; Hansen et al., 2007a; Perkins et al., 2007; Beveridge et al., 2010; Kim et al., 2010). In a mega analysis of genome-wide association study (GWAS) of schizophrenia (SZ) and bipolar disorders (BP), a polymorphism (rs1625579) located in the primary transcript of a miRNA gene, hsa-miR-137, was reported to be strongly associated with SZ. Four SZ loci (CACNA1C, TCF4, CSMD1, C10orf26) achieving genome-wide significance in the same study were predicted and later experimentally validated (Kwon et al., 2011) as hsa-miR-137 targets. Here, using in silico, cellular and luciferase based approaches we also provide evidence that another well replicated candidate schizophrenia gene, ZNF804A, is also target for hsa-miR-137.

Convergent functional genomics of schizophrenia: from comprehensive understanding to genetic risk prediction

We have used a translational convergent functional genomics (CFG) approach to identify and prioritize genes involved in schizophrenia, by gene-level integration of genome-wide association study data with other genetic and gene expression studies in humans and animal models. Using this polyevidence scoring and pathway analyses, we identify top genes (DISC1, TCF4, MBP, MOBP, NCAM1, NRCAM, NDUFV2, RAB18, as well as ADCYAP1, BDNF, CNR1, COMT, DRD2, DTNBP1, GAD1, GRIA1, GRIN2B, HTR2A, NRG1, RELN, SNAP-25, TNIK), brain development, myelination, cell adhesion, glutamate receptor signaling, G-protein-coupled receptor signaling and cAMP-mediated signaling as key to pathophysiology and as targets for therapeutic intervention. Overall, the data are consistent with a model of disrupted connectivity in schizophrenia, resulting from the effects of neurodevelopmental environmental stress on a background of genetic vulnerability. In addition, we show how the top candidate genes identified by CFG can be used to generate a genetic risk prediction score (GRPS) to aid schizophrenia diagnostics, with predictive ability in independent cohorts. The GRPS also differentiates classic age of onset schizophrenia from early onset and late-onset disease. We also show, in three independent cohorts, two European American and one African American, increasing overlap, reproducibility and consistency of findings from single-nucleotide polymorphisms to genes, then genes prioritized by CFG, and ultimately at the level of biological pathways and mechanisms. Finally, we compared our top candidate genes for schizophrenia from this analysis with top candidate genes for bipolar disorder and anxiety disorders from previous CFG analyses conducted by us, as well as findings from the fields of autism and Alzheimer. Overall, our work maps the genomic and biological landscape for schizophrenia, providing leads towards a better understanding of illness, diagnostics and therapeutics. It also reveals the significant genetic overlap with other major psychiatric disorder domains, suggesting the need for improved nosology.

Network-assisted investigation of combined causal signals from genome-wide association studies in schizophrenia

With the recent success of genome-wide association studies (GWAS), a wealth of association data has been accomplished for more than 200 complex diseases/traits, proposing a strong demand for data integration and interpretation. A combinatory analysis of multiple GWAS datasets, or an integrative analysis of GWAS data and other high-throughput data, has been particularly promising. In this study, we proposed an integrative analysis framework of multiple GWAS datasets by overlaying association signals onto the protein-protein interaction network, and demonstrated it using schizophrenia datasets. Building on a dense module search algorithm, we first searched for significantly enriched subnetworks for schizophrenia in each single GWAS dataset and then implemented a discovery-evaluation strategy to identify module genes with consistent association signals. We validated the module genes in an independent dataset, and also examined them through meta-analysis of the related SNPs using multiple GWAS datasets. As a result, we identified 205 module genes with a joint effect significantly associated with schizophrenia; these module genes included a number of well-studied candidate genes such as DISC1, GNA12, GNA13, GNAI1, GPR17, and GRIN2B. Further functional analysis suggested these genes are involved in neuronal related processes. Additionally, meta-analysis found that 18 SNPs in 9 module genes had P_meta<1×10⁻⁴, including the gene HLA-DQA1 located in the MHC region on chromosome 6, which was reported in previous studies using the largest cohort of schizophrenia patients to date. These results demonstrated our bi-directional network-based strategy is efficient for identifying disease-associated genes with modest signals in GWAS datasets. This approach can be applied to any other complex diseases/traits where multiple GWAS datasets are available.

The recent success of genome-wide association studies (GWAS) has generated a wealth of genotyping data critical to studies of genetic architectures of many complex diseases. In contrast to traditional single marker analysis, an integrative analysis of multiple genes and the assessment of their joint effects have been particularly promising, especially upon the availability of many GWAS datasets and other high-throughput datasets for numerous complex diseases. In this study, we developed an integrative analysis framework for multiple GWAS datasets and demonstrated it in schizophrenia. We first constructed a GWAS-weighted protein-protein interaction (PPI) network and then applied a dense module search algorithm to identify subnetworks with combinatory disease effects. We applied combinatorial criteria for module selection based on permutation tests to determine whether the modules are significantly different from random gene sets and whether the modules are associated with the disease in investigation. Importantly, considering there are many complex diseases with multiple GWAS datasets available, we proposed a discovery-evaluation strategy to search for modules with consistent combined effects from two or more GWAS datasets. This approach can be applied to any diseases or traits that have two or more GWAS datasets available.

Genetic overlap of schizophrenia and bipolar disorder in a high-density linkage survey in the Portuguese Island population

Recent family and genome-wide association studies strongly suggest shared genetic risk factors for schizophrenia (SZ) and bipolar disorder (BP). However, linkage studies have not been used to test for statistically significant genome-wide overlap between them. Forty-seven Portuguese families with sibpairs concordant for SZ, BP, or psychosis (PSY, which includes either SZ or psychotic BP) were genotyped for over 57,000 markers using the Affymetrix 50K Xba SNP array. NPL and Kong and Cox LOD scores were calculated in Merlin for all three phenotypes. Empirical significance was determined using 1,000 gene-dropping simulations. Significance of genome-wide genetic overlap between SZ and BP was determined by the number of simulated BP scans having the same number of loci jointly linked with the real SZ scan, and vice versa. For all three phenotypes, a number of regions previously linked in this sample remained so. For BP, chromosome 1p36 achieved significance (11.54-15.71 MB, LOD = 3.51), whereas it was not even suggestively linked at lower marker densities, as did chromosome 11q14.1 (89.32-90.15 MB, NPL = 4.15). Four chromosomes had loci at which both SZ and BP had NPL ≥ 1.98, which was more than would be expected by chance (empirical P = 0.01 using simulated SZ scans; 0.07 using simulated BP scans), although they did not necessarily meet criteria for suggestive linkage individually. These results suggest that high-density marker maps may provide greater power and precision in linkage studies than lower density maps. They also further support the hypothesis that SZ and BP share at least some risk alleles.

Comprehensive gene-based association study of a chromosome 20 linked region implicates novel risk loci for depressive symptoms in psychotic illness

BACKGROUND

Prior genomewide scans of schizophrenia support evidence of linkage to regions of chromosome 20. However, association analyses have yet to provide support for any etiologically relevant variants.

METHODS

We analyzed 2988 LD-tagging single nucleotide polymorphisms (SNPs) in 327 genes on chromosome 20, to test for association with schizophrenia in 270 Irish high-density families (ISHDSF, N = 270 families, 1408 subjects). These SNPs were genotyped using an Illumina iSelect genotyping array which employs the Infinium assay. Given a previous report of novel linkage with chromosome 20p using latent classes of psychotic illness in this sample, association analysis was also conducted for each of five factor-derived scores based on the Operational Criteria Checklist for Psychotic Illness (delusions, hallucinations, mania, depression, and negative symptoms). Tests of association were conducted using the PDTPHASE and QPDTPHASE packages of UNPHASED. Empirical estimates of gene-wise significance were obtained by adaptive permutation of a) the smallest observed P-value and b) the threshold-truncated product of P-values for each locus.

RESULTS

While no single variant was significant after LD-corrected Bonferroni-correction, our gene-dropping analyses identified loci which exceeded empirical significance criteria for both gene-based tests. Namely, R3HDML and C20orf39 are significantly associated with depressive symptoms of schizophrenia (P(emp)<2×10⁻⁵) based on the minimum P-value and truncated-product methods, respectively.

CONCLUSIONS

Using a gene-based approach to family-based association, R3HDML and C20orf39 were found to be significantly associated with clinical dimensions of schizophrenia. These findings demonstrate the efficacy of gene-based analysis and support previous evidence that chromosome 20 may harbor schizophrenia susceptibility or modifier loci.

A bias-reducing pathway enrichment analysis of genome-wide association data confirmed association of the MHC region with schizophrenia

BACKGROUND

After the recent successes of genome-wide association studies (GWAS), one key challenge is to identify genetic variants that might have a significant joint effect on complex diseases but have failed to be identified individually due to weak to moderate marginal effect. One popular and effective approach is gene set based analysis, which investigates the joint effect of multiple functionally related genes (eg, pathways). However, a typical gene set analysis method is biased towards long genes, a problem that is especially severe in psychiatric diseases.

METHODS

A novel approach was proposed, namely generalised additive model (GAM) for GWAS (gamGWAS), for gene set enrichment analysis of GWAS data, specifically adjusting the gene length bias or the number of single-nucleotide polymorphisms per gene. GAM is applied to estimate the probability of a gene to be selected as significant given its gene length, followed by weighted resampling and computation of empirical p values for the rank of pathways. We demonstrated gamGWAS in two schizophrenia GWAS datasets from the International Schizophrenia Consortium and the Genetic Association Information Network.

RESULTS

The gamGWAS results not only confirmed previous findings, but also highlighted several immune related pathways. Comparison with other methods indicated that gamGWAS could effectively reduce the correlation between pathway p values and its median gene length.

CONCLUSION

gamGWAS can effectively relieve the long gene bias and generate reliable results for GWAS data analysis. It does not require genotype data or permutation of sample labels in the original GWAS data; thus, it is computationally efficient.

Two non-synonymous markers in PTPN21, identified by genome-wide association study data-mining and replication, are associated with schizophrenia

We conducted data-mining analyses of genome wide association (GWA) studies of the CATIE and MGS-GAIN datasets, and found 13 markers in the two physically linked genes, PTPN21 and EML5, showing nominally significant association with schizophrenia. Linkage disequilibrium (LD) analysis indicated that all 7 markers from PTPN21 shared high LD (r(2)>0.8), including rs2274736 and rs2401751, the two non-synonymous markers with the most significant association signals (rs2401751, P=1.10 × 10(-3) and rs2274736, P=1.21 × 10(-3)). In a meta-analysis of all 13 replication datasets with a total of 13,940 subjects, we found that the two non-synonymous markers are significantly associated with schizophrenia (rs2274736, OR=0.92, 95% CI: 0.86-0.97, P=5.45 × 10(-3) and rs2401751, OR=0.92, 95% CI: 0.86-0.97, P=5.29 × 10(-3)). One SNP (rs7147796) in EML5 is also significantly associated with the disease (OR=1.08, 95% CI: 1.02-1.14, P=6.43 × 10(-3)). These 3 markers remain significant after Bonferroni correction. Furthermore, haplotype conditioned analyses indicated that the association signals observed between rs2274736/rs2401751 and rs7147796 are statistically independent. Given the results that 2 non-synonymous markers in PTPN21 are associated with schizophrenia, further investigation of this locus is warranted.