BioMed Research International

Volume 2014 (2014), Article ID 891945, 8 pages

http://dx.doi.org/10.1155/2014/891945

Research Article

A Graphic Method for Identification of Novel Glioma Related Genes

Yu-Fei Gao,¹ Yang Shu,² Lei Yang,³ Yi-Chun He,¹ Li-Peng Li,¹ GuaHua Huang,⁴ Hai-Peng Li,⁵ and Yang Jiang¹

¹Department of Surgery, China-Japan Union Hospital of Jilin University, Changchun 130033, China

²State Key Laboratory of Medical Genomics, Institute of Health Sciences, Shanghai Jiaotong University School of Medicine and Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200025, China

³Endoscopy Center, China-Japan Union Hospital of Jilin University, Changchun 130033, China

⁴Institute of Systems Biology, Shanghai University, Shanghai 200444, China

⁵CAS-MPG Partner Institute for Computational Biology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200031, China

Received 18 April 2014; Revised 25 May 2014; Accepted 28 May 2014; Published 23 June 2014

Academic Editor: Tao Huang

Copyright © 2014 Yu-Fei Gao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Glioma, as the most common and lethal intracranial tumor, is a serious disease that causes many deaths every year. Good comprehension of the mechanism underlying this disease is very helpful to design effective treatments. However, up to now, the knowledge of this disease is still limited. It is an important step to understand the mechanism underlying this disease by uncovering its related genes. In this study, a graphic method was proposed to identify novel glioma related genes based on known glioma related genes. A weighted graph was constructed according to the protein-protein interaction information retrieved from STRING and the well-known shortest path algorithm was employed to discover novel genes. The following analysis suggests that some of them are related to the biological process of glioma, proving that our method was effective in identifying novel glioma related genes. We hope that the proposed method would be applied to study other diseases and provide useful information to medical workers, thereby designing effective treatments of different diseases.

1. Introduction

Glioma is the most common and lethal intracranial tumor. It always revealed itself as malignant glioma which is usually divided into astrocytoma, oligodendroglioma, and oligoastrocytoma. Besides the classification based on histopathological features, glioma could also be graded on a WHO consensus-derived scale of I to IV by means of the degree of malignancy [1]. Clinically, most of the gliomas are high-grade gliomas (HGG). Glioblastoma (GBM), one of the HGG, accounts for more than half of gliomas [1]. Although the knowledge of glioma, especially HGG, has increased dramatically in recent years, many questions are still waiting for further elucidation. On the other hand, the overall 5-year survival rate of GBM remains less than 5% despite the advances in surgery, radiation, and chemotherapy [2].

In the previous reports, glioma always manifested itself with disordered pathways which regulated proliferation, survival, invasion, and angiogenesis. Among these biological processes, RB and p53 pathways are more inclined to be dysfunctional in GBM, and the disrepair could lead to the destruction of cell cycle by regulating the G1-to-S-phase transition [3, 4]. Furthermore, other pathways such as MAPK, PI3K/PTEN/AKT, and NF-κB pathway are also overactivated in glioma and contribute to the uncontrollable cellular proliferation [5–7]. As we know, the tumorigenesis of glioma is a complicated process which involved intricate pathways beyond the above ones. To widely understand the mechanism underlying this disease, identification of its related genes and uncovering the relationship of them and the biological process of glioma are very important. However, it is time-consuming and expensive to identify novel glioma related genes by conventional experiments. On the other hand, encouraged by the successful application of computational methods to deal with various biological problems such as drug design [8–13] and analysis of complicated biological pathway [14–17], computational methods may address this problem and provide some useful information for investigators.

In this study, we proposed a graphic method and attempted to apply this method to discover novel glioma related genes. The current known glioma related genes collected from various sources were the firsthand information. Based on these genes, some new discovered genes were obtained by the well-known shortest path algorithm. Furthermore, a permutation test was conducted to exclude false positives among them. The analysis of the final remaining genes suggests that some of them had direct or indirect relationship to the biological process of glioma, indicating that this method was effective and may give new insight to study other diseases.

2. Materials and Methods

2.1. Materials

The current known glioma related genes were retrieved from the following sources. All the 11 data sheets listed on the web page of COSMIC (http://cancer.sanger.ac.uk/cancergenome/projects/census/) were downloaded, from which we obtained 18 glioma related genes; search for human diseases in UniProt (http://www.uniprot.org/) with keywords “human glioma oncogene” and “human glioma suppressor gene,” thereby obtaining 49 and 32 genes (only reviewed genes were selected), respectively; select “Literature Search” in TSGene (http://bioinfo.mc.vanderbilt.edu/TSGene/search.cgi) and input “glioma” as keyword, obtaining 7 genes. After collecting all of the genes mentioned above, we finally obtained 77 glioma related genes, which were available in Supplementary Material I available online at http://dx.doi.org/10.1155/2014/891945.

2.2. Construction of a Weighted Graph from Protein-Protein Interactions

Protein-protein interaction (PPI) is useful information for investigating various biological problems [18–22]. Many computational methods were proposed based on the fact that proteins that can interact with each other always share similar functions. Since the known glioma related genes must have some common features related to glioma, it is reasonable to discover novel glioma related genes based on protein-protein interaction and known glioma related genes. The data concerning protein-protein interactions was downloaded from STRING (Search Tool for the Retrieval of Interacting Genes/Proteins, http://string.embl.de/) [23], a large database containing known and predicted protein interactions which are derived from genomic context, high-throughput experiments, (Conserved) Coexpression and Previous knowledge. In the obtained file, we extracted all protein-protein interactions of human. Each obtained interaction consists of two proteins and a score with range between 150 and 999, which can quantify the likelihood that an interaction may occur. For later formulation, let denote the score of the interaction between two proteins and . The constructed graph took proteins occurring at least one protein-protein interaction of human as nodes, while two nodes were adjacent if and only if the score of the interaction between the corresponding proteins was greater than zero. The obtained graph consisted of 18,600 nodes and 1,640,707 edges. Furthermore, to correctly reflect the strength of the interaction, each edge with nodes and was labeled by a weight, which can be computed by

where and were two corresponding proteins of nodes and , respectively.

2.3. Selection of Candidate Genes

It is obvious that the glioma related genes must have some common features which are related to glioma. On the other hand, as mentioned in Section 2.2, two proteins that can interact with each other, that is, they are adjacent in the constructed weighted graph, always share common features. The idea of our method was based on these facts. To clearly elaborate the idea of the method, we constructed a simple weighted graph which is shown in Figure 1, because the original graph was too large to exhibit in the paper. It is easy to observe from Figure 1 that the shortest path connecting and contains , and as the inner nodes. Based on the weights of edges on this path, we can obtain that genes and can share common functions with high probability, because the confidence score of the interaction between and is very high, which is . The similar results also hold for and , and , and . If genes and are two known glioma related genes, genes , , and are actual glioma related genes with high likelihood. In view of this, shortest paths between any pair of known glioma related genes obtained by Dijkstra’s algorithm, the most famous shortest path algorithm proposed by Dijkstra in 1956 [24], are useful information for further investigation.

Figure 1: A simple example of the weighted graph.

After obtaining the shortest paths connecting any pair of known glioma related genes, it can be seen that some nodes/genes occurred in many paths, while most of nodes/genes in the graph were not contained in any path. Thus, for each node/gene in the graph, we counted the number of paths containing the node/gene, termed as betweenness which is defined as the number of shortest paths containing the node/gene as an inner node. The concept of betweenness has been employed in some studies of natural and man-made networks [25–29]. In fact, the betweenness of some node/gene reflects the direct and indirect relationship of the gene and known glioma related genes. Thus, the likelihood of genes with high betweenness to be related to glioma was higher than those with low betweenness. In view of this, we selected genes with betweenness greater than 0 as the candidate genes, which may be the novel glioma related genes with high probability. It is necessary to point out that the known glioma related genes were not included in the set of candidate genes.

2.4. Filtering Candidate Genes by Permutation Test

As described in Section 2.3, some candidate genes can be obtained by researching the shortest paths connecting any two known glioma related genes. However, some of them may be false positives, because some nodes/genes may easily receive a high betweenness due to their location in the weighted graph even if we randomly select genes in STRING as the known glioma related genes. To exclude these false discoveries, a permutation test should be executed as follows.

(I)1,000 node/gene sets, denoted by , were randomly selected in the weighted graph such that each of them had the same size of known glioma related gene set.(II)For each candidate gene discovered in Section 2.3, calculate its betweenness on each set .(III)Calculate the permutation FDR of each candidate gene by

where was set to be 1 if the betweenness of on was larger than that of on the known glioma related gene set; otherwise, it was set to be 0.

Obviously, small permutation FDR of one candidate gene implies that it is the true positive with high probability.

3. Results and Discussions

3.1. Candidate Genes

For the 77 genes mentioned in Section 2.1, we searched the shortest path connecting any two of them by Dijkstra algorithm. After calculating the betweenness of each node/gene in the weighted graph, 215 candidate genes with betweenness larger than zero were obtained. These 215 genes and their betweenness were available in Supplementary Material II. To exclude the false positives, the permutation test was conducted as described in Section 2.4. By (2), we can calculate the permutation FDR of each candidate gene. These values were also provided in Supplementary Material II. Since the likelihood of gene with small permutation FDR to be the actual glioma related gene is high, we set the threshold to be 0.05, that is, selecting genes with permutation FDRs lower than 0.05 among 215 candidate genes, thereby obtaining 67 genes, listed in Table 1. These genes were deemed to have strong relationship with glioma and further discussions were based on these genes.

Table 1: Candidate genes with permutation FDR lower than 0.05.

Ensemble ID of candidate gene
Gene name
Betweenness
Permutation FDR

ENSP00000227638
PANX1
72
0

ENSP00000235310
MAD2L2
72
0

ENSP00000245323
EFNB2
335
0

ENSP00000258428
REV1
72
0

ENSP00000265727
ADAM22
72
0

ENSP00000281821
EPHA4
339
0

ENSP00000293831
EIF4A1
72
0

ENSP00000302719
KCNAB3
72
0

ENSP00000312697
DMAP1
88
0

ENSP00000329797
CADM1
72
0

ENSP00000335434
WDR20
72
0

ENSP00000341138
EPB41L3
72
0

ENSP00000351697
REV3L
72
0

ENSP00000354778
CNTNAP2
72
0

ENSP00000356150
MDM4
72
0

ENSP00000357177
ARHGEF11
142
0

ENSP00000361366
SFTPD
72
0

ENSP00000369218
RBM17
72
0

ENSP00000370119
SMN2
72
0

ENSP00000245304
RAP2A
72
0.001

ENSP00000275815
EPHA1
72
0.001

ENSP00000328511
KCNA4
72
0.001

ENSP00000377446
SUCLG1
210
0.001

ENSP00000399511
TNIK
72
0.001

ENSP00000229595
ASF1A
72
0.002

ENSP00000252699
ACTN4
72
0.002

ENSP00000263208
HIRA
72
0.002

ENSP00000263923
KDR
104
0.002

ENSP00000304169
PITX2
210
0.002

ENSP00000330633
CNTN2
72
0.002

ENSP00000276072
TAF1
72
0.003

ENSP00000295600
MITF
138
0.003

ENSP00000360157
FOXD3
6
0.003

ENSP00000264010
CTCF
72
0.004

ENSP00000271628
SF3B4
75
0.004

ENSP00000350941
SRC
421
0.004

ENSP00000361125
VEGFA
164
0.004

ENSP00000226091
EFNB3
67
0.006

ENSP00000358309
EPHA7
2
0.006

ENSP00000358918
SUFU
72
0.007

ENSP00000316879
EIF4G1
72
0.011

ENSP00000260653
SIX3
2
0.012

ENSP00000352516
DNMT1
90
0.016

ENSP00000358716
DDX20
72
0.017

ENSP00000357393
EFNA3
50
0.018

ENSP00000341680
DTNBP1
66
0.02

ENSP00000344456
CTNNB1
607
0.02

ENSP00000297904
FIGF
2
0.021

ENSP00000265165
LEF1
134
0.022

ENSP00000347948
TNFRSF14
68
0.023

ENSP00000288986
NCK1
78
0.027

ENSP00000261937
FLT4
2
0.029

ENSP00000333919
BTLA
68
0.03

ENSP00000332549
GRIN2A
58
0.032

ENSP00000376765
PIAS3
4
0.033

ENSP00000361818
SDC4
1
0.035

ENSP00000386165
CEBPD
38
0.035

ENSP00000348307
SIRPA
34
0.036

ENSP00000344666
NF2
1
0.037

ENSP00000219255
PARD6A
72
0.038

ENSP00000204961
EFNB1
5
0.039

ENSP00000172229
NGFR
72
0.043

ENSP00000344115
CDH5
24
0.043

ENSP00000405041
POU5F1
6
0.045

ENSP00000360532
CDC5L
6
0.046

ENSP00000295897
ALB
72
0.047

ENSP00000340944
PTPN11
112
0.047

3.2. Analysis of Enriched KEGG Pathways of Candidate Genes

As mentioned in Section 3.1, 67 candidate genes were obtained. To analyze the relationship of them and glioma, we employed DAVID (Database for Annotation, Visualization and Integrated Discovery) [30], a functional annotation tool to understand biological meaning behind large list of genes. 67 candidate genes comprised the input gene list of DAVID, thereby obtaining 9 KEGG pathways that were enriched by these 67 candidate genes. The detailed output of DAVID for KEGG pathway enrichment analysis was available as Supplementary Material III.

The top 5 pathways have the value less than 0.05 which were discussed below. The most enriched pathway is hsa04360: axon guidance (“count” = 8). Among the 8 genes, 7 ephrins-related genes are enriched whose corresponding proteins include 4 members (EFNA3, EFNB1, EFNB2, and EFNB3) of the ephrins family and 3 members (EPHA1, EPHA4, and EPHA7) of the ephrins receptor subfamily. Eph receptor tyrosine kinases (Ephs) and ephrins (EPH) could navigate cells by controlling cell-cell adhesion and segregation [31]. In other words, with the function of axon guidance Ephs/EPH could regulate the invasion, neoangiogenesis, and metastasis of gliomas [32, 33]. Ding and his colleagues have identified several somatic mutations of Ephs especially EphA7 in lung cancer [34]. Although the close connection between Ephs/EPH and cancer has been reported, its pathogenic mechanism in gliomas is still unknown. The second pathway is hsa04510: focal adhesion (“count” = 7). As we know, infiltration of tumor cells and angiogenesis are critical for the growth of tumor. Zagzag et al. reported that focal adhesion kinase (FAK), highly associated with these biological processes, plays an important role in tumorigenesis of gliomas via enhancing the ability of infiltration and angiogenesis [35]. The FAK-related genes, like CTNNB1, VEGFA, KDR, and FLT4 enriched in this pathway, are always mutated or aberrantly expressed in various types of cancers [36, 37]. The third pathway is hsa04530: tight junction (“count” = 5). In the brain, the expression of the tight junction proteins is important for blood-brain tumor-barrier (BTB) permeability. Hence destruction of the tight junction could facilitate the development of gliomas by increasing BTB permeability [38, 39]. The fourth pathway is hsa04520: adherens junction (“count” = 4). Adherens junction is reported to be disordered in the glioblastoma and to affect the invasive behavior of GBM [40, 41]. The last significantly enriched pathway is hsa05200: pathways in cancer (“count” = 7). The result shows that a common mechanism is shared by the gliomas and other types of cancers. Although these significant enriched pathways have been reported to be related to gliomas more or less, our results might expand the avenues to explore new mechanisms in the tumorigenesis of gliomas.

3.3. Analysis of Enriched GO Terms Candidate Genes

In addition to KEGG pathway enrichment analysis, DAVID also provided the GO terms enrichment analysis of the 67 candidate genes, which were available in Supplementary Material IV.

It can be seen that 227 GO terms were enriched by these 67 genes. Top 10 Go terms sorted by value are investigated and discussed as below. Among the top 10, 4 GO terms are biological process (BP) which included GO: 0007169: transmembrane receptor protein tyrosine kinase signaling pathway, GO: 0007167: enzyme linked receptor protein signaling pathway, GO: 0042127: regulation of cell proliferation, and GO: 0000904: cell morphogenesis involved in differentiation. From the results, we found that all these processes are connected with receptor-dependent signaling pathways. The cancer genome atlas (TCGA) group has revealed that the receptor tyrosine kinase (RTK) pathway was deregulated in 88% of the patients with glioblastoma [42]. After deregulation of the RTK, its downstream genes could function uncontrollably in the cellular proliferation and morphogenesis which are very pivotal for the growth of gliomas. In the top 10 GO terms, we also find 4 molecular function (MF) GO terms: GO: 0004714: transmembrane receptor protein tyrosine kinase activity, GO: 0005003: ephrin receptor activity, GO: 0046875: ephrin receptor binding, and GO: 0004713: protein tyrosine kinase activity. The MF classification also suggests the importance of RTK signaling pathways especially the Ephs/EPH pathway in the tumorigenesis of gliomas. As the previous reports, RTK pathways could regulate cell proliferation and migration which were indispensable for the development of gliomas [42, 43]. Besides BP and MF GO terms, 2 cellular component (CC) GO terms are also enriched in the top 10 GO terms: GO: 0044459: plasma membrane part and GO: 0005887: integral to plasma membrane. As we know, the transformation of cell membrane is necessary for the migration and invasion process during tumorigenesis of gliomas. Our results pave the way for understanding potential pathogenic mechanism of gliomas.

3.4. Analysis of Some Candidate Genes

Among the 67 genes, several genes are intriguing which may play pivotal role in tumorigenesis of glioma. This section gave the detailed discussion of some candidate genes.

VEGFA, also known as the vascular endothelial growth factor (VEGF), is a member of a large family of growth factors that also includes VEGFB, VEGFC, VEGFD, and placental growth factor (PLGF). VEGF is the only mitogen that specifically acts on endothelial cells and also a tumor angiogenesis factor in human glioma in vivo [44]. Knizetova et al. have demonstrated that the autocrine VEGF signaling is mediated via VEGFR2 (KDR), another gene in our list. They found that blockade of VEGFR2 would abrogate the VEFG-mediated enhancement of astrocytoma cell growth and viability [37]. In the in vivo level, Millauer et al. found that disrepair of VEGFR2/VEGF system in angiogenesis could prevent tumor growth in nude mice [45]. Another VEGF receptor found in our list is VEGFR3 (also known as FLT4). In contrast to VEGFR1/2, VEGFR3 does not bind VEGFA and mainly functions in lymphangiogenesis as a receptor of VEGFC and VEGFD [46, 47]. Jenny et al. reported that VEFGR3 was expressed in some tumor types such as haemangioblastoma and glioblastoma, despite their lack of lymphatic vessels [48]. Although the roles of VEGF signaling pathway in the tumorigenesis of glioma have been well studied, new findings have been explored in succession recently.

CTNNB1, with more famous name of catenin beta 1, encodes β-catenin protein which plays important roles in cellular morphogenesis, differentiation, and proliferation via regulating the Wnt signaling [49]. Yano et al. induced glioma in rat using N-ethyl-N-nitrosourea display aberrant nuclear accumulation of β-catenin in contrast to normal brains [50]. Moreover, Pu et al. found that the downregulation of β-catenin by siRNA could suppress malignant glioma cell growth [36]. To elucidate the connection between β-catenin and glioma, Liu and his colleagues performed a systemic research. They found a higher expression level of β-catenin in astrocytic glioma patients with high grade in comparison with the normal controls. Furthermore, they also illustrated that the overexpression of β-catenin may be an important contributing factor to glioma progression by means of facilitating proliferation and inhibiting apoptosis [51]. After Wnt pathway is activated, β-catenin accumulates and enters the nucleus where it can act as a coactivator for TCF/LEF-mediated transcription [52]. As the downstream of beta-catenin, LEF1, another gene in our list, also plays a crucial role in the Wnt signaling pathway. LEF1, with full name of lymphoid enhancer-binding factor 1, tends to be mutated in the tumors. Liu et al. have investigated that MiR-218 could reduce the invasiveness of glioblastoma cells by targeting LEF1 [53].

SRC, whose corresponding protein is a tyrosine-protein kinase, could play a pivotal role in the regulation of embryonic development and cell growth [54]. Besides functioning in the embryonic development, SRC could also regulate the tumorigenesis of various types of cancers like breast cancer, colon cancer, and brain cancer [55, 56]. Src protein always maintains an inactive state until its Y530 residue is dephosphorylated by protein tyrosine phosphatase-α [57]. Src could also be activated by direct binding of its SH2 and SH3 domains to intracellular proteins or activated tyrosine kinase growth factor receptors [58]. Stettner et al. have found elevated SRC activity in GBM compared with normal brain [59]. On the other hand, Lund et al. found that the infiltration of glioma reduced in Src-deficient mice [60]. It is reported that the increased SRC activity in GBM may be due to increased activation of cell surface growth factor receptors and integrins that activate SRC-family kinases (SFKs) rather than the amplification or mutation of SFK genes [42, 61].

4. Conclusion

In biomedicine and genomics, identification of disease genes is an important topic. This contribution proposed a graphic method to identify novel disease genes and the method was applied to glioma, one kind of cancers. The findings indicate that this method is quite effective. It is hopeful that the contribution can provide help for medical workers to discover effective treatments of glioma and give new insight to study various diseases.

Conflict of Interests

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

Authors’ Contribution

Yu-Fei Gao and Yang Shu contributed equally to this work.

Acknowledgments

This paper is supported by the National Science Foundation of China (81372696), the Natural Science Fund Projects of Jilin province (201215059), the Development of Science and Technology Plan Projects of Jilin province (20100733, 201101074), and SRF for ROCS, SEM (2009-36), Scientific Research Foundation (Jilin Department of Science & Technology; 200705314, 20090175, and 20100733), Scientific Research Foundation (Jilin Department of Health, 2010Z068), and SRF for ROCS (Jilin Department of Human Resource & Social Security, 2012–2014).

References

D. N. Louis, H. Ohgaki, O. D. Wiestler et al., “The 2007 WHO classification of tumours of the central nervous system,” Acta Neuropathologica, vol. 114, pp. 97–109, 2007.

T. A. Dolecek, J. M. Propp, N. E. Stroup, and C. Kruchko, “CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2005–2009,” Neuro-Oncology, vol. 14, supplement 5, pp. v1–v49, 2012.

D. N. Louis, “The p53 gene and protein in human brain tumors,” Journal of Neuropathology and Experimental Neurology, vol. 53, no. 1, pp. 11–21, 1994.

J. W. Henson, B. L. Schnitker, K. M. Correa et al., “The retinoblastoma gene is involved in malignant progression of astrocytomas,” Annals of Neurology, vol. 36, no. 5, pp. 714–721, 1994.

A. Guha, M. M. Feldkamp, N. Lau, G. Boss, and A. Pawson, “Proliferation of human malignant astrocytomas is dependent on Ras activation,” Oncogene, vol. 15, no. 23, pp. 2755–2765, 1997.

C. B. Knobbe, A. Trampe-Kieslich, and G. Reifenberger, “Genetic alteration and expression of the phosphoinositol-3-kinase/Akt pathway genes PIK3CA and PIKE in human glioblastomas,” Neuropathology and Applied Neurobiology, vol. 31, pp. 486–490, 2005.

A. L. Rinkenbaugh and A. S. Baldwin, “Monoallelic deletion of NFKBIA in glioblastoma: when less is more,” Cancer Cell, vol. 19, no. 2, pp. 163–165, 2011.

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.

Y. Yamanishi, M. Araki, A. Gutteridge, W. Honda, and M. Kanehisa, “Prediction of drug-target interaction networks from the integration of chemical and genomic spaces,” Bioinformatics, vol. 24, no. 13, pp. i232–i240, 2008.

Y. Yamanishi, M. Kotera, M. Kanehisa, and S. Goto, “Drug-target interaction prediction from chemical, genomic and pharmacological data in an integrated framework,” Bioinformatics, vol. 26, no. 12, pp. i246–i254, 2010.

L. Chen, J. Lu, X. Luo, and K.-Y. Feng, “Prediction of drug target groups based on chemical-chemical similarities and chemical-chemical/protein connections,” Biochimica et Biophysica Acta: Proteins and Proteomics, vol. 1844, no. 1, pp. 207–213, 2014.

B. Padhy and Y. Gupta, “Drug repositioning: re-investigating existing drugs for new therapeutic indications,” Journal of Postgraduate Medicine, vol. 57, no. 2, pp. 153–160, 2011.

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.

L. Chen, W.-M. Zeng, Y.-D. Cai, and T. Huang, “Prediction of metabolic pathway using graph property, chemical functional group and chemical structural set,” Current Bioinformatics, vol. 8, no. 2, pp. 200–207, 2013.

H.-W. Ma and A.-P. Zeng, “The connectivity structure, giant strong component and centrality of metabolic networks,” Bioinformatics, vol. 19, no. 11, pp. 1423–1430, 2003.

J. M. Dale, L. Popescu, and P. D. Karp, “Machine learning methods for metabolic pathway prediction,” BMC Bioinformatics, vol. 11, article 15, 2010.

L. Chen, B.-Q. Li, and K.-Y. Feng, “Predicting biological functions of protein complexes using graphic and functional features,” Current Bioinformatics, vol. 8, no. 5, pp. 545–551, 2013.

Y. F. Gao, L. Chen, Y. D. Cai, K. Y. Feng, T. Huang, and Y. Jiang, “Predicting metabolic pathways of small molecules and enzymes based on interaction information of chemicals and proteins,” PLoS ONE, vol. 7, no. 9, Article ID e45944, 2012.

R. Sharan, I. Ulitsky, and R. Shamir, “Network-based prediction of protein function,” Molecular Systems Biology, vol. 3, p. 88, 2007.

K. L. Ng, J. S. Ciou, and C. H. Huang, “Prediction of protein functions based on function-function correlation relations,” Computers in Biology and Medicine, vol. 40, no. 3, pp. 300–305, 2010.

P. Bogdanov and A. K. Singh, “Molecular function prediction using neighborhood features,” IEEE-ACM Transactions on Computational Biology and Bioinformatics, vol. 7, no. 2, pp. 208–217, 2010.

P. Gao, Q. P. Wang, L. Chen, and T. Huang, “Prediction of human genes' regulatory functions based on proteinprotein interaction network,” Protein and Peptide Letters, vol. 19, no. 9, pp. 910–916, 2012.

L. J. Jensen, M. Kuhn, M. Stark et al., “STRING 8—a global view on proteins and their functional interactions in 630 organisms,” Nucleic Acids Research, vol. 37, no. 1, pp. D412–D416, 2009.

T. H. Gormen, C. E. Leiserson, R. L. Rivest, and C. Stein, Eds., Introduction to Algorithms, MIT Press, Cambridge, Mass, USA, 1990.

J. Davis and M. Goadrich, “The relationship between precision-recall and ROC curves,” in Proceedings of the 23rd International Conference on Machine Learning (ICML '06), pp. 233–240, New York, NY, USA, June 2006.

R. Bunescu, R. Ge, R. J. Kate et al., “Comparative experiments on learning information extractors for proteins and their interactions,” Artificial Intelligence in Medicine, vol. 33, no. 2, pp. 139–155, 2005.

D. E. Johnson and G. H. I. Wolfgang, “Predicting human safety: screening and computational approaches,” Drug Discovery Today, vol. 5, no. 10, pp. 445–454, 2000.

B.-Q. Li, B. Niu, L. Chen et al., “Identifying chemicals with potential therapy of HIV based on protein-protein and protein-chemical interaction network,” PLoS ONE, vol. 8, no. 6, Article ID e65207, 2013.

J. Zhang, M. Jiang, F. Yuan, K. Y. Feng, Y. D. Cai, et al., “Identification of age-related macular degeneration related genes by applying shortest path algorithm in protein-protein interaction network,” BioMed Research International, vol. 2013, Article ID 523415, 8 pages, 2013.

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.

R. Klein, “Eph/ephrin signaling in morphogenesis, neural development and plasticity,” Current Opinion in Cell Biology, vol. 16, no. 5, pp. 580–589, 2004.

P. W. Janes, S. Adikari, and M. Lackmann, “Eph/ephrin signalling and function in oncogenesis: Lessons from embryonic development,” Current Cancer Drug Targets, vol. 8, no. 6, pp. 473–489, 2008.

E. B. Pasquale, “Eph receptors and ephrins in cancer: bidirectional signalling and beyond,” Nature Reviews Cancer, vol. 10, no. 3, pp. 165–180, 2010.

L. Ding, G. Getz, D. A. Wheeler et al., “Somatic mutations affect key pathways in lung adenocarcinoma,” Nature, vol. 455, no. 7216, pp. 1069–1075, 2008.

D. Zagzag, D. R. Friedlander, B. Margolis et al., “Molecular events implicated in brain tumor angiogenesis and invasion,” Pediatric Neurosurgery, vol. 33, no. 1, pp. 49–55, 2000.

P. Pu, Z. Zhang, C. Kang et al., “Downregulation of Wnt2 and β-catenin by siRNA suppresses malignant glioma cell growth,” Cancer Gene Therapy, vol. 16, no. 4, pp. 351–361, 2009.

P. Knizetova, J. Ehrmann, A. Hlobilkova et al., “Autocrine regulation of glioblastoma cell cycle progression, viability and radioresistance through the VEGF-VEGFR2 (KDR) interplay,” Cell Cycle, vol. 7, no. 16, pp. 2553–2561, 2008.

Y. T. Gu, Y. X. Xue, X. Y. Wei, H. Zhang, and Y. Li, “Calcium-activated potassium channel activator down-regulated the expression of tight junction protein in brain tumor model in rats,” Neuroscience Letters, vol. 493, no. 3, pp. 140–144, 2011.

H. Xie, Y. X. Xue, L. B. Liu, and Y. H. Liu, “Endothelial-monocyte-activating polypeptide II increases blood-tumor barrier permeability by down-regulating the expression levels of tight junction associated proteins,” Brain Research, vol. 1319, pp. 13–20, 2010.

T. Nikuseva-Martic, V. Beros, N. Pecina-Slaus, H. I. Pecina, and F. Bulic-Jakus, “Genetic changes of CDH1, APC, and CTNNB1 found in human brain tumors,” Pathology: Research and Practice, vol. 203, pp. 779–787, 2007.

C. Perego, C. Vanoni, S. Massari et al., “Invasive behaviour of glioblastoma cell lines is associated with altered organisation of the cadherin-catenin adhesion system,” Journal of Cell Science, vol. 115, no. 16, pp. 3331–3340, 2002.

The Cancer Genome Atlas Research Network, “Comprehensive genomic characterization defines human glioblastoma genes and core pathways,” Nature, vol. 455, pp. 1061–1068, 2008.

S. A. Rao, A. Arimappamagan, P. Pandey et al., “miR-219-5p inhibits receptor tyrosine kinase pathway by targeting EGFR in glioblastoma,” PLoS ONE, vol. 8, no. 5, Article ID e63164, 2013.

K. H. Plate, G. Breier, H. A. Weich, and W. Risau, “Vascular endothelial growth factor is a potential tumour angiogenssis factor in human gliomas in vivo,” Nature, vol. 359, no. 6398, pp. 845–848, 1992.

B. Millauer, L. K. Shawver, K. H. Plate, W. Risau, and A. Ullrich, “Glioblastoma growth inhibited in vivo by a dominant-negative Flk-1 mutant,” Nature, vol. 367, no. 6463, pp. 576–579, 1994.

S. Koch, S. Tugues, X. Li, L. Gualandi, and L. Claesson-Welsh, “Signal transduction by vascular endothelial growth factor receptors,” The Biochemical Journal, vol. 437, no. 2, pp. 169–183, 2011.

M. J. Karkkainen and T. V. Petrova, “Vascular endothelial growth factor receptors in the regulation of angiogenesis and lymphangiogenesis,” Oncogene, vol. 19, no. 49, pp. 5598–5605, 2000.

B. Jenny, J. A. Harrison, D. Baetens et al., “Expression and localization of VEGF-C and VEGFR-3 in glioblastomas and haemangioblastomas,” Journal of Pathology, vol. 209, no. 1, pp. 34–43, 2006.

J. N. Anastas and R. T. Moon, “WNT signalling pathways as therapeutic targets in cancer,” Nature Reviews Cancer, vol. 13, no. 1, pp. 11–26, 2013.

H. Yano, A. Hara, J. Shinoda et al., “Immunohistochemical analysis of β-catenin in N-ethyl-N-nitrosourea- induced rat gliomas: implications in regulation of angiogenesis,” Neurological Research, vol. 22, no. 5, pp. 527–532, 2000.

X. Liu, L. Wang, S. Zhao, X. Ji, Y. Luo, and F. Ling, “β-catenin overexpression in malignant glioma and its role in proliferation and apoptosis in glioblastma cells,” Medical Oncology, vol. 28, no. 2, pp. 608–614, 2011.

P. Polakis, “Wnt signaling and cancer,” Genes & Development, vol. 14, no. 15, pp. 1837–1851, 2000.

Y. Liu, W. Yan, W. Zhang et al., “MiR-218 reverses high invasiveness of glioblastoma cells by targeting the oncogenic transcription factor LEF1,” Oncology Reports, vol. 28, no. 3, pp. 1013–1021, 2012.

S. Nada, T. Yagi, H. Takeda et al., “Constitutive activation of Src family kinases in mouse embryos that lack Csk,” Cell, vol. 73, no. 6, pp. 1125–1135, 1993.

S. Zhang, W. C. Huang, L. Zhang, C. Zhang, F. J. Lowery, et al., “SRC family kinases as novel therapeutic targets to treat breast cancer brain metastases,” Cancer Research, vol. 73, pp. 5764–5774, 2013.

J. Chen, A. Elfiky, M. Han, C. Chen, and M. W. Saif, “The Role of Src in colon cancer and its therapeutic implications,” Clinical Colorectal Cancer, vol. 13, pp. 5–13, 2014.

C. Egan, A. Pang, D. Durda, H. Cheng, J. H. Wang, and D. J. Fujita, “Activation of Src in human breast tumor cell lines: elevated levels of phosphotyrosine phosphatase activity that preferentially recognizes the Src carboxy terminal negative regulatory tyrosine 530,” Oncogene, vol. 18, no. 5, pp. 1227–1237, 1999.

T. J. Yeatman, “A renaissance for SRC,” Nature Reviews Cancer, vol. 4, no. 6, pp. 470–480, 2004.

M. R. Stettner, W. Wang, L. B. Nabors et al., “Lyn kinase activity is the predominant cellular Src kinase activity in glioblastoma tumor cells,” Cancer Research, vol. 65, no. 13, pp. 5535–5543, 2005.

C. V. Lund, M. T. Nguyen, G. C. Owens et al., “Reduced glioma infiltration in Src-deficient mice,” Journal of Neuro-Oncology, vol. 78, no. 1, pp. 19–29, 2006.

M. S. Ahluwalia, J. de Groot, W. M. Liu, and C. L. Gladson, “Targeting SRC in glioblastoma tumors and brain metastases: rationale and preclinical studies,” Cancer Letters, vol. 298, no. 2, pp. 139–149, 2010.