Abstract
Schizophrenia is a highly heritable disorder. Genetic risk is conferred by a large number of alleles, including common alleles of small effect that might be detected by genome-wide association studies. Here we report a multi-stage schizophrenia genome-wide association study of up to 36,989 cases and 113,075 controls. We identify 128 independent associations spanning 108 conservatively defined loci that meet genome-wide significance, 83 of which have not been previously reported. Associations were enriched among genes expressed in brain, providing biological plausibility for the findings. Many findings have the potential to provide entirely new insights into aetiology, but associations at DRD2 and several genes involved in glutamatergic neurotransmission highlight molecules of known and potential therapeutic relevance to schizophrenia, and are consistent with leading pathophysiological hypotheses. Independent of genes expressed in brain, associations were enriched among genes expressed in tissues that have important roles in immunity, providing support for the speculated link between the immune system and schizophrenia.
Similar content being viewed by others
Main
Schizophrenia has a lifetime risk of around 1%, and is associated with substantial morbidity and mortality as well as personal and societal costs1,2,3. Although pharmacological treatments are available for schizophrenia, their efficacy is poor for many patients4. All available antipsychotic drugs are thought to exert their main therapeutic effects through blockade of the type 2 dopaminergic receptor5,6 but, since the discovery of this mechanism over 60 years ago, no new antipsychotic drug of proven efficacy has been developed based on other target molecules. Therapeutic stasis is in large part a consequence of the fact that the pathophysiology of schizophrenia is unknown. Identifying the causes of schizophrenia is therefore a critical step towards improving treatments and outcomes for those with the disorder.
High heritability points to a major role for inherited genetic variants in the aetiology of schizophrenia7,8. Although risk variants range in frequency from common to extremely rare9, estimates10,11 suggest half to a third of the genetic risk of schizophrenia is indexed by common alleles genotyped by current genome-wide association study (GWAS) arrays. Thus, GWAS is potentially an important tool for understanding the biological underpinnings of schizophrenia.
To date, around 30 schizophrenia-associated loci10,11,12,13,14,15,16,17,18,19,20,21,22,23 have been identified through GWAS. Postulating that sample size is one of the most important limiting factors in applying GWAS to schizophrenia, we created the Schizophrenia Working Group of the Psychiatric Genomics Consortium (PGC). Our primary aim was to combine all available schizophrenia samples with published or unpublished GWAS genotypes into a single, systematic analysis24. Here we report the results of that analysis, including at least 108 independent genomic loci that exceed genome-wide significance. Some of the findings support leading pathophysiological hypotheses of schizophrenia or targets of therapeutic relevance, but most of the findings provide new insights.
108 independent associated loci
We obtained genome-wide genotype data from which we constructed 49 ancestry matched, non-overlapping case-control samples (46 of European and three of east Asian ancestry, 34,241 cases and 45,604 controls) and 3 family-based samples of European ancestry (1,235 parent affected-offspring trios) (Supplementary Table 1 and Supplementary Methods). These comprise the primary PGC GWAS data set. We processed the genotypes from all studies using unified quality control procedures followed by imputation of SNPs and insertion-deletions using the 1000 Genomes Project reference panel25. In each sample, association testing was conducted using imputed marker dosages and principal components (PCs) to control for population stratification. The results were combined using an inverse-variance weighted fixed effects model26. After quality control (imputation INFO score ≥ 0.6, MAF ≥ 0.01, and successfully imputed in ≥ 20 samples), we considered around 9.5 million variants. The results are summarized in Fig. 1. To enable acquisition of large samples, some groups ascertained cases via clinician diagnosis rather than a research-based assessment and provided evidence of the validity of this approach (Supplementary Information)11,13. Post hoc analyses revealed the pattern of effect sizes for associated loci was similar across different assessment methods and modes of ascertainment (Extended Data Fig. 1), supporting our a priori decision to include samples of this nature.
Manhattan plot of the discovery genome-wide association meta-analysis of 49 case control samples (34,241 cases and 45,604 controls) and 3 family based association studies (1,235 parent affected-offspring trios). The x axis is chromosomal position and the y axis is the significance (–log10 P; 2-tailed) of association derived by logistic regression. The red line shows the genome-wide significance level (5 × 10−8). SNPs in green are in linkage disequilibrium with the index SNPs (diamonds) which represent independent genome-wide significant associations.
For the subset of linkage-disequilibrium-independent single nucleotide polymorphisms (SNPs) with P < 1 × 10−6 in the meta-analysis, we next obtained results from deCODE genetics (1,513 cases and 66,236 controls of European ancestry). We define linkage-disequilibrium-independent SNPs as those with low linkage disequilibrium (r2 < 0.1) to a more significantly associated SNP within a 500-kb window. Given high linkage disequilibrium in the extended major histocompatibility complex (MHC) region spans ∼8 Mb, we conservatively include only a single MHC SNP to represent this locus. The deCODE data were then combined with those from the primary GWAS to give a data set of 36,989 cases and 113,075 controls. In this final analysis, 128 linkage-disequilibrium-independent SNPs exceeded genome-wide significance (P ≤ 5 × 10−8) (Supplementary Table 2).
As in meta-analyses of other complex traits which identified large numbers of common risk variants27,28, the test statistic distribution from our GWAS deviates from the null (Extended Data Fig. 2). This is consistent with the previously documented polygenic contribution to schizophrenia10,11. The deviation in the test statistics from the null (λGC = 1.47, λ1000 = 1.01) is only slightly less than expected (λGC = 1.56) under a polygenic model given fully informative genotypes, the current sample size, and the lifetime risk and heritability of schizophrenia29.
Additional lines of evidence allow us to conclude the deviation between the observed and null distributions in our primary GWAS indicates a true polygenic contribution to schizophrenia. First, applying a novel method30 that uses linkage disequilibrium information to distinguish between the major potential sources of test statistic inflation, we found our results are consistent with polygenic architecture but not population stratification (Extended Data Fig. 3). Second, the schizophrenia-associated alleles at 78% of 234 linkage-disequilibrium-independent SNPs exceeding P < 1 × 10−6 in the case-control GWAS were again overrepresented in cases in the independent samples from deCODE. This degree of consistency between the case-control GWAS and the replication data is highly unlikely to occur by chance (P = 6 × 10−19). The tested alleles surpassed the P < 10−6 threshold in our GWAS before we added either the trios or deCODE data to the meta-analysis. This trend test is therefore independent of the primary case-control GWAS. Third, analysing the 1,235 parent-proband trios, we again found excess transmission of the schizophrenia-associated allele at 69% of the 263 linkage-disequilibrium-independent SNPs with P < 1 × 10−6 in the case-control GWAS. This is again unlikely to occur by chance (P = 1 × 10−9) and additionally excludes population stratification as fully explaining the associations reaching our threshold for seeking replication. Fourth, we used the trios trend data to estimate the expected proportion of true associations at P < 1 × 10−6 in the discovery GWAS, allowing for the fact that half of the index SNPs are expected to show the same allelic trend in the trios by chance, and that some true associations will show opposite trends given the limited number of trio samples (Supplementary Methods). Given the observed trend test results, around 67% (95% confidence interval: 64–73%) or n = 176 of the associations in the scan at P < 1 × 10−6 are expected to be true, and therefore the number of associations that will ultimately be validated from this set of SNPs will be considerably more than those that now meet genome-wide significance. Taken together, these analyses indicate that the observed deviation of test statistics from the null primarily represents polygenic association signal and the considerable excess of associations at the tail of extreme significance largely correspond to true associations.
Independently associated SNPs do not translate to well-bounded chromosomal regions. Nevertheless, it is useful to define physical boundaries for the SNP associations to identify candidate risk genes. We defined an associated locus as the physical region containing all SNPs correlated at r2 > 0.6 with each of the 128 index SNPs. Associated loci within 250 kb of each other were merged. This resulted in 108 physically distinct associated loci, 83 of which have not been previously implicated in schizophrenia and therefore harbour potential new biological insights into disease aetiology (Supplementary Table 3; regional plots in Supplementary Fig. 1). The significant regions include all but 5 loci previously reported to be genome-wide significant in large samples (Supplementary Table 3).
Characterization of associated loci
Of the 108 loci, 75% include protein-coding genes (40%, a single gene) and a further 8% are within 20 kb of a gene (Supplementary Table 3). Notable associations relevant to major hypotheses of the aetiology and treatment of schizophrenia include DRD2 (the target of all effective antipsychotic drugs) and many genes (for example, GRM3, GRIN2A, SRR, GRIA1) involved in glutamatergic neurotransmission and synaptic plasticity. In addition, associations at CACNA1C, CACNB2 and CACNA1I, which encode voltage-gated calcium channel subunits, extend previous findings implicating members of this family of proteins in schizophrenia and other psychiatric disorders11,13,31,32. Genes encoding calcium channels, and proteins involved in glutamatergic neurotransmission and synaptic plasticity have been independently implicated in schizophrenia by studies of rare genetic variation33,34,35, suggesting convergence at a broad functional level between studies of common and rare genetic variation. We highlight in the Supplementary Discussion genes of particular interest within associated loci with respect to current hypotheses of schizophrenia aetiology or treatment (although we do not imply that these genes are necessarily the causal elements).
For each of the schizophrenia-associated loci, we identified a credible causal set of SNPs (for definition, see Supplementary Methods)36. In only 10 instances (Supplementary Table 4) was the association signal credibly attributable to a known non-synonymous exonic polymorphism. The apparently limited role of protein-coding variants is consistent both with exome sequencing findings33 and with the hypothesis that most associated variants detected by GWAS exert their effects through altering gene expression rather than protein structure37,38 and with the observation that schizophrenia risk loci are enriched for expression quantitative trait loci (eQTL)39.
To try to identify eQTLs that could explain associations with schizophrenia, we merged the credible causal set of SNPs defined above with eQTLs from a meta-analysis of human brain cortex eQTL studies (n = 550) and an eQTL study of peripheral venous blood (n = 3,754)40 (Supplementary Methods). Multiple schizophrenia loci contained at least one eQTL for a gene within 1 Mb of the locus (Supplementary Table 4). However, in only 12 instances was the eQTL plausibly causal (two in brain, and nine in peripheral blood, one in both). This low proportion suggests that if most risk variants are regulatory, available eQTL catalogues do not yet provide power, cellular specificity, or developmental diversity to provide clear mechanistic hypotheses for follow-up experiments.
The brain and immunity
To further explore the regulatory nature of the schizophrenia associations, we mapped the credible sets (n = 108) of causal variants onto sequences with epigenetic markers characteristic of active enhancers in 56 different tissues and cell lines (Supplementary Methods). Schizophrenia associations were significantly enriched at enhancers active in brain (Fig. 2) but not in tissues unlikely to be relevant to schizophrenia (for example, bone, cartilage, kidney and fibroblasts). Brain tissues used to define enhancers consist of heterogeneous populations of cells. Seeking greater specificity, we contrasted genes enriched for expression in neurons and glia using mouse ribotagged lines41. Genes with strong expression in multiple cortical and striatal neuronal lineages were enriched for associations, providing support for an important neuronal pathology in schizophrenia (Extended Data Fig. 4) but this is not statistically more significant than, or exclusionary of, contributions from other lineages42.
Cell and tissue type specific enhancers were identified using ChIP-seq data sets (H3K27ac signal) from 56 cell line and tissue samples (y axis). We defined cell and tissue type enhancers as the top 10% of enhancers with the highest ratio of reads in that cell or tissue type divided by the total number of reads. Enrichment of credible causal associated SNPs from the schizophrenia GWAS was compared with frequency matched sets of 1000 Genomes SNPs (Supplementary Methods). The x axis is the –log10 P for enrichment. P values are uncorrected for the number of tissues or cells tested. A –log10 P of roughly 3 can be considered significant after Bonferroni correction. Descriptions of cell and tissue types at the Roadmap Epigenome website (http://www.roadmapepigenomics.org).
Schizophrenia associations were also strongly enriched at enhancers that are active in tissues with important immune functions, particularly B-lymphocyte lineages involved in acquired immunity (CD19 and CD20 lines, Fig. 2). These enrichments remain significant even after excluding the extended MHC region and regions containing brain enhancers (enrichment P for CD20 < 10−6), demonstrating that this finding is not an artefact of correlation between enhancer elements in different tissues and not driven by the strong and diffuse association at the extended MHC. Epidemiological studies have long hinted at a role for immune dysregulation in schizophrenia, the present findings provide genetic support for this hypothesis43.
To develop additional biological hypotheses beyond those that emerge from inspection of the individual loci, we further undertook a limited mining of the data through gene-set analysis. However, as there is no consensus methodology by which such analyses should be conducted, nor an established optimal significance threshold for including loci, we sought to be conservative, using only two of the many available approaches44,45 and restricting analyses to genes within genome-wide significant loci. Neither approach identified gene-sets that were significantly enriched for associations after correction for the number of pathways tested (Supplementary Table 5) although nominally significantly enrichments were observed among several predefined candidate pathways (Extended Data Table 1). A fuller exploratory analysis of the data will be presented elsewhere.
Overlap with rare mutations
CNVs associated with schizophrenia overlap with those associated with autism spectrum disorder (ASD) and intellectual disability9, as do genes with deleterious de novo mutations34. Here we find significant overlap between genes in the schizophrenia GWAS associated intervals and those with de novo non-synonymous mutations in schizophrenia (P = 0.0061) (Extended Data Table 2), suggesting that mechanistic studies of rare genetic variation in schizophrenia will be informative for schizophrenia more widely. We also find evidence for overlap between genes in schizophrenia GWAS regions and those with de novo non-synonymous mutations in intellectual disability (P = 0.00024) and ASD (P = 0.035), providing further support for the hypothesis that these disorders have partly overlapping pathophysiologies9,34.
Polygenic risk score profiling
Previous studies have shown that risk profile scores (RPS) constructed from alleles showing modest association with schizophrenia in a discovery GWAS can predict case-control status in independent samples, albeit with low sensitivity and specificity10,11,16. This finding was robustly confirmed in the present study. The estimate of Nagelkerke R2 (a measure of variance in case-control status explained) depends on the specific target data set and threshold (PT) for selecting risk alleles for RPS analysis (Extended Data Fig. 5 and 6a). However, using the same target sample as earlier studies and PT = 0.05, R2 is now increased from 0.03 (ref. 10) to 0.184 (Extended Data Fig. 5). Assuming a liability-threshold model, a lifetime risk of 1%, independent SNP effects, and adjusting for case-control ascertainment, RPS now explains about 7% of variation on the liability scale46 to schizophrenia across the samples (Extended Data Fig. 6b), about half of which (3.4%) is explained by genome-wide significant loci.
We also evaluated the capacity of RPS to predict case-control status using a standard epidemiological approach to a continuous risk factor. We illustrate this in three samples, each with different ascertainment schemes (Fig. 3). The Danish sample is population-based (that is, inpatient and outpatient facilities), the Swedish sample is based on all cases hospitalized for schizophrenia in Sweden, and the Molecular Genetics of Schizophrenia (MGS) sample was ascertained specially for genetic studies from clinical sources in the US and Australia. We grouped individuals into RPS deciles and estimated the odds ratios for affected status for each decile with reference to the lowest risk decile. The odds ratios increased with greater number of schizophrenia risk alleles in each sample, maximizing for the tenth decile in all samples: Denmark 7.8 (95% confidence interval (CI): 4.4–13.9), Sweden 15.0 (95% CI: 12.1–18.7) and MGS 20.3 (95% CI: 14.7–28.2). Given the need for measures that index liability to schizophrenia47,48, the ability to stratify individuals by RPS offers new opportunities for clinical and epidemiological research. Nevertheless, we stress that the sensitivity and specificity of RPS do not support its use as a predictive test. For example, in the Danish epidemiological sample, the area under the receiver operating curve is only 0.62 (Extended Data Fig. 6c, Supplementary Table 6).
Odds ratio for schizophrenia by risk score profile (RPS) decile in the Sweden (Sw1-6), Denmark (Aarhus), and Molecular Genetics of Schizophrenia studies (Supplementary Methods). Risk alleles and weights were derived from ‘leave one out’ analyses in which those samples were excluded from the GWAS meta-analysis (Supplementary Methods). The threshold for selecting risk alleles was PT < 0.05. The RPS were converted to deciles (1 = lowest, 10 = highest RPS), and nine dummy variables created to contrast deciles 2-10 to decile 1 as the reference. Odds ratios and 95% confidence intervals (bars) were estimated using logistic regression with PCs to control for population stratification.
Finally, seeking evidence for non-additive effects on risk, we tested for statistical interaction between all pairs of 125 autosomal SNPs that reached genome-wide significance. P values for the interaction terms were distributed according to the null, and no interaction was significant after correction for multiple comparisons. Thus, we find no evidence for epistatic or non-additive effects between the significant loci (Extended Data Fig. 7). It is possible that such effects could be present between other loci, or occur in the form of higher-order interactions.
Discussion
In the largest (to our knowledge) molecular genetic study of schizophrenia, or indeed of any neuropsychiatric disorder, ever conducted, we demonstrate the power of GWAS to identify large numbers of risk loci. We show that the use of alternative ascertainment and diagnostic schemes designed to rapidly increase sample size does not inevitably introduce a crippling degree of heterogeneity. That this is true for a phenotype like schizophrenia, in which there are no biomarkers or supportive diagnostic tests, provides grounds to be optimistic that this approach can be successfully applied to GWAS of other clinically defined disorders.
We further show that the associations are not randomly distributed across genes of all classes and function; rather they converge upon genes that are expressed in certain tissues and cellular types. The findings include molecules that are the current, or the most promising, targets for therapeutics, and point to systems that align with the predominant aetiological hypotheses of the disorder. This suggests that the many novel findings we report also provide an aetiologically relevant foundation for mechanistic and treatment development studies. We also find overlap between genes affected by rare variants in schizophrenia and those within GWAS loci, and broad convergence in the functions of some of the clusters of genes implicated by both sets of genetic variants, particularly genes related to abnormal glutamatergic synaptic and calcium channel function. How variation in these genes impact function to increase risk for schizophrenia cannot be answered by genetics, but the overlap strongly suggests that common and rare variant studies are complementary rather than antagonistic, and that mechanistic studies driven by rare genetic variation will be informative for schizophrenia.
References
Saha, S., Chant, D. & McGrath, J. A systematic review of mortality in schizophrenia: is the differential mortality gap worsening over time? Arch. Gen. Psychiatry 64, 1123–1131 (2007)
World Health Organization. The Global Burden of Disease: 2004 Update (WHO Press, 2008)
Knapp, M., Mangalore, R. & Simon, J. The global costs of schizophrenia. Schizophr. Bull. 30, 279–293 (2004)
Lieberman, J. A. et al. Effectiveness of antipsychotic drugs in patients with chronic schizophrenia. N. Engl. J. Med. 353, 1209–1223 (2005)
Carlsson, A. & Lindqvist, M. Effect of chlorpromazine or haloperidol on formation of 3-methoxytyramine and normetanephrine in mouse brain. Acta Pharmacol. Toxicol. 20, 140–144 (1963)
van Rossum, J. M. The significance of dopamine-receptor blockade for the mechanism of action of neuroleptic drugs. Arch. Int. Pharmacodyn. Ther. 160, 492–494 (1966)
Lichtenstein, P. et al. Recurrence risks for schizophrenia in a Swedish national cohort. Psychol. Med. 36, 1417–1425 (2006)
Sullivan, P. F., Kendler, K. S. & Neale, M. C. Schizophrenia as a complex trait: evidence from a meta-analysis of twin studies. Arch. Gen. Psychiatry 60, 1187–1192 (2003)
Sullivan, P. F., Daly, M. J. & O’Donovan, M. Genetic architectures of psychiatric disorders: the emerging picture and its implications. Nature Rev. Genet. 13, 537–551 (2012)
International Schizophrenia Consortium. Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature 460, 748–752 (2009)
Ripke, S. et al. Genome-wide association analysis identifies 13 new risk loci for schizophrenia. Nature Genet. 45, 1150–1159 (2013)
Ikeda, M. et al. Genome-wide association study of schizophrenia in a Japanese population. Biol. Psychiatry 69, 472–478 (2011)
Hamshere, M. L. et al. Genome-wide significant associations in schizophrenia to ITIH3/4, CACNA1C and SDCCAG8, and extensive replication of associations reported by the Schizophrenia PGC. Mol. Psychiatry 18, 708–712 (2013)
O’Donovan, M. C. et al. Identification of novel schizophrenia loci by genome-wide association and follow-up. Nature Genet. 40, 1053–1055 (2008)
Rietschel, M. et al. Association between genetic variation in a region on chromosome 11 and schizophrenia in large samples from Europe. Mol. Psychiatry 17, 906–917 (2012)
Schizophrenia Psychiatric Genome-Wide Association Study Consortium. Genome-wide association study identifies five new schizophrenia loci. Nature Genet. 43, 969–976 (2011)
Irish Schizophrenia Genomics Consortium & Wellcome Trust Case Control Consortium. Genome-wide association study implicates HLA-C*01:02 as a risk factor at the major histocompatibility complex locus in schizophrenia. Biol. Psychiatry 72, 620–628 (2012)
Shi, J. et al. Common variants on chromosome 6p22.1 are associated with schizophrenia. Nature 460, 753–757 (2009)
Shi, Y. et al. Common variants on 8p12 and 1q24.2 confer risk of schizophrenia. Nature Genet. 43, 1224–1227 (2011)
Stefansson, H. et al. Common variants conferring risk of schizophrenia. Nature 460, 744–747 (2009)
Steinberg, S. et al. Common variants at VRK2 and TCF4 conferring risk of schizophrenia. Hum. Mol. Genet. 20, 4076–4081 (2011)
Yue, W. H. et al. Genome-wide association study identifies a susceptibility locus for schizophrenia in Han Chinese at 11p11.2. Nature Genet. 43, 1228–1231 (2011)
Lencz, T. et al. Genome-wide association study implicates NDST3 in schizophrenia and bipolar disorder. Nature Commun. 4, 2739 (2013)
Psychiatric GWAS Consortium. A framework for interpreting genomewide association studies of psychiatric disorders. Mol. Psychiatry 14, 10–17 (2009)
The 1000 Genomes Project Consortium. A map of human genome variation from population-scale sequencing. Nature 467, 1061–1073 (2010)
Begum, F., Ghosh, D., Tseng, G. C. & Feingold, E. Comprehensive literature review and statistical considerations for GWAS meta-analysis. Nucleic Acids Res. 40, 3777–3784 (2012)
Lango Allen, H. et al. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature 467, 832–838 (2010)
Jostins, L. et al. Host–microbe interactions have shaped the genetic architecture of inflammatory bowel disease. Nature 491, 119–124 (2012)
Yang, J. et al. Genomic inflation factors under polygenic inheritance. Eur. J. Hum. Genet. 19, 807–812 (2011)
Bulik-Sullivan, B. K. et al. LD score regression distinguishes confounding from polygenicity in genome-wide association studies. Preprint at http://dx.doi.org/10.1101/002931 (2014)
Ferreira, M. A. et al. Collaborative genome-wide association supports a role for ANK3 and CACNA1C in bipolar disorder. Nature Genet. 40, 1056–1058 (2008)
Cross-Disorder Group of the Psychiatric Genomics Consortium. Identification of risk loci with shared effects on five major psychiatric disorders: a genome-wide analysis. Lancet 381, 1371–1379 (2013)
Purcell, S. M. et al. A polygenic burden of rare disruptive mutations in schizophrenia. Nature. 506, 185–190 (2014)
Fromer, M. et al. De novo mutations in schizophrenia implicate synaptic networks. Nature. 506, 179–184 (2014)
Kirov, G. et al. De novo CNV analysis implicates specific abnormalities of postsynaptic signalling complexes in the pathogenesis of schizophrenia. Mol. Psychiatry 17, 142–153 (2012)
Wellcome Trust Case Control Consortium Bayesian refinement of association signals for 14 loci in 3 common diseases. Nature Genet. 44, 1294–1301 (2012)
Nicolae, D. L. et al. Trait-associated SNPs are more likely to be eQTLs: annotation to enhance discovery from GWAS. PLoS Genet. 6, e1000888 (2010)
Maurano, M. T. et al. Systematic localization of common disease-associated variation in regulatory DNA. Science 337, 1190–1195 (2012)
Richards, A. L. et al. Schizophrenia susceptibility alleles are enriched for alleles that affect gene expression in adult human brain. Mol. Psychiatry 17, 193–201 (2012)
Wright, F. A. et al. Heritability and genomics of gene expression in peripheral blood. Nature Genet. 46, 430–437 (2014)
Doyle, J. P. et al. Application of a translational profiling approach for the comparative analysis of CNS cell types. Cell 135, 749–762 (2008)
Tkachev, D. et al. Oligodendrocyte dysfunction in schizophrenia and bipolar disorder. Lancet 362, 798–805 (2003)
Benros, M. E., Mortensen, P. B. & Eaton, W. W. Autoimmune diseases and infections as risk factors for schizophrenia. Ann. NY Acad. Sci. 1262, 56–66 (2012)
Holmans, P. et al. Gene ontology analysis of GWA study data sets provides insights into the biology of bipolar disorder. Am. J. Hum. Genet. 85, 13–24 (2009)
Lee, P. H., O’Dushlaine, C., Thomas, B. & Purcell, S. InRich: interval-based enrichment analysis for genome-wide association studies. Bioinformatics 28, 1797–1799 (2012)
Lee, S. H., Goddard, M. E., Wray, N. R. & Visscher, P. M. A better coefficient of determination for genetic profile analysis. Genet. Epidemiol. 36, 214–224 (2012)
Gottesman, I. I. & Gould, T. D. The endophenotype concept in psychiatry: etymology and strategic intentions. Am. J. Psychiatry 160, 636–645 (2003)
Insel, T. et al. Research domain criteria (RDoC): toward a new classification framework for research on mental disorders. Am. J. Psychiatry 167, 748–751 (2010)
Lips, E. S. et al. Functional gene group analysis identifies synaptic gene groups as risk factor for schizophrenia. Mol. Psychiatry 17, 996–1006 (2012)
Lewis, B. P., Burge, C. B. & Bartel, D. P. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 120, 15–20 (2005)
Acknowledgements
Core funding for the Psychiatric Genomics Consortium is from the US National Institute of Mental Health (U01 MH094421). We thank T. Lehner (NIMH). The work of the contributing groups was supported by numerous grants from governmental and charitable bodies as well as philanthropic donation. Details are provided in the Supplementary Notes. Membership of the Wellcome Trust Case Control Consortium and of the Psychosis Endophenotype International Consortium are provided in the Supplementary Notes.
Author information
Authors and Affiliations
Consortia
Contributions
The individual studies or consortia contributing to the GWAS meta-analysis were led by R.A., O.A.A., D.H.R.B., A.D.B., E. Bramon, J.D.B., A.C., D.A.C., S.C., A.D., E. Domenici, H.E., T.E., P.V.G., M.G., H.G., C.M.H., N.I., A.V.J., E.G.J., K.S.K., G.K., J. Knight, T. Lencz, D.F.L., Q.S.L., J. Liu, A.K.M., S.A.M., A. McQuillin, J.L.M., P.B.M., B.J.M., M.M.N., M.C.O’D., R.A.O., M.J.O., A. Palotie, C.N.P., T.L.P., M.R., B.P.R., D.R., P.C.S, P. Sklar. D.St.C., P.F.S., D.R.W., J.R.W., J.T.R.W. and T.W. Together with the core statistical analysis group led by M.J.D. comprising S.R., B.M.N. and P.A.H., this group comprised the management group led by M.C.O’D. who were responsible for the management of the study and the overall content of the manuscript. Additional analyses and interpretations were contributed by E.A., B.B.-S., D.K., K.-H.F., M. Fromer, H.H., P.L., P.B.M., S.M.P., T.H.P., N.R.W. and P.M.V. The phenotype supervisory group comprised A.C., A.H.F., P.V.G., K.K.K. and B.J.M. D.A.C. led the candidate selected genes subgroup comprised of M.J.D., E. Dominici, J.A.K., A.M.H., M.C.O’D, B.P.R., D.R., E.M.S. and P. Sklar. Replication results were provided by S.S., H.S. and K.S. The remaining authors contributed to the recruitment, genotyping, or data processing for the contributing components of the meta-analysis. A.C., M.J.D., B.M.N., S.R., P.F.S. and M.C.O’D. took responsibility for the primary drafting of the manuscript which was shaped by the management group. All other authors saw, had the opportunity to comment on, and approved the final draft.
Corresponding author
Ethics declarations
Competing interests
CFI statement–Several of the authors are employees of the following pharmaceutical companies; Pfizer (C.R.S., J.R.W., H.S.X.), F.Hoffman-La Roche (E.D., L.E.), Eli Lilly (D.A.C., Y.M., L.N.) and Janssen (S.G., D.W., Q.S.L.; also N.C. an ex-employee). Others are employees of deCODE genetics (S.S, H.S., K.S.). None of these companies influenced the design of the study, the interpretation of the data, or the amount of data reported, or financially profit by publication of the results which are pre-competitive. The other authors declare no competing interests.
Additional information
Results can be downloaded from the Psychiatric Genomics Consortium website (http://pgc.unc.edu) and visualized using Ricopili (http://www.broadinstitute.org/mpg/ricopili). Genotype data for the samples where the ethics permit deposition are available upon application from the NIMH Genetics Repository (https://www.nimhgenetics.org).
A list of authors and affiliations appear in the Supplementary Information.
Extended data figures and tables
Extended Data Figure 1 Homogeneity of effects across studies.
Plot of the first two principal components (PCs) from principal components analysis (PCA) of the logistic regression β coefficients for autosomal genome-wide significant associations. The input data were the β coefficients from 52 samples for 112 independent SNP associations (excluding 3 chrX SNPs and 13 SNPs with missing values in Asian samples). PCAs were weighted by the number of cases. Each circle shows the location of a study on PC1 and PC2. Circle size and colour are proportional to the number of cases in each sample (larger and darker red circles correspond to more cases). Most samples cluster. Outliers had either small numbers of cases (‘small’) or were genotyped on older arrays. Abbreviations: a500 (Affymetrix 500K); a5 (Affymetrix 5.0). Studies that did not use conventional research interviews are in the central cluster (CLOZUK, Sweden, and Denmark-Aarhus studies, see Supplementary Methods for sample descriptions).
Extended Data Figure 2 Quantile-quantile plot.
Quantile-quantile plot of the discovery genome-wide association meta-analysis of 49 case control samples (34,241 cases and 45,604 controls) and 3 family based association studies (1,235 parent affected-offspring trios). Expected –log10 P values are those expected under the null hypothesis. Observed are the GWAS association results derived by logistic regression (2-tailed) as in Fig. 1. For clarity, we avoided expansion of the y axis by setting the smallest association P values to 10−12. The shaded area surrounded by a red line indicates the 95% confidence interval under the null. λGC is the observed median χ2 test statistic divided by the median expected χ2 test statistic under the null hypothesis.
Extended Data Figure 3 Linkage disequilibrium score regression consistent with polygenic inheritance.
The relationship between marker χ2 association statistics and linkage disequilibrium (LD) as measured by the linkage disequilibrium score. Linkage disequilibrium score is the sum of the r2 values between a variant and all other known variants within a 1 cM window, and quantifies the amount of genetic variation tagged by that variant. Variants were grouped into 50 equal-sized bins based on linkage disequilibrium score rank. Linkage disequilibrium score bin and mean χ2 denotes mean linkage disequilibrium score and test statistic for markers each bin. a, b, We simulated (Supplementary Methods) test statistics under two scenarios: a, no true association, inflation due to population stratification; and b, polygenic inheritance (λ = 1.32), in which we assigned independent and identically distributed per-normalized-genotype effects to a randomly selected subset of variants. c, Results from the PGC schizophrenia GWAS (λ = 1.48). The real data are strikingly similar to the simulated data summarized in b but not a. The intercept estimates the inflation in the mean χ2 that results from confounding biases, such as cryptic relatedness or population stratification. Thus, the intercept of 1.066 for the schizophrenia GWAS suggests that ∼90% of the inflation in the mean χ2 results from polygenic signal. The results of the simulations are also consistent with theoretical expectation (see Supplementary Methods). λ is the median χ2 test statistic from the simulations (a, b) or the observed data (c) divided by the median expected χ2 test statistic under the null hypothesis.
Extended Data Figure 4 Enrichment of associations in tissues and cells.
Genes whose transcriptional start is nearest to the most associated SNP at each schizophrenia-associated locus were tested for enriched expression in purified brain cell subsets obtained from mouse ribotagged lines41 using enrichment analysis described in the Supplementary Methods. The red dotted line indicates P = 0.05.
Extended Data Figure 5 MGS risk profile score analysis.
Polygenic risk profile score (RPS) analyses using the MGS18 sample as target, and deriving risk alleles from three published schizophrenia data sets (x axis): ISC (2,615 cases and 3,338 controls)10, PGC1 (excluding MGS, 9,320 cases and 10,228 controls)16, and the current meta-analysis (excluding MGS) with 32,838 cases and 44,357 controls. Samples sizes differ slightly from the original publications due to different analytical procedures. This shows the increasing RPS prediction with increasing training data set size reflecting improved precision of estimates of the SNP effect sizes. The proportion of variance explained (y axis; Nagelkerke’s R2) was computed by comparison of a full model (covariates + RPS) score to a reduced model (covariates only). Ten different P value thresholds (PT) for selecting risk alleles are denoted by the colour of each bar (legend above plot). For significance testing, see the bottom legend which denotes the P value for the test that R2 is different from zero. All numerical data and methods used to generate these plots are available in Supplementary Table 6 and Supplementary Methods.
Extended Data Figure 6 Risk profile score analysis.
We defined 40 target subgroups of the primary GWAS data set and performed 40 leave-one-out GWAS analyses (see Supplementary Methods and Supplementary Table 7) from which we derived risk alleles for RPS analysis (x axis) for each target subgroup. a, The proportion of variance explained (y axis; Nagelkerke’s R2) was computed for each target by comparison of a full model (covariates + RPS) score to a reduced model (covariates only). For clarity, 3 different P value thresholds (PT) are presented denoted by the colour of each bar (legend above plot) as for Extended Data Fig. 5, but for clarity we restrict to fewer P value thresholds (PT of 5 × 10−8, 1 × 10−4 and 0.05) and removed the significance values. b, The proportion of variance on the liability scale from risk scores calculated at the PT 0.05 with 95% CI bar assuming baseline population disease risk of 1%. c, Area under the receiver operating curve (AUC). All numerical data and methods used to generate these plots are available in Supplementary Table 7 and Supplementary Methods.
Extended Data Figure 7 Pairwise epistasis analysis of significant SNPs.
Quantile-quantile plot for all pair-wise (n = 7,750) combinations of the 125 independent autosomal genome-wide significant SNPs tested for non-additive effects on risk using case-control data sets of European ancestry (32,405 cases and 42,221 controls). We included as covariates the principal components from the main analysis as well as a study indicator. The interaction model is described by: and
are genotypes at the two loci,
is the interaction between the two genotypes modelled in a multiplicative fashion,
is the vector of principal components,
is the vector of study indicator variables. Each
is the regression coefficient in the generalized linear model using logistic regression. The overall distribution of P values did not deviate from the null and the smallest P value (4.28 × 10−4) did not surpass the Bonferroni correction threshold (P = 0.05/7750 = 6.45 × 10−6). The line x = y indicates the expected null distribution with the grey area bounded by red lines indicating the expected 95% confidence interval for the null.
Supplementary information
Supplementary Information
This file contains Supplementary Text, Supplementary Tables 1-3, Supplementary References and Supplementary Notes (including a list of consortium members and acknowledgements) – see contents page for details. (PDF 1772 kb)
Supplementary Figure
This file contains Supplementary Figure 1. (PDF 2168 kb)
Supplementary Table 4
Credible causal schizophrenia SNPs, coding variants, and eQTLs. Worksheet 1: Coding variants: Index SNP is the schizophrenia associated SNP defining the schizophrenia associated region. Coding variant, R2, and gene denotes a coding credible SNP and the R2 with the index SNP, and the gene containing the coding variant. CHR (chromosome), BP (base position), A1A2 (alleles 1 and 2), frequencies of allele 1 (FRQ_A1), INFO (imputation quality) and P (P-value) refer to the index SNP in the discovery GWAS. P (incl rep) refers to replication P value for index SNP. Worksheets 2 and 3: Brain and blood eQTL: Credible SNP denotes a SNP within the schizophrenia credible set (defined in supplementary material) that is also a cis eQTL (transcript within 1Mb, PeQTL<1x10-4). P(cSCZ) is the schizophrenia (discovery) GWAS association P-value for the credible SNP. The Prob(cSCZ) is the normalized probability of the credible variant being causal for schizophrenia. N(cSCZ) is the number of variants in the credible set of schizophrenia variants within a region spanned by eQTLs at P<10-4. eQTL SNP is the most significant expression associated SNP in the region for the gene in next column (N.B., many regions have an eQTL for more than 1 gene). eQTLgene is the gene that is linked to the eQTL SNP. P(eQTL) is the association P-value between the eQTL SNP and the eQTLgene in the previous two columns. Prob(eQTL) is the normalized probability that the eQTL SNP is also the causal SNP for schizophrenia (high values mean higher probability of being causal). eQTLcumsum is the cumulative sum of the probability of all SNPs into the region, up to the inclusion of the max eQTL in locus ordered by probability of being the functional SNP. PeQTL(SCZ) is the schizophrenia association P-value for the eQTL SNP. R2 (cSCZ/ eQTL) is the R2 between the credibleSNP and eQTL SNP. Associations to schizophrenia that are plausibly explained by an eQTL are in bold. Separate worksheets provide information on brain and blood eQTL analyses. Distinct loci are alternately shaded/unshaded. (XLSX 106 kb)
Supplementary Table 5
Pathway analyses by ALIGATOR and INRICH. Enrichment analyses using ALIGATOR and INRICH were performed as described in Supplementary Text. Pathway ID denotes the pathway source: GO (Gene ontology; http://www.geneontology.org), KEGG (Kyoto Encyclopaedia of Genes and Genomes; http://www.genome.jp/kegg), PAN-PW (PANTHER; http://www.pantherdb.org/pathway), Reactome (http://www.reactome.org/download), BioCarta (downloaded from the Molecular Signatures Database v4.0 http://www.broadinstitute.org/gsea/msigdb/index.jsp), MGI (Mouse Genome Informatics; http://www.informatics.jax.org), and NCI pathways (NCI: http://pid.nci.nih.gov). (XLS 173 kb)
Supplementary Table 7
Risk Profile Score Analyses. Risk Profile Score (RPS) analysis was performed as described in supplementary text. RPS datasets tab provides the name given for sample in which RPS was performed (target label) and the datasets included (defined in Supplementary Table 1). The GWAS data used to define the risk alleles for RPS analysis represents the remaining GWAS samples. For various GWAS P-value thresholds (denoted PT), we calculated: 1) the significance of the case-control score difference was analyzed (P tab), 2) the proportion of variance explained (Nagelkerke’s R2, R2 tab), 3) the proportion of variance on the liability scale explained by RPS (h2I tab) with standard error in brackets, 4) area under the receiver operator characteristic curve (AUC tab), and 5) odds ratio for the 10th RPS decile group compared with lowest decile with confidence interval in brackets. Ncases tab denotes number of cases in each target set. (XLS 96 kb)
Supplementary Table 6
RPS analysis of MGS sample. Risk Profile Score (RPS) analyses was performed using the MGS dataset as target, using three distinct published results for SCZ GWAS, from the (1) ISC (2009) study of 2615 cases and 3338 controls11 (denoted ISC columns) (2) PGC1 (excluding MGS, denoted PGC1 columns) with 9320 cases and 10228 controls22, (3) current meta analysis (excluding MGS, denoted Current columns) with 32838 cases and 44357 controls. For various GWAS P value thresholds (denoted PT), we calculated 1) the significance of the case-control score difference was analyzed (P tab) 2) The proportion of variance explained (Nagelkerke’s R2, R2 tab) 3) The proportion of variance on the liability scale explained by RPS (h2I tab) with standard error in brackets 4) Area under the receiver operator characteristic curve (AUC tab) and 5) Odds ratio for 10th RPS decile group compared with lowest decile with confidence interval in brackets. Ncases tab denotes number of cases in each target set. (XLS 32 kb)
Rights and permissions
About this article
Cite this article
Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature 511, 421–427 (2014). https://doi.org/10.1038/nature13595
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nature13595