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.

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.

Association at SYNE1 in both bipolar disorder and recurrent major depression.

Genome-wide association studies (GWAS) have identified a number of loci that have strong support for their association with bipolar disorder (BD). The Psychiatric Genome-Wide Association Study (GWAS) Consortium Bipolar Disorder Working Group (PGC-BD) meta-analysis of BD GWAS data sets and replication samples identified evidence (P=6.7 × 10⁻7, odds ratio (OR)=1.147) of association with the risk of BD at the polymorphism rs9371601 within SYNE1, a gene which encodes nesprin-1. Here we have tested this polymorphism in an independent BD case (n=1527) and control (n=1579) samples, and find evidence for association (P=0.0095) with similar effect sizes to those previously observed in BD (allelic OR=1.148). In a combined (meta) analysis of PGC-BD data (both primary and replication data) and our independent BD samples, we found genome-wide significant evidence for association (P=2.9 × 10⁻8, OR=1.104). We have also examined the polymorphism in our recurrent unipolar depression cases (n=1159) and control (n=2592) sample, and found that the risk allele was associated with risk for recurrent major depression (P=0.032, OR=1.118). Our findings add to the evidence that association at this locus influences susceptibility to bipolar and unipolar mood disorders.

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.

Genome-wide significant associations in schizophrenia to ITIH3/4, CACNA1C and SDCCAG8, and extensive replication of associations reported by the Schizophrenia PGC.

The Schizophrenia Psychiatric Genome-Wide Association Study Consortium (PGC) highlighted 81 single-nucleotide polymorphisms (SNPs) with moderate evidence for association to schizophrenia. After follow-up in independent samples, seven loci attained genome-wide significance (GWS), but multi-locus tests suggested some SNPs that did not do so represented true associations. We tested 78 of the 81 SNPs in 2640 individuals with a clinical diagnosis of schizophrenia attending a clozapine clinic (CLOZUK), 2504 cases with a research diagnosis of bipolar disorder, and 2878 controls. In CLOZUK, we obtained significant replication to the PGC-associated allele for no fewer than 37 (47%) of the SNPs, including many prior GWS major histocompatibility complex (MHC) SNPs as well as 3/6 non-MHC SNPs for which we had data that were reported as GWS by the PGC. After combining the new schizophrenia data with those of the PGC, variants at three loci (ITIH3/4, CACNA1C and SDCCAG8) that had not previously been GWS in schizophrenia attained that level of support. In bipolar disorder, we also obtained significant evidence for association for 21% of the alleles that had been associated with schizophrenia in the PGC. Our study independently confirms association to three loci previously reported to be GWS in schizophrenia, and identifies the first GWS evidence in schizophrenia for a further three loci. Given the number of independent replications and the power of our sample, we estimate 98% (confidence interval (CI) 78-100%) of the original set of 78 SNPs represent true associations. We also provide strong evidence for overlap in genetic risk between schizophrenia and bipolar disorder.

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.

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.

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.

Analysis of 10 independent samples provides evidence for association between schizophrenia and a SNP flanking fibroblast growth factor receptor 2.

We and others have previously reported linkage to schizophrenia on chromosome 10q25-q26 but, to date, a susceptibility gene in the region has not been identified. We examined data from 3606 single-nucleotide polymorphisms (SNPs) mapping to 10q25-q26 that had been typed in a genome-wide association study (GWAS) of schizophrenia (479 UK cases/2937 controls). SNPs with P<0.01 (n=40) were genotyped in an additional 163 UK cases and those markers that remained nominally significant at P<0.01 (n=22) were genotyped in replication samples from Ireland, Germany and Bulgaria consisting of a total of 1664 cases with schizophrenia and 3541 controls. Only one SNP, rs17101921, was nominally significant after meta-analyses across the replication samples and this was genotyped in an additional six samples from the United States/Australia, Germany, China, Japan, Israel and Sweden (n=5142 cases/6561 controls). Across all replication samples, the allele at rs17101921 that was associated in the GWAS showed evidence for association independent of the original data (OR 1.17 (95% CI 1.06-1.29), P=0.0009). The SNP maps 85 kb from the nearest gene encoding fibroblast growth factor receptor 2 (FGFR2) making this a potential susceptibility gene for schizophrenia.

Genetic utility of broadly defined bipolar schizoaffective disorder as a diagnostic concept.

BACKGROUND: Psychiatric phenotypes are currently defined according to sets of descriptive criteria. Although many of these phenotypes are heritable, it would be useful to know whether any of the various diagnostic categories in current use identify cases that are particularly helpful for biological-genetic research. AIMS: To use genome-wide genetic association data to explore the relative genetic utility of seven different descriptive operational diagnostic categories relevant to bipolar illness within a large UK case-control bipolar disorder sample. METHOD: We analysed our previously published Wellcome Trust Case Control Consortium (WTCCC) bipolar disorder genome-wide association data-set, comprising 1868 individuals with bipolar disorder and 2938 controls genotyped for 276 122 single nucleotide polymorphisms (SNPs) that met stringent criteria for genotype quality. For each SNP we performed a test of association (bipolar disorder group v. control group) and used the number of associated independent SNPs statistically significant at P<0.00001 as a metric for the overall genetic signal in the sample. We next compared this metric with that obtained using each of seven diagnostic subsets of the group with bipolar disorder: Research Diagnostic Criteria (RDC): bipolar I disorder; manic disorder; bipolar II disorder; schizoaffective disorder, bipolar type; DSM-IV: bipolar I disorder; bipolar II disorder; schizoaffective disorder, bipolar type. RESULTS: The RDC schizoaffective disorder, bipolar type (v. controls) stood out from the other diagnostic subsets as having a significant excess of independent association signals (P<0.003) compared with that expected in samples of the same size selected randomly from the total bipolar disorder group data-set. The strongest association in this subset of participants with bipolar disorder was at rs4818065 (P = 2.42 x 10(-7)). Biological systems implicated included gamma amniobutyric acid (GABA)(A) receptors. Genes having at least one associated polymorphism at P<10(-4) included B3GALTS, A2BP1, GABRB1, AUTS2, BSN, PTPRG, GIRK2 and CDH12. CONCLUSIONS: Our findings show that individuals with broadly defined bipolar schizoaffective features have either a particularly strong genetic contribution or that, as a group, are genetically more homogeneous than the other phenotypes tested. The results point to the importance of using diagnostic approaches that recognise this group of individuals. Our approach can be applied to similar data-sets for other psychiatric and non-psychiatric phenotypes.

Susceptibility locus for Alzheimer's disease on chromosome 10.

The apolipoprotein E (APOE) gene is the only genetic risk factor that has so far been linked to risk for late-onset Alzheimer's disease (LOAD). However, 50 percent of Alzheimer's disease cases do not carry an APOE4 allele, suggesting that other risk factors must exist. We performed a two-stage genome-wide screen in sibling pairs with LOAD to detect other susceptibility loci. Here we report evidence for an Alzheimer's disease locus on chromosome 10. Our stage one multipoint lod score (logarithm of the odds ratio for linkage/no linkage) of 2.48 (266 sibling pairs) increased to 3.83 in stage 2 (429 sibling pairs) close to D10S1225 (79 centimorgans). This locus modifies risk for Alzheimer's disease independent of APOE genotype.

Alternative splicing of a human alpha-tropomyosin muscle-specific exon: identification of determining sequences.

The human alpha-tropomyosin gene hTMnm has two mutually exclusive versions of exon 5 (NM and SK), one of which is expressed specifically in skeletal muscle (exon SK). A minigene construct expresses only the nonmuscle (NM) isoform when transfected into COS-1 cells and both forms when transfected into myoblasts. Twenty-four mutants were produced to determine why the SK exon is not expressed in COS cells. The results showed that exons NM and SK are not in competition for splicing to the flanking exons and that there is no intrinsic barrier to splicing between the exons. Instead, exon SK is skipped whenever there are flanking introns. Splicing of exon SK was induced when the branch site sequence 70 nucleotides upstream of the exon was mutated to resemble the consensus and when the extremities of the exon itself were changed to the corresponding NM sequence. Precise swaps of the NM and SK exon sequences showed that the exon sequence effect was dominant to that of intron sequences. The mechanism of regulation appears to be unlike that of other tropomyosin genes. We propose that exclusion of exon SK arises because its 3' splicing signals are weak and are prevented by an exon-specific repressor from competing for splice site recognition.

Genome wide significant linkage in schizophrenia conditioning on occurrence of depressive episodes.

BACKGROUND: Schizophrenia shows substantial clinical heterogeneity. One common important clinical variable in presentation is the occurrence of episodes of major depression. METHODS: We undertook analyses in an attempt to detect loci that influence susceptibility to, or modify the clinical expression of, schizophrenia according to the occurrence of episodes of major depression. We used a logistic regression framework in which lifetime presence/absence of major depression was entered as a covariate in the linkage analysis of our UK schizophrenia affected sibling pair series (168 affected sibling pairs typed for a 10 cM map of microsatellite markers). RESULTS: Inclusion of presence/absence of depression as a covariate detected a genome wide significant linkage signal on chromosome 4q28.3 at 130.7 cM (LOD = 4.59; p = 0.038; increase in maximum LOD over univariate analysis (ILOD) = 3.62). Inclusion of the depression covariate also showed suggestive evidence of linkage on 20q11.21 (LOD = 4.10; expected to occur by chance 0.093 times per genome scan, ILOD = 2.83). CONCLUSIONS: Our findings identify loci that may harbour genes that play a role in susceptibility to, or modify the risk of, episodes of major depression in people with schizophrenia.

In vivo co-localisation of MBNL protein with DMPK expanded-repeat transcripts.

Myotonic dystrophy (DM1) is the most common form of adult muscular dystrophy and is inherited as an autosomal dominant trait. The genetic basis of DM1 is the expansion of a CTG repeat in the 3' untranslated region of a protein kinase gene (DMPK). The molecular mechanism by which this expanded repeat produces the pathophysiology of DM1 remains unknown. Transcripts from the expanded allele accumulate as foci in the nucleus of DM1 cells and it has been suggested that these transcript foci sequester cellular proteins that are required for normal nuclear function. We have investigated the role of three RNA-binding proteins, CUG-BP, hnRNP C and MBNL, as possible sequestered factors. Using a combination of indirect immunofluorescence to detect endogenous proteins and overexpression of proteins with green fluorescent protein (GFP) tags we have shown that CUG-BP and hnRNP C do not co-localise with expanded repeat foci in DM1 cell lines. However, GFP-tagged MBNL does itself form foci in DM1 cell lines and co-localises with the foci of expanded repeat transcripts. GFP-tagged MBNL does not appear as foci in non-DM1 cell lines. This work provides further support for the involvement of MBNL in DM1.

A genome-wide map showing common regions of loss of heterozygosity/allelic imbalance in breast cancer.

We report here a new framework map that integrates data from 143 studies on the loss of heterozygosity/allelic imbalance in breast cancer. Full details of these loss of heterozygosity maps can be found at the web site (http://www.nottingham.ac.uk/pdzmgh/loh/) that accompanies this report. By combining results from these data, we have also been able to highlight and identify minimum commonly deleted regions on each chromosome. In addition to finding commonly deleted regions at both the BRCA1 and BRCA2 loci, which confirmed the power of the technique, 24 other regions were identified on 16 different chromosomes.

A genomewide linkage study of age at onset in schizophrenia.

There is strong evidence for a genetic contribution to age at onset of schizophrenia, which probably involves both susceptibility loci for schizophrenia and modifying loci acting independent of disease risk. We sought evidence of linkage to loci that influence age at onset of schizophrenia in a sample of 94 affected sibling pairs with DSM-IV schizophrenia or schizoaffective disorder, and age at first psychiatric contact of 45 years or less. Individuals were genotyped for 229 microsatellite markers spaced at approximately 20 cM intervals throughout the genome. Loci contributing to age at onset were sought by a quantitative maximum-likelihood multipoint linkage analysis using MAPMAKER/SIBS. A nonparametric multipoint analysis was also performed. The genomewide significance of linkage results was assessed by simulation studies. There were six maximum-likelihood LOD score peaks of 1.5 or greater, the highest being on chromosome 17q (LOD = 2.54; genomewide P = 0.27). This fulfils Lander and Kruglyak's [1995: Nat Genet 11:241-247] criteria for suggestive linkage in that it would be expected to occur once or less (0.3 times) per genome scan. However, this finding should be treated with caution because the LOD score appeared to be almost solely accounted for by the pattern of ibd sharing at one marker (D17S787), with virtually no evidence of linkage over flanking markers. None of the linkage results achieved genomewide statistical significance, but the LOD score peak on chromosome 13q (LOD = 1.68) coincided with the region showing maximum evidence for linkage in the study by Blouin et al. [1998: Nat Genet 20:70-73] of categorical schizophrenia.