- Open Access
DisSetSim: an online system for calculating similarity between disease sets
Journal of Biomedical Semantics volume 8, Article number: 28 (2017)
Functional similarity between molecules results in similar phenotypes, such as diseases. Therefore, it is an effective way to reveal the function of molecules based on their induced diseases. However, the lack of a tool for obtaining the similarity score of pair-wise disease sets (SSDS) limits this type of application.
Here, we introduce DisSetSim, an online system to solve this problem in this article. Five state-of-the-art methods involving Resnik’s, Lin’s, Wang’s, PSB, and SemFunSim methods were implemented to measure the similarity score of pair-wise diseases (SSD) first. And then “pair-wise-best pairs-average” (PWBPA) method was implemented to calculated the SSDS by the SSD. The system was applied for calculating the functional similarity of miRNAs based on their induced disease sets. The results were further used to predict potential disease-miRNA relationships.
The high area under the receiver operating characteristic curve AUC (0.9296) based on leave-one-out cross validation shows that the PWBPA method achieves a high true positive rate and a low false positive rate. The system can be accessed from http://www.bio-annotation.cn:8080/DisSetSim/.
The similarity of pair-wise disease sets (SDS) has drawn more and more attention in identifying functional similarity of the disease-caused molecules , predicting potential relationships between diseases and molecules [2,3,4,5,6,7,8], and so on. In previous studies, Wang et al. utilized the SDS to construct a human miRNA functional similarity network (MFSN) . And Sun et al. used the SDS to predict novel disease lncRNA relationships .
The performance of calculating the SDS is mainly based on the method for computing the similarity of pair-wise diseases (SD). Currently, seven state-of-art methods involving Resnik’s , Lin’s , Wang’s , process-similarity based (PSB) , SemFunSim , ILNCSIM , and FMLNCSIM  methods were frequently used for computing the SD. Among these methods, Resnik’s , Lin’s , and Wang’s methods  are designed earlier for Gene Ontology (GO) [8, 17]. And these methods were introduced for calculating the SD by DOSim  and DisSim . Resnik’s and Lin’s methods [10, 11] are based on information content (IC) for computing similarity between terms of ontology. And Wang’s method  is based on the hierarchical structure of the ontology. PSB and SemFunSim methods are newly developed for Disease Ontology (DO) . PSB method  utilized the association of biological process between genes to calculate disease similarity. In comparison, SemFunSim method  considered more types of the functional associations including protein-protein interaction , human mRNA co-expression , and so on.
Resources for calculating the similarity score of pair-wise diseases (SSD) mainly includes the vocabularies of diseases and disease-related genes. The frequently used disease vocabularies contain Online Mendelian Inheritance in Man (OMIM) , Medical Subject Headings (MeSH) , and DO . OMIM records the names of genetic disorders without providing semantic associations between them. MeSH provides a hierarchy of terms in biomedical domain. It contains 16 categories, of which only C and F03 involve disease names. In comparison with OMIM and MeSH, DO has been established around the concept of disease, and it aims to provide a clear definition for each disease. The disease-related genes are scattered in the databases, such as Gene Reference into Function (GeneRIF) , OMIM , Genetic Association Database (GAD)  and Comparative Toxicogenomics Database (CTD) . It is better to use relationships of all of these databases.
“pair-wise-all pairs-maximum” (PWAPM) method and “pair-wise-best pairs-average” (PWBPA) method are optional for calculating similarity of pair-wise term sets . For comparing multiple aspects, the best measure is the PWBPA method, which is widely utilized in calculating similarity of DO and GO term sets [1, 7, 9, 12].
Although DOSim  and DisSim  implemented the disease similarity methods in R package and web interface, no tools provided the function to calculate the similarity score of pair-wise disease sets (SSDS) currently. In this article, we designed and implemented an online tool DisSetSim to calculate the SSDS. Five state-of-art disease similarity methods (Resnik’s, Lin’s, Wang’s, PSB, and SemFunSim) and the PWBPA method was implemented in the tool. The system is freely available at http://www.bio-annotation.cn:8080/DisSetSim/.
Data sets of DisSetSim are from open source databases, and they are listed in Table 1. DO  records disease names. It provides terms for calculating disease similarity. GeneRIF , OMIM , GAD  and CTD  are manually curated databases of disease-related genes. All of diseases in these databases are mapped to terms in DO based on SIDD . GO annotation (GOA)  includes functional annotation of genes. HumanNet is the gene functional network of human. In addition, HMDD v2.0  contains disease-related miRNAs, diseases of which were manually mapped to terms in DO by OAHG .
Methods for calculating similarity score of pair-wise diseases
Five state-of-art methods involving Resnik’s , Lin’s , Wang’s , PSB , and SemFunSim methods  have been implemented for calculating the SSD.
Resnik’s and Lin’s methods are based on IC. The IC of a disease t is described as Eq. 1:
where N is the total number of genes annotated by diseases, and n t is the number of genes annotated by t. Assuming t 1 and t 2 are two diseases, the similarity of them is defined by Resnik as following :
where t MICA is the most informative common ancestor (MICA) of t 1 and t 2 . Lin defines the similarity of t 1 and t 2 as Eq. 3 :
Assuming T 1 is the set involving t 1 and all of its ancestor terms of ontology. Semantic contribution of term t to t 1 is represented as following:
where w is the contribution factor of each semantic relationship. According to Wang et al. , w is defined as 0.5 for ‘IS_A’ relationship of DO . Then, all the semantic contributions of T 1 to t 1 is SV(t 1 ), which is defined as following:
Assuming T 2 is the set involving t 2 and all of its ancestor terms, the similarity between t 1 and t 2 is defined as following by Wang’s method :
Assuming t 1 and t 2 can be related with m and n biological processes of GO based on hypergeometric test, respectively, the similarity of t 1 and t 2 is defined by the PSB method as following:
where p 1i and p 2j is the ith and jth significant related biological process terms of t 1 and t 2 , respectively. Sim(p 1i , p 2j ) represents similarity between two processes p 1i and p 2j , which is defined as Eq. 8:
where IC GO and IC DO represent IC based on GO and DO, respectively. n(p 1 ∩ p 2 ) and n(p 1 ∪ p 2 ) denote the number of common genes of p 1 and p 2 , and the number of total genes of p 1 and p 2 , respectively.
Assuming G 1 and G 2 represent related gene sets of t 1 and t 2 , respectively. Then, the similarity of t 1 and t 2 by the SemFunSim method can be described as following:
where |GMICA| represents the number of genes in GMICA. m and n denote the number of genes in G 1 and G 2 , respectively. Sim(g 1i , g 2j ) is the functional similarity score between genes g 1i and g 2j , which could be obtained from HumanNet .
Method for calculating similarity score of pair-wise disease sets
The PWBPA method was utilized for calculating the SSDS. The similarity of two disease sets T 1 and T 2 is defined as following:
where T 1 and T 2 contains N and M diseases, respectively. t i and t j represents ith and jth terms of T 1 and T 2 , respectively.
Predicting potential association between diseases and miRNAs
Functional similarity between miRNAs could be calculated based on their related disease sets. Similarities of each pair-wise miRNAs are utilized to establish a MFSN. Node of the network represents miRNA. Weight of edge is the functional similarity score. Then, disease-related miRNAs were prioritized using the network ranking algorithm named random walk with restart (RWR) .
The random walker starts on one or several seed nodes and then randomly transits to neighboring nodes considering the probabilities of the edges between the two nodes. And the probability to return to the seed nodes is supposed as γ. Then, RWR algorithm can be defined as following:
where P 0 denotes the initial probability vector, P t is a vector in which the ith element represents the probability of finding the walker at node i and step t, A is the column-normalized adjacency matrix of the network. The algorithm was performed until the difference between P t and P t+1 falling below 10−10, which means all the nodes become stable.
In this study, the known miRNAs of a disease were considered as seed nodes. The unknown miRNAs of it could be scored based on RWR on the MFSN. After ranking the miRNAs based on the scores, disease-related miRNAs could be prioritized.
DisSetSim has been implemented on a JavaEE framework and run on the web server (2-core (2.26 GHz) processors) of Ucloud . The four-layer architecture involving DATABASE, ALGORITHM, TOOLS, and VIEW layer is shown in Fig. 1 The detailed description of the architecture is fixed as following.
(1) DATABASE layer. This layer stores DO, disease-related genes, and functional associations between genes. These are exploited by ALGORITHM layer for calculating the similarity between disease sets.
(2) ALGORITHM layer. Five algorithms of measuring the similarity between DO terms have been implemented, which include Resnik’s, Lin’s, Wang’s, PSB, and SemFunSim methods. And the method named PWBPA for calculating the SSDS were also implemented.
(3) TOOL layer. Two tools including PairSim and BatchSim have been provided for exploring the SSDS. PairSim calculates the similarity for a given pair of disease sets, and BatchSim computes similarity between each pair of multiple disease sets.
(4) VIEW layer. Web pages are provided for viewing the results. It shows the similarity of pair-wise disease sets.
DisSetSim provides two tools PairSim and BatchSim for querying the SSDS. The details about the usage of these two tools are described as follows.
The usage of PairSim
Figure 2a shows a case for searching the similarity score of a given pair of disease sets. The web page for inputting disease sets is http://www.bio-annotation.cn:8080/DisSetSim/ basic-init. Each of these disease sets could be inputted in a textbox. A disease set is comprised by several diseases. And each disease is represented by the identifier of term in DO. All the term identifiers could be downloaded from the hyperlink ‘disease terms’ in the inputting page. Here, we click the ‘example’ button to use our example. Then, we choose one of the five methods (Resnik’s, Lin’s, Wang’s, PSB, and SemFunSim) for calculating the SSD. After submitting this pair of disease sets, the system could return the similarity score based on the PWBPA method.
The usage of BatchSim
Figure 2b shows a case for searching the similarity score of all the pairs based on the selected files. The web page for inputting disease sets is http://www.bio-annotation.cn:8080/DisSetSim/batch-init. Two files including disease sets should be selected before submitting. The file should be a plain text which contains several disease sets. Each disease set must be in a newline, and each disease set contains several disease IDs which are separated by commas. The size of uploaded file must be <2 Mb. Here, we selected our example file in this page. Then, we choose one of the five methods for calculating the SSD. After clicking the ‘submit’ button, the system could return the similarity score of all the pairs of the selected files based on the PWBPA method.
miRNA functional similarity network
By applying DisSetSim to the inputted disease sets of miRNAs, the similarity score of each miRNA pair could be obtained. Using miRNA as node and similar miRNAs as edge, the MFSN was constructed based on various similarity cutoffs. As shown in Fig. 3a, the number of links dramatically decreases when the cutoff increases from low value to high value. When the cutoff is equal to or bigger than 0.7, the link numbers remain relatively stable. Therefore, we use 0.7 as cutoff for the MFSN. In total, 1042 miRNA-miRNA functional associations between 346 miRNAs were obtained as MFSN (Fig. 3c). Similar to the most of the reported biological networks, the degree of this MFSN also shows a scale-free distribution [5, 9, 35,36,37]. It means that most of the miRNAs only have a few functionally similar miRNAs, and a few of miRNAs have a numerous functional similar miRNA (Fig. 3b).
Here, the PWBPA method was utilized for calculating similarity between disease sets, and SemFunSim was used as computing the similarity of pair-wise diseases. This is because that the SemFunSim method was proven to obtain the best performance . Alternatively, other state-of-art methods could also be chosen to construct MFSN.
By applying the above similarity scores of miRNAs, novel disease-related miRNAs were predicted based on RWR algorithm (See ‘Methods’ section). To evaluate the performance of the similarity scores of miRNAs, leave-one-out cross validation of 5710 known experimentally confirmed miRNA-disease associations, including 265 diseases with at least two miRNAs, were used for this assessment. For a disease of interest, each known miRNA of this disease was left out as the testing case, and the remaining miRNAs of this disease were used as seed nodes. All the miRNAs except the miRNAs of this disease were considered as candidate miRNAs. We then examined how well the testing miRNA ranked relative to the candidate miRNAs. If the ranking of this testing miRNA exceeded a given cutoff, we regarded this miRNA-disease association as successfully predicted. As a result, an area under the ROC curve (AUC) of 0.9296 was achieved (Fig. 4), which demonstrated that our miRNA functional similarity was effective in recovering known experimentally confirmed disease-related miRNAs.
As the best of our knowledge, non-coding RNAs (ncRNAs) attract more and more attentions because of their important regulation roles in molecular level. However, the lack of protein limits the identification of their function. Here the application of our tool in constructing MFSN and predicting miRNA-disease associations provides a novel way to help for exploring the function of miRNAs especially for prioritizing miRNA-disease associations. This application can be extended to other ncRNAs, such as lncRNAs and circRNAs. Although methods for calculating the SDS have been implemented by previous methods, it is not easy to calculate the SSDS. Therefore, DisSetSim benefits researchers for exploring the function of disease-related molecular.
In this article, we designed and developed a web system DisSetSim to calculate the SSDS. Five state-of-art methods were implemented (see ‘METHODS’ section) for calculating disease similarity. And the PWBPA method was implemented for calculating the SSDS. Two tools involving PairSim and BatchSim provide the function to obtain the SSDS by inputting a pair-wise disease sets and multiple disease sets, respectively.
The functional similarity of miRNAs could be calculated based on our system. Here, the similarity of each pair-wise miRNAs was calculated. And then a MFSN was constructed based on miRNA similarity. The network was further utilized to predicate disease-related miRNAs based on RWR. The high AUC (0.9296) shows the MFSN is very suitable for predicting potential relationships between diseases and miRNAs.
Comparative Toxicogenomics Database
Genetic Association Database
Gene Reference into Function
miRNA functional similarity network
the most informative common ancestor
Online Mendelian Inheritance in Man
random walk with restart
the similarity score of pair-wise diseases
the similarity score of pair-wise disease sets.
Wang D, Wang J, Lu M, Song F, Cui Q. Inferring the human microRNA functional similarity and functional network based on microRNA-associated diseases. Bioinformatics. 2010;26(13):1644–50.
Chen X, Yan CC, Luo C, Ji W, Zhang Y, Dai Q. Constructing lncRNA functional similarity network based on lncRNA-disease associations and disease semantic similarity. Sci Rep. 2015;5:11338.
Chen X, Yan CC, Zhang X, You Z-H. Long non-coding RNAs and complex diseases: from experimental results to computational models. Brief Bioinform. 2016:bbw060.
Chen X, You Z, Yan G, Gong D. IRWRLDA: improved random walk with restart for LncRNA-disease association prediction. Oncotarget. 2016;7(36):57919–31.
Peng J, Wang T, Hu J, Wang Y, Chen J. Constructing networks of organelle functional modules in Arabidopsis. Current Genomics. 2016;17(5):427–38.
J. Peng, Bai, K., Shang, X., Wang, G., Xue, H., Jin, S., Cheng, L., Wang, Y., & Chen, J., “Predicting Disease-related Genes using Integrated Biomedical Networks,” BMC Genomics, vol. 18, no. 1, pp. 1043, 2017.
J. Peng, H. Li, Y. Liu, L. Juan, Q. Jiang, Y. Wang, and J. Chen, “InteGO2: a web tool for measuring and visualizing gene semantic similarities using Gene Ontology,” BMC Genomics, vol. 17 Suppl 5, pp. 530, 2016.
J. Peng, T. Wang, J. Wang, Y. Wang, and J. Chen, “Extending gene ontology with gene association networks,” Bioinformatics, vol. 32, no. 8, pp. 1185–94, Apr 15, 2016.
Sun J, Shi H, Wang Z, Zhang C, Liu L, Wang L, He W, Hao D, Liu S, Zhou M. Inferring novel lncRNA-disease associations based on a random walk model of a lncRNA functional similarity network. Mol BioSyst. 2014;10(8):2074–81.
P. Resnik, “Using information content to evaluate semantic similarity in a taxonomy,” arXiv preprint cmp-lg/9511007, 1995.
Lin D. An information-theoretic definition of similarity. 1998. pp. 296–304.
J. Z. Wang, Z. Du, R. Payattakool, P. S. Yu, and C. F. Chen, “A new method to measure the semantic similarity of GO terms,” Bioinformatics, vol. 23, no. 10, pp. 1274–81, May 15, 2007.
Mathur S, Dinakarpandian D. Finding disease similarity based on implicit semantic similarity. J Biomed Inform. 2012;45(2):363–71.
L. Cheng, J. Li, P. Ju, J. Peng, and Y. Wang, “SemFunSim: a new method for measuring disease similarity by integrating semantic and gene functional association,” PLoS One, vol. 9, no. 6, pp. e99415, 2014.
Huang Y, Chen X, You Z, Huang D, Chan K. ILNCSIM: improved lncRNA functional similarity calculation model. Oncotarget. 2016;7(18):25902–14.
Chen X, Huang Y, Wang X, You Z, Chan K. FMLNCSIM: fuzzy measure-based lncRNA functional similarity calculation model. Oncotarget. 2016;7(29):45948–58.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25(1):25–9.
Li J, Gong B, Chen X, Liu T, Wu C, Zhang F, Li C, Li X, Rao S, Li X. DOSim: an R package for similarity between diseases based on disease ontology. BMC Bioinformatics. 2011;12:266.
L. Cheng, Y. Jiang, Z. Wang, H. Shi, J. Sun, H. Yang, S. Zhang, Y. Hu, and M. Zhou, “DisSim: an online system for exploring significant similar diseases and exhibiting potential therapeutic drugs,” Sci Rep, vol. 6, pp. 30024, Jul 26, 2016.
Schriml LM, Arze C, Nadendla S, Chang YW, Mazaitis M, Felix V, Feng G, Kibbe WA. Disease ontology: a backbone for disease semantic integration. Nucleic Acids Res. 2012;40(Database issue):D940–6.
Ortutay C, Vihinen M. Identification of candidate disease genes by integrating gene Ontologies and protein-interaction networks: case study of primary immunodeficiencies. Nucleic Acids Res. Feb, 2009;37(2):622–8.
J. M. Stuart, E. Segal, D. Koller, and S. K. Kim, “A gene-coexpression network for global discovery of conserved genetic modules,” Science, vol. 302, no. 5643, pp. 249–55, Oct 10, 2003.
Amberger J, Bocchini C, Hamosh A. A new face and new challenges for online Mendelian inheritance in man (OMIM®). Hum Mutat. 2011;32(5):564–7.
Dhammi IK, Kumar S. Medical subject headings (MeSH) terms. Indian J Orthop. 2014;48(5):443–4.
Mitchell JA, Aronson AR, Mork JG, Folk LC, Humphrey SM, Ward JM. Gene indexing: characterization and analysis of NLM's GeneRIFs. AMIA Annu Symp Proc. 2003:460–4.
Becker KG, Barnes KC, Bright TJ, Wang SA. The genetic association database. Nat Genet. May, 2004;36(5):431–2.
Davis AP, Murphy CG, Johnson R, Lay JM, Lennon-Hopkins K, Saraceni-Richards C, Sciaky D, King BL, Rosenstein MC, Wiegers TC, Mattingly CJ. The comparative Toxicogenomics database: update 2013. Nucleic Acids Res. Jan, 2013;41(Database issue):D1104–14.
C. Pesquita, D. Faria, A. O. Falcao, P. Lord, and F. M. Couto, “Semantic similarity in biomedical ontologies,” PLoS Comput Biol, vol. 5, no. 7, pp. e1000443, 2009.
L. Cheng, G. Wang, J. Li, T. Zhang, P. Xu, and Y. Wang, “SIDD: A Semantically Integrated Database towards a Global View of Human Disease,” PLOS ONE, vol. 8, no. 10, pp. e75504, 2013.
Camon E, Magrane M, Barrell D, Lee V, Dimmer E, Maslen J, Binns D, Harte N, Lopez R, Apweiler R. The gene ontology annotation (goa) database: sharing knowledge in uniprot with gene ontology. Nucleic Acids Res. 2004;32(suppl 1):D262–6.
Li Y, Qiu C, Tu J, Geng B, Yang J, Jiang T, Cui Q. HMDD v2.0: a database for experimentally supported human microRNA and disease associations. Nucleic Acids Res. 2014;42(Database issue):D1070–4.
Cheng L, Sun J, Xu W, Dong L, Hu Y, Zhou M. OAHG: an integrated resource for annotating human genes with multi-level ontologies. Sci Rep. 2016;10:34820.
Lee I, Blom UM, Wang PI, Shim JE, Marcotte EM. Prioritizing candidate disease genes by network-based boosting of genome-wide association data. Genome Res. 2011;21(7):1109–21.
Sqalli MH, Al-Saeedi M, Binbeshr F, Siddiqui M. UCloud: A simulated Hybrid Cloud for a university environment. 2012. pp. 170–172.
K. I. Goh, M. E. Cusick, D. Valle, B. Childs, M. Vidal, and A. L. Barabasi, “The human disease network,” Proc Natl Acad Sci U S A, vol. 104, no. 21, pp. 8685–90, May 22, 2007.
Zhang F, Ren C, Lau KK, Zheng Z, Lu G, Yi Z, Zhao Y, Su F, Zhang S, Zhang B, Sobie EA, Zhang W, Walsh MJ. A network medicine approach to build a comprehensive atlas for the prognosis of human cancer. Brief Bioinform, Aug. 2016;24
F. Zhang, B. Gao, L. Xu, C. Li, D. Hao, S. Zhang, M. Zhou, F. Su, X. Chen, H. Zhi, and X. Li, “Allele-specific behavior of molecular networks: understanding small-molecule drug response in yeast,” PLoS One, vol. 8, no. 1, pp. e53581, 2013.
Yadong Wang and Liang Cheng are the corresponding authors. Lingling Zhao and Zhiyan Liu are the co-first authors.
This work was supported by the Major State Research Development Program of China [No. 2016YFC1202302], the National Natural Science Foundation of China (Grant No. 61502125, $2000), Heilongjiang Postdoctoral Fund (Grant No. LBH-Z15179, $800), and China Postdoctoral Science Foundation (Grant No. 2016 M590291, $1000).
Availability of data and materials
The system can be accessed from http://www.bio-annotation.cn:8080/DisSetSim/.
About this supplement
This article has been published as part of Journal of Biomedical Semantics Volume 8 Supplement 1, 2017: Selected articles from the Biological Ontologies and Knowledge bases workshop. The full contents of the supplement are available online at https://jbiomedsem.biomedcentral.com/articles/supplements/volume-8-supplement-1.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Hu, Y., Zhao, L., Liu, Z. et al. DisSetSim: an online system for calculating similarity between disease sets. J Biomed Semant 8 (Suppl 1), 28 (2017). https://doi.org/10.1186/s13326-017-0140-2
- Functional similarity
- Similarity score
- Disease sets
- Disease-miRNA relationships