Skip to main content

Multiple signatures of a disease in potential biomarker space: Getting the signatures consensus and identification of novel biomarkers

Abstract

Background

The lack of consensus among reported gene signature subsets (GSSs) in multi-gene biomarker discovery studies is often a concern for researchers and clinicians. Subsequently, it discourages larger scale prospective studies, prevents the translation of such knowledge into a practical clinical setting and ultimately hinders the progress of the field of biomarker-based disease classification, prognosis and prediction.

Methods

We define all "gene identificators" (gIDs) as constituents of the entire potential disease biomarker space. For each gID in a GSS of interest ("tested GSS"/tGSS), our method counts the empirical frequency of gID co-occurrences/overlaps in other reference GSSs (rGSSs) and compares it with the expected frequency generated via implementation of a randomized sampling procedure. Comparison of the empirical frequency distribution (EFD) with the expected background frequency distribution (BFD) allows dichotomization of statistically novel (SN) and common (SC) gIDs within the tGSS.

Results

We identify SN or SC biomarkers for tGSSs obtained from previous studies of high-grade serous ovarian cancer (HG-SOC) and breast cancer (BC). For each tGSS, the EFD of gID co-occurrences/overlaps with other rGSSs is characterized by scale and context-dependent Pareto-like frequency distribution function. Our results indicate that while independently there is little overlap between our tGSS with individual rGSSs, comparison of the EFD with BFD suggests that beyond a confidence threshold, tested gIDs become more common in rGSSs than expected. This validates the use of our tGSS as individual or combined prognostic factors. Our method identifies SN and SC genes of a 36-gene prognostic signature that stratify HG-SOC patients into subgroups with low, intermediate or high-risk of the disease outcome. Using 70 BC rGSSs, the method also predicted SN and SC BC prognostic genes from the tested obesity and IGF1 pathway GSSs.

Conclusions

Our method provides a strategy that identify/predict within a tGSS of interest, gID subsets that are either SN or SC when compared to other rGSSs. Practically, our results suggest that there is a stronger association of the IGF1 signature genes with the 70 BC rGSSs, than for the obesity-associated signature. Furthermore, both SC and SN genes, in both signatures could be considered as perspective prognostic biomarkers of BCs that stratify the patients onto low or high risks of cancer development.

Background

Current technology encourages the study of biological phenomena on a genome-wide scale. Technological platforms such as microarrays, next-generation sequencing, and mass spectrometry have resulted in generation of data on an unprecedented scale [1, 2]. Inadvertently, the field of bioinformatics which includes high-performance cloud computing, adaptation of statistical methods, design of novel algorithms and generation of databases, play critical roles in the analysis of these massive and diverse datasets [3]. The variation in the type and amount of biological data, coupled with the fact that investigators may sometimes be confronted with a question that cannot be answered using current statistical techniques or algorithms [4], means that the field of statistical methods and algorithms is under constant refinement, adaptation and improvement [5].

Today, analysis of data from high-throughput experiments often yields a set of high-dimensional variable (HDV) list which typically represent a particular phenotype with respect to another. Such HDV lists commonly include signature lists of expressed genes, loci or proteins. Subsequent types of analysis to be performed on the gene list, depend greatly on the biological question an investigator is interested in. The work has been greatly simplified, partly due to the presence of many databases that were created, mostly in recent years [6–8]. The wealth of raw or curated, but nonetheless collated information in these databases is often critical in the subsequent analysis of gene (or other HDV) lists derived from these high-throughput experiments.

One of the most common analyses one could perform with a set of gene lists is an enrichment study of biological functions, processes or pathways with respect to a well-annotated reference gene list which commonly includes all the annotated genes in the genome. This analysis is commonly termed gene ontology analysis [9] which is based on simple statistical tests such as hypergeometric, binomial, or Chi-square tests [10]. These statistical tests could also be used if one is merely interested in whether one list of genes is similar to another, e.g. whether the gene products differentially expressed in human breast cancer (BC) are similar to the gene products differentially expressed in human ovarian cancer [8]. In addition, complementary methods such as Gene Set Enrichment Analysis (GSEA) allow the assessment of the relative relevance of one gene list of interest with reference to the expression differences of ranked genes between two phenotype cell classes of an organism [11].

Despite improvements in experimental technology and techniques, poor reproducibility and stability of results from independent but similar experiments can often hinder scientific discovery. These issues can arise due to many reasons such as small sample size [12, 13], high-noise data [12–14], use of different technological platforms as well as poorly reported clinical or research protocols (different cohort classifications, treatment differences) [8, 14–16]. In particular, small sample sizes, in combination with technical and biological noises, often complicate efforts to identify statistical differences of expression signals between many functionally important genes of distinct tumor subtypes or clinical groups. These limitations lead to bias in signature predictions and poor consistency. Inconsistency, divergence and poor overlap of many dozen of reported signatures suggest that our knowledge of nature and space dimensionality of tumor-associated genes and potential biomarkers is essentially incomplete [13, 14, 16]. Identification of potential biomarker space can reveal specific genetic patterns of cancer cells, for example, cell junctions in non-small cell lung cancer subtypes [17]. Our recent integrative studies of microarray gene expression profiles in lung adenocarcinoma (AC) versus normal adjacent lung tissue suggested that space dimensionality of potential biomarkers of lung adenocarcinoma (AC) is at least 2300 known genes [13]. Among these, hundreds of genes could be essential for disease-driving because they encode proteins containing mutagenesis sites, implying that these genes could be considered as relevant for diagnostic or prognostic applications [13]. In BC, the number of the genes that grade primary tumor by its aggressiveness is even larger; it consists of ~4000 microarray U133A detected genes [16].

To the best of our knowledge, current statistical analysis of matched lists mostly centres on comparison of two lists [18, 19]. For example, independent studies of prostate cancer from two different countries may each yield a set of gene lists where genes are ranked by biological relevance, either by the magnitude of fold change, statistical significance of differential expression or other statistical measures. The comparison of the two gene lists based on robustness and stability can then be evaluated using a recently published method, which incorporates permutation studies and uses Canberra distance as the measure of dissimilarity [18]. However, methods which evaluate elements of one (new) list with reference to many other (known) lists is required in certain disciplines, including bioinformatics, genome-associated disease studies, medical statistics, epidemiology, ecology and authorship identification etc [20].

In cancer research, advances in high-throughput experimental techniques as well as statistical algorithms have resulted in the discovery of many gene signature sets (GSSs) which are representative of a particular cancer phenotype based on patient prognosis, tumor subtypes or other molecular features. However, it was reported that many of the GSSs generally do not show strong consensus for a given disease, even for more clinically homogeneous sub-groups of a disease (e.g. stage, tumor subtype) when independent patient cohorts are compared [12, 13, 21]. In BC, two well-known signatures which predict BC recurrence, comprising of 79 and 76 genes were derived independently in Amsterdam and Rotterdam cohorts respectively [15, 22]. However, only three genes were found to be common. Nevertheless, it was stated that the rest of the unique genes could share similar pathways and be associated with similar mechanisms leading to the disease.

Recently, a systematic collection of gene expression signatures in cancers have been identified and collected in several databases [6–8]. For example, Abba et al. have collected 42 BC gene expression signatures in an effort to identify the most relevant BC biomarkers [23]. Comparison of these 42 signatures revealed limited or zero overlap between signatures. Specifically, comparison of the 3427 distinct gene symbols revealed that only 15 genes (RRM2, MELK, MAD2L1, MYBL2, BIRC5, PTTG1, AURKA, PRC1, CKS2, CDCA8, MKI67, UBE2C, DUSP4, CENPF and CDC2) are found in at least 10 signatures, which indicates the great disparity across gene signatures. The reasons for the disparity have been attributed to differences in clinical attributes of patients analysed which include ER status, stages, histological patterns, disease subtypes and treatment received by the patients.

Also, it is likely that each of the signatures represented only a partial picture of the heterogeneous and complex BC biology, which subsequently limits its potential for clinical implementation. The authors demonstrated that selection of several of these signatures having better consensus may provide more relevant and robust cancer biomarkers [23]. However, their strategy has a bias towards the larger size signatures and over-representation of cell cycle related genes. Other authors published a systematic review of gene expression signatures in colorectal cancer and identified 31 prognostic gene signatures [24]. It was reported that these gene lists, comprising a total of 1530 genes do not show great overlap, as there were only two common genes in four signatures, 10 common genes in three signatures, and 102 common genes in two signatures. It was stated by the authors that "the lack of gene overlap is generally interpreted as if each signature is a random sampling of a small subset of genes from a larger signature that represent the involved pathways [24]." However, without strong agreement and reproducibility of the gene signatures from independent studies, future large-scale prospective validation studies are unlikely to proceed, and may prove to be a major hindrance in achieving the desired goals of biomarker-based disease diagnosis, prognosis and prediction of therapeutic efficacy.

Many of the available biomarkers were generated in connection with biological functions of a given medical condition/disease. However, identity of the biomarkers for the same medical condition/disease could be quite different due to differences in the study design or analytical methods. Generally, it is very difficult to compare the GSSs due to such poorly-controlled variations in the process of discovery. However at times, it would be interesting to know how many times a particular gID in a newly derived signature has been reported in other known and reference GSSs (rGSSs) and whether the presence of the gID in these rGSSs occurs more frequently than expected by chance.

In our work, we address medical bioinformatics issues via computational intelligence, statistical analysis and computational simulation. Here, we proposed that a randomized sampling approach can provide a confidence indication of whether a gene is likely to be a common (and perhaps relatively more reliable) potential biomarker present in many other GSSs, or whether it is a relatively unique biomarker specific to a particular biological or disease state.

Results

Definition of novel or common biomarkers

Traditional definition of novel biomarkers typically involves an absolute threshold value of 0, where only biomarkers which are not present in other published reference gene signature subsets (rGSSs) are defined as novel biomarkers (Figure 1A). In contrast, biomarkers which are present in at least one published rGSS are defined as "known" biomarkers. However, this definition does not account for the number of rGSSs under comparison, as well as ignore the variation in the number of biomarkers across rGSSs. For example, a biomarker present in only one of hundreds of rGSSs might be considered as a "statistically novel (SN) biomarker".

Figure 1
figure 1

Definition of novel or common biomarkers. (A) Traditional definition of novel or common biomarkers. (B) Statistical definition of novel or common biomarkers. A further vertical dimension is extended which provides a statistical measure of whether the signature gene is considered "novel".

To address the above-mentioned issues, we first propose the construction of the observed frequency distribution function that describes the number of published rGSSs that contain each of the genes from the gene signature of interest. This function has a skewed shape with long right side tail. More precisely, the observed frequencies have the following characteristics in common: there are few frequent, and many rare events (clusters, interactions, co-occurrence etc). Such skewed functions are often observed in many natural and technological processes (the birth-death processes, biological evolution, interaction events in genome, transcriptome and proteome scales, artificial complex systems, physics phenomena, biological and social networks, industry incidences [25, 26]. Sampling from such populations could be commonly fitted by the Pareto-like frequency distribution function, which is sample-size and context-dependent [25, 26]. In practical applications, the left side of the observed skewed distribution could be enriched with 'admixture' events which consist of 'null' or 'background' additive noise events due to error measurements. In comparison with the Pareto-like frequency distribution function, such additional ('admixture') null or background frequency distribution (BFD) has a relatively shorter right-side tail. Such function could be described by the exponential distribution function or Poison distribution function. The identification of BFD and 'de-noising' of the empirical Pareto-like distribution function is a great challenge, specifically when sample sizes are relatively small [25, 26].

We propose a simulation-based approach of random sampling to generate an expected BFD of the number of published rGSSs expected by chance to contain each of the genes from the GSS of interest (or "tested GSS"/tGSS) (Figure 1B). The method of random sampling to generate the BFD provides a fair basis of evaluating a significance measure of the empirical frequency of the number of the published rGSSs that contain each of the genes from the signature of interest. Implicitly, it takes into account the number of rGSSs under comparison, as well as taking into account the size variation across rGSSs. Effectively, it extends a vertical dimension (of expected and actual frequencies of the published rGSSs) to facilitate identification of statistically novel (SN) or common (SC) biomarkers (Figure 1B).

To discriminate between SN and SC biomarkers within tGSS derived from genome-scale data, we assume that the identity of individual element (potential biomarker represented by gene and/or probesets ID) is sufficiently well recorded as a potential biomarker, if that element appeared in the other lists more than r times (Figure 1B). With this assumption, we can compare the null BFD with actual empirical frequency distribution (EFD) and at the given confidence level, estimate a critical cut-off value of signature co-occurrences which subsequently allows us to discriminate biomarkers as SN or SC.

Description of method

An illustration of our proposed methodology is shown in Figure 2. We first define a background gene list which is a superset of all gene signatures. Appropriately, a background gene list could be defined as all human gene symbols or a subset of genes such as only those represented on a particular microarray platform. Genes in signatures which are not in the background gene list are removed. For each gene in a newly derived gene signature (AS0), the number of co-occurrences in M other published rGSSs is counted (ASi = 1,2,3...M) (Figure 2A). The EFD of number of genes for each co-occurrence events with the published rGSSs can be counted and plotted.

Figure 2
figure 2

Schema of gene list comparison with other defined sets. (A) Actual observations of gene lists overlap between single list of interest (AS 0 ) with other defined sets. (B) Observations of gene lists overlap in a simulation where other defined gene sets are randomly and independently sampled without replacement. AS and RS denote actual and random set respectively. O m and RO m denotes overlap segments and random overlap segments respectively. Blue solid circle represents our gene list of interest (AS 0 ). Green oval, red rectangle and yellow triangle represent 3 other defined set of genes with sizes |AS i = 1 |, |AS i = 2 |, |AS i = 3 | respectively.

To generate the null BFD, several rounds of simulation are performed. In each simulation, random gene lists RSi = 1,2,3...M of equal sizes (where |RSj,i| = |ASi| for gene lists i = 1,2,3...M) are generated from the background gene list, without replacement within each simulation and with replacement between two simulations (Figure 2B). The number of co-occurrences of overlap between genes in AS0 with RSi = 1,2,3...M are counted in each simulation, and summed across all simulations. Subsequently, a null BFD of the number of gID matches with the random gene lists can be generated, which represents the expected frequency distribution of gID matches.

Effect of background gene list

We first performed studies on simulated data to establish the family of null distributions that may result from this analysis. Specifically, the effect of the size of the background gene list, i.e. the effect of biomarker space on the null BFD is of interest. The background gene list, D, is first assumed to comprise of 20,000 genes. We assume that our gene signature, AS0 contains 100 genes, i.e. |AS0| = 100. Next, we assume that 10 other well-defined gene signatures are reported in the literature, each containing an unequal number of genes, e.g. |AS1|=20, |AS2|=40, |AS3|=60, |AS4|=80, |AS5|=100, |AS6|=120, |AS7|=140, |AS8|=160, |AS9|=180 and |AS10|=200.

Subsequently, we performed 100 simulations where in each simulation, the actual gene sets (ASi) are simulated by sampling the same number of genes independently from D without replacement. The expected frequency of observed co-occurrences of each gene in our gene signature can be determined.

Similar analyses are performed by repeating the procedures with reduced background gene lists (|D| = 10000, 5000, 1000, 500 genes). The effect of the size of the background gene lists can be observed in Figure 3. In addition, the expected null BFD can be approximated and fitted using Weibull or Sigmoid functions.

Figure 3
figure 3

Family of null frequency distribution of expected co-occurrences of our signature genes with other signatures. The horizontal axis represents the number of samples that contain the gene from our signature of interest. The dotted lines represent the fitted curves of Weibull function whereas the dashed lines represent the fitted curves of Sigmoid function

Analysis of 36-genes prognostic signature for epithelial ovarian cancer

In our previous work [27], we studied expression and clinical data from patients diagnosed with high-grade serous ovarian cancer, and subsequently identified a 36-mRNA signature (assigned as the tGSS in this analysis) which stratifies patients into subgroups with very distinct and varied overall survival rates: low, intermediate or high-risk, with 5 year overall survival rates of 65%, 20% and 10% respectively. In addition to its prognostic significance, the 36-mRNA signature is also predictive of patients' response to chemo-therapy.

We compared our 36-gene tGSS with other published rGSSs of ovarian cancer. From the literature, we collected 63 rGSSs which were previously reported to show associations with survival, disease subtype, chemo-sensitivity, disease detection, development, progression or recurrence (Additional file 1). We restrict our analysis to official gene symbols present in the background gene set (RefGene version 30th November 2012). After pre-processing of the gene symbols, each of the 63 rGSSs contain between 1 to 966 gene symbols. Subsequently, via 100 independent simulations, we generated a null BFD of overlap between our tGSS (comprising of 36 genes) with the randomly generated lists of signatures. Additionally, to understand the effect of the number of simulations on the null BFDs, we performed 1000 and 10000 independent simulations.

The results of the comparison of our tGSS with other rGSSs are shown in Table 1 and Figure 4A (see also Additional file 3).

Table 1 Analysis of occurrence events for genes in the ovarian cancer prognostic gene signature set.
Figure 4
figure 4

Actual and expected frequency distribution of gene overlap from a query signature with other reference signatures. Comparison of genes from (A) 36-gene ovarian cancer prognostic gene signature, (B) tumor breast obesity gene signature and (C) tumor breast IGF1 gene signature, with other reference gene signatures for that disease. (D) Comparison of the actual frequency distribution generated from tumor breast obesity (From B) and tumor breast IGF1 (From C). The expected frequency distributions were generated via performing N simulations, where N is 100, 1000 or 10000. The y-axis is log10 transformed. p1 denotes the two-sided p-value from Kolmogorov-Smirnov statistic which tests if the actual and expected (for N = 100) distribution are similar. p2 denotes the p-value that represents the significance of that threshold in dichotomizing statistically novel or common biomarkers from a GSS of interest.

Analysis of obesity and IGF1 signatures in breast cancer

Recently, the link between obesity and its contribution to poorer disease outcome in BC has been reported. Creighton and Sada et al. studied the effects of obesity on primary breast tumor gene expression and their results revealed an obesity-associated cancer transcriptional signature of 662 genes (assigned as the tGSS in this analysis) [28]. After preprocessing, this tGSS contains 683 RefSeq gene symbols (Additional file 2). Subsequently, we investigate whether individual genes in this obesity-associated tGSS derived from BC tumors are significantly enriched among the 70 BC rGSSs which were previously reported to be associated with clinical observations such as response to chemotherapy, distant metastasis, ER-alpha status, tumor subtypes and grades, as well as clinical outcomes such as patient prognosis. The 70 rGSSs comprise 42 gene signatures reviewed by Abba et al. [23] as well as those manually curated by us (Additional file 2). The background gene list consists of all RefGene symbols downloaded from UCSC Genome Browser on 30th November 2012 and any gene that is not found within the background gene list is excluded from subsequent analysis.

As described in the methods, we performed 100 independent simulations and the enrichment of the genes in rGSSs relative to random expectation are shown in Figure 4B (Additional file 4). Overall, our results indicated that a large proportion of the obesity-associated tGSS genes have been reported as biomarkers in more BC rGSSs than expected from random simulations. For instance, based on random simulations, only 13.0% of the obesity-associated signature is expected to be found in at least 2 other rGSSs. However, 161 of 683 (23.6%) genes from the obesity-associated tGSS are actually found in at least 2 other rGSSs.

Also, it has been suggested that the link between obesity and BC outcome could be due to increased endocrine signaling involving insulin and insulin-like growth factors (IGFs) [28]. Specifically IGF1, whose expression is elevated in human breast cancer [29], is known to increase breast cancer cell growth and invasion [30]. Furthermore, activation of IGF1 pathway is correlated with early recurrence and decreased relapse-free survival [31].

Therefore, in subsequent analysis, we studied an IGF1 pathway gene signature (assigned as the tGSS in this analysis) whose genes were differentially expressed in MCF-7 BC cell line after IGF1 stimulation and which were correlated with several poor prognostic factors and disease outcome in patients with BC [32]. This IGF1 pathway tGSS contains 925 genes after gene identifier conversion and preprocessing. Similar analyses were performed for this tGSS with respect to the 70 other BC rGSSs and the results are shown in Figure 4C (Additional file 5). Similar to the obesity-associated tGSS, our results indicated that genes belonging to the IGF1-tGSS are generally found to co-occur in more rGSSs than expected from random simulations. Specifically, 43.6% of the IGF1-signature genes are found in at least 2 rGSSs, which is 3.3 times more than expected.

When we compared the EFDs for both the obesity-associated tGSS and IGF1 pathway tGSS, results suggested that there is a stronger association of the IGF1 pathway tGSS's genes with the 70 BC rGSSs than observed for the obesity-associated tGSS's genes with the 70 BC rGSSs (Figure 4D).

Next, we studied both the obesity and IGF1 tGSSs of BC with respect to the MAPK signalling pathway. From the KEGG database, the MAPK signalling pathway comprises of 267 genes, which included a list of 31 genes (ATF2, CHUK, DDIT3, DUSP4, DUSP5, DUSP6, DUSP8, FGF13, FGFR2, FGFR4, FLNA, FOS, GADD45G, HRAS, IKBKB, IL1R1, MAP2K4, MAP2K6, MAP3K8, MAP4K3, MAPK9, MAPKAPK5, MAX, PAK2, PPM1B, PPP3R1, RAPGEF2, RASA1, RPS6KA3, RPS6KA4 and SOS1) which are present in either or both of obesity-associated tGSS and the IGF1 pathway tGSS (Table 2). Subject to a minimum threshold of four rGSSs, genes present in at least one of the obesity or IGF1 tGSSs could be considered as SC gIDs. On the other hand, genes in either the obesity or IGF1-associated tGSSs that occur in three or less rGSSs are termed SN gIDs (Table 2). Typically, such genes which individually occur less frequently in other rGSSs are not easily interpreted with respect to their association with disease pathways or phenotypes. On the other hand, our method identifies many dozen genes of MAPK signalling pathway which have not been considered as a common BC signature in context of their functional association of obesity, IGF1 pathway and BC.

Table 2 Occurrence of MAP kinases signalling pathway genes in the studied breast tissue/cell signatures.

Discussion

In this work, we developed a method for the identification of gene subsets that are statistically novel (SN) or common (SC) in a newly defined signature of interest when compared to other (known) reference gene signature sets (rGSSs) that are representative of a medical condition (e.g., breast cancer).

We first provide an application of our methodology to one of our gene signature previously identified for the prognosis of patients diagnosed with high-grade serous ovarian cancer [27]. Our prognostic signature comprise of 36 mRNA genes (Signature#1 of Additional file 1). For this analysis, we manually curate 63 gene sets from the literature which were reported to show associations with survival, disease subtype, chemo-sensitivity, disease detection, development, progression or recurrence (Additional file 1). Our analysis revealed that more than 50% of the genes in our 36-gene prognostic signature were not previously considered as potential biomarkers in any of the curated publications (Table 1). These genes include ARPC1B, CCT2, CD93, CDK4, DNMT1, FGFR1, FZD1, HGF, MIS12, NCAPD2, NCAPG2, NCAPH, PIK3R1, PLAUR, POLA2, POLR2D, POLR2J, TGFBR2 and VCL. Despite the absence of these genes in currently known gene sets associated with the different aspects of ovarian cancers, some of these were nonetheless shown to play critical roles in cancer development and progression. For instance, it was shown that the majority of epithelial ovarian cancer tumors exhibited positive staining for FZD1 [33], and it was associated with chemo-resistance via the Wnt/Beta-catenin pathway [34, 35]. Similarly for HGF, deregulated HGF/MET signaling is a common hallmark of many tumors and is associated with various aspects of tumor progression [36]. Furthermore, elevated HGF serum levels could predict poor prognosis in advanced ovarian cancers [37]. To the best of our knowledge, although FZD1 and HGF have not been reported as high-confidence biomarkers in ovarian cancer, their values as prognostic markers should not be neglected as is often the case. In fact in our previous analysis of high-grade serous ovarian cancer patients, both FZD1 and HGF exhibited significant prognostic properties (Figure 5A-B).

Figure 5
figure 5

Classification of high-grade serous ovarian cancer patients. The patients diagnosed with high-grade serous ovarian cancer were classified using a data-driven method for statistically novel biomarkers (A) FZD1and (B) HGFand common biomarkers (C) COL3A1and (D) EDNRA. Log-rank tests were used to assess the survival statistical significance of the two patient subgroups. Expr: expression.

On the other hand, genes such as CCL2, PDGFRA, TUBB, COL3A1 and EDNRA have previously been identified as relevant biomarkers in at least 3 other published gene sets (Table 1). For instance, EDNRA, either independently or in combination with other biomarkers, was reported to be able to predict benign and malignant tumors from borderline tumors [38], tumor subtype classification [39, 40] and classification of tumors related to cell plasticity [41]. COL3A1 is also a commonly studied gene in ovarian cancer, where it was revealed that it is one of the most expressed proteins in advanced relative to local ovarian adenocarcinoma [42], and that its expression was observed to be higher in platinum-resistant relative to platinum-sensitive cells [43]. Its use as a biomarker in prediction of chemo-sensitivity or platinum-sensitivity [43, 44], determination of tumor subtype [39] or molecular subtype [40], and classification of tumors related to cell plasticity [41] was previously reported. Our data also indicated that both EDNRA and COL3A1 showed significant prognostic properties in patients diagnosed with high-grade serous ovarian cancer (Figure 5C-D).

In fact, our signature was originally derived for patient prognosis and prediction of chemo-sensitivity. In this analysis, we further found that certain genes in our 36-gene signature could be also relevant in tumor subtype classification. Also, we show that survival curves of several of our biomarkers show significance in patient stratification, regardless of whether they have been or not been reported in previous publications of ovarian cancer (Figure 5).

Next, we also studied two BC signatures from two published research reports (Signatures#1-2 of Additional file 2). The first signature was derived from differential expression analysis of obese versus non-obese BC patients [28]. The second signature was derived from differentially regulated genes in MCF-7 cells after IGF1 stimulation [32]. Independent comparison of these two signatures with the 70 BC rGSSs revealed that a larger proportion of genes from the IGF1 pathway gene signature were more likely to be represented in the rGSSs when compared to expectations from random simulation (Figure 4D, Additional files 3, 4). This might suggest a tighter association between the IGF1 signaling pathway with clinical observations such as response to chemotherapy, distant metastasis, ER-alpha status, tumor subtypes and grades, as well as clinical outcomes such as patient prognosis. In contrast, while we showed slight association between the obesity-associated gene signature with the rGSSs (Figure 4B,Additional file 4), this association appear to be weaker when compared to that of the IGF1 signaling pathway (Figure 4C, Additional file 5). This seems to suggest that the probable effects of obesity on cancer association could be less specific or directed, when compared to that of the IGF1 signaling pathway. Our method allows us to quantify a measure of associations between obesity, IGF1 pathway and BC signatures, as well as predict potentially SC and SN therapeutic target genes.

Additionally, using our published data-driven prognostic analytical method [45], we studied the survival significance of SN gIDs (PIK3C3 and APPBP2) and SC gIDs (IL6ST and DUSP6) of BC (Figure 6). Expression levels of SC gIDs such as IL6ST and DUSP6, which were found in 7 and 8 BC rGSSs respectively (Additional files 4 and 5), were found to be able to stratify both Stockholm and Uppsala BC patient cohorts into two survival significant subgroups via the data-driven grouping method (Figures 6C-D). Genes that were less studied and less represented in other rGSSs were traditionally neglected in terms of their prognostic capability. However, we show that expression levels of genes such as PIK3C3 and APPBP2, which were not found in any of the 70 BC rGSSs (Additional files 4 and 5), could also stratify both Stockholm and Uppsala patient cohort into two survival significant subgroups (Figures 6A-B).

Figure 6
figure 6

Classification of breast cancer patients for both Stockholm and Uppsala patient cohort. The patients diagnosed with breast cancer were classified using data-driven method for statistically novel biomarkers (A) PIK3C3and (B) APPBP2and common biomarkers (C) IL6STand (D) DUSP6. Top panel: Stockholm breast cancer patient cohort, Bottom panel: Uppsala breast cancer patient cohort. Log-rank tests were used to assess the survival statistical significance of the two patient subgroups. Expr: expression.

While our general schema allows us to uncover the null background frequency distribution (BFD) of expected co-occurrences of our tSS's genes with other well-defined rGSSs, the specific use of this approach requires a careful formulation of the null hypothesis. Therefore, the background gene set D needs to be judiciously defined. In our example studies of BC and HG-SOC, we have defined the background gene set D as all annotated RefSeq gene symbols. However, sometimes it may be necessary to restrict the set of background gene list D to the list of genes that are tested, e.g. probesets on a microarray platform. The effect of a smaller background gene set D is non-trivial. Our simulated experiments on reduced background gene sets (|D| = 20000, 10000, 5000, 1000, 500) clearly revealed the rightward shift of the null BFD as the background gene set becomes more restricted (Figure 3). The interpretation of the actual empirical frequency distribution (EFD) and null BFD is therefore context-dependent and depends on the selection of the background gene list and other rGSSs. Therefore, the selection of the background set is non-trivial, and should generally be guided by the design of the experiments or formulation of the hypotheses.

We also studied the effect of the number of simulations on the null BFD. For each of the three tGSSs analysed, we generated null BFD from 100, 1000 or 10000 random and independent simulations (see Methods). Our results show that generally, 100 simulations are sufficient to uncover the structure and shape of the null BFD, which is representative of the percentage of tGSS that is expected to be found in the rGSSs at various levels of overlap counts. Increasing the number of simulations (1000 or 10000) might increase the sensitivity of our method in detecting random genes that likely appear in more rGSSs, but this is at the expense of higher computational cost.

Finally, in order to obtain a meaningful expected null BFD, the use of the method requires a considerable size of the signature of interest (|AS0|) and the number of other well-defined rGSSs (M). Therefore, the rapidly growing number of databases as well as the accumulating wealth at databases including but not limited to Gene Expression Omnibus (GEO) [6], ArrayExpress [7], GeneSigDB [8], meant that analysis of a newly generated gene lists with these previously defined gene sets could become more feasible. Certainly, the selection of the previously defined gene sets from the public repositories or literatures requires careful consideration as the eventual null BFD is representative of a scenario subject to the assumptions and conditions.

Conclusions

The accumulating wealth of molecular signatures in publications and public repositories for a wide range of disease types and experimental conditions is a useful resource for any researcher. Importantly, it allows projection of a newly-derived GSS onto published, defined and well-annotated GSSs and subsequently identifies individual gIDs that are also independently associated with another feature. Our randomized sampling approach allows the generation of a null background frequency distribution, which could be used to identify a preliminary co-occurrence threshold and subsequently select subsets of the signature genes as statistically unique or significantly associated with another molecular feature. We have successfully applied our method to expression microarray and clinical data. In particular, our results suggest that there is a stronger association of the IGF1 signature genes with 70 BC rGSSs, than for the obesity-associated signature. Also, both SN and SC gIDs could be considered as perspective prognostic biomarkers of BCs. Our method identifies many dozen genes which have not been considered in context of their functional association of obesity and BC. In particular, our analysis suggests a close association of many genes expressed in MAPK pathways with obesity, IGF1 and BC pathways. We propose that the method reported in this study could also be applicable and of interest to researchers in other fields including but not limited to social science, epidemiology, linguistics, forensic analysis and internet networks.

Methods

Collection of gene sets

In our previous work, we identified 36 mRNA genes that are associated with overall survival and chemotherapeutic response of patients diagnosed with high-grade serous ovarian carcinoma [27]. This 36-gene signature is used as the testing signature of interest (or tGSS) for ovarian cancer (Signature#1 of Additional file 1). For breast cancer, we downloaded and analyzed two independent gene sets associated with obesity and IGF1 pathway respectively [28, 32]. These two gene sets are defined as our testing signatures (or tGSS) of interest for BC (Signatures#1-2 of Additional file 2).

For the reference gene signature subsets (rGSSs), we compiled gene sets associated with BC and ovarian cancer from the literature. A total of 70 and 63 gene signature sets were collected for breast and ovarian cancers respectively. For BC gene sets, they were reported to predict clinical observations such as response to chemotherapy, distant metastasis, ER-alpha status, tumor subtypes and grades, as well as clinical outcomes such as patient prognosis (Additional file 2). The ovarian cancer gene sets were known to show associations with survival, disease subtype, chemo-sensitivity, disease detection, development, progression or recurrence (Additional file 1).

All gene accession identifiers are converted to official gene symbols.

Definition of biomarker space

In this work, we defined the biomarker space as the official gene symbols represented on the RefGene database (downloaded from the UCSC genome browser on 30th November 2012).

Frequency distribution of actual gene occurrences in reference gene subsets

The visual representation as well as the notations of our methodology is shown in Figure 2. Assume that the entire biomarker space, D is well-defined and contains |D| genes. For example, a probable biomarker space could be defined by the number of reliable microarray probe-sets used for expression signal detection. Let ASi denote the actual set of genes, where i = 0 for a newly derived GSS (e.g. our gene signature or "testing GSS"/tGSS), and i = 1,2,3...M for M other reference GSSs (rGSSs), each containing a list of genes with size |ASi|>0. We assume that |ASi|<<|D| for i = 1,2,3...M.

Within our actual gene set AS0, the number of genes that are found in m other gene sets are denoted by observations, Om = 1,2,3...M, where

O m = 1 , 2 , 3 … M = ∑ s = s e g m e n t a l l s e g m e n t s w i t h m o b s e r v a t i o n s O m , s
(1)

Also, the number of genes in our gene set (AS0) can be denoted by |AS0| where

A S 0 = ∑ m = 1 M O m = O 1 + O 2 + O 3 +…+ O M
(2)

The observed normalized frequency Om = 1,2,3...M is denoted by Fm = 1,2,3...m where

F m = 1 , 2 , 3 … m = O m A S 0
(3)

Generation of random reference gene subsets via random sampling

For jth of N simulations, RSi,j denotes a randomly and independently generated set of genes with size equal to ASi, i.e. |RSj,i| = |ASi| for gene lists i = 1,2,3...M. Each RSj,i is sampled randomly without replacement from the biomarker space D. Thus, the lists are produced by randomly sampling without replacement from biomarker space D, but the intensity of sampling might differ from list to list. In statistical literature such random collection was called a multiple record system [46].

The procedure is repeated N times, for N simulations.

Frequency distribution of expected gene occurrences in reference gene subsets

Next, the background model of the expected biomarker co-occurrence across other gene signatures is developed via sampling and simulations.

For jth simulation, the number of genes in our gene set (AS0) that are also found in m gene sets RSj,i = 1,2,3...M is denoted by ROj,m = 1,2,3,...M, where

R O j , m = 1 , 2 , 3 … M = ∑ s = s e g m e n t a l l s e g m e n t s w i t h m o b s e r v a t i o n s R O j , m , s
(4)

The expected normalized frequency of ROm = 1,2,3,...M after N simulations is

E F m = 1 , 2 , 3 … M = ∑ j = 1 N R O j , m N × A S 0
(5)

Comparison of the expected and actual frequency distribution of gene overlaps with other signatures

Subsequently, the observed and expected frequency of gene overlaps with other signatures can be seen, by comparing Fm and EFm for m = 1,2,3,...M.

Statistics of differences and identification of threshold stratifying statistically novel and common biomarkers of a signature

Initially, to assess the statistical differences between the observed normalized frequency (F) and expected normalized frequency (EF), we use Kolmogorov-Smirnov statistics to assess the statistical differences between the reverse cumulative frequency of F and EF.

When the difference between F and EF is significant, we aim to discriminate between statistically novel (SN) or common (SC) gIDs within a GSS of interest, by assuming that the biomarker is sufficiently well recorded as a potential marker if that biomarker appears in other rGSSs more than r times (Figure 1B).

For each discrete value r, where r ≥ 0, the p-value representing the significance of that threshold in stratifying SN or SC gIDs from a GSS of interest can be calculated from the ratio of the cumulative frequencies (expected with respect to actual) at that discrete value.

Abbreviations

AC:

Adenocarcinoma

BC:

Breast cancer

BFD:

Background frequency distribution

EFD:

Empirical frequency distribution

gID:

Gene identificator

HDV:

High-dimensional variable

HG-SOC:

High-grade serous ovarian carcinoma

ID:

Identifier

IGF1:

Insulin-like growth factor 1 (somatomedin C)

GEO:

Gene Expression Omnibus

GSS:

Gene signature subset

ncRNA:

Non-coding RNA

RefSeq:

Reference sequence

SC:

Statistically common

SN:

Statistically novel

rGSS:

Reference GSS

References

  1. Chin L, Hahn WC, Getz G, Meyerson M: Making sense of cancer genomic data. Genes & development. 2011, 25 (6): 534-555. 10.1101/gad.2017311.

    Article  CAS  Google Scholar 

  2. Lizardi PM, Forloni M, Wajapeyee N: Genome-wide approaches for cancer gene discovery. Trends Biotechnol. 2011, 29 (11): 558-568. 10.1016/j.tibtech.2011.06.003.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Fortney K, Jurisica I: Integrative computational biology for cancer research. Hum Genet. 2011, 130 (4): 465-481. 10.1007/s00439-011-0983-z.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Li Y, Chen L: Big Biological Data: Challenges and Opportunities. Genomics, Proteomics & Bioinformatics. 2014, 10.1016/j.gpb.2014.10.001.

    Google Scholar 

  5. Wang Y, Zhang XS, Chen L: Computational systems biology in the big data era. BMC Syst Biol. 2013, 7 (Suppl 2): S1-10.1186/1752-0509-7-S2-S1.

    Article  PubMed Central  PubMed  Google Scholar 

  6. Barrett T, Suzek TO, Troup DB, Wilhite SE, Ngau WC, Ledoux P, Rudnev D, Lash AE, Fujibuchi W, Edgar R: NCBI GEO: mining millions of expression profiles--database and tools. Nucleic Acids Res. 2005, 33 (Database): D562-566.

    PubMed Central  CAS  PubMed  Google Scholar 

  7. Parkinson H, Sarkans U, Shojatalab M, Abeygunawardena N, Contrino S, Coulson R, Farne A, Lara GG, Holloway E, Kapushesky M, et al: ArrayExpress--a public repository for microarray gene expression data at the EBI. Nucleic Acids Res. 2005, 33 (Database): D553-555.

    PubMed Central  CAS  PubMed  Google Scholar 

  8. Culhane AC, Schroder MS, Sultana R, Picard SC, Martinelli EN, Kelly C, Haibe-Kains B, Kapushesky M, St Pierre AA, Flahive W, et al: GeneSigDB: a manually curated database and resource for analysis of gene expression signatures. Nucleic Acids Res. 2011, 40 (Database): D1060-1066.

    PubMed Central  PubMed  Google Scholar 

  9. 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.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Rivals I, Personnaz L, Taing L, Potier MC: Enrichment or depletion of a GO category within a class of genes: which test?. Bioinformatics. 2007, 23 (4): 401-407. 10.1093/bioinformatics/btl633.

    Article  CAS  PubMed  Google Scholar 

  11. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Ein-Dor L, Zuk O, Domany E: Thousands of samples are needed to generate a robust gene list for predicting outcome in cancer. Proc Natl Acad Sci USA. 2006, 103 (15): 5923-5928. 10.1073/pnas.0601231103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Toh SH, Prathipati P, Motakis E, Kwoh CK, Yenamandra SP, Kuznetsov VA: A robust tool for discriminative analysis and feature selection in paired samples impacts the identification of the genes essential for reprogramming lung tissue to adenocarcinoma. BMC Genomics. 2011, 12 (Suppl 3): S24-10.1186/1471-2164-12-S3-S24.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Clarke R, Ressom HW, Wang A, Xuan J, Liu MC, Gehan EA, Wang Y: The properties of high-dimensional data spaces: implications for exploring gene and protein expression data. Nat Rev Cancer. 2008, 8 (1): 37-49. 10.1038/nrc2294.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Wang Y, Klijn JG, Zhang Y, Sieuwerts AM, Look MP, Yang F, Talantov D, Timmermans M, Meijer-van Gelder ME, Yu J, et al: Gene-expression profiles to predict distant metastasis of lymph-node-negative primary breast cancer. Lancet. 2005, 365 (9460): 671-679. 10.1016/S0140-6736(05)70933-8.

    Article  CAS  PubMed  Google Scholar 

  16. Chua ALS, Ivshina AV, Kuznetsov VA: Pareto-Gamma Statistics reveals global rescaling in transcriptomes of low and high aggressive breast cancer phenotypes. Pattern Recognition in Bioinformatics (PRIB-2006). Edited by: Ragapakese LW, R. Acharya. 2006, LNCS 4146: Springer-Verlag Berlin-Heidelberg, 49-59. 10.1007/11818564_7.

    Chapter  Google Scholar 

  17. Kuner R, Muley T, Meister M, Ruschhaupt M, Buness A, Xu EC, Schnabel P, Warth A, Poustka A, Sultmann H, et al: Global gene expression analysis reveals specific patterns of cell junctions in non-small cell lung cancer subtypes. Lung Cancer. 2009, 63 (1): 32-38. 10.1016/j.lungcan.2008.03.033.

    Article  PubMed  Google Scholar 

  18. Jurman G, Riccadonna S, Visintainer R, Furlanello C: Algebraic comparison of partial lists in bioinformatics. PLoS One. 2012, 7 (5): e36540-10.1371/journal.pone.0036540.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Damavandi B: Estimating the Overlap of Top Instances in Lists Ranked by Correlation to Label. MS thesis. 2012, University of Alberta, [http://hdl.handle.net/10402/era.24985]

    Google Scholar 

  20. Ow GS, Jenjaroenpun P, Thiery JP, Kuznetsov VA: How to discriminate between potentially novel and considered biomarkers within molecular signature?. 2013 IEEE Symposium on Computational Intelligence in Bioinformatics and Computational Biology (CIBCB): 2013; Singapore. 2013, IEEE Publishing, 176-182. 10.1109/CIBCB.2013.6595405.

    Chapter  Google Scholar 

  21. Gormley M, Dampier W, Ertel A, Karacali B, Tozeren A: Prediction potential of candidate biomarker sets identified and validated on gene expression data from multiple datasets. BMC Bioinformatics. 2007, 8: 415-10.1186/1471-2105-8-415.

    Article  PubMed Central  PubMed  Google Scholar 

  22. van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, et al: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.

    Article  PubMed  Google Scholar 

  23. Abba MC, Lacunza E, Butti M, Aldaz CM: Breast cancer biomarker discovery in the functional genomic age: a systematic review of 42 gene expression signatures. Biomark Insights. 2010, 5: 103-118.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Sanz-Pamplona R, Berenguer A, Cordero D, Riccadonna S, Sole X, Crous-Bou M, Guino E, Sanjuan X, Biondo S, Soriano A, et al: Clinical value of prognosis gene expression signatures in colorectal cancer: a systematic review. PLoS One. 2012, 7 (11): e48877-10.1371/journal.pone.0048877.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Kuznetsov V: Scale-Dependent Statistics of the Numbers of Transcripts and Protein Sequences Encoded in the Genome. Computational and Statistical Approaches to Genomics. Edited by: Zhang W, Shmulevich I. 2006, Springer US, 163-208. 10.1007/0-387-26288-1_10. 2

    Chapter  Google Scholar 

  26. Kuznetsov VA, Singh O, Jenjaroenpun P: Statistics of protein-DNA binding and the total number of binding sites for a transcription factor in the mammalian genome. BMC Genomics. 2010, 11 (Suppl 1): S12-10.1186/1471-2164-11-S1-S12.

    Article  PubMed Central  PubMed  Google Scholar 

  27. Tang Z, Ow GS, Thiery JP, Ivshina AV, Kuznetsov VA: Meta-analysis of transcriptome reveals let-7b as an unfavorable prognostic biomarker and predicts molecular and clinical sub-classes in high-grade serous ovarian carcinoma. Int J Cancer. 2013, 134 (2): 306-318.

    Article  CAS  PubMed  Google Scholar 

  28. Creighton CJ, Sada YH, Zhang Y, Tsimelzon A, Wong H, Dave B, Landis MD, Bear HD, Rodriguez A, Chang JC: A gene transcription signature of obesity in breast cancer. Breast Cancer Res Treat. 2012, 132 (3): 993-1000. 10.1007/s10549-011-1595-y.

    Article  CAS  PubMed  Google Scholar 

  29. Papa V, Pezzino V, Costantino A, Belfiore A, Giuffrida D, Frittitta L, Vannelli GB, Brand R, Goldfine ID, Vigneri R: Elevated insulin receptor content in human breast cancer. J Clin Invest. 1990, 86 (5): 1503-1510. 10.1172/JCI114868.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  30. Sarkissyan S, Sarkissyan M, Wu Y, Cardenas J, Koeffler HP, Vadgama JV: IGF-1 regulates Cyr61 induced breast cancer cell proliferation and invasion. PLoS One. 2014, 16 (7): e103534-

    Article  Google Scholar 

  31. Champ CE, Volek JS, Siglin J, Jin L, Simone NL: Weight gain, metabolic syndrome, and breast cancer recurrence: are dietary recommendations supported by the data?. Int J Breast Cancer. 2012, 2012: 506868-

    Article  PubMed Central  PubMed  Google Scholar 

  32. Creighton CJ, Casa A, Lazard Z, Huang S, Tsimelzon A, Hilsenbeck SG, Osborne CK, Lee AV: Insulin-like growth factor-I activates gene transcription programs strongly associated with poor breast cancer prognosis. J Clin Oncol. 2008, 26 (25): 4078-4085. 10.1200/JCO.2007.13.4429.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Badiglian Filho L, Oshima CT, De Oliveira Lima F, De Oliveira Costa H, De Sousa Damiao R, Gomes TS, Goncalves WJ: Canonical and noncanonical Wnt pathway: a comparison among normal ovary, benign ovarian tumor and ovarian cancer. Oncol Rep. 2009, 21 (2): 313-320.

    PubMed  Google Scholar 

  34. Flahaut M, Meier R, Coulon A, Nardou KA, Niggli FK, Martinet D, Beckmann JS, Joseph JM, Muhlethaler-Mottet A, Gross N: The Wnt receptor FZD1 mediates chemoresistance in neuroblastoma through activation of the Wnt/beta-catenin pathway. Oncogene. 2009, 28 (23): 2245-2256. 10.1038/onc.2009.80.

    Article  CAS  PubMed  Google Scholar 

  35. Zhang H, Zhang X, Wu X, Li W, Su P, Cheng H, Xiang L, Gao P, Zhou G: Interference of Frizzled 1 (FZD1) reverses multidrug resistance in breast cancer cells through the Wnt/beta-catenin pathway. Cancer Lett. 2012, 323 (1): 106-113. 10.1016/j.canlet.2012.03.039.

    Article  CAS  PubMed  Google Scholar 

  36. Zhou HY, Pon YL, Wong AS: HGF/MET signaling in ovarian cancer. Curr Mol Med. 2008, 8 (6): 469-480. 10.2174/156652408785747933.

    Article  CAS  PubMed  Google Scholar 

  37. Aune G, Lian AM, Tingulstad S, Torp SH, Forsmo S, Reseland JE, Stunes AK, Syversen U: Increased circulating hepatocyte growth factor (HGF): a marker of epithelial ovarian cancer and an indicator of poor prognosis. Gynecol Oncol. 2011, 121 (2): 402-406. 10.1016/j.ygyno.2010.12.355.

    Article  CAS  PubMed  Google Scholar 

  38. Curry EW, Stronach EA, Rama NR, Wang YY, Gabra H, El-Bahrawy MA: Molecular subtypes of serous borderline ovarian tumor show distinct expression patterns of benign tumor and malignant tumor-associated signatures. Mod Pathol. 2013

    Google Scholar 

  39. Bentink S, Haibe-Kains B, Risch T, Fan JB, Hirsch MS, Holton K, Rubio R, April C, Chen J, Wickham-Garcia E, et al: Angiogenic mRNA and microRNA gene expression signature predicts a novel subtype of serous ovarian cancer. PLoS One. 2012, 7 (2): e30269-10.1371/journal.pone.0030269.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Tan TZ, Miow QH, Huang RY, Wong MK, Ye J, Lau JA, Wu MC, Bin Abdul Hadi LH, Soong R, Choolani M, et al: Functional genomics identifies five distinct molecular subtypes with clinical relevance and pathways for growth control in epithelial ovarian cancer. EMBO Mol Med. 2013, 5 (7): 983-998.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. De Cecco L, Marchionni L, Gariboldi M, Reid JF, Lagonigro MS, Caramuta S, Ferrario C, Bussani E, Mezzanzanica D, Turatti F, et al: Gene expression profiling of advanced ovarian cancer: characterization of a molecular signature involving fibroblast growth factor 2. Oncogene. 2004, 23 (49): 8171-8183. 10.1038/sj.onc.1207979.

    Article  CAS  PubMed  Google Scholar 

  42. Gagne JP, Ethier C, Gagne P, Mercier G, Bonicalzi ME, Mes-Masson AM, Droit A, Winstall E, Isabelle M, Poirier GG: Comparative proteome analysis of human epithelial ovarian cancer. Proteome Sci. 2007, 5: 16-10.1186/1477-5956-5-16.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Helleman J, Jansen MP, Span PN, van Staveren IL, Massuger LF, Meijer-van Gelder ME, Sweep FC, Ewing PC, van der Burg ME, Stoter G, et al: Molecular profiling of platinum resistant ovarian cancer. Int J Cancer. 2006, 118 (8): 1963-1971. 10.1002/ijc.21599.

    Article  CAS  PubMed  Google Scholar 

  44. Cheng L, Lu W, Kulkarni B, Pejovic T, Yan X, Chiang JH, Hood L, Odunsi K, Lin B: Analysis of chemotherapy response programs in ovarian cancers by the next-generation sequencing technologies. Gynecol Oncol. 2010, 117 (2): 159-169. 10.1016/j.ygyno.2010.01.041.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Motakis E, Ivshina AV, Kuznetsov VA: Data-driven approach to predict survival of cancer patients: estimation of microarray genes' prediction significance by Cox proportional hazard regression model. IEEE Eng Med Biol Mag. 2009, 28 (4): 58-66. 10.1109/MEMB.2009.932937.

    Article  PubMed  Google Scholar 

  46. Lloyd CJ: Statistical Analysis of Categorical Data. Wiley Series in Probability and Statistics. 1999, John Wiley & Sons, Inc, 367-

    Google Scholar 

Download references

Acknowledgements

We thank the reviewers for their valuable comments. We acknowledge the funding of this work by the Agency of Science, Technology and Research (A*STAR), Singapore. The funding agency had no role in the study design, data analysis, and decision to publish.

Declarations

The publication costs for this article were funded by A*STAR, Singapore.

This article has been published as part of BMC Genomics Volume 16 Supplement 7, 2015: Selected articles from The International Conference on Intelligent Biology and Medicine (ICIBM) 2014: Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/16/S7.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Vladimir A Kuznetsov.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

VAK conceived and designed the study and together with OGS developed the algorithm of proposed methods. OGS implemented the methods. All authors analyzed the data, interpreted the results and wrote the manuscript.

Electronic supplementary material

Additional file 1: Gene sets compiled from published studies of ovarian cancer. (XLSX 38 KB)

Additional file 2: Gene sets compiled from published studies of breast cancer. (XLSX 75 KB)

12864_2015_7196_MOESM3_ESM.xlsx

Additional file 3: Analysis of occurrence events for genes in the ovarian cancer gene signature set. Expected and actual frequency of co-occurrences of 36 genes in the ovarian cancer signature with other 60 gene signatures. All signatures are ovarian cancer-related. (XLSX 10 KB)

12864_2015_7196_MOESM4_ESM.xlsx

Additional file 4: Analysis of occurrence events for genes in the tumor breast obesity gene signature set. Expected and actual frequency of co-occurrences of 683 genes in the tumor breast obesity signature with other 70 gene signatures. All signatures are breast cancer-related. (XLSX 13 KB)

12864_2015_7196_MOESM5_ESM.xlsx

Additional file 5: Analysis of occurrence events for genes in the tumor breast IGF1 gene signature set. Expected and actual frequency of co-occurrences of 925 genes in the tumor breast IGF1 signature with other 70 gene signatures. All signatures are breast cancer-related. (XLSX 14 KB)

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ow, G.S., Kuznetsov, V.A. Multiple signatures of a disease in potential biomarker space: Getting the signatures consensus and identification of novel biomarkers. BMC Genomics 16 (Suppl 7), S2 (2015). https://doi.org/10.1186/1471-2164-16-S7-S2

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-16-S7-S2

Keywords