Document Detail

Two gene co-expression modules differentiate psychotics and controls.
Jump to Full Text
MedLine Citation:
PMID:  23147385     Owner:  NLM     Status:  Publisher    
Abstract/OtherAbstract:
Schizophrenia (SCZ) and bipolar disorder (BD) are highly heritable psychiatric disorders. Associated genetic and gene expression changes have been identified, but many have not been replicated and have unknown functions. We identified groups of genes whose expressions varied together, that is co-expression modules, then tested them for association with SCZ. Using weighted gene co-expression network analysis, we show that two modules were differentially expressed in patients versus controls. One, upregulated in cerebral cortex, was enriched with neuron differentiation and neuron development genes, as well as disease genome-wide association study genetic signals; the second, altered in cerebral cortex and cerebellum, was enriched with genes involved in neuron protection functions. The findings were preserved in five expression data sets, including sets from three brain regions, from a different microarray platform, and from BD patients. From those observations, we propose neuron differentiation and development pathways may be involved in etiologies of both SCZ and BD, and neuron protection function participates in pathological process of the diseases.Molecular Psychiatry advance online publication, 13 November 2012; doi:10.1038/mp.2012.146.
Authors:
C Chen; L Cheng; K Grennan; F Pibiri; C Zhang; J A Badner; ; E S Gershon; C Liu
Related Documents :
15758235 - Diverse roles for hspr in campylobacter jejuni revealed by the proteome, transcriptome ...
3277955 - Roles of the escherichia coli heat shock sigma factor 32 in early and late gene express...
10455235 - Schizosaccharomyces pombe homologue of glutathione peroxidase, which does not contain s...
10567555 - Yap1p activates gene transcription in an oxidant-specific fashion.
18820255 - Evolution of soldier-specific venomous protease in social aphids.
11287605 - Identification of liver x receptor-retinoid x receptor as an activator of the sterol re...
Publication Detail:
Type:  JOURNAL ARTICLE     Date:  2012-11-13
Journal Detail:
Title:  Molecular psychiatry     Volume:  -     ISSN:  1476-5578     ISO Abbreviation:  Mol. Psychiatry     Publication Date:  2012 Nov 
Date Detail:
Created Date:  2012-11-13     Completed Date:  -     Revised Date:  -    
Medline Journal Info:
Nlm Unique ID:  9607835     Medline TA:  Mol Psychiatry     Country:  -    
Other Details:
Languages:  ENG     Pagination:  -     Citation Subset:  -    
Affiliation:
Department of Psychiatry, The University of Illinois at Chicago, Chicago, IL, USA.
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-journal-id): 9607835
Journal ID (pubmed-jr-id): 20545
Journal ID (nlm-ta): Mol Psychiatry
Journal ID (iso-abbrev): Mol. Psychiatry
ISSN: 1359-4184
ISSN: 1476-5578
Article Information
Download PDF

nihms-submitted publication date: Day: 14 Month: 4 Year: 2014
Electronic publication date: Day: 13 Month: 11 Year: 2012
Print publication date: Month: 12 Year: 2013
pmc-release publication date: Day: 01 Month: 6 Year: 2014
Volume: 18 Issue: 12
First Page: 1308 Last Page: 1314
PubMed Id: 23147385
ID: 4018461
DOI: 10.1038/mp.2012.146
ID: NIHMS560305

Two Gene Co-expression Modules Differentiate Psychotics and Controls
Chao Chen, Ph.D Affiliation: Psychiatry Department, the University of Illinois at Chicago, Chicago, IL, US 60607
Lijun Cheng, Ph.D Affiliation: Psychiatry Department, the University of Illinois at Chicago, Chicago, IL, US 60607
Kay Grennan, Ph.D Affiliation: Department of Psychiatry and Behavioral Neuroscience, The University of Chicago, Chicago, IL, US 60637
Fabio Pibiri, Ph.D Affiliation: Department of Pediatrics, the University of Illinois at Chicago, Chicago, IL, US 60607
Chunling Zhang, Ph.D Affiliation: Psychiatry Department, the University of Illinois at Chicago, Chicago, IL, US 60607
Judith A. Badner, M.D. Ph.D Affiliation: Department of Psychiatry and Behavioral Neuroscience, The University of Chicago, Chicago, IL, US 60637
Elliot S. Gershon, M.D. Affiliation: Department of Psychiatry and Behavioral Neuroscience, The University of Chicago, Chicago, IL, US 60637
Chunyu Liu, Ph.D* Affiliation: Psychiatry Department, the University of Illinois at Chicago, Chicago, IL, US 60607
*Correspondence: Dr Chunyu Liu, Associate Professor, Psychiatry Department, The University of Illinois at Chicago, 900 S Ashland Ave, M/C915 Chicago, IL 60607 liucy@uic.edu Tel: 312-413-2599; Fax: 312-355-1741

Introduction

Schizophrenia (SCZ) and bipolar disorder (BD) are psychiatric disorders with world-wide lifetime prevalences of around 1%1-4. These disorders have been shown to be highly heritable by monozygotic and dizygotic twin studies and by adoption studies: the heritability estimates for SCZ range from 70 to 85% and, for BD, from 60 to 85%5. Over the past two decades, genetic studies of thousands of samples have identified hundreds of candidate genes6-8. Most of these findings have not been supported by genome-wide studies. Moreover, functional genomic roles have not yet been determined for most of the associated genes.

Gene expression can bridge the gap between genetic variation and disease susceptibility as an intermediate phenotype that is regulated by a combination of genetic and epigenetic factors. Gene expression transcription profiling is widely used and has been thoroughly validated by the MicroArray Quality Control (MAQC) project9,10. Expression microarray studies test thousands of gene transcripts for differential expression simultaneously. When all these genes are tested for association with disease, it creates a multiple testing problem: to minimize false positives, the multiple test correction sets the significance threshold so high that true positives might be missed as well. Additionally, testing individual genes for disease association also ignores the interaction between genes. So far, hundreds of gene expression microarray studies of psychiatric diseases have been reported, including studies of SCZ, BD and autism. However, the significant gene changes detected by one study are seldom replicated in another, let alone across three or more11.

Gene expression network analyses are an alternative approach for analyzing expression data that reduce the sample space tested (that is, the number of hypotheses to be tested). Gene expression networks are constructed from expression data from thousands of genes, and describe the interactions among groups of transcripts. They can be used to observe systematic alterations in expression associated with complex diseases such as psychiatric disorders. Horvath and colleagues have developed an algorithm for creating networks, called Weighted Gene Co-expression Network Analysis (WGCNA), which is widely used12,13. This method identifies groups of genes within a network whose expressions are highly correlated. These groups, called modules, can then be compared between cases and controls, among different tissues, species, or other phenotypes or clinical traits14-17.

Oldham et al. were the first to apply this method to expression in brain. They compared the gene expression across different brain regions and demonstrated that the modules reflect the underlying cell-type composition of the regions18.

The method was first applied to a psychiatric disease by Torkamani et al., who detected SCZ-associated gene co-expression modules in one microarray data analysis, and found that aging affected gene expression in normal controls, but not in SCZ patients19. They made the interesting distinction between constructing networks from case and control data separately and constructing a network from a combination of case and control data. Constructing modules from case and control networks separately allows comparison of basic module structure between the two groups. If modules are not substantially preserved in cases, a conclusion can be drawn that fundamental relationships among genes have been disrupted. If modules are substantially preserved in cases, the case and control data can be combined for network construction, and the detected modules assessed for differential expression between cases and controls. Torkamani found that modules were substantially preserved in cases, and identified five SCZ-associated modules; Gene Ontogeny (GO) enrichment analysis showed those five modules were associated with oxidative phosphorylation, angiogenesis, neuron differentiation, chromatin and nucleosome assembly, and inositol phosphate metabolism, separately.

Voineagu et al. applied WGCNA to another psychiatric disease, autism. In addition to a GO enrichment analysis, they found one of the differentially expressed modules was enriched in disease genome-wide association study (GWAS) signals. They interpreted this as evidence that the module's member genes were causally associated with autism. They also found another module with altered expression but not enriched in GWAS signals and took it as an indication of a non-genetic etiology20.

Note that all three of these network studies were based on analysis of a single microarray data set; their findings have not yet been replicated in other data sets at the network level.

To demonstrate the robustness and reproducibility of WGCNA findings, we used multiple microarray data sets to study co-expression networks in brains of psychiatric patients. We constructed gene expression networks, and then identified gene co-expression modules within the networks. The modules were tested for association with SCZ. We assessed whether the disease-associated modules were detected in independent data sets, including sets from several different brain collections and different brain regions. Reproducible, or preserved, modules were then evaluated for case-control differences in the independent data sets.

Since it is widely accepted that SCZ and BD have genetic factors in common21,22, we tested whether the gene modules perturbed in SCZ were similarly perturbed in BD. We also tested whether the disease-associated modules were enriched with genetic variants that had been previously associated with disease by GWASs23,24.


Subjects and Methods
Samples and quality control

Cerebellum (CB) and parietal cortex (PCX) brain tissues were obtained from Stanley Medical Research Institute (SMRI)25,. They came from the SMRI's Neuropathology Consortium and Array collections, and included 50 SCZ samples, 50 BD samples, and 50 unaffected control samples. Expression data for these samples came from the NCBI Gene Expression Omnibus (GEO) database (GSE35978). One of the two prefrontal cortex (PFC) brain expression data sets came from Dobrin's group using SMRI samples (PFC-SMRI)25, while the second came from the Victorian Brain Bank Network (PFC-VBBN)26 and was obtained from at GEO (GSE21138, sample information and preparation in Supplementary File).

ComBat, a batch effects adjustment program that we have previously shown to be the best available, was used to remove batch effects from the data sets (see Supplementary File for detailed preprocessing steps)27,28.

Gene network construction and module detection

We used weighted gene co-expression network analysis (WGCNA)13 to identify modules of co-expressed genes within gene expression networks. To construct the network, the absolute values of Pearson correlation coefficients were calculated for all possible gene pairs. Values were entered into a matrix, and the data were transformed so that the matrix followed an approximate scale-free topology (see Supplementary File for detailed information). A dynamic tree cut algorithm was used to detect network modules29. WGCNA and the dynamic tree cut algorithm were implemented in R12,29.

We ran singular value decomposition (SVD) on each module's expression matrix and used the resulting module eigengene, which is equivalent to the first principal component 30, to represent the overall expression profiles of the modules.

Module preservation statistics

We utilized the module preservation statistic Zsummary, to assess the module preservation from different expression data sets33 (see Supplementary File for formal definition). Unlike the cross-tabulation test, Zsummary not only takes into account the overlap in module membership, but also the density and connectivity patterns of modules. In addition, for our study, network-based preservation statistics only require that module membership be identified in the original data sets, reducing the variation coming from various parameters setting to identify new modules in validation data sets.

We converted the probe-level measurements into gene-level measurements to make data from different platforms comparable. The probe within a gene that had the highest coefficient of variation was used to represent that gene. Overall, 8497 genes were retained in our preservation calculation.

Differential expression test for modules

We used multiple linear regressions on the modules’ eigengenes to remove the effects of sex, age, pH and PMI from the SMRI samples and VBBN samples. The residual eigengene values were then used to test the disease association using Pearson's correlation test. We used false discovery rate (FDR) for multiple testing correction49.

Module-based GWAS signal enrichment test

GAIN-SCZ was the genome-wide association study of the GAIN SCZ data set, which comprised 4591 cases and controls (1217 European-American cases, 1442 European-American controls, 953 African-American cases, 979 African-American controls). The data was downloaded from dbGaP (http://www.ncbi.nlm.nih.gov/projects/gap/cgibin/study.cgi?study_id=phs000021.v3.p2). GAIN-BD was a genome-wide association study of the GAIN BD data set, which comprised 3261 cases and controls (1079 European-American cases, 1081 European-American controls, 415 African-American cases, 686 African-American controls). The data was downloaded from dbGaP (http://www.ncbi.nlm.nih.gov/projects/gap/cgibin/study.cgi?study_id=phs000017.v3.p1). Whole genome genotyping of GAIN data was done with the Affymetrix Genome-Wide Human SNP Array 6.0.

The TGen-BD GWAS data came from the genome-wide association study of the Translational Genomics Research Institute's (TGen) BD data (http://www.tgen.org/), which comprised 1,190 BD cases and 401 controls. Sample genotyping was conducted using the Affymetrix GeneChip Mapping 5.0K Array. We performed imputation using MaCH v1.0 to increase the density of interrogated SNPs31, with HapMap data as reference. Overall, 2,593,107 SNPs in GAIN-SCZ, 3,281,319 SNPs in GAIN-BD, and 2,542,706 SNPs in TGen-BD were included after imputation.

Imputed GWAS data was used to test whether the two disease-associated modules were enriched in SCZ/BD association signals32. We used a previously-reported procedure to run the enrichment test32. First, the max −log(P-value) of a SNP located between 20kb upstream and downstream of a gene was assigned to represent the gene, then the gene set's enrichment scores (ES) were calculated based the gene's rank. SNP level permutation was applied to generate the distribution of the ES and then the distribution was normalized. For multiple gene sets, FDR was calculated by joining all the distributions of ESs, each for one gene set, generated by permutation. As the difference of gene length distribution between genes in one module was significant (p=5.47e-27), we also applied a permutation procedure to verify whether there was a gene length bias in the genetic signals enrichment test43, 50, and there was no bias (see Supplementary Methods section 8).


Results
Gene modules in parietal cortex

We first analyzed the parietal cortex (PCX) brain gene expression in 50 SCZ patients and 50 normal controls from the Stanley Medical Research Institute (SMRI) using the Affymetrix Human Gene 1.0 ST Array25. After a series of sample and array-level quality control measures (see Supplementary methods for details), we retained 45 SCZ patients and 46 normal controls with measures of 19,884 transcripts.

Networks were constructed in two different ways: first, we constructed one network from control data and one network from case data, then identified modules in the control network and assessed their preservation in the case network. Second, we constructed a network from the combined case and control data, and identified modules within it.

Structure of co-expression modules in schizophrenia cases and controls

Case and control sample gene networks were constructed separately to detect whether there was any gene co-regulation disruption or creation in cases relative to controls. We assessed module preservation using a permutation-based preservation statistic, Zsummary, developed by Langfelder et al., which assesses whether the connectivity level and pattern of a module in one data set is preserved in another33. The authors suggest the following significance thresholds: Zsummary < 2 implies no evidence for module preservation, 2 < Zsummary < 10 implies weak to moderate evidence, and Zsummary > 10 implies strong evidence for module preservation. In our study, all modules detected in the control data had Zsummary greater than 10 in the case data (Fig. 1), suggesting well-preserved membership and connectivity of the control modules in the SCZ cases. This is consistent with two previous gene network-based psychosis studies20,34, where the case modules had no obvious perturbations relative to control modules.

To test the reliability of the module construction results, we compared the modules identified from control samples to published gene expression networks. Cross-tabulation showed that our control modules have similar module membership to modules previously reported by Oldham et al. 18(Supplementary Fig. 1). Slight differences, including some modules in our controls collapsing into one module in the Oldham data, might have been due to the differences in samples and platforms, data-filtering steps and/or parameters used for network construction and module detection.

Differential expression of modules in schizophrenia patients

Since there was no significant difference in module structure between cases and controls, modules identified in a gene network constructed from both case and control data were analyzed for differential expression in SCZ patients. Twenty-four modules were detected (Supplementary Fig. 2). We used an eigengene to summarize each module's expression profile30 (see Methods for formal eigengene definition). After multiple linear regressions on the eigengenes to remove the effects of sex, pH and post-mortem interval (PMI), three covariates that can affect measures of gene expression11, the residual module eigengenes for each individual were tested for disease association. The eigengenes of two PCX modules, referred to as M1A and M3A, were significantly associated with SCZ after multiple testing correction (Fig. 2A). M1A comprised 490 genes, while M3A comprised 106. NOTCH2 and MT1X, respectively, were the hub genes of M1A and M3A, meaning that, of all the genes in their modules, these two genes had the strongest correlation with the module's eigengene, as well as being highly connected to the other genes in their modules (Supplementary Tables 5 and 6).


Replication of findings in different data sets
Replication of module structure

We tested co-expression module preservation across different microarray data sets. We again used the module preservation statistic, Zsummary, to compare modules. Data sets included prefrontal cortex (PFC) tissues and cerebellum (CB) tissues, also from the Stanley Medical Research Institute (SMRI) but performed by a different research group25, and another PFC data set, where samples came from the Victorian Brain Bank Network (VBBN)26.

As measured by the Zsummary statistic, six of the 18 modules identified in the SMRI PCX data were strongly preserved in the SMRI PFC data set, including M1A. Six modules were moderately preserved, including M3A, and four modules were not preserved (Supplementary Fig. 3). Similarly, M1A was strongly preserved and M3A moderately preserved in the SMRI CB and VBBN PFC data sets (Table 1; Supplementary Fig. 4 and 5).

Replication of disease association

We tested the modules’ hub genes for disease association, since hub genes by definition are highly correlated with their modules’ eigengenes35. After controlling for the effects of covariates, NOTCH2, the hub gene of M1A, was significantly associated with SCZ in VBBN and SMRI PFC data (p=0.049 and 0.022, respectively), but not in SMRI CB data (p=0.583). MT1X, the hub gene of M3A, was consistently upregulated in SCZ in VBBN PFC, SMRI PFC and SMRI CB data (p=3.5E-03, 0.026 and 1.3E-04, respectively) (Fig. 2B, 2C). These results suggest that M1A gene network perturbation between SCZ and controls occurs in cerebral cortex but not in cerebellum, while M3A perturbation happens across all brain regions.

Modules shared between schizophrenia and bipolar disorder

We were interested in whether the M1A and M3A modules associated with SCZ were also associated with BD. One BD data set with two regions studied was tested, PCX and CB from SMRI. The control samples here overlapped with the control samples used to construct the SCZ case and control network and modules. We first confirmed that the M1A and M3A in PCX were well-preserved in the BD case-control data sets (Table 1; Supplementary Fig. 6 and 7), then we tested whether NOTCH2 and MT1X gene expression were associated with BD. NOTCH2 showed an association trend in PCX (p=0.054) but not in CB (p=0.117), while MT1X showed significant association in both regions (p=0.015 in PCX and 1.5e-3 in CB) (Fig. 2B, 2C), indicating M1A in PCX and M3A in two brain regions were altered in BD, and were shared by two diseases.

Characterization of two disease-associated modules

NOTCH2, a member of the Notch gene family, had the highest intramodular connectivity in the SMRI PCX M1A module. The Notch signaling pathway plays a key role in cell to cell communication, and was reported to play multiple roles in central nervous system development36-38. Functional enrichment analyses were performed with The Database for Annotation, Visualization and Integrated Discovery (DAVID) (Fig. 3, Supplementary Table 7)39. Consistent with the hub gene's function, M1A was highly significantly enriched with gene ontology terms related to neuron differentiation and neuron development genes (p=7.50e-08 and 3.63e-06, FDR q=1.27e-04 and 6.14e-03, respectively).

We also found that 46 genes in M1A overlapped with a synaptic gene group associated with SCZ identified by Lips et al.40. These overlapping genes are interesting because ten of them have previously been reported to be associated with either SCZ or BD, which may suggest a relationship between synaptic dysfunction and these two psychiatric diseases. These genes include ABLIM1, APOE, AQP4, SLC1A2, SLC1A3, SLC4A4, GABRG1, NTRK2, ADD3 and MAOA.

Functional enrichment analysis showed M3A to be highly enriched in gene ontology categories related to metallothioneins (MT) and metal binding site functions (p=2.41e-05, FDR q=0.031) (Fig. 4, Supplementary Table 8). MT genes respond to several stimuli, including metals, oxidative agents, inflammation and stress; they influence cognition, protect against neurotoxicity and play a role in the astrocytic response to CNS injuries41. The expressions of our MT gene family members, including the hub gene MT1X, and MT1E, MT1F and MT2A have previously been reported to be differentially expressed in SCZ42.

Enrichment of genetic association signals

To determine whether the expression alterations to modules M1A and M3A had genetic bases, we tested the modules for enrichment with SCZ and BD genetic association signals. M1A showed significant enrichment of signals from both SCZ and BD GWASs (GAIN, SCZ, p<0.001, FDR q=0.009; GAIN-BD, p<0.001, FDR q=0.0015; TGen, BD, p<0.001, FDR q=0.002) (Supplementary Fig. 8, 9 and 10) (see Supplementary materials for further descriptions of GWAS data sets). Those significances were retained after controlling for gene length bias43. M3A was not enriched in signals from either disease. We also used other GWAS data, from non-psychiatric disease type 2 diabetes (T2D), as a negative control for our enrichment analysis. Neither M1A nor M3A were enriched with signals from a T2D GWAS.


Discussion

We looked for case-control differences at the level of gene co-expression regulation, rather than the level of individual genes. We found strong evidence of such differences in one expression data set, and then demonstrated robust gene network preservation and disease association across five additional different data sets. These data sets came from different brain banks and different brain regions, and were produced on different expression microarray platforms; they contained data from patients with SCZ and BD, as well as normal controls. To the best of our knowledge, this is the first time gene expression alterations associated with disease have been replicated across so many gene expression data sets. This reproducibility suggests that our results are reliable and that the WGCNA network-based method increased the statistical power of our expression study.

M1A was enriched with neuron development and neuron differentiation genes, among others. Using module preservation statistics, we determined that the M1A module was highly preserved in the five other data sets, and was consistently up-regulated in SCZ and BD cerebral cortex relative to controls, but not in cerebellum. This result is consistent with Torkamani et al's finding that modules enriched in neuron differentiation are differentially expressed between SCZ patients and normal controls19. A test for enrichment with GWAS signals suggested that M1A expression changes may have a genetic basis. Based on those observations, we speculate that the etiology of psychiatric diseases could be related to developmental dysfunction in the cerebral cortex, extending the neurodevelopmental hypothesis of SCZ44, which has recently received further support from imaging studies45 and deep sequencing analysis46.

M3A was enriched with metallothionein genes. It was moderately conserved in other data sets and consistently showed differential expression in both cerebral cortex and cerebellum between SCZ and BD patients, relative to controls. M3A was not significantly enriched in GWAS signals. Expression changes in these MT genes in SCZ and BD has been previously reported42,47, but, so far, no positive genetic associations of these genes have been reported. Considering that the major physiological function of MT genes is binding of heavy metals to protect the central nervous system (CNS) against neurotoxicity or other injuries41, the change in the expression of M3A may be a downstream effect of disease or may be related to environmental insult. Alternatively, loci regulating M3A or its member genes may exist, but not have been detected by GWAS.

One advantage of module-based association studies is that they substantially reduce the multiple testing correction burden inherent in genome-wide association studies51 and microarray studies, particularly in brain expression studies where sample size is always limited by tissue availability. Microarray studies test tens of thousands of genes at one time, which require that p-values be very small to reach genome-wide significance. However, WGCNA produced tens of modules of highly correlated genes to be tested for association. In our case, we found M1A to be significant in most of our data sets, but when we ran an individual gene-based analysis, the hub gene of M1A, NOTCH2, was nominally significant but did not survive the multiple testing correction (data not shown). WGCNA may also reveal the function of an otherwise poorly characterized gene, if that gene is a member of a module highly enriched in a particular function. In addition, by comparing results from different data sets, we demonstrated that findings at the network/pathway level are consistent.

Our results suggested that the observed similarities between SCZ and BD are manifested at the transcriptional level. Shared genetic elements have previously been demonstrated21,22,48, and some expression studies reported overlapping aberrant individual genes at the transcriptional level47. Our study explored the shared pathology defect at a higher level, revealing a larger scale gene pathway-based perturbation. Further functional work can focus on those genes, especially the hub genes, which may reveal the molecular bases of the symptoms shared between these two diseases.

In summary, our study revealed that 1) SCZ and BD shared altered expression of neuron differentiation and neuron development genes (M1A) and metallothionein genes (M3A), which may indicate common etiology or pathology; 2) M1A's up-regulation in cerebral cortex could be regulated by genetic variants already found to be associated with SCZ and BD risk; 3) M3A expression change was found in both cerebral cortex and cerebellum, and was not detectably regulated by genetic variants.



Notes

FN1Conflict of interest:

The authors declare no conflict of interest.

FN2Supplementary information is available at Molecular Psychiatry's website.

Acknowledgement

This work was funded by 5R01MH080425 (to C.L.), and supported by Eklund Family. We thank Elizabeth Thomas and Ali Torkamani for sharing Victorian Brain Bank Network (VBBN) brain expression data with us. We thank investigators of Wellcome Trust Case Control Consortium (WTCCC) for sharing non-psychiatric disease type 2 diabetes (T2D) GWAS data. We thank the investigators of Bipolar Genome Study (BiGS) for generating those bipolar GAIN-BD and Tgen-BD GWAS data, thank investigators of GAIN for generating SZ GWAS data. We thank Seth E. Dobrin, Maree Webster and other collaborators at the Stanley Medical Research Institute for providing the Stanley Online Genomics Database and brain samples. We also thank Lorenzo Pesce, Alex Rodriguez and collaborators at Computation Institute, University of Chicago for supplying their supercomputing devices and services.


References
1. Jablensky A. The 100-year epidemiology of schizophrenia.Schizophrenia ResearchDay: 19 Month: 12 Year: 1997282-31111259468347
2. McGuffin P,Rijsdijk F,Andrew M,Sham P,Katz R,Cardno A. The heritability of bipolar affective disorder and the genetic relationship to unipolar depression.Archives of General PsychiatryMonth: 5 Year: 200360549750212742871
3. Merikangas KR,Akiskal HS,Angst J,Greenberg PE,Hirschfeld RMA,Petukhova M,et al. Lifetime and 12-month prevalence of bipolar spectrum disorder in the national comorbidity survey replication.Archives of General PsychiatryMonth: 5 Year: 200764554355217485606
4. McGrath J,Saha S,Chant D,Welham J. Schizophrenia: A Concise Overview of Incidence, Prevalence, and Mortality.Epidemiologic ReviewsDay: 1 Month: 11 Year: 2008301677618480098
5. Burmeister M,McInnis MG,Zollner S. Psychiatric genetics: progress amid controversy.Nature Reviews GeneticsMonth: 7 Year: 200897527540
6. Sanders J,Gill M. Unravelling the genome: a review of molecular genetic research in schizophrenia.Irish Journal of Medical ScienceMonth: 3 Year: 200717615917849516
7. Hattori E,Liu CY,Badner JA,Bonner TI,Christian SL,Maheshwari M,et al. Polymorphisms at the G72/G30 gene locus, on 13q33, are associated with bipolar disorder in two independent pedigree series.American Journal of Human GeneticsMonth: 5 Year: 20037251131114012647258
8. O'Donovan MC,Williams HJ,Owen MJ. New findings from genetic association studies of schizophrenia.Journal of Human GeneticsMonth: 1 Year: 200954191419158819
9. Shi LM,Campbell G,Jones WD,Campagne F,Wen ZN,Walker SJ,et al. The MicroArray Quality Control (MAQC)-II study of common practices for the development and validation of microarray-based predictive models.Nature BiotechnologyMonth: 10 Year: 2010S5S16
10. Shi LM,Reid LH,Jones WD,Shippy R,Warrington JA,Baker SC,et al. The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements.Nature BiotechnologyMonth: 9 Year: 200624911511161
11. Sequeira PA,Martin MV,Vawter MP. The first decade and beyond of transcriptional profiling in schizophrenia.Neurobiol DisMonth: 1 Year: 2011451233621396449
12. Horvath S,Langfelder P. WGCNA: an R package for weighted correlation network analysis.BMC BioinformaticsDay: 29 Month: 12 Year: 20089559
13. Zhang B,Horvath S. A general framework for weighted gene co-expression network analysis.Stat Appl Genet Mol BiolYear: 20054 Article17.
14. Cai CC,Langfelder P,Fuller TF,Oldham MC,Luo R,van den Berg LH,et al. Is human blood a good surrogate for brain tissue in transcriptional studies?BMC GenomicsDay: 20 Month: 10 Year: 20101120053297
15. Fuller TF,Ghazalpour A,Aten JE,Drake TA,Lusis AJ,Horvath S. Weighted gene coexpression network analysis strategies applied to mouse weight.Mammalian GenomeMonth: 7 Year: 2007186-746347217668265
16. Kang HJ,Kawasawa YI,Cheng F,Zhu Y,Xu XM,Li MF,et al. Spatio-temporal transcriptome of the human brain.NatureDay: 27 Month: 10 Year: 2011478737048348922031440
17. Oldham MC,Horvath S,Geschwind DH. Conservation and evolution of gene co-expression networks in human and chimpanzee brains.Proceedings of the National Academy of Sciences of the United States of AmericaDay: 21 Month: 11 Year: 200610347179731797817101986
18. Oldham MC,Konopka G,Iwamoto K,Langfelder P,Kato T,Horvath S,et al. Functional organization of the transcriptome in human brain.Nature NeuroscienceMonth: 11 Year: 2008111112711282
19. Torkamani A,Dean B,Schork NJ,Thomas EA. Coexpression network analysis of neural tissue reveals perturbations in developmental processes in schizophrenia.Genome ResearchMonth: 4 Year: 201020440341220197298
20. Geschwind DH,Voineagu I,Wang XC,Johnston P,Lowe JK,Tian Y,et al. Transcriptomic analysis of autistic brain reveals convergent molecular pathology.NatureDay: 16 Month: 6 Year: 2011474735138038421614001
21. Ripke S,Sanders AR,Kendler KS,Levinson DF,Sklar P,Holmans PA,et al. Genome-wide association study identifies five new schizophrenia loci.Nature GeneticsMonth: 10 Year: 2011431096997621926974
22. 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.NatureDay: 6 Month: 8 Year: 2009460725674875219571811
23. Shi JX,Levinson DF,Duan JB,Sanders AR,Zheng YL,Pe'er I,et al. Common variants on chromosome 6p22.1 are associated with schizophrenia.NatureDay: 6 Month: 8 Year: 2009460725675375719571809
24. Goes FS,Zandi PP,Miao K,McMahon FJ,Steele J,Willour VL. Mood-incongruent psychotic features in bipolar disorder: Familial aggregation and suggestive linkage to 2p11-q14 and 13q21-33.American Journal of PsychiatryMonth: 2 Year: 2007164223624717267786
25. Torrey EF,Webster M,Knable M,Johnston N,Yolken RH. The Stanley Foundation brain collection and Neuropathology Consortium.Schizophrenia ResearchDay: 3 Month: 8 Year: 200044215115510913747
26. Narayan S,Tang B,Head SR,Gilmartin TJ,Sutcliffe JG,Dean B,et al. Molecular profiles of schizophrenia in the CNS at different stages of illness.Brain ResearchDay: 6 Month: 11 Year: 2008123923524818778695
27. Chen C,Grennan K,Badner J,Zhang DD,Gershon E,Jin L,et al. Removing Batch Effects in Analysis of Expression Microarray Data: An Evaluation of Six Batch Adjustment Methods.PLoS OneDay: 28 Month: 2 Year: 201162
28. Johnson WE,Li C,Rabinovic A. Adjusting batch effects in microarray expression data using empirical Bayes methods.BiostatisticsMonth: 1 Year: 20078111812716632515
29. Langfelder P,Zhang B,Horvath S. Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R.BioinformaticsDay: 1 Month: 3 Year: 200824571972018024473
30. Langfelder P,Horvath S. Eigengene networks for studying the relationships between co-expression modules.BMC Systems BiologyDay: 21 Month: 11 Year: 2007117408505
31. Li Y,Willer CJ,Ding J,Scheet P,Abecasis GR. MaCH: Using Sequence and Genotype Data to Estimate Haplotypes and Unobserved Genotypes.Genetic EpidemiologyMonth: 12 Year: 201034881683421058334
32. Wang J,Zhang KL,Cui SJ,Chang SH,Zhang LY. i-GSEA4GWAS: a web server for identification of pathways/gene sets associated with traits by applying an improved gene set enrichment analysis to genome-wide association study.Nucleic Acids ResearchMonth: 7 Year: 201038W90W9520435672
33. Langfelder P,Luo R,Oldham MC,Horvath S. Is My Network Module Preserved and Reproducible?PLoS Computational BiologyMonth: 1 Year: 201171
34. Thomas EA,Torkamani A,Dean B,Schork NJ. Coexpression network analysis of neural tissue reveals perturbations in developmental processes in schizophrenia.Genome ResearchMonth: 4 Year: 201020440341220197298
35. Saris CGJ,Horvath S,van Vught PWJ,van Es MA,Blauw HM,Fuller TF,et al. Weighted gene co-expression network analysis of the peripheral blood from Amyotrophic Lateral Sclerosis patients.BMC GenomicsDay: 27 Month: 8 Year: 20091019128516
36. Alexson TO,Hitoshi S,Coles BL,Bernstein A,van der Kooy D. Notch signaling is required to maintain all neural stem cell populations - Irrespective of spatial or temporal niche.Developmental NeuroscienceYear: 2006281-2344816508302
37. Breunig JJ,Silbereis J,Vaccarino FM,Sestan N,Rakic P. Notch regulates cell fate and dendrite morphology of newborn neurons in the postnatal dentate gyrus.Proceedings of the National Academy of Sciences of the United States of AmericaDay: 18 Month: 12 Year: 200710451205582056318077357
38. Lathia JD,Mattson MP,Cheng A. Notch: from neural development to neurological disorders.Journal of NeurochemistryMonth: 12 Year: 200810761471148119094054
39. Huang DW,Sherman BT,Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources.Nature ProtocolsYear: 2009414457
40. Lips ES,Cornelisse LN,Toonen RF,Min JL,Hultman CM,Consortium tIS,et al. Functional gene group analysis identifies synaptic gene groups as risk factor for schizophrenia.Mol PsychiatryYear: 2011
41. West AK,Hidalgo J,Eddins D,Levin ED,Aschner M. Metallothionein in the central nervous system: Roles in protection, regeneration and cognition.NeurotoxicologyMonth: 5 Year: 200829348950318313142
42. Choi KH,Elashoff M,Higgs BW,Song J,Kim S,Sabunciyan S,et al. Putative psychosis genes in the prefrontal cortex: combined analysis of gene expression microarrays.BMC PsychiatryDay: 7 Month: 11 Year: 2008818271956
43. Peilin Jia Jian Tian,Zhao Zhongming. Assessing gene length biases in gene set analysis of Genome-Wide Association Studies.International Journal of Computational Biology and Drug DesignYear: 20103429731021297229
44. Weinberger DR. Implications of Normal Brain-Development for the Pathogenesis of Schizophrenia.Archives of General PsychiatryMonth: 7 Year: 19874476606693606332
45. Pantelis C,Yucel M,Wood SJ,Velakoulis D,Sun DQ,Berger G,et al. Structural brain imaging evidence for multiple pathological processes at different stages of brain development in schizophrenia.Schizophrenia BulletinMonth: 7 Year: 200531367269616020551
46. Myers RA,Casals F,Gauthier J,Hamdan FF,Keebler J,Boyko AR,et al. A Population Genetic Approach to Mapping Neurological Disorder Genes Using Deep Resequencing.PLoS GeneticsMonth: 2 Year: 201172
47. Shao L,Vawter MP. Shared gene expression alterations in schizophrenia and bipolar disorder.Biological PsychiatryDay: 15 Month: 7 Year: 2008642899718191109
48. Craddock N O'Donovan MC,Owen MJ. Genes for schizophrenia and bipolar disorder? Implications for psychiatric nosology.Schizophrenia BulletinMonth: 1 Year: 200632191616319375
49. Storey JD,Tibshirani R. Statistical significance for genome-wide studies.Proc Natl Acad Sci USADay: 1 Month: 4 Year: 2003100738899412660367
50. Jia P,Wang L,Fanous AH,Chen X,Kendler KS,International Schizophrenia Consortiumet al. A bias-reducing pathway enrichment analysis of genome-wide association data confirmed association of the MHC region with schizophrenia.J Med GenetMonth: 2 Year: 2012492961022187495
51. Jia P,Wang L,Fanous AH,Pato CN,Edwards TL,The Internatio nal Schizophrenia Consortiumet al. Network-Assisted Investigation of Combined Causal Signals from Genome-Wide Association Studies in Schizophrenia.PLoS Computational BiologyYear: 201287e100258722792057

Article Categories:
  • Article

Keywords: Gene Expression, Schizophrenia, WGCNA, Neuron differentiation, Neuron protection.

Previous Document:  Effects of the mu-opioid receptor antagonist GSK1521498 on hedonic and consummatory eating behaviour...
Next Document:  CRHR1 genotypes, neural circuits and the diathesis for anxiety and depression.