3.1. Identification of DEGs in Epilepsy and IBD
In the GSE59071 dataset, 1243 DEGs (UC tissues vs. control) and 697 DEGs (CD tissues vs. control) were identified. In the GSE134697 dataset, 1846 DEGs (hippocampus vs. control) were identified. The shared DEGs of datasets GSE134697 and GSE59071 were 78 genes, including 56 upregulated and 22 downregulated genes (Figure 1, Table 1).
3.2. GO Enrichment Analysis of DEGs
GO enrichment analysis showed various upregulated and downregulated DEGs in biological processes (BP), molecular functions (MF), and cellular components (CC), respectively. In BP, upregulated DEGs were significantly enriched in regulation of vasculature development, response to molecule of bacterial origin, response to lipopolysaccharide, positive regulation of cytokine production, and positive regulation of cell-cell adhesion (Figure 2a, Table 2). In MF, DEGs were enriched in glycosaminoglycan binding, cytokine activity, receptor ligand activity, and signaling receptor activator activity (Figure 2b, Table 2). In CC, upregulated DEGs were significantly enriched in secretory granule membrane and collagen-containing extracellular (Figure 2c, Table 2). Downregulated DEGs were only enriched in MF (Figure 2d, Table 2). More detailed GO enrichment results are shown in Table 2. The amounts of upregulated genes involved in the different processes of BP, CC, and MF were much higher than that of downregulated genes. These results suggest that comorbidity of epilepsy and IBD may undergo complex metabolic activities involving cytokine associated inflammatory response and other complex biological process.
3.3. KEGG Enrichment Analysis of DEGs
KEGG enrichment analysis showed that upregulated genes were enriched in 26 pathways, while the downregulated genes were not enriched in any pathway (Table 3). Upregulated DEGs were enriched in pathways of rheumatoid arthritis, graft-versus-host disease, type I diabetes mellitus, hematopoietic cell lineage, leishmaniasis, intestinal immune network for IgA production, toll-like receptor signaling pathway, inflammatory bowel disease, and cytokine-cytokine receptor interaction pathway and so on (Figure 3). These pathways are mainly involved in inflammatory or proinflammatory response.
3.4. GSEA
GSEA showed that the dataset GSE134697 is enriched in these pathways: antigen processing and presentation, cytokine-cytokine receptor interaction, graft-versus-host disease, hematopoietic cell lineage, and leishmaniasis infection. As for dataset GSE59071, the expression profiles of both UC and CD samples were enriched in cell adhesion molecules cams, chemokine signaling pathway, complement and coagulation cascades, cytokine-cytokine receptor interaction, and ECM receptor interaction. The results indicated that the cytokine-cytokine receptor interaction pathway is the communal pathway between epilepsy and IBD (Figure 4a, 4b).
3.5. PPI Network Construction, Significant Module Analysis and Candidate Hub Gene Identification
We developed a PPI network containing 51 nodes and 162 edges and visualized using Cytoscape software (Figure 5a). Nodes are drawn in different colors, representing the node degree and control (up or down), there were 56 upregulated genes and 22 downregulated genes. The most significant module with the highest score (module score: 9.111) was shown in Figure 5b, containing 10 genes and 41 edges. GO analysis revealed that these genes are involved in chemokine-mediated signaling pathway, inflammatory response, and the other biological process are showed in Table 4. Genes in the most significant module were selected as candidate hub genes. These genes were all upregulated genes, as follows: CXCL8, CXCR4, ITGAX, CD80, SELL, CD44, CD69, IL1B, CD38, and CXCL2.
3.6. Validation of Candidate Hub genes
We finally identified three hub genes: CXCL8, CXCR4 and ITGAX. CXCL8 is significantly downregulated (adjusted P <0.05 and log2FC =-1.2911696) and the other two genes are slightly downregulated in epilepsy patients with response to VPA treatment, while these three genes are significantly upregulated in GSE134697 and GSE59071 datasets.