Skip to main content

Analysis of schizophrenia and hepatocellular carcinoma genetic network with corresponding modularity and pathways: novel insights to the immune system



Schizophrenic patients show lower incidences of cancer, implicating schizophrenia may be a protective factor against cancer. To study the genetic correlation between the two diseases, a specific PPI network was constructed with candidate genes of both schizophrenia and hepatocellular carcinoma. The network, designated schizophrenia-hepatocellular carcinoma network (SHCN), was analysed and cliques were identified as potential functional modules or complexes. The findings were compared with information from pathway databases such as KEGG, Reactome, PID and ConsensusPathDB.


The functions of mediator genes from SHCN show immune system and cell cycle regulation have important roles in the eitology mechanism of schizophrenia. For example, the over-expressing schizophrenia candidate genes, SIRPB1, SYK and LCK, are responsible for signal transduction in cytokine production; immune responses involving IL-2 and TREM-1/DAP12 pathways are relevant for the etiology mechanism of schizophrenia. Novel treatments were proposed by searching the target genes of FDA approved drugs with genes in potential protein complexes and pathways. It was found that Vitamin A, retinoid acid and a few other immune response agents modulated by RARA and LCK genes may be potential treatments for both schizophrenia and hepatocellular carcinoma.


This is the first study showing specific mediator genes in the SHCN which may suppress tumors. We also show that the schizophrenic protein interactions and modulation with cancer implicates the importance of immune system for etiology of schizophrenia.


Recent studies suggest that schizophrenia may result from neuropathological abnormalities and imbalanced immune systems. Signal transduction dysfunction of the neuroendocrine system are responsible for schizophrenia, especially the dopamine, serotonin and glutamate system in the temporal and frontal lobe of the brain area [1, 2]. Although an increasing number of studies show that the immune-mediated mechanism for inflammation responses are the pathogenesis of schizophrenia [3], the corresponding specific complexes, pathways and candidate genes are not well-documented for the etiological model of schizophrenia.

In recent years, there have been many studies focusing on the discovery of schizophrenic candidate genes and the construction of PPI networks and related pathways for the hope of a better understanding of schizophrenia. However, genetic association researches have been published with largely inconsistent results [4]. It was generally believed that a protein sub-network, rather than a single gene or genetic variants, accounts for the susceptibility of schizophrenia. Sun J. et al. (2008) surveyed the increased association studies from the SchizophreniaGene database in ethnic populations [5], in which candidate genes are selected and ranked by the combined odds ratio method as an important index of the candidate genes [6]. It provides a basis for the investigation of molecular and cellular mechanisms of schizophrenia by the analysis of gene features for a genetic network. A regularly updated online database of genetic association studies for schizophrenia (SZGene) was collected from Allen NC. et al. (2008)[4]. Sun J. et al. (2010) [7] selected a list of schizophrenia candidate genes by a multi-dimensional evidence-based approach to provide a comprehensive review of the schizophrenia molecular networks. The identified pathway characteristics of schizophrenic candidate genes have important implications of molecular features for schizophrenia. Another gene risk prediction study used the translational convergent functional genomics approach introduced by Ayalew M. et al. (2012) to prioritize schizophrenia genes by gene-level integration of genome-wide association study data to identify top candidate genes [8]. These candidate gene studies conclude the specific genetic variants or patterns contributing to the schizophrenic model by integrating functional and genotypic data. The previous literatures provide different databases and integration of formulated reliability analysis, ranking and scoring for important candidate genes of schizophrenia.

Schizophrenic patients have less chance to develop cancer than the general population [9]. Lower incidence of cancers, especially in lung, prostate and bladder cancer, was found in schizophrenic patients [1012]. Research suggests that cancer risk decreases as the duration and age of onset of schizophrenia increases [13]. Cancer protective factors in schizophrenic patients are genetic predisposition [14, 15]. These literature reviews have implication of sharing common disease genes or pathways between schizophrenia and cancer, and that schizophrenia is a protective factor for cancer [16].

To demonstrate the genetic relationship between schizophrenia and cancer, network biology and systemic bioinformatics data such as protein-protein interactions (PPIs) and related pathways were introduced. The data of human PPIs brought insights to the network biology of diseases and explained the interrelationships among disease-related genes and proteins. Through the development of modulation interaction networks of schizophrenic candidate genes, the related resources of molecular biology were integrated to explore the molecular biological information of disease mechanism and related drug targets or complexes.

Efforts on the exploration of schizophrenic common pathways from corresponding candidate gene analysis are gaining more attention and represent for novel treatment approaches in schizophrenia. Postulated disease networks are analyzed by tools or algorithms such as modularity, centrality (closeness and degree) and clique analysis derived from network biology, which the functional relevance of different gene sets and related biological significance were analyzed. In functional genomics, there are available integrative protein interaction databases developed to identify gene sets of interest which involve similar disorders. These gene sets are commonly presented as gene modules, protein complexes or pathways such as in the Database for Annotation, Visualization and Integrated Discovery (DAVID) [17], Kyoto Encyclopedia of Genes and Genomes (KEGG) [18] and ConsensusPathDB [19]. In these integrative databases, candidate gene sets from disease-related network to gene ontology classification were mapped to the related molecular pathways and PPI networks.

This study integrates comparative analysis of different genetic research results. From the RNA extraction of microarray data, the expression level of each gene was acquired from BA22-derived brain cells and hepatocellular carcinoma cells. Generated from two group sets of candidate genes, the corresponding PPI networks were constituted and analyzed. The over- and under- expression level of genetic interactions between schizophrenia and hepatocellular carcinoma are not only found by the direct effect of inhibition of candidate genes for cancer, but also through an indirect modulation of protein-protein interactions in the cancer genetic network which have potential effects on tumor suppression by analysis of the core schizophrenia-cancer genetic network. The differences in gene expression and PPI sub-networks between schizophrenia and hepatocellular carcinoma were analyzed to discover protein complexes and possible drug targets.


By analyzing microarray data of Brodman Area 22 (BA22), susceptible genes for schizophrenia were proposed. Through the analysis of the PPI network of schizophrenia and cancer, the potential complexes or possible drugs were proposed. The research flow is shown in Figure 1. The related tools for candidate gene resources, PPI networks and pathway databases with analytic tools or algorithms are described in the following sections.

Figure 1
figure 1

Flow chart of research protocol. By analysis of BA22 microarray data and Roessler liver 2 sample, the significant candidate genes were screened and constructed the PPI network for schizophrenia and cancer. The potential complexes or drugs are discovered by their corresponding pathways. The right part shows the corresponding tools used for analysis of candidate gene resources, PPI networks and pathway databases, the potential pathways and drugs.

Schizophrenia related genetic information

Schizophrenic candidate genes from literature studies such as Sun J. et al. and Ayalew M. et al. were selected based on the chromosome classification, mapping genetic literature data and statistic measures, and the highly relevant genes were sorted with different ranking systems [7, 8]. Ayalew M. et al. analyzed the schizophrenic candidate genes by ranking and scoring the relevant candidate genes from NCBI literature [8].

Selection of schizophrenic candidate genes by microarray data

Human BA22 of the prefrontal cortex is believed to be responsible for many positive symptoms and cognitive dysfunction in patients with psychiatric illness. RNA was extracted from post-mortem BA22 tissue from schizophrenic and control patients. The RNA samples were analyzed by Affymetrix GeneChip HG-U133 Plus2.0. We downloaded the microarray data (GSE21935) from the NCBI GEO database [20]. This dataset consists of 19 control and 23 schizophrenia samples.

The over- and under-expression genes in the BA22 samples were selected using the Student's t-test between the schizophrenia and control samples. The genes of the corresponding probes with p-value < 0.05 were defined as abnormally expressed and proposed as the candidate genes for schizophrenia.

Cancer-related genes by microarray

The expression data for human cancer including breast cancer, leukemia, colon cancer and prostate cancer is collected in the ONCOMINE database ( This database currently contains 674 datasets and information of 73327 tissue samples (ONCOMINE version 4.4.3).

The Roessler liver 2 sample [21] includes 445 samples which contains 225 hepatocellular carcinoma and 220 liver samples, a total of 12624 mRNA expressions are measured by the Human Genome U133A 2.0 Array and the data is released on 2009/11/1 by ONCOMINE [22]. 126 over- and 126 under-expression genes are selected which respectively account for 1% of top candidate genes from the Roessler liver 2 sample of ONCOMINE. The over- and under- expression candidate genes from BA22 samples and the Roessler liver 2 samples are listed in Additional file 1.

Construction of schizophrenia and cancer network

In order to construct a PPI network, the fundamental basis of human PPI network was formulated by the integration of interaction databases including BIND [23], HPRD [24], MINT [25], BioGrid [26], DIP [27] and ConsensusPathDB [19]. ConsensusPathDB currently contains the most comprehensive publicly available repository including genes, proteins and complexes interaction for Homo sapiens.

In PPI networks, each node represents an encoded gene and each edge represents a protein interaction by literature reviews or experiments. The interaction network with significant functionalities generates genetic network through the selection of different query genes, such as Level-One PPI (L1PPI) and Query-Query PPI (QQPPI). QQPPI networks include only the query marker genes as the nodes and show direct interactions among these queries. L1PPI networks also show other non-query nodes directly connected to the queries. L1PPI network allowed analysis of an extended network and indicated indirect interactions [28].

The over- and under-expression candidate genes from schizophrenia and hepatocellular carcinoma were combined into a single gene list, defined as SHCGene. The SHCN is defined as the L1PPI and QQPPI networks using the SHCGene as the query genes. The QQPPI of SHCN in Figure 2 illustrates the direct relationship of protein-protein interaction between the two diseases. In order to find the mediator genes between diseases, the L1PPI of SHCN was analyzed to find the mediators of query genes. By the comparison of different mediators, the potential core modulation network for schizophrenia was effectively extracted.

Figure 2
figure 2

The over- and under-expression PPI network for schizophrenia and hepatocellular carcinoma. Red nodes denote over-expression genes; green nodes represent under-expression genes, circle nodes indicate hepatocellular carcinoma genes, triangular nodes indicate schizophrenia genes. The under-expression gene sets between schizophrenia and hepatocellular carcinoma genes are noted including CD4-CXCL12-DPP4, HLA-DPA1-CD4-FCGR3-PTPRC, MYO1C-ESR1-RARA, FOS-JUN-EGR1, etc.

Selection of candidate genes, complexes and modularity from SHCN

GenRev [29] was used to construct the functional modularity of schizophrenia by calculating the betweeness and closeness centrality of genes. The gene included by the QQPPI of SHCN were used as the input for GenRev to construct a reference network and define the importance of genes using sub-network analysis by calculating MCL [30]. Combined with the SHCGene, the candidate genes collectively represented highly modulating functionality or inhibitory relationship between the diseases [31]. Using the clique analysis network algorithm, the involving complexes were extracted from the PPI network.

Analysis of schizophrenic pathways and drugs

By analysis and verification of known pathways, image information resource ( provided by pathway central presents the pathway information related to schizophrenia. The current available pathway databases for analyzing modulation network mechanisms include KEGG [18], Reactome [32], PID [33] and ConsensuspathDB[19]. ConsensuspathDB enables analysis of different types of functional interactions between genes and regulatory pathways [34] through integration of meta-databases such as KEGG, Reactome and PID.

The interactors of SHCN were mapped to the respective pathways of ConsensusPathDB which provide corresponding information from original pathway databases. The corresponding detailed pathway information from each pathway database were searched for the relationship between pathways and SHCN. The significant pathways from the pathway enrichment analysis test involving genes from the L1PPI of SHCN. The corresponding genes from the pathway databases were then prioritized. The importance of the respective pathways was evaluated by p-values [35].

By integration of schizophrenic related genes and pathways, novel drugs were discovered for further investigation. The target genes for specific drugs or complexes for query gene functions and interactions were searched against DrugBank [36] which lists FDA-approved agents and target complexes.

Potential drugs or complexes related to tumor suppression

In order to analyze the potential complexes or targets which might have tumor suppression effects, the clique complexes were searched against CORUM [37] to find significant protein complexes. Tumor suppression genes were collected for comparison from TSGene [38].


Schizophrenia candidate genes related to tumor suppression

Despite the increased risk factors for the schizophrenic patients such as heavy smoking, poor diet habit and inadequate physical activities, the protective factor for cancer incidence in schizophrenic patients [16] such as TP53 and APC, which plays a key role in the susceptibility of schizophrenia and the reduced cancer risk by apoptosis [3942], have implication of the explanation for less incidence of cancer for schizophrenic patients. Therefore, to see the differential interaction between diseases in the PPI network, schizophrenic candidate genes were examined from three different literature results and compared with the BA22 over- and under-expression genes with tumor suppression. Of the 716 tumor suppression genes, 26 genes were contained in Sun J. et al.(2008), 38 genes were contained in Allen NC. et al. (2008) and 47 genes were contained in Ayalew M. et al.(2012), this implicates that schizophrenia candidate genes are related to the tumor suppression genes and share a common genetic biological regulation. A total of 7.6% tumor suppression genes also appeared in schizophrenia candidate gene list. ZFHX3, RND3, KLF5, ERBB4, EGR1 and APC genes are both schizophrenia candidate genes and tumor suppression genes (Additional file 2). From the genetic aspect, the genes which were contained in both the schizophrenia candidate gene list and tumor suppression gene list have implication of protective genetic interaction relationship between schizophrenia and cancer.

The overlapped genes in schizophrenic genetic network

The schizophrenic candidate genes reported by the literature reviews of Sun J. et al.(2008), Allen NC. et al.(2008), Ayalew M. et al.(2012) and BA22 microarray samples are inconsistent, which might be attributed to the ethnic difference in allele and haplotype frequencies [43]. Furthermore, when comparing the 781 genes of the SZGene database from Allen NC. et al. (2008) [4], the 183 candidate genes from Ayalew M. et al.(2012)[8] and the 75 candidate genes from Sun J. et al.(2008) [5], only 1 overlapped gene (BDNF) was found, and only 36 overlapped genes appeared in any two of three sets. However, when we extended the previous three gene lists using L1PPI analysis, we obtained 10679, 4109, and 1821 PPIs and 4386, 2611, and 1438 interactors for the 781, 183, and 75 genes, respectively. As a result, there are now 663 interactors co-existing in the three L1PPI extended gene lists, and 1520 interactors co-existing in any two of three L1PPI extended gene lists. The overlapping candidate genes in the three original sets range from 0.12% to 1.4%; however, the overlapped interactors range from 15% to 46%, which indicate the importance of mediator genes (Additional file 3) in the candidate genetic network of schizophrenia rather than query genes as key roles of disease susceptibility.

Analysis of the schizophrenic genetic network with different expression levels by human protein-protein interactions

By analyzing the microarrays of human BA22 samples, the over- and under-expression network for schizophrenia and hepatocellular carcinoma reflect the interaction of both diseases by annotation of each different node with over- and under-expression features. A total of 472 genes including 138 over-expression genes and 334 under-expression genes are selected from the BA22 samples with a p-value less than 0.05, derived 3247 interacts of direct interactors of L1PPI for schizophrenia.

Direct interaction genes observed from QQPPI interaction formulates potential common functional modularity between schizophrenia and cancer. The genes which are contained in both the schizophrenic candidate gene list and the hepatocellular carcinoma candidate gene list from the Roessler liver 2 sample constitutes 197 over- and under-expression level genes and 264 PPIs in the QQPPI of SHCN network (Figure 2).

Many under-expression gene sets between schizophrenia and hepatocellular carcinoma genes are noted, including CD4-CXCL12-DPP4, HLA-DPA1-CD4-FCGR3-PTPRC, MYO1C-ESR1-RARA, FOS-JUN-EGR1, etc. Moreover, the combination of over- and under-expression level genes of both diseases such as CASC5-ZWINT, MARCO-SCGB3A2, PLD1-PEA15, GTF3C3-POLR3C, TBCE-TBCD-XRCC6, etc.(pink region) highlight the potential biological significance of gene set combination implicating the protective factor for both diseases. In the schizophrenia-hepatocellular carcinoma network (SHCN), the numbers of under-expression genes exceed the over-expression genes for schizophrenia. Furthermore, the under-expression schizophrenia genes interact with the over-expression genes of hepatocellular carcinoma (Figure 2) which have the implication of a genetic modulation mechanism for both diseases. These gene sets might have important roles in potential cellular modulation or neurodevelopmental regulation functions of disease pathophysiological mechanism by their involvement in the molecular pathways of related complexes.

Modularity and complex analysis of SHCN

Direct interaction genes could be easily observed from the QQPPI of SHCN and proposed as significant common functional modularity from the potential modulation process of both diseases which may have a lot of hidden biological significance. GenRev analysis showed that the significant genes and functional modularity were discovered from the QQPPI of SHCN by the expression level of each gene input as the score of GenRev. SYK and LCK genes were calculated with highest betweeness and ranked as top genes for the SHCN. The discovery of functional modularity by the MCL algorithm is listed in Table 1. With the use of systemic and integrated analysis of functional annotation information by using the DAVID, novel enriched gene list with functional modularity was analyzed. The functional modularity from GenRev can be mapped to DAVID to annotate the potential functions involving T cell and lymphocyte activation, cell mediated cytotoxicity, DNA and calcium binding and inflammatory response. The over-expression schizophrenia genes involving the SYK and LCK genes play a key role in functional modularity in T cell, leukemia and lymphocyte activation which is crucial in immune-related responses.

Table 1 Modularity analysis from SHCN

Common mediator genes of both diseases could be easily observed from the L1PPI of SHCN. Clique analysis formulates common functional modules and enables easy screening of the co-expression or functional units of gene sets which maps to related target complexes or common pathways of biological function [28, 31]. As shown in Figure 2, three clique-4 gene sets were derived from SHCN: CD4-FCGR3A-LCK-PTPRC (sky blue region), ESR1-NMP1-PARP1-RARA (green region) and HOXB7-PARP1-PRKDC-XRCC6 (purple region). The corresponding complexes were retrieved from CORUM[37] and are shown in Table 2, which mediate genetic functions such as cell cycle, transcriptional activation and immune responses. These complexes involve potential modulating mechanism for the discovery of novel treatment agents.

Table 2 Clique-4 complexes in SHCN

Furthermore, the L1PPI of SHCN formulates the critical functional modules and mediators by clique analyses. The most extended sub-network is clique-5. Figure 3 presents the 34 genes with 118 PPIs in the clique-5 sub-network. Two gene groups are identified. One is the clique CD4-FCGR3A-ZAP70-LCK-PTPRC which formulates the functional module of cytokine production related to immune system response; the other is the clique CKS1B-CCNB1-CDK1-CKS2-UBC which is responsible for cell cycle regulation. The schizophrenia-related SNP study supports the significant result of the involving pathway of translocation of ZAP-70 to immunological synapse by Reactome [44].

Figure 3
figure 3

Clique-5s for potential complexes in SHCN. Blue star nodes indicate mediators; red nodes represent up-expression genes; green nodes represent under-expression genes; circle nodes indicate hepatocellular carcinoma genes; triangular nodes indicate schizophrenia genes: (a) The CD4-FCGR3A-ZAP70-LCK-PTPRC clique formulates the functional module of cytokine production related to the immune system response. (b) The CKS1B-CCNB1-CDK1-CKS2-UBC clique is responsible for cell cycle regulation.

UBC and TP53 are mediator genes which are potential targets involved in the disease mechanism for schizophrenia and cancer[4548] which appears in the clique-5 network. PRKDC, PARP1, NPM1 and XRCC6 are hepatocellular carcinoma over-expression genes which modulate cell cycle regulation network through the modulation of UBC[49]. HOXB7 and RARA are schizophrenic genes with different gene expressions, modulated through the retinoid signaling pathway by the hepatocellular carcinoma over-expression gene PARP1 [50]. Furthermore, these genes formulate an important genetic functionality of the immune system, which illustrate the relationship between schizophrenia and autoimmune diseases.

The immune-related pathway responsible for pathological mechanism of schizophrenia

In order to prioritize the potential pathways in which schizophrenia and hepatocellular carcinoma candidate genes are involved, the L1PPI extended over- and under-expression genes from BA22 sample and Roessler liver 2 sample were used to search for significant pathways in PID. The crucial pathways are listed in Additional file 4 which is ranked by the p-value and FDR-adjusted p-value by the Benjamini-Hochberg procedure. The top ranked pathways with significant p-value associated with BA22 and Roessler liver 2 sample include the PDGFR-beta signaling pathway, the ErbB1 downstream signaling and BCR signaling pathway which indicates common pathways for schizophrenia and cancer. It is appealing that another group of significant pathways including the Fc-epsilon receptor I signaling in mast cells, the TCR signaling in CD4+ T cells and IL2-mediated signaling events highlight the importance of immune system mediated pathways in the key role of schizophrenia susceptibility.

A specific gene set of schizophrenic over-expression genes includes SIRPB1, LCK and SYK (Figure 2). SYK and LCK were the over-expression genes with high betweenness, however, SIRPB1 was an over-expressed gene for schizophrenia and directly linked to SYK from the QQPPI of SHCN, and it may have an implication of a potential disease mechanism for schizophrenia. Another responsible pathway involves the IL-2 pathway which Interleukin-2 binds to the IL-2 receptor to activate LCK and SYK, which helps explain that the over-expression genes: SIRPB1, LCK and SYK might be responsible for one of the possible disease mechanisms for schizophrenia. The corresponding pathways containing SIRPB1, LCK and SYK were retrieved from KEGG, PID and Reactome. Strong correlations with immune system pathways were found which involve the IL-2 and TREM-1/DAP12 pathway which are responsible for the etiology mechanism for schizophrenia in PID (Figure 4). The SIRPB1 protein interacts with the glycoprotein DAP12 to induce signal transduction and formation of the SIRPB1/DAP12 complex, and then transduces signals to the nucleus with the activation of SYK. Activation of DAP12 and triggering of the TREM-1receptor result in the production of pro-inflammatory cytokines and receptor expression on natural killer (NK) cells, monocytes and neutrophils with degranulation of neutrophilic granules and phagocytosis.

Figure 4
figure 4

IL-2 and TREM-1/DAP12 pathway. A specific gene set of schizophrenic over-expression genes includes SIRPB1, LCK and SYK. Strong correlations with immune system pathways were found which involve the IL-2 and TREM-1/DAP12 pathway which are responsible for etiology mechanism for schizophrenia in PID. The SIRPB1 protein interacts with the glycoprotein DAP12 to induce the signal transduction and the formation of the SIRPB1/DAP12 complex, and then transduces signals to the nucleus with the activation of SYK. Activation of DAP12 and triggering of the TREM-1 receptor result in the production of pro-inflammatory cytokines including interleukin family cytokines, TNF-a, and GM-CSF, etc.

Another responsible pathway involves the IL-2 pathway in which Interleukin-2 binds to the IL-2 receptor to activate LCK and SYK, then induces cascade signal transduction through Rho, PI3K and Akt/PKB signaling pathways. The results help to explain that the over-expression genes SIRPB1, LCK and SYK are responsible for one of the possible disease mechanisms for schizophrenia. The immune-related pathways reveal the novel discovery for the treatment and pathophysiological mechanism for schizophrenia.

Discovery of candidate drugs or treatments for both schizophrenia and cancer

The FDA-approved potential candidate drugs mapped from L1PPI of SHCN by the clique analysis and prioritized from DrugBank are listed as Table 3. The relationship between candidate drugs and encoded genes are illustrated in Figure 5. Cetuximab, Immunoglobulin, Rituxan and Thymoglobulin are immune-related drugs encoded by FCGR3A and CD4 genes in which Cetuximab is an epidermal growth factor receptor (EGFR) inhibitor used for the treatment of metastatic colorectal cancer and head and neck cancer [51, 52]. The use of Rituxan(Rituximab) has been proposed for patients with B-cell non-Hodgkin's lymphoma (NHL) and B-cell chronic lymphocytic leukaemia (CLL) [53, 54]. It shows that the immune-related drugs with the potential of antibody therapy are associated with several mechanisms of potential treatment for cancer including interference of vital signaling pathways by tumor-bound antibody through the Fc portion of the antibody [55].

Table 3 Potential drugs discovered from Clique-5 network
Figure 5
figure 5

Drugs and related genes. Potential candidate drugs are mapped to the L1PPI of SHCN using clique analysis and prioritized using DrugBank data. (A) The drugs with green background are immune-related drugs. (B) The drugs with red background are anti-cancer drugs. (C) The drugs with yellow background are sex-hormone drugs. (D) The drugs with purple background are retinoid drugs.

The sex hormones, estrogen and thyroid hormone (T3), are candidate drugs involving THRA, THRB and ESR1 in the SHCN. Dasatinib and Indolocarbazole are protein kinase inhibitors involving the LCK and ZAP70 genes. Vorinostat and Flavopiridol are anti-cancer agents involving the HDAC1, HDAC2 and CDK1 genes. Lovastatin involving HDAC2 is a cholesterol-lowering agent. These candidate drugs mostly act as cancer treatment agents but have little evidence for treatment of schizophrenia.

Adapalene, Tazarotene, Tamibarotene are retinoids which involved RARA gene with multiple functions including eye vision, immune function, and activation of tumor suppressor genes. The retinoid acid has been reported for cancer prevention and treatment of leukemia. Promising results by Fenretinide (retinoid acid) in breast cancer prevention provides a strong rationale for cancer treatment especially in combination with chemotherapy in non-small cell lung cancer [56].

The retinoids with the promising role in chemoprevention of premalignant lesions in the head and neck have been the focus of cancer intervention treatment [57] which Stra6 unregulated RA-responsive genes unregulated by DNA damage with important role in cell death responses. Novel findings between the retinoid acid and TP53 pathways provide a new insight which enhances the tumor suppression functions. It implicates the significant role of vitamin A metabolites in cancer prevention and treatment [58]. Moreover, Vitamin A (retinol), the biologically active form of retinoic acid, has been proposed to be involved in the pathogenesis of schizophrenia by the genetic basis of encoding retinoid acid metabolism enzymes. Seven genes were investigated and involved in the synthesis, degradation and transportation of RA, ALDH1A1, ALDH1A2, ALDH1A3, CYP26A1, CYP26B1, CYP26C1 and Transthyretin (TTR), for their roles in the development of schizophrenia [59]. The expression of the transthyretin (TTR) tetramer, which is a retinoid transporter, is increased significantly in the plasma of schizophrenic patients. Retinoid dysfunction might be involved in the pathology of schizophrenia[60].


Few microarray studies of mental disorders have used post-mortem brain samples of human species from schizophrenic patients. Researchers did not have convenient access to brain samples of psychiatric patients until 1994 when the Stanley Brain Collection started. However, the result of using a single expression dataset could be biased. A more diverse microarray dataset of BA22, BA10, BA46 samples and other tissue specific samples from schizophrenic patients could be compared and analyzed for the future study.

In order to determine the extent of which the BA22 genetic network is the result of chance, we introduced 5000 randomly generated control networks by randomly selecting 472 genes from 32560 official human genome gene symbols. The 472 query genes generated from the BA22 sample formulate a genetic network with 15 subgraphs and 36 QQPPIs. However, the mean subgraphs and QQPPIs of 5000 randomly generated networks was 7.54 and 18.54 respectively. The number of QQPPIs of the BA22 genetic network ranked top 4.52% of all randomly generated networks, which shows that the BA22 genetic network unlikely to be the result of chance.

In the search of schizophrenia specific pathways, there are consistent results compared with Sun J. et al. (2010) that there are 4 among 8 pathways involving in or related to the immune system [7] including the glucocorticoid receptor regulatory network, the Fc-epsilon receptor I signaling in mast cells, the NF-kappaB pathway and IL-10 signaling. The interleukin family (IL1, IL2, IL3, IL4, IL5, IL6, IL8, IL10, IL12, IL23 and IL27) pathways also implicate significant schizophrenic pathways. These immune-related pathways with significant p-value (<0.01) supports the autoimmune hypothesis of schizophrenia [61].

Schizophrenia and the immune system

A higher prevalence of several autoimmune disorders has been reported in schizophrenic patients. A growing evidence of researches suggests that the human immune system is associated with the susceptibility and increased risk for schizophrenia, which alterations in the inflammation process and cytokine production have been focused as important mediators in the inflammatory process. Alteration of the immune system and increased level of cytokine are also associated with schizophrenia [62, 63]. However, evidence for common genetic susceptibility between schizophrenia and autoimmune disorders is mostly indirect and not intuitive. On the molecular level, schizophrenia and autoimmune disorders seem to share specific genes with family predisposition [64].

Accumulated evidence has identified abnormalities of the immune system in schizophrenia patients. Neuroinfalmmatory and arachidonic acid cascade markers are increased in schizophrenic patients [65]. Dysregulation of the alternative complement pathway in schizophrenia patients provides evidence that the imbalance of immune system contribute to schizophrenia [66].

Putative association of SIRPB1-LCK-SYK genes in SHCN

The LCK gene encodes a 56-kDa protein-tyrosine kinase, predominantly expressed in T lymphocytes, crucial for initiating T cell antigen receptor (TCR) signal transduction pathways is associated with phosphorylation of the T cell antigen receptor(TCR) by tyrosine kinase which is an essential step in the activation of T cell [67]. Isothiazolinones is a kind of fungicidal and bactericidal effect with properties of broad spectrum, which can quickly inhibit microbe growth, leading to the death of microbes. It is also a novel inhibitors of p56(LCK), which is identified to inhibited kinase activity [68]. TCR-induced stimulation of T cells led to simultaneous phosphorylation of p56(LCK) residues at Y505 and Y394 [69]. Serial activation of the tyrosine kinases LCK and ZAP-70 initiates signaling downstream of the T-cell receptor. ZAP-70 and SYK which is essential for B-cell receptor signaling, share a unique domain structure for protein kinases and undergo conformational change on binding to doubly phosphorylated ITAM peptide [70].

In summary, the SIRPB1 mediated LCK and SYK gene activation are associated with schizophrenia related to BA22 tissue specific gene. The figure illustrates the use of genetic network analysis as an explanation for potential mechanisms of schizophrenic pathway.

Schizophrenia and IL-2/TREM-1 pathway

The cellular and molecular module for immune system involving IL-2 pathway and TREM-1/DAP12 pathway were proposed for potential susceptibility for schizophrenia. Recent study proposed strong evidence of the association between schizophrenia and immune functions, elevated levels of inflammation in the dorsolateral prefrontal cortex has been found. To find specific immune patterns in schizophrenia raises the possibility of developing a disease mechanism. Based on the finding of overactive immune system in the brains of schizophrenia, suppression treatment targets in the immune system would cast the future of novel research which introduced a whole new range of treatment possibilities [71].

BDNF and schizophrenia

Brain-derived neurotropic factor (BDNF) plays an important role in the susceptibility of schizophrenia which is involved in the neurodevelopmental abnormalities of the brain and influences the neuroplasticity in schizophrenia [7274]. Low BDNF is associated with schizophrenia [75, 76]. N-methyl-d-aspartate(NMDA) receptor dysfunction mediated glutamatergic system underlying abnormalities in serum BDNF level and NMDA receptor hypofunction contribute to one of the etiology of schizophrenia [77, 78]. Although BDNF does not appear in either the SHCN or the clique network, a simple BDNF network (Figure 6) is proposed and postulated the important mechanism for schizophrenia.

Figure 6
figure 6

BDNF-related PPI network. The PPI of BDNF and NARG2 through EVAVL1 interacts with UBC to cascade the molecular function of cell cycle. Genome-wide association studies indicate significance for the SNP in the ELAVL2 gene associated with schizophrenia. The ELAVL family of RNA-binding proteins regulates gene expression with the implication of biological processes of cancer, and plays an important role in schizophrenia. The KIAA0101 gene interacts with NFkappaB which is important in the BDNF pathway. KIAA0101 is critical in modulating schizophrenic pathway.

BDNF binds to the TrkB receptor with the presynaptic glutamate to the NMDA receptor activate cascades of PI3K, Akt, and Ras pathway, which formulates the PPI of BDNF and NARG2 through EVAVL1 [78, 79] then interacts with UBC to cascade the molecular function of cell cycle. Genome-wide association studies indicate significance for the SNP in the ELAVL2 gene associated with schizophrenia [80]. The extent to which the ELAVL family of RNA-binding proteins regulates gene expression with the implication of biological processes of cancer plays an important role in schizophrenia.

KIAA0101 and schizophrenia

The KIAA0101 gene has the highest centrality and closeness. KIAA0101 has been observed in a variety of human malignancies and plays a key factor in DNA repair and apoptosis in cell cycle regulation. High-level KIAA0101 expression was also identified as an independent prognostic factor for determining postoperative adjuvant treatments for non-small cell lung carcinoma[81].

The over-expression of KIAA0101 was involved in tumor progression through inhibiting the transcriptional activity of the TP53 gene [82]. KIAA0101 functions as a regulator, promoting cell survival in hepatocellular carcinoma through the regulation of TP53. Suppression of the KIAA0101 function is likely to develop novel cancer therapeutic drugs. In SHCN, KIAA0101 interacts with RUVBL2 which is over-expressed in schizophrenic genes alone with NME2 (yellow region in Figure 2), which indicates the important role in the modulation of disease genes. RUVBL2 is a novel repressor of ARF transcription, ARF is the second most commonly inactivated tumor suppressor gene behind TP53. The genes including KIAA0101, RUVBL2, ARF and TP53 are crucial for schizophrenia.

The KIAA0101 gene is an important cancer gene. In fact, it has PPIs with many other schizophrenic candidate genes. It indirectly interacts with the TP53 gene through the interaction with RUVBL2 and ARF genes [83]. The KIAA0101 gene also interacts with NFkappaB which is important in the BDNF pathway [84]. Through the observation of the PPI network, it is postulated that the KIAA0101 is critical in the modulation of schizophrenic pathways.


It is not clear that cross-talk among various schizophrenic candidate genes is essential for the explanation of the etiology of schizophrenia. The aim of this research is to evaluate the candidate genes chosen from significant over- and under-expression genes of schizophrenia and hepatocellular carcinoma. The SHCGene formulates the SHCN, including the QQPPI, L1PPI and clique network as a major approach for the discovery of potential complexes and pathways. Investigation of potential schizophrenic pathways with the IL-2/TREM-1 pathway reveals possible complexes or drugs responsible for novel treatment of schizophrenia and hepatocellular carcinoma.


  1. Muller N, Schwarz MJ: Neuroimmune-endocrine crosstalk in schizophrenia and mood disorders. Expert Rev Neurother. 2006, 6: 1017-1038. 10.1586/14737175.6.7.1017.

    Article  PubMed  Google Scholar 

  2. Muller N, Schwarz MJ: [Immunology in schizophrenic disorders]. Nervenarzt. 2007, 78: 253-256. 10.1007/s00115-006-2108-9. 258-260, 262-253

    Article  PubMed  CAS  Google Scholar 

  3. Richard MD, Brahm NC: Schizophrenia and the immune system: pathophysiology, prevention, and treatment. Am J Health Syst Pharm. 2012, 69: 757-766. 10.2146/ajhp110271.

    Article  PubMed  CAS  Google Scholar 

  4. 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.171.

    Article  PubMed  CAS  Google Scholar 

  5. Sun J, Kuo PH, Riley BP, Kendler KS, Zhao Z: Candidate genes for schizophrenia: a survey of association studies and gene ranking. Am J Med Genet B Neuropsychiatr Genet. 2008, 147B: 1173-1181. 10.1002/ajmg.b.30743.

    Article  PubMed  Google Scholar 

  6. Sun J, Han L, Zhao Z: Gene- and evidence-based candidate gene selection for schizophrenia and gene feature analysis. Artif Intell Med. 2010, 48: 99-106. 10.1016/j.artmed.2009.07.009.

    Article  PubMed  PubMed Central  Google Scholar 

  7. 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.0011351.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Ayalew M, Le-Niculescu H, Levey DF, Jain N, Changala B, Patel SD, Winiger E, Breier A, Shekhar A, Amdur R, et al: Convergent functional genomics of schizophrenia: from comprehensive understanding to genetic risk prediction. Mol Psychiatry. 2012, 17: 887-905. 10.1038/mp.2012.37.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  9. Fond G, Macgregor A, Attal J, Larue A, Brittner M, Ducasse D, Capdevielle D: Antipsychotic drugs: pro-cancer or anti-cancer? A systematic review. Med Hypotheses. 2012, 79: 38-42. 10.1016/j.mehy.2012.03.026.

    Article  PubMed  CAS  Google Scholar 

  10. Jablensky A, Lawrence D: Schizophrenia and cancer: is there a need to invoke a protective gene?. Arch Gen Psychiatry. 2001, 58: 579-580. 10.1001/archpsyc.58.6.579.

    Article  PubMed  CAS  Google Scholar 

  11. Barak Y, Achiron A, Mandel M, Mirecki I, Aizenberg D: Reduced cancer incidence among patients with schizophrenia. Cancer. 2005, 104: 2817-2821. 10.1002/cncr.21574.

    Article  PubMed  Google Scholar 

  12. Mortensen PB: The incidence of cancer in schizophrenic patients. J Epidemiol Community Health. 1989, 43: 43-47. 10.1136/jech.43.1.43.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  13. Lin GM, Chen YJ, Kuo DJ, Jaiteh LE, Wu YC, Lo TS, Li YH: Cancer Incidence in Patients With Schizophrenia or Bipolar Disorder: A Nationwide Population-Based Study in Taiwan, 1997-2009. Schizophr Bull. 2011

    Google Scholar 

  14. Ji J, Sundquist K, Ning Y, Kendler KS, Sundquist J, Chen X: Incidence of Cancer in Patients With Schizophrenia and Their First-Degree Relatives: A Population-Based Study in Sweden. Schizophr Bull. 2012

    Google Scholar 

  15. Gal G, Goral A, Murad H, Gross R, Pugachova I, Barchana M, Kohn R, Levav I: Cancer in parents of persons with schizophrenia: is there a genetic protection?. Schizophr Res. 2012, 139: 189-193. 10.1016/j.schres.2012.04.018.

    Article  PubMed  Google Scholar 

  16. Hodgson R, Wildgust HJ, Bushe CJ: Cancer and schizophrenia: is there a paradox?. J Psychopharmacol. 2010, 24: 51-60.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Huang da W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.

    Article  PubMed  Google Scholar 

  18. Kanehisa M: The KEGG database. Novartis Found Symp. 2002, 247: 91-101. discussion 101-103, 119-128, 244-152

    Article  PubMed  CAS  Google Scholar 

  19. Kamburov A, Stelzl U, Lehrach H, Herwig R: The ConsensusPathDB interaction database: 2013 update. Nucleic Acids Res. 2013, 41: D793-800. 10.1093/nar/gks1055.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  20. Barnes MR, Huxley-Jones J, Maycox PR, Lennon M, Thornber A, Kelly F, Bates S, Taylor A, Reid J, Jones N, et al: Transcription and pathway analysis of the superior temporal cortex and anterior prefrontal cortex in schizophrenia. J Neurosci Res. 2011, 89: 1218-1227. 10.1002/jnr.22647.

    Article  PubMed  CAS  Google Scholar 

  21. Roessler S, Jia HL, Budhu A, Forgues M, Ye QH, Lee JS, Thorgeirsson SS, Sun Z, Tang ZY, Qin LX, Wang XW: A unique metastasis gene signature enables prediction of tumor relapse in early-stage hepatocellular carcinoma patients. Cancer Res. 2010, 70: 10202-10212. 10.1158/0008-5472.CAN-10-2607.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, Barrette TR, Anstet MJ, Kincead-Beal C, Kulkarni P, et al: Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007, 9: 166-180. 10.1593/neo.07112.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  23. Bader GD, Donaldson I, Wolting C, Ouellette BF, Pawson T, Hogue CW: BIND--The Biomolecular Interaction Network Database. Nucleic Acids Res. 2001, 29: 242-245. 10.1093/nar/29.1.242.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  24. Goel R, Harsha HC, Pandey A, Prasad TS: Human Protein Reference Database and Human Proteinpedia as resources for phosphoproteome analysis. Mol Biosyst. 2012, 8: 453-463. 10.1039/c1mb05340j.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  25. Cesareni G, Chatr-aryamontri A, Licata L, Ceol A: Searching the MINT database for protein interaction information. Curr Protoc Bioinformatics. 2008, Chapter 8:Unit 8 5

    Google Scholar 

  26. Chatr-Aryamontri A, Breitkreutz BJ, Heinicke S, Boucher L, Winter A, Stark C, Nixon J, Ramage L, Kolas N, O'Donnell L, et al: The BioGRID interaction database: 2013 update. Nucleic Acids Res. 2013, 41: D816-823. 10.1093/nar/gks1158.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  27. Xenarios I, Salwinski L, Duan XJ, Higney P, Kim SM, Eisenberg D: DIP, the Database of Interacting Proteins: a research tool for studying cellular networks of protein interactions. Nucleic Acids Res. 2002, 30: 303-305. 10.1093/nar/30.1.303.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  28. Lee SA, Chan CH, Chen TC, Yang CY, Huang KC, Tsai CH, Lai JM, Wang FS, Kao CY, Huang CY: POINeT: protein interactome with sub-network analysis and hub prioritization. BMC Bioinformatics. 2009, 10: 114-10.1186/1471-2105-10-114.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Zheng S, Zhao Z: GenRev: exploring functional relevance of genes in molecular networks. Genomics. 2012, 99: 183-188. 10.1016/j.ygeno.2011.12.005.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  30. Ma Q, Chirn GW, Cai R, Szustakowski JD, Nirmala NR: Clustering protein sequences with a novel metric transformed from sequence similarity scores and sequence alignments with neural networks. BMC Bioinformatics. 2005, 10: 114-

    Google Scholar 

  31. Chen TC, Lee SA, Chan CH, Juang YL, Hong YR, Huang YH, Lai JM, Kao CY, Huang CY: Cliques in mitotic spindle network bring kinetochore-associated complexes to form dependence pathway. Proteomics. 2009, 9: 4048-4062. 10.1002/pmic.200900231.

    Article  PubMed  CAS  Google Scholar 

  32. Jupe S, Akkerman JW, Soranzo N, Ouwehand WH: Reactome - a curated knowledgebase of biological pathways: megakaryocytes and platelets. J Thromb Haemost. 2012

    Google Scholar 

  33. Schaefer CF, Anthony K, Krupa S, Buchoff J, Day M, Hannay T, Buetow KH: PID: the Pathway Interaction Database. Nucleic Acids Res. 2009, 37: D674-679. 10.1093/nar/gkn653.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Kamburov A, Pentchev K, Galicka H, Wierling C, Lehrach H, Herwig R: ConsensusPathDB: toward a more complete picture of cell biology. Nucleic Acids Res. 2011, 39: D712-717. 10.1093/nar/gkq1156.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  35. Sreenivasaiah PK, Rani S, Cayetano J, Arul N, Kim do H: IPAVS: Integrated Pathway Resources, Analysis and Visualization System. Nucleic Acids Res. 2012, 40: D803-808. 10.1093/nar/gkr1208.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  36. Knox C, Law V, Jewison T, Liu P, Ly S, Frolkis A, Pon A, Banco K, Mak C, Neveu V, et al: DrugBank 3.0: a comprehensive resource for 'omics' research on drugs. Nucleic Acids Res. 2011, 39: D1035-1041. 10.1093/nar/gkq1126.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  37. Ruepp A, Waegele B, Lechner M, Brauner B, Dunger-Kaltenbach I, Fobo G, Frishman G, Montrone C, Mewes HW: CORUM: the comprehensive resource of mammalian protein complexes--2009. Nucleic Acids Res. 2010, 38: D497-501. 10.1093/nar/gkp914.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  38. Zhao M, Sun J, Zhao Z: TSGene: a web resource for tumor suppressor genes. Nucleic Acids Res. 2013, 41: D970-976. 10.1093/nar/gks937.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  39. Catts VS, Catts SV: Apoptosis and schizophrenia: is the tumour suppressor gene, p53, a candidate susceptibility gene?. Schizophr Res. 2000, 41: 405-415. 10.1016/S0920-9964(99)00077-8.

    Article  PubMed  CAS  Google Scholar 

  40. Cui DH, Jiang KD, Jiang SD, Xu YF, Yao H: The tumor suppressor adenomatous polyposis coli gene is associated with susceptibility to schizophrenia. Mol Psychiatry. 2005, 10: 669-677. 10.1038/

    Article  PubMed  CAS  Google Scholar 

  41. Aoki K, Taketo MM: Adenomatous polyposis coli (APC): a multi-functional tumor suppressor gene. J Cell Sci. 2007, 120: 3327-3335. 10.1242/jcs.03485.

    Article  PubMed  CAS  Google Scholar 

  42. Li T, Kon N, Jiang L, Tan M, Ludwig T, Zhao Y, Baer R, Gu W: Tumor suppression in the absence of p53-mediated cell-cycle arrest, apoptosis, and senescence. Cell. 2012, 149: 1269-1283. 10.1016/j.cell.2012.04.026.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  43. Shiota S, Tochigi M, Shimada H, Ohashi J, Kasai K, Kato N, Tokunaga K, Sasaki T: Association and interaction analyses of NRG1 and ERBB4 genes with schizophrenia in a Japanese population. J Hum Genet. 2008, 53: 929-935. 10.1007/s10038-008-0332-9.

    Article  PubMed  Google Scholar 

  44. Aberg KA, Liu Y, Bukszar J, McClay JL, Khachane AN, Andreassen OA, Blackwood D, Corvin A, Djurovic S, Gurling H, et al: A comprehensive family-based replication study of schizophrenia genes. JAMA Psychiatry. 2013, 70: 1-9. 10.1001/jamapsychiatry.2013.1322.

    Article  Google Scholar 

  45. Yang Y, Xiao Z, Chen W, Sang H, Guan Y, Peng Y, Zhang D, Gu Z, Qian M, He G, et al: Tumor suppressor gene TP53 is genetically associated with schizophrenia in the Chinese population. Neurosci Lett. 2004, 369: 126-131. 10.1016/j.neulet.2004.07.068.

    Article  PubMed  CAS  Google Scholar 

  46. Lee SA, Tsao TT, Yang KC, Lin H, Kuo YL, Hsu CH, Lee WK, Huang KC, Kao CY: Construction and analysis of the protein-protein interaction networks for schizophrenia, bipolar disorder, and major depression. BMC Bioinformatics. 2011, 12 (Suppl 13): S20-10.1186/1471-2105-12-S13-S20.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  47. Yin F, Liu X, Li D, Wang Q, Zhang W, Li L: Tumor suppressor genes associated with drug resistance in ovarian cancer (Review). Oncol Rep. 2013, 30: 3-10.

    PubMed  CAS  Google Scholar 

  48. Guan X, Wang LE, Liu Z, Sturgis EM, Wei Q: Association between a rare novel TP53 variant (rs78378222) and melanoma, squamous cell carcinoma of head and neck and lung cancer susceptibility in non-Hispanic Whites. J Cell Mol Med. 2013

    Google Scholar 

  49. Pagano M: Cell cycle regulation by the ubiquitin pathway. FASEB J. 1997, 11: 1067-1075.

    PubMed  CAS  Google Scholar 

  50. Pavri R, Lewis B, Kim TK, Dilworth FJ, Erdjument-Bromage H, Tempst P, de Murcia G, Evans R, Chambon P, Reinberg D: PARP-1 determines specificity in a retinoid signaling pathway via direct modulation of mediator. Mol Cell. 2005, 18: 83-96. 10.1016/j.molcel.2005.02.034.

    Article  PubMed  CAS  Google Scholar 

  51. Specenier P, Vermorken JB: Cetuximab: its unique place in head and neck cancer treatment. Biologics. 2013, 7: 77-90.

    PubMed  CAS  PubMed Central  Google Scholar 

  52. Arnold D, Stein A: New developments in the second-line treatment of metastatic colorectal cancer: potential place in therapy. Drugs. 2013, 73: 883-891. 10.1007/s40265-013-0076-5.

    Article  PubMed  CAS  Google Scholar 

  53. Huang YC, Liu CY, Lu HJ, Liu HT, Hung MH, Hong YC, Hsiao LT, Gau JP, Liu JH, Hsu HC, et al: Comparison of prognostic models for patients with diffuse large B-cell lymphoma in the rituximab era. Ann Hematol. 2013

    Google Scholar 

  54. Plosker GL, Figgitt DP: Rituximab: a review of its use in non-Hodgkin's lymphoma and chronic lymphocytic leukaemia. Drugs. 2003, 63: 803-843. 10.2165/00003495-200363080-00005.

    Article  PubMed  CAS  Google Scholar 

  55. Noguchi T, Ritter G, Nishikawa H: Antibody-based therapy in colorectal cancer. Immunotherapy. 2013, 5: 533-545. 10.2217/imt.13.35.

    Article  PubMed  CAS  Google Scholar 

  56. Connolly R, Nguyen NK, Sukumar S: Molecular Pathways: Current Role and Future Directions of the Retinoic Acid Pathway In Cancer Prevention and Treatment. Clin Cancer Res. 2013

    Google Scholar 

  57. Smith W, Saba N: Retinoids as chemoprevention for head and neck cancer: where do we go from here?. Crit Rev Oncol Hematol. 2005, 55: 143-152. 10.1016/j.critrevonc.2005.02.003.

    Article  PubMed  Google Scholar 

  58. Carrera S, Cuadrado-Castano S, Samuel J, Jones GD, Villar E, Lee SW, Macip S: Stra6, a retinoic acid-responsive gene, participates in p53-induced apoptosis after DNA damage. Cell Death Differ. 2013

    Google Scholar 

  59. Wan C, Shi Y, Zhao X, Tang W, Zhang M, Ji B, Zhu H, Xu Y, Li H, Feng G, He L: Positive association between ALDH1A2 and schizophrenia in the Chinese population. Prog Neuropsychopharmacol Biol Psychiatry. 2009, 33: 1491-1495. 10.1016/j.pnpbp.2009.08.008.

    Article  PubMed  CAS  Google Scholar 

  60. Wan C, Yang Y, Li H, La Y, Zhu H, Jiang L, Chen Y, Feng G, He L: Dysregulation of retinoid transporters expression in body fluids of schizophrenia patients. J Proteome Res. 2006, 5: 3213-3216. 10.1021/pr060176l.

    Article  PubMed  CAS  Google Scholar 

  61. Jones AL, Mowry BJ, Pender MP, Greer JM: Immune dysregulation and self-reactivity in schizophrenia: do some cases of schizophrenia have an autoimmune basis?. Immunol Cell Biol. 2005, 83: 9-17. 10.1111/j.1440-1711.2005.01305.x.

    Article  PubMed  CAS  Google Scholar 

  62. Fineberg AM, Ellman LM: Inflammatory Cytokines and Neurological and Neurocognitive Alterations in the Course of Schizophrenia. Biol Psychiatry. 2013

    Google Scholar 

  63. Frydecka D, Beszlej A, Karabon L, Pawlak-Adamska E, Tomkiewicz A, Partyka A, Jonkisz A, Monika SB, Kiejna A: The role of genetic variations of immune system regulatory molecules CD28 and CTLA-4 in schizophrenia. Psychiatry Res. 2013

    Google Scholar 

  64. Ferentinos P, Dikeos D: Genetic correlates of medical comorbidity associated with schizophrenia and treatment with antipsychotics. Curr Opin Psychiatry. 2012, 25: 381-390. 10.1097/YCO.0b013e3283568537.

    Article  PubMed  Google Scholar 

  65. Rao JS, Kim HW, Harry GJ, Rapoport SI, Reese EA: Increased neuroinflammatory and arachidonic acid cascade markers, and reduced synaptic proteins, in the postmortem frontal cortex from schizophrenia patients. Schizophr Res. 2013

    Google Scholar 

  66. Martins-de-Souza D, Gattaz WF, Schmitt A, Rewerts C, Maccarrone G, Dias-Neto E, Turck CW: Prefrontal cortex shotgun proteome analysis reveals altered calcium homeostasis and immune system imbalance in schizophrenia. Eur Arch Psychiatry Clin Neurosci. 2009, 259: 151-163. 10.1007/s00406-008-0847-2.

    Article  PubMed  Google Scholar 

  67. Rossy J, Owen DM, Williamson DJ, Yang Z, Gaus K: Conformational states of the kinase Lck regulate clustering in early T cell signaling. Nat Immunol. 2013, 14: 82-89.

    Article  PubMed  CAS  Google Scholar 

  68. Trevillyan JM, Chiou XG, Ballaron SJ, Tang QM, Buko A, Sheets MP, Smith ML, Putman CB, Wiedeman P, Tu N, et al: Inhibition of p56(lck) tyrosine kinase by isothiazolones. Arch Biochem Biophys. 1999, 364: 19-29. 10.1006/abbi.1999.1099.

    Article  PubMed  CAS  Google Scholar 

  69. Nyakeriga AM, Garg H, Joshi A: TCR-induced T cell activation leads to simultaneous phosphorylation at Y505 and Y394 of p56(lck) residues. Cytometry A. 2012, 81: 797-805.

    Article  PubMed  Google Scholar 

  70. Yan Q, Barros T, Visperas PR, Deindl S, Kadlecek TA, Weiss A, Kuriyan J: Structural basis for activation of ZAP-70 by phosphorylation of the SH2-kinase linker. Mol Cell Biol. 2013

    Google Scholar 

  71. Fillman SG, Cloonan N, Catts VS, Miller LC, Wong J, McCrossin T, Cairns M, Weickert CS: Increased inflammatory markers identified in the dorsolateral prefrontal cortex of individuals with schizophrenia. Mol Psychiatry. 2013, 18: 206-214. 10.1038/mp.2012.110.

    Article  PubMed  CAS  Google Scholar 

  72. Huang TL: Effects of antipsychotics on the BDNF in schizophrenia. Curr Med Chem. 2013, 20: 345-350.

    PubMed  CAS  Google Scholar 

  73. Xu MQ, St Clair D, Feng GY, Lin ZG, He G, Li X, He L: BDNF gene is a genetic risk factor for schizophrenia and is related to the chlorpromazine-induced extrapyramidal syndrome in the Chinese population. Pharmacogenet Genomics. 2008, 18: 449-457. 10.1097/FPC.0b013e3282f85e26.

    Article  PubMed  CAS  Google Scholar 

  74. Muglia P, Vicente AM, Verga M, King N, Macciardi F, Kennedy JL: Association between the BDNF gene and schizophrenia. Mol Psychiatry. 2003, 8: 146-147. 10.1038/

    Article  PubMed  CAS  Google Scholar 

  75. Zhang XY, Liang J, Chen da C, Xiu MH, Yang FD, Kosten TA, Kosten TR: Low BDNF is associated with cognitive impairment in chronic patients with schizophrenia. Psychopharmacology (Berl). 2012, 222: 277-284. 10.1007/s00213-012-2643-y.

    Article  CAS  Google Scholar 

  76. Green MJ, Matheson SL, Shepherd A, Weickert CS, Carr VJ: Brain-derived neurotrophic factor levels in schizophrenia: a systematic review with meta-analysis. Mol Psychiatry. 2011, 16: 960-972. 10.1038/mp.2010.88.

    Article  PubMed  CAS  Google Scholar 

  77. Kantrowitz JT, Javitt DC: N-methyl-d-aspartate (NMDA) receptor dysfunction or dysregulation: the final common pathway on the road to schizophrenia?. Brain Res Bull. 2010, 83: 108-121. 10.1016/j.brainresbull.2010.04.006.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  78. Snyder MA, Gao WJ: NMDA hypofunction as a convergence point for progression and symptoms of schizophrenia. Front Cell Neurosci. 2013, 10: 114-

    Google Scholar 

  79. Favalli G, Li J, Belmonte-de-Abreu P, Wong AH, Daskalakis ZJ: The role of BDNF in the pathophysiology and treatment of schizophrenia. J Psychiatr Res. 2012, 46: 1-11. 10.1016/j.jpsychires.2011.09.022.

    Article  PubMed  Google Scholar 

  80. Yamada K, Iwayama Y, Hattori E, Iwamoto K, Toyota T, Ohnishi T, Ohba H, Maekawa M, Kato T, Yoshikawa T: Genome-wide association study of schizophrenia in Japanese population. PLoS One. 2011, 6: e20468-10.1371/journal.pone.0020468.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  81. Kato T, Daigo Y, Aragaki M, Ishikawa K, Sato M, Kaji M: Overexpression of KIAA0101 predicts poor prognosis in primary lung cancer patients. Lung Cancer. 2012, 75: 110-118. 10.1016/j.lungcan.2011.05.024.

    Article  PubMed  Google Scholar 

  82. Liu L, Chen X, Xie S, Zhang C, Qiu Z, Zhu F: Variant 1 of KIAA0101, overexpressed in hepatocellular carcinoma, prevents doxorubicin-induced apoptosis by inhibiting p53 activation. Hepatology. 2012, 56: 1760-1769. 10.1002/hep.25834.

    Article  PubMed  CAS  Google Scholar 

  83. Xie C, Wang W, Yang F, Wu M, Mei Y: RUVBL2 is a novel repressor of ARF transcription. FEBS Lett. 2012, 586: 435-441. 10.1016/j.febslet.2012.01.026.

    Article  PubMed  CAS  Google Scholar 

  84. Li K, Ma Q, Shi L, Dang C, Hong Y, Wang Q, Li Y, Fan W, Zhang L, Cheng J: NS5ATP9 gene regulated by NF-kappaB signal pathway. Arch Biochem Biophys. 2008, 479: 15-19. 10.1016/

    Article  PubMed  CAS  Google Scholar 

Download references


We would like to give special thanks to Beitou Branch, Tri-Service General Hospital for providing space and facilities for necessary computing machine and discussion. We were also grateful that Information Technology Center, Kainan University provided high-performance computing resources. We also thanked DAVID, GenRev, CORUM and DrugBank for helpful webtools available online for data analysis and retrieving.


Publication of this research was funded by Project NSC100-2218-E-424-001 of National Science Council, Taiwan.

This article has been published as part of BMC Genomics Volume 14 Supplement 5, 2013: Twelfth International Conference on Bioinformatics (InCoB2013): Computational biology. The full contents of the supplement are available online at

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Sheng-An Lee or Cheng-Yan Kao.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

KCH interpreted the results, drafted the manuscript, and contributed to the design of the bioinformatics analysis tools. SAL programmed the bioinformatics analysis tools and carried out the data analysis. KCY, HL and WKL assisted in the interpretation of results. SAL, TTHT and CYK conceived the study and participate in coordination and management of the research project.

Electronic supplementary material


Additional file 1: The SHCGene contains the over- and under- expression candidate genes from the BA22 sample and the Roessler liver 2 sample(XLS 84 KB)


Additional file 2: The schizophrenia candidate genes and tumor suppression genes are compared from different literature reviews. Sun J. et al. (2008), Allen NC. et al. (2008), Ayalew M. et al. (2012), BA22 microarray and hepatocellular carcinoma samples overlap with 716 genes in tumor suppression genes which implicate that schizophrenia candidate genes related to the tumor suppression genes share a common biological regulation. (XLS 51 KB)


Additional file 3: The co-expressed mediator genes generated from SHCN L1PPI. The common mediators of the L1PPI networks of Sun J. et al.(2008), Allen NC. et al. (2008), Ayalew M. et al. (2012) and BA22 microarray samples. The results implicate the importance of mediator genes rather than query genes as key roles of disease susceptibility. (XLS 191 KB)


Additional file 4: The pathways of BA22 and Roessler liver 2 samples are ranked by the p-value and FDR-adjusted p-value by the Benjamini-Hochberg procedure. (XLS 76 KB)

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( ) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Cite this article

Huang, KC., Yang, KC., Lin, H. et al. Analysis of schizophrenia and hepatocellular carcinoma genetic network with corresponding modularity and pathways: novel insights to the immune system. BMC Genomics 14 (Suppl 5), S10 (2013).

Download citation

  • Published:

  • DOI:


  • Schizophrenia
  • Schizophrenic Patient
  • Tazarotene
  • Fenretinide
  • Retinoid Acid