A comprehensive network and pathway analysis of candidate genes in major depressive disorder
- Peilin Jia†1, 2,
- Chung-Feng Kao†3,
- Po-Hsiu Kuo3, 4Email author and
- Zhongming Zhao1, 2, 5Email author
© Jia et al. 2011
Published: 23 December 2011
Numerous genetic and genomic datasets related to complex diseases have been made available during the last decade. It is now a great challenge to assess such heterogeneous datasets to prioritize disease genes and perform follow up functional analysis and validation. Among complex disease studies, psychiatric disorders such as major depressive disorder (MDD) are especially in need of robust integrative analysis because these diseases are more complex than others, with weak genetic factors at various levels, including genetic markers, transcription (gene expression), epigenetics (methylation), protein, pathways and networks.
In this study, we proposed a comprehensive analysis framework at the systems level and demonstrated it in MDD using a set of candidate genes that have recently been prioritized based on multiple lines of evidence including association, linkage, gene expression (both human and animal studies), regulatory pathway, and literature search. In the network analysis, we explored the topological characteristics of these genes in the context of the human interactome and compared them with two other complex diseases. The network topological features indicated that MDD is similar to schizophrenia compared to cancer. In the functional analysis, we performed the gene set enrichment analysis for both Gene Ontology categories and canonical pathways. Moreover, we proposed a unique pathway crosstalk approach to examine the dynamic interactions among biological pathways. Our pathway enrichment and crosstalk analyses revealed two unique pathway interaction modules that were significantly enriched with MDD genes. These two modules are neuro-transmission and immune system related, supporting the neuropathology hypothesis of MDD. Finally, we constructed a MDD-specific subnetwork, which recruited novel candidate genes with association signals from a major MDD GWAS dataset.
This study is the first systematic network and pathway analysis of candidate genes in MDD, providing abundant important information about gene interaction and regulation in a major psychiatric disease. The results suggest potential functional components underlying the molecular mechanisms of MDD and, thus, facilitate generation of novel hypotheses in this disease. The systems biology based strategy in this study can be applied to many other complex diseases.
During the past decade, rapid advances in high throughput technologies have helped investigators generate numerous genetic and genomic datasets, aiming to uncover disease causal genes and their actions in complex diseases. These datasets are often heterogeneous and multi-dimensional; thus, it is difficult to find consistent genetic signals for the connection to the corresponding disease. Specifically in psychiatric genetics, there have been numerous datasets from different platforms or sources such as association studies, including genome-wide association studies (GWAS), genome-wide linkage scans, microarray gene expression, and copy number variation, among others. Analyses of these datasets have led to many exciting discoveries, including disease susceptibility genes or loci, providing important insights into the underlying molecular mechanisms of the diseases. However, the results based on single domain data analysis are often inconsistent, with a very low replication rate in psychiatric disorders [1, 2]. It has now been commonly accepted that psychiatric disorders, such as schizophrenia and major depressive disorder (MDD), have been caused by many genes, each of which has a weak or moderate risk to the disease [3, 4]. Thus, a convergent analysis of multi-dimensional datasets to prioritize disease candidate genes is urgently needed. Such an approach may overcome the limitation of each single data type and provide a systematic view of the evidence at the genomic, transcriptomic, proteomic, metabolomic, and regulatory levels [5, 6].
Recently, pathway and network-assisted analyses of genomic and transcriptomic datasets have been emerging as powerful approaches to analyze disease genes and their biological implications [7–11]. According to the observation of "guilt by association", genes with similar functions have been demonstrated to interact with each other more closely in the protein-protein interaction (PPI) networks than those functionally unrelated genes . Similarly, we have seen accumulating evidence that complex diseases are caused by functional related genes (e.g., in pathways or protein complex) through their dynamic interaction and regulation rather than action by single gene alone. Taken together, a systematic analysis and comparison of disease genes in the PPI network would provide additional insights into the diseases that otherwise could not be identified by single gene or single marker analysis. It is important to note that, although network-based analysis has been widely applied in major complex diseases such as cancer, its application in psychiatric diseases has been limited so far.
MDD is a complex mental disorder with a lifetime prevalence of 9-19% [13–15] and moderate heritability (37-43%) . Previous studies have suggested the involvement of polygenic and mutifactorial features in the pathology of MDD, as well as complex interactions among genes (G×G) and environmental factors (G×E) [17, 18]. Recently, we have performed the first gene prioritization using multi-dimensional evidence-based datasets in MDD, including association, linkage, gene expression (both human and animal studies), regulatory pathway, and literature search (both human and animal studies) . A list of depression candidate genes (which we named DEPgenes) with high reliability has been generated based on this strategy . However, several characteristics remain unclear: the functional relationships among these DEPgenes, how they interact and regulate with each other, and how they act in the MDD. Such investigations are warranted for a deeper understanding of the molecular mechanisms of MDD but require comprehensive analysis at the systems biology level.
In this study, we first explored DEPgenes in the context of the PPI network for their topological characteristics and compared them with two representative complex diseases: schizophrenia and cancer. We performed the functional enrichment analyses using annotations from both Gene Ontology (GO)  and canonical pathways. More importantly, we examined crosstalk among the significantly enriched pathways by quantitatively measuring the shared protein components between each pair of pathways. Finally, we constructed a MDD-specific subnetwork using the DEPgenes and validated them using the association data from an independent GWAS dataset for MDD. Our work demonstrated a practical framework for complex disease candidate gene analysis at the functional level, which can be applied to other complex diseases.
Materials and methods
Depression candidate genes
We modified the scoring scheme in the gene prioritization system proposed by Kao et al  and reprioritized a list of 151 DEPgenes for MDD using the updated data information. Briefly, several lines of evidence-based datasets were collected for MDD, including association studies, linkage scans, gene expression (both human and animal studies), literature search (both human and animal studies), and biological regulatory pathways. A dataset-specific score was assigned for each gene in each data source, and all data types were combined by an optimized weighting matrix to indicate the priority of a gene's association with MDD. The final gene list was selected based on a set of previously implicated core genes for MDD and validated by the GWAS dataset. Detailed information of this gene prioritization procedure can be found in Kao et al . Of note, the number of genes we used here is slightly different from that in Kao et al  due to the data and annotation updates, but the two lists were very similar.
Other data sources and process
For the purpose of comparison, we collected schizophrenia candidate genes and cancer genes. Schizophrenia is a severe psychiatric disorder and has been suggested to share certain comorbidity with MDD clinically and genetically . We included this disorder here to represent other psychiatric disorders for the purpose of comparison. We retrieved 160 schizophrenia candidate genes prioritized in our recent work using a similar multi-dimensional evidence-based strategy . Cancer has been the most studied among all complex disease and is expected to have substantially different pathological features from MDD. Thus, it would be interesting to see how those genes act differently at the network and pathway levels. Cancer genes were downloaded from the Cancer Gene Census database  (CGC, July 2011).
The human PPI data was downloaded from the Protein Interaction Network Analysis (PINA) platform (downloaded in March 2010) , which collected and annotated data from six public PPI databases (MINT, IntAct, DIP, BioGRID, HPRD, and MIPS/MPact). Only proteins that could be successfully mapped to NCBI protein-coding genes were included in our analysis (see below). After removing self-interaction and duplicates, the final network included a total of 10,377 nodes and 50,109 interactions.
The GWAS dataset for major depression (dbGaP Study Accession: phs000020.v2.p1) was retrieved through our approved access to dbGaP . We developed a pipeline for quality controls (QC) to the dataset. Detailed information can be found in our previous studies [19, 26–28]. As a brief summary, there were 1,738 depression patients and 1,802 matched normal controls, and 424,861 markers after QC, covering a total of 16,758 genes. This dataset was used to evaluate the genes identified in this work.
To coordinate these heterozygous datasets in this study, we downloaded several key annotation files from the National Center for Biotechnology Information (NCBI)  for the ease of integration. These included the annotation files of Homo_sapiens.gene_info.zip, gene_refseq_uniprotkb_collab.zip, and gene2refseq.zip (as of November 24, 2010). DEPgenes, schizophrenia candidate genes, cancer genes, PPI data, and GWAS data were all mapped to human protein-coding genes from NCBI. Those genes that could not be mapped appropriately were discarded from the subsequent analysis.
Network topological properties
In network analysis, there are several key topological indicators that have been defined to describe the behaviors or characteristics of the nodes in a network. The most representative ones are degree, betweenness, and shortest path. Degree is defined as the number of adjacent edges of a given node (protein) or the number of neighbor nodes interacting with it. Betweenness of a node is defined as the number of shortest paths going through the node; shortest path measures the nearest distance traveling from one node to another. We chose to examine the distribution of degree and betweenness of DEPgenes for exploration of their topological behaviors, and compared them with those of schizophrenia candidate genes  and cancer genes .
Functional enrichment tests
To perform functional enrichment tests of the candidate genes, we used WebGestalt  for Gene Ontology (GO) term analysis and used the Ingenuity Pathway Analysis (IPA) system  for both canonical pathways and molecular networks. Although WebGestalt can perform enrichment tests for the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways , the IPA system provides a more comprehensive pathway resource based on manual collection and curation. The rich information returned by IPA is also suitable for pathway crosstalk analysis (see below), as it has more molecules and their connections included. Briefly, WebGestalt implements the hypergeometric test for the enrichment of GO terms in the candidate genes, followed by the correction of multiple testing using the Benjamini & Hochberg (BH) method . The IPA system implements Fisher's exact test to determine whether a canonical pathway is enriched with genes of interest. Furthermore, the network analysis in the IPA system searches for significant molecular networks in a commercial knowledge base, including integrative information from literature, gene expression, and gene annotation.
We performed pathway crosstalk analysis using the pathways that were significantly enriched with DEPgenes after multiple testing correction. Two pathways are considered to crosstalk if they share a proportion of DEPgenes. We introduced two measurements to computationally indicate the overlap of a pair of pathways: the Jaccard Coefficient and the Overlap Coefficient , where A and B denote the number of candidate genes in the two pathways, respectively. To avoid non-specific inclusion of crosstalk, we further implemented the following rules: (1) only pathways with at least 5 DEPgenes were used; (2) only pathways with adjusted P values < 0.01 were used; and (3) two pathways in crosstalk were required to share at least 3 DEPgenes. These criteria were introduced to ensure that each of the pathways, as well as its crosstalk pair, have not only statistical significance but also a biologically meaningful number of genes, as some pathways may be too small. Finally, we found many significant pathways were identified by IPA; thus, they generated thousands of crosstalk events when all the pathway combinations were compared. In practice, we chose only those crosstalk events that had scores within the top 10% of the score distribution. Although these criteria were arbitrary, we found it worked efficiently to balance an appropriate number of pathways and crosstalk events.
Construction of MDD-specific subnetwork
To construct a MDD-specific subnetwork, we applied the Steiner minimum tree algorithm that is implemented in our software framework GenRev  to the 151 DEPgenes. Solving the Steiner minimum tree algorithm was proposed by Klein and Ravi , which can be used for constructing a connected subnetwork given a list of query nodes. In our case, the query nodes are those encoded by DEPgenes, and the whole network is the human interactome extracted from the PINA database (see above). This algorithm aims to connect a maximum proportion of the query nodes. To accomplish this, additional nodes in the network, but not in the query list, would be recruited in order to make the target subnetwork interconnected, while the algorithm is optimized towards a minimum list of the additional nodes. GenRev is a recently developed software tool which implements the Steiner minimum tree algorithm, as well as two other popular algorithms for subnetwork construction. It has been successfully applied in our previous work [6, 22, 37]. In the work discussed here, we used it for DEPgenes to construct MDD-specific subnetwork.
Network topological properties of depression genes
For the measurement of betweenness, the average value was 5.02 × 104 for DEPgenes, 4.01 × 104 for the schizophrenia genes, and 5.61 × 104 for cancer genes, while their median values were 5.12 × 103, 3.54 × 103, and 1.02 × 104, respectively. Similar to the measurement of degree, there was no significant difference in the betweenness values between the MDD and schizophrenia candidate genes (P = 0.21), but cancer genes had significantly larger betweenness values than DEPgenes (P = 0.03). These results indicated that the candidate genes for the two major psychiatric disorders, MDD and schizophrenia, shared similar topological features in the human interactome, while both had substantially different features when compared to cancer genes.
Gene Ontology enrichment analysis by WebGestalt
Gene Ontology (GO) terms enriched with module genes (GO level ≥ 4)
P BH $
GO:0007268: synaptic transmission
7.18 × 10-38
1.18 × 10-34
GO:0007267: cell-cell signaling
8.78 × 10-37
7.21 × 10-34
GO:0019226: transmission of nerve impulse
2.84 × 10-36
1.55 × 10-33
GO:0044057: regulation of system process
1.23 × 10-29
2.89 × 10-27
GO:0051239: regulation of multicellular organismal process
1.07 × 10-29
2.89 × 10-27
GO:0050877: neurological system process
4.72 × 10-26
9.69 × 10-24
GO:0007154: cell communication
6.14 × 10-26
1.12 × 10-23
GO:0008227: G-protein coupled amine receptor activity
5.13 × 10-21
7.18 × 10-19
GO:0035240: dopamine binding
2.25 × 10-13
1.58 × 10-11
GO:0005230: extracellular ligand-gated ion channel activity
4.17 × 10-12
2.34 × 10-10
GO:0004888: transmembrane receptor activity
6.55 × 10-12
2.78 × 10-10
GO:0005102: receptor binding
9.87 × 10-11
3.07 × 10-9
GO:0022834: ligand-gated channel activity
2.46 × 10-10
5.74 × 10-9
GO:0043005: neuron projection
2.43 × 10-22
3.91 × 10-20
GO:0044459: plasma membrane part
3.12 × 10-20
2.51 × 10-18
GO:0000267: cell fraction
1.45 × 10-19
7.78 × 10-18
GO:0005887: integral to plasma membrane
5.84 × 10-19
2.35 × 10-17
GO:0031226: intrinsic to plasma membrane
1.27 × 10-18
4.09 × 10-17
GO:0042995: cell projection
5.50 × 10-18
1.48 × 10-16
GO:0005886: plasma membrane
1.64 × 10-17
3.77 × 10-16
1.15 × 10-15
2.31 × 10-14
2.76 × 10-14
4.94 × 10-13
GO:0005626: insoluble fraction
7.11 × 10-13
1.14 × 10-11
Pathway enrichment by Ingenuity Pathway Analysis
Canonical pathways enriched with module genes by Ingenuity Pathway Analysis (IPA) (PBH < 10-6)
Ingenuity canonical pathways
P BH $
6.31 × 10-16
G-protein coupled receptor signaling
5.01 × 10-15
Serotonin receptor signaling
6.31 × 10-15
Corticotropin releasing hormone signaling
7.94 × 10-11
Dopamine receptor signaling
3.80 × 10-9
Glucocorticoid receptor signaling
1.35 × 10-8
2.29 × 10-8
Amyotrophic lateral sclerosis signaling
5.89 × 10-8
Neuropathic pain signaling on dorsal horn neurons
5.89 × 10-8
1.02 × 10-7
CREB signaling in neurons
1.62 × 10-7
6.61 × 10-7
Crosstalk among significantly enriched pathways
Since many genes and pathways might be involved in MDD, to more deeply understand how these pathways are related, we performed a pathway crosstalk analysis. We first selected the significantly enriched pathways from the IPA results. Specifically, we selected those pathways having PBH < 0.01 and ≥ 5 DEPgenes. There were 71 pathways that met these criteria. Among them, 69 pathways shared at least 3 genes with other pathways. A total of 571 edges (links) connected between any two of these pathways, and these edges were ranked according to the average scores of the Jaccard Coefficient and the Overlap Coefficient (see the Materials and methods section). We selected the top 10% edges, which resulted in 57 pairs of pathway crosstalk, and constructed the pathway crosstalk network for MDD. This pathway crosstalk was the first of its kind in MDD.
The second large pathway crosstalk module mainly consisted of immune-related pathways, such as IL-6 signaling and LXR/RXR activation (Figure 2). This strongly supports recent discoveries of immunity and inflammation related processes in psychiatric disorders [42, 43], including MDD . Many genes that drove the crosstalk in the figure were also found to function in both neuro- and immune-related processes like APOE , TNF , and IL6 .
Although there have been numerous reports of susceptibility genes or loci to psychiatric disorders such as major depressive disorder and schizophrenia, no disease causal genes have been confirmed [48–50]. One important task now is to reduce the data noise and prioritize the candidate genes from multiple dimensional genetic and genomic datasets that have been made available during the last decade and then explore their functional relationships for further validation. To our knowledge, this is the first systematic network and pathway analysis for MDD using candidate genes prioritized from comprehensive evidence-based data sources. By overlaying the MDD candidate genes in the context of the human interactome, we examined the topological characteristics of these genes by comparing them with those of schizophrenia and cancer candidate genes. We further performed pathway enrichment analysis to better understand the biological implications of these genes in the context of the regulatory system. Building on our observation of the large number of pathways enriched with DEPgenes, we developed novel approaches to measure pathway crosstalk so that complex gene action and regulation could be explored, thus providing us new insights into the interpretation of the underlying molecular mechanisms in MDD.
Our network topological analysis revealed that DEPgenes showed similar topological characteristics to schizophrenia, supporting previous reports that depression and schizophrenia might share comorbidity both clinically and genetically . For example, clinical symptoms such as psychosis and neuro-cognitive impairments have been observed in both depression and schizophrenia patients , and shared genetic variance has been reported between major depression and schizophrenia [51, 52]. Although similar network topological features are expected by many investigators, our study was the first to confirm, and provided further evidence, that the topological features of depression genes are different from cancer genes. It is worth noting that, although depression and schizophrenia genes had similar degree distributions (Figure 1), depression genes had moderately stronger connectivity and betweenness than schizophrenia genes.
Of significance, our pathway crosstalk analysis revealed two large clustered modules, both of which had important implications to MDD (Figure 2). The first cluster included 17 pathways, and it was dominated by neuro-signaling pathways. Among these pathways, neuropathic pain signaling in dorsal horn neurons (PBH = 5.89 × 10-8), CREB signaling in neurons (PBH = 1.62 × 10-7), synaptic long term potentiation (PBH = 6.17 × 10-5), and axonal guidance signaling (PBH = 1.55 × 10-4) are involved in neuron/brain tissues and have been reported to be involved in MDD [53, 54]. Our further examination of the genes contributing to the crosstalk revealed that the most frequently shared genes in this cluster were PRKACA (functioning in n = 15 pathways in this cluster), GNAS (n = 14), GNB3 (n = 13), ADCY7 (n = 10), GNAL (n = 9), AKT1 (n = 9), CREB1 (n = 8), CAMK2A (n = 6), GRIN2B (n = 5), GRIN2A (n = 5), and GRIN1 (n = 5), among others.
The second cluster is primarily related to immunity and inflammation, including the IL-6 signaling pathway (PBH = 6.17 × 10-3), differential regulation of cytokine production in macrophages and T helper cells by IL-17A and IL-17F (PBH = 8.13 × 10-6), and LXR/RXR activation (PBH = 4.57 × 10-4). For example, the LXR/RXR pathway may play a role in the prevention of programmed cell death and a role in immune responses to inhibit inflammatory gene expression . The most frequently shared genes in this cluster included TNF (functioning in n = 14 pathways), IL6 (n = 13), IL1B (n = 13), IL10 (n = 9), CCL2 (n = 8), NGFR (n = 7), and AKT1 (n = 7), among others. These genes further support the observation that immune- and inflammation-related functions are involved in this cluster. During recent years, evidence of immune and inflammation systems in psychiatric disorders has accumulated quickly [3, 4, 56].
In addition to the two major clusters, there are other crosstalk pairs that are noteworthy. The most interesting one is the pathway pair of cAMP-mediated signaling and G-protein coupled receptor signaling. The evidence linking these two pathways is strong, as its edge had a score 0.87. Moreover, these two pathways had the most significant enrichment test P values (6.31 × 10-16 and 5.01 × 10-15, respectively) in the IPA canonical pathway analysis (Table 2). The interaction between these two pathways involved 23 DEPgenes, including several serotonin receptor genes like HTR1A, HTR1B, and HTR5A. The cAMP-mediated signaling and G-protein coupled receptor signaling pathways have long been studied for their roles in the nervous system. Of note, there were several crosstalk links between one of these two pathways and other pathways that were enriched with the DEPgenes. Those pathway crosstalk connections were not shown in Figure 2 because they did not meet our stringent criteria for pathway inclusion (at least 3 DEPgenes shared between the pair of pathways or not within the top 10% crosstalk score, see the Materials and methods section). One example is the link between the cAMP-mediated signaling pathway and the serotonin receptor signaling, both of which were significantly enriched with DEPgenes, but their crosstalk score fell outside of the top 10% in the score distribution.
Our aim of the depression-specific subnetwork construction was to explore functional interactions of DEPgenes in a local protein-protein interaction environment. Our follow-up evaluation of the disease association of both DEPgenes and the additionally recruited genes using a major GWAS dataset for depression found that these genes tended to have small P-values (i.e., at the nominal significance level). Since the GWAS data we used here was an independent dataset, and GWAS was designed to be hypothesis free in genome-wide association studies, our survey of MDD-specific subnetwork genes demonstrated that this approach is efficient to find a set of genes that are both functionally interactive and enriched with the association signals of the corresponding disease. Therefore, this approach is not only promising to find novel disease candidate genes for future validation but also useful to study the disease at the systems biology level.
This work has a few limitations. First, our DEPgenes and the follow up pathway/network analyses were conducted based on computational strategies. Although informative, this approach generally requires extensive experimental validation. Thus, although we validated subnetwork genes at the genome-wide level using the GWAS dataset, further validation of specific novel genes using more samples is urgently needed. Second, the pathway crosstalk analysis was based on the scores measured by Jaccard Coefficient (JC) and Overlap Coefficient (OC). In this study, we selected the pathway pairs empirically, that is, those ranked in the top 10%. P values from a statistic test would be better applied to select significant crosstalk. We did not apply this method because the Ingenuity Pathway Analysis system is a commercial software tool, and the information needed to conduct such a statistic test is not publically available. Accordingly, we could only use the limited information for pathway crosstalk analysis. Third, the MDD-specific subnetwork was built on available human interactome data. Although the number and quality of protein interactions has recently improved greatly, the human interactome is still incomplete with many false positives . Additionally, subnetwork extraction relies on specific algorithms and corresponding parameters. Several algorithms exist for subnetwork extraction. In this study, we applied the Steiner minimum tree algorithm, which can effectively reduce unrelated nodes (genes) to be included, but it may also miss some important functional links. Our analysis, along with our recent application of this algorithm in other complex diseases (schizophrenia [5, 57], hepatocellular carcinoma , and epilepsy ), has demonstrated this strategy is practical and could provide valuable information of the interactions among DEPgenes.
We developed a systems biology framework for advanced and functional analyses of major depressive disorder candidate genes. The network topological analysis revealed similar network characteristics between depression and schizophrenia, but network characteristics of both depression and schizophrenia differed from cancer, consistent with previous clinical and genetic studies. However, the depression genes interacted moderately stronger than schizophrenia genes in the context of the protein-protein interaction network. Our pathway enrichment tests followed by pathway crosstalk analysis revealed that neurotransmission and immune systems might play key roles in the etiology of depression, assuming that our evidence-based DEPgenes were representative of depression. Notably, we found two major functional clusters in the pathway crosstalk network. We further constructed a depression-specific subnetwork, in which additional candidate genes were identified with enriched association signals using the depression GWAS dataset. These findings present a wealth of information for future validation. The framework we presented in this work can be applied to many other complex diseases, such as addiction and bipolar disorder.
This work was partially supported by NIH grants, 2009 NARSAD Maltz Investigator Award (to Z.Z.), 2010 NARSAD Young Investigator Award (to P. J.), and NSC and NHRI grants in Taiwan (NSC 99-2314-B-002-140-MY3, EX99-9918NC to P.H.K.). The GWA data sets were accessed through the Genetic Association Information Network (GAIN), database of Genotypes and Phenotypes (dbGaP) accession number phs000020.v2.p1 .
This article has been published as part of BMC Systems Biology Volume 5 Supplement 3, 2011: BIOCOMP 2010 - The 2010 International Conference on Bioinformatics & Computational Biology: Systems Biology. The full contents of the supplement are available online at http://www.biomedcentral.com/1752-0509/5?issue=S3.
- Kendler KS: Schizophrenia genetics and dysbindin: a corner turned? Am J Psychiatry 2004, 161: 1533-1536. 10.1176/appi.ajp.161.9.1533View ArticlePubMedGoogle Scholar
- Bosker FJ, Hartman CA, Nolte IM, Prins BP, Terpstra P, Posthuma D, van Veen T, Willemsen G, DeRijk RH, de Geus EJ, et al.: Poor replication of candidate genes for major depressive disorder using genome-wide association data. Mol Psychiatry 2011, 16: 516-532. 10.1038/mp.2010.38View ArticlePubMedGoogle Scholar
- Purcell SM, Wray NR, Stone JL, Visscher PM, O'Donovan MC, Sullivan PF, Sklar P: Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature 2009, 460: 748-752.PubMedGoogle Scholar
- Jia P, Wang L, Meltzer HY, Zhao Z: Common variants conferring risk of schizophrenia: A pathway analysis of GWAS data. Schizophr Res 2010, 12: 38-42.View ArticleGoogle Scholar
- Sun J, Jia P, Fanous AH, van den Oord E, Chen X, Riley BP, Amdur RL, Kendler KS, Zhao Z: Schizophrenia gene networks and pathways and their applications for novel candidate gene selection. PLoS ONE 2010, 5: e11351. 10.1371/journal.pone.0011351PubMed CentralView ArticlePubMedGoogle Scholar
- Jia P, Ewers JM, Zhao Z: Prioritization of epilepsy associated candidate genes by convergent analysis. PLoS One 2011, 6: e17162. 10.1371/journal.pone.0017162PubMed CentralView ArticlePubMedGoogle Scholar
- Baranzini SE, Galwey NW, Wang J, Khankhanian P, Lindberg R, Pelletier D, Wu W, Uitdehaag BM, Kappos L, Polman CH, et al.: Pathway and network-based analysis of genome-wide association studies in multiple sclerosis. Hum Mol Genet 2009, 18: 2078-2090. 10.1093/hmg/ddp120PubMed CentralView ArticlePubMedGoogle Scholar
- Wang K, Li M, Hakonarson H: Analysing biological pathways in genome-wide association studies. Nat Rev Genet 2010, 11: 843-854. 10.1038/nrg2884View ArticlePubMedGoogle Scholar
- Wang L, Jia P, Wolfinger RD, Chen X, Zhao Z: Pathway analysis of genome-wide association studies: methodological issues and perspectives. Genomics 2011, 98: 1-8.View ArticlePubMedGoogle Scholar
- Guo AY, Sun J, Jia P, Zhao Z: A novel microRNA and transcription factor mediated regulatory network in schizophrenia. BMC Syst Biol 2010, 4: 10. 10.1186/1752-0509-4-10PubMed CentralView ArticlePubMedGoogle Scholar
- Guo AY, Webb BT, Miles MF, Zimmerman MP, Kendler KS, Zhao Z: ERGR: an ethanol-related gene resource. Nucleic Acids Res 2009, 37: D840-845. 10.1093/nar/gkn816PubMed CentralView ArticlePubMedGoogle Scholar
- Song J, Singh M: How and when should interactome-derived clusters be used to predict functional modules and protein function? Bioinformatics 2009, 25: 3143-3150. 10.1093/bioinformatics/btp551PubMed CentralView ArticlePubMedGoogle Scholar
- Kessler RC, Berglund P, Demler O, Jin R, Merikangas KR, Walters EE: Lifetime prevalence and age-of-onset distributions of DSM-IV disorders in the national comorbidity survey replication. Arch Gen Psychiatry 2005, 62: 593-602. 10.1001/archpsyc.62.6.593View ArticlePubMedGoogle Scholar
- Bijl RV, A R, G vZ: Prevalence of psychiatric disorder in the general population: results of The Netherlands Mental Health Survey and Incidence Study (NEMESIS). Soc Psychiatry Psychiatr Epidemiol 1998, 33: 587-595. 10.1007/s001270050098View ArticlePubMedGoogle Scholar
- Vicente B, Kohn R, Rioseco P, Saldivia S, Levav I, Torres S: Lifetime and 12-month prevalence of DSM-III-R disorders in the Chile psychiatric prevalence study. Am J Psychiatry 2006, 163: 1632-1370.View ArticleGoogle Scholar
- Sullivan PF, Neale MC, Kendler KS: Genetic epidemiology of major depression: review and meta-analysis. Am J Psychiatry 2000, 157: 1552-1562. 10.1176/appi.ajp.157.10.1552View ArticlePubMedGoogle Scholar
- Farmer A, Breen G, Brewster S, Craddock N, Gill M, Korszun A, Maier W, Middleton L, Mors O, Owen M, et al.: The Depression Network (DeNT) Study: methodology and sociodemographic characteristics of the first 470 affected sibling pairs from a large multi-site linkage genetic study. BMC Psychiatry 2004, 4: 1-7. 10.1186/1471-244X-4-1View ArticleGoogle Scholar
- Belmaker RH, Agam G: Major depressive disorder. N Engl J Med 2008, 358: 55-68. 10.1056/NEJMra073096View ArticlePubMedGoogle Scholar
- Kao CF, Fang YS, Zhao Z, Kuo PH: Prioritization and evaluation of depression candidate genes by combining multidimensional data resources. PLoS One 2011, 6: e18696. 10.1371/journal.pone.0018696PubMed CentralView ArticlePubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000, 25: 25-29. 10.1038/75556PubMed CentralView ArticlePubMedGoogle Scholar
- Huang J, Perlis RH, Lee PH, Rush AJ, Fava M, Sachs GS, Lieberman J, Hamilton SP, Sullivan P, Sklar P, et al.: Cross-disorder genomewide analysis of schizophrenia, bipolar disorder, and depression. Am J Psychiatry 2010, 167: 1254-1263. 10.1176/appi.ajp.2010.09091335View ArticlePubMedGoogle Scholar
- Sun J, Jia P, Fanous AH, Webb BT, van den Oord EJ, Chen X, Bukszar J, Kendler KS, Zhao Z: A multi-dimensional evidence-based candidate gene prioritization approach for complex diseases-schizophrenia as a case. Bioinformatics 2009, 25: 2595-2602. 10.1093/bioinformatics/btp428PubMed CentralView ArticlePubMedGoogle Scholar
- The Cancer Gene Census database[http://www.sanger.ac.uk/genetics/CGP/Census/]
- Wu J, Vallenius T, Ovaska K, Westermarck J, Makela TP, Hautaniemi S: Integrated network analysis platform for protein-protein interactions. Nat Methods 2009, 6: 75-77. 10.1038/nmeth.1282View ArticlePubMedGoogle Scholar
- MDD GWAS dataset through dbGaP[http://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000020.v2.p1]
- Jia P, Wang L, Meltzer HY, Zhao Z: Pathway-based analysis of GWAS datasets: effective but caution required. Int J Neuropsychopharmacol 2011, 14: 567-572. 10.1017/S1461145710001446View ArticlePubMedGoogle Scholar
- Jia P, Zheng S, Long J, Zheng W, Zhao Z: dmGWAS: dense module searching for genome-wide association studies in protein-protein interaction networks. Bioinformatics 2011, 27: 95-102. 10.1093/bioinformatics/btq615PubMed CentralView ArticlePubMedGoogle Scholar
- Kao CF, Jia P, Zhao Z, Kuo PH: Enriched pathways for major depressive disorder identified from a genome-wide association study. Int J Neuropsychopharmacol 2011, in press.Google Scholar
- National Center for Biotechnology Information (NCBI)[http://www.ncbi.nlm.nih.gov/]
- Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, Rahman N, Stratton MR: A census of human cancer genes. Nat Rev Cancer 2004, 4: 177-183. 10.1038/nrc1299PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang B, Kirov S, Snoddy J: WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res 2005, 33: W741-W748. 10.1093/nar/gki475PubMed CentralView ArticlePubMedGoogle Scholar
- Ingenuity Pathway Analysis system[http://www.ingenuity.com/]
- The Kyoto Encyclopedia of Genes and Genomes (KEGG) Database[http://www.genome.jp/kegg/]
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Statist Soc Ser B 1995, 57: 289-300.Google Scholar
- Klein P, Ravi R: A nearly best-possible approximation algorithm for node-weighted Steiner trees. J Algorithms 1995, 19: 104-115. 10.1006/jagm.1995.1029View ArticleGoogle Scholar
- Zheng S, Tansey WP, Hiebert SW, Zhao Z: Integrative network analysis identifies key genes and pathways in the progression of hepatitis C virus induced hepatocellular carcinoma. BMC Med Genomics 2011, 4: 62. 10.1186/1755-8794-4-62PubMed CentralView ArticlePubMedGoogle Scholar
- Catena-Dell'Osso M, Bellantuono C, Consoli G, Baroni S, Rotella F, Marazziti D: Inflammatory and neurodegenerative pathways in depression: a new avenue for antidepressant development? Curr Med Chem 2011, 18: 245-255. 10.2174/092986711794088353View ArticlePubMedGoogle Scholar
- Opmeer EM, Kortekaas R, Aleman A: Depression and the role of genes involved in dopamine metabolism and signalling. Prog Neurobiol 2010, 92: 112-133. 10.1016/j.pneurobio.2010.06.003View ArticlePubMedGoogle Scholar
- Hettema JM, An SS, van den Oord EJ, Neale MC, Kendler KS, Chen X: Association study of CREB1 with Major Depressive Disorder and related phenotypes. Am J Med Genet B Neuropsychiatr Genet 2009, 150B: 1128-1132. 10.1002/ajmg.b.30935PubMed CentralView ArticlePubMedGoogle Scholar
- Juhasz G, Dunham JS, McKie S, Thomas E, Downey D, Chase D, Lloyd-Williams K, Toth ZG, Platt H, Mekli K, et al.: The CREB1-BDNF-NTRK2 pathway in depression: multiple gene-cognition-environment interactions. Biol Psychiatry 2011, 69: 762-771. 10.1016/j.biopsych.2010.11.019View ArticlePubMedGoogle Scholar
- Buehler MR: A proposed mechanism for autism: an aberrant neuroimmune response manifested as a psychiatric disorder. Med Hypotheses 2011, 76: 863-870. 10.1016/j.mehy.2011.02.038View ArticlePubMedGoogle Scholar
- Muller N, Myint AM, Schwarz MJ: The impact of neuroimmune dysregulation on neuroprotection and neurotoxicity in psychiatric disorders--relation to drug treatment. Dialogues Clin Neurosci 2009, 11: 319-332.PubMed CentralPubMedGoogle Scholar
- Christmas DM, Potokar J, Davies SJ: A biological pathway linking inflammation and depression: activation of indoleamine 2,3-dioxygenase. Neuropsychiatr Dis Treat 2011, 7: 431-439.PubMed CentralPubMedGoogle Scholar
- Laskowitz DT, Vitek MP: Apolipoprotein E and neurological disease: therapeutic potential and pharmacogenomic interactions. Pharmacogenomics 2007, 8: 959-969. 10.2217/14622422.214.171.1249View ArticlePubMedGoogle Scholar
- Gallowitsch-Puerta M, Pavlov VA: Neuro-immune interactions via the cholinergic anti-inflammatory pathway. Life Sci 2007, 80: 2325-2329. 10.1016/j.lfs.2007.01.002PubMed CentralView ArticlePubMedGoogle Scholar
- Papassotiropoulos A, Hock C, Nitsch RM: Genetics of interleukin 6: implications for Alzheimer's disease. Neurobiol Aging 2001, 22: 863-871. 10.1016/S0197-4580(01)00294-9View ArticlePubMedGoogle Scholar
- Allen NC, Bagade S, McQueen MB, Ioannidis JP, Kavvoura FK, Khoury MJ, Tanzi RE, Bertram L: Systematic meta-analyses and field synopsis of genetic association studies in schizophrenia: the SzGene database. Nat Genet 2008, 40: 827-834. 10.1038/ng.171View ArticlePubMedGoogle Scholar
- Jia P, Sun J, Guo A, Zhao Z: SZGR: a comprehensive schizophrenia gene resource. Mol Psychiatry 2010, 15: 453-462. 10.1038/mp.2009.93PubMed CentralView ArticlePubMedGoogle Scholar
- Piletz JE, Zhang X, Ranade R, Liu C: Database of genetic studies of bipolar disorder. Psychiatr Genet 2011, 21: 57-68. 10.1097/YPG.0b013e328341a346PubMed CentralView ArticlePubMedGoogle Scholar
- Bassett AS, Chow EW, AbdelMalik P, Gheorghiu M, Husted J, Weksberg R: The schizophrenia phenotype in 22q11 deletion syndrome. Am J Psychiatry 2003, 160: 1580-1586. 10.1176/appi.ajp.160.9.1580PubMed CentralView ArticlePubMedGoogle Scholar
- Bassett AS, Chow EW, Husted J, Weksberg R, Caluseriu O, Webb GD, Gatzoulis MA: Clinical features of 78 adults with 22q11 Deletion Syndrome. Am J Med Genet A 2005, 138: 307-313.PubMed CentralView ArticlePubMedGoogle Scholar
- Gass P, Riva MA: CREB, neurogenesis and depression. Bioessays 2007, 29: 957-961. 10.1002/bies.20658View ArticlePubMedGoogle Scholar
- Alme MN, Wibrand K, Dagestad G, Bramham CR: Chronic fluoxetine treatment induces brain region-specific upregulation of genes associated with BDNF-induced long-term potentiation. Neural Plast 2007, 2007: 26496.PubMed CentralPubMedGoogle Scholar
- Valledor AF: The innate immune response under the control of the LXR pathway. Immunobiology 2005, 210: 127-132. 10.1016/j.imbio.2005.05.007View ArticlePubMedGoogle Scholar
- Shi J, Levinson DF, Duan J, Sanders AR, Zheng Y, Pe'er I, Dudbridge F, Holmans PA, Whittemore AS, Mowry BJ, et al.: Common variants on chromosome 6p22.1 are associated with schizophrenia. Nature 2009, 460: 753-757.PubMed CentralPubMedGoogle Scholar
- Sun J, Wan C, Jia P, Fanous AH, Kendler KS, Riley BP, Zhao Z: Application of systems biology approach identifies and validates GRB2 as a risk gene for schizophrenia in the Irish Case Control Study of Schizophrenia (ICCSS) sample. Schizophr Res 2011, 125: 201-208. 10.1016/j.schres.2010.12.002PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.