Email updates

Keep up to date with the latest news and content from BMC Systems Biology and BioMed Central.

Open Access Research article

Mining breast cancer genes with a network based noise-tolerant approach

Yaling Nie and Jingkai Yu*

Author Affiliations

National Key Laboratory of Biochemical Engineering, Institute of Process Engineering, Chinese Academy of Sciences, Beijing 100190, China

For all author emails, please log on.

BMC Systems Biology 2013, 7:49  doi:10.1186/1752-0509-7-49

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1752-0509/7/49


Received:28 November 2012
Accepted:21 June 2013
Published:25 June 2013

© 2013 Nie and Yu; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Mining novel breast cancer genes is an important task in breast cancer research. Many approaches prioritize candidate genes based on their similarity to known cancer genes, usually by integrating multiple data sources. However, different types of data often contain varying degrees of noise. For effective data integration, it’s important to design methods that work robustly with respect to noise.

Results

Gene Ontology (GO) annotations were often utilized in cancer gene mining works. However, the vast majority of GO annotations were computationally derived, thus not completely accurate. A set of genes annotated with breast cancer enriched GO terms was adopted here as a set of source data with realistic noise. A novel noise tolerant approach was proposed to rank candidate breast cancer genes using noisy source data within the framework of a comprehensive human Protein-Protein Interaction (PPI) network. Performance of the proposed method was quantitatively evaluated by comparing it with the more established random walk approach. Results showed that the proposed method exhibited better performance in ranking known breast cancer genes and higher robustness against data noise than the random walk approach. When noise started to increase, the proposed method was able to maintained relatively stable performance, while the random walk approach showed drastic performance decline; when noise increased to a large extent, the proposed method was still able to achieve better performance than random walk did.

Conclusions

A novel noise tolerant method was proposed to mine breast cancer genes. Compared to the well established random walk approach, it showed better performance in correctly ranking cancer genes and worked robustly with respect to noise within source data. To the best of our knowledge, it’s the first such effort to quantitatively analyze noise tolerance between different breast cancer gene mining methods. The sorted gene list can be valuable for breast cancer research. The proposed quantitative noise analysis method may also prove useful for other data integration efforts. It is hoped that the current work can lead to more discussions about influence of data noise on different computational methods for mining disease genes.

Keywords:
Network; Breast cancer; Data noise; Noise tolerance

Background

Novel disease genes remain difficult to identify in most genetic diseases, and in particular, in highly polygenic disorders. Currently, not all genes have yet been detected even for those diseases whose molecular mechanisms are partially known [1], for instance, breast cancer [2]. Breast cancer is a common cancer and a major cause of cancer death among females around the world, which makes up 23% of total cancer cases and 14% of cancer deaths [3]. Mining breast cancer genes is conducive to understand its pathogenic mechanism and search for effective treatments. With rapid growth of disease-related genomic and functional data, computational approaches can be utilized to mine for new cancer genes [4].

In the past two decades, a number of computational methods had been developed to mine potential disease related genes. Most of those methods rank candidate genes based on the idea that proteins similar to each other tend to cause similar or same diseases [5]. They involve setting up a candidate gene set to be compared with a known disease gene set on their physical or functional attributes [6]. On one hand, physical attribute-based methods include screening direct neighbors of known disease genes in the PPI network [7,8], comparing shortest path length [9] between candidate genes and known disease genes, clustering or graph partitioning to uncover disease modules in the interaction network [10-12]. Some approaches also used global network features to find genes similar with known disease genes [13,14]. On the other hand, several methods rely on functional similarities between candidate and disease genes [15], for example, some methods measured similarity between genes by their functional annotations [16] (e.g., Gene Ontology (GO) [17]). Methods using other data sources had also been developed, such as gene expression [18,19], biological pathways and sequence features [20].

Cancers such as breast cancer are complex and heterogeneous in nature, cancer-related genes often do not function in isolation but interact with one another [5]. Integrating multiple data types was found to be effective for gene mining in alleviating problems caused by incomplete information [21-23]. For instance, ENDEAVOUR [24] is an online tool based on using multiple data sources. It integrated candidate gene rankings from different data sources into a final ranking with the order statistic algorithm. However, different data categories usually contain inherent noise or systematic errors [25]. For instance, data from computational predictions will no doubt contain some amount of uncertainty. Experimental data obtained from different labs or experimental platforms can contain appreciable amount of noise. Noise in source data can push computed results away from their true values, lead to erroneous reporting.

A better method must be able to tolerate certain amount of noise, which makes the integration of different data sources more applicable to real-life scenarios. Despite the fact that some approaches can work with precision when presented with highly accurate data, few studies have shown that those methods worked robustly when faced with increasingly noisy data. A number of papers had discussed the task of balancing noise and precision when using multiple data sources for cancer gene mining, however, hardly anyone had analyzed the noise problem quantitatively [26-29]. It is important to calibrate how robust a method works with respect to noise, namely, how fast a method deteriorates when percentage of noise in source data goes up. With that knowledge, users can then be confident about the method’s effectiveness when it is applied to real life data sets.

To tackle the data noise problem, a novel noise tolerant data fusion approach was proposed here for breast cancer gene mining (Figure 1), which integrated information from PPI network with gene expression data to rank genes based on their probabilities of being breast cancer related. Satisfactory results were obtained even when noise level was high. To demonstrate advantages of the proposed method, its performance was compared with that of the random walk method [13], which utilized a node’s global neighborhood in a network to rank genes. Random walk based methods had been shown to produce good performance [18,23,29-31] in gene ranking. Results showed the proposed method exhibited better robustness when faced with increasingly noisy data, as compared to the random walk approach.

thumbnailFigure 1. Schematic chart for mining breast cancer genes. Four different types of data were used as input: PPI data, human gene expression data, known breast cancer genes and GO annotations. Gene expression data (GDSes) from the GEO database were clustered. Known breast cancer genes and their enriched GO annotations were used to rank genes in those clusters. From the PPI network, three network topological attributes were computed to rank genes in the network. Finally, all individual rankings were combined into a final rank, which represents a gene’s overall probability of being involved in breast cancer.

Results and discussion

An effective data integration method was developed to mine breast cancer genes from four major data sources: Protein-Protein interactions, gene expression data, GO annotations, and known breast cancer genes (Table 1).

Table 1. Data source

After removing redundancy, a comprehensive human PPI network was constructed with data obtained from multiple interactions databases. The resultant network contained a total of 156,459 PPIs with 15,494 genes. A noise tolerant method was designed to rank potential breast cancer genes.

Rationale for data integration

A network-based score (SN) and an expression-based score (SE) were respectively derived for each gene, which were then integrated into a final score (S) by weighting them with a coefficient λ. A P-score was computed to represent performance of the proposed method when λ was changed from 0 to 1. P-score was the average ranking of known breast cancer genes in top 10% of the final gene ranking list (see Methods). Smaller P-score (ranked higher) meant better performance (Figure 2). As shown in Figure 2, better performance was achieved when λ=0 (i.e., using only network-based ranking) compared with λ=1 (i.e., using only expression-based ranking); the best performance was obtained when λ=0.2, which suggested that utilizing genes’ complex relations in the PPI network can help cancer gene mining tasks. It also confirmed that multiple data sources can complement each other in ranking cancer genes. Final ranking results were listed in Additional file 1: Table S1.

thumbnailFigure 2. Performance of our method for different λ values. P-score was the average ranking of KnownSet in top 10% of the final sorted list. A smaller P-score (ranked higher) meant better capability to correctly rank known breast cancer genes.

Additional file 1: Table S1. Supplementary Table 1.

Format: XLSX Size: 369KB Download fileOpen Data

Evaluation of performance in ranking known cancer genes

Known cancer genes were derived from the OMIM (Online Mendelian Inheritance in Man) and CGC (Cancer Gene Census) databases (see Methods). Known breast cancer genes derived from OMIM was called OMIMSet, those from CGC called CGCSet. OMIMSet was used to train the proposed method, and CGCSet was used for evaluation. Table 2 showed that the proposed method achieved better performance in ranking known breast cancer genes. Counting only genes ranked in top 10%, the proposed method achieved an average ranking of 279, compared with 545 by random walk. Counting all 11 genes in CGCSet, the proposed method achieved an average ranking of 1801, compared with 2207 by random walk.

Table 2. Ranking performance comparison

Robustness with respect to realistic data noise

In general, a method based on multiple types of knowledge is more objective than those utilizing singular information; combining independent data sources can alleviate effects of biases inherent in single data types [32]. Most data sets, especially genome wide data, tend to contain appreciable amount of noise. For instance, GO is a powerful tool which provides a controlled vocabulary to describe biological functions on multiple levels [33]. It was also widely used on cancer research (there are about 39,000 citations for GO in Google Scholar by Oct. 16, 2012). However, not all GO annotations are equally credible [34]. As of October 2012, there were over 3 million GO annotations for Homo sapiens genes. Each GO annotation includes an evidence code to indicate how the annotation was inferred. All evidence codes (except IEA, which is automatically derived) are manually assigned. They can be divided into four categories: experimental (EXP, IDA, IPI, IMP, IGI, and IEP), computational analysis (ISS, ISO, ISA, ISM, IGC, IBA, IBD, IKR, IRD, and RCA), author statements (TAS and NAS), and curatorial statements (IC and ND). For Homo sapiens, 16 evidence codes were used to describe GO annotations [35]. The vast majority of GO annotations were found to be computationally derived and not manually curated (IEA evidence code) (Figure 3). GO annotations therefore contain appreciable amount of noise within themselves.

thumbnailFigure 3. Evidence code distribution of GO annotation (for Homo sapiens). IEA means that the denoted GO annotations were inferred from electronic annotation and not curated manually. IEA counted for the majority of GO annotations.

Approaches that exhibit robust performance with regard to noise are needed if they are to prove useful in cancer gene hunting endeavours. Nevertheless, as mentioned before, few projects had specifically analysed data noise effects quantitatively. A network based noise tolerant method was proposed here to mine breast cancer genes. Its performance was compared with that of the well performing random walk approach by five-fold cross-validation. The results confirmed the proposed method’s robust performance with respect to data noise.

The set of known breast cancer genes (KnownSet, see Methods) was enlarged by including genes sharing GO annotations with those in the KnownSet. The enlarged set was called the GOSet (GO enriched gene set, see Methods), which was adopted as a noisy set of likely breast cancer genes. The GOSet was utilized to check an algorithm’s robustness with respect to data noise. Data were sampled from the GOSet, and combined with the KnowSet to generate a noisy set of training data. This way of synthesizing noisy data set is unique in that it doesn’t simply using random data as noise, which is too artificial. The GOSet contains enriched but still imperfect data, which can better mimic data noise in real life scenarios. An algorithm’s ability to retain its performance was checked when fraction of noisy data in the training set went up.

The proposed method did not work quite as well as random walk approach when input data was 100% accurate, however, when noise level in input data increased, we observed the following phenomena (Figure 4).

1. Random walk approach exhibited a sharp decrease in its performance, while our method was able to maintain a relatively stable performance.

2. When noise increased to a large extent, the proposed method was able to perform about twice as well as random walk approach did.

thumbnailFigure 4. Five-fold cross-validation performance evaluation. In cross validation, different ratios of noise data were added to input data, with the ratio changing from 0 to 50, where 0 meant only known breast cancer genes were used as input data. Performance was evaluated in terms of the F-score. k was a ranking threshold to judge a ranked gene as a true breast cancer gene by the proposed method. (a)k=300; (b)k=400.

It can thus be stated that the proposed method was more robust with respect to noise in input data, compared to the state-of-art random walk based approach. The results also confirmed the power of data integration, which was able to let different data sets complementing each other [22,23].

Robustness with respect to completely random noise

GOSet tried to simulate realistic data noise. However, it might be suspected that GOSet was biased toward the proposed method in one way or another. To make sure the comparison between the proposed method and the random walk approach was not done unfairly, randomly picked genes were added to the KnownSet (Figure 1, Box B, see Methods), and performances of the proposed method and random walk approach were then compared. Figure 5 again showed that random walk approach showed linear performance decline when noise increased in the source data, while the proposed method was able to maintain stable performance.

thumbnailFigure 5. Ranking comparison performance. Randomly picked genes were added to the KnownSet, and performances of the proposed method and random walk approach were then compared.

Conclusions

Cancers are highly complex processes, the majority of cancer genes are yet to be mapped. Currently available data (known breast cancer genes) are too limited to be really effective for cancer gene searching purpose. Broadening the scope of input data (both volume and type) should enable better use of available data to mine for new cancer genes. Approaches that work robustly against data noise are needed.

A novel noise tolerant breast cancer gene mining method was presented here, which integrated a comprehensive PPI network, gene expression data, prior knowledge of breast cancer and GO annotations to rank potential breast cancer genes. From each data source, a ranked list for each candidate gene was computed, and they were then combined into a final ranking order. Influence of data noise was quantitatively evaluated. Random walk approach performed better than the proposed method using 100% accurate input data (known breast cancer genes). However, the proposed method showed much greater noise tolerance. To our best knowledge, this is the first effort to quantitatively analyse noise tolerance between different cancer gene mining methods. The framework of the proposed mining method and the quantitative way of appraising noise effects are flexible enough to be useful for other data sources, and hopefully, lead to more discussions on data noise issue for different computational methods in cancer gene mining field.

Methods

Figure 1 presented a schematic view of our approach. A comprehensive PPI network was obtained by integrating data from different interactions databases [36] (Box A). A set of known breast cancer genes (KnownSet) was extracted from the OMIM and CGC databases (Box B). Candidate genes were first ranked by three network topological attributes: node degree, node betweenness and by their closeness to known cancer genes in the network (Box C). GO term enrichment analyses were performed for KnownSet, producing a GO term set enriched with breast cancer related terms, into which a group of cancer-hallmark GO terms were also added [37]. A set of genes which were annotated with terms in the obtained GO set were generated, which was called the GO enriched gene set (GOSet) (Box D). A batch of breast cancer-related expression data was extracted from the GEO database [38] on April 7, 2011 and expression profiles in those data files were clustered based on their similarity with each other (Box E). Expression clusters were intersected with GOSet. Overlap significance was represented by a p-value computed with the normal distribution. The p-value was utilized to rank genes in expression clusters (Box F). All individual rankings from different data sources were finally combined into a final ranking, which represented a gene’s overall probability of being involved in breast cancer (Box G).

Deriving known breast cancer gene set

Thirty known breast cancer genes were extracted from the OMIM database [39] and 19 from the CGC database [40] (Table 3). For a gene to be usable, it was required to be covered by both the PPI network and expression data sets. With that requirement, 26 genes derived from OMIM (OMIMSet) and 11 additional genes from CGC (CGCSet) were obtained, this set of 37 genes was called the KnownSet.

Table 3. The known breast cancer genes

Ranking by PPI network

The human PPI data were derived from five sources: HPRD [41], BioGRID [42], homoMINT [43], IntAct [44] and a manually curated human signalling network [45]. Protein identifiers were mapped to uniform coding gene identifiers. Official gene symbols were used as identifier. Redundant interactions were removed, along with interactions with identifiers that could not be mapped to gene symbols (Table 1). The final PPI network was represented by an undirected graph where nodes representing genes and edges representing interactions. The graph contained 156,459 interactions connecting 15,494 genes.

Similarities between proteins were found to be correlated with their proximity in the PPI network [46]. It was assumed that when a gene in the PPI network exhibited topological features similar to known breast cancer genes, it’s more likely to be involved in breast cancer processes. Several papers had shown that cancer genes could be effectively distinguished from others by their topological attributes in the PPI network, such as node degree [47], betweenness centrality [48] and shortest path length [10]. The above three network topological indices were computed and used to assess gene similarity in the PPI network. Genes were then sorted according to values of the topological indices.

Let G(V,E) be the PPI network, where V is the set of genes, and E the set of interactions in the network.

For a node vV, degreecd(v) is the number of direct neighbours of v in the network. Betweenness centrality is the sum of the fraction of all-pair shortest paths that pass through v[49]. It can be expressed as following,

<a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M1">View MathML</a>

(1)

where <a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M2">View MathML</a>jk is the number of shortest paths from a source jV to a target kV, and <a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M3">View MathML</a>jk(v) is the number of those paths passing through some node v other than j,k. If j=k, <a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M4">View MathML</a>jk =1, and if vj,k, <a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M5">View MathML</a>jk(v)=0.

The shortest path length was defined as the average shortest path distance from all known breast cancer genes to node v, it can be denoted as following,

<a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M6">View MathML</a>

(2)

where G is the KnownSet, d(v,t) is the shortest path length between node v and t. n is the number of known breast cancer genes which can be reached by v.

The above topological attributes were computed with the Python package networkx [50]. To facilitate later integration step, they were transformed into rankings; that is, each list was sorted, and a gene was assigned a positive integer number for a specific attribute according to that attribute’s value in the sorted list (Figure 6). In the end, three rankings for each gene were obtained based on its topological attributes in the network. SN(v) was the topology-based ranking score of a gene v, which was calculated from the three network topological attributes: node degree, shortest path length and node betweenness, SN(v) = (cd(v) + cspl(v) + cb(v))/3.

thumbnailFigure 6. Ranking genes by network attributes. Genes of the PPI network were sorted according to network attributes (node degree was used as an example here). The scores were then converted into ranking values.

Deriving GO enriched gene set

GO provides a controlled vocabulary of terms for describing genes and gene products [17]. GO enrichment analysis assesses whether certain GO annotations are significantly over represented among a set of genes [51]. The level of enrichment can be represented by a p-value based on specific probability distribution. The p-value is calculated by randomly picking sets of genes from the genome and computing the probability of obtaining more genes with GO terms annotated to those in the study set [52]. The smaller the p-value, the more significant the GO term is enriched in the gene set. For the KnownSet, we used three tools (DAVID [53], GOstats [54], GOEAST [55]) to perform enrichment analyses in the BP (Biological Process) sub-ontology. DAVID and GOEAST are online tools and GOstats is an R package from Bioconductor [56]. The top 50 enriched terms were picked from results obtained by each of those tools. The three enriched GO term sets from the three tools were combined into one GO term set by taking their union. In addition, cancer-hallmark related GO terms are those characteristically related to cancers, they should also be included. The set of cancer-hallmark GO terms listed in Table 1 of [37] were added to the above obtained GO term set, which was then remapped to a set of corresponding genes based on human GO annotations. Genes not covered by our PPI network were removed. The obtained set of genes was enriched with BP terms annotated to known breast cancer genes; they thus were more likely to be involved in breast cancer than randomly selected genes. This set of genes was called the GO enriched gene set (GOSet) (Figure 7).

thumbnailFigure 7. Getting GO enriched gene set by GO enrichment analysis. Three different tools were used to perform GO term enrichment analyse for known breast cancer genes. Top 50 enriched GO terms were picked from results obtained by each tool, and their union were generated. Nine cancer-hallmark GO terms from [37] were added into the enriched GO term set. The enriched GO term set were re-mapped back to a set of human genes based on Homo sapiens GO annotations, called the GO enriched gene set (GOSet).

Ranking by gene expression and GO

All breast cancer-related gene expression datasets (keywords: Homo sapiens & breast cancer) were download from the GEO database [38]. Data sets with fewer than five samples or conditions were deleted. Data sets of normal versus cancer samples were used so those containing recurred versus non-recurred samples were deleted. 53 GDSes (GEO data sets) were thus obtained.

For each GDS, records with “null” information and genes which didn’t exist in the PPI network were removed, and if a gene had more than one expression profiles, its expression was defined as the profile which had the largest mean value [57], defined as:

<a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M7">View MathML</a>

(3)

where n is the set of expression profiles for gene i in a GDS, m is the number of samples/conditions in one of those profiles, and ek(j) is the corresponding expression value of sample j.

After the above mentioned preprocessing steps, genes in each GDS were clustered by the APCluster algorithm according to their expression profiles. APCluster is an algorithm based on affinity propagation which works by considering all data points as potential cluster centers at the same time and setting up messages of similarity between any two data points, messages are exchanged among data points until all clusters are determined. APCluster had been shown to perform well compared to other clustering approaches [58,59]. Pearson correlation coefficient between gene expression profiles was used as the similarity metric for APCluster. It was assumed that genes within a cluster would have higher probability of being involved in certain biological processes than those across clusters.

Overlaps of expression clusters with the GOSet were next computed. To evaluate significance of the overlaps, the same number of genes was randomly sampled as those in the cluster from GDS, and its overlap with the GOSet was computed; the procedure was repeated 1000 times for each cluster. A p-value was then computed for the clusters (Figure 8). Rankings were assigned to clusters according to their p-values; that is, the lower the p-value, the higher the position in the ranking list. Genes in the same clusters were assigned the same ranking. A ranking score SE(i) was thus obtained for genes in each GDS, where i represented a specific GDS. A score SE was assigned for each gene by computing the average of SE(i) from all relevant GDSes.

thumbnailFigure 8. Ranking gene clusters from GEO expression profiles. For a cluster i, random samples of same size were drawn from the same GDS and their overlaps with GOSet were computed. P-value was used to represent significance of a cluster’s overlap with GOSet. The smaller the p-value, the higher the ranking.

Ranking integration

For a gene v, the ranking scores based on network topology and expression clustering were combined into an overall ranking as following,

<a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M8">View MathML</a>

(4)

SE(v) was the expression-based ranking score of gene v, which was computed from breast cancer-related gene expression data and GO annotations. λ (0≤λ≤1) is a coefficient to weigh the contribution of topological attributes and expression information in ranking breast cancer genes. The average ranking of genes in the KnownSet that sorted into top 10% was computed as the P-score. A smaller P-score meant better performance, that is, it was more likely to find true breast cancer genes from the top of the sorted list. S(v) is the final ranking for a gene v, which reflected the belief that a specific gene was a potential breast cancer gene. The higher a gene was ranked, the more likely it was involved in breast cancer related processes.

Random walk approach

The steps in [13] were followed to perform random walk. It iteratively transitions a walker from its current node to a randomly selected neighbour. Assuming W is the adjacency matrix of the PPI network and pt is a vector whose i-th element holds the probability of arriving at node i at step t. Random walk was computed by

<a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M9">View MathML</a>

(5)

where p0 is the initial probability, which is 1/37 for the 37 genes in the KnownSet and 0 for all others; r represents the probability of remaining at the same node at the next step. [30] showed that random walk worked robustly against different r values, which was also confirmed by our computation (data not shown). r was taken to be 0.7 in the current work. For details of random walk approach, see [13].

Ranking performance comparison

OMIMSet contained 26 known breast cancer genes, CGCSet contained additional 11. The 26 known breast cancer genes in the OMIMSet were used as the KnownSet. Procedure in Figure 1 was followed and the model built, which was then used to rank the 11 known breast cancer genes in the CGCSet. Ranking values in italic meant those genes were ranked in the top 10% of the final list. The row of “top 10% average” represented average rankings of those known breast cancer genes in CGCSet that ranked in top 10%, while “all 11 test genes average” represented the average rankings of all 11 genes in the CGCSet (Table 2). In later computation, OMIMSet and CGCSet were combined into a KnownSet of 37 genes.

Performance evaluation against realistic data noise

GOSet was a set of genes enriched with probable breast cancer genes. Its quality was obviously lower than the KnownSet, but higher than a set of random genes nonetheless. To evaluate the proposed method’s robustness with respect to noise in real life data sources, multiples of 37 (the number of known breast cancer genes in the PPI network) genes were sampled from the GOSet (ranging from 0 to 50, with 0 meaning no noise added). Each sample was combined with the KnownSet into an ntegrated source data. Stratified cross-validation was utilized for performance comparison between the proposed approach and the random walk method [13]. The sampled genes were randomly partitioned into five equal subsets, one of the five subsets was retained for testing, and the remaining four subsets were used as training data. The KnownSet was also randomly divided into five equal parts. One of them was combined with one subset of the sampled genes as testing data, the other four with the four remaining subsets of the sample as training data. This procedure was then repeated five times, with each of the five subsets used once as testing data. All results from the five folds were averaged to generate the final result. Performance of the proposed approach and random walk method with respect to data noise were evaluated in terms of the F-score, which was computed from precision and recall (Figure 9). Precision was the fraction of genes ranked within top k in the test data that were true known cancer genes; recall was the fraction of known breast cancer genes ranking within top k. F-score was then the harmonic mean of precision and recall. k was the ranking threshold that was used to decide whether a ranked gene was considered as a predicted positive, that is, genes ranked higher than k were judged as breast cancer genes. The F-score was defined as following,

<a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M10">View MathML</a>

(6)

where

<a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M11">View MathML</a>

(7)

<a onClick="popup('http://www.biomedcentral.com/1752-0509/7/49/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1752-0509/7/49/mathml/M12">View MathML</a>

(8)

where A is the number of genes in the KnownSet, B is the number of genes ranked within top k.

thumbnailFigure 9. Computing precision and recall. A was the set of known breast cancer genes, and B was the set of breast cancer genes which had been predicted as breast cancer genes by the proposed method. Precision represented ability to reject unrelated genes, and recall represented ability to obtain true breast cancer genes.

One F-score was computed for each fold, averaging five F-scores (for five-fold cross validation) produced the final F-score.

Performance evaluation against completely random noise

Random genes were first sampled from the PPI network and added to the KnownSet. Procedure in Figure 1 and random walk computation were then performed. Figure 5 plotted the ratio of random genes to the number of genes in KnownSet (37) and the proportion of added random genes that ranked within top 10%.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

YN and JY jointly developed the framework. YN collected data. YN and JY performed analyses. YN and JY jointly wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This work was partly supported by the National Natural Science Foundation of China (No. 61179008).

References

  1. Wang X, Gulbahce N, Yu H: Network-based methods for human disease gene prediction.

    Brief Funct Genomics 2011, 10:280-293. PubMed Abstract | Publisher Full Text OpenURL

  2. Wu X, Li S: Cancer gene prediction using a network approach.

    Cancer Systems Biology 2010, 191-212. OpenURL

  3. Siegel R, Naishadham D, Jemal A: Cancer statistics, 2012.

    CA Cancer J Clin 2012, 62:10-29. PubMed Abstract | Publisher Full Text OpenURL

  4. Materi W, Wishart DS: Computational systems biology in cancer: modeling methods and applications.

    Gene Regul Syst Bio 2007, 1:91-110. PubMed Abstract | PubMed Central Full Text OpenURL

  5. Ideker T, Sharan R: Protein networks in disease.

    Genome Res 2008, 18:644-652. PubMed Abstract | Publisher Full Text OpenURL

  6. Chuang HY, Hofree M, Ideker T: A decade of systems biology.

    Annu Rev Cell Dev Biol 2010, 26:721-744. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Oti M, Snel B, Huynen MA, Brunner HG: Predicting disease genes using protein–protein interactions.

    J Med Genet 2006, 43:691-698. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Östlund G, Lindskog M, Sonnhammer ELL: Network-based Identification of Novel Cancer Genes.

    Mol Cell Proteomics 2010, 9:648-655. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Franke L, 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:1011-1025. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Wu X, Jiang R, Zhang MQ, Li S: Network-based global inference of human disease genes.

    Mol Syst Biol 2008, 4:189. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Navlakha S, Kingsford C: The power of protein interaction networks for associating genes with diseases.

    Bioinformatics 2010, 26:1057-1063. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Navlakha S, Schatz MC, Kingsford C: Revealing biological modules via graph summarization.

    J Comput Biol 2009, 16:253-264. PubMed Abstract | Publisher Full Text OpenURL

  13. Köhler S, Bauer S, Horn D, Robinson PN: Walking the interactome for prioritization of candidate disease genes.

    Am J Hum Genet 2008, 82:949-958. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Vanunu O, Magger O, Ruppin E, Shlomi T, Sharan R: Associating genes and protein complexes with disease via network propagation.

    PLoS Comput Biol 2010, 6:e1000641. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Tiffin N, Andrade-Navarro M, Perez-Iratxeta C: Linking genes to diseases: it's all in the data.

    Genome Med 2009, 1:77. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  16. Perez-Iratxeta C, Bork P, Andrade MA: Association of genes to genetically inherited diseases using data mining.

    Nat Genet 2002, 31:316-319. PubMed Abstract | Publisher Full Text OpenURL

  17. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT: Gene Ontology: tool for the unification of biology.

    Nat Genet 2000, 25:25-29. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Nitsch D, Gonçalves J, Ojeda F, De Moor B, Moreau Y: Candidate gene prioritization by network analysis of differential expression using machine learning approaches.

    BMC Bioinforma 2010, 11:460. BioMed Central Full Text OpenURL

  19. Linh T, Bin Z, Zhan Z, Chunsheng Z, Tao X, John L, Hongyue D, Eric S, Jun Z: Inferring causal genomic alterations in breast cancer using gene expression data.

    BMC Syst Biol 2011, 5:121. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  20. George RA, Liu JY, Feng LL, Bryson-Richardson RJ, Fatkin D, Wouters MA: Analysis of protein sequence and interaction data for candidate disease gene prediction.

    Nucleic Acids Res 2006, 34:e130-e130. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. Gortzak-Uzan L, Ignatchenko A, Evangelou AI, Agochiya M, Brown KA, St. Onge P, Kireeva I, Schmitt-Ulms G, Brown TJ, Murphy J, Rosen B, Shaw P, Jurisica I, Kislinger T: A proteome resource of ovarian cancer ascites: integrated proteomic and bioinformatic analyses to identify putative biomarkers.

    J Proteome Res 2007, 7:339-351. PubMed Abstract | Publisher Full Text OpenURL

  22. Fortney K, Jurisica I: Integrative computational biology for cancer research.

    Hum Genet 2011, 4:465-481. OpenURL

  23. Nibbe RK, Koyutürk M, Chance MR: An integrative-omics approach to identify functional sub-networks in human colorectal cancer.

    PLoS Comp Biol 2010, 6:e1000639. Publisher Full Text OpenURL

  24. Aerts S, Lambrechts D, Maity S, Van Loo P, Coessens B, De Smet F, Tranchevent LC, De Moor B, Marynen P, Hassan B: Gene prioritization through genomic data fusion.

    Nat Biotechnol 2006, 24:537-544. PubMed Abstract | Publisher Full Text OpenURL

  25. Xia Z, Wu L-Y, Zhou X, Wong ST: Semi-supervised drug-protein interaction prediction from heterogeneous biological spaces.

    BMC Syst Biol 2010, 4:S6. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. De Bie T, Tranchevent LC, Van Oeffelen LMM, Moreau Y: Kernel-based data fusion for gene prioritization.

    Bioinformatics 2007, 23:i125-i132. PubMed Abstract | Publisher Full Text OpenURL

  27. Chen Y, Wang W, Zhou Y, Shields R, Chanda SK, Elston RC, Li J: In silico gene prioritization by integrating multiple data sources.

    PLoS One 2011, 6:e21137. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  28. Barabási AL, Gulbahce N, Loscalzo J: Network medicine: a network-based approach to human disease.

    Nat Rev Genet 2011, 12:56-68. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Huan T, Wu X, Bai Z, Chen JY: Seed-weighted random walks ranking method and its application to leukemia cancer biomarker prioritizations.

    Proceedings of the 2009 IEEE International Conference on Bioinformatics and Biomedicine Workshop 2009, 220-220. OpenURL

  30. Chen X, Yan G, Ren W, Qu JB: Modularized random walk with restart for candidate disease genes prioritization.

    Syst Biol (Stevenage) 2009, 353-360. OpenURL

  31. Li Y, Patra JC: Genome-wide inferring gene–phenotype relationship by walking on the heterogeneous network.

    Bioinformatics 2010, 26:1219-1224. PubMed Abstract | Publisher Full Text OpenURL

  32. Yu J, Finley RL Jr: Combining multiple positive training sets to generate confidence scores for protein–protein interactions.

    Bioinformatics 2009, 25:105-111. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Berardini TZ, Khodiyar VK, Lovering RC, Talmud P: The Gene Ontology in 2010: extensions and refinements.

    Nucleic Acids Res 2010, 38:D331-D335. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Rhee SY, Wood V, Dolinski K, Draghici S: Use and misuse of the gene ontology annotations.

    Nat Rev Genet 2008, 9:509-515. PubMed Abstract | Publisher Full Text OpenURL

  35. Harris M, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C: The Gene Ontology (GO) database and informatics resource.

    Nucleic Acids Res 2004, 32:D258-D261. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Murali T, Pacifico S, Yu J, Guest S, Roberts GG, Finley RL: DroID 2011: a comprehensive, integrated resource for protein, transcription factor, RNA and gene interactions for Drosophila.

    Nucleic Acids Res 2011, 39:D736-D743. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Li J, Lenferink AE, Deng Y, Collins C, Cui Q, Purisima EO, O'Connor-McCourt MD, Wang E: Identification of high-quality cancer prognostic markers and metastasis network modules.

    Nat Commun 2010, 1:34. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository.

    Nucleic Acids Res 2002, 30:207-210. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. 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:D514-D517. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, Rahman N, Stratton MR: A census of human cancer genes.

    Nat Rev Cancer 2004, 4:177-183. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. 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:D767-D772. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  42. Stark C, Breitkreutz BJ, Chatr-Aryamontri A, Boucher L, Oughtred R, Livstone MS, Nixon J, Van Auken K, Wang X, Shi X: The BioGRID interaction database: 2011 update.

    Nucleic Acids Res 2011, 39:D698-D704. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  43. Persico M, Ceol A, Gavrila C, Hoffmann R, Florio A, Cesareni G: HomoMINT: an inferred human network based on orthology mapping of protein interactions discovered in model organisms.

    BMC Bioinforma 2005, 6:S21. OpenURL

  44. Hermjakob H, Montecchi‒Palazzi L, Lewington C, Mudali S, Kerrien S, Orchard S, Vingron M, Roechert B, Roepstorff P, Valencia A, Margalit H, Armstrong J, Bairoch A, Cesareni G, Sherman D, Apweiler R: IntAct: an open source molecular interaction database.

    Nucleic Acids Res 2004, 32:D452-D455. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  45. Cui Q, Ma Y, Jaramillo M, Bari H, Awan A, Yang S, Zhang S, Liu L, Lu M, O'Connor-McCourt M, Purisima EO, Wang E: A map of human cancer signaling.

    Mol Syst Biol 2007, 3:152. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  46. Wang PI, Marcotte EM: It's the machine that matters: predicting gene function and phenotype from protein networks.

    J Proteomics 2010, 73:2277-2289. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  47. Xu J, Li Y: Discovering disease-genes by topological features in human protein-protein interaction network.

    Bioinformatics 2006, 22:2800-2805. PubMed Abstract | Publisher Full Text OpenURL

  48. Özgür A, Vu T, Erkan G, Radev DR: Identifying gene-disease associations using centrality on a literature mined gene-interaction network.

    Bioinformatics 2008, 24:i277-i285. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  49. Brandes U: On variants of shortest-path betweenness centrality and their generic computation.

    Social Netwks 2008, 30:136-145. Publisher Full Text OpenURL

  50. Hagberg AA, Schult DA, Swart PJ: Exploring network structure, dynamics, and function using NetworkX. Proceedings of the 7th Python in Science Conference (SciPy2008). Pasadena CA USA; 2008:11-15. OpenURL

  51. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles.

    PNAS 2005, 102:15545-15550. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Rivals I, Personnaz L, Taing L, Potier M-C: Enrichment or depletion of a GO category within a class of genes: which test?

    Bioinformatics 2007, 23:401-407. PubMed Abstract | Publisher Full Text OpenURL

  53. Sherman BT, Tan Q, Kir J, Liu D, Bryant D, Guo Y, Stephens R, Baseler MW, Lane HC, Lempicki RA: DAVID bioinformatics resources: expanded annotation database and novel algorithms to better extract biology from large gene lists.

    Nucleic Acids Res 2007, 35:W169-W175. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Falcon S, Gentleman R: Using GOstats to test gene lists for GO term association.

    Bioinformatics 2007, 23:257-258. PubMed Abstract | Publisher Full Text OpenURL

  55. Zheng Q, Wang XJ: GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis.

    Nucleic Acids Res 2008, 36:W358-W363. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  56. Davis JW: Bioinformatics and computational biology solutions using R and Bioconductor.

    J Amer Statistical Assoc 2007, 102:388-389. OpenURL

  57. Yu J, Murali T, Finley RL Jr: Assigning confidence scores to protein–protein interactions. In Two Hybrid Technologies. Springer; 2012:161-174. OpenURL

  58. Frey BJ, Dueck D: Clustering by passing messages between data points.

    Science 2007, 315:972-976. PubMed Abstract | Publisher Full Text OpenURL

  59. Bodenhofer U, Kothmeier A, Hochreiter S: APCluster: an R package for affinity propagation clustering.

    Bioinformatics 2011, 27:2463-2464. PubMed Abstract | Publisher Full Text OpenURL