A computational analysis of protein-protein interaction networks in neurodegenerative diseases
© Goñi et al; licensee BioMed Central Ltd. 2008
Received: 09 December 2007
Accepted: 20 June 2008
Published: 20 June 2008
Recent developments have meant that network theory is making an important contribution to the topological study of biological networks, such as protein-protein interaction (PPI) networks. The identification of differentially expressed genes in DNA array experiments is a source of information regarding the molecular pathways involved in disease. Thus, considering PPI analysis and gene expression studies together may provide a better understanding of multifactorial neurodegenerative diseases such as Multiple Sclerosis (MS) and Alzheimer disease (AD). The aim of this study was to assess whether the parameters of degree and betweenness, two fundamental measures in network theory, are properties that differentiate between implicated (seed-proteins) and non-implicated nodes (neighbors) in MS and AD. We used experimentally validated PPI information to obtain the neighbors for each seed group and we studied these parameters in four networks: MS-blood network; MS-brain network; AD-blood network; and AD-brain network.
Specific features of seed-proteins were revealed, whereby they displayed a lower average degree in both diseases and tissues, and a higher betweenness in AD-brain and MS-blood networks. Additionally, the heterogeneity of the processes involved indicate that these findings are not pathway specific but rather that they are spread over different pathways.
Our findings show differential centrality properties of proteins whose gene expression is impaired in neurodegenerative diseases.
The structural and functional relationships underlying the organization of living systems imply the need to coordinate molecular interactions, principally those involving gene expression and protein activity. Although the genome in each cell of a given organism is virtually the same, dynamic changes in gene expression and thus in the protein content depend on the functional state of the cell . Genome-wide expression profiles using DNA arrays, together with the development of bioinformatic approaches , enable both genetic and protein-protein interaction (PPI) networks to be modeled, thereby helping to understand how biological networks operate .
From a systems point of view, the arrangement of molecular networks from gene expression data based on known interactions permits the understanding of the basic mechanisms upon which the complexity and adaptability of a living cell is founded . This information also helps to decipher processes involved in illness, for instance the molecular heterogeneity of cancer . However, and consistent with the model of multifactorial diseases, it is difficult to find genes that account for direct genotype-phenotype correlations . Thus, network modeling and topological analysis at the meso-scale (intermediate level between local and global features of networks) may provide additional knowledge about the particular properties of genes and proteins involved in diseases of multifactorial nature, where the pathogenesis does not depend on the malfunction of a single gene or protein . In this case, the analysis of gene, protein and pathway interactions might indicate common properties of good candidates to be targeted by therapy. In addition, understanding the emergent properties of a system might help identify new targets that would not be captured by a molecular approach .
Multiple Sclerosis (MS) is a chronic inflammatory and neurodegenerative disease of the central nervous system (CNS) . Although its etiology remains elusive, the interplay between environmental and genetic factors is ultimately thought to be critical for the development of the disease. MS is considered as an autoimmune disease due to the presence of inflammatory infiltrates in the brain in the absence of infection, and through its association with HLA alleles, among other factors . The chronic inflammatory activity within the CNS is the main mediator of tissue damage, even in the late neurodegenerative stage of the disease that involves widespread demyelination and axon loss . In addition to the autoimmune processes, MS also has a neurodegenerative component whereby axons and neurons are lost through unknown processes in the late chronic stages of the disease. Several lines of evidence suggest that the degeneration of demyelinated axons is the most important factor in MS neurodegeneration . Thus, MS is a multifactorial disease in which many immune system and CNS pathways are involved . Current therapies partially ameliorate the inflammatory process, but more effective therapeutic approaches are required to stop disease progression and prevent neurodegeneration.
Alzheimer's Disease (AD) is the most common neurodegenerative disease and it represents one of the biggest unmet needs in modern medicine. AD is characterized by the loss of neurons in conjunction with the presence of oxidative stress, axonal dystrophy, mature senile plaques and neurofibrillary tangles . A set of gene mutations involved in the amyloid beta and tau pathways have been associated with hereditary AD and, in conjunction with neuropathological findings, it has been demonstrated that amyloid and tau are involved in the pathogenesis of AD. However, current evidence suggests that sporadic AD is a multifactorial disease in which many pathways are involved [15, 16]. Indeed, recent studies have also identified molecular abnormalities in the blood of patients with AD . Because the AD therapies available are symptomatic, and considering the epidemic proportions of this disease in western countries, the development of new therapies to stop its progress is an important health priority.
To better understand the basis of neurodegenerative diseases, we set out to study the centrality related features of proteins whose genes were differentially expressed (seed proteins) in MS and AD with respect to their protein neighbors. The main features examined were the degree and the betweenness of these seed proteins and its comparison to their neighbors.
Network measurements for the four disease networks.
number of nodes
mean shortest path length
Connectivity analysis of the disease networks.
As indicated in Tables 1 and 2, we found seed-proteins displayed a lower average degree with respect to the degree of their PPI neighbors in both diseases and in both tissues. In addition, direct interactions between seed proteins were very low: MS-blood: 1 (total links: 386); MS-brain: 4 (total links: 231); AD-blood: 2 (total links: 245); AD-brain: 2 (total links: 191). There were not big differences in the network measurements (<k>, <C>, D and mspl) among the four disease networks studied, indicating certain homogeneity in the architecture of the PPI subnetworks analyzed in this study (Table 1). With regards the centrality of seed proteins, our study shows a low correspondence between their degree and betweenness (Table 2), indicating that critical proteins in disease pathogenesis are not highly connected, but tend to be located in bottleneck regions.
Network theory provides a useful tool to study the complexity of neurodegenerative diseases. Here we report a novel approach to study PPI networks at the meso-scale based on the products of genes differentially expressed in MS and AD. Our approach was to analyze PPI networks based on seed-protein neighborhoods from the genes that were differentially expressed in DNA array studies. The method for growing networks from seed-proteins is critical for determining their topological properties . For this reason, the network growth in our study was carried out by expanding it through experimentally validated protein interactions. The stability, dynamics and functioning of networks are generally characterized by determining the topology of the map, i.e., the configuration of its nodes and the connecting edges . For example, networks with a scale-free topology are resistant to random failure but they are vulnerable to targeted attack, specifically against the most connected nodes. In terms of identifying common properties among the genes involved in neurodegenerative disorders, very interesting results were obtained by carrying out a topological analysis. There were multiple pathways affected by proteins with a low degree, and half the time with high betweenness.
During the last decade, network studies have been applied to biological data bearing in mind that the degree of connectivity is a key property of any network, as demonstrated in yeast . The most common approach to identify key nodes consists of obtaining networks from high throughput data and having obtained the network, searching for the most connected nodes (hubs). The underlying assumption was that these hubs could be critical to explain the pathogenesis of diseases. However, betweenness is another key indicator of centrality that demonstrates how nodes with a low degree of centrality may be relevant in a network (i.e. bottleneck effect) [22, 23]. Our study was performed from a novel viewpoint, since we analyzed whether degree was any different respect to the PPI neighbors starting from critical nodes (in terms of differentially expressed genes). Accordingly, we found that the degree of seed-proteins was lower than that of the PPI neighbors, situating seed proteins in peripheral regions of the network. According to our results of the GO analysis, such peripheral regions are distributed among several pathways that could be involved in disease. Indeed, our results are in agreement with a recent study in asthma showing that hubs exhibit small changes in gene expression . Therefore our results support the application of strategies other than those previously applied, whereby only hubs that might compromise the robustness of networks were generally searched [25, 26].
The fact that we obtained similar results with regards the low average degree of seed proteins in two diseases and two different tissues suggests that this might be a common property in complex diseases, more relevant than the issues associated with the techniques applied such as DNA array technology. However, our approach relies on the current knowledge of interactions, which depends more on how much the gene/protein has been studied than on how many real interactions it participates in. Although it is more difficult to relate gene expression data from hubs with that of other genes, this would not bias our analysis since we focused on whether genes that are differentially expressed (but not necessarily correlated) have a particular distribution with regards their neighbors (neighbors found in a database that includes structural and experimental evidence and not correlation profiles).
We can consider complex diseases as an evolutionary stage in which the pathogenesis process hijacks the robustness of the biological pathways. Such an event may be followed by a cascade of failures in these pathways [8, 27]. In this sense and from a therapeutic point of view, it may be necessary to target many of the pathways involved following a systems biology rationale, and based on the dynamics and topology of the networks involved. The aim of this therapy would be to drive those pathways to a non-pathological state or at least, to a less deleterious state.
The topological implications of the observed scale-free properties in biological networks would indicate that the best therapeutic targets to modify network behavior would be the genes (or proteins) corresponding to the hubs in the network. However, our findings suggest that less extensively connected proteins might be more appropriate therapeutic targets than hyper-connected ones, at least in neurodegenerative diseases. The fact that in both diseases (MS and AD) and in two different tissues analyzed (blood and cerebral tissue), seed-proteins are weakly connected nodes taking part in many different pathways, strengthens the concept of the multifactorial pathogenesis of neurodegenerative diseases. Thus, our results suggest that to modify the disease course we need to target many genes or proteins in several pathways. In a previous network analysis in MS we demonstrated that therapies act on different regions of the gene network that control T-cell activation, suggesting that a pleiotropic activity is required in order to modulate the immune response . In addition, recent network studies in neurodegenerative diseases suggest that several common pathways are involved in their pathogenesis, reinforcing the need to interact with several regions of the PPI network [29, 30]. Another reason why hubs might not be good therapeutic targets is because their critical role in the network modules might prevent them from fluctuating substantially. For the same reason, we can speculate that networks would poorly tolerate modifications in hub behavior without spreading such changes across the network and thereby, inducing significant side effects.
The results we present here indicate that both neurodegenerative diseases (MS and AD) share common characteristics, such as the low degree of seed-proteins and in two of the four disease networks, a high degree of betweenness. These findings mainly situate seed-proteins in peripheral regions of the PPI map (in terms of degree), involved in different pathways as indicated by the associated GO terms and the direct interactions, and integrated into subnetworks of the complete Human proteome network.
Some definitions are introduced to better explain the development of our topological studies:
– Seed-proteins: proteins whose genes were differentially expressed in DNA array studies focused on the specific disease and on a particular tissue. In this study, the diseases considered are Multiple Sclerosis (MS) and Alzheimer Disease (AD) and the tissues are blood and brain.
– MS-blood seed-proteins: proteins whose genes were differentially expressed in DNA array studies of blood from MS patients .
– MS-blood neighbors: nodes selected as a consequence of adding experimentally validated interactions starting from seed-proteins.
– MS-blood network: network that includes MS seed-proteins, MS-neighbors and their interactions. Only seed-proteins linked to neighboring proteins were included in the network analysis (isolated seed-proteins were not included in the analysis shown in Table 1).
– MS-brain seed-proteins: proteins whose genes were differentially expressed in DNA array studies of brain tissue from MS patients .
– MS-brain neighbors: nodes selected as a consequence of adding experimentally validated interactions starting from seed-proteins.
– MS-brain network: network that includes MS seed-proteins, MS-neighbors and their interactions.
– AD-blood seed-proteins: proteins whose genes were differentially expressed in DNA array studies of blood from AD patients .
– AD-blood neighbors: nodes selected as a consequence of adding experimentally validated interactions starting from AD seed-proteins.
– AD-blood network: network that includes AD seed-proteins, AD-neighbors and their interactions.
– AD-brain seed-proteins: proteins whose genes were differentially expressed in DNA array studies of brain tissue from AD patients .
– AD-brain neighbors: nodes selected as a consequence of adding experimentally validated interactions starting from AD seed-proteins.
– AD-brain network: network that includes AD seed-proteins, AD-neighbors and their interactions.
– Disease-networks: the term used to refer to the networks obtained from MS or AD patients that contain seed-proteins and their neighbors. It is important to note that we did not consider neighbors as newly proposed proteins implicated in the disease but rather, they were taken simply to capture the network context in which seed-proteins are located.
– Giant component: term used to refer the largest part of a network whose nodes are connected either directly or indirectly.
Gene expression data
Starting from seed-proteins involved in either MS or AD, we obtained a PPI network through the interaction of these proteins with their direct neighbors. A general scheme of the approach adopted here is presented in Figure 1. The growth of each network was carried out using the STRING database  and the parameters used to generate the network in the STRING database were: active prediction method – experiments; confidence score – 0.7-high confidence; network depth – 2 (only direct neighbors); and an edge scaling factor of 80%. This configuration implies that only the experimental evidence of interactions with a high level of confidence were extracted from the database as valid links for each PPI network. A detailed description of each parameter can be found elsewhere . We did not consider either the direction of each protein interaction or self-interactions. Network files in Pajek format (.net) of each network are available as additional files [see additional files 5, 6, 7, 8].
Topological analysis and measurements of centrality
In order to characterize the disease networks (all of them undirected graphs) and assess the centrality of seed-proteins we applied the following measurements [25, 35] using MATLAB (The Mathworks, MA, United States):
-Degree (k): in an undirected graph, the degree of a vertex is the number of adjacent links. In this study, it represents the number of experimentally validated interactions (links) that connect one protein (node) to its neighbors.
– Average degree (<k>): it represents the mean of all degree values of nodes in a network.
– Clustering coefficient (<C>): is the average clustering coefficient of nodes where the clustering coefficient of a node i (C i ) is the proportion of links between the nodes within the i- neighborhood divided by the number of links that could possibly exist between them.
– Mean shortest path length (mspl): is the average of the steps (number of links) needed to connect every pair of nodes through their shortest path.
– Diameter (D): is the longest among all shortest paths, i.e. the minimum number of links that separate the two most distant nodes in a network.
– Betweenness centrality: measures how often nodes occur on the shortest paths between other nodes. When combined with the degree, it is a key measure to assess the relevance of the location of nodes within a network (vertices within a graph).
Gene symbol identities corresponding to the four different lists of seed-proteins were loaded into the ExPlainTm 2.3 Tool , where functional groups of Gene Ontology Biological Processes (GO-BP) were detected using a p-value threshold of 0.05 as the classification criteria and one as the minimal number of genes assigned to a group (i.e.: number of hits). A FDR multiple hypothesis test adjustment was further carried out using the Benjamini-Hochberg (BH) procedure  and taking the total number of GO-BP as those in which at least one protein of the seed-protein list is included.
We used the Kolmogorov-Smirnov test to compare the distributions of degree and betweenness between seed-proteins and neighbors for each disease. The level of significance was set at p < 0.05.
Central nervous system
mean shortest path length
We thank Cristian von Mering and Ivan Martinez-Forero for helping us with the use of STRING database, Mark Sefton for linguistic assistance, and Ricard V. Solé and Bartolomé Bejarano for their comments and suggestions. J.G, N.VdM. and J.S. are fellows supported by the Government of Navarra (Spain), the Government of the Basque Country (Spain) and the Spanish Ministry of Health (CM05/00222), respectively. This work was partially supported by Junta de Andalucía (BIO-302; to F.J.E.), Spanish Ministry of Science and Technology (BFM2003-05878 to S.A.T.) and Spanish Ministry of Health (FIS PI060117), the CIMA UTE project and the European Commission (6FP: ComplexDis) to PV.
- Lodish H, Berk A, Zipursky SL, Matsudaira P, Baltimore D, Darnell J: Molecular Cell Biology . 2000, New York , W.H. Freeman and Company, 4thGoogle Scholar
- Villoslada P, Oksenberg J: Neuroinformatics in clinical practice: are computers going to help neurological patients and their physicians? . Future Neurology. 2006, 1 (2): 1-12. 10.2217/14796708.1.2.159.View ArticleGoogle Scholar
- Xia Y, Yu H, Jansen R, Seringhaus M, Baxter S, Greenbaum D, Zhao H, Gerstein M: Analyzing cellular biochemistry in terms of molecular networks. Ann Rev Biochem. 2004, 73 (1051-1087):View ArticlePubMedGoogle Scholar
- Hiesinger PR, Hassan BA: Genetics in the age of systems biology. Cell. 2005, 123 (1173-1174):View ArticlePubMedGoogle Scholar
- Rhodes DR, Chinnaiyan AM: Integrative analysis of the cancer transcriptome. Nat Genet. 2005, 37 (Suppl): S31-S37. 10.1038/ng1570View ArticlePubMedGoogle Scholar
- Gunsalus KC, Ge H, Schetter AJ, Goldberg DS, Han JD, Hao T, Berriz GH, Bertin N, Huang J, Chuang LS, Li N, Mani R, Hyman AA, Sonnichsen B, Echeverri CJ, Roth FP, Vidal M, F P: Predictive models of molecular machines involved in Caenorhabditis elegans early embryogenesis. Nature. 2005, 436 (7052): 861-865. 10.1038/nature03876View ArticlePubMedGoogle Scholar
- Goh KI, Cusick ME, Valle D, Childs B, Vidal M, Barabasi AL: The human disease network. Proc Natl Acad Sci U S A. 2007, 104 (21): 8685-8690. 10.1073/pnas.0701361104PubMed CentralView ArticlePubMedGoogle Scholar
- Kitano H: A robustness-based approach to systems-oriented drug design. Nat Rev Drug Discov. 2007, 6 (3): 202-210. 10.1038/nrd2195View ArticlePubMedGoogle Scholar
- Hauser SL, Oksenberg JR: The neurobiology of multiple sclerosis: genes, inflammation, and neurodegeneration. Neuron. 2006, 52 (1): 61-76. 10.1016/j.neuron.2006.09.011View ArticlePubMedGoogle Scholar
- Oksenberg JR, Hauser SL: Genetics of multiple sclerosis. NeurolClin. 2005, 23 (1): 61-75, vi.Google Scholar
- Bruck W: The pathology of multiple sclerosis is the result of focal inflammatory demyelination with axonal damage. J Neurol. 2005, 252 Suppl 5: v3-v9. 10.1007/s00415-005-5002-7View ArticlePubMedGoogle Scholar
- Imitola J, Chitnis T, Khoury SJ: Insights into the molecular pathogenesis of progression in multiple sclerosis: potential implications for future therapies. Arch Neurol. 2006, 63 (1): 25-33. 10.1001/archneur.63.1.25View ArticlePubMedGoogle Scholar
- Fernald GH, Yeh RF, Hauser SL, Oksenberg JR, Baranzini SE: Mapping gene activity in complex disorders: Integration of expression and genomic scans for multiple sclerosis. J Neuroimmunol. 2005, 167 (1-2): 157-169. 10.1016/j.jneuroim.2005.06.032View ArticlePubMedGoogle Scholar
- Mattson MP: Pathways towards and away from Alzheimer's disease. Nature. 2004, 430 (7000): 631-639. 10.1038/nature02621PubMed CentralView ArticlePubMedGoogle Scholar
- Cummings JL: Alzheimer's disease. N Engl J Med. 2004, 351 (1): 56-67. 10.1056/NEJMra040223View ArticlePubMedGoogle Scholar
- Palop JJ, Chin J, Mucke L: A network dysfunction perspective on neurodegenerative diseases. Nature. 2006, 443 (7113): 768-773. 10.1038/nature05289View ArticlePubMedGoogle Scholar
- Ray S, Britschgi M, Herbert C, Takeda-Uchimura Y, Boxer A, Blennow K, Friedman LF, Galasko DR, Jutel M, Karydas A, Kaye JA, Leszek J, Miller BL, Minthon L, Quinn JF, Rabinovici GD, Robinson WH, Sabbagh MN, So YT, Sparks DL, Tabaton M, Tinklenberg J, Yesavage JA, Tibshirani R, Wyss-Coray T: Classification and prediction of clinical Alzheimer's diagnosis based on plasma signaling proteins. Nat Med. 2007, 13 (11): 1359-1362. 10.1038/nm1653View ArticlePubMedGoogle Scholar
- von Mering C, Jensen LJ, Snel B, Hooper SD, Krupp M, Foglierini M, Jouffre N, Huynen MA, Bork P: STRING: known and predicted protein-protein associations, integrated and transferred across organisms. Nucleic Acids Res. 2005, 1 (33 (Database issue)): D433-D437.Google Scholar
- Hormozdiari F, Berenbrink P, Pr Ulj N, Sahinalp SC: Not All Scale-Free Networks Are Born Equal: The Role of the Seed Graph in PPI Network Evolution. PLoS Comput Biol. 2007, 3 (7): e118- 10.1371/journal.pcbi.0030118PubMed CentralView ArticlePubMedGoogle Scholar
- Han JD, Dupuy D, Bertin N, Cusick ME, Vidal M: Effect of sampling on topology predictions of protein-protein interaction networks. Nat Biotechnol. 2005, 23 (7): 839-844. 10.1038/nbt1116View ArticlePubMedGoogle Scholar
- Jeong H, Mason SP, Barabasi AL, Oltvai ZN: Lethality and centrality in protein networks. Nature. 2001, 411 (6833): 41-42. 10.1038/35075138View ArticlePubMedGoogle Scholar
- Yu H, Kim PM, Sprecher E, Trifonov V, Gerstein M: The importance of bottlenecks in protein networks: correlation with gene essentiality and expression dynamics. PLoS Comput Biol. 2007, 3 (4): e59- 10.1371/journal.pcbi.0030059PubMed CentralView ArticlePubMedGoogle Scholar
- Joy MP, Brock A, Ingber DE, Huang S: High-betweenness proteins in the yeast protein interaction network. J Biomed Biotechnol. 2005, 2005 (2): 96-103. 10.1155/JBB.2005.96PubMed CentralView ArticlePubMedGoogle Scholar
- Lu X, Jain VV, Finn PW, Perkins DL: Hubs in biological interaction networks exhibit low changes in expression in experimental asthma. Mol Syst Biol. 2007, 3: 98- 10.1038/msb4100138PubMed CentralView ArticlePubMedGoogle Scholar
- Barabasi AL, Oltvai ZN: Network biology: understanding the cell's functional organization. NatRevGenet. 2004, 5 (2): 101-113.Google Scholar
- He X, Zhang J: Why do hubs tend to be essential in protein networks?. PLoS Genet. 2006, 2 (6): e88- 10.1371/journal.pgen.0020088PubMed CentralView ArticlePubMedGoogle Scholar
- Kitano H: Biological robustness. Nat Rev Genet. 2004, 5 (11): 826-837. 10.1038/nrg1471View ArticlePubMedGoogle Scholar
- Palacios R, Goni J, Martinez-Forero I, Iranzo J, Sepulcre J, Melero I, Villoslada P: A network analysis of the human T-cell activation gene network identifies JAGGED1 as a therapeutic target for autoimmune diseases. PLoS ONE. 2007, 2 (11): e1222- 10.1371/journal.pone.0001222PubMed CentralView ArticlePubMedGoogle Scholar
- Lim J, Hao T, Shaw C, Patel AJ, Szabo G, Rual JF, Fisk CJ, Li N, Smolyar A, Hill DE, Barabasi AL, Vidal M, Zoghbi HY: A protein-protein interaction network for human inherited ataxias and disorders of Purkinje cell degeneration. Cell. 2006, 125 (4): 801-814. 10.1016/j.cell.2006.03.032View ArticlePubMedGoogle Scholar
- Limviphuvadh V, Tanaka S, Goto S, Ueda K, Kanehisa M: The commonality of protein interaction networks determined in neurodegenerative disorders (NDDs). Bioinformatics. 2007, 23 (16): 2129-2138. 10.1093/bioinformatics/btm307View ArticlePubMedGoogle Scholar
- Bomprezzi R, Ringner M, Kim S, Bittner ML, Khan J, Chen Y, Elkahloun A, Yu AJ, Bielekova B, Meltzer PS, Martin R, McFarland HF, Trent JM: Gene expression profile in multiple sclerosis patients and healthy controls: identifying pathways relevant to disease. Hum Mol Genet. 2003, 12 (17): 2191-2199. 10.1093/hmg/ddg221View ArticlePubMedGoogle Scholar
- Lindberg RL, De Groot CJ, Certa U, Ravid R, Hoffmann F, Kappos L, Leppert D: Multiple sclerosis as a generalized CNS disease--comparative microarray analysis of normal appearing white matter and lesions in secondary progressive MS. J Neuroimmunol. 2004, 152 (1-2): 154-167. 10.1016/j.jneuroim.2004.03.011View ArticlePubMedGoogle Scholar
- Maes OC, Xu S, Yu B, Chertkow HM, Wang E, Schipper HM: Transcriptional profiling of Alzheimer blood mononuclear cells by microarray. Neurobiol Aging. 2007, 28 (12): 1795-1809. 10.1016/j.neurobiolaging.2006.08.004View ArticlePubMedGoogle Scholar
- Walker PR, Smith B, Liu QY, Famili AF, Valdes JJ, Liu Z, Lach B: Data mining of gene expression changes in Alzheimer brain. Artificial Intelligence in Medicine. 2004, 31 (2): 137-154. 10.1016/j.artmed.2004.01.008View ArticlePubMedGoogle Scholar
- Hsu CW, Juan HF, Huang HC: Characterization of microRNA-Regulated Protein-Protein Interaction Network . Proc ICSB. 2007Google Scholar
- ExPlain 2.3 Tool ., http://www.biobase-international.com
- Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I: Controlling the false discovery rate in behavior genetics research. Behav Brain Res. 2001, 125 (1-2): 279-284. 10.1016/S0166-4328(01)00297-2View ArticlePubMedGoogle Scholar