- Research article
- Open Access
Revealing and avoiding bias in semantic similarity scores for protein pairs
BMC Bioinformatics volume 11, Article number: 290 (2010)
Semantic similarity scores for protein pairs are widely applied in functional genomic researches for finding functional clusters of proteins, predicting protein functions and protein-protein interactions, and for identifying putative disease genes. However, because some proteins, such as those related to diseases, tend to be studied more intensively, annotations are likely to be biased, which may affect applications based on semantic similarity measures. Thus, it is necessary to evaluate the effects of the bias on semantic similarity scores between proteins and then find a method to avoid them.
First, we evaluated 14 commonly used semantic similarity scores for protein pairs and demonstrated that they significantly correlated with the numbers of annotation terms for the proteins (also known as the protein annotation length). These results suggested that current applications of the semantic similarity scores between proteins might be unreliable. Then, to reduce this annotation bias effect, we proposed normalizing the semantic similarity scores between proteins using the power transformation of the scores. We provide evidence that this improves performance in some applications.
Current semantic similarity measures for protein pairs are highly dependent on protein annotation lengths, which are subject to biological research bias. This affects applications that are based on these semantic similarity scores, especially in clustering studies that rely on score magnitudes. The normalized scores proposed in this paper can reduce the effects of this bias to some extent.
Many scores for measuring semantic similarity (also termed functional similarity) between proteins have been proposed, based on the Gene Ontology (GO) terms  used to annotate the proteins. Some semantic similarity scores for a protein pair [2, 3] are calculated by combining the similarity scores for the term pairs [4–7] describing the two proteins. Other scores between proteins that do not use pairwise similarity scores between terms have also been proposed [8, 7, 9–12]. Similarity scores for protein pairs have been widely applied in functional genomic research . These scores are commonly used to analyze the correlation between functional similarity and similarities on other aspects, such as amino acid sequence similarity [2, 8, 14–16], or expression similarity [17–19]. Another type of applications is finding functional clusters of proteins [7, 20–22], or functional modules in physical or genetic protein-protein interaction networks [23–28]. Similarity scores are also used to predict protein functions [29–35], protein-protein interactions [36–41] and putative disease genes [42–45].
GO protein annotations are known to be incomplete , and suffer from a large research bias, because certain proteins, such as those related to diseases, tend to be studied more intensively [43, 47, 48]. Such an annotation bias may affect protein semantic similarity scores. In this paper, we evaluated 14 common semantic similarity scores for protein pairs, and demonstrated that the scores significantly correlated with the numbers of annotation terms for the proteins (i.e., the annotation length). Thus, we proposed normalizing the scores based on their power transformation to reduce annotation bias effects, and we provide evidence that this improves performance in some applications.
Gene Ontology (GO)
The GO annotation data for human was downloaded from the UniProt database http://www.ebi.ac.uk/GOA/index.html, including versions in November (Nov) 2001, Nov 2002, Nov 2003, Nov 2004, Nov 2005, Nov 2006, Nov 2007 and August 2008. The GO vocabulary data were downloaded from the GO website http://www.geneontology.org in August 2008. Here, we considered only IS-A links in GO [5, 6], and we mainly present the results based on the "Biological Process" (BP) sub-ontology. We also observed that all the semantic similarity scores for pairs of term groups are dependent on the annotation lengths when using "Molecular Function" and "Cellular Component" (data not shown).
Online Mendelian Inheritance in Man (OMIM) database and disease classification
The data for 1996 genes associated with 2192 diseases were downloaded from the OMIM database ftp://ftp.ncbi.nih.gov/repository/OMIM in August 2008, of which 1752 genes were annotated to GO BP terms. According to Goh et al. , the 2192 diseases were classified into 20 primary disorder classes based on the affected physiological systems. Diseases with multiple clinical features were assigned to the "multiple" class, and disease assigned to "Unclassified" class were not analyzed.
Similarity scores for term pairs
Many semantic similarity scores for two proteins are based on combinations of the similarity scores for term pairs between two groups of protein annotation terms. We evaluated four semantic similarity scores for term pairs based on the information contents: Resnik score , Lin score , Relevance score (RS)  and Jiang score . The information content of a term c was defined as IC(c) = -log(p(c)), where p(c) is the number of gene products annotated to the term and its descendants, divided by the number of all gene products annotated to the GO BP ontology. Let P(m, n) represent the set of common ancestor terms of m and n, then the four scores between terms m and n were calculated as:
Similarity scores for protein pairs based on pairwise similarity scores between term groups
In some methods, the similarity scores for term pairs describing two proteins are combined to calculate the semantic similarity scores of the two proteins. Here, two combination methods were evaluated: the arithmetic average (AVG) of the pairwise semantic similarity scores between two groups of GO terms describing the two proteins  and the best-match average (BMA) approach .
A 1 and A 2 were the groups of annotation terms for proteins P 1 and P 2 , and #P 1 and #P 2 were the number of terms included in A 1 and A 2 . Then the two combined scores between the two proteins were defined as:
where , .
In total, eight semantic similarity measures for protein pairs were evaluated, using the four semantic similarity scores for term pairs (Resnik, Lin, RS and Jiang) combined with the AVG and BMA methods (see Table 1).
Similarity scores for protein pairs based on groupwise similarity scores between term groups
We also evaluated six protein semantic similarity scores that do not use pairwise similarity scores between two term groups. These similarity scores are briefly described as below (please see details in the original papers).
The TO (Term Overlap) score  simply counts the number of overlapped terms for two proteins P 1 and P 2 as follows:
where GA1 and GA2 include the terms directly annotated with proteins P 1 and P 2 and all their ancestral terms, respectively. #(GA 1 ∩ GA 2 ) is the number of the overlapped terms between GA 1 and GA 2 .
The NTO (Normalized Term Overlap) score  is defined as:
The Dice score  is defined as:
The Kappa score  is defined as:
where and represent the observed and random co-occurrence of GO annotation terms for the two proteins, respectively.
The Graph Information Content (GIC) score  is defined as:
The Vector Space Model (VSM) score  is defined as follow:
where n is the number of all the GO BP terms and is the information content of term k if it is annotated with protein P 1 (P 2 ) while is 0 if the term k is not an annotation of the protein P 1 (P 2 ).
In total, we evaluated 14 semantic similarity scores for protein pairs (see Table 1). We note that some other semantic similarity scores for protein pairs [13, 51] were not evaluated in this paper. For example, the score proposed by Wang et al. , which weights the IS-A and PART-OF links of GO, was not analyzed, because we considered only IS-A links in this study.
Using randomly selected pairs of term groups, we evaluated the increase in protein semantic similarity score that resulted from only the increased annotation length, regardless of other biological factors. First, we randomly selected 10,000 pairs of term groups with the same sizes (corresponding to the annotation lengths of proteins) ranging from 1 to 10, since only 1.5% of proteins had more than 10 annotations in GO BP ontology. Then, using each of the 14 semantic similarity scores described above, we calculated the semantic similarity scores for random term group pairs, and analyzed whether these scores increased as the group size increased using the Spearman rank correlation coefficient .
Normalization based on power transformation
As demonstrated in the Results section, a similarity score for two groups of terms is dependent on the lengths of the term groups. To reduce the effect of the lengths on the scores, we took two steps to make the scores for pairs of term groups with given length combinations follow the standard normal distribution.
Firstly, we applied the commonly used power transformation approach to transform the scores to achieve normality [53, 54]. Suppose SS(TG1, TG2) is the score for term groups TG 1 with length L 1 and TG 2 with length L 2 , we power-transformed it to TSS(TG1, TG2) = . Here, the power was estimated as follow [53, 54]:
where M SS is the median of the random SS(TG1, TG2) distribution which was estimated by the similarity scores for 10,000 pairs of random term groups (with lengths L 1 and L 2 ). T q and T1-qare the lower and upper q th quantiles of this distribution (). By the one-sample Kolmogorov-Smirnov test for distribution goodness-of-fit , at the significance level of 0.1, we tested whether the power-transformed scores for pairs of term groups with given length combinations fit normal distributions.
Secondly, we normalized the above power-transformed scores to make them fit the standard normal distribution as follow:
In the above normalization formula, M TSS and STD TSS are the median and standard deviation of the power-transformed scores respectively. Here, we used the median rather than the mean in the normalization formula because it might be more appropriate for measuring the location parameter of a distribution when the distribution might be skewed [56–58]. As shown in the Results section, most of the normalized scores for pairs of term groups with given length combinations follow normal distributions. In this situation, the means and the medians are equal.
Sequence similarity scores for protein pairs
Amino acid sequence data for human proteins was downloaded from UniProt ftp://ftp.uniprot.org in August 2008. The sequence similarity between two proteins was measured by the ln(bit score), and calculated by the NCBI "blastall" program . Sequence similarity scores were obtained for a total of 499,878 protein pairs with GO BP annotations.
Clustering algorithm and enrichment analysis
Suppose the original and normalized similarity scores for two proteins (P 1 and P 2 ) annotated with two groups of terms are SS(P 1 , P 2 ) and NSS(P 1 , P 2 ) respectively, we firstly transformed both SS(P 1 , P 2 ) and NSS(P 1 , P 2 ) to the range [0, 1] by the Min-Max normalization method [59, 60] as follows
where Max SS and Min SS are the maximum and minimum values of the original similarity scores for all protein pairs from a protein set (e.g., a set of proteins encoded by a set of disease genes). Max NSS and Min NSS are the maximum and minimum values of the normalized similarity scores for all these protein pairs.
Then, we calculated the distance between two proteins as D(P 1 , P 2 ) = 1-MM(P 1 , P 2 ) based on the original score. Similarly, based on the normalized score, the distance was calculated as ND(P 1 , P 2 ) = 1-NM(P 1 , P 2 ). Both D(P 1 , P 2 ) and ND(P 1 , P 2 ) take values ranging from 0 to 1 and satisfy three main properties of distance metrics : (i) symmetry, D(P 1 , P 2 ) = D(P 2 , P 1 ) (ND(P 1 , P 2 ) = ND(P 2 , P 1 )); (ii) non-negative, D(P 1 , P 2 ) ≥ 0 (ND(P 1 , P 2 ) ≥ 0); (iii) triangle inequality, D(P 1 , P 3 ) ≤ D(P 1 , P 2 )+D(P 2 , P 3 ) (ND(P 1 , P 3 ) ≤ ND(P 1 , P 2 )+ND(P 2 , P 3 )). Using D(P 1 , P 2 ) and ND(P 1 , P 2 ) respectively, we clustered disease genes by the complete linkage clustering algorithm .
To evaluate the clustering results, using the hypergeometric distribution model [63, 64], we calculated the probability p of detecting at least the observed number of genes related to a disease category proposed by Goh et al.  in a cluster of disease genes by random chance. The p values were corrected by the false discovery rate (FDR) by the Benjamini-Hochberg (BH) procedure . With FDR of 1%, we found the disease categories enriched in a cluster of disease genes found by the clustering algorithm.
The dependence of the semantic similarity scores on annotation lengths
From 2001 to 2008, the average number of GO BP terms annotated with disease genes increased from 2.6 to 5.1, as shown in Figure 1(A). In contrast, the average annotation length of "non-disease" genes increased slightly from 1.7 to 2.1 (Figure 1(B)). These results indicated that disease genes tend to be studied more extensively, and are biased to have more annotations.
As shown in Figure 2, for each of the 14 protein semantic similarity scores analyzed, the median score for 10,000 random pairs of term groups increased significantly as the annotation lengths (the group sizes) increased (Spearman r ≥ 0.99, p < 1E-07). Based on the Resnik , Lin , RS  and Jiang  semantic similarity scores for term pairs, all four AVG combined scores for protein pairs were relatively stable when the annotation length was greater than three, especially for the Jiang(AVG). In contrast, the combined scores based on Resnik(BMA), Lin(BMA) and RS(BMA) increased rapidly with the annotation length. We evaluated six other semantic similarity scores for protein pairs, which do not use pairwise similarity scores between terms. As shown in Figure 2, the TO score was linearly dependent on the annotation length, while other scores increased slightly but significantly as the annotation length increased (Spearman r = 1, p < 1E-09). Notably, as shown in Figure 2, as the annotation lengths increased, the standard deviation of TO scores increased but it decreased for other similarity scores, which could be explained statistically. For example, because TO scores follow the hypergeometric probability model [63, 64], we can derive that its standard deviation increases with the annotation lengths.
Applications of the normalized scores
As shown in Table 2, based on each of the 14 similarity measures for term groups, most of the original scores (SS(TG1, TG2)) for pairs of term groups with given length combinations did not fit normal distributions (p ≥ 0.1, one-sample Kolmogorov-Smirnov test). For nine similarity measures, namely the Resnik(AVG), Resnik(BMA), Lin(AVG), Lin(BMA), RS(AVG), RS(BMA), Jiang(AVG), Dice and Kappa scores, over 80% of the power-transformed scores for pairs of term groups with given length combinations followed normal distributions. Then, these power-transformed scores were normalized to the standard normal distribution. Thus, for these nine similarity measures, the normalization method based on the power transformation is largely suitable for comparing scores for pairs of term groups with different length combinations. However, based on each of the other five similarity measures, less than 60% of the power-transformed scores fitted normal distributions. We also analyzed another five simple transformation methods and the results (see Table 2) showed that all these simple transformation methods performed worse than the power-transformation method using the estimated .
Then, for two types of applications, by comparing the results based on the original scores and their corresponding normalized scores, we showed that the bias affects certain analysis more than others. One type of applications based on semantic similarity scores for protein pairs study the correlation between functional similarity and similarities on other aspects [2, 8, 14–19]. Based on the normalized RS(BMA) score (the corresponding distribution for this measure was shown in Figure 3), we analyzed the correlation between protein semantic similarity scores and their amino acid sequence similarity scores (ln(bit score)). As shown in Figure 4, the correlation was significant (p < 2.20E-16), supporting the model that proteins with similar sequences tend to be functionally similar [2, 14]. Based on the RS(BMA) scores, similar results were observed, because of significant correlation between the ranks of the RS(BMA) scores and the normalized RS(BMA) scores for protein pairs (Spearman r = 0.88, p < 2.20E-16).
Another type of applications is clustering of functionally similar proteins [7, 20–22] or finding functional modules in physical or genetic protein-protein interaction networks [23–28]. Using the RS(BMA) and the normalized RS(BMA) distance, we applied a hierarchy clustering algorithm to cluster the disease genes into 21 categories, and compared the results with the categories determined by Goh et al. . As evaluated by the hypergenomic distribution test, using FDR of 1%, 16 clusters based on the normalized distance were enriched with disease genes with the same or similar phenotypes while only 6 clusters were enriched based on the original distance. To analyze more clearly the effect of annotation length on the cluster results, we clustered only the genes determined to the "Hematological" and "Immunological" categories. As shown in Figure 5(A), based on the normalized RS(BMA) distance, 73.5% of "Hematological" genes (red) were clustered into one group (p = 7.3E-13), while 78.4% of "Immunological" genes (blue) were in another (p = 3.96E-13). In contrast, as shown in Figure 5(B), when clustering these two classes of disease genes into two groups based on the RS(BMA) distance, no group was significantly enriched with a class of disease genes (p > 0.10). As shown in Figure 6, after normalization, the ranks of some disease gene pairs with different annotation lengths changed, improving the clustering results based on the normalized RS(BMA) distance. In general, based on the normalized RS(BMA) scores, our results suggested that disease genes related to the same or similar diseases tend to work together in the same disease-related functional gene modules [50, 66–78].
In this paper, we found that most semantic similarity scores for protein pairs increased as protein annotation lengths increased. Because protein annotations are likely to be subject to biological research bias, most applications based on current semantic similarity scores for protein pairs will be biased. Without the annotation bias, one could argue that over-annotated proteins might be more likely to be similar than under-annotated proteins, when considering only shared functions, and disregarding differences. However, currently, most semantic similarity scores for protein pairs evaluate the overall functional similarity between proteins. Depending on the available knowledge about domains, and the final aim of the application, different criteria could be used to define similarity between proteins. We note that protein annotations in GO for most model organisms (e.g., Saccharomyces cerevisiae) are also incomplete and suffer from the research bias because important genes such as the homologues of human disease genes tend to be studied more intensively [79, 80]. By analyzing the Saccharomyces cerevisiae data, we also found that the similarity scores between two groups of terms increased significantly with the annotation lengths (data not shown). Thus, our conclusion on the bias of semantic similarity scores for proteins would be applicable to other organisms.
A protein is usually annotated to a group of GO terms. Often, the semantic similarity scores between two proteins are calculated using some combination methods [2, 3] based on the semantic similarity scores for pairs of terms annotated with the two proteins. Many semantic similarity scores for term pairs such as the Resnik , Lin , Relevance (RS)  and Jiang  are based on the information content (related to the annotation specificity) of the terms. Based on these similarity scores for term pairs, the similarity scores for two proteins might not always increase, if the proteins have many annotations with low-specificity. However, as shown here, all the AVG and BMA scores for protein pairs based on the Resnik, Lin, RS and Jiang scores for term pairs still significantly correlated with the protein annotation lengths.
To reduce the effects of protein annotation bias, we normalized the scores based on the power transformation by estimating power . The normalization method based on the power transformation can transform most scores based on nine of the similarity measures to fit normal distributions but it performs poorly for the other five similarity measures. Thus, future works are needed to further improve the data transformation and normalization method.
The feasibility of the normalized scores was analyzed for two types of applications and the results showed that the normalized scores were useful in these applications. Analysis of the correlation between functional similarities and similarities on other aspects [2, 8, 14–19] might be less affected by the annotation bias, because the ranks of semantic similarity scores for protein pairs and their corresponding normalized scores were highly correlated. Our results also showed that clustering analysis [7, 20–22] using the magnitude of the semantic similarity scores might be more seriously affected by biased protein annotations, and the results could be improved by using the normalized scores.
A third type of applications that uses protein semantic similarity scores is predicting protein functions [29–35], protein-protein interactions [36–41] and disease genes [42–45]. However, because many other factors, such as the selection of algorithms and the definition of positive and negative sets  can affect the prediction results, we did not evaluate the effect of the annotation bias on these uses. Nevertheless, because this type of applications also uses the similarity score magnitudes, we recommend also using normalized scores in prediction studies, to reduce the effects of the annotation bias.
To avoid the influence of annotation bias, other approaches may be attempted. For example, the statistical p-value of a semantic similarity score for a protein pair could be evaluated by comparing this score with the scores of random protein pairs with the same annotation lengths. If the semantic similarity score of the two proteins was significantly larger than the score expected by random chance, at a given significance level (p-value), we could determine that the two proteins are functionally similar . Functional modules could be found by linking functionally related proteins. To analyze the functional relationships of proteins more comprehensively, the semantic similarity scores should be combined with other functional data, such as protein-protein interaction, co-expression and co-conservation of proteins [83–85].
Current protein semantic similarity scores are highly dependent on protein annotation lengths, which are subject to biological research bias. This bias may affect many current applications based on these scores. The proposed normalization method might solve this problem to some extend.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al.: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet 2000, 25(1):25–29. 10.1038/75556
Lord PW, Stevens RD, Brass A, Goble CA: Investigating semantic similarity measures across the Gene Ontology: the relationship between sequence and annotation. Bioinformatics 2003, 19(10):1275–1283. 10.1093/bioinformatics/btg153
Marino-Ramirez L, Bodenreider O, Kantz N, Jordan IK: Co-evolutionary Rates of Functionally Related Yeast Genes. Evol Bioinform Online 2006, 2: 295–300.
Schlicker A, Domingues FS, Rahnenfuhrer J, Lengauer T: A new measure for functional similarity of gene products based on Gene Ontology. BMC Bioinformatics 2006, 7: 302. 10.1186/1471-2105-7-302
Lin D: An information-theoretic definition of similarity. Proc 15th International Conf on Machine Learning: 1998 1998, 296–304.
Resnik P: Using information content to evaluate semantic similarity in a taxonomy. Proceedings of the 14th International Joint Conference on Artificial Intelligence: 1995 1995, 448–453.
Ovaska K, Laakso M, Hautaniemi S: Fast Gene Ontology based clustering for microarray experiments. BioData Min 2008, 1(1):11. 10.1186/1756-0381-1-11
Pesquita C, Faria D, Bastos H, Ferreira AE, Falcao AO, Couto FM: Metrics for GO based protein semantic similarity: a systematic evaluation. BMC Bioinformatics 2008, 9(Suppl 5):S4. 10.1186/1471-2105-9-S5-S4
Mistry M, Pavlidis P: Gene Ontology term overlap as a measure of gene functional similarity. BMC Bioinformatics 2008, 9: 327. 10.1186/1471-2105-9-327
Chabalier J, Mosser J, Burgun A: A transversal approach to predict gene product networks from ontology-based similarity. BMC Bioinformatics 2007, 8: 235. 10.1186/1471-2105-8-235
Huang da W, Sherman BT, Tan Q, Collins JR, Alvord WG, Roayaei J, Stephens R, Baseler MW, Lane HC, Lempicki RA: The DAVID Gene Functional Classification Tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol 2007, 8(9):R183. 10.1186/gb-2007-8-9-r183
Martin D, Brun C, Remy E, Mouren P, Thieffry D, Jacq B: GOToolBox: functional analysis of gene datasets based on Gene Ontology. Genome Biol 2004, 5(12):R101. 10.1186/gb-2004-5-12-r101
Pesquita C, Faria D, Falcao AO, Lord P, Couto FM: Semantic similarity in biomedical ontologies. PLoS Comput Biol 2009, 5(7):e1000443. 10.1371/journal.pcbi.1000443
Joshi T, Xu D: Quantitative assessment of relationship between sequence similarity and function similarity. BMC Genomics 2007, 8: 222. 10.1186/1471-2164-8-222
Yang L, Yu J: A comparative analysis of divergently-paired genes (DPGs) among Drosophila and vertebrate genomes. BMC Evol Biol 2009, 9: 55. 10.1186/1471-2148-9-55
Altenhoff AM, Dessimoz C: Phylogenetic and functional assessment of orthologs inference projects and methods. PLoS Comput Biol 2009, 5(1):e1000262. 10.1371/journal.pcbi.1000262
Elo LL, Jarvenpaa H, Oresic M, Lahesmaa R, Aittokallio T: Systematic construction of gene coexpression networks with applications to human T helper cell differentiation process. Bioinformatics 2007, 23(16):2096–2103. 10.1093/bioinformatics/btm309
Lee HK, Hsu AK, Sajdak J, Qin J, Pavlidis P: Coexpression analysis of human genes across many microarray data sets. Genome Res 2004, 14(6):1085–1094. 10.1101/gr.1910904
Wang H, Azuaje F, Bodenreider O, Dopazo J: Gene expression correlation and gene ontology-based similarity: an assessment of quantitative relationships. Proceedings of the 2004 IEEE Symposium on Computational Intelligence in Bioinformatics and Computational Biology, CIBCB'04: 2004 2004, 25–31.
Chen JL, Liu Y, Sam LT, Li J, Lussier YA: Evaluation of high-throughput functional categorization of human disease genes. BMC Bioinformatics 2007, 8(Suppl 3):S7. 10.1186/1471-2105-8-S3-S7
Du Z, Li L, Chen CF, Yu PS, Wang JZ: G-SESAME: web tools for GO-term-based gene similarity analysis and knowledge discovery. Nucleic Acids Res 2009, (37 Web Server):W345–349. 10.1093/nar/gkp463
Wang JZ, Du Z, Payattakool R, Yu PS, Chen CF: A new method to measure the semantic similarity of GO terms. Bioinformatics 2007, 23(10):1274–1281. 10.1093/bioinformatics/btm087
Ulitsky I, Shlomi T, Kupiec M, Shamir R: From E-MAPs to module maps: dissecting quantitative genetic interactions using physical interactions. Mol Syst Biol 2008, 4: 209. 10.1038/msb.2008.42
Bandyopadhyay S, Kelley R, Krogan NJ, Ideker T: Functional maps of protein complexes from quantitative genetic interaction data. PLoS Comput Biol 2008, 4(4):e1000065. 10.1371/journal.pcbi.1000065
Sen TZ, Kloczkowski A, Jernigan RL: Functional clustering of yeast proteins from the protein-protein interaction network. BMC Bioinformatics 2006, 7: 355. 10.1186/1471-2105-7-355
Lubovac Z, Gamalielsson J, Olsson B: Combining functional and topological properties to identify core modules in protein interaction networks. Proteins 2006, 64(4):948–959. 10.1002/prot.21071
Krogan NJ, Cagney G, Yu H, Zhong G, Guo X, Ignatchenko A, Li J, Pu S, Datta N, Tikuisis AP, et al.: Global landscape of protein complexes in the yeast Saccharomyces cerevisiae. Nature 2006, 440(7084):637–643. 10.1038/nature04670
Taylor IW, Linding R, Warde-Farley D, Liu Y, Pesquita C, Faria D, Bull S, Pawson T, Morris Q, Wrana JL: Dynamic modularity in protein interaction networks predicts breast cancer outcome. Nat Biotechnol 2009, 27(2):199–204. 10.1038/nbt.1522
Chen XW, Liu M, Ward R: Protein function assignment through mining cross-species protein-protein interactions. PLoS ONE 2008, 3(2):e1562. 10.1371/journal.pone.0001562
Zhu M, Gao L, Guo Z, Li Y, Wang D, Wang J, Wang C: Globally predicting protein functions based on co-expressed protein-protein interaction networks and ontology taxonomy similarities. Gene 2007, 391(1–2):113–119. 10.1016/j.gene.2006.12.008
Tao Y, Sam L, Li J, Friedman C, Lussier YA: Information theory applied to the sparse gene ontology annotation network to predict novel gene function. Bioinformatics 2007, 23(13):i529–538. 10.1093/bioinformatics/btm195
Tu K, Yu H, Guo Z, Li X: Learnability-based further prediction of gene functions in Gene Ontology. Genomics 2004, 84(6):922–928. 10.1016/j.ygeno.2004.08.005
Cakmak A, Ozsoyoglu G: Discovering gene annotations in biomedical text databases. BMC Bioinformatics 2008, 9: 143. 10.1186/1471-2105-9-143
Cho YR, Shi L, Ramanathan M, Zhang A: A probabilistic framework to predict protein function from interaction data integrated with semantic knowledge. BMC Bioinformatics 2008, 9: 382. 10.1186/1471-2105-9-382
Fontana P, Cestaro A, Velasco R, Formentin E, Toppo S: Rapid annotation of anonymous sequences from genome projects using semantic similarities and a weighting scheme in gene ontology. PLoS ONE 2009, 4(2):e4619. 10.1371/journal.pone.0004619
Futschik ME, Chaurasia G, Herzel H: Comparison of human protein-protein interaction maps. Bioinformatics 2007, 23(5):605–611. 10.1093/bioinformatics/btl683
Wu X, Zhu L, Guo J, Zhang DY, Lin K: Prediction of yeast protein-protein interaction network: insights from the Gene Ontology and annotations. Nucleic Acids Res 2006, 34(7):2137–2150. 10.1093/nar/gkl219
Ofran Y, Yachdav G, Mozes E, Soong TT, Nair R, Rost B: Create and assess protein networks through molecular characteristics of individual proteins. Bioinformatics 2006, 22(14):e402–407. 10.1093/bioinformatics/btl258
Tarassov K, Messier V, Landry CR, Radinovic S, Serna Molina MM, Shames I, Malitskaya Y, Vogel J, Bussey H, Michnick SW: An in vivo map of the yeast protein interactome. Science 2008, 320(5882):1465–1470. 10.1126/science.1153878
Xu T, Du L, Zhou Y: Evaluation of GO-based functional similarity measures using S. cerevisiae protein interaction and expression profile data. BMC Bioinformatics 2008, 9: 472. 10.1186/1471-2105-9-472
Soong TT, Wrzeszczynski KO, Rost B: Physical protein-protein interactions predicted from microarrays. Bioinformatics 2008, 24(22):2608–2614. 10.1093/bioinformatics/btn498
Gaulton KJ, Mohlke KL, Vision TJ: A computational system to select candidate genes for complex human traits. Bioinformatics 2007, 23(9):1132–1140. 10.1093/bioinformatics/btm001
Adie EA, Adams RR, Evans KL, Porteous DJ, Pickard BS: SUSPECTS: enabling fast and effective prioritization of positional candidates. Bioinformatics 2006, 22(6):773–774. 10.1093/bioinformatics/btk031
Chen J, Bardes EE, Aronow BJ, Jegga AG: ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res 2009, (37 Web Server):W305–311. 10.1093/nar/gkp427
Freudenberg J, Propping P: A similarity-based method for genome-wide prediction of disease-relevant human genes. Bioinformatics 2002, 18(Suppl 2):S110–115.
Rhee SY, Wood V, Dolinski K, Draghici S: Use and misuse of the gene ontology annotations. Nat Rev Genet 2008, 9(7):509–515. 10.1038/nrg2363
Verver O, Ridder Jd, Reinders MJT, Wessels LFA: Prioritization of Candidate Disease Genes using Microarray Data and Functional Relations. Bioinformatics 2007, 00: 1–12.
Adie EA, Adams RR, Evans KL, Porteous DJ, Pickard BS: Speeding disease gene discovery by sequence based candidate prioritization. BMC Bioinformatics 2005, 6: 55. 10.1186/1471-2105-6-55
Hamosh A, Scott AF, Amberger JS, Bocchini CA, McKusick VA: Online Mendelian Inheritance in Man (OMIM), a knowledgebase of human genes and genetic disorders. Nucleic Acids Res 2005, (33 Database):D514–517.
Goh KI, Cusick ME, Valle D, Childs B, Vidal M, Barabasi AL: The human disease network. Proc Natl Acad Sci USA 2007, 104(21):8685–8690. 10.1073/pnas.0701361104
Chagoyen M, Carazo J, Pascual-Montano A: Pairwise similarity scores using functional annotations: review and comparison. 8th Spanish Symposium on Bioinformatics and Computational Biology: 2008 2008.
Spearman C: The Proof and Measurement of Association Between Two Things. American Journal of Psychology 1904, 15: 72–101. 10.2307/1412159
Tan W, Gan F, Chang T: Using normal quantile plot to select an appropriate transformation to achieve normality. Computational Statistics & Data Analysis 2004, 45(3):609–619. 10.1016/S0167-9473(03)00009-4
Emerson J, Stoto M: Exploratory Methods for Choosing Power Transformations. Journal of the American Statistical Association 1982, 77: 103–108. 10.2307/2287775
Massey J, Frank J: The Kolmogorov-Smirnov test for goodness of fit. Journal of the American Statistical Association 1951, 46: 68–78. 10.2307/2280095
Kuczmarski RJ, Ogden CL, Grummer-Strawn LM, Flegal KM, Guo SS, Wei R, Mei Z, Curtin LR, Roche AF, Johnson CL: CDC growth charts: United States. Adv Data 2000, 314: 1–27.
Kuk A, Mak T: Median estimation in the presence of auxiliary information. Journal of the Royal Statistical Society Series B Methodological 1989, 51: 261–269.
Waterlow JC, Buzina R, Keller W, Lane JM, Nichaman MZ, Tanner JM: The presentation and use of height and weight data for comparing the nutritional status of groups of children under the age of 10 years. Bull World Health Organ 1977, 55(4):489–498.
Indovina M, Uludag U, Snelick R, Mink A, Jain A: Multimodal biometric authentication methods: a COTS approach. Proceedings of the MMUA: 2003 2003, 99–106.
Sorace JM, Zhan M: A data review and re-assessment of ovarian cancer serum proteomic profiling. BMC Bioinformatics 2003, 4: 24. 10.1186/1471-2105-4-24
Wang J, Wang X, Lin K, Shasha D, Shapiro B, Zhang K: Evaluating A Class of Distance-Mapping Algorithms for Data Mining and Clustering. Proceedings of the fifth ACM SIGKDD international conference on Knowledge discovery and data mining 1999 1999, 307–311. full_text
Yamada T, Kanehisa M, Goto S: Extraction of phylogenetic network modules from the metabolic network. BMC Bioinformatics 2006, 7: 130. 10.1186/1471-2105-7-130
Fury W, Batliwalla F, Gregersen PK, Li W: Overlapping probabilities of top ranking gene lists, hypergeometric distribution, and stringency of gene selection criterion. Conf Proc IEEE Eng Med Biol Soc 2006, 1: 5531–5534. full_text
Gonin H: The use of factorial moments in the treatment of the hypergeometric distribution and in tests for regression. Philosophical Magazine Series 7 1936, 21(139):215–226.
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B Methodological 1995, 57: 289–330.
Oti M, Brunner HG: The modular nature of genetic diseases. Clin Genet 2007, 71(1):1–11. 10.1111/j.1399-0004.2006.00708.x
Brunner HG, van Driel MA: From syndrome families to functional genomics. Nat Rev Genet 2004, 5(7):545–551. 10.1038/nrg1383
Oti M, Huynen MA, Brunner HG: Phenome connections. Trends Genet 2008, 24(3):103–106. 10.1016/j.tig.2007.12.005
Lage K, Karlberg EO, Storling ZM, Olason PI, Pedersen AG, Rigina O, Hinsby AM, Tumer Z, Pociot F, Tommerup N, et al.: A human phenome-interactome network of protein complexes implicated in genetic disorders. Nat Biotechnol 2007, 25(3):309–316. 10.1038/nbt1295
Lage K, Hansen NT, Karlberg EO, Eklund AC, Roque FS, Donahoe PK, Szallasi Z, Jensen TS, Brunak S: A large-scale analysis of tissue-specific pathology and gene expression of human disease genes and complexes. Proc Natl Acad Sci USA 2008, 105(52):20870–20875. 10.1073/pnas.0810772105
Girirajan S, Truong HT, Blanchard CL, Elsea SH: A functional network module for Smith-Magenis syndrome. Clin Genet 2009, 75(4):364–374. 10.1111/j.1399-0004.2008.01135.x
Moran LB, Graeber MB: Towards a pathway definition of Parkinson's disease: a complex disorder with links to cancer, diabetes and inflammation. Neurogenetics 2008, 9(1):1–13. 10.1007/s10048-007-0116-y
Li Y, Agarwal P: A pathway-based view of human diseases and disease relationships. PLoS One 2009, 4(2):e4346. 10.1371/journal.pone.0004346
Lim J, Hao T, Shaw C, Patel AJ, Szabo G, Rual JF, Fisk CJ, Li N, Smolyar A, Hill DE, et al.: 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.032
Goehler H, Lalowski M, Stelzl U, Waelter S, Stroedicke M, Worm U, Droege A, Lindenberg KS, Knoblich M, Haenig C, et al.: A protein interaction network links GIT1, an enhancer of huntingtin aggregation, to Huntington's disease. Mol Cell 2004, 15(6):853–865. 10.1016/j.molcel.2004.09.016
Bergholdt R, Storling ZM, Lage K, Karlberg EO, Olason PI, Aalund M, Nerup J, Brunak S, Workman CT, Pociot F: Integrative analysis for finding genes and networks involved in diabetes and other complex diseases. Genome Biol 2007, 8(11):R253. 10.1186/gb-2007-8-11-r253
Ergun A, Lawrence CA, Kohanski MA, Brennan TA, Collins JJ: A network biology approach to prostate cancer. Mol Syst Biol 2007, 3: 82. 10.1038/msb4100125
Jiang X, Liu B, Jiang J, Zhao H, Fan M, Zhang J, Fan Z, Jiang T: Modularity in the genetic disease-phenotype network. FEBS Lett 2008, 582(17):2549–2554. 10.1016/j.febslet.2008.06.023
Kaletta T, Hengartner MO: Finding function in novel targets: C. elegans as a model organism. Nat Rev Drug Discov 2006, 5(5):387–398. 10.1038/nrd2031
Langenau DM, Jette C, Berghmans S, Palomero T, Kanki JP, Kutok JL, Look AT: Suppression of apoptosis by bcl-2 overexpression in lymphoid cells of transgenic zebrafish. Blood 2005, 105(8):3278–3285. 10.1182/blood-2004-08-3073
Yu S, Van Vooren S, Tranchevent LC, De Moor B, Moreau Y: Comparison of vocabularies, representations and ranking algorithms for gene prioritization by text mining. Bioinformatics 2008, 24(16):i119–125. 10.1093/bioinformatics/btn291
Yang D, Li Y, Xiao H, Liu Q, Zhang M, Zhu J, Ma W, Yao C, Wang J, Wang D, et al.: Gaining confidence in biological interpretation of the microarray data: the functional consistence of the significant GO categories. Bioinformatics 2008, 24(2):265–271. 10.1093/bioinformatics/btm558
Lee I, Date SV, Adai AT, Marcotte EM: A probabilistic functional network of yeast genes. Science 2004, 306(5701):1555–1558. 10.1126/science.1099511
Franke L, van Bakel H, Fokkens L, de Jong ED, Egmont-Petersen M, Wijmenga C: Reconstruction of a functional human gene network, with an application for prioritizing positional candidate genes. Am J Hum Genet 2006, 78(6):1011–1025. 10.1086/504300
Segal E, Friedman N, Koller D, Regev A: A module map showing conditional activity of expression modules in cancer. Nat Genet 2004, 36(10):1090–1098. 10.1038/ng1434
This work was supported by the National Natural Science Foundation of China (grant nos. 30571034, 30970668, 30670539).
Conceived and designed the experiments: ZG, JW. Performed the experiments: JW, XZ and CZ. Analyzed the data: ZG, JW and JZ. Wrote the paper: JW, JZ and ZG. All authors read and approved the final manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Wang, J., Zhou, X., Zhu, J. et al. Revealing and avoiding bias in semantic similarity scores for protein pairs. BMC Bioinformatics 11, 290 (2010). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2105-11-290
- Gene Ontology
- Similarity Score
- Protein Pair
- Relevance Score
- Term Group