Printer Friendly

Identifying and Analyzing Novel Epilepsy-Related Genes Using Random Walk with Restart Algorithm.

1. Introduction

As a classical neurological condition that may suddenly interrupt normal life activities and result in physical injury, epilepsy has been widely used to describe a group of epileptic seizure associated diseases [1, 2]. Epileptic seizures are the typical symptoms of the disease, which is the consequence of a disruption of the electrical communications between neurons [3]. As a common neurological disease, epilepsies are found all over the world and affect people of all ages [4]. Only in America, more than one hundred thousand incident cases are diagnosed as epilepsy per year, seriously threatening their mental and physical health [5,6]. Nowadays, the clinical study on epilepsy has been progressively deepening and typical diagnosis routines of epilepsy have been set up. Considering that epileptic seizures induced by disruption of electrical communications between neurons are the typical symptoms of epilepsy, the long-term paroxysmal epileptic seizures might suggest the initiation and progression of such neurological disease [7, 8]. Generally, it is quite necessary to turn to the doctor for help after someone has more than twice abnormal seizures excluding those with known medical conditions. Considering epilepsy has many subtypes induced by different pathogenic factors, resulting in different complications with similar early seizures symptoms, the diagnosis of epilepsy contributes not only to a confirmation of epilepsy but also to classification of the epilepsy from which the patients suffer into its subtype [9, 10]. The diagnosis of epilepsy can be divided into two main procedures: medical history taking and instrumental inspections. Previous experiments have confirmed that typical family history and multiple medical conditions may lead to neurological abnormalities, which may contribute to the initiation and progression of epilepsy [11-13]. Therefore, the first step of the diagnosis of epilepsy is to inquire about the medical history of the patients and their respective family. However, a definitive diagnosis of epilepsy is performed by the following instrumental inspections. As we have mentioned above, epilepsy is referred to as a group of neurological diseases induced by abnormal electrical communication between neurons [4]. Therefore, the measurement of electrical impulses in brain by an electroencephalogram (EEG) test has been regarded as one of the golden standards for epilepsy diagnosis [14, 15]. Apart from EGG, magnetic resonance spectroscopy (MRS), positron emission tomography (PET), and magnetic resonance imaging (MRI) have also been widely used to diagnose epilepsy [16-18].

Although great progresses have been made to diagnose epilepsy, the therapeutic methods to treat epilepsy are still quite limited and they mainly contribute to the symptomatic relief. There are two main functional methods to relieve the seizure symptoms of epilepsy: certain nutrient intakes and vagus nerve stimulations by surgery [19-22]. The therapeutic nutrients that have been confirmed to contribute to the relief of epilepsy include folic acid, melatonin, and vitamins (large doses) [19]. However, such a treatment cannot provide a permanent cure but tries to temporarily relieve the symptoms. The vagus nerve stimulation, as the most effective treatment for epilepsy, necessitates implantation of a pacemaker-like device in the patient's body to stimulate the vagus nerve, relieving the symptoms with few side-effects [23]. These treatments do not have an in-depth consideration of the pathogenic factors that cause the disease but mainly concentrate on the relief of the seizure symptoms. To develop more effective curing methods, it is quite fundamental to understand and reveal the pathogenesis of epilepsy.

As it is known, traditionally epilepsy is referred to a group of diseases characterized by similar symptoms (i.e., epileptic seizures), but not by its pathogenesis. And the underlying pathogenesis of epilepsy may be quite complicated. In the past, due to technological constraints, the pathogenesis of epilepsy is largely unknown. It has been reported that the occurrence of some epilepsy cases turned out to exhibit some degree of familial aggregation, not only implicating the significance of history taking, but also suggesting that the genetic background contributes to the disease [24, 25]. According to some clinical data, it is without any doubt that if a sibling suffers from epilepsy, the brothers and sisters who have similar genetic background inherited from their patients are at higher risk of epilepsy comparing to those who do not [26, 27]. However, the detailed pathogenesis cannot be clearly revealed. Recently, with the development of sequencing technologies, some epilepsy associated genes with either pathological mutations or copy number variants have been identified [28, 29]. Among these genes, a group of functional genes, the sodium channel protein family, encode the core sodium channel in the nerve system [30]; for example, genes like SCN1A and SCN8A encoding core component of the sodium channel have been confirmed to be associated with the progression of hereditary epilepsies [31-33]. Therefore, specific genes may play a definitive role during the initiation and progression of the epilepsy, as hereditary epilepsies are deemed to have a certain genetic background.

Although specific genes have been strongly suggested to be associated with epilepsy, however, it is quite hard to identify the core regulatory genes related to epilepsy by time-consuming experimental methods such as the western blot [34, 35]. Here, based on some known epilepsy-related genes, we presented a new computational workflow to search out potential genes of interest. The random walk with restart (RWR) algorithm was employed in our workflow to search possible novel genes in a protein-protein interaction (PPI) network. Compared to the network methods based on guilt-by-association [36] which only consider the neighbors of known genes [37-39], the RWR algorithm can inspect the whole network to make extensive decisions; that is, methods based on guilt-by-association used part of the network, while the RWR algorithm can utilize the whole network. The brief procedures of our workflow was described as follows. Firstly, the RWR algorithm was executed on a PPI network using validated epilepsy associated genes as seed nodes and genes receiving high probabilities were selected as possible candidate genes. Then, these possible genes were screened by a permutation test, followed by functional association tests, resulting in thirty-three novel epilepsy-related genes. Further analysis indicates that all genes obtained may directly or indirectly contribute to the initiation and progression of epilepsy. To the best of our knowledge, this is the first study attempting to identify core regulatory factors of epilepsy using computational methods. These newly found genes may reveal the underlying mechanisms of epilepsy, and the approach may be extended to solve the similar problems of other complex diseases.

2. Materials and Methods

2.1. Epilepsy Related Genes. 499 genes related to epilepsy were retrieved from EpilepsyGene ( [40], a genetic resource for genes and their mutations related to epilepsy. The epilepsy genes in EpilepsyGene database were collected by searching the PubMed database ( Because our method was executed on a PPI network (referred to in Section 2.2), all 499 genes were linked to their Ensembl IDs. Those without Ensembl IDs or those Ensembl IDs do not occur in the PPI network were excluded. Finally, 470 genes with their Ensembl IDs were obtained for investigation in this study. The detailed information of these genes is provided in Supplementary Material S1 in Supplementary Material available online at

2.2. PPI Network. The interactions between proteins within and outside the cells provide useful information about their activities, properties, and functions. Two proteins that can interact with each other produce a PPI (protein-protein interaction), which often share similar functions or involve in the same biological processes. The PPI network comprised of large amounts of PPIs representing proteins' complicated interaction relationships and remote functional relationships in signaling pathways, such as the proteins involved in regulation and catalysis activity in glycolysis and tricarboxylic acid cycle [41-43]. Some computational predictors and workflows employed PPIs to predict protein functions [44-46] and search for novel genes related to a variety of diseases [47-51]. Therefore, PPIs could be useful to infer novel epilepsy-related genes based on the validated 470 epilepsy genes mentioned in Section 2.1.

To obtain the PPI information and construct a PPI network, the human PPI information was retrieved from STRING (Version 10.0, [52], a well-known public database collecting known and predicted protein-protein interactions. In the current version, it covers 9,643,763 proteins from 2,031 organisms. Interactions reported in STRING are derived from the following five sources: (I) Genomic Context Predictions; (II) High-throughput Lab Experiments; (III) (Conserved) CoExpression; (IV) Automated Text mining; (V) Previous Knowledge in Databases. The human PPIs are collected in the file "9606.protein.links.v10.txt.gz" that can be accessed from the download page of STRING using "Homo sapiens" as a restriction to the data. Accordingly, we obtained 4,274,001 human PPIs covering 19,247 proteins. Because the 4,274,001 human PPIs include not only direct (physical) but also indirect (functional) interactions between proteins, these PPIs can offer relatively more information about the novel genes related to epilepsy.

For each PPI, there are two Ensembl IDs representing two proteins and a score ranging from 150 to 999 that indicates the strength of the interaction. A larger score assigned to a PPI indicates that the two proteins are more likely to interact with each other. For proteins [p.sub.a] and [p.sub.a], their interaction score was denoted as S([p.sub.a], [p.sub.b]). In the network, the 19,247 proteins were denoted as the nodes and two proteins were connected by an edge if and only if they can form a PPI. Thus, there were 4,274,001 edges in the network, and each edge represented a PPI. In addition, the interaction score was added to the network as the weight of the corresponding edge. For convenience, the PPI network was denoted as G in the following sections.

2.3. RWR Algorithm. As a ranking algorithm, the RWR algorithm simulated a walker starting from a seed node or several seed nodes and randomly moved on the network G [53]. In this study, 470 Ensembl IDs of epilepsy genes were set as the seed nodes. Based on them, we aim to mine some potential genes functionally related to epilepsy. In the beginning of the algorithm, an initialization vector P0 was constructed with 19,247 components in it and each component was a score rating the probability of each node being a potential epilepsy-related gene. The probability scores of 470 Ensembl IDs that represented validated epilepsy genes in [p.sub.0] were set to 1/470 (0.0021) and other components were set to zeros. If the vector [p.sub.i] was the probability vector after the RWR algorithm was executed ith round, then the iteration equation can be formulated as follows:

[p.sub.i+1] = (1- r) [A.sup.T] [P.sub.i] + r[P.sub.0], (1)

where A was the column-wise normalized adjacency matrix and r was the probability that it returned to the start nodes, which was set to 0.8 in this study. When probability vector [P.sub.i+1] and [p.sub.i] satisfy the inequality [mathematical expression not reproducible] < 1E - 06, the iteration stopped and [P.sub.i+1] RWR algorithm.

According to the probability vector yielded by RWR algorithm, each node (gene) in the network was given a number representing the probability of it being a novel epilepsy gene. Genes with larger values are more likely to be epilepsy-related genes. Threshold 1E - 05 was adopted in this study; that is, genes receiving scores larger than 1E - 05 were selected from the network G, because it filtered out a large portion of genes and there remained enough genes for further analysis. For convenience, the obtained genes were called RWR genes.

2.4. Filtering Methods. After RWR algorithm was executed on the network, many RWR genes could be selected. However, there are likely many false positive genes among them as elaborated in our previous study [44]. These genes are not special to the epilepsy and should be excluded. In this section, a two-step filtering method was proposed to screen out false positive genes.

2.4.1. Permutation Test. The structure of network G can influence the output of RWR algorithm, which may lead to the false selection of some RWR genes. For example, a node with a degree higher than average degree of the network G may receive a larger score by RWR algorithm even if it was not related to epilepsy. To mine this type of nodes in the network, a permutation test was applied on the network. Firstly, 1,000 Ensembl ID sets, denoted as [S.sub.1], [S.sub.2], ..., [S.sub.1000], were randomly produced and each set contained 470 random Ensembl IDs. Secondly, for each set, the RWR algorithm was executed on the network G using the 470 Ensembl IDs in this set as seed nodes, thereby yielding a probability for each RWR gene. Thirdly, for each RWR gene g, a measurement, namely, permutation FDR, was calculated based on the following equation:

FDR(g) = [THETA]/1000, (2)

where [THETA] was the number of randomly produced sets in which the score of gene g was larger than the score computed by the validated epilepsy related genes. According to (2), the higher permutation FDR an RWR gene had, the less possible the gene was an epilepsy related gene. Because 0.05 was widely used as a common cutoff in statistical test, it was also set to be the threshold of permutation FDR in this study. Therefore, the RWR genes with permutation FDRs less than 0.05 were selected and called candidate genes for further analysis.

2.4.2. Functional Association Test. Among the candidate genes, some of them were functionally highly associated with epilepsy while others weakly associated with it. To select essential genes among them, a functional association test that consisted of two selection schemes was proposed.

It is known that proteins in a PPI with a higher interaction score are more likely to share similar functions. Among the candidate genes, those having strong associations with validated epilepsy-related genes could be the most likely novel epilepsy genes. If a candidate gene has strong associations with exact one validated epilepsy-related gene and has weak or no associations with other epilepsy-related genes, it may still be a novel epilepsy-related gene. In view of this, we believe that using the associations between a candidate gene and its most related epilepsy-related gene is more proper to indicate its associations with epilepsy. Accordingly, an interaction measurement called maximum interaction score (MIS) was calculated for each candidate gene g, which can be defined as

MIS (g)

= max {S (g, g') : g' is an epilepsy-related gene}. (3)

Candidate genes with large MISs mean that it is highly possible that they can directly interact with at least one validated epilepsy gene and may cause the symptoms of epilepsy. In STRING, the value 900 is set to be the cutoff to achieve a highest confidence. Thus, it was also set to be the threshold of MIS; that is, candidate genes with MISs less than 900 were discarded.

Besides, genes related to epilepsy may share some common gene ontology (GO) [54] terms and often occurred in the same Kyoto Encyclopedia of Genes and Genomes (KEGG) [55] pathways. Thus, candidate genes sharing same or similar GO terms and KEGG pathways with validated epilepsy genes are more likely to be the genes related to epilepsy. The enrichment theory of GO terms (KEGG pathways) [56-58] was used to quantitatively measure the relationship between a gene and GO terms (KEGG pathways). For a gene g, let us denote the set containing g and its direct neighbor genes in the PPI network reported in STRING by H(g). Then, the relationship between g and one GO term or KEGG pathway can be encoded into a numeric value as follows:

[mathematical expression not reproducible] (4)

where G represented a GO term or a KEGG pathway, N was the total number of genes in humans, M was the number of genes annotated to G, n was the number of genes in H(g), and m was the number of genes that are in H(g) and annotated by G. The values for all GO terms and KEGG pathways can be collected into a vector ES(g). The similarity of two genes g and g' on GO terms and KEGG pathways can be measured by the proximity of the two vectors ES(g) and ES(g') as follows:

[GAMMA](g, g') = ES(g) x ES(g')/[parallel] ES(g) [parallel] x [parallel] ES(g')[parallel]. (5)

It is clear that if the resultant number of (5) is large, g and g are similar on GO terms and KEGG pathways, implicating a strong relationship between them. With similar arguments on MIS, for each candidate gene g, the maximum function score (MFS) was calculated as follows:

MFS (g)

= max {[GAMMA] (g, g') : g' is an epilepsy-related gene}. (6)

A candidate gene receiving a high MFS means it shares relatively more GO terms and KEGG pathways with at least one validated epilepsy gene. In this study, we tried 0.9 as the threshold of MFS; that is, candidate genes with MFSs larger than 0.9 were selected.

In short, candidate genes resulting from permutation test with MISs larger than or equal to 900 and MFSs larger than 0.9 were selected. For convenience, they were named as core candidate genes.

3. Results

An outline for the procedure of the method, including RWR algorithm and filtering methods described in Sections 2.3 and

2.4, by a flowchart is illustrated in Figure 1. This section would show the detailed results yielded by the method.

As described in Section 2.3, the RWR algorithm was executed on the PPI network G, inwhich the470 Ensembl IDs were used as seed nodes. A probability vector can be obtained, in which each composition represents the probability score of the corresponding node (gene) being a novel epilepsy-related gene. Genes with probabilities larger than 1E - 05 were selected, producing 6,886 RWR genes, which are listed in Supplementary Material S2.

For the 6,886 RWR genes derived from RWR algorithm, a permutation test was applied on them to screen out RWR genes that are not special for epilepsy. The permutation FDR was calculated for each RWR gene, which is provided in Supplementary Material S2. Value 0.05 was set to be the threshold of permutation FDR, thereby producing 980 candidate genes, which are listed in Supplementary Material S3.

To further select genes that are functionally related to epilepsy, a functional association test was applied to the 980 candidate genes. As described in Section 2.4.2, for each candidate gene, we calculated its MIS (cf. (3)) and MFS (cf. (6)). Values 900 and 0.9 were used as the threshold of MIS and MFS, respectively. And finally thirty-three core candidate genes were obtained. These genes were deemed to be closely related to epilepsy and are listed in Table 1. According to some recent published literature as discussed in Section 4, these core candidate genes, which had similar functions with the validated genes, are highly likely to be the novel epilepsy genes.

4. Discussions

For a long time, epilepsy has been regarded as complicated neurological diseases with various pathogenesis. Based on clinical data, the occurrence of epilepsy has shown conspicuous familial aggregation characteristics, implying that genetic background features (such as mutations, copy number variants of genes) may play an irreplaceable role for epilepsy [59-61]. Recent publications have also confirmed such implication. Various epilepsy associated genes have been identified [62-66]. However, it is quite expensive and time-consuming to identify epilepsy associated genes with experiments. Based on our computational method, we identified thirty-three candidate epilepsy associated genes, listed in Table 1. According to some recent publications, all these core candidates show specific relationship with the initiation and progression of epilepsy, validating the effectiveness of our computational method. According to the gene families of these thirty-three core candidate genes, we classified them into seven clusters, shown in Figure 2, and analyzed them accordingly.

Ankyrins. In our prediction list, two of the candidate genes, ANK1 and ANK2, turned out to be the functional members of the ankyrins. As we all know, ankyrins are a group of connexin that link the integral membrane proteins to the cytoskeleton, which have been widely reported that they contribute to cell proliferation, motility, and the maintenance of specialized membrane domains [67-69]. In human bodies, ankyrins have been confirmed to bind to the voltage-gated potassium channel subunits KCNQ2 and KCNQ3, regulating their normal functions [70]. Considering that KCNQ2 and KCNQ3 turned out to directly contributing to the initiation and progression of epilepsy, our predicted genes ANK1 and ANK2 as the functional components of ankyrins may very probably be epilepsy associated genes, validating our prediction [71].

EPH Subfamily. Apart from the Ankyrin protein family, another group of proteins, the EPH subfamily, have also been identified to contribute to epilepsy. In our prediction list, five genes can be classified in such subfamily: EPHA3, EPHA4, EPHA5, EPHA7, and EPHB2. Such five genes all encode the receptors for the erythropoietin-producing hepatoma amplified sequences (EPH), acting as the tyrosine-protein kinase receptor [72,73]. In mouse model, it has been confirmed that the activation of EPH receptor associated genes, like EPHB3, contributes to the onset of epilepsy [74]. Considering the functional similarity and underlying correlations, it is quite reasonable that our predicted genes of EPH receptor family may also contribute to such processes [75]. Apart from that, another publication confirmed that, by stimulating NMDA receptor activity, ERK activates the progression of epilepsy [76]. During the activation, various genes of our predicted EPH family have been identified to promote such biological processes, validating the crucial role of EPH family including our predicted genes EPHA3, EPHA4, EPHA5, EPHA7, and EPHB2 during epilepsy.

Protein Kinases. Two of our predicted candidates, PRKCA and PRKCG, can be clustered into another functional family, the family of serine- and threonine-specific protein kinases. Such two genes turn out to be functional components of the protein kinase C, a core member of the protein family we mentioned above [77-79]. The protein kinase C associated signaling pathway has been widely reported to be associated with epilepsy and may be a candidate therapeutic target for such disease [80]. As two major components for such pathway, our predicted genes PRKCA and PRKCG may definitely contribute to such disease. Apart from such evidence, a specific mutation of PRKCG (SCA-14) has been reported to be associated with a typical movement disorder, which can be called Ramsay Hunt phenotype [81]. Considering that such disorder has been widely identified in epilepsy patients, such mutation may be functionally related to the progression of epilepsy, validating our prediction [82, 83]. As for MAPK7, such gene has been widely regarded as a multifunctional gene that involves various biological processes including proliferation, differentiation, transcription regulation, and development [84]. MAPK7 has been reported to interact with a specific protein Aquaporin 4 (AQP4) in human beings [85]. Since AQP4 has been confirmed to accumulate in neuron cell during epilepsy and contribute to the pathological processes, it is quite reasonable to summarize that, as a functional related protein of AQP4, our predicted gene MAPK7 very probably contributes to epilepsy, validating the accuracy and efficacy of our prediction [85]. The next gene is also a crucial kinase for human beings, the PTK7. Although, different from other proteins from protein tyrosine kinase family, PTK7 lacks detectable catalytic tyrosine kinase activity, it has still been reported to contribute to the functional Wnt signaling pathway and regulate the cellular polarity and adhesion [86]. Though no direct relationship between PTK7 and epilepsy has been confirmed, recent publications reported that PTK7 may participate in the metabolism of antiepileptic drugs (AED), suggesting that there may remain uncovered interactions between PTK7 and epilepsy, validating our prediction [87]. GSK3A, as a multifunctional Ser/Thr protein kinase, has been reported to contribute to glycogen synthesis and transcriptional regulation [88, 89]. Such gene has been reported to be quite essential for the development and maturation of cortical neurons [90]. Considering that cortical neurons, especially the migration of neurons, are quite significant for epilepsy, our predicted gene GSK3A may very probably be epilepsy associated gene [91, 92].

Small GTPases. Apart from PKC associated genes, there are also six genes (rAlA, RAPA, RAP1B, RAP2A, RRAS, and MRAS) that can be clustered into the famous Ras family of small GTPases. Based on recent publications, various small GTPases have been identified to contribute to the progression of epilepsy, including Cdc42, RAB39B [93, 94]. As for our predicted candidates, it has been confirmed that, during the pathological processes of epilepsy, the normal function of RAP1A and its related Ras signaling pathway has been regulated and altered by microRNAs, implying the potential role of Ras signaling pathway during epilepsy [95]. As for RAP1B and RAP2A, RAP1B has been confirmed to be specifically activated in nerve system and contributes to the regeneration of neuronal connectivity, the dysfunction of which turns out to be one of the pathological factors for epilepsy, validating the regulatory role of RAP1B during such disease [96]. RAP2A has been validated and confirmed to contribute to the childhood absence epilepsy, a specific subtype of epilepsy, and may be related with a specific glioma inducing epilepsy associated symptoms, validating our prediction of epilepsy associated genes [97, 98]. As for RRAS and MRAS, considering the functional similarity of MRAS and MAPK, the detailed analysis of such genes can be seen below, while the inner relationship between RRAS and epilepsy has also been revealed in mouse model, validating our prediction [99,100]. RALA, encoding a functional small GTPase belonging to Ras family, has been confirmed to mediate the transmembrane signaling by the occupancy of functional receptors [101]. Such gene has been definitely confirmed to be associated with epilepsy by regulating the drug resistance of such disease [102]. Another gene, RND1, encodes a small GTPase, which does not belong to Ras family but to Rho GTPase family. In response to various extracellular signaling, the protein encoded by such gene turns out to regulate the actin cytoskeleton [103, 104]. In intractable epilepsy, a clinical subtype of epilepsy which is hard to cure, recent publications confirmed the expression of such gene in the central nerve system of the patients, implying that such gene may definitely contribute to the progression and prognosis of such disease [105,106]. MAPK7 and MRAS are two proliferation-associated genes in our candidate epilepsy associated gene list. Among them, MRAS turn out to contribute to Ras signaling pathway, which has been confirmed above to be associated with the progression of epilepsy [95]. What is more, as for MRAS itself, it has been reported that such gene may contribute to the development of brain in early stage and the abnormal activation of such gene may induce epilepsy-like syndrome, validating our prediction [107].

Calmodulin (CALM) Family. Apart from that, such genes may also contribute to the specific seizure like features during the progression of epilepsy, suggesting its core regulatory role [108]. Four genes (CALM1, CALM2, CALM3, and CALM6) of the functional calmodulin (CALM) family have also been predicted to contribute to epilepsy. Genes of calmodulin family mainly act as a calcium binding protein that participate in cell cycle and proliferation associated biological processes [109, 110]. Considering that epilepsy has been confirmed to be associated with abnormal calcium ion transportation, it is quite reasonable to speculate that our predicted genes of CALM family may contribute to epilepsy which has also been confirmed by recent publications [111,112].

14-3-3 Family of Proteins. The remaining group of functional proteins, including YWHAB, YWHAE, YWHAQ, and YWHAZ, that contribute to epilepsy turn out to be encoded by the so-called 14-3-3 family of proteins. Such family of proteins contribute to the signaling transduction by binding to phosphoserine-containing proteins [113]. As we have mentioned above, epilepsy has been confirmed to be associated with the protein kinase C signaling pathway [81]. Recent publications identified that, during the progression of epilepsy, our predicted candidates, proteins of 14-3-3 family, may interact with protein kinase C and further promote the progression, implying the functional role of such genes [114]. Apart from that, such five genes that we sorted out have also been directly identified in epilepsy cases. Take YWHAB as an example. Such gene has been identified to contribute to the regeneration of neurons after physical or chemical injury. The dysfunction of such gene maybe related to the initiation of epilepsy in certain pathological conditions [115]. Therefore, such four genes in our prediction list have all been confirmed to participate in epilepsy associated biological processes, validating the accuracy and efficacy of our prediction.

Other Crucial Genes. Apart from such genes, there still remain six genes with no clear family enrichment in our prediction list that may also contribute to epilepsy in their respective ways. Among such genes, ATP2A2 turns out to encode a significant intracellular pump located in the sarcoplasmic or endoplasmic reticula [116,117]. Epilepsy has been confirmed to be a specific complication of various diseases, including Darier's disease [118, 119]. Our predicted gene ATP2A2 has been identified in clinical cases of Darier's disease and may directly contribute to epilepsy associated symptoms, validating the efficacy and accuracy of our prediction [120]. Another gene INSR turns out to be the functional receptor for a core endogenous hormone insulin, which further activates the downstream of insulin signaling pathway [121, 122]. As for the underlying relationship between INSR and epilepsy, it has been reported that a group of specific mutations in INSR: INSR H1085H C>T, G972R has been confirmed to specifically occur in the epilepsy patients in Han Chinese, validating the specific role of INSR during the progression of epilepsy [123]. Another gene FYN is also a candidate oncogene just like genes from Ras family as we have mentioned above [124, 125]. Recent publications reveal the potential relationship between our predicted gene FYN and amygdala kindling, a specific phenotype associated with epileptogenesis [126]. Participating in mTOR signaling pathway, though the detailed mechanism of the pathology of FYN initiated epilepsyhas not been fully revealed, our predicted gene FYN may definitely be an epilepsy associated gene [126-128].

Phosphodiesterase 6C, encoded by our predicted gene PDE6C, has been confirmed to contribute to pyrimidine metabolism and phototransduction [129, 130]. Autosomaldominant cerebellar ataxia is a specific symptom of epilepsy in human beings [131]. It has been reported that a specific mutation, p.Arg95His, of our predicted gene, PDE6C, may be associated with the autosomal-dominant cerebellar [132]. Considering the inner linkage between autosomal-dominant cerebellar ataxia and epilepsy in human beings, our predicted gene PDE6C may also contribute to the progression of epilepsy, validating our prediction [131]. As we have mentioned above, quite a lot of kinases have been reported to contribute to the abnormal biological process during epilepsy. The remaining two genes SGK1 and YES1 have also been suggested to contribute to the progression of epilepsy. SGK1 turns out to encode a serum/glucocorticoid regulated kinase, which further contributes to the regulation of cellular stress response [133,134]. It has been confirmed that our predicted gene SGK1 is upregulated by a functional cellular component, aldosterone [134]. It has been confirmed that aldosterone has strong chemical and biological effects on epileptic seizures, implying that our predicted gene SGK1 may contribute to the regulation of the typical symptoms of epilepsy, the seizures [135]. Therefore, SGK1 very probably is a functional epilepsy associated gene. The last gene, YES1, also encodes a small GTPase that has been widely regarded as a tumor associated gene [136]. Although no direct report confirms the relationship between YES1 and epilepsy, considering the core regulatory role of small GTPases for epilepsy and that it has been confirmed that YES1 is expressed in brain and central nerve system, it is quite reasonable for us to believe that YES may be a functional epilepsy associated gene [137].

5. Conclusion

Based on our newly developed computational method, we have identified thirty-three novel genes that may contribute to the initiation and progression of epilepsy. According to the comprehensive analyses on these genes, they are strongly suspected to be either directly or indirectly related to epilepsy, validating the effectiveness of our method. In summary, this method can not only contribute to the identification of potential epilepsy associated genes but also provide a new tool to investigate the underlying mechanisms of the pathological processes of epilepsy.

Competing Interests

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

Authors' Contributions

Wei Guo and Dong-Mei Shang contributed equally to this work.


This paper is supported by Norman Bethune Program of Jilin University (2015218), Science and Technology Department of Jilin Province (20160414007GH, 20160414047GH), Education Department of Jilin Province (2015509, 2016449), and Development and Reform Commission of Jilin Province (2015Y032).


[1] S. Wiebe, "Brain surgery for epilepsy," The Lancet, vol. 362, pp. s48-s49, 2003.

[2] S. Baxendale, "Epilepsy at the movies: possession to presidential assassination," Lancet Neurology, vol. 2, no. 12, pp. 764-770, 2003.

[3] F. M. Liu, S. Dai, A. Napoli et al., "Epileptic seizures are induced by intracerebral ablation of astrocytes in the brain, a novel model for dissecting the interaction of neurons with glial cells," Journal of NeuroVirology, vol. 21, pp. S41-S42, 2015.

[4] N. K. Sethi, "Psychogenic non-epileptic seizures--the age matters," Clinical Neurology and Neurosurgery, vol. 120, p. 142, 2014.

[5] R. Zepeda, K. A. Gleason, E. J. Bubrick, D. J. Pallin, and B. A. Dworetzky, "Disparities of epilepsy care in the emergency department," Epilepsia, vol. 50, pp. 307-307, 2009.

[6] J. G. Burneo, N. Jette, W. Theodore et al., "Disparities in epilepsy: report of a systematic review by the North American Commission of the international league against epilepsy," Epilepsia, vol. 50, no. 10, pp. 2285-2295, 2009.

[7] K. Sugai, E. Nakagawa, H. Komaki, H. Sakuma, Y. Saito, and M. Sasaki, "Pharmacotherapy for childhood nonidio-pathic partial epilepsies based on seizure symptoms: retrospective and prospective studies," Epilepsia, vol. 50, pp. 113-113, 2009.

[8] H. Choi, M. R. Winawer, S. Kalachikov, T. A. Pedley, W. A. Hauser, and R. Ottman, "Classification of partial seizure symptoms in genetic studies of the epilepsies," Neurology, vol. 66, no. 11, pp. 1648-1653, 2006.

[9] B. Aktekin, "Up-to-date critical review of the classification of epilepsies and epileptic seizures," Noropsikiyatri Arsivi, vol. 52, no. 2, pp. 109-110, 2015.

[10] S.-H. Lee, J. S. Lim, J.-K. Kim, J. Yang, and Y. Lee, "Classification of normal and epileptic seizure EEG signals using wavelet transform, phase-space reconstruction, and Euclidean distance," Computer Methods and Programs in Biomedicine, vol. 116, no. 1, pp. 10-25, 2014.

[11] M. Jose and S. V. Thomas, "Family history of congenital malformations does not increase the risk of fetal malformations in women with epilepsy," Epilepsia, vol. 56, pp. 31-31, 2015.

[12] C. Alonso-Cerezo, I. Herrera-Peco, V. Fernandez-Millares et al., "Family history of epilepsy resistant to treatment," Revista de Neurologia, vol. 52, no. 9, pp. 522-526, 2011.

[13] U. C. Wieshmann, "Family history of epilepsy in epilepsy and other neurological conditions," Journal of Neurology, Neurosurgery & Psychiatry, vol. 75, pp. 799-800, 2004.

[14] Y. Bahbiti, F. Moutaouakil, A. Ouichou, A. El Hessni, B. Benazzouz, and A. Mesfioui, "Epilepsy: electroencephalogram and brain maturation," Epilepsia, vol. 54, pp. 156-156, 2013.

[15] J. Liao, L. Song, and Y. Chen, "Seizures captured with video-electroencephalogram in infants with epilepsy," Epilepsia, vol. 53, p. 132, 2012.

[16] G. D. Wang, Z. Y. Dai, W. G. Song et al., "Grey matter anomalies in drug-naive childhood absence epilepsy: a voxel-based morphometry study with MRI at 3.0 T," Epilepsy Research, vol. 124, pp. 63-66, 2016.

[17] L. Lipatova and T. Kapustina, "Functional neuroimaging using the method 1H MRS in epilepsy," European Journal of Neurology, vol. 22, supplement 1, p. 632, 2015.

[18] J. Peter, S. Houshmand, T. J. Werner, D. Rubello, and A. Alavi, "Novel assessment of global metabolism by 18F-FDG-PET for localizing affected lobe in temporal lobe epilepsy," Nuclear Medicine Communications, vol. 37, no. 8, pp. 882-887, 2016.

[19] M. Kinney and J. Morrow, "Vitamin K is important for epilepsy in pregnancy," British Medical Journal, vol. 354, article i3929, 2016.

[20] F. M. Snoeijen-Schouwenaars, K. C. Van Deursen, I. Y. Tan, P. Verschuure, and M. H. Majoie, "Vitamin D supplementation in children with epilepsy and intellectual disability," Pediatric Neurology, vol. 52, no. 2, pp. 160-164, 2015.

[21] C. Bodin, S. Aubert, G. Daquin et al., "Responders to vagus nerve stimulation (VNS) in refractory epilepsy have reduced interictal cortical synchronicity on scalp EEG," Epilepsy Research, vol. 113, pp. 98-103, 2015.

[22] S. Mannino, G. Colicchio, R. Di Bonaventura et al., "Patients/caregivers satisfaction following vagal nerve stimulation (Vns) for drug-resistant epilepsies," Epilepsia, vol. 55, no. 1, pp. 103-104, 2014.

[23] A. Cukiert, J. Burattini, and C. Cukiert, "Vagus nerve stimulation (Vns) in refractory epilepsy," Epilepsia, vol. 54, p. 84, 2013.

[24] P. Fabera, H. Krijtova, M. Tomasek et al., "Familial temporal lobe epilepsy due to focal cortical dysplasia type IIIa," Seizure, vol. 31, pp. 120-123, 2015.

[25] A. Chentouf, A. Dahdouh, M. Guipponi et al., "Familial epilepsy in Algeria: clinical features and inheritance profiles," Seizure, vol. 31, pp. 12-18, 2015.

[26] A. Hames and R. Appleton, "Living with a brother or sister with epilepsy: siblings' experiences," Seizure, vol. 18, no. 10, pp. 699-701, 2009.

[27] E. Kurca, M. Grofik, P. Kucera, and P. Varsik, "Familial occurrence of adrenocortical insufficiency in two brothers with all-grove syndrome. A case report of 4A (Allgrove) syndrome with epilepsy and a new AAAs gene mutation," Neuroendocrinology Letters, vol. 26, no. 5, pp. 499-502, 2005.

[28] R. H. Purcell, L. A. Papale, C. D. Makinson et al., "Effects of an epilepsy-causing mutation in the SCN1A sodium channel gene on cocaine-induced seizure susceptibility in mice," Psychopharmacology, vol. 228, no. 2, pp. 263-270, 2013.

[29] A. Escayg and A. L. Goldin, "Sodium channel SCN1A and epilepsy: mutations and mechanisms," Epilepsia, vol. 51, no. 9, pp. 1650-1658, 2010.

[30] T. H. Rhodes, C. G. Vanoye, and A. L. George, "Functional characterization of SCN1A sodium channel mutations associated with familial epilepsy," Biophysical Journal, vol. 88, p. 378a, 2005.

[31] L. Baum, B. S. Haerian, H.-K. Ng et al., "Case-control association study of polymorphisms in the voltage-gated sodium channel genes SCN1A, SCN2A, SCN3A, SCN1B, and SCN2B and epilepsy," Human Genetics, vol. 133, no. 5, pp. 651-659, 2014.

[32] A. J. Barela, S. P. Waddy, J. G. Lickfett et al., "An epilepsy mutation in the sodium channel SCN1A that decreases channel excitability," Journal of Neuroscience, vol. 26, no. 10, pp. 2714-2723, 2006.

[33] C. Lossin, T. H. Rhodes, R. R. Desai et al., "Epilepsy-associated dysfunction in the voltage-gated neuronal sodium channel SCN1A," Journal of Neuroscience, vol. 23, no. 36, pp. 11289-11295, 2003.

[34] A. B. Holt and T. I. Netoff, "Computational modeling of epilepsy for an experimental neurologist," Experimental Neurology, vol. 244, pp. 75-86, 2013.

[35] R. A. Stefanescu, R. G. Shivakeshavan, and S. S. Talathi, "Computational models of epilepsy," Seizure, vol. 21, no. 10, pp. 748-759, 2012.

[36] S. Oliver, "Guilt-by-association goes global," Nature, vol. 403, no. 6770, pp. 601-603, 2000.

[37] 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.

[38] 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 disease," Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 42, pp. 15148-15153, 2004.

[39] 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," American Journal of Human Genetics, vol. 78, no. 6, pp. 1011-1025, 2006.

[40] X. Ran, J. Li, Q. Shao et al., "EpilepsyGene: a genetic resource for genes and mutations related to epilepsy," Nucleic Acids Research, vol. 43, no. 1, pp. D893-D899, 2015.

[41] C. Depre, M. H. Rider, and L. Hue, "Mechanisms of control of heart glycolysis," European Journal of Biochemistry, vol. 258, no. 2, pp. 277-290, 1998.

[42] L. Hue and M. H. Rider, "Role of fructose 2,6-bisphosphate in the control of glycolysis in mammalian tissues," The Biochemical Journal, vol. 245, no. 2, pp. 313-324, 1987.

[43] R. G. Hansford and D. Zorov, "Role of mitochondrial calcium transport in the control of substrate oxidation," Molecular and Cellular Biochemistry, vol. 184, no. 1-2, pp. 359-369, 1998.

[44] L. Chen, Y.-H. Zhang, T. Huang, and Y.-D. Cai, "Identifying novel protein phenotype annotations by hybridizing protein-protein interactions and protein sequence similarities," Molecular genetics and genomics: MGG, vol. 291, no. 2, pp. 913-934, 2016.

[45] L. Hu, T. Huang, X. Shi, W.-C. Lu, Y.-D. Cai, and K.-C. Chou, "Predicting functions of proteins in mouse based on weighted protein-protein interaction network and protein hybrid properties," PLoS ONE, vol. 6, no. 1, Article ID e14556, 2011.

[46] 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.

[47] 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.

[48] J. Zhang, J. Yang, T. Huang, Y. Shu, and L. Chen, "Identification of novel proliferative diabetic retinopathy related genes on protein-protein interaction network," Neurocomputing, vol. 217, pp. 63-72, 2016.

[49] L. Chen, T. Huang, Y.-H. Zhang, Y. Jiang, M. Zheng, and Y.D. Cai, "Identification of novel candidate drivers connecting different dysfunctional levels for lung adenocarcinoma using protein-protein interactions and a shortest path approach," Scientific Reports, vol. 6, Article ID 29849,


[50] L. Chen, J. Yang, T. Huang, X. Kong, L. Lu, and Y.-D. Cai, "Mining for novel tumor suppressor genes using a shortest path approach," Journal of Biomolecular Structure and Dynamics, vol. 34, no. 3, pp. 664-675, 2016.

[51] L. Chen, Z. 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.

[52] D. Szklarczyk, A. Franceschini, S. Wyder et al., "STRING v10: protein-protein interaction networks, integrated over the tree of life," Nucleic Acids Research, vol. 43, no. 1, pp. D447-D452, 2015.

[53] 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. 949958, 2008.

[54] "Gene ontology consortium: going forward," Nucleic Acids Research, vol. 43, no. D1, pp. D1049-D1056, 2015.

[55] M. Kanehisa and S. Goto, "KEGG: kyoto encyclopedia of genes and genomes," Nucleic Acids Research, vol. 28, no. 1, pp. 27-30, 2000.

[56] J. Yang, L. Chen, X. Kong, T. Huang, and Y.-D. Cai, "Analysis of tumor suppressor genes based on gene ontology and the KEGG pathway," PLoS ONE, vol. 9, no. 9, Article ID e107202, 2014.

[57] J. Zhang, Z. Xing, M. Ma et al., "Gene ontology and KEGG enrichment analyses of genes related to age-related macular degeneration," BioMed Research International, vol. 2014, Article ID 450386,10 pages, 2014.

[58] L. Chen, Y.-H. Zhang, M. Zheng, T. Huang, and Y.-D. Cai, "Identification of compound-protein interactions through the analysis of gene ontology, KEGG enrichment for proteins and molecular fragments of compounds," Molecular Genetics and Genomics, vol. 291, no. 6, pp. 2065-2079, 2016.

[59] R. Fjaer, E. Brodtkorb, A.-M. 0ye et al., "Generalized epilepsy in a family with basal ganglia calcifications and mutations in SLC20A2 and CHRNB2," European Journal of Medical Genetics, vol. 58, no. 11, pp. 624-628, 2015.

[60] S. Partemi, S. Cestele, M. Pezzella et al., "Loss-of-function KCNH2 mutation in a family with long QT syndrome, epilepsy, and sudden death," Epilepsia, vol. 54, no. 8, pp. el12-el16, 2013.

[61] B. Berghuis, E. H. Brilstra, D. Lindhout, S. Baulac, G. J. de Haan, and M. van Kempen, "Hyperactive behavior in a family with autosomal dominant lateral temporal lobe epilepsy caused by a mutation in the LGI1/epitempin gene," Epilepsy and Behavior, vol. 28, no. 1, pp. 41-46, 2013.

[62] W. Bi, I. A. Glass, D. M. Muzny et al., "Whole exome sequencing identifies the first STRADA point mutation in a patient with polyhydramnios, megalencephaly, and symptomatic epilepsy syndrome (PMSE)," American Journal of Medical Genetics A, vol. 170, no. 8, pp. 2181-2185, 2016.

[63] M. Gal, D. Magen, Y. Zahran et al., "A novel homozygous splice site mutation in NALCN identified in siblings with cachexia, strabismus, severe intellectual disability, epilepsy and abnormal respiratory rhythm," European Journal of Medical Genetics, vol. 59, no. 4, pp. 204-209, 2016.

[64] G. Li, R. Shi, J. Wuet al., "Association of the hERG mutation with long-QT syndrome type 2, syncope and epilepsy," Molecular Medicine Reports, vol. 13, no. 3, pp. 2467-2475, 2016.

[65] M. G. Sweeney, S. R. Hammans, L. W. Duchen et al., "Mitochondrial DNA mutation underlying Leigh's syndrome: clinical, pathological, biochemical, and genetic studies of a patient presenting with progressive myoclonic epilepsy," Journal of the Neurological Sciences, vol. 121, no. 1, pp. 57-65, 1994.

[66] S. R. Hammans, M. G. Sweeney, M. Brockington et al., "The mitochondrial-DNA transfer Rna(Lys) a-]G(8344) mutation and the syndrome of myoclonic epilepsy with ragged-red fibers (Merrf)--relationship of clinical phenotype to proportion of mutant mitochondrial-DNA," Brain, vol. 116, pp. 617-632, 1993.

[67] K. J. Chang, T. S. Ho, K. Susuki et al., "Paranodal ankyrins: enigmatic glial anchors?" Journal of Neurochemistry, vol. 125, p. 198, 2013.

[68] A. Armani, E. Giacomello, S. Galli, D. Rossi, and V. Sorrentino, "Muscle-specific ankyrins and the organization of the sarcoplasmic reticulum in striated muscle cells," Biophysical Journal 86(1): 222a, vol. 86, no. 1, p. 222a, 2004.

[69] P J. Mohler, A. O. Gramolini, and V. Bennett, "Ankyrins," Journal of Cell Science, vol. 115, no. 8, pp. 1565-1566, 2002.

[70] Z. Pan, T. Kao, Z. Horvath et al., "A common ankyrin-G-based mechanism retains KCNQ and Na V channels at electrically active domains of the axon," Journal of Neuroscience, vol. 26, no. 10, pp. 2599-2613, 2006.

[71] S. R. Cunha and P J. Mohler, "Ankyrin protein networks in membrane formation and stabilization," Journal of Cellular and Molecular Medicine, vol. 13, no. 11-12, pp. 4364-4376, 2009.

[72] J. Chen, W. Song, and K. Amato, "Eph receptor tyrosine kinases in cancer stem cells," Cytokine & Growth Factor Reviews, vol. 26, no. 1, pp. 1-6, 2015.

[73] O. Eriksson, M. Ramstrom, K. Hornaeus, J. Bergquist, D. Mokhtari, and A. Siegbahn, "The Eph tyrosine kinase receptors EphB2 and EphA2 are novel proteolytic substrates of tissue factor/coagulation factor Vila," Journal of Biological Chemistry, vol. 289, no. 47, pp. 32379-32391, 2014.

[74] H. Huang, R. H. Li, J. X. Yuan et al., "Up-regulated ephrinB3/ EphB3 expression in intractable temporal lobe epilepsy patients and pilocarpine induced experimental epilepsy rat model," Brain Research, vol. 1639, pp. 1-12, 2016.

[75] B. Hock, B. Bohme, T. Karn et al., "PDZ-domain-mediated interaction of the Eph-related receptor tyrosine kinase EphB3 and the ras-binding protein AF6 depends on the kinase activity of the receptor," Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 17, pp. 9779-9784, 1998.

[76] A. S. Nateri, G. Raivich, C. Gebhardt et al., "ERK activation causes epilepsy by stimulating NMDA receptor activity," EMBO Journal, vol. 26, no. 23, pp. 4891-4901, 2007

[77] G. Sanchez-Fernandez, S. Cabezudo, A. Caballero et al., "Protein kinase C [zeta] interacts with a novel binding region of Gaq to act as a functional effector," Journal of Biological Chemistry, vol. 291, no. 18, pp. 9513-9525, 2016.

[78] A. M. F. Liu and Y. H. Wong, "G16-mediated activation of nuclear factor kB by the adenosine A1 receptor involves cSrc, protein kinase C, and ERK signaling," Journal of Biological Chemistry, vol. 279, no. 51, pp. 53196-53204, 2004.

[79] C. J. Doering, A. E. Kisilevsky, Z.-P Feng et al., "A single Gj3 subunit locus controls cross-talk between protein kinase C and G protein regulation of N-type calcium channels," The Journal of Biological Chemistry, vol. 279, no. 28, pp. 29709-29717, 2004.

[80] Z. Gajda, R. Torok, Z. Horvath et al., "Protein kinase inhibitor as a potential candidate for epilepsy treatment," Epilepsia, vol. 52, no. 3, pp. 579-588, 2011.

[81] J. E. Visser, B. R. Bloem, and B. P. C. Van De Warrenburg, "PRKCG mutation (SCA-14) causing a Ramsay Hunt phenotype," Movement Disorders, vol. 22, no. 7, pp. 1024-1026, 2007

[82] M.-C. Hsiao, C.-Y. Liu, Y.-Y. Yang, C.-S. Lu, and E.-K. Yeh, "Progressive myoclonic epilepsies syndrome (Ramsay Hunt syndrome) with mental disorder: report of two cases," Psychiatry and Clinical Neurosciences, vol. 53, no. 5, pp. 575-578, 1999.

[83] T. D. Bird and C. M. Shaw, "Progressive myoclonus and epilepsy with dentatorubral degeneration: a clinicopathological study of the Ramsay Hunt syndrome," Journal of Neurology Neurosurgery and Psychiatry, vol. 41, no. 2, pp. 140-149, 1978.

[84] S. Javaid, J. Zhang, G. A. Smolen et al., "MAPK7 regulates emt features and modulates the generation of CTCs," Molecular Cancer Research, vol. 13, no. 5, pp. 934-943, 2015.

[85] T. S. Lee, T. Eid, S. Mane et al., "Aquaporin-4 is increased in the sclerotic hippocampus in human temporal lobe epilepsy," Acta Neuropathologica, vol. 108, no. 6, pp. 493-502, 2004.

[86] V. S. Golubkov, N. L. Prigozhina, Y. Zhang et al., "Protein-Tyrosine Pseudokinase 7 (PTK7) directs cancer cell motility and metastasis," The Journal of Biological Chemistry, vol. 289, no. 35, pp. 24238-24239, 2014.

[87] M. E. Ross, "Gene-environment interactions, folate metabolism and the embryonic nervous system," Wiley Interdisciplinary Reviews: Systems Biology and Medicine, vol. 2, no. 4, pp. 471-480, 2010.

[88] M. Bouskila, M. F. Hirshman, J. Jensen, L. J. Goodyear, and K. Sakamoto, "Insulin promotes glycogen synthesis in the absence of GSK3 phosphorylation in skeletal muscle," American Journal of Physiology--Endocrinology and Metabolism, vol. 294, no. 1, pp. E28-E35, 2008.

[89] V. Matys, O. V. Kel-Margoulis, E. Fricke et al., "TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes," Nucleic Acids Research, vol. 34, pp. D108-D110, 2006.

[90] M. Morgan-Smith, Y. Wu, X. Zhu, J. Pringle, and W. D. Snider, "GSK-3 signaling in developing cortical neurons is essential for radial migration and dendritic orientation," eLife, vol. 3, p. e02663, 2014.

[91] D. Sattarova, M. I. Sigatullina, S. S. Shamansurov, and N. A. Mirsaidova, "Outcome of epilepsy surgery in focal cortical dysplasia," European Journal of Neurology, vol. 19, p. 605, 2012.

[92] F. Tanaka, H. Otsubo, W. C. Gaetz et al., "FDG PET and MEG evaluation of focal cortical dysplasia: comparison with the results of intracranial invasive EEG and epilepsy surgery," Journal of Nuclear Medicine, vol. 41, p. 220, 2000.

[93] M. Giannandrea, V. Bianchi, M. L. Mignogna et al., "Mutations in the small GTPase gene RAB39B are responsible for X-linked mental retardation associated with autism, epilepsy, and macrocephaly," American Journal of Human Genetics, vol. 86, no. 2, pp. 185-195, 2010.

[94] Y. Zhang, J. Liu, G. Luan, and X. Wang, "Inhibition of the small GTPase Cdc42 in regulation of epileptic-seizure in rats," Neuroscience, vol. 289, pp. 381-391, 2015.

[95] A. Kretschmann, B. Danis, L. Andonovic et al., "Different MicroRNA profiles in chronic epilepsy versus acute seizure mouse models," Journal of Molecular Neuroscience, vol. 55, no. 2, pp. 466-479, 2015.

[96] T. Nakamura, S. Yasuda, H. Nagai et al., "Longest neurite-specific activation of Rap1B in hippocampal neurons contributes to polarity formation through RalA and Nore1A in addition to PI3-kinase," Genes to Cells, vol. 18, no. 11, pp. 1020-1031, 2013.

[97] H. Kano, T. Takayama, Y. Midorikawa, and H. Nagase, "Promoter hypomethylation of RAR-related orphan receptor a 1 is correlated with unfavorable clinicopathological features in patients with colorectal cancer," BioScience Trends, vol. 10, no. 3, pp. 202-209, 2016.

[98] L. Wang, W. Zhan, S. Xie et al., "Over-expression of Rap2a inhibits glioma migration and invasion by down-regulating pAKT," Cell Biology International, vol. 38, no. 3, pp. 326-334, 2014.

[99] Q. Zhu, L. Wang, Z. Xiao et al., "Decreased expression of RasGRF1 in the brain tissue of the intractable epilepsy patients and experimental rats," Brain Research, vol. 1493, pp. 99-109, 2013.

[100] M. Adachi, Y. Abe, Y. Aoki, and Y. Matsubara, "Epilepsy in RAS/MAPK syndrome: two cases of cardio-facio-cutaneous syndrome with epileptic encephalopathy and a literature review," Seizure, vol. 21, no. 1, pp. 55-60, 2012.

[101] Z. Lu, A. Hornia, T. Joseph et al., "Phospholipase D and RalA cooperate with the epidermal growth factor receptor to transform 3Y1 rat fibroblasts," Molecular and Cellular Biology, vol. 20, no. 2, pp. 462-467, 2000.

[102] E. Manguoglu, S. Akdeniz, N. Dundar et al., "RLIP76 gene variants are not associated with drug response in turkish epilepsy patients," Balkan Journal of Medical Genetics, vol. 14, no. 1, pp. 25-30, 2011.

[103] C. D. Nobes, I. Lauritzen, M.-G. Mattei, S. Paris, A. Hall, and P Chardin, "A new member of the Rho family, Rnd1, promotes disassembly of actin filament structures and loss of cell adhesion," Journal of Cell Biology, vol. 141, no. 1, pp. 187-197, 1998.

[104] Y. Ishikawa, H. Katoh, and M. Negishi, "A role of Rnd1 GTPase in dendritic spine formation in hippocampal neurons," Journal of Neuroscience, vol. 23, no. 35, pp. 11065-11072, 2003.

[105] C. J. Wingard, V. Chintalgattu, G. Harris, R. Nolan, J. Narron, and L. C. Katwa, "Testosterone-dependent expression of RhoA, ROCK I, ROCK II and Rnd1 in rat corpus cavernosum," The FASEB Journal, vol. 19, p. A123, 2005.

[106] S. M. Zanata, I. Hovatta, B. Rohm, and A. W. Piischel, "Antagonistic effects of Rnd1 and RhoD GTPases regulate receptor activity in semaphorin 3A-induced cytoskeletal collapse," Journal of Neuroscience, vol. 22, no. 2, pp. 471-477, 2002.

[107] R. A. Teixeira, V. A. Zanardi, L. M. Li, S. L. M. Santos, and F. Cendes, "Epilepsy and destructive brain insults in early life: a topographical classification on the basis of MRI findings," Seizure, vol. 13, no. 6, pp. 383-391, 2004.

[108] Y.-S. Bae, W. Chung, K. Han et al., "Down-regulation of RalBP1 expression reduces seizure threshold and synaptic inhibition in mice," Biochemical and Biophysical Research Communications, vol. 433, no. 2, pp. 175-180, 2013.

[109] M. C. Lee, S. S. Ban, Y.-J. Woo, and S. U. Kim, "Calcium/calmodulin kinase II activity of hippocampus in kainite-induced epilepsy," Journal of Korean Medical Science, vol. 16, no. 5, pp. 643-648, 2001.

[110] L. S. Butler, A. J. Silva, A. Abeliovich, Y. Watanabe, S. Tonegawa, and J. O. Mcnamara, "Limbic epilepsy in transgenic mice carrying a Ca2+/calmodulin-dependent kinase II a-subunit mutation," Proceedings of the National Academy of Sciences of the United States of America, vol. 92, no. 15, pp. 6852-6855, 1995.

[111] P. Ambrosino, A. Alaimo, S. Bartollino et al., "Epilepsy-causing mutations in Kv72 C-terminus affect binding and functional modulation by calmodulin," Biochimica et Biophysica Acta-Molecular Basis of Disease, vol. 1852, no. 9, pp. 1856-1866, 2015.

[112] S. B. Churn, L. D. Kochan, and R. J. Delorenzo, "Chronic inhibition of Ca2+/calmodulin kinase II activity in the pilocarpine model of epilepsy," Brain Research, vol. 875, no. 1-2, pp. 66-77, 2000.

[113] A. J. Muslin, J. W. Tanner, P. M. Allen, and A. S. Shaw, "Interaction of 14-3-3 with signaling proteins is mediated by the recognition of phosphoserine," Cell, vol. 84, no. 6, pp. 889-897, 1996.

[114] Y. S. Kim, M. Y. Choi, Y. H. Kim et al., "Protein kinase Cdelta is associated with 14-3-3 phosphorylation in seizure-induced neuronal death," Epilepsy Research, vol. 92, no. 1, pp. 30-40, 2010.

[115] S. Shinoda, S. L. Skradski, T. Arakiet al., "Formation of a tumour necrosis factor receptor 1 molecular scaffolding complex and activation of apoptosis signal-regulating kinase 1 during seizure-induced neuronal death," European Journal of Neuroscience, vol. 17, no. 10, pp. 2065-2076, 2003.

[116] E. A. Knopp, C. Saraceni, J. Moss, J. M. McNiff, and K. A. Choate, "Somatic ATP2A2 mutation in a case of papular acantholytic dyskeratosis: mosaic Darier disease," Journal of Cutaneous Pathology, vol. 42, pp. 853-857, 2015.

[117] J. Dhitavat, L. Dode, N. Leslie, S. Burge, and A. Hovnanian, "Effects of mutations in ATP2A2 on calcium transport across sarco/endoplasmic reticulum (ER) membrane," Journal of Investigative Dermatology, vol. 117, pp. 774-774, 2001.

[118] A. Takagi, M. Kamijo, and S. Ikeda, "Darier disease," Journal of Dermatology, vol. 43, no. 3, pp. 275-279, 2016.

[119] R. P. Dodiuk-Gad, E. Cohen-Barak, M. Khayat et al., "Darier disease in Israel: combined evaluation of genetic and neuropsychiatric aspects," British Journal of Dermatology, vol. 174, no. 3, pp. 562-568, 2016.

[120] N. J. O. Jacobsen, I. Lyons, B. Hoogendoorn et al., "ATP2A2 mutations in Darier's disease and their relationship to neuropsychiatric phenotypes," Human Molecular Genetics, vol. 8, no. 9, pp. 1631-1636, 1999.

[121] K. Karimi, T. Mahmoudi, N. Karimi et al., "Is there an association between variants in candidate insulin pathway genes IGFI, IGFBP-3, INSR, and IRS2 and risk of colorectal cancer in the Iranian Population?" Asian Pacific Journal of Cancer Prevention, vol. 14, no. 9, pp. 5011-5016, 2013.

[122] S. Pechlivanis, B. Pardini, J. L. Bermejo et al., "Insulin pathway related genes and risk of colorectal cancer: INSR promoter polymorphism shows a protective effect," Endocrine-Related Cancer, vol. 14, no. 3, pp. 733-740, 2007.

[123] F. Che, Q. Fu, X. Li et al., "Association of insulin receptor H1085H C>T, insulin receptor substrate 1 G972R and insulin receptor substrate 2 1057G/A polymorphisms with refractory temporal lobe epilepsy in Han Chinese," Seizure, vol. 25, pp. 178-180, 2015.

[124] Q. Wang, J. Qian, F. Wang, and Z. Ma, "Cellular prion protein accelerates colorectal cancer metastasis via the Fyn-SP1-SATB1 axis," Oncology Reports, vol. 28, no. 6, pp. 2029-2034, 2012.

[125] A. Strom, S. Diecke, G. Hunsmann, and A. W. Stuke, "Cellular prion protein promotes glucose uptake through the Fyn-HIF-2 alpha-Glut1 pathway to support colorectal cancer cell survival," Cancer Science, vol. 103, pp. 606-606, 2011.

[126] D. P. Cain, S. G. N. Grant, D. Saucier, E. L. Hargreaves, and E. R. Kandel, "Fyn tyrosine kinase is required for normal amygdala kindling," Epilepsy Research, vol. 22, no. 2, pp. 107-114, 1995.

[127] X. Yang, C. Marshall, T. Dentchev et al., "A topical PI3K/mTOR inhibitor induces regression of squamous cell carcinomas in K14-Fyn Y528F mice," Journal of Investigative Dermatology, vol. 132, p. S25, 2012.

[128] Y. Bermudez, S. P. Stratton, G. T. Bowden et al., "Abstract 3673: expression profile of phosphorylated proteins from the mTOR and Fyn/RSK2 signaling pathways in solar UV-induced skin carcinogenesis," Cancer Research, vol. 71, no. 8 supplement, pp. 3673-3673, 2011.

[129] B. Chang, T. Grau, S. Dangel et al., "A homologous genetic basis of the murine cpfl1 mutant and human achromatopsia linked to mutations in the PDE6C gene," Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 46, pp. 19581-19586, 2009.

[130] S. E. Martinez, C. C. Heikaus, R. E. Klevit, and J. A. Beavo, "The structure of the GAF a domain from phosphodiesterase 6C reveals determinants of cGMP binding, a conserved binding surface, and a large cGMP-dependent conformational change," Journal of Biological Chemistry, vol. 283, no. 38, pp. 25913-25919, 2008.

[131] A. Filla, G. De Michele, S. Cocozza et al., "Early onset autosomal dominant dementia with ataxia, extrapyramidal features, and epilepsy," Neurology, vol. 58, no. 6, pp. 922-928, 2002.

[132] M. Coutelier, I. Blesneac, A. Monteil et al., "A recurrent mutation in CACNA1G alters Cav3.1 T-type calcium-channel conduction and causes autosomal-dominant cerebellar ataxia," American Journal of Human Genetics, vol. 97, no. 5, pp. 726-737, 2015.

[133] H. Wang, D. Xu, M. F. Toh, A. C. Pao, and G. You, "Serum- and glucocorticoid-inducible kinase SGK2 regulates human organic anion transporters 4 via ubiquitin ligase Nedd4-2," Biochemical Pharmacology, vol. 102, pp. 120-129, 2016.

[134] B. Friedrich, Y. Feng, P. Cohen et al., "The serine/threonine kinases SGK2 and SGK3 are potent stimulators of the epithelial Na+ channel [alpha], [beta], [gamma]-ENaC," Pflugers Archiv European Journal of Physiology, vol. 445, no. 6, pp. 693-696, 2003.

[135] M. J. Seller and R. G. Spector, "Effect of aldosterone and cortisol on leptazol-induced seizures in rats," British Journal of Pharmacology and Chemotherapy, vol. 19, no. 2, pp. 271-273, 1962.

[136] W. Tan, S.-G. Lim, and T. M. C. Tan, "Up-regulation of microRNA-210 inhibits proliferation of hepatocellular carcinoma cells by targeting Yes1," World Journal of Gastroenterology, vol. 21, no. 46, pp. 13030-13041, 2015.

[137] P. R. Patel, H. Sun, S. Q. Li et al., "Identification of potent Yes1 kinase inhibitors using a library screening approach," Bioorganic & Medicinal Chemistry Letters, vol. 23, no. 15, pp. 4398-4403, 2013.

Wei Guo, (1) Dong-Mei Shang, (1) Jing-Hui Cao, (2) Kaiyan Feng, (3) Yi-Chun He, (2) Yang Jiang, (4) ShaoPeng Wang, (5) and Yu-Fei Gao (2)

(1) Department of Outpatient, China-Japan Union Hospital of Jilin University, Changchun 130033, China

(2) Department of Neurosurgery, China-Japan Union Hospital of Jilin University, Changchun 130033, China

(3) Department of Computer Science, Guangdong AIB Polytechnic, Guangzhou 510507, China

(4) Department of Surgery, China-Japan Union Hospital of Jilin University, Changchun 130033, China

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

Correspondence should be addressed to Yu-Fei Gao;

Received 23 October 2016; Accepted 15 January 2017; Published 1 February 2017

Academic Editor: Ansgar Poetsch

Caption: Figure 1: The flowchart of RWR algorithm and filtering methods for identifying core candidate genes.

Caption: Figure 2: The distribution of the thirty-three core candidate genes according to their protein families.
Table 1: Thirty-three core candidate genes identified by our method.

Gene       Ensembl ID         Probability    Permutation
symbol                                       FDR

ANK2       ENSP00000349588    4.53E - 05     <0.001
ANK1       ENSP00000265709    3.52E - 05     0.034
EPHA7      ENSP00000358309    3.28E - 05     0.025
EPHA5      ENSP00000273854    3.46E - 05     0.032
PRIKCG     ENSP00000263431    4.75E - 05     0.025
PTK7       ENSP00000230419    3.10E - 05     0.023
EPHA3      ENSP00000337451    4.40E - 05     0.017
PDE6C      ENSP00000360502    3.61E - 05     0.032
EPHA4      ENSP00000281821    4.16E - 05     0.008
YWHAQ      ENSP00000238081    5.50E - 05     0.004
GSK3A      ENSP00000222330    6.40E - 05     0.017
CALM1      ENSP00000349467    6.55E - 05     0.023
EPHB2      ENSP00000363763    4.39E - 05     0.026
PRIKCA     ENSP00000408695    7.02E - 05     0.039
CALM2      ENSP00000272298    5.52E-05       0.048
YWHAE      ENSP00000264335    1.63E - 04     0.009
YWHAB      ENSP00000300161    5.56E - 05     0.047
ATP2A2     ENSP00000440045    5.03E - 05     0.005
YES1       ENSP00000324740    6.54E - 05     0.002
CALML3     ENSP00000315299    3.96E - 05     0.027
SGK1       ENSP00000356832    4.90E - 05     0.034
CALML6     ENSP00000304643    4.44E - 05     0.01
YWHAZ      ENSP00000309503    1.86E - 04     0.002
MAPK7      ENSP00000311005    7.15E - 05     0.031
RRAS       ENSP00000246792    5.94E - 05     0.004
RAP2A      ENSP00000245304    4.97E - 05     0.021
RALA       ENSP00000005257    6.57E - 05     0.006
RAP1B      ENSP00000250559    5.28E - 05     0.009
MRAS       ENSP00000289104    4.70E - 05     0.022
INSR       ENSP00000303830    7.27E - 05     0.024
PAP1A      ENSP00000348786    6.49E - 05     0.014
RND1       ENSP00000308461    5.70E - 05     0.002
FYN        ENSP00000346671    1.02E - 04     <0.001

Gene       MIS     MFS

ANK2       990     0.997
ANK1       995     0.995
EPHA7      906     0.988
EPHA5      906     0.988
PRIKCG     905     0.987
PTK7       943     0.987
EPHA3      912     0.986
PDE6C      900     0.980
EPHA4      990     0.980
YWHAQ      999     0.978
GSK3A      977     0.976
CALM1      966     0.976
EPHB2      908     0.975
PRIKCA     992     0.975
CALM2      985     0.974
YWHAE      999     0.973
YWHAB      999     0.971
ATP2A2     908     0.970
YES1       967     0.969
CALML3     906     0.968
SGK1       999     0.966
CALML6     909     0.965
YWHAZ      999     0.958
MAPK7      999     0.956
RRAS       951     0.941
RAP2A      940     0.937
RALA       981     0.932
RAP1B      972     0.931
MRAS       932     0.931
INSR       996     0.927
PAP1A      995     0.925
RND1       996     0.903
FYN        999     0.900
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:Guo, Wei; Shang, Dong-Mei; Cao, Jing-Hui; Feng, Kaiyan; He, Yi-Chun; Jiang, Yang; Wang, ShaoPeng; Ga
Publication:BioMed Research International
Article Type:Report
Geographic Code:9CHIN
Date:Jan 1, 2017
Previous Article:Rheumatoid Arthritis Affecting the Upper Cervical Spine: Biomechanical Assessment of the Stabilizing Ligaments.
Next Article:Autologous Stem Cells Combined Core Decompression for Treatment of Avascular Necrosis of the Femoral Head: A Systematic Meta-Analysis.

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