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

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

Exome arrays capture polygenic rare variant contributions to schizophrenia.

Schizophrenia is a highly heritable disorder. Genome-wide association studies based largely on common alleles have identified over 100 schizophrenia risk loci, but it is also evident from studies of copy number variants (CNVs) and from exome-sequencing studies that rare alleles are also involved. Full characterization of the contribution of rare alleles to the disorder awaits the deployment of sequencing technology in very large sample sizes, meanwhile, as an interim measure, exome arrays allow rare non-synonymous variants to be sampled at a fraction of the cost. In an analysis of exome array data from 13 688 individuals (5585 cases and 8103 controls) from the UK, we found that rare (minor allele frequency < 0.1%) variant association signal was enriched among genes that map to autosomal loci that are genome-wide significant (GWS) in common variant studies of schizophrenia genome-wide association study (PGWAS = 0.01) as well as gene sets known to be enriched for rare variants in sequencing studies (PRARE = 0.026). We also identified the gene-wise equivalent of GWS support for WDR88 (WD repeat-containing protein 88), a gene of unknown function (P = 6.5 × 10(-7)). Rare alleles represented on exome chip arrays contribute to the genetic architecture of schizophrenia, but as is the case for GWAS, very large studies are required to reveal additional susceptibility alleles for the disorder.

Functional analysis of rare variants found in schizophrenia implicates a critical role for GIT1-PAK3 signaling in neuroplasticity.

Although the pathogenesis of schizophrenia (SCZ) is proposed to involve alterations of neural circuits via synaptic dysfunction, the underlying molecular mechanisms remain poorly understood. Recent exome sequencing studies of SCZ have uncovered numerous single-nucleotide variants (SNVs); however, the majority of these SNVs have unknown functional consequences, leaving their disease relevance uncertain. Filling this knowledge gap requires systematic application of quantitative and scalable assays to assess known and novel biological functions of genes. Here we demonstrate loss-of-function effects of multiple rare coding SNVs found in SCZ subjects in the GIT1 (G protein-coupled receptor kinase interacting ArfGAP 1) gene using functional cell-based assays involving coexpression of GIT1 and PAK3 (p21 protein (Cdc42/Rac)-activated kinase 3). Most notably, a GIT1-R283W variant reported in four independent SCZ cases was defective in activating PAK3 as well as MAPK (mitogen-activated protein kinase). Similar functional deficits were found for a de novo SCZ variant GIT1-S601N. Additional assays revealed deficits in the capacity of GIT1-R283W to stimulate PAK phosphorylation in cultured hippocampal neurons. In addition, GIT1-R283W showed deficits in the induction of GAD1 (glutamate decarboxylase 1) protein expression. Extending these functional assays to 10 additional rare GIT1 variants revealed the existence of an allelic series with the majority of the SCZ case variants exhibiting loss of function toward MAPK activation in a manner correlated with loss of PAK3 activation. Taken together, we propose that rare variants in GIT1, along with other genetic and environmental factors, cause dysregulation of PAK3 leading to synaptic deficits in SCZ.

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.

Detecting large copy number variants using exome genotyping arrays in a large Swedish schizophrenia sample.

Although copy number variants (CNVs) are important in genomic medicine, CNVs have not been systematically assessed for many complex traits. Several large rare CNVs increase risk for schizophrenia (SCZ) and autism and often demonstrate pleiotropic effects; however, their frequencies in the general population and other complex traits are unknown. Genotyping large numbers of samples is essential for progress. Large cohorts from many different diseases are being genotyped using exome-focused arrays designed to detect uncommon or rare protein-altering sequence variation. Although these arrays were not designed for CNV detection, the hybridization intensity data generated in each experiment could, in principle, be used for gene-focused CNV analysis. Our goal was to evaluate the extent to which CNVs can be detected using data from one particular exome array (the Illumina Human Exome Bead Chip). We genotyped 9100 Swedish subjects (3962 cases with SCZ and 5138 controls) using both standard genome-wide association study (GWAS) and exome arrays. In comparison with CNVs detected using GWAS arrays, we observed high sensitivity and specificity for detecting genic CNVs 400 kb including known pathogenic CNVs along with replicating the literature finding that cases with SCZ had greater enrichment for genic CNVs. Our data confirm the association of SCZ with 16p11.2 duplications and 22q11.2 deletions, and suggest a novel association with deletions at 11q12.2. Our results suggest the utility of exome-focused arrays in surveying large genic CNVs in very large samples; and thereby open the door for new opportunities such as conducting well-powered CNV assessment and comparisons between different diseases. The use of a single platform also minimizes potential confounding factors that could impact accurate detection.

Thy-1 is a component common to multiple populations of synaptic vesicles.

Thy-1, a glycosylphosphatidylinositol-linked integral membrane protein of the immunoglobulin superfamily, is a component of both large dense-core and small clear vesicles in PC12 cells. A majority of this protein, formerly recognized only on the plasma membrane of neurons, is localized to regulated secretory vesicles. Thy-1 is also present in synaptic vesicles in rat central nervous system. Experiments on permeabilized PC12 cells demonstrate that antibodies against Thy-1 inhibit the regulated release of neurotransmitter; this inhibition appears to be independent of any effect on the Ca2+ channel. These findings suggest Thy-1 is an integral component of many types of regulated secretory vesicles, and plays an important role in the regulated vesicular release of neurotransmitter at the synapse.

Polygenic risk for schizophrenia and neurocognitive performance in patients with schizophrenia.

Both neurocognitive deficits and schizophrenia are highly heritable. Genetic overlap between neurocognitive deficits and schizophrenia has been observed in both the general population and in the clinical samples. This study aimed to examine if the polygenic architecture of susceptibility to schizophrenia modified neurocognitive performance in schizophrenia patients. Schizophrenia polygenic risk scores (PRSs) were first derived from the Psychiatric Genomics Consortium (PGC) on schizophrenia, and then the scores were calculated in our independent sample of 1130 schizophrenia trios, who had PsychChip data and were part of the Schizophrenia Families from Taiwan project. Pseudocontrols generated from the nontransmitted parental alleles of the parents in these trios were compared with alleles in schizophrenia patients in assessing the replicability of PGC-derived susceptibility variants. Schizophrenia PRS at the P-value threshold (PT) of 0.1 explained 0.2% in the variance of disease status in this Han-Taiwanese samples, and the score itself had a P-value 0.05 for the association test with the disorder. Each patient underwent neurocognitive evaluation on sustained attention using the continuous performance test and executive function using the Wisconsin Card Sorting Test. We applied a structural equation model to construct the neurocognitive latent variable estimated from multiple measured indices in these 2 tests, and then tested the association between the PRS and the neurocognitive latent variable. Higher schizophrenia PRS generated at the PT of 0.1 was significantly associated with poorer neurocognitive performance with explained variance 0.5%. Our findings indicated that schizophrenia susceptibility variants modify the neurocognitive performance in schizophrenia patients.

Genome-wide association studies of smooth pursuit and antisaccade eye movements in psychotic disorders: findings from the B-SNIP study.

Eye movement deviations, particularly deficits of initial sensorimotor processing and sustained pursuit maintenance, and antisaccade inhibition errors, are established intermediate phenotypes for psychotic disorders. We here studied eye movement measures of 849 participants from the Bipolar-Schizophrenia Network on Intermediate Phenotypes (B-SNIP) study (schizophrenia N=230, schizoaffective disorder N=155, psychotic bipolar disorder N=206 and healthy controls N=258) as quantitative phenotypes in relation to genetic data, while controlling for genetically derived ancestry measures, age and sex. A mixed-modeling genome-wide association studies approach was used including ~4.4 million genotypes (PsychChip and 1000 Genomes imputation). Across participants, sensorimotor processing at pursuit initiation was significantly associated with a single nucleotide polymorphism in IPO8 (12p11.21, P=8 × 10-11), whereas suggestive associations with sustained pursuit maintenance were identified with SNPs in SH3GL2 (9p22.2, P=3 × 10-8). In participants of predominantly African ancestry, sensorimotor processing was also significantly associated with SNPs in PCDH12 (5q31.3, P=1.6 × 10-10), and suggestive associations were observed with NRSN1 (6p22.3, P=5.4 × 10-8) and LMO7 (13q22.2, P=7.3x10-8), whereas antisaccade error rate was significantly associated with a non-coding region at chromosome 7 (P=6.5 × 10-9). Exploratory pathway analyses revealed associations with nervous system development and function for 40 top genes with sensorimotor processing and pursuit maintenance (P=4.9 × 10-2-9.8 × 10-4). Our findings suggest novel patterns of genetic variation relevant for brain systems subserving eye movement control known to be impaired in psychotic disorders. They include genes involved in nuclear trafficking and gene silencing (IPO8), fast axonal guidance and synaptic specificity (PCDH12), transduction of nerve signals (NRSN1), retinal degeneration (LMO7), synaptic glutamate release (SH3GL2), and broader nervous system development and function.

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.

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.