Genome-wide common and rare variant analysis provides novel insights into clozapine-associated neutropenia.

This corrects the article DOI: 10.1038/mp.2016.97.

Schizophrenia copy number variants and associative learning.

Large-scale genomic studies have made major progress in identifying genetic risk variants for schizophrenia. A key finding from these studies is that there is an increased burden of genomic copy number variants (CNVs) in schizophrenia cases compared with controls. The mechanism through which these CNVs confer risk for the symptoms of schizophrenia, however, remains unclear. One possibility is that schizophrenia risk CNVs impact basic associative learning processes, abnormalities of which have long been associated with the disorder. To investigate whether genes in schizophrenia CNVs impact on specific phases of associative learning we combined human genetics with experimental gene expression studies in animals. In a sample of 11 917 schizophrenia cases and 16 416 controls, we investigated whether CNVs from patients with schizophrenia are enriched for genes expressed during the consolidation, retrieval or extinction of associative memories. We show that CNVs from cases are enriched for genes expressed during fear extinction in the hippocampus, but not genes expressed following consolidation or retrieval. These results suggest that CNVs act to impair inhibitory learning in schizophrenia, potentially contributing to the development of core symptoms of the disorder.

Genome-wide common and rare variant analysis provides novel insights into clozapine-associated neutropenia.

The antipsychotic clozapine is uniquely effective in the management of schizophrenia; however, its use is limited by its potential to induce agranulocytosis. The causes of this, and of its precursor neutropenia, are largely unknown, although genetic factors have an important role. We sought risk alleles for clozapine-associated neutropenia in a sample of 66 cases and 5583 clozapine-treated controls, through a genome-wide association study (GWAS), imputed human leukocyte antigen (HLA) alleles, exome array and copy-number variation (CNV) analyses. We then combined associated variants in a meta-analysis with data from the Clozapine-Induced Agranulocytosis Consortium (up to 163 cases and 7970 controls). In the largest combined sample to date, we identified a novel association with rs149104283 (odds ratio (OR)=4.32, P=1.79 × 10-8), intronic to transcripts of SLCO1B3 and SLCO1B7, members of a family of hepatic transporter genes previously implicated in adverse drug reactions including simvastatin-induced myopathy and docetaxel-induced neutropenia. Exome array analysis identified gene-wide associations of uncommon non-synonymous variants within UBAP2 and STARD9. We additionally provide independent replication of a previously identified variant in HLA-DQB1 (OR=15.6, P=0.015, positive predictive value=35.1%). These results implicate biological pathways through which clozapine may act to cause this serious adverse effect.

Genome-wide common and rare variant analysis provides novel insights into clozapine-associated neutropenia.

This corrects the article DOI: 10.1038/mp.2016.97.

Common alleles contribute to schizophrenia in CNV carriers.

The genetic architecture of schizophrenia is complex, involving risk alleles ranging from common alleles of weak effect to rare alleles of large effect, the best exemplar of the latter being large copy number variants (CNVs). It is currently unknown whether pathophysiology in those with defined rare mutations overlaps with that in other individuals with the disorder who do not share the same rare mutation. Under an extreme heterogeneity model, carriers of specific high-penetrance mutations form distinct subgroups. In contrast, under a polygenic threshold model, high-penetrance rare allele carriers possess many risk factors, of which the rare allele is the only one, albeit an important, factor. Under the latter model, cases with rare mutations can be expected to share some common risk alleles, and therefore pathophysiological mechanisms, with cases without the same mutation. Here we show that, compared with controls, individuals with schizophrenia who have known pathogenic CNVs carry an excess burden of common risk alleles (P=2.25 × 10(-17)) defined from a genome-wide association study largely based on individuals without known CNVs. Our finding is not consistent with an extreme heterogeneity model for CNV carriers, but does offer support for the polygenic threshold model of schizophrenia. That this is so provides support for the notion that studies aiming to model the effects of rare variation may uncover pathophysiological mechanisms of relevance to those with the disorder more widely.

Copy number variation in bipolar disorder.

Large (>100 kb), rare (<1% in the population) copy number variants (CNVs) have been shown to confer risk for schizophrenia (SZ), but the findings for bipolar disorder (BD) are less clear. In a new BD sample from the United Kingdom (n=2591), we have examined the occurrence of CNVs and compared this with previously reported samples of 6882 SZ and 8842 control subjects. When combined with previous data, we find evidence for a contribution to BD for three SZ-associated CNV loci: duplications at 1q21.1 (P=0.022), deletions at 3q29 (P=0.03) and duplications at 16p11.2 (P=2.3 × 10(-4)). The latter survives multiple-testing correction for the number of recurrent large CNV loci in the genome. Genes in 20 regions (total of 55 genes) were enriched for rare exonic CNVs among BD cases, but none of these survives correction for multiple testing. Finally, our data provide strong support for the hypothesis of a lesser contribution of very large (>500 kb) CNVs in BD compared with SZ, most notably for deletions >1 Mb (P=9 × 10(-4)).

Copy number variation in schizophrenia in Sweden.

Schizophrenia (SCZ) is a highly heritable neuropsychiatric disorder of complex genetic etiology. Previous genome-wide surveys have revealed a greater burden of large, rare copy number variations (CNVs) in SCZ cases and identified multiple rare recurrent CNVs that increase risk of SCZ although with incomplete penetrance and pleiotropic effects. Identification of additional recurrent CNVs and biological pathways enriched for SCZ CNVs requires greater sample sizes. We conducted a genome-wide survey for CNVs associated with SCZ using a Swedish national sample (4719 cases and 5917 controls). High-confidence CNV calls were generated using genotyping array intensity data, and their effect on risk of SCZ was measured. Our data confirm increased burden of large, rare CNVs in SCZ cases as well as significant associations for recurrent 16p11.2 duplications, 22q11.2 deletions and 3q29 deletions. We report a novel association for 17q12 duplications (odds ratio=4.16, P=0.018), previously associated with autism and mental retardation but not SCZ. Intriguingly, gene set association analyses implicate biological pathways previously associated with SCZ through common variation and exome sequencing (calcium channel signaling and binding partners of the fragile X mental retardation protein). We found significantly increased burden of the largest CNVs (>500 kb) in genes present in the postsynaptic density, in genomic regions implicated via SCZ genome-wide association studies and in gene products localized to mitochondria and cytoplasm. Our findings suggest that multiple lines of genomic inquiry--genome-wide screens for CNVs, common variation and exonic variation--are converging on similar sets of pathways and/or genes.

Evidence that duplications of 22q11.2 protect against schizophrenia.

A number of large, rare copy number variants (CNVs) are deleterious for neurodevelopmental disorders, but large, rare, protective CNVs have not been reported for such phenotypes. Here we show in a CNV analysis of 47 005 individuals, the largest CNV analysis of schizophrenia to date, that large duplications (1.5-3.0 Mb) at 22q11.2--the reciprocal of the well-known, risk-inducing deletion of this locus--are substantially less common in schizophrenia cases than in the general population (0.014% vs 0.085%, OR=0.17, P=0.00086). 22q11.2 duplications represent the first putative protective mutation for schizophrenia.

Replication of bipolar disorder susceptibility alleles and identification of two novel genome-wide significant associations in a new bipolar disorder case-control sample.

We have conducted a genotyping study using a custom Illumina Infinium HD genotyping array, the ImmunoChip, in a new UK sample of 1218 bipolar disorder (BD) cases and 2913 controls that have not been used in any studies previously reported independently or in meta-analyses. The ImmunoChip was designed before the publication of the Psychiatric Genome-Wide Association Study Consortium Bipolar Disorder Working Group (PGC-BD) meta-analysis data. As such 3106 single-nucleotide polymorphisms (SNPs) with a P-value <1 × 10(-3) from the BD meta-analysis by Ferreira et al. were genotyped. We report support for two of the three most strongly associated chromosomal regions in the Ferreira study, CACNA1C (rs1006737, P=4.09 × 10(-4)) and 15q14 (rs2172835, P=0.043) but not ANK3 (rs10994336, P=0.912). We have combined our ImmunoChip data (569 quasi-independent SNPs from the 3016 SNPs genotyped) with the recently published PGC-BD meta-analysis data, using either the PGC-BD combined discovery and replication data where available or just the discovery data where the SNP was not typed in a replication sample in PGC-BD. Our data provide support for two regions, at ODZ4 and CACNA1C, with prior evidence for genome-wide significant (GWS) association in PGC-BD meta-analysis. In addition, the combined analysis shows two novel GWS associations. First, rs7296288 (P=8.97 × 10(-9), odds ratio (OR)=0.9), an intergenic polymorphism on chromosome 12 located between RHEBL1 and DHH. Second, rs3818253 (P=3.88 × 10(-8), OR=1.16), an intronic SNP on chromosome 20q11.2 in the gene TRPC4AP, which lies in a high linkage disequilibrium region along with the genes GSS and MYH7B.

Schizophrenia susceptibility alleles are enriched for alleles that affect gene expression in adult human brain.

It is widely thought that alleles that influence susceptibility to common diseases, including schizophrenia, will frequently do so through effects on gene expression. As only a small proportion of the genetic variance for schizophrenia has been attributed to specific loci, this remains an unproven hypothesis. The International Schizophrenia Consortium (ISC) recently reported a substantial polygenic contribution to that disorder, and that schizophrenia risk alleles are enriched among single-nucleotide polymorphisms (SNPs) selected for marginal evidence for association (P<0.5) from genome-wide association studies (GWAS). It follows that if schizophrenia susceptibility alleles are enriched for those that affect gene expression, those marginally associated SNPs, which are also expression quantitative trait loci (eQTLs), should carry more true association signals compared with SNPs that are not marginally associated. To test this, we identified marginally associated (P<0.5) SNPs from two of the largest available schizophrenia GWAS data sets. We assigned eQTL status to those SNPs based upon an eQTL data set derived from adult human brain. Using the polygenic score method of analysis reported by the ISC, we observed and replicated the observation that higher probability cis-eQTLs predicted schizophrenia better than those with a lower probability for being a cis-eQTL. Our data support the hypothesis that alleles conferring risk of schizophrenia are enriched among those that affect gene expression. Moreover, our data show that notwithstanding the likely developmental origin of schizophrenia, studies of adult brain tissue can, in principle, allow relevant susceptibility eQTLs to be identified.

Genome-wide linkage analysis of 972 bipolar pedigrees using single-nucleotide polymorphisms.

Because of the high costs associated with ascertainment of families, most linkage studies of Bipolar I disorder (BPI) have used relatively small samples. Moreover, the genetic information content reported in most studies has been less than 0.6. Although microsatellite markers spaced every 10 cM typically extract most of the genetic information content for larger multiplex families, they can be less informative for smaller pedigrees especially for affected sib pair kindreds. For these reasons we collaborated to pool family resources and carried out higher density genotyping. Approximately 1100 pedigrees of European ancestry were initially selected for study and were genotyped by the Center for Inherited Disease Research using the Illumina Linkage Panel 12 set of 6090 single-nucleotide polymorphisms. Of the ~1100 families, 972 were informative for further analyses, and mean information content was 0.86 after pruning for linkage disequilibrium. The 972 kindreds include 2284 cases of BPI disorder, 498 individuals with bipolar II disorder (BPII) and 702 subjects with recurrent major depression. Three affection status models (ASMs) were considered: ASM1 (BPI and schizoaffective disorder, BP cases (SABP) only), ASM2 (ASM1 cases plus BPII) and ASM3 (ASM2 cases plus recurrent major depression). Both parametric and non-parametric linkage methods were carried out. The strongest findings occurred at 6q21 (non-parametric pairs LOD 3.4 for rs1046943 at 119 cM) and 9q21 (non-parametric pairs logarithm of odds (LOD) 3.4 for rs722642 at 78 cM) using only BPI and schizoaffective (SA), BP cases. Both results met genome-wide significant criteria, although neither was significant after correction for multiple analyses. We also inspected parametric scores for the larger multiplex families to identify possible rare susceptibility loci. In this analysis, we observed 59 parametric LODs of 2 or greater, many of which are likely to be close to maximum possible scores. Although some linkage findings may be false positives, the results could help prioritize the search for rare variants using whole exome or genome sequencing.

De novo CNV analysis implicates specific abnormalities of postsynaptic signalling complexes in the pathogenesis of schizophrenia.

A small number of rare, recurrent genomic copy number variants (CNVs) are known to substantially increase susceptibility to schizophrenia. As a consequence of the low fecundity in people with schizophrenia and other neurodevelopmental phenotypes to which these CNVs contribute, CNVs with large effects on risk are likely to be rapidly removed from the population by natural selection. Accordingly, such CNVs must frequently occur as recurrent de novo mutations. In a sample of 662 schizophrenia proband-parent trios, we found that rare de novo CNV mutations were significantly more frequent in cases (5.1% all cases, 5.5% family history negative) compared with 2.2% among 2623 controls, confirming the involvement of de novo CNVs in the pathogenesis of schizophrenia. Eight de novo CNVs occurred at four known schizophrenia loci (3q29, 15q11.2, 15q13.3 and 16p11.2). De novo CNVs of known pathogenic significance in other genomic disorders were also observed, including deletion at the TAR (thrombocytopenia absent radius) region on 1q21.1 and duplication at the WBS (Williams-Beuren syndrome) region at 7q11.23. Multiple de novos spanned genes encoding members of the DLG (discs large) family of membrane-associated guanylate kinases (MAGUKs) that are components of the postsynaptic density (PSD). Two de novos also affected EHMT1, a histone methyl transferase known to directly regulate DLG family members. Using a systems biology approach and merging novel CNV and proteomics data sets, systematic analysis of synaptic protein complexes showed that, compared with control CNVs, case de novos were significantly enriched for the PSD proteome (P=1.72 × 10⁻6. This was largely explained by enrichment for members of the N-methyl-D-aspartate receptor (NMDAR) (P=4.24 × 10⁻6) and neuronal activity-regulated cytoskeleton-associated protein (ARC) (P=3.78 × 10⁻8) postsynaptic signalling complexes. In an analysis of 18 492 subjects (7907 cases and 10 585 controls), case CNVs were enriched for members of the NMDAR complex (P=0.0015) but not ARC (P=0.14). Our data indicate that defects in NMDAR postsynaptic signalling and, possibly, ARC complexes, which are known to be important in synaptic plasticity and cognition, play a significant role in the pathogenesis of schizophrenia.

Association analysis of NOTCH4 loci in schizophrenia using family and population-based controls.

A genetic association between NOTCH4 and schizophrenia has previously been proposed. Unsing all markers previously shown to be associated, we found no evidence for such in three independent family-based samples (n=519 parent-offspring trios), and a case-control sample derived from the same ethnic background as the original observation. These data strongly suggest that NOTCH4 is not a significant susceptibility allele for schizophrenia.

Fine mapping of ZNF804A and genome-wide significant evidence for its involvement in schizophrenia and bipolar disorder.

A recent genome-wide association study (GWAS) reported evidence for association between rs1344706 within ZNF804A (encoding zinc-finger protein 804A) and schizophrenia (P=1.61 × 10(-7)), and stronger evidence when the phenotype was broadened to include bipolar disorder (P=9.96 × 10(-9)). In this study we provide additional evidence for association through meta-analysis of a larger data set (schizophrenia/schizoaffective disorder N=18 945, schizophrenia plus bipolar disorder N=21 274 and controls N=38 675). We also sought to better localize the association signal using a combination of de novo polymorphism discovery in exons, pooled de novo polymorphism discovery spanning the genomic sequence of the locus and high-density linkage disequilibrium (LD) mapping. The meta-analysis provided evidence for association between rs1344706 that surpasses widely accepted benchmarks of significance by several orders of magnitude for both schizophrenia (P=2.5 × 10(-11), odds ratio (OR) 1.10, 95% confidence interval 1.07-1.14) and schizophrenia and bipolar disorder combined (P=4.1 × 10(-13), OR 1.11, 95% confidence interval 1.07-1.14). After de novo polymorphism discovery and detailed association analysis, rs1344706 remained the most strongly associated marker in the gene. The allelic association at the ZNF804A locus is now one of the most compelling in schizophrenia to date, and supports the accumulating data suggesting overlapping genetic risk between schizophrenia and bipolar disorder.

GWA study data mining and independent replication identify cardiomyopathy-associated 5 (CMYA5) as a risk gene for schizophrenia.

We conducted data-mining analyses using the Clinical Antipsychotic Trials of Intervention Effectiveness (CATIE) and molecular genetics of schizophrenia genome-wide association study supported by the genetic association information network (MGS-GAIN) schizophrenia data sets and performed bioinformatic prioritization for all the markers with P-values ≤0.05 in both data sets. In this process, we found that in the CMYA5 gene, there were two non-synonymous markers, rs3828611 and rs10043986, showing nominal significance in both the CATIE and MGS-GAIN samples. In a combined analysis of both the CATIE and MGS-GAIN samples, rs4704591 was identified as the most significant marker in the gene. Linkage disequilibrium analyses indicated that these markers were in low LD (3 828 611-rs10043986, r(2)=0.008; rs10043986-rs4704591, r(2)=0.204). In addition, CMYA5 was reported to be physically interacting with the DTNBP1 gene, a promising candidate for schizophrenia, suggesting that CMYA5 may be involved in the same biological pathway and process. On the basis of this information, we performed replication studies for these three single-nucleotide polymorphisms. The rs3828611 was found to have conflicting results in our Irish samples and was dropped out without further investigation. The other two markers were verified in 23 other independent data sets. In a meta-analysis of all 23 replication samples (family samples, 912 families with 4160 subjects; case-control samples, 11 380 cases and 15 021 controls), we found that both markers are significantly associated with schizophrenia (rs10043986, odds ratio (OR)=1.11, 95% confidence interval (CI)=1.04-1.18, P=8.2 × 10(-4) and rs4704591, OR=1.07, 95% CI=1.03-1.11, P=3.0 × 10(-4)). The results were also significant for the 22 Caucasian replication samples (rs10043986, OR=1.11, 95% CI=1.03-1.17, P=0.0026 and rs4704591, OR=1.07, 95% CI=1.02-1.11, P=0.0015). Furthermore, haplotype conditioned analyses indicated that the association signals observed at these two markers are independent. On the basis of these results, we concluded that CMYA5 is associated with schizophrenia and further investigation of the gene is warranted.

Polygenic dissection of the bipolar phenotype.

BACKGROUND: Recent data provide strong support for a substantial common polygenic contribution (i.e. many alleles each of small effect) to genetic susceptibility for schizophrenia and overlapping susceptibility for bipolar disorder. AIMS: To test hypotheses about the relationship between schizophrenia and psychotic types of bipolar disorder. METHOD: Using a polygenic score analysis to test whether schizophrenia polygenic risk alleles, en masse, significantly discriminate between individuals with bipolar disorder with and without psychotic features. The primary sample included 1829 participants with bipolar disorder and the replication sample comprised 506 people with bipolar disorder. RESULTS: The subset of participants with Research Diagnostic Criteria schizoaffective bipolar disorder (n = 277) were significantly discriminated from the remaining participants with bipolar disorder (n = 1552) in both the primary (P = 0.00059) and the replication data-sets (P = 0.0070). In contrast, those with psychotic bipolar disorder as a whole were not significantly different from those with non-psychotic bipolar disorder in either data-set. CONCLUSIONS: Genetic susceptibility influences at least two major domains of psychopathological variation in the schizophrenia-bipolar disorder clinical spectrum: one that relates to expression of a 'bipolar disorder-like' phenotype and one that is associated with expression of 'schizophrenia-like' psychotic symptoms.

The bipolar disorder risk allele at CACNA1C also confers risk of recurrent major depression and of schizophrenia.

Molecular genetic analysis offers opportunities to advance our understanding of the nosological relationship between psychiatric diagnostic categories in general, and the mood and psychotic disorders in particular. Strong evidence (P=7.0 × 10(-7)) of association at the polymorphism rs1006737 (within CACNA1C, the gene encoding the α-1C subunit of the L-type voltage-gated calcium channel) with the risk of bipolar disorder (BD) has recently been reported in a meta-analysis of three genome-wide association studies of BD, including our BD sample (N=1868) studied within the Wellcome Trust Case Control Consortium. Here, we have used our UK case samples of recurrent major depression (N=1196) and schizophrenia (N=479) and UK non-psychiatric comparison groups (N=15316) to examine the spectrum of phenotypic effect of the bipolar risk allele at rs1006737. We found that the risk allele conferred increased risk for schizophrenia (P=0.034) and recurrent major depression (P=0.013) with similar effect sizes to those previously observed in BD (allelic odds ratio ∼1.15). Our findings are evidence of some degree of overlap in the biological underpinnings of susceptibility to mental illness across the clinical spectrum of mood and psychotic disorders, and show that at least some loci can have a relatively general effect on susceptibility to diagnostic categories, as currently defined. Our findings will contribute to a better understanding of the pathogenesis of major psychiatric illness, and such knowledge should be useful in providing an etiological rationale for shaping psychiatric nosology, which is currently reliant entirely on descriptive clinical data.

Evidence for rare and common genetic risk variants for schizophrenia at protein kinase C, alpha.

We earlier reported a genome-wide significant linkage to schizophrenia at chromosome 17 that was identified in a single pedigree (C702) consisting of six affected, male siblings with DSM-IV schizophrenia and prominent mood symptoms. In this study, we adopted several approaches in an attempt to map the putative disease locus. First, mapping the source of linkage to chromosome 17 in pedigree C702. We refined the linkage region in family C702 to a 21-marker segment spanning 11.7 Mb at 17q23-q24 by genotyping a total of 50 microsatellites across chromosome 17 in the pedigree. Analysis of data from 1028 single nucleotide polymorphisms (SNPs) across the refined linkage region identified a single region of homozygosity present in pedigree C702 but not in 2938 UK controls. This spanned ~432 kb of the gene encoding protein kinase C, alpha (PRKCA), the encoded protein of which has been implicated in the pathogenesis of psychiatric disorders. Analysis of pedigree C702 by oligonucleotide-array comparative genome hybridization excluded the possibility that this region of homozygosity was because of a deletion. Mutation screening of PRKCA identified a rare, four-marker haplotype (C-HAP) in the 3' untranslated region of the gene, which was present in the homozygous state in all six affected members of pedigree C702. No other homozygotes were observed in genotype data for a total of 6597 unrelated Europeans (case N=1755, control N=3580 and parents of probands N=1262). Second, association analysis of C702 alleles at PRKCA. The low-frequency haplotype (C-HAP) showed a trend for association in a study of unrelated schizophrenia cases and controls from the UK (661 cases, 2824 controls, P=0.078 and odd ratio (OR)=1.9) and significant evidence for association when the sample was expanded to include cases with bipolar (N=710) and schizoaffective disorder (N=50) (psychosis sample: 1421 cases, 2824 controls, P=0.037 and OR=1.9). Given that all the affected members of C702 are male, we also undertook sex-specific analyses. This revealed that the association was strongest in males for both schizophrenia (446 male cases, 1421 male controls, P=0.008 and OR=3.9) and in the broader psychosis group (730 male cases, 1421 male controls, P=0.008 and OR=3.6). Analysis of C-HAP in follow-up samples from Ireland and Bulgaria revealed no evidence for association in either the whole sample or in males alone, and meta-analysis of all male psychosis samples yielded no significant evidence of association (969 male cases, 1939 male controls, 311 male probands P=0.304 and OR=1.4). Third, association mapping of the pedigree C702 linkage region. Independent of pedigree C702, genotype data from the Affymetrix 500k GeneChip set were available for 476 patients with schizophrenia and 2938 controls from the United Kingdom. SNPs in PRKCA showed evidence for association with schizophrenia that achieved gene-wide significance (P=0.027). Moreover, the same SNP was the most significantly associated marker out of the 1028 SNPs genotyped across the linkage region (rs873417, allelic P=0.0004). Follow-up genotyping in samples from Ireland, Bulgaria and Germany did not show consistent replication, but meta-analysis of all samples (4116 cases and 6491 controls) remained nominally significant (meta-analysis P=0.026, OR=1.1). We conclude that, although we have obtained convergent lines of evidence implicating both rare and common schizophrenia risk variants at PRKCA, none of these is individually compelling. However, the evidence across all approaches suggests that further study of this locus is warranted.

Strong genetic evidence for a selective influence of GABAA receptors on a component of the bipolar disorder phenotype.

Despite compelling evidence for a major genetic contribution to risk of bipolar mood disorder, conclusive evidence implicating specific genes or pathophysiological systems has proved elusive. In part this is likely to be related to the unknown validity of current phenotype definitions and consequent aetiological heterogeneity of samples. In the recent Wellcome Trust Case Control Consortium genome-wide association analysis of bipolar disorder (1868 cases, 2938 controls) one of the most strongly associated polymorphisms lay within the gene encoding the GABA(A) receptor beta1 subunit, GABRB1. Aiming to increase biological homogeneity, we sought the diagnostic subset that showed the strongest signal at this polymorphism and used this to test for independent evidence of association with other members of the GABA(A) receptor gene family. The index signal was significantly enriched in the 279 cases meeting Research Diagnostic Criteria for schizoaffective disorder, bipolar type (P=3.8 x 10(-6)). Independently, these cases showed strong evidence that variation in GABA(A) receptor genes influences risk for this phenotype (independent system-wide P=6.6 x 10(-5)) with association signals also at GABRA4, GABRB3, GABRA5 and GABRR3. [corrected] Our findings have the potential to inform understanding of presentation, pathogenesis and nosology of bipolar disorders. Our method of phenotype refinement may be useful in studies of other complex psychiatric and non-psychiatric disorders.