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

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

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.

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.

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

Analysis of exome sequence in 604 trios for recessive genotypes in schizophrenia.

Genetic associations involving both rare and common alleles have been reported for schizophrenia but there have been no systematic scans for rare recessive genotypes using fully phased trio data. Here, we use exome sequencing in 604 schizophrenia proband-parent trios to investigate the role of recessive (homozygous or compound heterozygous) nonsynonymous genotypes in the disorder. The burden of recessive genotypes was not significantly increased in probands at either a genome-wide level or in any individual gene after adjustment for multiple testing. At a system level, probands had an excess of nonsynonymous compound heterozygous genotypes (minor allele frequency, MAF ⩽ 1%) in voltage-gated sodium channels (VGSCs; eight in probands and none in parents, P = 1.5 × 10(-)(4)). Previous findings of multiple de novo loss-of-function mutations in this gene family, particularly SCN2A, in autism and intellectual disability provide biological and genetic plausibility for this finding. Pointing further to the involvement of VGSCs in schizophrenia, we found that these genes were enriched for nonsynonymous mutations (MAF ⩽ 0.1%) in cases genotyped using an exome array, (5585 schizophrenia cases and 8103 controls), and that in the trios data, synaptic proteins interacting with VGSCs were also enriched for both compound heterozygosity (P = 0.018) and de novo mutations (P = 0.04). However, we were unable to replicate the specific association with compound heterozygosity at VGSCs in an independent sample of Taiwanese schizophrenia trios (N = 614). We conclude that recessive genotypes do not appear to make a substantial contribution to schizophrenia at a genome-wide level. Although multiple lines of evidence, including several from this study, suggest that rare mutations in VGSCs contribute to the disorder, in the absence of replication of the original findings regarding compound heterozygosity, this conclusion requires evaluation in a larger sample of trios.

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.

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.

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.

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.

Mutation screening of the 3q29 microdeletion syndrome candidate genes DLG1 and PAK2 in schizophrenia.

Deletion of chromosome 3q29, which is associated with mental retardation and autism, was recently identified as being present in excess or occurring de novo in schizophrenia cases, being present in approximately 1/1,000 cases and 1/40,000 unscreened controls. Of the ∼20 genes in the commonly deleted region two are prominent candidates for involvement in the behavioral features of the microdeletion syndrome: DLG1 and PAK2. We report the result of mutation screening of the entire protein coding sequence of both genes in a sample of 234 unrelated cases and 272 unrelated controls from the UK. We find no evidence for any amino acid changing genetic variants in PAK2. We observe several rare and singleton non-synonymous genetic variations at DLG1, however there is no excess of these variants in cases when compared to controls. Our sample was underpowered to detect very rare or low-penetrance disease relevant alleles in the studied genes. Therefore very rare, low-to-moderate penetrance protein coding mutations or non-coding mutations at DLG1 and/or PAK2, or a nearby gene, may reproduce the behavioral characteristics of the 3q29 microdeletion.

Genome-wide association study on bipolar disorder in the Bulgarian population.

Bipolar disorder is a severe psychiatric disorder influenced by environmental and genetic factors. Genetic studies have implicated many variants in the disease's etiology but only few have been successfully replicated. We conducted a genome-wide association study (GWAS) on bipolar disorder in the Bulgarian population followed by a replication study of the top 100 single nucleotide polymorphisms (SNPs) showing the smallest P values. The GWAS was performed on 188 bipolar disorder patients and 376 control subjects genotyped on the Illumina 550 platform. The replication study was conducted on 122 patients and 328 controls. Although our study did not show any association P value that achieved genome-wide significance, and none of the top 100 SNPs reached the Bonferroni-corrected P value in the replication study, the plausible involvement of some variants cannot be entirely discarded. Three polymorphisms, rs8099939 [P = 2.12 × 10(-6), odds ratio (OR) = 1.95, 95% confidence interval (CI) = 1.43-2.67] in GRIK5, rs6122972 (P = 3.11 × 10(-6), OR = 2.02, 95% CI = 1.46-2.80) in PARD6B and rs2289700 (P = 9.14 × 10(-6), OR = 2.13, 95% CI = 1.53-2.95) in CTSH remained associated at a similar level after Mantel-Haenszel test for combining the results from the genome-wide and replication studies. A modest association was also detected for SNP rs1012053 (GWAS P = 4.50 × 10(-2)) in DGKH, which has already been reported as the most significant variant in a previous genome-wide scan on bipolar disorder. However, further studies using larger datasets are needed to identify variants with smaller effects that contribute to the risk of bipolar disorder.

DISC1 exon 11 rare variants found more commonly in schizoaffective spectrum cases than controls.

We previously performed a linkage study using families identified through probands meeting criteria for DSM-IV schizoaffective disorder, bipolar type (SABP) and observed a genome-wide significant signal (LOD = 3.54) at chromosome 1q42 close to DISC1. An initial sequencing study of DISC1 using 14 unrelated DSM-IV SABP samples from the linkage study identified 2 non-synonymous coding SNPs in exon 11 in 2 separate individuals. Here we provide evidence of additional rare coding SNPs within exon 11. In sequencing exon 11 in 506 cases and 1,211 controls for variants that occurred only once, 4 additional rare variants were found in cases (P-value = 0.008, Fisher's exact trend test).