eTable 1. Genes Impacted by Both HSS150612 and HSS150613 ZNF804A siRNA Conditions (Hill, 2012)
eTable 2. Patient Analyses: ZNF804A Pathway Polygenic Score Regression Analysis for Each Neuropsychological Variable
eTable 3. Patient Analyses: ZNF804A Pathway Polygenic Score Regression Analysis for Each Neuropsychological Variable, Excluding ZNF804A SNPs
eAppendix 1. Members of the Schizophrenia Psychiatric Genome-wide Association Study (GWAS) Consortium
eAppendix 2. Members of The Wellcome Trust Case Control Consortium 2
Nicodemus KK, Hargreaves A, Morris D, Anney R, Gill M, Corvin A, Donohoe G, for the Schizophrenia Psychiatric Genome-wide Association Study (GWAS) Consortium and The Wellcome Trust Case Control Consortium 2. Variability in Working Memory Performance Explained by Epistasis vs Polygenic Scores in the ZNF804A Pathway. JAMA Psychiatry. 2014;71(7):778-785. doi:10.1001/jamapsychiatry.2014.528
We investigated the variation in neuropsychological function explained by risk alleles at the psychosis susceptibility gene ZNF804A and its interacting partners using single nucleotide polymorphisms (SNPs), polygenic scores, and epistatic analyses. Of particular importance was the relative contribution of the polygenic score vs epistasis in variation explained.
To (1) assess the association between SNPs in ZNF804A and the ZNF804A polygenic score with measures of cognition in cases with psychosis and (2) assess whether epistasis within the ZNF804A pathway could explain additional variation above and beyond that explained by the polygenic score.
Design, Setting, and Participants
Patients with psychosis (n = 424) were assessed in areas of cognitive ability impaired in schizophrenia including IQ, memory, attention, and social cognition. We used the Psychiatric GWAS Consortium 1 schizophrenia genome-wide association study to calculate a polygenic score based on identified risk variants within this genetic pathway. Cognitive measures significantly associated with the polygenic score were tested for an epistatic component using a training set (n = 170), which was used to develop linear regression models containing the polygenic score and 2-SNP interactions. The best-fitting models were tested for replication in 2 independent test sets of cases: (1) 170 individuals with schizophrenia or schizoaffective disorder and (2) 84 patients with broad psychosis (including bipolar disorder, major depressive disorder, and other psychosis).
Main Outcomes and Measures
Participants completed a neuropsychological assessment battery designed to target the cognitive deficits of schizophrenia including general cognitive function, episodic memory, working memory, attentional control, and social cognition.
Higher polygenic scores were associated with poorer performance among patients on IQ, memory, and social cognition, explaining 1% to 3% of variation on these scores (range, P = .01 to .03). Using a narrow psychosis training set and independent test sets of narrow phenotype psychosis (schizophrenia and schizoaffective disorder), broad psychosis, and control participants (n = 89), the addition of 2 interaction terms containing 2 SNPs each increased the R2 for spatial working memory strategy in the independent psychosis test sets from 1.2% using the polygenic score only to 4.8% (P = .11 and .001, respectively) but did not explain additional variation in control participants.
Conclusions and Relevance
These data support a role for the ZNF804A pathway in IQ, memory, and social cognition in cases. Furthermore, we showed that epistasis increases the variation explained above the contribution of the polygenic score.
Genome-wide association studies (GWASs) have been at the forefront of identifying candidate genes for schizophrenia. One of the genetic variants achieving genome-wide significance for psychosis was rs1344706 in the zinc finger binding protein 804A (ZNF804A).1 Several independent replication studies2,3 and a recent meta-analysis4 have supported an association between schizophrenia and the risk allele of this single nucleotide polymorphism (SNP). The gene ZNF804A, which is expressed in the brain, is believed to encode a protein with a C2H2 zinc finger domain. This suggests a role in the regulation of gene expression through DNA and/or RNA binding.5 The gene ZNF804A has been reported to show an association with brain activity and structure.6,7 A recent study by Hill et al8 assessed the effects of its knockdown on the cellular transcriptome, which linked ZNF804A to cell adhesion molecules, suggesting a role in neural migration, neurite outgrowth, and synapse formation, which are hypothesized to be aberrant in schizophrenia.
Several studies have linked ZNF804A to cognition, based on imaging studies,6,9 traditional neuropsychological measures,10- 12 and measures of social cognition.13,14 The results of these studies deviate from what would be expected in that while the rs1344706 risk allele appears to convey impairments in cognition in control individuals based on behavioral and imaging studies,11,15- 17 the literature points to preserved cognition in patients.10,12,18,19 Although some studies suggest that ZNF804A is associated with impaired social cognition in control participants, it is uncertain whether ZNF804A might also confer a disadvantage to patients because only 1 study14 has assessed patients. If ZNF804A confers risk for psychosis, why does the risk allele–carrying patient population show preserved cognitive function? The answer may lie in the impact of ZNF804A embedded within its functional pathway. The development of complex traits, such as psychosis, is likely to involve the contribution of a large number of independent and/or interacting genetic variants, mostly of modest effect.20,21 In case-control analyses, polygenic scores (a simple or weighted summation of the top sets of SNPs) have been shown to predict case status in related disorders and explain a significant percentage of variability.22 Limiting this polygenic risk score to include variants within genes shown to be altered by ZNF804A knockdown,8 we investigated whether more of the variance in patients’ neuropsychological function can be explained than is explained by single variants. We used the P values and odds ratios from the Schizophrenia Psychiatric GWAS Consortium 1 (PGC1) schizophrenia case-control analysis23 to rank SNPs for inclusion in the polygenic score. As variation in complex traits is thought to include both polygenic and epistatic components, we examined whether epistasis between SNPs within the ZNF804A pathway could explain variation above that explained by the polygenic score alone. We used half of the narrow psychosis (schizophrenia and schizoaffective disorder) set as a training set to test for pairwise epistasis among all SNPs included in the polygenic score, then assessed whether adding these interactions to the regression model containing the polygenic score increased the R2 among 3 independent test sets including (1) additional narrow psychosis cases, (2) broad psychosis (bipolar disorder, major depressive disorder, and psychosis not otherwise specified), and (3) healthy control participants.
Four hundred twenty-four cases and 89 healthy participants who completed a full neuropsychological assessment battery and for whom GWAS data were available were included. Cases were clinically stable patients with a DSM-IV diagnosis of schizophrenia (n = 282), schizoaffective disorder (n = 58), bipolar disorder (n = 61), major depressive disorder with psychotic features (n = 11), or psychosis not otherwise specified (n = 12) recruited from 5 sites across Ireland (Table 1). Inclusion criteria required that participants were clinically stable at neuropsychological assessment, aged 18 to 65 years, had no history of comorbid psychiatric disorder, no substance abuse in the preceding 6 months, no prior head injury with loss of consciousness, and no history of seizures. Diagnosis was confirmed by trained psychiatrists using the Structured Clinical Interview for DSM-IV-TR.24 Additional diagnostic details and clinical sample characteristics, including symptom severity (Scale for the Assessment of Negative Symptoms/Scale for the Assessment of Positive Symptoms)25,26 and medication dosage, are detailed elsewhere.10 Healthy control participants were recruited via online and poster advertising. They were aged 18 to 65 years, with no history of substance abuse in the preceding 6 months, no prior head injury with loss of consciousness, no history of seizures, and no personal history of psychosis or in their first-degree relatives. All assessments were conducted in accordance with the relevant ethics committees’ approval from each participating site. All participants had 4 grandparents born in Ireland and provided written informed consent.
Participants completed a neuropsychological assessment battery designed to target the cognitive deficits of schizophrenia including general cognitive function, episodic memory, working memory, attentional control, and social cognition. General cognitive functioning (IQ) was measured using selected subtests (Vocabulary, Similarities, Block Design, and Matrix Reasoning) from the Wechsler Adult Intelligence Scale,27 yielding a full-scale, verbal, and performance IQ. Episodic memory was assessed using the logical memory subtest from the Wechsler Memory Scale III.28 Working memory was assessed using the spatial working memory task (SWM) from the Cambridge Neuropsychological Test Automated Battery29 and letter-number sequencing from the Wechsler Memory Scale III.28 Attentional control was assessed using the Continuous Performance Task–Identical Pairs version,30 the Intradimensional-Extradimensional shift task (IDED) (Cambridge Neuropsychological Test Automated Battery),29 and the Sustained Attention to Response Task.31 Social cognition was assessed using the Reading the Mind in the Eyes Task32 and the Internal, Personal, and Situational Attributions Questionnaire (IPSAQ),33 which yields 2 bias scores; externalizing bias (EB), which indicates a propensity to attribute positive events to oneself rather than to other people, and a personalizing bias, which indicates a propensity to attribute negative events to other people rather than to situational factors.
Genotyping was conducted on DNA extracted from blood in patients and saliva in control participants. Single nucleotide polymorphism data were obtained from a recent GWAS using the Affymetrix SNP Array 6.0 platform, conducted as part of the Wellcome Trust Case Control Consortium 2, described in detail elsewhere.34
Polygenic scores22 for variants located within the ZNF804A pathway were calculated starting with all available SNPs within 20 Kb of genes in the ZNF804A pathway8 (eTable 1 in Supplement). Target alleles at these SNPs were identified as risk associated based on the PGC1 GWAS. Two polygenic scores for each individual were calculated. The count polygenic score was based on the simple sum of the number of risk-associated alleles they carried averaged across the total number of valid genotypes for that individual. The weighted polygenic score was based on the count polygenic score, with the exception that each SNP was weighted by the log of the odds ratio from the PGC1. We used 3 P value thresholds from the PGC1 case-control analysis as arbitrary thresholds as in a previous polygenic analysis22 (thresholds and number of SNPs: P < 1.0e-05, n = 10; P < .05, n = 218; and P < .50, n = 1525). As individual ZNF804A SNPs have been associated with cognition, we examined whether an association observed between the pathway-based polygenic score remained after removing all ZNF804A SNPs (n = 27). A subset of the Wellcome Trust Case Control Consortium 2 sample was included in the PGC1 schizophrenia case-control analysis, thus the samples are not entirely independent. However, the outcome of interest in the PGC1 was case status, whereas in the present study the outcome of interest was cognition in cases with psychosis. As suggested in the supplementary materials of the original polygenic score report,22 we used all genotyped SNPs in genes within the ZNF804A pathway without pruning for linkage disequilibrium.
Associations between ZNF804A polygenic score and the phenotypes of IQ, episodic memory, working memory, attention, and social cognition were tested in multiple regression analyses implemented in SPSS version 1735 or the R Statistical Computing Environment.36 In each case, scores for each neuropsychological phenotype were entered as dependent variables, controlling for age and sex as necessary.
To test whether additional variation could be explained by epistasis beyond that explained by the polygenic score, we developed a novel approach within the context of the regression models previously described. To reduce multiple testing, we restricted our search for epistasis to the models where the polygenic score accounted for a small, but significant, amount of variation in the neuropsychological phenotype in broad or narrow psychosis including SWM strategy, IPSAQ-EB, and performance IQ. First, we created equal-sized training and test sets from the narrow psychosis cases (n = 170). Second, to obtain a stable and consistent estimate of the P value from 2-SNP interactions, we took 100 bootstrap samples with replacement from the narrow psychosis training sample and performed linear regression analysis for all possible pairs of SNPs falling under the threshold at hand (SWM threshold P ≤ .05, SNPs: n = 218; IPSAQ-EB and performance IQ threshold P < 1.0e−05, SNPs: n = 10). For performance IQ and IPSAQ-EB, the 10 SNPs were all on chromosome 10 and in tight linkage disequilibrium (range, r2 = 0.6-1.0), which led to collinearity in the interactions and were not tested further. For SWM, the linear regression model in the training data contained the unweighted polygenic score plus an epistatic term, which was the product of the risk-associated alleles at 2 SNPs. The unweighted polygenic score was used so the alleles comprising the score were on the same scale of measurement as those used in the interactions; however, the use of the weighted score did not change the results (discussed below). The average P value from the 100 replicates was used to determine which interactions would be evaluated using the 3 independent test sets, using an uncorrected P < .05 threshold estimated from the training set. To account for linkage disequilibrium, a further condition was that only the interaction with the smallest P value containing a particular SNP would be tested for replication in the independent test sets and all other interactions containing that SNP would not be considered for replication due to collinearity. This led to 3 significant 2-SNP interaction terms, independent of one another, being brought forward for replication in the 3 test sets. R2 values on the independent test sets were calculated as the square of the correlation between the fitted values for the test set based on the model estimated from the training data and the observed values from the test set.
Demographic and clinical characteristics for patients and healthy participants appear in Table 1. The characteristics of the broad psychosis group were compared with the narrow psychosis group and with the control group using t tests. No significant differences were observed between the narrow and broad psychosis groups for age, sex, age at onset, full-scale IQ, medication dosage as measured by chlorpromazine equivalents, or positive and negative symptoms, with the exception of the mania factor where the broad psychosis group scored significantly higher. The patient group contained significantly more men than the healthy group, was significantly older at the time of assessment, and had a significantly lower full-scale IQ.
ZNF804A weighted polygenic scores were associated (uncorrected P < .05) with measurements of both general and social cognition, including SWM, IPSAQ-EB, and performance IQ, although the effect size was moderate (1.2%-3%) (Table 2; results for all phenotypes are in eTable 2 in Supplement). Larger ZNF804A polygenic scores were predictive of poorer performance on performance IQ in the broad psychosis group at the threshold of P = 1.0e−05, but did not predict performance IQ in the narrow psychosis group (Table 2). The IPSAQ-EB demonstrated significant association with the ZNF804A polygenic score among broad and narrow psychosis at a threshold of P = 1.0e−05. Patients with psychosis who had greater polygenic scores demonstrated a decreased IPSAQ-EB score, suggesting that these patients were less likely to show what is known as a self-serving bias: the adaptive tendency to attribute causality for negative events to external factors and positive events to oneself. Among both narrow and broad psychosis groups, a higher ZNF804A polygenic score led to significantly poorer performance on SWM strategy at a threshold of P = .05.
To determine whether the polygenic score results were powered by ZNF804A SNPs, we excluded all variants in ZNF804A and reassessed association with the polygenic score; the results were virtually unchanged (Table 2), with the exception of SWM strategy, which showed a slightly reduced R2 and larger P value in the narrow psychosis set (results for all phenotypes are in eTable 3 in Supplement).
As just described, we tested for additional variation in SWM explained by 2-SNP epistasis in the narrow psychosis training set. Because of linkage disequilibrium, 112 average P values across 100 bootstrap samples of the training set were less than .05; however, of these, only 3 were completely independent of the other sets and thus only these 3 interaction terms were brought forward for replication in the 3 independent test sets. Beginning with SWM strategy, adding the most significant interaction (rs17186340:rs140512) to the model containing the polygenic count score led to an increase of 1.4% in the narrow psychosis test set (P = .05), 3.7% in the broad psychosis test set (P = .04), and a combined case test set increase of 2.3% (P = .006) (Table 3). Although adding this interaction increased the variation explained in 2 independent sets of cases, it did not increase variation explained in control participants (R2 = 0.0069, P = .45). Adding a second interaction term to the model (rs2295984:rs34138673) further increased the variation explained in cases: in the narrow psychosis test set, the R2 increased by 1.3% (P = .02) to a total of 4.0%; in the broad psychosis test set, the R2 increased further by 1.2% (P = .04) for a total variation explained of 6.2%; and in the combined set of cases, the R2 increase was 1.3% (P = .001) for a total of 4.8%. In control participants, adding the second interaction did not increase the variation explained (R2 = 0.001, P = .77). The third epistatic term did not increase the R2 values in any independent test set and thus was not considered further.
To test whether the increase in R2 was attributable to strongly associated SNPs contained in the interactions themselves, interaction SNPs were tested individually for association with SWM strategy in a model containing the polygenic count score. Two SNPs were associated with SWM strategy at an uncorrected P < .05 and both were in the second interaction term (rs17186340:rs34138673). The improvement in R2 values from the narrow psychosis training set model containing the polygenic count score plus rs17186340 on the narrow psychosis, broad psychosis, and total test cases set were 5.89e−06, 9.1e−04, and 8.0e−05, respectively. Results for the same analysis using rs34138673 were 0.0056, 0.013, and 0.0079, respectively, indicating that the interaction term explained more variation than either single SNP.
To see whether the type of polygenic score influenced the amount of variation explained, the weighted polygenic score was substituted for the count polygenic score and R2 values were calculated for the case test sets. In the narrow psychosis test set, the use of the weighted polygenic score in the 3 models (polygenic only, 1 and 2 interactions) increased the R2 value by 9.1e−04 − 0.001 vs models containing the count polygenic score, whereas the use of the weighted polygenic score in broad psychosis reduced the R2 value by 4.7e−04 − 0.0061 vs the count polygenic score, suggesting the choice of polygenic score was trivial in this instance. When strongly associated SNPs with very large or very small effect sizes (odds ratios) are present, we would expect to see differences between using the count or weighted score.
We used polygenic scores to investigate whether ZNF804A pathway schizophrenia risk–associated alleles were associated with neuropsychological function among 424 patients with psychosis. Higher ZNF804A polygenic scores were significantly associated with poorer performance in IQ, working memory, and biased social cognition and explained 1% to 3% of the variation in these measures, consistent with estimates previously reported in general intelligence in control individuals.37 Removal of the SNPs within ZNF804A reduced the R2 values only slightly, suggesting the combined contribution of genes within the pathway was driving the association. Furthermore, we showed that considering epistasis along with the polygenic score resulted in more than 3 times the amount of variation explained in 2 independent test sets of cases (range, total R2 = 4.0%-6.2%). Because the SNPs participating in interactions were not in ZNF804A, we provide further evidence that this gene was not the key contributor to our pathway-based results. Thus, our findings are not inconsistent with previous studies showing preserved cognitive function in cases is associated with ZNF804A. Although the polygenic score explained a similar amount of variation in control participants, improvements due to epistasis were specific to psychosis cases.
Although previous studies have shown that the risk allele of ZNF804A rs1344706 shows differential effects in cases and control participants,10- 19 we showed that, at the pathway level, the effect of the combination of schizophrenia risk alleles leads to poorer performance in patients with psychosis on measures of intelligence, working memory, and social cognition. The 4 SNPs participating in epistasis that increased the R2 values in our test sets were near STAC (rs17186340), MAPK8IP2 (rs140512), and flanking either side of FAM46A (rs2295984 and rs34138673). In mice, Stac is expressed in the brain, neurons, and postsynaptic densities; within the brain, the expression is highest in the hippocampus and cerebellum.38,39 A Stac knockout mouse model showed reduced social interaction, impaired learning, and deficits in exploration of novel environments.39MAPK8IP2 is located within the Chr22q13.3 deletion region associated with autism spectrum disorders and Phelan-McDermid syndrome, which is characterized by developmental delay. Also known as JIP2, it is a scaffold protein that is necessary for N-methyl-D-aspartic acid receptor function and modulates signal transduction.40FAM46A is expressed in adult human brain and shows higher expression in human fetal brain.41 The SNPs participating in the second interaction term, rs2295984 and rs34138673, are located on either side of FAM46A, approximately 19.5 kbp apart, possibly indicating a promoter and/or enhancer role.
How can we reconcile previous research that has shown the risk allele at ZNF804A rs1344706 is associated with less impaired cognition in patients with psychosis with the results of the present study, which showed that the polygenic score from the ZNF804A pathway was associated with poorer performance IQ, working memory, and social cognition? The P values from the PGC1 were used to select SNPs for inclusion in the polygenic score, and the smallest ZNF804AP value was .002 for rs1344706. Therefore, the set with the most stringent threshold did not include any ZNF804A SNPs, and this set was negatively associated with IQ and social cognition. For working memory, the removal of the ZNF804A SNPs at a threshold of P = .05 would have included 12 of the 27 SNPs with P values ranging between .002 and .05. We have shown that the removal of these SNPs did not lead to significant differences in the magnitude of association between the polygenic score and working memory. The use of the weighted polygenic score would have ensured a weak contribution of these SNPs because they were not strongly associated with schizophrenia in the PGC1 and they comprised only 5.5% of the total number of SNPs at that threshold. Interestingly, the PGC1 P values for the 4 SNPs participating in epistasis ranged between .007 (rs2295984) and .04 (rs17186340), showing that although they are marginally associated with schizophrenia, they would not have been considered for follow-up. As was the case with the previous use of the polygenic score,22 we showed that the polygenic score and epistatic models based on a narrow psychosis training sample were able to significantly account for variation in working memory in 2 independent psychosis samples, but they were not able to predict variation in control participants. Although the control sample size was modest (n = 89), the broad psychosis sample was of a similar size (n = 84) so a lack of statistical power cannot fully explain the inability to account for additional variation in control participants.
Interestingly, the 10 SNPs in linkage disequilibrium that comprised the most strongly associated polygenic score for IPSAQ-EB and performance IQ included 5 SNPs in CNNM2 and 5 SNPs either 3′ or 5′ of the gene. This gene is strongly associated with schizophrenia23,42,43 and we and other groups have shown variation within the gene is associated with gray matter volume in patients with schizophrenia44,45 and with attributional style.45 Our results support a polygenic contribution of variation within and around this gene weakly contributing to attributional style and performance IQ.
We have introduced a novel method to evaluate the combined effect of the polygenic score and epistasis, which is both simple and computationally tractable. The addition of the epistatic terms also increased the interpretability of the model, as it is difficult to determine which genes were contributing signal to the polygenic score. We showed that the results estimated on the training sample were generalizable to 2 independent test sets of patients—one with narrow psychosis and the other with nonschizophrenia psychosis—similar to previous studies’ use of the polygenic score.22 In both instances, the variation explained by our epistatic terms was much larger than that explained by the polygenic score itself: the polygenic score explained 1.2% to 1.3% of variation, whereas increases in R2 using our novel approach in the test sets were between 2.7% and 4.9%. Epistasis is thought to be a key element in complex phenotypes20,21 and has been shown to influence the risk for schizophrenia and inefficient dorsolateral prefrontal cortex processing during a working memory task in healthy control individuals.20,21,46,47 The potential limitations of our study included modest sample sizes and the fact that we may not have captured all variation (especially rare variation) in the genes in the ZNF804A pathway because of our reliance on SNPs from a GWAS.
To our knowledge, this study is the first to investigate the role of the ZNF804A pathway in the cognitive decline commonly observed among patients with psychosis. We have identified 3 new candidate genes for working memory in the ZNF804A pathway: STAC, MAPK8IP2, and FAM46A. Perhaps more critically, we introduced an improvement in the use of polygenic scores by adding an epistatic component that explained additional variation in working memory that was specific to cases with psychosis.
Corresponding Author: Kristin K. Nicodemus, PhD, MPH, Neuropsychiatric Genetics Group, Department of Psychiatry, Trinity College Dublin, St James Hospital, Dublin 8, Ireland (email@example.com).
Submitted for Publication: November 21, 2013; final revision received January 9, 2014; accepted February 18, 2014.
Published Online: May 14, 2014. doi:10.1001/jamapsychiatry.2014.528.
Author Contributions: Dr Nicodemus had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Dr Nicodemus and Ms Hargreaves contributed equally to this work.
Study concept and design: Nicodemus, Hargreaves, Anney, Gill, Corvin, Donohoe.
Acquisition, analysis, or interpretation of data: All authors.
Drafting of the manuscript: Nicodemus, Hargreaves, Donohoe.
Critical revision of the manuscript for important intellectual content: Nicodemus, Morris, Anney, Gill, Corvin, Donohoe.
Statistical analysis: Nicodemus, Hargreaves, Anney, Donohoe.
Obtained funding: Nicodemus, Gill, Corvin, Donohoe.
Administrative, technical, or material support: Morris, Anney, Corvin, Donohoe.
Study supervision: Corvin, Donohoe.
Conflict of Interest Disclosures: None reported.
Funding/Support: Recruitment and genotyping was supported by Science Foundation Ireland (SFI) (grant 08/IN.1/B1916) and the Wellcome Trust Case Control Consortium 2 project (grants 085475/B/08/Z and 085475/Z/08/Z) and the Wellcome Trust (grants 072894/Z/03/Z, 090532/Z/09/Z, and 075491/Z/04/B), respectively. This publication has emanated from research conducted with the financial support of SFI and the Marie-Curie Action COFUND under grant 11/SIRG/B2183 to Dr Nicodemus. Dr Donohoe’s work is generously supported by grant funding from the Health Research Board (HRA_POR/2012/54) and SFI (12.IP.1359).
Role of the Sponsor: The funding bodies had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Group Information: The Schizophrenia Psychiatric Genome-wide Association Study (GWAS) Consortium investigators include Stephan Ripke, Alan R. Sanders, Kenneth S. Kendler, Douglas F. Levinson, Pamela Sklar, Peter A. Holmans, Dan-Yu Lin, Jubao Duan, Roel A. Ophoff, Ole A. Andreassen, Edward Scolnick, Sven Cichon, David St. Clair, Aiden Corvin, Hugh Gurling, Thomas Werge, Dan Rujescu, Douglas H. R. Blackwood, Carlos N. Pato, Anil K. Malhotra, Shaun Purcell, Frank Dudbridge, Benjamin M. Neale, Lizzy Rossin, Peter M. Visscher, Danielle Posthuma, Douglas M. Ruderfer, Ayman Fanous, Hreinn Stefansson, Stacy Steinberg, Bryan J. Mowry, Vera Golimbet, Marc De Hert, Erik G. Jönsson, István Bitter, Olli P. H. Pietiläinen, David A. Collier, Sarah Tosato, Ingrid Agartz, Margot Albus, Madeline Alexander, Richard L. Amdur, Farooq Amin, Nicholas Bass, Sarah E. Bergen, Donald W. Black, Anders D. Børglum, Matthew A. Brown, Richard Bruggeman, Nancy G. Buccola, William F. Byerley, Wiepke Cahn, Rita M. Cantor, Vaughan J. Carr, Stanley V. Catts, Khalid Choudhury, C. Robert Cloninger, Paul Cormican, Nicholas Craddock, Patrick A. Danoy, Susmita Datta, Lieuwe de Haan, Ditte Demontis, Dimitris Dikeos, Srdjan Djurovic, Peter Donnelly, Gary Donohoe, Linh Duong, Sarah Dwyer, Anders Fink-Jensen, Robert Freedman, Nelson B. Freimer, Marion Friedl, Lyudmila Georgieva, Ina Giegling, Michael Gill, Birte Glenthøj, Stephanie Godard, Marian Hamshere, Mark Hansen, Thomas Hansen, Annette M. Hartmann, Frans A. Henskens, David M. Hougaard, Christina M. Hultman, Andrés Ingason, Assen V. Jablensky, Klaus D. Jakobsen, Maurice Jay, Gesche Jürgens, René S. Kahn, Matthew C. Keller, Gunter Kenis, Elaine Kenny, Yunjung Kim, George K. Kirov, Heike Konnerth, Bettina Konte, Lydia Krabbendam, Robert Krasucki, Virginia K. Lasseter, Claudine Laurent, Jacob Lawrence, Todd Lencz, F. Bernard Lerer, Kung-Yee Liang, Paul Lichtenstein, Jeffrey A. Lieberman, Don H. Linszen, Jouko Lönnqvist, Carmel M. Loughland, Alan W. Maclean, Brion S. Maher, Wolfgang Maier, Jacques Mallet, Pat Malloy, Manuel Mattheisen, Morten Mattingsdal, Kevin A. McGhee, John J. McGrath, Andrew McIntosh, Duncan E. McLean, Andrew McQuillin, Ingrid Melle, Patricia T. Michie, Vihra Milanova, Derek W. Morris, Ole Mors, Preben B. Mortensen, Valentina Moskvina, Pierandrea Muglia, Inez Myin-Germeys, Deborah A. Nertney, Gerald Nestadt, Jimmi Nielsen, Ivan Nikolov, Merete Nordentoft, Nadine Norton, Markus M. Nöthen, Colm T. O’Dushlaine, Ann Olincy, Line Olsen, F. Anthony O’Neill, Torben F. Ørntoft, Michael J. Owen, Christos Pantelis, George Papadimitriou, Michele T. Pato, Leena Peltonen, Hannes Petursson, Ben Pickard, Jonathan Pimm, Ann E. Pulver, Vinay Puri, Digby Quested, Emma M. Quinn, Henrik B. Rasmussen, János M. Réthelyi, Robert Ribble, Marcella Rietschel, Brien P. Riley, Mirella Ruggeri, Ulrich Schall, Thomas G. Schulze, Sibylle G. Schwab, Rodney J. Scott, Jianxin Shi, Engilbert Sigurdsson, Jeremy M. Silverman, Chris C. A. Spencer, Kari Stefansson, Amy Strange, Eric Strengman, T. Scott Stroup, Jaana Suvisaari, Lars Terenius, Srinivasa Thirumalai, Johan H. Thygesen, Sally Timm, Draga Toncheva, Edwin van den Oord, Jim van Os, Ruud van Winkel, Jan Veldink, Dermot Walsh, August G. Wang, Durk Wiersma, Dieter B. Wildenauer, Hywel J. Williams, Nigel M. Williams, Brandon Wormley, Stan Zammit, Patrick F. Sullivan, Michael C. O’Donovan, Mark J. Daly, and Pablo V Gejman.
The Wellcome Trust Case Control Consortium 2 investigators include Peter Donnelly, Ines Barroso, Jenefer M. Blackwell, Elvira Bramon, Matthew A. Brown, Juan P. Casas, Aiden Corvin, Panos Deloukas, Audrey Duncanson, Janusz Jankowski, Hugh S. Markus, Christopher G. Mathew, Colin N. A. Palmer, Robert Plomin, Anna Rautanen, Stephen J. Sawcer, Richard C. Trembath, Ananth C. Viswanathan, Nicholas W. Wood, Chris C. A. Spencer, Gavin Band, Céline Bellenguez, Colin Freeman, Garrett Hellenthal, Eleni Giannoulatou, Matti Pirinen, Richard Pearson, Amy Strange, Zhan Su, Damjan Vukcevic, Cordelia Langford, Sarah E. Hunt, Sarah Edkins, Rhian Gwilliam, Hannah Blackburn, Suzannah J. Bumpstead, Serge Dronov, Matthew Gillman, Emma Gray, Naomi Hammond, Alagurevathi Jayakumar, Owen T. McCann, Jennifer Liddle, Simon C. Potter, Radhi Ravindrarajah, Michelle Ricketts, Matthew Waller, Paul Weston, Sara Widaa, and Pamela Whittaker.
The affiliations for the group members can be found in the eAppendices in the Supplement.
Additional Contributions: We sincerely thank all patients who contributed to this study and all staff who facilitated their involvement.