- Open Access
Unraveling human protein interaction networks underlying co-occurrences of diseases and pathological conditions
Journal of Translational Medicine volume 12, Article number: 99 (2014)
Human diseases frequently cause complications such as obesity-induced diabetes and share numbers of pathological conditions, such as inflammation, by dysfunctions of common functional modules, such as protein–protein interactions (PPIs).
Our developed pipeline, ICod (Interaction analysis for disease Comorbidity), grades similarities between pairs of disease-related PPIs including comorbid diseases and pathological conditions. ICod displayed a disease similarity network consisting of nodes of disease PPIs and edges of similarity value. As a proof of concept, eight complex diseases and pathological conditions, such as type 2 diabetes, obesity, inflammation, and cancers, were examined to discover whether PPIs shared between diseases were associated with comorbidities.
By comparing Medicare reports of disease co-occurrences from 31 million patients, the disease similarity network shows that PPIs of pathological conditions, including insulin resistance, and inflammation, overlap significantly with PPIs of various comorbid diseases, including diabetes, obesity, and cancers (p < 0.05). Interestingly, maintaining connectivity between essential genes was more drastically perturbed by removing a node of a disease-related gene rather than a pathological condition-related gene, such as one related to inflammations.
Thus, PPIs of pathological symptoms are underlying functional modules across diseases accompanying comorbidity phenomena, whereas they contribute only marginally to maintaining interactions between essential genes.
Most diseases are the result of the collapse of cellular processes together with interaction networks among components of the genome, proteome, and metabolome, and these perturbed components are likely to be linked with other diseases . Indeed, disease comorbidities such that the onset of one disease increases the likelihood of the development of other diseases were correlated with the breakdown of common functional modules of disease pairs, such as metabolic and cellular networks [2, 3]. Therefore, exploring the biological network between diseases, such as protein–protein interactions (PPIs) of chronic diseases and complications, might give us a more detailed understanding of disease comorbidity and the functional differences between complex diseases.
A number of previous attempts at “network analysis” of diseases have revolutionized our knowledge about the relationships between human diseases and comorbidity [1, 2, 4–6]. For instance, disease-related genetic mutations of genes tend to be peripheral nodes of the essential network, while somatic mutations of genes related to cancers were central nodes . However, pathological phenotypes linked with comorbid diseases and complications remain unclear in the graph-theoretic frame. While distinct diseases share pathological symptoms and various comorbidity patterns, such as inflammations commonly associated with obesity and diabetes, a network model to depict sharing of conditions between diseases remains uncertain. In addition, network models to portray differences between diseases and pathological symptoms leading to severe (or minor) abnormalities of vital functions have been scarcely addressed, whereas distinct mortality issues have been highlighted among cancer-like diseases and pathological symptoms [7, 8].
Here, we designed a novel method, ICod, to build similarity networks among PPIs of disease and pathological conditions to address relationships between comorbid disease pairs and pathological symptoms. While there are various patterns of disease comorbidity, we focused obesity as one of leading risk factors contributing to the overall burden of disease worldwide . Among various obesity related complications, we selected seven diseases and pathological symptoms, which have been remarked as obesity related diseases and manifestations [9, 10]. Thus, the disease studied are obesity, type 2 diabetes mellitus (T2DM), breast cancer, colon cancer, and prostate cancer, and pathological symptoms are inflammation, insulin resistance, and immune response. The main assumption of ICod is that dysfunctions of common protein interactions between diseases might lead to disease comorbidities. To evaluate phenomic associations between network similarities of diseases and comorbidity patterns, disease co-occurrences in a human population were also interrogated using onset co-occurrence relationships based on 31 million patients  (http://hudine.neu.edu/). Furthermore, we address the structural importance of disease- and pathological condition-related genes in maintaining connectivity in the network of essential genes to suggest distinct network models for the dysfunction degrees under diseases or pathological symptoms including inflammation. The attack tolerance of the essential network was determined by measuring alterations of network diameter following removal of disease- and pathological symptom-related essential genes, respectively. The network diameter, defined as the average length of the shortest paths between any two nodes in a network, represents the ability to communicate between any two nodes within the network .
Materials and methods
ICod: Similarity of disease- and pathological condition-related PPIs
ICod used the following measures to grade the disease–disease similarities based on disease-related PPIs. The PPI network modules of each disease were explored using the disease-related genes in our datasets. A disease-related network was produced based on the first neighboring nodes of each disease-related gene. The distance between each pair of disease-related genes was calculated based on the shortest path in the integrated human PPI network. The distances were normalized by transforming them using the formula of Perlman et al.. The transformed distances of all the protein pairs in the two disease-related networks were used to compute the proportion of overlap by considering the size of each corresponding network. The detailed equation is:
where NET i and NET j denote the networks related to diseases i and j, respectively, D(p n , p m ) is the length of the shortest path between protein p n and p m , and S(p n , p m ) are the transformed distances between the networks based on the definition provided by Perlman et al.:
We used A = 0.9 and b = 1, as recommended by Perlman et al.. C is the threshold of D(p n , p m ), which indicates sufficient proximity between two proteins. We used C = 0 to consider directly overlapping proteins in two disease-related networks. Thus, μ(NET i , NET j ) represents the normalized proportion of the overlap based on the overall size of the networks. The statistical significance of μ(NET i , NET j ) was measured as the p value based on the background distribution of μ in 1000 randomly permuted tests. With identical manner, we also determined similarities between pathological conditions.
Preparation of datasets
List of seed genes
We used 317 genes related to five diseases (obesity, T2DM, breast cancer, prostate cancer, and colon cancer) and three pathological conditions (inflammation, immune response, and insulin resistance). These genes were collected from the public resource, GeneCards , which was searched by using related keywords such as “breast cancer”, “malignant neoplasm of breast”, “T2DM,” and “insulin resistance” (Additional file 1: Table S1). All of disease related keywords were manually selected from the results of concept ID search on the largest biomedical terminology database, UMLS (Unified Medical Language System) . In case of immune response, we combined results of immune disorder to comprising immune response related disease symptoms.
Protein–protein interaction (PPI) network
We integrated various well-known resources to prepare human PPI networks: the Human Protein Reference Database (HPRD) ; BioGrid (the Biological General Repository for Interaction Datasets) ; IntAct ; the Molecular INTeraction database (MINT) ; and the Database of Interacting Proteins (DIP) . To produce valid PPI networks, we only used protein interactions with physical evidence; i.e., those with Proteomics Standard Initiative —Molecular Interactions (PSI-MI) codes, such as physical interactions (MI: 0218), direct interactions (MI: 0407), and physical associations (MI: 0915).
Overview of ICod pipeline
The starting point of ICod is the generation of seed genes related to diseases or pathological conditions. A search of GeneCards was made using 33 keywords, and a total of 317 genes related to five diseases (obesity, T2DM, prostate cancer, breast cancer, and colon cancer) and three pathological symptoms (inflammation, insulin resistance, and immune response) were prepared  (Additional file 2: Table S2). Using these seed genes, subnetworks of diseases of interest and pathological conditions were prepared by exploring nearest neighbors among integrated human PPIs over 11 K of nodes and 113 K of edges from five different databases [16–20] (Figure 1A). Then, ICod computes disease–disease similarities considering degrees of overlap between disease- or pathological condition-related PPIs to produce a disease similarity network graphically, as depicted in Figure 1B. The disease similarity network consists of nodes for disease PPIs and edges for degree of similarity between the i- th and j- th PPIs (μ(d i , d j )). The detailed equations for μ(d i , d j ) are described in the Methods section. Figure 1C shows an example ICod pipeline to determine similarity, μ, between obesity- and T2DM-related PPIs after building subnetworks consisting of 3,969 obesity genes and 3,760 T2DM genes. For the statistical significance, p value of μ(d obesity , d T2DM ) was assessed using random permutation tests.
PPI similarity and comorbidity patterns among diseases and pathological conditions
Figure 2A shows the p values of the disease–disease PPI similarity (μ(d i , d j )) matrix calculated by ICod. Housekeeping and essential gene networks, as controls for nondisease networks, were built from 2,164 seed genes and their neighbor nodes (2064 for housekeeping and 115 for essential genes) from previous attempts [22, 23]. Figure 2B shows the background distribution of similarity values to compute the p values of the disease–disease similarity value, μ. Based on these permutation test results, inflammation, insulin resistance, and immune response show high similarity with multiple diseases (p < 0.05, Figure 2A).
As shown in Figure 2C, the similarity between obesity and T2DM is significantly high (p = 2.27E–03) and pathological conditions significantly overlapped with various diseases. Disease- and pathological condition - PPIs, except inflammation, are significantly similar to the essential network (p < 0.05). Thus, insulin resistance- and immune response-related PPIs were commonly incorporated in various disease and essential gene networks.
Figure 2D depicts onset co-occurrence of disease groups interrogated from previous attempts utilizing the medical records of 31 million patients . Using the network frame, we displayed statistical significances of disease co-occurrences (i.e., comorbidities) based on relative risk values . In the presented comorbidity network, nodes mean disease onsets and edges are relative risk values between nodes. As displayed in Figure 2D, inflammation-related symptoms (gray nodes) are closely associated with the onset of obesity (yellow nodes) and T2DM (green nodes).
As Figure 2C and D depict, comorbidities showed similar tendencies to PPI similarity networks of diseases. The network analysis presented supports the hypothesis that the collapse of common PPIs between disease and pathological symptoms were closely associated with comorbidity.
Topological role of disease- and pathological condition-related genes for essential interactions
Irrespective of mortality rate, pathological symptom-related PPIs overlap significantly with various diseases including cancers and even networks of essential genes. Using the measure of node degree (i.e., number of nearest neighbors), a previous network-based attempt suggested that disease-related genes were peripheral nodes in the essential network . Nevertheless, cancer-like diseases are major causes of death in the world , although models showing a severe impact on essential gene networks remain imprecise. Here, except for node degree, we compared topological importance between disease- and pathological condition-related nodes through measuring the collapsed connectivity of the essential network by elimination of disease- or symptom-related essential genes.
By the definition of network diameter (mean of shortest paths between all pairs of nodes in a network) , a pattern of increasing the diameter by removing a node denotes the breakage of network links and the vital role of the removed node in maintaining information flows in the network. According to Figure 3A and B, the network diameters of disease PPIs, pathological condition PPIs, housekeeping PPIs, and essential PPIs indicate a diverse density of connectivity within each network unrelated to the number of nodes or edges. Because the PPI network of human essential genes is a scale-free network, it is robust against network errors while it is vulnerable against attacks on high-degree nodes (Figure 3C). “Network attack” means sequential removal of nodes according to node rank (i.e., a hub-node attack), whereas “network error” (i.e., random error) means withdrawing a node randomly.
Interestingly, alteration of the essential network diameter by attacks on the pathological condition-related node was negligible (Figure 3D), whereas the connectivity of the essential network collapsed dramatically on removal of a disease-related essential gene (Figure 3E). The connectivity of the essential network (diameter 4.13) was dramatically perturbed even under attack of a small fraction of nodes related to diseases (0.1% of nodes in Figure 3E), such as T2DM, cancers, and obesity. However, attacks on a larger fraction of nodes related to pathological symptoms, including insulin resistance, inflammation, and immune response, showed subtle effects in truncating interactions among essential genes (Figure 3D). Based on these distinct topological roles of disease- and symptom-related nodes in the essential network, we suggest that disease-related nodes are vital nodes of information flow in the essential network, whereas nodes of pathological symptom play a less pivotal role.
In summary, using ICod, we determined relationships among five diseases (prostate cancer, breast cancer, colon cancer, T2DM, and obesity), three pathological conditions (inflammation, insulin resistance, and immune response), and the essential gene network. As expected from pathological symptoms in complex diseases sharing common phenotypic signals including inflammation, the results of ICod support our knowledge at the network level. The pathological condition network is closely associated with various disease networks. Our findings are the first attempt at uncovering the differences in topological role between each disease- and symptom-related network within the essential gene network using analysis of attack tolerance. Although PPIs of pathological conditions significantly overlapped with disease PPIs, the patterns of collapsing the essential network by removing the condition-related nodes were clearly distinct from attacks on disease-related nodes. While our network analysis covered partial sets of human diseases and symptoms, our conceptual approach successfully modeled functional roles of pathological states in disease etiology and maintenance of the essential network. Typical pathological symptoms, such as inflammation and immune responses, are widely spread mechanisms behind complex diseases with subtle impacts that can cause severe dysfunctions of the essential network.
Since our network model focused topological similarity, our method suggested network relationships between disease pairs, or disease-pathological symptoms without causal understandings and functional significance. To address network related functional impact (i.e., complete node removal and partial mutation), Zhong et al. attempted computational and experimental validation using Yeast-Two-Hybrid system (Y2H) . In our previous study, we analyzed gene expression patterns in diet induced obese mice . Interestingly, diet-induced obese mice displayed differentially expressed genes, which were related inflammation, immune response and insulin resistance as we suggested our network similarity analysis. While our previous work suggested enriched functional signatures under induced obesity condition without node-removal effect, significantly depict obesity derived pathological phenotype in time-resolving frame. Based on theses attempts, we suggest an approach combining Zhong et al’ s Y2H and time-resolving frame of ours for further functional understanding. Owing to the utilization of model organisms, Zhang et al. and our mouse data analysis give us limited understandings to depict underlying mechanisms of human diseases. Thus, as we conducted in our previous attempt , large-scale human cohort based analysis might shed light shared genetic and functional features, which lead disease comorbidity.
While cancers have shown high mortality rates , obesity and T2DM have low attributes for viability issues. In stark contrast with our expectation, topological roles for the interactions within the essential gene network are homogeneous between lethal diseases (cancers) and other chronic diseases. Therefore, further study is necessary on the associations between disease mortality and aspects of network structure, such as “bottleneckness” . In addition, our keyword-based approach to preparing disease and pathological symptom related genes were introduced for the proof-of-concept. Thus, it is necessary for advanced validation of disease related genes using various approaches, such as scrutinizing gene expression databases .
As shown in our network analysis of disease PPI similarity and US Medicare data, disease onset and comorbidity are closely associated with the breakage of common functional modules. Complex diseases and pathological conditions share molecular mechanisms such as PPIs, whereas mortalities are heterogeneous. We distinguished network models of complex diseases and pathological conditions using our analysis of attack tolerance of the essential network to find the different impacts on mortality issues.
One of our contributions is computing quantitative degree of overlapping between disease PPIs involving across similarity among diseases and related clinical manifestations considering connectivity of compared PPI pairs (Figure 2). Since ICod utilized public repository of PPI networks and list of disease related genes, our method can be a streamlined route to visualize similarity between diseases and pathological phenotypes, which are associated disease comorbidity . In addition, our network-based disease similarity might present drug targets related various diseases as presented by Suthram et al.. For example, ICod remarks a probability for the repositioning of drugs related pathological symptoms, such as inflammation, for the therapy of PPI overlapped diseases including obesity; The anti-inflammation drug, amlexanox, elevate energy expenditure and produce weigh loss in mice .
Therefore, analysis of ICod network similarity and attack tolerance has successfully modeled existing knowledge for disease comorbidities and the co-occurrence of pathological symptoms, and we have identified perturbation models for disease-related essential genes through the network frame.
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.0701361104.
Lee DS, Park J, Kay KA, Christakis NA, Oltvai ZN, Barabasi AL: The implications of human metabolic network topology for disease comorbidity. Proc Natl Acad Sci U S A. 2008, 105 (29): 9880-9885. 10.1073/pnas.0802208105.
Park J, Lee DS, Christakis NA, Barabasi AL: The impact of cellular networks on disease comorbidity. Mol Syst Biol. 2009, 5: 262-
Lage K, Karlberg EO, Storling ZM, Olason PI, Pedersen AG, Rigina O, Hinsby AM, Tumer Z, Pociot F, Tommerup N, Moreau Y, Brunak S: A human phenome-interactome network of protein complexes implicated in genetic disorders. Nat Biotechnol. 2007, 25 (3): 309-316. 10.1038/nbt1295.
Park S, Yang JS, Shin YE, Park J, Jang SK, Kim S: Protein localization as a principal feature of the etiology and comorbidity of genetic diseases. Mol Syst Biol. 2011, 7: 494-
Park S, Yang JS, Kim J, Shin YE, Hwang J, Park J, Jang SK, Kim S: Evolutionary history of human disease genes reveals phenotypic connections and comorbidity among genetic diseases. Sci Rep. 2012, 2: 757-
Jemal A, Siegel R, Ward E, Hao Y, Xu J, Murray T, Thun MJ: Cancer statistics, 2008. CA Cancer J Clin. 2008, 58 (2): 71-96. 10.3322/CA.2007.0010.
Fuller JH, Elford J, Goldblatt P, Adelstein AM: Diabetes mortality: new light on an underestimated public health problem. Diabetologia. 1983, 24 (5): 336-341.
Haslam DW, James WP: Obesity. Lancet. 2005, 366 (9492): 1197-1209. 10.1016/S0140-6736(05)67483-1.
De Pergola G, Silvestris F: Obesity as a major risk factor for cancer. J Obes. 2013, 2013: 291546-
Hidalgo CA, Blumm N, Barabasi AL, Christakis NA: A dynamic network approach for the study of human phenotypes. PLoS Comput Biol. 2009, 5 (4): e1000353-10.1371/journal.pcbi.1000353.
Albert R, Jeong H, Barabasi AL: Error and attack tolerance of complex networks. Nature. 2000, 406 (6794): 378-382. 10.1038/35019019.
Perlman L, Gottlieb A, Atias N, Ruppin E, Sharan R: Combining drug and gene similarity measures for drug-target elucidation. J Comput Biol. 2011, 18 (2): 133-145. 10.1089/cmb.2010.0213.
Safran M, Solomon I, Shmueli O, Lapidot M, Shen-Orr S, Adato A, Ben-Dor U, Esterman N, Rosen N, Peter I, Olender T, Chalifa-Caspi V, Lancet D: GeneCards 2002: towards a complete, object-oriented, human gene compendium. Bioinformatics. 2002, 18 (11): 1542-1543. 10.1093/bioinformatics/18.11.1542.
Bodenreider O: The Unified Medical Language System (UMLS): integrating biomedical terminology. Nucleic Acids Res. 2004, 32 (Database issue): D267-D270.
Keshava Prasad TS, Goel R, Kandasamy K, Keerthikumar S, Kumar S, Mathivanan S, Telikicherla D, Raju R, Shafreen B, Venugopal A, Balakrishnan L, Marimuthu A, Banerjee S, Somanathan DS, Sebastian A, Rani S, Ray S, Harrys Kishore CJ, Kanth S, Ahmed M, Kashyap MK, Mohmood R, Ramachandra YL, Krishna V, Rahiman BA, Mohan S, Ranganathan P, Ramabadran S, Chaerkady R, Pandey A: Human Protein Reference Database–2009 update. Nucleic Acids Res. 2009, 37 (Database issue): D767-D772.
Chatr-Aryamontri A, Breitkreutz BJ, Heinicke S, Boucher L, Winter A, Stark C, Nixon J, Ramage L, Kolas N, O’Donnell L, Reguly T, Breitkreutz A, Sellam A, Chen D, Chang C, Rust J, Livstone M, Oughtred R, Dolinski K, Tyers M: The BioGRID interaction database: 2013 update. Nucleic Acids Res. 2013, 41 (Database issue): D816-D823.
Kerrien S, Alam-Faruque Y, Aranda B, Bancarz I, Bridge A, Derow C, Dimmer E, Feuermann M, Friedrichsen A, Huntley R, Kohler C, Khadake J, Leroy C, Liban A, Lieftink C, Montecchi-Palazzi L, Orchard S, Risse J, Robbe K, Roechert B, Thorneycroft D, Zhang Y, Apweiler R, Hermjakob H: IntAct–open source resource for molecular interaction data. Nucleic Acids Res. 2007, 35 (Database issue): D561-D565.
Licata L, Briganti L, Peluso D, Perfetto L, Iannuccelli M, Galeota E, Sacco F, Palma A, Nardozza AP, Santonico E, Castagnoli L, Cesareni G: MINT, the molecular interaction database: 2012 update. Nucleic Acids Res. 2012, 40 (Database issue): D857-D861.
Xenarios I, Salwinski L, Duan XJ, Higney P, Kim SM, Eisenberg D: DIP, the Database of Interacting Proteins: a research tool for studying cellular networks of protein interactions. Nucleic Acids Res. 2002, 30 (1): 303-305. 10.1093/nar/30.1.303.
Safran M, Dalah I, Alexander J, Rosen N, Iny Stein T, Shmoish M, Nativ N, Bahir I, Doniger T, Krug H, Sirota-Madi A, Olender T, Golan Y, Stelzer G, Harel A, Lancet D: GeneCards Version 3: the human gene integrator. Database (Oxford). 2010, 2010: baq020-
Zhang CT, Zhang R: Gene essentiality analysis based on DEG, a database of essential genes. Methods Mol Biol. 2008, 416: 391-400. 10.1007/978-1-59745-321-9_27.
Chang CW, Cheng WC, Chen CR, Shu WY, Tsai ML, Huang CL, Hsu IC: Identification of human housekeeping genes and tissue-selective genes by microarray meta-analysis. PLoS One. 2011, 6 (7): e22859-10.1371/journal.pone.0022859.
Zhong Q, Simonis N, Li QR, Charloteaux B, Heuze F, Klitgord N, Tam S, Yu H, Venkatesan K, Mou D, Swearingen V, Yildirim MA, Yan H, Dricot A, Szeto D, Lin C, Hao T, Fan C, Milstein S, Dupuy D, Brasseur R, Hill DE, Cusick ME, Vidal M: Edgetic perturbation models of human inherited disorders. Mol Syst Biol. 2009, 5: 321-
Heo HS, Kim E, Jeon SM, Kwon EY, Shin SK, Paik H, Hur CG, Choi MS: A nutrigenomic framework to identify time-resolving responses of hepatic genes in diet-induced obese mice. Mol Cells. 2013, 36 (1): 25-38. 10.1007/s10059-013-2336-3.
Ban HJ, Kim SC, Seo J, Kang HB, Choi JK: Genetic and metabolic characterization of insomnia. PLoS One. 2011, 6 (4): e18455-10.1371/journal.pone.0018455.
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.0030059.
Sirota M, Dudley JT, Kim J, Chiang AP, Morgan AA, Sweet-Cordero A, Sage J, Butte AJ: Discovery and preclinical validation of drug indications using compendia of public gene expression data. Sci Transl Med. 2011, 3 (96): 96ra77-
Suthram S, Dudley JT, Chiang AP, Chen R, Hastie TJ, Butte AJ: Network-based elucidation of human disease similarities reveals common functional modules enriched for pluripotent drug targets. PLoS Comput Biol. 2010, 6 (2): e1000662-10.1371/journal.pcbi.1000662.
Reilly SM, Chiang SH, Decker SJ, Chang L, Uhm M, Larsen MJ, Rubin JR, Mowers J, White NM, Hochberg I, Downes M, Yu RT, Liddle C, Evans RM, Oh D, Li P, Olefsky JM, Saltiel AR: An inhibitor of the protein kinases TBK1 and IKK-varepsilon improves obesity-related metabolic dysfunctions in mice. Nat Med. 2013, 19 (3): 313-321. 10.1038/nm.3082.
This work was supported by the Research Program funded by the Korea Centers for Disease Control and Prevention (2012-NG72001-00). HP was also supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2012R1A1A3019523) and Lucile Packard Foundation for Children’s Health and National Institute of General Medical Sciences (R01 GM079719) of USA.
There are no conflicts of interest.
HP built the main algorithm, designed the analysis methods, prepared artworks of figures and helped write the manuscript. HSH prepared datasets and helped write the manuscript. HB analyzed co-occurrence patterns of human diseases. SBC oversaw the overall organization of the manuscript as a corresponding author. All authors read and approved the final manuscript.
Hyojung Paik, Hyoung-Sam Heo contributed equally to this work.
About this article
Cite this article
Paik, H., Heo, H., Ban, H. et al. Unraveling human protein interaction networks underlying co-occurrences of diseases and pathological conditions. J Transl Med 12, 99 (2014) doi:10.1186/1479-5876-12-99
- Protein–protein interaction
- Attack tolerance