Printer Friendly

Identification of Candidate Genes Related to Inflammatory Bowel Disease Using Minimum Redundancy Maximum Relevance, Incremental Feature Selection, and the Shortest-Path Approach.

1. Introduction

Inflammatory bowel disease (IBD) is a common systemic disease that involves the intestinal tissue [1]. It usually refers to chronic conditions that lead to intestinal inflammation and lesions. With the gradual development of inflammation, the intestinal walls become swollen, inflamed, and ulcerogenic [2]. Due to such lesions, several classical symptoms have been considered to be diagnostic indicators. Abdominal pain or cramping, diarrhea multiple times per day, and bloody stools are all classical symptoms of IBD [3]. Such severe symptoms are induced by violent unhealthy inflammation reactions and lesions in the intestinal tissue. Additionally, several complications outside the digestive tract may also be induced by IBD. Mouth sores and skin problems have both been reported in IBD patients [4]. Furthermore, arthritis is also related to IBD, as well as eye problems [5, 6].

As we have mentioned above, several complications have been identified in IBD patients. Such severe complications and related chronic characteristics strongly increase the risk of death [4-6]. In 2013 alone, thousands of people in the world died from IBD [7]. Additionally, IBD has been proven to be associated with colorectal cancer, with a high mortality. Apart from the risk of death, IBD is a lifetime disease, and life with IBD can be quite challenging. The complications associated with IBD and disease relapse severely impact the quality of life [8]. Therefore, the prevention, diagnosis, and treatment of IBD are quite crucial. It is known that IBD is a widespread disease that can develop at any stage of life. However, the disease usually initiates during the teenage years or the early adulthood of the patients [8]. As we mentioned above, genetic factors participate in the initiation and progression of IBD [9,10]. Therefore, people with a family history of IBD are at least ten times more likely to suffer from it. Racial factors also contribute to the morbidity of IBD [11].

Although IBD is a very severe and widespread disease, the essential mechanism behind the disease has not been demonstrated clearly. Most people believe that some types of exogenous materials trigger the initiation of inflammation [12, 13]. However, genetic factors may also contribute to the progression of such disease. Several specific genes have been linked to IBD. Pathogenic genes such as IL23R and IL12B play a crucial role in the intestinal immune system, which may induce the initiation of IBD [14, 15]. Several transcriptional factors also contribute to disease progression. The transcription factor NKX2-3 regulates the correct localization of lymphocytes and may further contribute to the immune response in intestinal tissue that induces IBD [16]. Several genes such as ZNF365 and PTGER4 show diversity in different subtypes of IBD and contribute to IBD through their respective methods and pathways [17,18].

As mentioned above, IBD has several subtypes. Basically, there are two main clinical classifications of IBD: Crohn's disease and ulcerative colitis [17]. Both classifications share the basic symptoms of IBD. However, Crohn's disease can occur anywhere along the digestive tract and typically appears as "skip lesions" between healthy areas [19], while the other type, ulcerative colitis, only involves the colon and rectum. Inflammation and ulcers typically affect only the innermost lining in these areas, with more superficial lesions than those with Crohn's disease [20]. Apart from the differences in clinical symptoms, genetic diversity is also observed between Crohn's disease and ulcerative colitis. Although they share most of the disease-causing genes, genes like ATG16L1, PTGER4, IRGM, and NOD2 have been proven to be specifically related to Crohn's disease but independent with ulcerative colitis [18, 21-23]. The roles of genes such as SLC22A5, ZNF365, and PTPN2 in ulcerative colitis are still unclear, even though they have been proven to be strongly related to Crohn's disease [24, 25].

Because genetic factors have been shown to be related to IBD and its specific subtypes, we developed a new computational method to screen differential expressing genes among different clusters based on a database for Crohn's disease and ulcerative colitis. From the Gene Expression Omnibus (GEO), we obtained the gene expression profiles (information often used to deduce and understand gene functions) for 59 Crohn's disease, 26 ulcerative colitis, and 42 normal samples. Each sample was represented using the expression levels of 12,754 genes. Two feature selection methods, minimum redundancy maximum relevance (mRMR) and incremental feature selection (IFS) [26], and a basic machine learning algorithm, sequential minimal optimization (SMO) [27, 28], were adopted to analyze the gene expression profiles and extract 21 promising candidate genes that could be used to distinguish the samples from the two subtypes of IBD and the normal samples; that is, they may be related to IBD. Furthermore, based on these 21 genes, the shortest-path (SP) approach was employed to identify additional 20 genes in a network constructed using protein-protein interaction (PPI) information. It was concluded that the 41 (21 + 20) genes obtained are closely associated with IBD and can be used to clearly distinguish healthy people from those who have IBD and to identify the subtypes of IBD.

2. Materials and Methods

2.1. Dataset. We downloaded the gene expression profiles of 59 Crohn's disease, 26 ulcerative colitis, and 42 normal samples from GEO under accession number GSE3365 [29]. The expression levels of 12,754 genes were measured using an Affymetrix Human Genome U133A Array. The gene expression profiles were quantile normalized. Each sample was represented using the expression levels of 12,754 genes; that is, each sample was encoded into a 12754-D vector. These features/genes were analyzed to identify the genes that can best discriminate the samples from these three different classes.

2.2. mRMR Method. It is known that some genes can effectively help us discriminate the samples from the three different classes mentioned in Section 2.1, while others offer few or no contributions. To identify these genes, the mRMR method, proposed by Peng et al. [26], was adopted to analyze the gene expression data. The mRMR method employed two criteria, Max-Relevance and Min-Redundancy, to analyze the features. Using the Max-Relevance criterion, the MaxRel feature list can be obtained, in which features are sorted by measuring the relevance between them and sample class labels. Features with high relevance receive high ranks, whereas those with low relevance receive low ranks. It is clearly seen that the rank of a feature in the MaxRel feature list indicates its single contribution to classification. Furthermore, another list, namely, the mRMR feature list, was created using both Max-Relevance and Min-Redundancy criteria. The rank of a feature in this list is determined using the relevance between it and sample class labels and the redundancies between it and the features listed before it. The MaxRel feature list and mRMR feature list in this study were formulated as follows:

MaxRel features list is

[F.sub.MaxRel] = [[f.sub.1.sup.M], [f.sub.2.sup.M], ..., [f.sub.N.sup.M]; (1)

mRMR features list is

[F.sub.mRMR] = [[f.sub.1.sup.m], [f.sub.2.sup.m], ..., [f.sub.N.sup.m]], (2)

where N represents the total number of features. Many investigators have used the mRMR method to analyze various complicated biological systems [30, 31], and it is deemed to be a useful tool for extracting important information from a complicated system. Readers can refer to Peng et al.'s paper [26] for the detailed procedures and principle of this method.

2.3. Prediction Engine. SMO is a type of support vector machine that uses Platt's sequential minimal optimization algorithm to train and optimize the support vector classifier. The kernels can be polynomial or Gaussian [27, 28]. For implementing our method, we employed the classifier SMO implemented in Weka [32] as the prediction engine.

2.4. Tenfold Cross-Validation. Tenfold cross-validation [33] is a type of cross-validation method that is widely used to examine the performance of a classifier on a given dataset. The given dataset is randomly and equally divided into ten partitions. Samples in each partition are singled out in turn as the test data, while other samples are used to train the classifier. Compared to the jackknife test [34, 35], another popular cross-validation method, this method involves a lower amount of computational time and always yields similar results. Thus, it was used in this study for evaluating the performance of the current prediction engine.

2.5. IFS Method. Using the mRMR method, features/genes were sorted and listed in the MaxRel feature list and mRMR feature list. Because the MaxRel feature list sorted features/genes by only measuring their own contributions to classification, the combination of some features/genes with high ranks in this list is not always an optimal combination for classification. The mRMR feature list is more appropriate for this purpose because it further considers the redundancies between features. The IFS method uses the mRMR feature list and the SMO prediction engine to extract the optimal combination of features/genes as biomarkers. First, according to the mRMR feature list [F.sub.mRMR] = [[f.sub.1.sup.m], [f.sub.2.sup.m], ..., [f.sub.N.sup.m]], we constructed N feature set, denoted by [F.sub.1], [F.sub.2], ..., [F.sub.N], where [F.sub.i] = {[f.sub.1.sup.m], [f.sub.2.sup.m], ..., [f.sub.i.sup.m]}; that is, [F.sub.i] contained the top i features in the mRMR feature list. Second, for each [F.sub.i], SMO was executed on the dataset, in which samples were represented using features in [F.sub.i], with its performance evaluated by tenfold cross-validation. Finally, we counted the total prediction accuracy and accuracies for each class. The feature set yielding the highest total prediction accuracy was deemed to be the optimal gene set ([G.sub.optimal]) for IBD, as features in this set may be significant for IBD.

2.6. Network Construction from PPI Information. The optimal gene set [G.sub.optimal] containing some genes closely related to IBD can be obtained using the mRMR and IFS methods. To further mine for other related genes, we constructed a large network from the PPI data and searched for additional candidate genes in the network.

To construct the network, we downloaded the file "protein.links.v9.1.txt.gz" containing the PPI information from STRING (Search Tool for the Retrieval of Interacting Genes/Proteins, version 9.1,, from which the human PPI data were extracted by identifying lines starting with "9606." A total of 2,425,314 human PPIs involving 20,770 proteins represented using Ensembl IDs were obtained. According to STRING ( [36, 37], these PPIs are derived from the following sources: (i) genomic context, (ii) high-throughput experiments, (iii) (conserved) coexpression, and (iv) previous knowledge. Thus, the obtained PPIs contained actual PPIs validated using experiments and predicted PPIs, suggesting that they can be used to widely measure the physical and functional relationships between proteins. Each PPI contained two proteins represented using Ensembl IDs and one score that indicates the strength of the interaction with a range between 150 and 999. The constructed network had 20,770 proteins as nodes. Two nodes were adjacent if and only if the corresponding proteins comprise an interaction that is contained in the 2,425,314 human PPIs. Furthermore, the interaction score was also added to the network. Each edge was assigned a weight defined to be 1,000 minus the corresponding interaction score.

2.7. SP Approach for Searching for Additional Candidates. Network method is an important type of approaches for investigation of disease genes, such as methods based on guilt-by-association (GBA) [38-40] and Random Walk with Restart (RWR) [41-43]. This section proposed another network method for identifying novel disease genes.

It has been elaborated in some previous studies [44-46] that two proteins in an interaction are more likely to share similar functions. It can be induced that the interactive proteins of the proteins encoded by genes in Goptimal are also related to IBD. Furthermore, if we consider a series of proteins [p.sub.1], [p.sub.2], ..., [p.sub.s] such that the consecutive proteins comprise a PPI with a high score and [p.sub.1], [p.sub.s] are proteins encoded by genes in [G.sub.optimal], [p.sub.2], [p.sub.3], ..., [p.sub.s-1] may also be related to IBD. From the construction of the network mentioned in Section 2.6, the corresponding nodes of [p.sub.1], [p.sub.2], ..., [p.sub.s] may comprise a shortest path connecting [p.sub.1] and [p.sub.s]. Therefore, for any two genes in [G.sub.optimal], we searched the shortest path connecting these two genes, thereby collecting a number of shortest paths. Because the endpoints of these paths represented proteins encoded by genes in [G.sub.optimal], genes on these paths may be related to IBD. Thus, we extracted inner nodes on the obtained shortest paths and their corresponding genes can be obtained. To identify novel genes related to IBD, genes in [G.sub.optimal] were excluded from the obtained genes. The remaining genes were called shortest-path genes for convenience. To identify these shortest-path genes, a measurement, namely, the betweenness [47], was recorded for each shortest-path gene, and it was defined to be the number of shortest paths containing the shortest-path gene.

Because some nodes occupied general hubs in the constructed network, the corresponding genes may always be selected even if we searched for the shortest path connecting any pair of randomly selected genes; some of these genes may be selected as the shortest-path genes obtained as described above. In fact, they have few or no associations with IBD. Thus, a permutation test is necessary to control for this type of gene. The procedures used are as follows:

(1) Randomly produce 1,000 gene sets, say [G.sub.1], [G.sub.2], ..., [G.sub.1000], where the size of each set is the same as that of [G.sub.optimal].

(2) For each [G.sub.i], search for all the shortest paths connecting any pair of genes in [G.sub.i] and count the betweenness of the shortest-path gene based on these paths.

(3) A total of 1,000 betweenness scores on 1,000 randomly produced gene sets can be obtained for each shortest-path gene. After comparing the betweenness on [G.sub.optimal], we calculate another measurement, the permutation FDR, for each shortest-path gene, which is defined to be "the number of betweenness scores on randomly produced gene sets that was larger than that on [G.sub.optimal]"/100[degrees].

(4) Because it is implied that shortest-path genes with high permutation FDRs are general hubs in the network and not specific to IBD, those with permutation FDRs larger than or equal to 0.05 are excluded. The remaining genes are termed candidate genes.

To select genes with core relationships with IBD from the candidate genes, the human PPIs and their interaction scores were directly used. For each candidate gene g, we checked the scores of the interactions between g and genes in [G.sub.optimal] and selected the maximum value among them as the maximum interaction score of g. If a candidate gene has a high maximum interaction score, this suggests that it is highly related to at least one gene in [G.sub.optimal], indicating that it is more likely to be related to IBD. As 900 is set to be the threshold of the highest confidence cutoff in STRING, we also set 900 as the threshold for the maximum interaction score; that is, genes with maximum interaction scores no less than 900 were finally selected as the candidate genes in this study.

3. Results and Discussion

3.1. Results of the mRMR and IFS Methods. The mRMR method was executed on the dataset containing 59 Crohn's disease, 26 ulcerative colitis, and 42 normal samples, and each sample was represented using the expression levels of 12,754 genes, thereby yielding the MaxRel feature list and mRMR feature list, which are provided in Supplementary Material I in Supplementary Material available online at

To extract the optimal gene sets for discriminating the samples from two subtypes of IBD and normal samples, the IFS method was used with the mRMR feature list obtained using the mRMR method and SMO as the prediction engine. To reduce computational time and account for the fact that genes with important contributions for discriminating samples from two subtypes of IBD and normal samples are few in number, we only investigated the first 2,000 feature sets. According to the procedures of the IFS method, each feature set can yield four accuracies: three accuracies for three classes and the total prediction accuracy. All of these are provided in Supplementary Material II. Furthermore, an IFS curve was plotted by representing the total prediction accuracy along y-axis and the size of the feature set, that is, the number of features participating in the classification, along x-axis, as shown in Figure 1. It can be seen that the highest total prediction accuracy was 97.64% using the 1170th feature set. The corresponding accuracies for the three classes were 100%, 92.31%, and 97.62%, respectively. Although the accuracies were quite good, the involved features/genes were too many in number, which is not realistic. By carefully checking the IFS curve shown in Figure 1, we observe a sharp increasing trend with more and more features participating in the classification at the beginning of the curve with a rather high total prediction accuracy (93.70%) using the 21st feature set. Then the curve is unstable; increasing trends and decreasing trend occur in succession. Thus, we believe that the first 21 features in the mRMR feature list are more important for discriminating the samples from two subtypes of IBD and normal samples than others and set the optimal gene set [G.sub.optimal] to be the 21st feature set. These 21 genes are listed in Table 1. The associations between these 21 genes and IBD are elaborated in Section 3.4. However, some important IBD-related genes may not be omitted using the mRMR and IFS methods. Based on these genes, the SP approach was applied to discover additional genes related to IBD, which is described in the following sections.

3.2. Shortest-Path Genes. As mentioned in Section 3.1, 21 genes were obtained and deemed to be important for discriminating the samples from the two subtypes of IBD and the normal samples. To further identify more candidate genes, we constructed a large network, as described in Section 2.6. These 21 genes were mapped to 20 genes in the network. We searched for all shortest paths connecting any pair of 20 genes, resulting in 190 paths. The graph of these 190 paths is shown in Figure 2, where we can see that there are 110 Ensembl genes on these paths other than the 21 genes obtained in Section 3.1. By mapping to their gene symbols, we obtained 107 shortest-path genes. These genes and their betweenness are listed in Supplementary Material III.

3.3. Additional Candidate Genes. According to Section 2.7, a permutation test was executed to exclude general genes in the network. The obtained permutation FDRs of 107 shortest-path genes are also provided in Supplementary Material III. By setting the threshold of the permutation FDR to be 0.05, 57 candidate genes were obtained, which are listed in Supplementary Material IV.

To select the core genes among the 57 candidate genes, the maximum interaction score of each candidate gene was calculated. These values are also provided in Supplementary Material IV. The threshold of the maximum interaction score was set to 900, resulting in 20 candidate genes that are listed in Table 2.

3.4. Analysis of Candidate Genes. Based on feature analysis of 59 Crohn's disease, 26 ulcerative colitis, and 42 normal samples, we obtained 21 genes, listed in Table 1, which maybe related to IBD and can help distinguish healthy people from those who have two subtypes of IBD. Furthermore, according to the above 21 genes and the SP approach, we obtained additional 20 candidate genes, listed in Table 2. These genes are also thought to be related to IBD. This section provides some evidence for this claim.

We combined two candidate gene sets and analyzed the biological meaning behind them using Functional Annotation Bioinformatics Microarray Analysis (DAVID) (version 6.7, [48]. The obtained results are provided in Supplementary Material V. According to the results yielded by DAVID, crucial gene ontology (GO) terms and KEGG pathways like hsa04660 (T cell receptor signaling pathway), GO: 0001775 (cell activation), and GO: 0045449 (regulation oftranscription) were screened out to be enriched by 41 candidate genes. In addition, the results also gave clues for clustering 41 candidate genes into some groups, which provided convenience for analyzing candidate genes.

3.4.1. Candidate Genes Contributing to T Cell Receptor Signaling Pathway (hsa04660). As mentioned above, IBD is a severe disease induced by inflammation reactions [1]. Considering the core regulatory role of T cells in immune system, it is quite reasonable that various candidate genes contribute to such pathway. Based on SP approach, we identified a specific gene FOS. It is also a tumor-associated gene, which encodes a leucine zipper protein that can dimerize with proteins of the JUN family, thereby forming the transcription factor complex AP-1 [49]. Related to crucial pathways such as NF-[kappa]B and MAPK, FOS is quite significant in inflammation initiation, especially in the digestive tract [50, 51]. Another calcium-associated gene PLCG1 was also discovered. PLC1 participates in the intracellular transduction of receptor-mediated tyrosine kinase activators and may participate in the inflammation reaction through a specific function [52, 53]. LCK (also known as p56lck) is another predicted IBD-related gene that encodes a functional tyrosine kinase. Similar to ZAP70, LCK also regulates the metabolism and maturation of T cells and may further regulate the inflammation process [54,55]. In terms of IBD, LCK has been reported to be associated with ulcerative colitis but not with Crohn's disease [56]. Our predicted gene ZAP70 is a protein tyrosine kinase participating in the development and activation of T cells [57, 58]. ZAP70 has been reported to be associated with a specific subtype of IBD, Crohn's disease, but not ulcerative colitis [59]. Therefore, the expression level of ZAP70 can be a useful biomarker for distinguishing different subtypes of IBD.

While based on the mRMR and IFS method, we also identified a group of candidate genes. Among them, CD247 (rank 3 in the mRMR feature list) and CD4 (rank 19 in the mRMR feature list) are both crucial genes for T cells and have been confirmed to further regulate the inflammation reaction [60,61]. Our predicted gene CD4 is characteristically expressed in IBD [62]. However, CD4 has also been reported as a differentially expressed gene in Crohn's disease and ulcerative colitis, and it may further serve as a new biomarker for distinguishing these two diseases [63]. Based on our functional clustering, various screened and predicted genes also are enriched in a similar GO term, GO: 0042101 which describes T cell receptor complex as a cellular component, validating the enrichment of T cell receptor signaling pathway of our screened out IBD associated genes.

3.4.2. Candidate Genes Contributing to Cell Activation (GO: 0001775). In Table 2, a highly conserved monooxygenase-associated protein YWHAZ was identified as a functional protein in intestinal bowel disease. Such gene is a crucial housekeeping gene that has been proven to be a suitable normalizer for bowel inflammation and cancer [60]. As a functional factor of innate immune response, which is also crucial in intestinal tissues, TLR4 is predicted to be associated with IBD. TLR4 has been reported as a crucial factor in the innate immune barrier of the intestine [61]. Such factors can be activated by specific factors (FFA, etc.) and further induce the initiation of IBD [64, 65]. Another thrombin-associated gene, F2 (coagulation factor II) was also identified by the SP approach. F2 and THBD are both coagulation-associated genes. The coagulation process is reported to be associated with Crohn's disease but not with ulcerative colitis, which reflects the differences between various subtypes of IBD [66]. There are two major subtypes of IBD: Crohn's disease and ulcerative colitis. Some candidate genes yielded by mRMR and IFS methods may distinguish these two subtypes. PF4 (rank 5 in the mRMR feature list), a crucial diagnostic biomarker for IBD, has been clearly reported to be overexpressed in Crohn's disease and thought that it does not play a clear role in ulcerative colitis [67, 68]. PF4 can also separate IBD from normal inflammation, which is crucial for diagnosis [69].

3.4.3. Candidate Genes Contributing to Regulation of Transcription (GO: 0045449). Among the 41 candidate genes, quite a lot of genes contribute to the regulation of transcription, implying the complicated endogenous pathological factors of IBD on multiple levels. Based on mRMR and IFS methods, the candidate gene ZNF207 (rank 1 in the mRMR feature list), which is a specific microtubule-associated zinc finger protein, may regulate the inflammation of IBD [70]. As a regulator of mitotic chromosome alignment, ZNF207 has been reported to be related to another type of inflammation disorder, chronic obstructive pulmonary disease (COPD). Since both COPD and IBD are localized inflammation involving the mucosal tissue, ZNF207 as our candidate gene may also contribute to inflammatory bowel disease [71]. As a T cell regulator, EGR3 (rank 8 in the mRMR feature list) was also identified. As a member of the EGR family, EGR3 maybe a crucial transcriptional factor for T cells, with high similarity with EGR2 [72]. SLTM (rank 4 in the mRMR feature list) acts as a general inhibitor of transcription that eventually leads to apoptosis via the regulation of telomere [73]. Because IBD is associated with abnormal cell death, SLTM may participate in IBD through the regulation of the apoptosis of intestinal cells [74]. CNOT8 (rank 14 in the mRMR feature list) is a significant predicted gene that interacts with BTG, the regulator of the cell cycle, especially in B cells [75]. Therefore, CNOT8 may indirectly participate in the intestinal inflammation reaction [75, 76]. TH1L (rank 13 in the mRMR feature list), as a negative elongation factor complex member C/D (NELFCD), promotes the proliferation of intestinal cells and has been proved to induce carcinoma progression [77]. As a regulator of B cells, HMGB1 (rank 9 in the mRMR feature list) and its homolog HMGB2 constitute a complex that is differentially expressed in Crohn's disease and ulcerative colitis [78, 79]. Such a complex has also been reported as a new marker of IBD and may be a sensitive marker of mucosal inflammation [80]. As we have mentioned above, our predicted gene CD4 is characteristically expressed in IBD [62]. However, CD4 has also been reported as a differentially expressed gene in Crohn's disease and ulcerative colitis, and it may further serve as a new biomarker for distinguishing these two diseases [63]. UBE2I (rank 12 in the mRMR feature list) also regulates the proliferation of intestinal cells [81]. Unlike FOLR1, which we will analyze below, UBE2I is a major part of the SUMO ligases and further promotes the proliferation of intestinal cells via multiple means even under pathological conditions [81, 82].

For the candidate genes obtained by the SP approach, HCFC1 is a functional nuclear activator. As a unique cleavage signal, it has been reported to be associated with cell cycle regulation and may have a specific function in tumorigenesis [83, 84]. As a part of the CCR4-NOT complex, CNOT1 is a crucial immune associated gene that is a major cellular mRNA deadenylase and has been reported to participate in several processes related to immune reactions [85]. Regulated by the CCR4-NOT complex, a crucial microRNA, miR155, has been reported to be directly associated with inflammation, which may further reveal the tight connection between CNOT1 and the inflammation reaction [86, 87]. Such functional genes may also participate in the initiation of inflammation and tumors. CNOT4 is also a part of the CCRNOT complex, and CNOT4 may act similarly to CNOT1 and contribute to the regulation of the immune reaction [85]. As a functional factor of innate immune response, which is also crucial in intestinal tissues, TRAK1 is a regulatory gene that maybe related to endosome-to-lysosome trafficking and EGF-EGFR interaction [88]. Such an EGF-EGFR interaction is definitely associated with the initiation of bowel inflammation [89]. The candidate gene HDAC1 regulates the acetylation of specific genes and further participates in the regulation of corresponding functions [90]. Gene acetylation and deacetylation are functional regulatory methods for cell metabolism, which have been identified in IBD [9193]. Therefore, HDAC1 may play a regulatory role in the initiation and progression of intestinal bowel diseases. BTG1 is a functional regulatory gene associated with cell growth and differentiation. Similar to FASLG, it also regulates the apoptosis of specific target cells and may further regulate specific cytokines associated with inflammation such as IFN[gamma] [94]. Histone deacetylase is commonly used to modify the epigenetic status and regulate gene expression [95]. RUNX1, known as runt-related transcription factor 1, is quite crucial in the development of normal hematopoiesis as a part of CBF (core binding factor). Associated with T cell function and TGF-[beta], RUNX1 has been proven to be quite crucial in inflammation initiation [96, 97]. Considering the strong relationship between IBD and immune reaction, RUNX1, which regulates the function of T cells, may also participate in the initiation of IBD [98].

3.4.4. Candidate Genes Contributing to Protein Kinase Cascade (GO: 0007243). Four functional genes have been clustered into such group. Genes like F2, ZAP70, and TLR4 have already been analyzed above. The gene MARK2 (rank 2 in the mRMR feature list) may also contribute to the initiation and progression of IBD by interfering with the protein kinase cascade. Inflammation is a basic pathological process regulated by the immune system [99]. Therefore, the immune system plays an irreplaceable role in IBD [100]. Several predicted genes have been confirmed to be associated with the immune system and participate in the immune reaction. MARK2 is a serine/threonine-protein kinase that is the major regulator of cell polarity in epithelial cells, including intestinal epithelial cells. Since immune cells in intestinal system have been proven to be regulated by such gene, the abnormal expression and effect of MARK2 may contribute to the unusual activation of focal inflammatory reaction in the digestive system, which may further promote IBD [101].

3.4.5. Candidate Genes Contributing to Intrinsic to Plasma Membrane (GO: 0031226). Among the candidate genes obtained by the SP approach, THBD is an endothelial-specific type I membrane receptor that binds thrombin [102]. As a specific protein in coagulation mechanisms, this receptor has also been reported as a potential inflammation mediator and may have a specific function in IBD [103,104]. We also predicted a specific member of the TNF family, FASLG, as a candidate gene. FASLG has been proven to be involved in the induction of apoptosis triggered by binding to FAS [105]. Members of the TNF family have been widely reported to participate in IBDs by regulating the apoptosis of specific local cells [106,107].

For candidate genes listed in Table 1, FOLR1 (rank 6 in the mRMR feature list), the folate receptor, participates in intestinal inflammation via the regulation of folate. Folate is associated with cell apoptosis in bowel tissues and has been reported to be crucial in colonic epithelial cell proliferation implying its potential role in inflammatory bowel diseases [108, 109]. SLC22A4 (rank 18 in the mRMR feature list) is a homolog of SLC22A5, which has been reported to be crucial in Crohn's disease and is also overexpressed in this disease [110]. However, just like SLC22A5, SLC22A4 has not been confirmed to be overexpressed in ulcerative colitis [111, 112]. The genes mentioned above can distinguish IBD subtypes at the genetic level and may serve as new markers for the classification of inflammation in intestinal tissues. As a receptor of significant biological signals, LEPROT (rank 10 in the mRMR feature list) encodes a crucial receptor of GH and has been reported to be associated with the initiation of inflammation in the intestine in mice [113]. IBD has been regarded to be the result of immune systematic disorders and autoimmune reactions [1, 100]. Another gene, CLEC1B (rank 16 in the mRMR feature list), also participates in the development of IBD via the regulation of the intestinal immune system, especially the proliferation of NK cells and the formation of lymph nodes [114]. Apart from NK cells, activated cell is also a major part of the immune system and has been shown to be related to IBD [115,116].

3.4.6. Candidate Genes Contributing to Apoptosis (GO: 0006915). Some candidate genes have been confirmed to participate in the apoptosis processes during the pathological processes of IBD. Apart from genes like SLTM, LCK, F2, FASLG, and BLCAP which we have just analyzed above, the candidate gene RHOT2 (rank 21 in the mRMR feature list), a mitochondrial GTPase involved in mitochondrial trafficking, has been proven to be crucial regulator of [Ca.sup.2+] in T cells. Thus, RHOT2 may also contribute to IBD [117]. IBD is a common disease involving the digestive system, especially the intestinal tissue [1]. However, IBD has also been shown to be associated with carcinoma in the digestive system, especially colorectal cancer [100]. Several of our predicted genes are also involved in tumor initiation, where cells may have mutated in precancerous lesions, including severe IBD. Most of these genes are related to cell proliferation. BLCAP (rank 7 in the mRMR feature list), which was first reported in bladder cancer, regulates the proliferation of cells that are quite common in intestinal tissue of IBD patients [118].

3.4.7. Candidate Genes Contributing to Regulation of Cell Proliferation (GO: 0042127). Among the candidate genes listed in Table 2, several have been confirmed to contribute to cell proliferation, implying the potential role of that during IBD initiation and progression. STK11, a functional serine/threonine kinase, regulates the polarity of cells and may participate in tumor suppression [119]. NF-[kappa]B is a crucial transcriptional factor that participates in the inflammation process [120]. STK11 (also known as LKB1) directly regulates the function of NF-[kappa]B and is definitely associated with inflammation [121]. STK11 also regulates the proliferation and maturation of intestinal cells, which indirectly reflects the regulatory function of STK11 in intestinal tissues. The calcium binding protein S100A6 is also on our predicted list, and it is located in the cytoplasm and nucleus of a wide range of cells. S100A6 regulates the progression of the cell cycle and the differentiation of specific cells [122]. Considering the tight relationship between IBD and cancer, some of our predicted genes are also associated with tumor initiation [123]. We also predict as a candidate gene a serine proteinase inhibitor SERPINE1, which encodes the principal inhibitor of tissue plasminogen activator (tPA) and urokinase (uPA). Tissue plasminogen activator and urokinase are both associated with inflammation and the process of wound healing [124]. SERPINE1 and proteins in the downstream of its specific pathway have also been reported to be directly associated with IBD as a functional regulator [125,126]. As a candidate gene, we also predicted an angiogenesis-associated gene VEGFC, which regulates angiogenesis and endothelial cell growth [127,128]. VEGFC has been reported to participate in several intestinal disorders including IBD and some specific digestive tract cancers [129,130].

3.4.8. Other Candidate Genes. Four candidate genes obtained by mRMR and IFS methods were not clustered into any above group. The candidate gene ANXA11 (rank 20 in the mRMR feature list) is a predicted gene that regulates the autoimmune reaction. Such a gene has been reported to be related to several autoimmune disorders and may further participate in intestinal inflammation [131]. As a part of the MLL complex, OGT (rank 15 in the mRMR feature list) regulates the cell cycle of intestinal cells, including immune cells [132]. Therefore, the abnormity of the OGT gene may induce IBD in various downstream pathways. Another candidate gene, USPL1 (rank 17 in the mRMR feature list), also participates in intestinal inflammation reaction via the SUMO complex [133]. Large-scale mapping of human protein-protein interactions by mass spectrometry revealed several genes associated with inflammation, especially in the intestine [76]. The last gene, HIST1H2AC (rank 11 in the mRMR feature list), is also a candidate gene for cancer. Such gene was first reported in breast cancer and regulates the proliferation of tissue cells, similar to BLCAP [134].

3.5. Comparison of Other Methods. To indicate the effectiveness of the proposed method and the reliability of the obtained genes, we compared our method with other methods. Before making the comparison, 77 validated IBD-related genes were retrieved from [135], which are provided in Supplementary Material VI. These genes were used to test the results yielded by our method and other methods.

DisGeNET (Verison4.0) [136] is a discovery platform that collects gene-disease associations from several public data sources and the literature. Here, it was used to search IBD-related genes. The obtained material is provided in Supplementary Material VII, from which we extracted 100 genes with high confidence (score > 0.1) as the predicted genes of this method. DAVID 6.7 ( [48] was employed again to analyze the biological meanings behind the validated genes, predicted genes by our method, and predicted genes by DisGeNET. The enriched gene ontology (GO) terms and KEGG pathways for three gene lists are listed in Supplementary Material V. It can be observed that 209 GO terms and KEGG pathways were enriched by 77 validated genes, while, for predicted genes by our method and DisGeNET, we obtained 154 and 314 GO terms and KEGG pathways, respectively. For the 154 GO terms and KEGG pathways enriched by 41 predicted genes of our method, 51 (51/154 = 33.12%) were also enriched by 77 validated genes, while there were 117 (117/314 = 37.26%) GO terms and KEGG pathways enriched by both 77 validated genes and 100 predicted genes of DisGeNET.

At a first glance, the performance of the DisGeNET is superior to our method. However, our method still has its advantages. According to our method, 21 genes were extracted by analyzing the gene expression profiles using mRMR, IFS, and SMO methods. In fact, these genes can only help us to distinguish two subtypes of IBD (rather than all subtypes of IBD) and normal samples. Thus, they are parts of IBD-related genes even if they are really IBD-related genes. 20 additional candidate genes were further obtained based on these genes, thereby accessing 41 predicted genes. These 41 predicted genes, in fact, are deemed to be related to two subtypes of IBD rather other all IBD subtypes. On the other hand, 100 predicted genes yielded by DisGeNET considered all subtypes of IBD. It is an important reason why DisGeNET gave the better performance. However, the performance of our method is only slightly lower than that of DisGeNET. Therefore, we believe that the proposed method is still quite effective and the obtained genes can be important and reliable materials for the investigation of IBD.

4. Conclusions

This contribution provides a novel computational method to identify genes related to IBD, which consists of two main steps: (1) analyzing the gene expression profiles and extracting important genes for IBD and (2) applying the shortestpath approach to the network constructed using proteinprotein interactions and identifying additional related genes. By analyzing the obtained genes, it is concluded that they have special relationships with IBD, implying that our method is effective. It is also believed that our method has potential applicability for the investigation of other diseases.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this article.


This paper is supported by National Natural Science Foundation of China (31371335).


[1] D. J. Mulder, A. J. Noble, C. J. Justinich, and J. M. Duffin, "A tale of two diseases: the history of inflammatory bowel disease," Journal of Crohn's and Colitis, vol. 8, no. 5, pp. 341-348, 2014.

[2] G. Barbara, C. Cremon, and V Stanghellini, "Inflammatory bowel disease and irritable bowel syndrome: similarities and differences," Current Opinion in Gastroenterology, vol. 30, no. 4, pp. 352-358, 2014.

[3] G. Bassotti, E. Antonelli, V. Villanacci, M. Salemme, M. Coppola, and V. Annese, "Gastrointestinal motility disorders in inflammatory bowel diseases," World Journal of Gastroenterology, vol. 20, no. 1, pp. 37-44, 2014.

[4] M. Zippi, C. Corrado, R. Pica et al., "Extraintestinal manifestations in a large series of Italian inflammatory bowel disease patients," World Journal of Gastroenterology, vol. 20, no. 46, pp. 17463-17467, 2014.

[5] A. Marineata, E. Rezus, C. Mihai, and C. C. Prelipcean, "Extra intestinal manifestations and complications in inflammatory bowel disease," Revista Medico-Chirurgicala a Societatii de Medici si Naturalisti din Iasi, vol. 118, no. 2, pp. 279-288, 2014.

[6] P. L. Lakatos, L. Lakatos, L. S. Kiss, L. Peyrin-Biroulet, A. Schoepfer, and S. Vavricka, "Treatment of extraintestinal manifestations in inflammatory bowel disease," Digestion, vol. 86, no. I, pp. 28-35, 2012.

[7] S. Singh, I. J. Kullo, D. S. Pardi, and E. V. Loftus Jr., "Epidemiology, risk factors and management of cardiovascular diseases in IBD," Nature Reviews Gastroenterology & Hepatology, vol. 12, no. 1, pp. 26-35, 2015.

[8] J. Ruel, D. Ruane, S. Mehandru, C. Gower-Rousseau, and J.F. Colombel, "IBD across the age spectrum--is it the same disease?" Nature Reviews Gastroenterology and Hepatology, vol. II, no. 2, pp. 88-98, 2014.

[9] E. JaZwinska-Tarnawska, I. Jeskowiak, E. Waszczuk et al., "Genetic polymorphism of ABCB1 gene (C3435T) in patients with inflammatory bowel diseases. Is there any gender dependency?" Pharmacological Reports, vol. 67, no. 2, pp. 294-298, 2015.

[10] C. Jakobsen, I. Cleynen, P. S. Andersen et al., "Genetic susceptibility and genotype-phenotype association in 588 Danish children with inflammatory bowel disease," Journal of Crohn's & Colitis, vol. 8, no. 7, pp. 678-685, 2014.

[11] G. C. Nguyen, C. A. Chong, and R. Y. Chong, "National estimates of the burden of inflammatory bowel disease among racial and ethnic groups in the United States," Journal of Crohn's and Colitis, vol. 8, no. 4, pp. 288-295, 2014.

[12] Y. Zhang and Y. Y. Li, "Inflammatory bowel disease: pathogenesis," World Journal of Gastroenterology, vol. 20, no. 1, pp. 91-99, 2014.

[13] P. Flanagan, B. J. Campbell, and J. M. Rhodes, "Bacteria in the pathogenesis of inflammatory bowel disease," Biochemical Society Transactions, vol. 39, no. 4, pp. 1067-1072, 2011.

[14] C. W. Lees, J. C. Barrett, M. Parkes, and J. Satsangi, "New IBD genetics: common pathways with other diseases," Gut, vol. 60, no. 12, pp. 1739-1753, 2011.

[15] J. Glas, J. Seiderer, J. Wagner et al., "Analysis of IL12B gene variants in inflammatory bowel disease," PLOS ONE, vol. 7, no. 3, Article ID e34349, 2012.

[16] G. John, J. P. Hegarty, W. Yu et al., "NKX2-3 variant rs11190140 is associated with IBD and alters binding of NFAT," Molecular Genetics and Metabolism, vol. 104, no. 1-2, pp. 174-179, 2011.

[17] K. L. VanDussen, T. C. Liu, D. Li et al., "Genetic variants synthesize to produce paneth cell phenotypes that define subtypes of Crohn's disease," Gastroenterology, vol. 146, no. 1, pp. 200-209, 2014.

[18] K. Kabashima, T. Saji, T. Murata et al., "The prostaglandin receptor EP4 suppresses colitis, mucosal damage and CD4 cell activation in the gut," Journal of Clinical Investigation, vol. 109, no. 7, pp. 883-893, 2002.

[19] G. Patman, "Crohn's disease: suppression of p21Rac1 signalling contributes to skip-lesion phenotype in Crohn's disease," Nature Reviews Gastroenterology & Hepatology, vol. 11, no. 6, article 332, 2014.

[20] J. J. Rumessen, "Ultrastructure of interstitial cells of Cajal at the colonic submuscular border in patients with ulcerative colitis," Gastroenterology, vol. 111, no. 6, pp. 1447-1455, 1996.

[21] J. Hampe, A. Franke, P. Rosenstiel et al., "A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1," Nature Genetics, vol. 39, no. 2, pp. 207-211, 2007

[22] M. Parkes, J. C. Barrett, N. J. Prescott et al., "Sequence variants in the autophagy gene IRGM and multiple other replicating loci contribute to Crohn's disease susceptibility," Nature Genetics, vol. 39, no. 7, pp. 830-832, 2007.

[23] T. Hisamatsu, M. Suzuki, H.-C. Reinecker, W. J. Nadeau, B. A. McCormick, and D. K. Podolsky, "CARD15/NOD2 functions as an antibacterial factor in human intestinal epithelial cells," Gastroenterology, vol. 124, no. 4, pp. 993-1000, 2003.

[24] E. Leung, J. Hong, A. G. Fraser, T. R. Merriman, P. Vishnu, and G. W. Krissansen, "Polymorphisms in the organic cation transporter genes SLC22A4 and SLC22A5 and Crohn's disease in a New Zealand Caucasian cohort," Immunology and Cell Biology, vol. 84, no. 2, pp. 233-236, 2006.

[25] D. F. McCole, "Regulation of epithelial barrier function by the inflammatory bowel disease candidate gene, PTPN2," Annals of the New York Academy of Sciences, vol. 1257, no. 1, pp. 108-114, 2012.

[26] H. Peng, F. Long, and C. Ding, "Feature selection based on mutual information: criteria of max-dependency, max-relevance, and min-redundancy," IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 27, no. 8, pp. 1226-1238, 2005.

[27] J. Platt, Fast Training of Support Vector Machines Using Sequential Minimal Optimization, MIT Press, Cambridge, Mass, USA, 1998.

[28] S. S. Keerthi, S. K. Shevade, C. Bhattacharyya, and K. R. K. Murthy, "Improvements to Platt's SMO algorithm for SVM classifier design," Neural Computation, vol. 13, no. 3, pp. 637-649, 2001.

[29] M. E. Burczynski, R. L. Peterson, N. C. Twine et al., "Molecular classification of Crohn's disease and ulcerative colitis patients using transcriptional profiles in peripheral blood mononuclear cells," The Journal of Molecular Diagnostics, vol. 8, no. 1, pp. 51-61, 2006.

[30] L. Chen, C. Chu, and K. Feng, "Predicting the types of metabolic pathway of compounds using molecular fragments and sequential minimal optimization," Combinatorial Chemistry & High Throughput Screening, vol. 19, no. 2, pp. 136-143, 2016.

[31] H. Mohabatkar, M. Mohammad Beigi, and A. Esmaeili, "Prediction of GABAA receptor proteins using the concept of Chou's pseudo-amino acid composition and support vector machine," Journal of Theoretical Biology, vol. 281, no. 1, pp. 18-23, 2011.

[32] I. H. Witten and E. Frank, Data Mining: Practical Machine Learning Tools and Techniques, Morgan Kaufmann, San Francisco, Calif, USA, 2005.

[33] R. Kohavi, "A study of cross-validation and bootstrap for accuracy estimation and model selection," in Proceedings of the International Joint Conference on Artificial Intelligence, Lawrence Erlbaum Associates Ltd, Quebec, Canada, 1995.

[34] L. Chen, J. Lu, N. Zhang, T. Huang, and Y.-D. Cai, "A hybrid method for prediction and repositioning of drug Anatomical Therapeutic Chemical classes," Molecular BioSystems, vol. 10, no. 4, pp. 868-877, 2014.

[35] L. Chen, W.-M. Zeng, Y.-D. Cai, K.-Y. Feng, and K.-C. Chou, "Predicting anatomical therapeutic chemical (ATC) classification of drugs by integrating chemical-chemical interactions and similarities," PLOS ONE, vol. 7, no. 4, Article ID e35254, 2012.

[36] C. von Mering, M. Huynen, D. Jaeggi, S. Schmidt, P. Bork, and B. Snel, "STRING: a database of predicted functional associations between proteins," Nucleic Acids Research, vol. 31, no. 1, pp. 258-261, 2003.

[37] A. Franceschini, D. Szklarczyk, S. Frankild et al., "STRING v9.1: protein-protein interaction networks, with increased coverage and integration," Nucleic Acids Research, vol. 41, no. 1, pp. D808-D815, 2013.

[38] M. Oti, B. Snel, M. A. Huynen, and H. G. Brunner, "Predicting disease genes using protein-protein interactions," Journal of Medical Genetics, vol. 43, no. 8, pp. 691-698, 2006.

[39] M. Krauthammer, C. A. Kaufmann, T. C. Gilliam, and A. Rzhetsky, "Molecular triangulation: bridging linkage and molecular-network information for identifying candidate genes in Alzheimer's desease," Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 42, pp. 15148-15153, 2004.

[40] L. Franke, H. van Bakel, L. Fokkens, E. D. de Jong, M. Egmont-Petersen, and C. Wijmenga, "Reconstruction of a functional human gene network, with an application for prioritizing positional candidate genes," The American Journal of Human Genetics, vol. 78, no. 6, pp. 1011-1025, 2006.

[41] S. Kohler, S. Bauer, D. Horn, and P. N. Robinson, "Walking the Interactome for Prioritization of Candidate Disease Genes," The American Journal of Human Genetics, vol. 82, no. 4, pp. 949-958, 2008.

[42] R. Jiang, M. Gan, and P. He, "Constructing a gene semantic similarity network for the inference of disease genes," BMC Systems Biology, vol. 5, supplement 2, article no. S2, 2011.

[43] H. Shi, J. Xu, G. Zhang et al., "Walking the interactome to identify human miRNA-disease associations through the functional link between miRNA targets and disease genes," BMC Systems Biology, vol. 7, no. 1, article 101, 2013.

[44] M. Jiang, Y. Chen, Y. Zhang et al., "Identification of hepatocellular carcinoma related genes with k-th shortest paths in a protein-protein interaction network," Molecular BioSystems, vol. 9, no. 11, pp. 2720-2728, 2013.

[45] L. Chen, Z. H. Xing, T. Huang, Y. Shu, G. Huang, and H.-P. Li, "Application of the shortest path algorithm for the discovery of breast cancer-related genes," Current Bioinformatics, vol. 11, no. 1, pp. 51-58, 2016.

[46] T. Gui, X. Dong, R. Li, Y. Li, and Z. Wang, "Identification of hepatocellular carcinoma-related genes with a machine learning and network analysis," Journal of Computational Biology, vol. 22, no. 1, pp. 63-71, 2015.

[47] M. Kitsak, S. Havlin, G. Paul, M. Riccaboni, F. Pammolli, and H. E. Stanley, "Betweenness centrality of fractal and nonfractal scale-free model networks and tests on real networks," Physical Review E--Statistical, Nonlinear, and Soft Matter Physics, vol. 75, no. 5, part 2, Article ID 056115, 2007

[48] D. W. Huang, B. T. Sherman, and R. A. Lempicki, "Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources," Nature Protocols, vol. 4, no. 1, pp. 44-57, 2009.

[49] C. Li, H. Li, S. Wang et al., "The c-Fos and c-Jun from Litopenaeus vannamei play opposite roles in Vibrio parahaemolyticus and white spot syndrome virus infection," Developmental and Comparative Immunology, vol. 52, no. 1, pp. 26-36, 2015.

[50] D. Thummuri, M. K. Jeengar, S. Shrivastava et al., "Thymoquinone prevents RANKL-induced osteoclastogenesis activation and osteolysis in an in vivo model of inflammation by suppressing NF-[kappa]B and MAPK Signalling," Pharmacological Research, vol. 99, pp. 63-73, 2015.

[51] M. G. Welch, K. G. Margolis, Z. Li, and M. D. Gershon, "Oxytocin regulates gastrointestinal motility, inflammation, macromolecular permeability, and mucosal maintenance in mice," American Journal of Physiology--Gastrointestinal and Liver Physiology, vol. 307, no. 8, pp. G848-G862, 2014.

[52] C. Auesukaree, H. Tochio, M. Shirakawa, Y. Kaneko, and S. Harashima, "Plc1p, Arg82p, and Kcs1p, enzymes involved in inositol pyrophosphate synthesis, are essential for phosphate regulation and polyphosphate accumulation in Saccharomyces cerevisiae," Journal of Biological Chemistry, vol. 280, no. 26, pp. 25127-25133, 2005.

[53] D. Engelberg, R. Perlman, and A. Levitzki, "Transmembrane signaling in Saccharomyces cerevisiae as a model for signaling in metazoans: state of the art after 25 years," Cellular Signalling, vol. 26, no. 12, pp. 2865-2878, 2014.

[54] Y. J. Chiang and R. J. Hodes, "Regulation of T cell development by c-Cbl: essential role of Lck," International Immunology, vol. 27, no. 5, pp. 245-251, 2015.

[55] M. D. Perron, S. Chowdhury, I. Aubry, E. Purisima, M. L. Tremblay, and H. U. Saragovi, "Allosteric noncompetitive small molecule selective inhibitors of CD45 tyrosine phosphatase suppress T-cell receptor signals and inflammation in vivo," Molecular Pharmacology, vol. 85, no. 4, pp. 553-563, 2014.

[56] L. S. Toy, X. Y. Yio, A. Lin, S. Honig, and L. Mayer, "Defective expression of gp180, a novel CD8 ligand on intestinal epithelial cells, in inflammatory bowel disease," Journal of Clinical Investigation, vol. 100, no. 8, pp. 2062-2071,1997

[57] Z. Liao, L. Zhou, C. Wang et al., "Characteristics of TCR[zeta], ZAP70, and Fc[epsilon]RI[gamma] Gene Expression in Patients with T- and NK/T-Cell Lymphoma," DNA and Cell Biology, vol. 34, no. 3, pp. 201-207, 2015.

[58] C. Sinclair, M. Ono, and B. Seddon, "A Zap70-dependent feedback circuit is essential for efficient selection of CD4 lineage thymocytes," Immunology and Cell Biology, vol. 93, no. 4, pp. 406-416, 2015.

[59] D. Bouzid, H. Fourati, A. Amouri et al., "Association of ZAP70 and PTPN6, but not BANK1 or CLEC2D, with inflammatory bowel disease in the tunisian population," Genetic Testing and Molecular Biomarkers, vol. 17, no. 4, pp. 321-326, 2013.

[60] M. Krzystek-Korpacka, D. Diakowska, J. Bania, and A. Gamian, "Expression stability of common housekeeping genes is differently affected by bowel inflammation and cancer: implications for finding suitable normalizers for inflammatory bowel disease studies," Inflammatory Bowel Diseases, vol. 20, no. 7, pp. 1147-1156, 2014.

[61] W. Wang, T. Xia, and X. Yu, "Wogonin suppresses inflammatory response and maintains intestinal barrier function via TLR4-MyD88-TAK1-mediated NF-[kappa]B pathway in vitro," Inflammation Research, vol. 64, no. 6, pp. 423-431, 2015.

[62] C. S. De Almeida, V. Andrade-Oliveira, N. O. S. Camara, J. F. Jacysyn, and E. L. Faquim-Mauro, "Crotoxin from Crotalus durissus terrificus is able to down-modulate the acute intestinal inflammation in mice," PLoS ONE, vol. 10, no. 4, Article ID e0121427, 2015.

[63] G. Brandhorst, S. Weigand, C. Eberle et al., "CD4+ immune response as a potential biomarker of patient reported inflammatory bowel disease (IBD) activity," Clinica Chimica Acta, vol. 421, pp. 31-33, 2013.

[64] R. A. Gupta, M. N. Motiwala, N. G. Dumore, K. R. Danao, and A. B. Ganjare, "Effect of piperine on inhibition of FFA induced TLR4 mediated inflammation and amelioration of acetic acid induced ulcerative colitis in mice," Journal of Ethnopharmacology, vol. 164, pp. 239-246, 2015.

[65] A. T. Cao, S. Yao, A. T. Stefka et al., "TLR4 regulates IFN-[gamma] and IL-17 production by both thymic and induced Foxp3+ Tregs during intestinal inflammation," Journal of Leukocyte Biology, vol. 96, no. 5, pp. 895-905, 2014.

[66] D. Kohoutova, M. Pecka, M. Cihak, J. Cyrany, J. Maly, and J. Bures, "Prevalence of hypercoagulable disorders in inflammatory bowel disease," Scandinavian Journal of Gastroenterology, vol. 49, no. 3, pp. 287-294, 2014.

[67] T. Bennike, S. Birkelund, A. Stensballe, and V. Andersen, "Biomarkers in inflammatory bowel diseases: current status and proteomics identification strategies," World Journal of Gastroenterology, vol. 20, no. 12, pp. 3231-3244, 2014.

[68] M.-A. Meuwis, M. Fillet, L. Lutteri et al., "Proteomics for prediction and characterization of response to infliximab in Crohn's disease: a pilot study," Clinical Biochemistry, vol. 41, no. 12, pp. 960-967, 2008.

[69] M.-A. Meuwis, M. Fillet, P. Geurts et al., "Biomarker discovery for inflammatory bowel disease, using proteomic serum profiling," Biochemical Pharmacology, vol. 73, no. 9, pp. 1422-1433, 2007.

[70] H. Jiang, X. He, S. Wang et al., "A microtubule-associated zinc finger protein, BuGZ, regulates mitotic chromosome alignment by ensuring Bub3 stability and kinetochore targeting," Developmental Cell, vol. 28, no. 3, pp. 268-281, 2014.

[71] S. Bhattacharya, S. Srisuma, D. L. DeMeo et al., "Molecular biomarkers for quantitative and discrete COPD phenotypes," American Journal of Respiratory Cell and Molecular Biology, vol. 40, no. 3, pp. 359-367, 2009.

[72] T. Okamura, K. Fujio, M. Shibuya et al., "[CD4.sup.+][CD25.sup.-][LAG3.sup.+] regulatory T cells controlled by the transcription factor Egr-2," Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 33, pp. 13974-13979, 2009.

[73] R. J. Giannone, H. W. McDonald, G. B. Hurst, R.-F. Shen, Y. Wang, and Y. Liu, "The protein network surrounding the human telomere repeat binding factors TRF1, TRF2, and POT1," PLoS ONE, vol. 5, no. 8, Article ID e12407, 2010.

[74] M. Zemljic, B. Pejkovic, I. Krajnc, and S. Lipovsek, "Biological pathways involved in the development of infammatory bowel disease," Wiener Klinische Wochenschrift, vol. 126, no. 19-20, pp. 626-633, 2014.

[75] Y. Du, P. Liu, W. Zang et al., "BTG3 upregulation induces cell apoptosis and suppresses invasion in esophageal adenocarcinoma," Molecular and Cellular Biochemistry, vol. 404, no. 1-2, pp. 31-38, 2015.

[76] R. M. Ewing, P. Chu, F. Elisma et al., "Large-scale mapping of human protein-protein interactions by mass spectrometry," Molecular Systems Biology, vol. 3, article no. 89, 2007

[77] B. Carvalho, C. Postma, S. Mongera et al., "Multiple putative oncogenes at the chromosome 20q amplicon contribute to colorectal adenoma to carcinoma progression," Gut, vol. 58, no. 1, pp. 79-89, 2009.

[78] M. McDonnell, Y. Liang, A. Noronha et al., "Systemic toll-like receptor ligands modify B-cell responses in human inflammatory bowel disease," Inflammatory Bowel Diseases, vol. 17, no. 1, pp. 298-307, 2011.

[79] R. Vitali, L. Stronati, A. Negroni et al., "Fecal HMGB1 is a novel marker of intestinal mucosal inflammation in pediatric inflammatory bowel disease," The American Journal of Gastroenterology, vol. 106, no. 11, pp. 2029-2040, 2011.

[80] H. Takaishi, T. Kanai, A. Nakazawa et al., "Anti-high mobility group box 1 and box 2 non-histone chromosomal proteins (HMGB1/HMGB2) antibodies and anti-Saccharomyces cerevisiae antibodies (ASCA): accuracy in differentially diagnosing UC and CD and correlation with inflammatory bowel disease phenotype," Journal of Gastroenterology, vol. 47, no. 9, pp. 969-977, 2012.

[81] M. D. Demarque, K. Nacerddine, H. Neyretkahn et al., "Sumoylation by Ubc9 regulates the stem cell compartment and structure and function of the intestinal epithelium in mice," Gastroenterology, vol. 140, no. 1, pp. 286-296, 2011.

[82] N. S. Belaguli, M. Zhang, A.-H. Garcia, and D. H. Berger, "PIAS1 Is a GATA4 SUMO ligase that regulates GATA4-dependent intestinal promoters independent of SUMO ligase activity and GATA4 sumoylation," PLoS ONE, vol. 7, no. 4, Article ID e35717, 2012.

[83] P. Zhou, Z. Wang, X. Yuan et al., "Mixed Lineage Leukemia 5 (MLL5) protein regulates cell cycle progression and E2F1-responsive gene expression via association with Host Cell Factor-1 (HCF-1)," The Journal of Biological Chemistry, vol. 288, no. 24, pp. 17532-17543, 2013.

[84] Y. J. Machida, Y. Machida, A. A. Vashisht, J. A. Wohlschlegel, and A. Dutta, "The deubiquitinating enzyme BAP1 regulates cell growth via interaction with HCF-1," The Journal of Biological Chemistry, vol. 284, no. 49, pp. 34179-34188, 2009.

[85] C. Chapat and L. Corbo, "Novel roles of the CCR4-NOT complex," Wiley Interdisciplinary Reviews: RNA, vol. 5, no. 6, pp. 883-901, 2014.

[86] A. S. Prabowo, J. van Scheppingen, A. M. Iyer et al., "Differential expression and clinical significance of three inflammation-related microRNAs in gangliogliomas," Journal of Neuroinflammation, vol. 12, no. 1, article 97, 2015.

[87] T. S. Elton, H. Selemon, S. M. Elton, and N. L. Parinandi, "Regulation of the MIR155 host gene in physiological and pathological processes," Gene, vol. 532, no. 1, pp. 1-12, 2013.

[88] O. Loss and F. A. Stephenson, "Localization of the kinesin adaptor proteins trafficking kinesin proteins 1 and 2 in primary cultures of hippocampal pyramidal and cortical neurons," Journal of Neuroscience Research, vol. 93, no. 7, pp. 1056-1066, 2015.

[89] R. A. Isidro, M. L. Cruz, A. A. Isidro et al., "Immunohistochemical expression of SP-NK-1R-EGFR pathway and VDR in colonic inflammation and neoplasia," World Journal of Gastroenterology, vol. 21, no. 6, pp. 1749-1758, 2015.

[90] P.-J. Chen, C. Huang, X.-M. Meng, and J. Li, "Epigenetic modifications by histone deacetylases: biological implications and therapeutic potential in liver fibrosis," Biochimie, vol. 116, pp. 61-69, 2015.

[91] C. Felice, A. Lewis, A. Armuzzi, J. O. Lindsay, and A. Silver, "Review article: selective histone deacetylase isoforms as potential therapeutic targets in inflammatory bowel diseases," Alimentary Pharmacology and Therapeutics, vol. 41, no. 1, pp. 26-38, 2015.

[92] I. A. Lee, A. Kamba, D. Low, and E. Mizoguchi, "Novel methylxanthine derivative-mediated anti-inflammatory effects in inflammatory bowel disease," World Journal of Gastroenterology, vol. 20, no. 5, pp. 1127-1138, 2014.

[93] S. Garcia-Maurino, A. Alcaide, and C. Dominguez, "Pharmacological control of autophagy: therapeutic perspectives in inflammatory bowel disease and colorectal cancer," Current Pharmaceutical Design, vol. 18, no. 26, pp. 3853-3873, 2012.

[94] H. Lee, S. Cha, M.-S. Lee, G. J. Cho, W. S. Choi, and K. Suk, "Role of antiproliferative B cell translocation gene-1 as an apoptotic sensitizer in activation-induced cell death of brain microglia," Journal of Immunology, vol. 171, no. 11, pp. 5802-5811, 2003.

[95] C. A. Hamm and F. F. Costa, "Epigenomes as therapeutic targets," Pharmacology & Therapeutics, vol. 151, pp. 72-86, 2015.

[96] H. Liu, A. T. Cao, T. Feng et al., "TGF-[beta] converts Th1 cells into Th17 cells through stimulation of Runx1 expression," European Journal of Immunology, vol. 45, no. 4, pp. 1010-1018, 2015.

[97] W. F. Wong, K. Kohu, A. Nakamura et al., "Runx1 deficiency in [CD4.sup.+] T cells causes fatal autoimmune inflammatory lung disease due to spontaneous hyperactivation of cells," Journal of Immunology, vol. 188, no. 11, pp. 5408-5420, 2012.

[98] G. P Christophi, R. Rong, P G. Holtzapple, P. T. Massa, and S. K. Landas, "Immune markers and differential signaling networks in ulcerative colitis and Crohn's disease," Inflammatory Bowel Diseases, vol. 18, no. 12, pp. 2342-2356, 2012.

[99] D. Salisbury and U. Bronas, "Inflammation and immune system contribution to the etiology of atherosclerosis: mechanisms and methods of assessment," Nursing Research, vol. 63, no. 5, pp. 375-385, 2014.

[100] D. Dunkin, S. Mehandru, and J.-F. Colombel, "Immune cell therapy in IBD," Digestive Diseases, vol. 32, supplement 1, pp. 61-66, 2014.

[101] J. B. Hurov, T. S. Stappenbeck, C. M. Zmasek et al., "Immune system dysfunction and autoimmune disease in mice lacking Emk (Par-1) protein kinase," Molecular and Cellular Biology, vol. 21, no. 9, pp. 3206-3219, 2001.

[102] Y. Miwa, S. Yazaki, M. Iwamoto et al., "Functional difference between membrane-bound and soluble human thrombomodulin," Transplantation, vol. 99, no. 4, pp. 702-709, 2015.

[103] M. C. Soult, Y. Dobrydneva, K. H. Wahab, L. D. Britt, and C. J. Sullivan, "Outer membrane vesicles alter inflammation and coagulation mediators," Journal of Surgical Research, vol. 192, no. 1, pp. 134-142, 2014.

[104] J. Pekow, U. Dougherty, Y. Huang et al., "Gene signature distinguishes patients with chronic ulcerative colitis harboring remote neoplastic lesions," Inflammatory Bowel Diseases, vol. 19, no. 3, pp. 461-470, 2013.

[105] M. Lettau, M. Paulsen, D. Kabelitz, and O. Janssen, "FasL expression and reverse signalling," Results and Problems in Cell Differentiation, vol. 49, pp. 49-61, 2009.

[106] T. J. Slebioda and Z. Kmiec, "Tumour necrosis factor superfamily members in the pathogenesis of inflammatory bowel disease," Mediators of Inflammation, vol. 2014, Article ID 325129, 15 pages, 2014.

[107] W. Ben Aleya, I. Sfar, L. Mouelhi et al., "Association of Fas/Apo1 gene promoter (-670 A/G) polymorphism in Tunisian patients with IBD," World Journal of Gastroenterology, vol. 15, no. 29, pp. 3643-3648, 2009.

[108] C. V. Antunes, A. E. Hallack Neto, C. R. Nascimento et al., "Anemia in Inflammatory bowel disease outpatients: prevalence, risk factors, and etiology," BioMed Research International, vol. 2015, Article ID 728925, 7 pages, 2015.

[109] J. W. Crott, Z. Liu, M. K. Keyes et al., "Moderate folate depletion modulates the expression of selected genes involved in cell cycle, intracellular signaling and folate uptake in human colonic epithelial cell lines," Journal of Nutritional Biochemistry, vol. 19, no. 5, pp. 328-335, 2008.

[110] L. Pochini, M. Scalise, M. Galluccio, G. Pani, K. A. Siminovitch, and C. Indiveri, "The human OCTN1 (SLC22A4) reconstituted in liposomes catalyzes acetylcholine transport which is defective in the mutant L503F associated to the Crohn's disease," Biochimica et Biophysica Acta--Biomembranes, vol. 1818, no. 3, pp. 559-565, 2012.

[111] K. Repnik and U. Potocnik, "Haplotype in the IBD5 region is associated with refractory Crohn's disease in Slovenian patients and modulates expression of the SLC22A5 gene," Journal of Gastroenterology, vol. 46, no. 9, pp. 1081-1091, 2011.

[112] P Sarlos, D. Varszegi, V. Csongei et al., "Susceptibility to ulcerative colitis in Hungarian patients determined by gene-gene interactions," World Journal of Gastroenterology, vol. 20, no. 1, pp. 219-227, 2014.

[113] M. E. Gove, D. H. Rhodes, M. Pini et al., "Role of leptin receptor-induced STAT3 signaling in modulation of intestinal and hepatic inflammation in mice," Journal of Leukocyte Biology, vol. 85, no. 3, pp. 491-496, 2008.

[114] C. Benezech, S. Nayar, B. A. Finney et al., "CLEC-2 is required for development and maintenance of lymph nodes," Blood, vol. 123, no. 20, pp. 3200-3207, 2014.

[115] J. Shin, I. Yoon, J. Lim et al., "CD4+VEGFR1HIGH T cell as a novel Treg subset regulates inflammatory bowel disease in lymphopenic mice," Cellular and Molecular Immunology, vol. 12, no. 5, pp. 592-603, 2015.

[116] S. Maeda, K. Ohno, A. Fujiwara-Igarashi, K. Uchida, and H. Tsujimoto, "Changes in Foxp3-positive regulatory T cell number in the intestine of dogs with idiopathic inflammatory bowel disease and intestinal lymphoma," Veterinary Pathology, vol. 53, no. 1, pp. 102-112, 2016.

[117] A. Di Sabatino, L. Rovedatti, R. Kaur et al., "Targeting gut T cell Ca2+ release-activated Ca2+ channels inhibits T cell cytokine production and T-box transcription factor T-bet in inflammatory bowel disease," Journal of Immunology, vol. 183, no. 5, pp. 3454-3462, 2009.

[118] I. Gromova, P Gromov, N. Kroman et al., "Immunoexpression analysis and prognostic value of BLCAP in breast cancer," PLoS ONE, vol. 7, no. 9, Article ID e45967, 2012.

[119] A. K.-F. Lo, K.-W. Lo, C.-W. Ko, L. S. Young, and C. W. Dawson, "Inhibition of the LKB1-AMPK pathway by the Epstein-Barr virus-encoded LMP1 promotes proliferation and transformation of human nasopharyngeal epithelial cells," The Journal of Pathology, vol. 230, no. 3, pp. 336-346, 2013.

[120] L. Verstrepen and R. Beyaert, "Receptor proximal kinases in NF-[kappa]B signaling as potential therapeutic targets in cancer and inflammation," Biochemical Pharmacology, vol. 92, no. 4, pp. 519-529, 2014.

[121] Z. Liu, W. Zhang, M. Zhang, H. Zhu, C. Moriasi, and M. Zou, "Liver kinase B1 suppresses lipopolysaccharide-induced nuclear factor [kappa]B (NF-[kappa]B) activation in macrophages," Journal of Biological Chemistry, vol. 290, no. 4, pp. 2312-2320, 2015.

[122] F Q. Calvo, M. Fillet, D. De Seny et al., "Biomarker discovery in asthma-related inflammation and remodeling," Proteomics, vol. 9, no. 8, pp. 2163-2170, 2009.

[123] L. Beaugerie, "IBD and increased risk of cancer: what is the reality?" La Revue de l'Infirmiere, vol. 63, no. 199, p. 28, 2014.

[124] M. C. Montesinos, A. Desai-Merchant, and B. N. Cronstein, "Promotion of wound healing by an agonist of adenosine A2A receptor is dependent on tissue plasminogen activator," Inflammation, vol. 38, no. 6, pp. 2036-2041, 2015.

[125] Z. Shaghaghi, M. Bonyadi, M. H. Somi, and M. Khoshbaten, "Association of plasminogen activator inhibitor-1 gene polymorphism with inflammatory bowel disease in Iranian Azeri Turkish patients," Saudi Journal of Gastroenterology, vol. 20, no. 1, pp. 54-58, 2014.

[126] I. E. Koutroubakis, A. Sfiridaki, G. Tsiolakidou, C. Coucoutsi, A. Theodoropoulou, and E. A. Kouroumalis, "Plasma thrombinactivatable fibrinolysis inhibitor and plasminogen activator inhibitor-1 levels in inflammatory bowel disease," European Journal of Gastroenterology and Hepatology, vol. 20, no. 9, pp. 912-916, 2008.

[127] E. Balboa-Beltran, M. J. Fernandez-Seara, A. Perez-Munuzuri et al., "A novel stop mutation in the vascular endothelial growth factor-C gene (VEGFC) results in Milroy-like disease," Journal of Medical Genetics, vol. 51, no. 7, pp. 475-478, 2014.

[128] L. Le Guen, T. Karpanen, D. Schulte et al., "Ccbe1 regulates Vegfc-mediated induction of Vegfr3 signaling during embryonic lymphangiogenesis," Development, vol. 141, no. 6, pp. 1239-1249, 2014.

[129] C. Tacconi, C. Correale, A. Gandelli et al., "Vascular endothelial growth factor C disrupts the endothelial lymphatic barrier to promote colorectal cancer invasion," Gastroenterology, vol. 148, no. 7, pp. 1438.e8-1451.e8, 2015.

[130] S. D'Alessio, C. Correale, C. Tacconi et al., "VEGF-C-dependent stimulation of lymphatic function ameliorates experimental inflammatory bowel disease," Journal of Clinical Investigation, vol. 124, no. 9, pp. 3863-3878, 2014.

[131] C. S. Jorgensen, G. Levantino, G. Houen et al., "Determination of autoantibodies to annexin XI in systemic autoimmune diseases," Lupus, vol. 9, no. 7, pp. 515-520, 2000.

[132] M. Heuser, D. B. Yap, M. Leung et al., "Loss of MII5 results in pleiotropic hematopoietic defects, reduced neutrophil immune function, and extreme sensitivity to DNA demethylation," Blood, vol. 113, no. 7, pp. 1432-1443, 2009.

[133] S. Schulz, G. Chachami, L. Kozaczkiewicz et al., "Ubiquitin-specific protease-like 1 (USPL1) is a SUMO isopeptidase with essential, non-catalytic functions," EMBO Reports, vol. 13, no. 10, pp. 930-938, 2012.

[134] J. Parssinen, E.-L. Alarmo, S. Khan, R. Karhu, M. Vihinen, and A. Kallioniemi, "Identification of differentially expressed genes after PPM1D silencing in breast cancer," Cancer Letters, vol. 259, no. 1, pp. 61-70, 2008.

[135] J. Z. Liu, S. van Sommeren, H. Huang et al., "Association analyses identify 38 susceptibility loci for inflammatory bowel disease and highlight shared genetic risk across populations," Nature Genetics, vol. 47, no. 9, pp. 979-986, 2015.

[136] J. Pinero, N. Queralt-Rosinach, A. Bravo et al., "DisGeNET: a discovery platform for the dynamical exploration of human diseases and their genes," Database, vol. 2015, Article ID bav028, 2015.

Fei Yuan, (1) Yu-Hang Zhang, (2) Xiang-Yin Kong, (2) and Yu-Dong Cai (3)

(1) Department of Science & Technology, Binzhou Medical University Hospital, Binzhou 256603, Shandong, China

(2) Institute of Health Sciences, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200031, China

(3) School of Life Sciences, Shanghai University, Shanghai 200444, China

Correspondence should be addressed to Fei Yuan;

Received 26 August 2016; Accepted 11 January 2017; Published 14 February 2017

Academic Editor: Mikihiro Fujiya

Caption: Figure 1: IFS curve. y-axis represents the total prediction accuracy, and x-axis represents the number of features participating in the classification.

Caption: Figure 2: The graph consisting of 190 shortest paths connecting any two genes in the optimal gene set. The yellow diamonds represent genes in the optimal gene set. The blue diamonds represent shortest-path genes. The numbers on the edges represent the edge weights in the network.
Table 1: Twenty-one important genes for IBD obtained using the
mRMR and IFS methods.

GO term/KEGG    Description        Rank (a)    Gene symbol
pathway ID

hsa04660        T cell receptor       3           CD247
                signaling             19           CD4

GO: 0001775     Cell activation       5            PF4
                                      1          ZNF207
                                      8           EGR3
                                      4           SLTM

GO: 0045449     Regulation of         14          CNOT8
                transcription         13      TH1L (NELFCD)

                                      9           HMGB1
                                      12          UBE2I

GO: 0007243     Protein kinase        2           MARK2
                cascade               6           FOLR1

GO: 0031226     Intrinsic             18         SLC22A4
                to plasma

                                      10         LEPROT
                                      16         CLEC1B

GO: 0006915     Apoptosis             21          RHOT2
                                      7           BLCAP
                                      20         ANXA11
                                      15           OGT

Non-grouped                           17          USPL1
genes                                 11        HIST1H2AC

GO term/KEGG                     Description
pathway ID

hsa04660                       CD247 molecule
                                CD4 molecule

GO: 0001775                   Platelet factor 4
                           Zinc finger protein 207
                           Early growth response 3
                     SAFB-like, transcription modulator

GO: 0045449       CCR4-NOT transcription complex, subunit 8
                  Negative elongation factor complex member
                          High mobility group box 1
                      Ubiquitin-conjugating enzyme E2I

GO: 0007243     MAP/microtubule affinity-regulating kinase 2
                          Folate receptor 1 (adult)

GO: 0031226           Solute carrier family 22 (organic

                  cation/zwitterion transporter), member 4

                   Leptin receptor overlapping transcript
                   C-type lectin domain family 1, member B

GO: 0006915             ras homolog family member T2
                      Bladder cancer associated protein
                                 Annexin A11
                    O-linked N-acetylglucosamine (GlcNAc)

Non-grouped          Ubiquitin specific peptidase like 1
genes                      Histone cluster 1, H2ac

(a): this column indicates the ranks of related features in the mRMR
feature list.

Table 2: Twenty candidate genes obtained by SP approach.

GO term/       Description     Gene        Ensembl ID
KEGG                           symbol

                                FOS      ENSP00000306245

hsa04660         T cell        PLCG1     ENSP00000244007
               signaling        LCK      ENSP00000337825
                               ZAP70     ENSP00000264972

                               YWHAZ     ENSP00000309503
GO: 0001775       Cell
                                TLR4     ENSP00000363089
                                 F2      ENSP00000308541

                               HCFC1     ENSP00000309555
                               CNOT1     ENSP00000320949

GO:           Regulation of    CNOT4     ENSP00000354673
45449         transcription
                               TRAK1     ENSP00000328998

                               HDAC1     ENSP00000362649
                                BTG1     ENSP00000256015

                               RUNX1     ENSP00000300305

GO:            Intrinsic        THBD     ENSP00000366307
31226           to plasma
                membrane       FAS LG    ENSP00000356694

                               STK11     ENSP00000324856
                               S100A6    ENSP00000357708
GO:            Regulation
42127            of cell      SERPINE1   ENSP00000223095

                               VEGFC     ENSP00000280193

GO term/                 Description               Betweenness

                   FBJ murine osteosarcoma             20
                    viral oncogene homolog
                      Phospholipase C,
hsa04660                   gamma 1                     19
                     LCK proto oncogene,               21
                          Src family
                       tyrosine kinase
                      Zeta-chain (TCR)                 16
                     associated protein
                        kinase 70 kDa

                  3-monooxygenase/tryptophan           19
GO: 0001775       5-monooxygenase activation
                        protein, zeta
                     Toll-like receptor 4              19
               coagulation factor II (thrombin)        19

                     Host cell factor Cl               36
                    CCR4-NOT transcription              6
                      complex, subunit 1
GO:                 CCR4-NOT transcription              6
45449                 complex, subunit 4
                 Trafficking protein, kinesin          19

                          binding 1
                    Histone deacetylase 1              19
                 B-cell translocation gene 1,          13
              Runt-related transcription factor        19

GO:                     Thrombomodulin                 19
31226            Fas ligand (TNF superfamily,          19

                          member 6)

                  Serine/threonine kinase 11           19
               S100 calcium binding protein A6         19
GO:           Serpin peptidase inhibitor, clade
42127          E (nexin, plasminogen activator         18

                 inhibitor type 1), member 1
                 Vascular endothelial growth           19
                           factor C

GO term/      Permutation     Maximum      Most related
KEGG              FDR       interaction       gene in
pathway                        score       the optimal
ID                                           gene set

                 0.035          950             CD4

hsa04660         0.022          927             CD4
                 0.02           999             CD4

                 0.016          999            CD247

                 0.009          962            MARK2
GO: 0001775

                <0.001          970            HMGB1
                <0.001          953             PF4

                <0.001          997             OCT
                 0.004          998            CNOT8

GO:              0.008          987            CNOT8
                <0.001          946             OCT

                 0.031          967            UBE2I
                 0.006          981            CNOT8

                <0.001          927           SLC22A4

GO:             <0.001          985             PF4
31226           <0.001          985            EGR3

                <0.001          986            MARK2
                <0.001          987           ANXA11
42127            0.034          933             PF4

                <0.001          919             PF4
COPYRIGHT 2017 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Yuan, Fei; Zhang, Yu-Hang; Kong, Xiang-Yin; Cai, Yu-Dong
Publication:BioMed Research International
Article Type:Report
Date:Jan 1, 2017
Previous Article:Growth Profile and Its Association with Nutrient Intake and Dietary Patterns among Children and Adolescents in Hail Region of Saudi Arabia.
Next Article:Shock, Cardiac Arrest, and Resuscitation.

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters