PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Mol Psychiatry. Author manuscript; available in PMC 2013 May 1.
Published in final edited form as:
PMCID: PMC3235542
NIHMSID: NIHMS315559

Statistical Epistasis and Progressive Brain Change in Schizophrenia: An Approach for Examining the Relationships Between Multiple Genes

Nancy C. Andreasen, MD, PhD, Marsha A. Wilcox, EdD, ScD,*,** Beng-Choon Ho, MD,* Eric Epping, MD, PhD,** Steven Ziebell, BS,* Eugene Zeien, BS,* Brett Weiss, BS,* and Thomas Wassink, MD**

Abstract

Although schizophrenia is generally considered to occur as a consequence of multiple genes that interact with one another, very few methods have been developed to model epistasis. Phenotype definition has also been a major challenge for research on the genetics of schizophrenia. In this report we use novel statistical techniques to address the high dimensionality of genomic data, and we apply a refinement in phenotype definition by basing it on the occurrence of brain changes during the early course of the illness, as measured by repeated MR scans (i.e., an “intermediate phenotype.” The method combines a machine learning algorithm, the ensemble method using stochastic gradient boosting, with traditional general linear model statistics. We began with fourteen genes that are relevant to schizophrenia based on association studies or their role in neurodevelopment and then used statistical techniques to reduce them to five genes and 17 SNPs that had a significant statistical interaction: 5 for PDE4B, 4 for RELN, 4 for ERBB4, 3 for DISC1, and one for NRG1. Five of the SNPs involved in these interactions replicate previous research, in that these five SNPs have previously been identified as schizophrenia vulnerability markers or implicate cognitive processes relevant to schizophrenia. This ability to replicate previous work suggests that our method has potential for detecting a meaningful epistatic relationships among the genes that influence brain abnormalities in schizophrenia.

Introduction

Schizophrenia is conceptualized as a genetically complex disorder that is probably the result of multiple alleles of small effect in many different genes, at least some of which are likely to interact with one another.1 Despite the widely-recognized plausibility of this conceptualization, progress in delineating the genetic and neural mechanisms of schizophrenia has been limited by the difficulty of developing adequate models that can account for genetic complexity, particularly epistasis. Epistasis has several connotations: one in which a gene suppresses the effect of another; and a number of genes interacting together to influence the expression of the phenotype. Here, we use epistasis in the latter sense.

Phenotype definition has also been a major challenge for research on the genetics of schizophrenia. Schizophrenia is characterized by a variable mixture of psychotic, disorganized, and negative symptoms; it also has a broad variety of clinical outcomes, ranging from prolonged periods of remission to a disabling treatment-refractory course that is sometimes accompanied by cognitive decline.2 Most current genetic studies do not take the complexity of the phenotype into account. Genome-wide association studies (GWAS) have been disappointing because of poor replicability, a problem likely to be due in part to weak or imprecise phenotype definition.315 GWAS only require that subjects meet specified diagnostic criteria, such as those put forth in DSM IV or ICD 10. This limited approach to phenotype definition runs the risk of grouping clinically (and probably genetically) diverse patients under a single rubric, thereby reducing the likelihood that replicable and accurate associations will be found. In this report we describe several approaches that we have implemented to address the problems of examining epistasis and improving phenotype definition.

An inherent problem in the effort to link improved phenotype definition with complex genetic concepts such as epistasis is the difficulty inherent in detecting meaningful signals within the noisy vastness of the human genome. Various approaches have been developed to reduce the high dimensionality of genomic data, such as the use of gene ontology or pathway analysis.16 Another approach for addressing high dimensionality is the use of statistical techniques such as machine learning algorithms (MLA). The utility of these approaches has been demonstrated recently in studies that have used techniques such as recursive partitioning to model epistatic relationships.17 We have developed a strategy for using similar algorithms to identify potentially meaningful relationships between schizophrenia vulnerability genes. Our approach employs a MLA, stochastic gradient boosting using regression trees,18 and performs a series of steps to winnow through a large genomics data set and to identify genes and SNPs that may interact with one another to produce brain tissue loss in schizophrenia.

In this report we also address the problem of phenotype definition by basing it on the organ that produces the clinical picture of the illness: the brain. Numerous well-replicated studies have documented that patients with schizophrenia have measurable brain tissue reduction when compared with healthy volunteers.1933 The tissue loss affects multiple brain regions, particularly overall cerebral volume and frontal and temporal lobe gray matter (GM) and white matter (WM); tissue loss in these regions is accompanied by a corresponding increase in cerebrospinal fluid (CSF). These brain volume reductions are present at the onset of the illness, but they also continue to progress over time, with the greatest severity occurring during the first few years after onset.33 These tissue changes are considered to be due to an aberration in a neurodevelopmental process for several reasons: the typical onset of illness during adolescence, when the brain is undergoing major developmental changes; and the absence of an neuropathologic signature suggesting a degenerative process such as neuronal loss or gliosis.

Because we have conducted a prospective longitudinal study of patients with schizophrenia ascertained at the time of onset ("first episode patients") and followed them with repeated Magnetic Resonance (MR) scans for a time period up to eighteen years, we have been able to develop a refinement in phenotype definition by basing it on the occurrence of brain changes during the early course of the illness, as measured by repeated MR scans.33

Therefore, in this report we describe a data analytic strategy for evaluating epistatic relationships between multiple genes and their alleles with a biologically-based phenotype of schizophrenia. This strategy is summarized in Figure 1 and consists of seven steps:

  1. Measure brain change over time in multiple brain regions, using structural Magnetic Resonance (sMR) volumetric measures in a sample of first episode schizophrenia patients
  2. Identify candidate genes for analysis, based on their relevance to neurodevelopmental processes that regulate brain change over time or their association with schizophrenia, and genotype patients for these genes (See Table 1)
    Table 1
    Fourteen Genes Evaluated for Epistasis
  3. Use stochastic gradient boosting (SGB) to reduce the number of genes to those that have the strongest relationship with multiple measures of brain change, as indicated by variable importance measures for genes and their associated SNPs (See Table 2 and Figure 2)
    Figure 2
    Variable Importance of Genes for Predicting Brain Volume Change Based on Initial Stochastic Gradient Boosting
    Table 2
    Variable Importance* of Genes for Predicting Brain Volume Change Based on Initial Stochastic Gradient Boosting
  4. Repeat SGB for the genes that survive the initial SGB to identify the top 20 SNPs with the highest VI; reduce data to those genes/SNPs that may have a significant epistatic relationship by determining which ones share an influence on the same specific brain regions
  5. Reduce data for those genes/SNPs by using a GLM analysis to find the genes/SNPs that have a significant gene/SNP by gene/SNP interaction
  6. Further reduce data by requiring that gene/SNP interactions affect at least two brain regions in a biologically plausible pattern (See Table 3)
    Table 3
    Epistatic Relationships Between Genes in Specific Brain Regions
  7. Examine relationships between specific genotypes and severity of brain tissue change over time in order to identify potential epistatic relationships that affect disease progression in the brain

Figure 1
Steps Involved in Stochastic Gradient Boosting

Materials and methods

Study sample

Subjects for this study were recruited through the Iowa Longitudinal Study (ILS). The ILS was initiated nineteen years ago and includes a total cohort of 542 first episode schizophrenia patients. These patients were recruited from consecutive admissions to the University of Iowa Psychiatry inpatient service at a rate of 25–30/year; recruitment ended in 2007. They have been followed at six month intervals after initial intake, with assessment of clinical symptoms, psychosocial function, and treatment received. Intensive assessments (sMR and cognitive testing) are done at intake and at two, five, nine, twelve, fifteen, and eighteen years.33 In this report we focus on those subjects for whom we have both genomic data and adequate sMR data to provide a relatively definitive determination of whether progressive brain change occurs over a time interval that is up to 12 years after intake. These comprise a total of 144 patients for whom we have a minimum of 2 scans and a maximum of six. (See Table 1 in Supplementary Materials for demographics and clinical characteristics.)

sMR data acquisition and analysis

Acquisition and analysis of sMR data have been previously described.33,3944 Details are available in supplementary materials.

Genotyping

To genotype study participants, DNA was prepared by high-salt extraction from whole blood and assayed using Illumina Infininum II array BeadChips which were designed, manufactured and completed by Illumina (San Diego, CA). The locus success rate and genotype call rate on the 26,122-SNP BeadChips were 97.2% and 99.8% respectively. We used a customized Illumina chip provided by Ortho-McNeill Janssen that contains SNPs for 1204 genes; the SNPs are primarily tag SNPs. The genes were selected for inclusion on the chip based on their association with major mental illnesses, their role in neurotransmission, their role in regulating major metabolic or developmental pathways, and their role in drug metabolism. For this study we selected fourteen genes for examination of statistical epistasis in relation to brain change measures. Because we were interested in the specific examination of the neurodevelopmental mechanisms influencing brain change in schizophrenia, we selected the genes based on the following criteria: established association with schizophrenia; influence on neurodevelopmental pathways; neurotransmitters implicated in schizophrenia. The rationale for the selection of these fourteen genes (810 SNPs) is summarized in Table 1.

Statistical analysis

To create a measure of percent change for use in our analyses, the MRI volumes were divided by the intracranial volume to correct for variability in head size. We have previously found that the greatest brain changes occur early in the course of the illness.33 Therefore, our measures are a percent change score based on the earliest available scan and its subsequent one. The value for the first scan was subtracted from the value for the second scan; that value was then divided by the number of years between the scans; the result was then divided by the value for the first scan, and finally that number was multiplied by 100 to provide an index of percent change.

The goal of our statistical approach was to identify a logical and empirically-driven method for detecting the signal (potential epistatic relationships between genes that may affect brain volume changes) within our vast array of genomic data. (See Figure 1)

Our first step was designed to determine which of the fourteen candidate genes had the strongest relationship with brain tissue change over time. We used regression trees to identify the relative importance of genes (SNPs) for predicting the amount of change in the size of each of the brain regions. This analysis employed an ensemble approach with stochastic gradient boosting (SGB) (TreeNet).18 Variable importance (VI) scores, generated by the SGB analysis, were used to identify the SNPs and genes that were retained during the initial steps of variable reduction. Variable importance scores show the relative contribution of each of the variables (SNPs) to predicting the outcome (a measure of brain volume change).

VI scores are relative; the largest value is assigned a value of 100 and the remaining variables are scaled accordingly. The variable importance measure for each gene, as a predictor of each brain volume measure, is calculated using the variable importance of the top 20 SNPs for predicting each brain volume measure. It is a percent of variable importance, calculated by summing the VI scores for the top 20 SNPs and dividing each SNP VI score by the total of the 20. A genewise VI score is then calculated by adding up the VIs for each SNP to get a gene summary.

The steps involved are summarized as follows.

  1. Begin with the entire data set and sample rows and columns of the data matrix (SNPs and brain measures for each individual); randomly sample a subset of individuals/SNPs
  2. Build a small tree (maximum of 8 nodes) and compute variable importance for SNPs
  3. Randomly resample individuals and SNPS and repeat, building an ensemble of 1000 trees
  4. Aggregate VI across the ensemble

From SNPs identified in the second SGB, we used general linear models (GLM) to further assess SNP*SNP interaction effects on brain volume changes over time. Analyses were performed using SAS (version 9.2; SAS, Cary, North Carolina). Percent brain volume change measure was entered as the dependent variable. Genotype of SNP pairs and SNP*SNP interaction terms were entered as independent measures in the statistical models. For all analyses, intracranial volume at initial MR scan, gender, imaging protocol (MR5 versus MR6), age at initial MR scan, and amount of neuroleptic exposure measured using dose-years were included as covariates. Because our technique is exploratory, we did not correct for multiple comparisons; we report only findings with P<.02.

Results

Table 2 and Figure 2 show the VI of the fourteen genes for predicting severity of brain volume changes, broken down into fourteen specific brain regions. Three genes stand out by virtue of having a high variable importance for all the brain regions that we examined: DISC1, ERBB4, and RELN. Two others, NRG1 and PDE4B, had high variable importance for twelve regions. The remaining genes had lesser (e.g., BDNF, COMT, and GRM3) or very little importance (e.g., AKT1, DAOA, GDNF, DTNBP1, KCNH2, and RGS4). We therefore retained the top five genes for further examination.

We then repeated SGB using these five genes and their associated SNPs (N=735). We again identified the top 20 SNPs for each brain region based on VI. The analysis further reduced the “candidate SNPs” to a total of 267 (Table 2, Supplementary Materials). A prominent role continued to be played by the SNPs belonging to DISC1 (58 discrete SNPs), ERBB4 (87 discrete SNPs), and RELN (61 discrete SNPs). NRG1 and PDE4B contribute fewer (23 and 38 respectively). We then reduced these data by identifying groups of SNPs that may have an epistatic relationship because they share an association with a given specific brain region. These potential interactions were then tested for statistical significance by being placed in a general linear model, using relevant covariates, and determining whether there was a significant SNP*SNP interaction. Eighteen SNP pairs were found to pass this filter.

The final filter in the search consisted of examining these pairs in relation to the brain regions with which they were associated. We limited this final group to those SNP pairs that had a significant relationship with more than one brain region and that was biologically plausible (e.g., reduced cerebral volume accompanied by increased CSF, as an indicator of a generalized tissue loss; reduced WM in two regions as indicator of tissue specificity). Eleven interactions were selected for analysis using these criteria. They are shown in Table 3. They include all five genes retained up to this point and a total of 17 SNPs: 5 for PDE4B, 4 for RELN, 4 for ERBB4, 3 for DISC, and one for NRG1.

Five of these SNPs have been previously identified as increasing schizophrenia risk or are in linkage disequilibrium with known schizophrenia risk alleles or alleles affecting a cognitive function relevant to schizophrenia.4552 Many of them involve RELN. Two of the SNPs involve an interaction between RELN and PDE4B alleles.

One is an interaction between RELN rs2229860 and PDE4B rs54402. This interaction is associated with a decrease in cerebral tissue and an increase in surface CSF. Rs2229860 is an identified schizophrenia risk allele on the RELN gene; it is a rare A/G mutation that results in a proline to arginine change at position 1703.46 A second interaction is between RELN rs580884 and DISC1 rs3738401. This interaction is associated with decreases in cerebral tissue and an increase in VBR. Rs3738401 is an A/G polymorphism that is a missense mutation; the presence of the minor A allele leads to an arginine to glycine change at position 844. The significant epistatic relationship appears to derive from an AA/TG combination, which occurred in 8 subjects. When these two genotypes co-occur, they have the largest amount of cerebral tissue loss (−.769%) and the largest amount of VBR increase (14.65%). Both are statistically significant (p<.01). A third interaction is between RELN rs499953 and PDE4B rs11576970. This interaction is associated with decreases in frontal and cerebral white matter. The PDE4B SNP is located in or near block 5, 2.8 kilobases away from rs7412571, with which it is in linkage disequilibrium; rs7412571 is significantly associated with schizophrenia.53 The epistatic relationship derives from a TC/GG combination, resulting in a significant decrease in cerebral (−6.10%, p<.01) and frontal (−6.62%, p<.03) WM.

A fourth interaction is between DISC1 rs11578905 and PDE4B rs2455012. This interaction is associated with a decrease in surface and parietal GM. The PDE4B SNP is in LD with two identified alleles, rs1354064 and rs2503222, both of which increase risk for schizophrenia.48,53 A final interaction involving a previously identified allele is between NRG1 and RELN. This interaction is associated with an increase in the volume of the caudate and putamen. The RELN SNP, rs2237628, is in LD with rs2711870, which has been implicated with problems in shifting response set on a card sorting task, an impairment in executive function that has repeatedly been shown to be abnormal in schizophrenia.54 The interaction of the RELN and the NRG1 heterozygotes leads to an increase of 2.22% in the caudate and 3.90% in the putamen.

Discussion

We have described a method for identifying potential epistatic relationships between genes that may contribute to schizophrenia and affect disease progression in the brain in patients suffering from schizophrenia. The method combines a machine learning algorithm, the ensemble method using stochastic gradient boosting, with traditional general linear model statistics. The method was used to identify genes/SNPs that were interacting with one another and predicting a continuous outcome measure that is a biologically meaningful phenotype (“intermediate phenotype”) for schizophrenia: changes in brain structure occurring after the onset of the illness. The method began with fourteen genes and 810 SNPs and reduced them to five genes and seventeen SNPs, identifying eleven interactions. Five of the SNPs involved in these interactions replicate previous research, in that these five SNPs have previously been identified as schizophrenia vulnerability markers or implicate cognitive processes relevant to schizophrenia. This ability to replicate previous work suggests that our method has potential for detecting a meaningful signal within the human genome.

RELN was the gene that was found to have the highest number of interactions. It had a total of six, involving PDE4B, DISC1, and NRG1. RELN’s relationship with schizophrenia is well-supported by linkage or association studies,45,46 and it has also been shown to be associated with treatment resistance.57 RELN expression is markedly reduced (by 50%) in mRNA and protein in post mortem brain tissue from schizophrenia patients.58 Epigenetic abnormalities have also been found; hypermethylation in the promotor region may be an explanation for the reduced expression of RELN in the brain.55,56 The known functions of RELN are highly consistent with current views that schizophrenia may be due to a disruption in neurodevelopmental processes. RELN encodes a large glycoprotein that affects both early and late neurodevelopment. It is responsible for neuronal migration, cell-cell interactions, and positioning of proliferating neurons early in development, and it also modulates neuronal and synaptic plasticity throughout life.57

DISC1 was also found to have significant epistatic relationships with other genes, including PDE4B and ERBB4 in addition to RELN. DISC1 is a well-established schizophrenia vulnerability gene.5962 It has multiple SNPs that have been extensively studied, including one identified in the present study (rs3738401).6369 The prevailing model for the role of DISC1 in schizophrenia emphasizes haploinsufficiency. DISC1 has decreased expression in mRNA and protein in both post mortem tissue and lymphoblastoid cell lines from patients with schizophrenia.66 DISC1 functions are consistent with both the neurodevelopmental hypothesis and current thinking about the mechanisms of schizophrenia. It is a hub protein whose interactors modulate expression of neurodevelopmental, synaptogenic, and sensory perception genes.6162 Its established interactome includes PDE4B, PDE4D, NDE1, and NDEL1; however, pathway analysis suggests that it may in fact be much larger.70 The DISC1 SNP identified in this analysis, rs3738401, is among those identified as a psychoactive drug target for the treatment of psychosis and/or schizophrenia through an Ingenuity pathway analysis.70

PDE4B was identified as having an epistatic relationship with DISC1 in our analyses, as well as with RELN and ERBB4. PDE4B has also been identified as a risk factor for schizophrenia.4849,53 It is a large gene (580 kb) with 17 exons; at least five isoforms occur as a result of alternative splicing. PDE4B2 and PDE4B4 isoforms have been found to be reduced in post mortem tissue from schizophrenia patients.48 It is orthologous to Drosophila dunce, a gene involved in learning and memory; dunce mutants have altered axonal grown cone motility, synaptic plasticity, and neuronal function.47 It also influences myelination, consistent with the findings in this study of significant WM loss associated with PDE4B genotypes. Its interaction with DISC1 occurs when the N-terminus of DISC1 binds to the UCR2 regulatory domain of PDE4B.71 This interactome plays a critical role in CNS signaling and homeostasis. DISC1 sequesters PDE4B in an inactive state and releases it when cyclicAMP is elevated; it thereby regulates intracellular signaling by inactivating cAMP.

NRG1 was found to have an epistatic relationship with RELN in our study. It has been implicated in multiple studies as a potential susceptibility gene for schizophrenia.7276 NRG1 is a large gene (1200 kb) that encodes at least 15 isoforms through alternative splicing.75 It is responsible for a host of functions that affect brain development and cell-cell communications and that are relevant to schizophrenia; they include axon guidance, myelination, glial differentiation, synaptogenesis, synaptic plasticity, and neurotransmission.75 mRNA expression studies implicate variation in isoform expression, due to alternative splicing, as a the most likely molecular mechanism for the genetic association of NRG1 with schizophrenia.77

ERBB4 would be predicted to have an epistatic relationship with NRG1, based on the fact that ERBB4 encodes a receptor for NRG1,75,78 but none was found in our analyses. It did have relationships with PDE4B and DISC 1, which are consistent with the hub protein nature of DISC1. ERBB4 is part of a family of genes for protein tyrosine kinases.75 We found an interaction with a known RELN SNP that is related to cognitive function. This relationship is limited exclusively to two basal ganglia structures, the caudate and putamen. Unlike other components of brain parenchyma, these two structures have increased in size. The epistatic relationship may provide some clue as to the mechanisms of antipsychotic action and, perhaps, toxicity.

This study suggests that exploratory analyses using MLAs may be a fruitful approach for examining the relationships between multiple genes that may be implicated in schizophrenia. It is noteworthy that five of the SNPs identified through this approach have been previously identified as risk alleles, in LD with risk alleles, or relevant to functions that are disrupted in schizophrenia. This replication of prior findings is of particular interest because the prior studies have generally been based on case-control comparisons, while the current study is a case only study, and brain measures rather than diagnosis were the outcome variable. This study also suggests that examination of genomic factors in schizophrenia may be enhanced by using a novel phenotypic outcome measure: brain change over time as measured by MR imaging.

This study has several limitations. It is a case-only sample; an examination of the genetic associations with brain change in a case-control setting would strengthen these results. It is based on a relatively small sample and is primarily useful for hypothesis generation. Because it is an exploratory study, we did not correct for multiple comparisons. The findings need replication in an independent sample. Finally, from a genetic perspective, we know that copy number variations (CNVs) and epigenetic factors such as methylation play an important role in conferring liability to disease. We have not examined any of these factors in this study.

Supplementary Material

Acknowledgments

This paper was written with support from the following grants: MHCRC: Neurobiology and Phenomenology of the Major Psychoses (MH43271); Phenomenology and the Classification of Schizophrenia (5R01MH031593); MR Imaging in the Major Psychoses (5R01MH040856); and BRAINS Morphology and Image Analysis (5R01NS050568); Investigator Initiated Research grant from Ortho-McNeil Janssen Scientific Affairs: Genetic Predictors of Long-term Course and Outcome in Schizophrenia.

Drs. Andreasen, Wassink, and Ho report receiving research support from Ortho-McNeil Janssen Scientific Affairs. Marsha Wilcox is an employee of Ortho-McNeil Janssen.

Footnotes

The other authors have no conflicts to report.

Supplementary information is available at Molecular Psychiatry’s website.

References

1. Cordell HJ, Clayton DG. Genetic association studies. Lancet. 2005 Sep 24–30;366(9491):1121–1131. [PubMed]
2. Andreasen NC. Pieces of the schizophrenia puzzle fall into place. Neuron. 1996;16(4):697–700. [PubMed]
3. Lencz T, Morgan TV, Athanasiou M, Dain B, Reed CR, Kane JM, et al. Converging evidence for a pseudoautosomal cytokine receptor gene locus in schizophrenia. Mol Psychiatry. 2007 Jun;12(6):572–580. [PubMed]
4. Sun S, Wang F, Wei J, Cao LY, Wu GY, Lu L, et al. Association between interleukin-3 receptor alpha polymorphism and schizophrenia in the Chinese population. Neurosci Lett. 2008 Jul 25;440(1):35–37. [PubMed]
5. Sun S, Wei J, Li H, Jin S, Li P, Ju G, et al. A family-based study of the IL3RA gene on susceptibility to schizophrenia in a Chinese Han population. Brain Res. 2009 May 1;1268:13–16. [PubMed]
6. Sullivan PF, Lin D, Tzeng JY, van den Oord E, Perkins D, Stroup TS, et al. Genomewide association for schizophrenia in the CATIE study: results of stage 1. Mol Psychiatry. 2008 Jun;13(6):570–584. [PubMed]
7. O'Donovan MC, Craddock N, Norton N, Williams H, Peirce T, Moskvina V, et al. Identification of loci associated with schizophrenia by genome-wide association and follow-up. Nat Genet. 2008 Sep;40(9):1053–1055. [PubMed]
8. Esslinger C, Walter H, Kirsch P, Erk S, Schnell K, Arnold C, et al. Neural mechanisms of a genome-wide supported psychosis variant. Science. 2009 May 1;324(5927):605. [PubMed]
9. Riley B, Thiselton D, Maher BS, Bigdeli T, Wormley B, McMichael GO, et al. Replication of association between schizophrenia and ZNF804A in the Irish Case-Control Study of Schizophrenia sample. Mol Psychiatry. 2010 Jan;15(1):29–37. [PMC free article] [PubMed]
10. Purcell SM, Wray NR, Stone JL, Visscher PM, O'Donovan MC, Sullivan PF, et al. Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature. 2009 Aug 6;460(7256):748–752. [PubMed]
11. Steinberg S, Mors O, Borglum AD, Gustafsson O, Werge T, Mortensen PB, et al. Expanding the range of ZNF804A variants conferring risk of psychosis. Mol Psychiatry. 2010 Jan 5; [PMC free article] [PubMed]
12. Need AC, Ge D, Weale ME, Maia J, Feng S, Heinzen EL, et al. A genome-wide investigation of SNPs and CNVs in schizophrenia. PLoS Genet. 2009 Feb;5(2) e1000373. [PMC free article] [PubMed]
13. Shi J, Levinson DF, Duan J, Sanders AR, Zheng Y, Pe'er I, et al. Common variants on chromosome 6p22.1 are associated with schizophrenia. Nature. 2009 Aug 6;460(7256):753–757. [PMC free article] [PubMed]
14. Stefansson H, Ophoff RA, Steinberg S, Andreassen OA, Cichon S, Rujescu D, et al. Common variants conferring risk of schizophrenia. Nature. 2009 Aug 6;460(7256):744–747. [PMC free article] [PubMed]
15. Athanasiu L, Mattingsdal M, Kahler AK, Brown A, Gustafsson O, Agartz I, et al. Gene variants associated with schizophrenia in a Norwegian genome-wide study are replicated in a large European cohort. J Psychiatr Res. 2010 Feb 23; [PMC free article] [PubMed]
16. Holmans P, Green EK, Pahwa JS, Ferreira MAR, Purcell SM, et al. Gene ontology analysis of GWA study data sets provides insights into the biology of bipolar disorder. Am J Med Genet. 2009 Jul;85:13–24. [PubMed]
17. Nicodemus KK, Callicott JH, Higier RG, Luna A, Nixon DC, Lipska BK, et al. Evidence of statistical epistasis between DISC1, CIT and NDEL1 impacting risk for schizophrenia: biological validation with functional neuroimaging. Hum Genet. 2010 Jan 19; [PubMed]
18. Breiman L, Friedman JH, Olshen RA, Stone CJ. Classification and regression trees. Boca Raton: Chapman and Hall; 1984.
19. Ellison-Wright I, Glahn DC, Laird AR, Thelen SM, Bullmore E. The anatomy of first-episode and chronic schizophrenia: an anatomical likelihood estimation meta-analysis. Amer J Psychiat. 2008;165:1015–1023. [PMC free article] [PubMed]
20. Nopoulos P, Torres I, Flaum M, Andreasen NC, Ehrhardt J, Yuh WTC. Brain morphology in first episode schizophrenia. Amer J Psychiatry. 1995;152:1721–1723. [PubMed]
21. Ho BC, Andreasen NC, Nopoulos P, Arndt S, Magnotta V, Flaum M. Progressive structural brain abnormalities and their relationship to clinical outcome. Arch Gen Psychiatry. 2003;60:585–594. [PubMed]
22. Lim KO, Tew W, Kushner M, Chow K, Matsumoto B, DeLisi L. Cortical gray matter volume deficit in first episode schizophrenia. Am J Psychiatry. 1996;153:1548–1553. [PubMed]
23. Kasai H, Shenton ME, Salisbury DF, et al. Progressive decrease of left Heschl gyrus and planum temporale gray matter volume in first-episode schizophrenia: a longtitudinal magnetic resonance imaging study. Arch Gen Psychiatry. 2003;60:766–775. [PMC free article] [PubMed]
24. Arango C, Moreno C, Martinez S, et al. Longitudinal brain changes in early-onset psychosis. Schizophr Bull. 2008 Mar;34(2):341–353. [PMC free article] [PubMed]
25. Degreef G, Ashtari M, Wu HW, Borenstein M, Geisler S, Lieberman J. Follow up MRI study in first episode schizophrenia. Schizophr Res. 1991 Oct;5(3):204–206. [PubMed]
26. DeLisi LE, Sakuma M, Maurizio AM, Relja M, Hoff AL. Cerebral ventricular change over the first 10 years after the onset of schizophrenia. Psychiatry Res. 2004 Jan 15;130(1):57–70. [PubMed]
27. DeLisi LE, Sakuma M, Tew W, Kushner M, Hoff AL, Grimson R. Schizophrenia as a chronic active brain process: a study of progressive brain structural change subsequent to the onset of schizophrenia. Psychiatry Res. 1997 Jul 4;74(3):129–140. [PubMed]
28. Gur RE, Cowell P, Turetsky BI, et al. A follow-up magnetic resonance imaging study of schizophrenia. Relationship of neuroanatomical changes to clinical and neurobehavioral measures. Arch Gen Psychiatry. 1998 Feb;55(2):145–152. [PubMed]
29. Lieberman J, Chakos M, Wu H, et al. Longitudinal study of brain morphology in first episode schizophrenia. Biol Psychiatry. 2001 Mar 15;49(6):487–499. [PubMed]
30. Whitford TJ, Grieve SM, Farrow TF, et al. Progressive grey matter atrophy over the first 2–3 years of illness in first-episode schizophrenia: a tensor-based morphometry study. Neuroimage. 2006 Aug 15;32(2):511–519. [PubMed]
31. Cahn W, Hulshoff Pol HE, Lems EB, et al. Brain volume changes in first-episode schizophrenia: a 1-year follow-up study. Arch Gen Psychiatry. 2002;59(11):1002–1010. [PubMed]
32. Ho BC, Andreasen NC, Dawson JD, Wassink TH. Association between brain-derived neurotrophic factor Val66Met gene polymorphism and progressive brain volume changes in schizophrenia. Am J Psychiatry. 2007;164(12):1890–1899. [PMC free article] [PubMed]
33. Andreasen NC, Nopoulos P, Magnotta V, Pierson R, Ziebell S, et al. Progressive brain change in schizophrenia: a prospective longitudinal study of first episode patients. Biol Psych. In Press. [PMC free article] [PubMed]
34. Weinberger DR. Implications of normal brain development for the pathogenesis of schizophrenia. Arch Gen Psychiatry. 1987 Jul;44(7):660–669. [PubMed]
35. Woods BT. Is schizophrenia a progressive neurodevelopmental disorder? Toward a unitary pathogenetic mechanism. Am J Psychiatry. 1998 Dec;155(12):1661–1670. [PubMed]
36. Feinberg I. Schizophrenia: caused by a fault in programmed synaptic elimination during adolescence? J Psychiatr Res. 1982;17(4):319–334. [PubMed]
37. Selemon LD, Kleinman JE, Herman MM, Goldman-Rakic PS. Smaller frontal gray matter volume in postmortem schizophrenic brains. Am J Psychiatry. 2002 Dec;159(12):1983–1991. [PubMed]
38. Benes FM, Turtle M, Khan Y, Farol P. Myelination of a key relay zone in the hippocampal formation occurs in the human brain during childhood, adolescence, and adulthood. Arch Gen Psychiatry. 1994;51:874–887. [PubMed]
39. Andreasen NC, Cizadlo T, Harris G, et al. Voxel processing techniques for the antemortem study of neuroanatomy and neuropathology using magnetic resonance imaging. J Neuropsychiatry Clin Neurosci. 1993;5(2):121–130. [PubMed]
40. Andreasen NC, Rajarethinam R, Cizadlo T, et al. Automatic atlas-based volume estimation of human brain regions from MR images. J Comput Assist Tomogr. 1996;20(1):98–106. [PubMed]
41. Harris G, Andreasen NC, Cizadlo T, et al. Improving tissue classification in MRI: a three-dimensional multispectral discriminant analysis method with automated training class selection. J Comput Assist Tomogr. 1999 Jan-Feb;23(1):144–154. [PubMed]
42. Magnotta VA, Andreasen NC, Schultz SK, et al. Quantitative in vivo measurement of gyrification in the human brain: changes associated with aging. Cereb Cortex. 1999 Mar;9(2):151–160. [PubMed]
43. Magnotta VA, Heckel D, Andreasen NC, et al. Measurement of brain structures with artificial neural networks: two- and three-dimensional applications. Radiology. 1999 Jun;211(3):781–790. [PubMed]
44. Pierson R, Johnson H, Harris G, Keefe H, Paulsen JS, Andreasen NC, Magnotta VA. Fully automated analysis using BRAINS. Neuroimage. 2010 [PubMed]
45. Wedenoja J, Tuulio-Henriksson A, Suvisaari J, Loukola A, Paunio T, Partonen T, et al. Replication of association between working memory and Reelin, a potential modifier gene in schizophrenia. Biol Psychiatry. 2010 May 15;67(10):983–991. [PMC free article] [PubMed]
46. van Schijndel JE, van Loo KM, van Zweeden M, Djurovic S, Andreassen OA, Hansen T, et al. Three-cohort targeted gene screening reveals a non-synonymous TRKA polymorphism associated with schizophrenia. J Psychiatr Res. 2009 Oct;43(15):1195–1199. [PubMed]
47. Chubb JE, Bradshaw NJ, Soares DC, Porteous DJ, Millar JK. The DISC locus in psychiatric illness. Mol Psychiatry. 2008 Jan;13(1):36–64. [PubMed]
48. Fatemi SH, King DP, Reutiman TJ, Folsom TD, Laurence JA, Lee S, et al. PDE4B polymorphisms and decreased PDE4B expression are associated with schizophrenia. Schizophr Res. 2008 Apr;101(1–3):36–49. [PubMed]
49. Millar JK, Wilson JC, Anderson S, Christie S, Taylor MS, et al. Disruption of two novel genes by a translocation co-segregating with schizophrenia. Hum Mol Genet. 2000;9:1415–1423. [PubMed]
50. Law AJ, Lipska BK, Weickert CS, Hyde TM, Straub RE, Hashimoto R, et al. Neuregulin 1 transcripts are differentially expressed in schizophrenia and regulated by 5' SNPs associated with the disease. Proc Natl Acad Sci U S A. 2006 Apr 25;103(17):6747–6752. [PubMed]
51. Benzel I, Bansal A, Browning BL, Galwey NW, Maycox PR, McGinnis R, et al. Interactions among genes in the ErbB-Neuregulin signalling network are associated with increased susceptibility to schizophrenia. Behav Brain Funct. 2007;3:31. [PMC free article] [PubMed]
52. Alaerts M, Ceulemans S, Forero D, Moens LN, De Zutter S, Heyrman L, et al. Support for NRG1 as a susceptibility factor for schizophrenia in a northern Swedish isolated population. Arch Gen Psychiatry. 2009 Aug;66(8):828–837. [PubMed]
53. Tomppo L, Hennah W, Laherma P, Loukola A, Tuulio-Henriksson A, et al. Association between genes of disrupted in schizophrenia 1 (DISC1) interactors and schizophrenia supports the role of the DISC1 pathway in the etiology of major mental illnesses. Biol Psychiatry. 2009;65:1055–1062. [PMC free article] [PubMed]
54. Baune BT, Konrad C, Suslow T, Domschke K, Birasova E, et al. The reelin (RELN) gene is associated with executive function in healthy individuals. Neurol Biol Learning and Memory. Epub 2010. [PubMed]
55. Grayson DR, Jia X, Chen Y, Sharma RP, Mitchell CP, Guidotti A, et al. Reelin promoter hypermethylation in schizophrenia. Proc Natl Acad Sci U S A. 2005 Jun 28;102(26):9341–9346. [PubMed]
56. Lintas C, Persico AM. Neocortical RELN promoter methylation increases significantly after puberty. Neuroreport. 2010 Jan 27;21(2):114–118. [PubMed]
57. Goldberger C, Gourion D, Leroy S, Schurhoff F, Bourdel MC, Leboyer M, et al. Population-based and family-based association study of 5'UTR polymorphism of the reelin gene and schizophrenia. Am J Med Genet B Neuropsychiatr Genet. 2005 Aug 5;137B(1):51–55. [PubMed]
58. Impagnatiello F, Guidotti AR, Pesold C, Dwivedi Y, Caruncho H, Pisu MG, et al. A decrease of reelin expression as a putative vulnerability factor in schizophrenia. Proc Natl Acad Sci U S A. 1998 Dec 22;95(26):15718–15723. [PubMed]
59. Millar JK, Wilson-Annan JC, Anderson S, Christie S, Taylor MS, Semple CA, et al. Disruption of two novel genes by a translocation co-segregating with schizophrenia. Hum Mol Genet. 2000 May 22;9(9):1415–1423. [PubMed]
60. Ekelund J, Hennah W, Hiekkalinna T, Parker A, Meyer J, Lonnqvist J, et al. Replication of 1q42 linkage in Finnish schizophrenia pedigrees. Mol Psychiatry. 2004 Nov;9(11):1037–1041. [PubMed]
61. Chubb JE, Bradshaw NJ, Soares DC, Porteous DJ, Millar JK. The DISC locus in psychiatric illness. Mol Psychiatry. 2008 Jan;13(1):36–64. [PubMed]
62. Hennah W, Thomson P, Peltonen L, Porteous D. Genes and schizophrenia: beyond schizophrenia: the role of DISC1 in major mental illness. Schizophr Bull. 2006 Jul;32(3):409–416. [PMC free article] [PubMed]
63. Eastwood SL, Hodgkinson CA, Harrison PJ. DISC-1 Leu607Phe alleles differentially affect centrosomal PCM1 localization and neurotransmitter release. Mol Psychiatry. 2009 Jun;14(6):556–557. [PubMed]
64. Nakata K, Lipska BK, Hyde TM, Ye T, Newburn EN, Morita Y, et al. DISC1 splice variants are upregulated in schizophrenia and associated with risk polymorphisms. Proc Natl Acad Sci U S A. 2009 Sep 15;106(37):15873–15878. [PubMed]
65. Hodgkinson CA, Goldman D, Jaeger J, Persaud S, Kane JM, Lipsky RH, et al. Disrupted in schizophrenia 1 (DISC1): association with schizophrenia, schizoaffective disorder, and bipolar disorder. Am J Hum Genet. 2004 Nov;75(5):862–872. [PubMed]
66. Callicott JH, Straub RE, Pezawas L, Egan MF, Mattay VS, Hariri AR, et al. Variation in DISC1 affects hippocampal structure and function and increases risk for schizophrenia. Proc Natl Acad Sci U S A. 2005 Jun 14;102(24):8627–8632. [PubMed]
67. Saetre P, Agartz I, De Franciscis A, Lundmark P, Djurovic S, Kahler A, et al. Association between a disrupted-in-schizophrenia 1 (DISC1) single nucleotide polymorphism and schizophrenia in a combined Scandinavian case-control sample. Schizophr Res. 2008 Dec;106(2–3):237–241. [PubMed]
68. Liu YL, Fann CS, Liu CM, Chen WJ, Wu JY, Hung SI, et al. A single nucleotide polymorphism fine mapping study of chromosome 1q42.1 reveals the vulnerability genes for schizophrenia, GNPAT and DISC1: Association with impairment of sustained attention. Biol Psychiatry. 2006 Sep 15;60(6):554–562. [PubMed]
69. Szeszko PR, Hodgkinson CA, Robinson DG, Derosse P, Bilder RM, Lencz T, et al. DISC1 is associated with prefrontal cortical gray matter and positive symptoms in schizophrenia. Biol Psychol. 2008 Sep;79(1):103–110. [PMC free article] [PubMed]
70. Hennah W, Porteous D. The DISC1 pathway modulates expression of neurodevelopmental, synaptogenic and sensory perception genes. PLoS One. 2009;4(3):e4906. [PMC free article] [PubMed]
71. Millar JK, Pickard BS, Kackie S, James R, Christie S, et al. DISC1 and PDE4B are interacting genetic factors in schizophrenia that regulate cAMP signaling. Science. 2005;310:1187–1191. [PubMed]
72. Stefansson H, Sarginson J, Kong A, Yates P, Steinthorsdottir V, Gudfinnsson E, et al. Association of neuregulin 1 with schizophrenia confirmed in a Scottish population. Am J Hum Genet. 2003 Jan;72(1):83–87. [PubMed]
73. Stefansson H, Sigurdsson E, Steinthorsdottir V, Bjornsdottir S, Sigmundsson T, Ghosh S, et al. Neuregulin 1 and susceptibility to schizophrenia. Am J Hum Genet. 2002 Oct;71(4):877–892. [PubMed]
74. Harrison PJ, Law AJ. Neuregulin 1 and schizophrenia: genetics, gene expression, and neurobiology. Biol Psychiatry. 2006 Jul 15;60(2):132–140. [PubMed]
75. Mei L, Xiong WC. Neuregulin 1 in neural development, synaptic plasticity and schizophrenia. Nat Rev Neurosci. 2008 Jun;9(6):437–452. [PMC free article] [PubMed]
76. Hall J, Whalley HC, Job DE, Baig BJ, McIntosh AM, Evans KL, et al. A neuregulin 1 variant associated with abnormal cortical function and psychotic symptoms. Nat Neurosci. 2006 Dec;9(12):1477–1478. [PubMed]
77. Law AJ, Lipska BK, Weickert CS, Hyde TM, Straub RE, Hashimoto R, et al. Neuregulin 1 transcripts are differentially expressed in schizophrenia and regulated by 50027; SNPs associated with the disease. Proc Natl Acad Sci U S A. 2006 Apr 25;103(17):6747–6752. [PubMed]
78. Benzel I, Bansal A, Browning BL, Galwey NW, Maycox PR, McGinnis R, et al. Interactions among genes in the ErbB-Neuregulin signalling network are associated with increased susceptibility to schizophrenia. Behav Brain Funct. 2007;3:31. [PMC free article] [PubMed]