<title>Background</title> <p>Many studies have demonstrated the ability of using gene-expression "signatures" derived from DNA microarray data to define cancer subtypes, predict disease recurrence, and guide treatment decisions. In breast cancer, van't Veer <italic>et al</italic>. [##REF##11823860##1##] derived a 70-gene profile to predict a patient's risk of developing distant metastases. Perou <italic>et al</italic>. [##REF##10963602##2##] and Sorlie <italic>et al</italic>. [##REF##11553815##3##] developed an intrinsic-subtype signature that classifies breast tumors into molecular subtypes showing distinct differences in prognosis. From a cancer biology perspective, Chang <italic>et al</italic>. [##UREF##0##4##] studied the links between the wound healing process and cancer progression. Based on the expression pattern of a wound-response signature of 512 genes, they classify a tumor to have either activated or quiescent response and found this to be a significant prognostic predictor of tumor metastasis. These are promising results and a few of these signatures have begun to be assessed in clinical settings. Two questions have often been asked: (1) are these signatures identifying the same set of genes and (2) will they generate similar prediction performance when tested in new data sets?</p> <p>The answer to the first question has been discouraging. Any pair of these signatures share only a few common genes. Possible reasons have been suggested including the differences in patient cohort characteristics (such as the distribution of age or stage of the disease), lack of comparability and reproducibility of the data generated using different microarray platforms, and varying statistical procedures used to generate the gene list. Nevertheless, Ein-Dor <italic>et al</italic>. [##REF##15308542##5##] showed that the inconsistency still exists when eliminating all three differences. In particular, the authors repeated the same analysis in a single data set and identified many lists of genes equally predictive of the outcome. Any two of these gene lists share only a small number of genes. In another study by Son <italic>et al</italic>. [##REF##15741514##6##], the authors reported that any randomly selected subgroup of around 100 differentially expressed genes generates similar hierarchical clustering results in the same data set.</p> <p>Ein-Dor <italic>et al</italic>. [##REF##16585533##7##] further suggested perhaps the main source of the problem lies in the small sample size and large number of genes the signatures were derived from. For several published breast cancer data sets, the authors estimated that several thousands of samples would be needed to achieve a typical gene overlap of 50%. On the other hand, the problem is compounded by analyzing and interpreting genes in isolation. A common approach to gene selection involves selecting a handful of top-ranking genes that best differentiate sample classes (such as tumor vs. normal tissue) or are most predictive of clinical outcome. The univariate selection procedure ignores correlation between genes. The biological and statistical validity of such assumption seems tenuous. As a result, gene-set based approaches have emerged in recent years to identify sets of biologically related genes that are deregulated as a group. Examples of gene-set analysis include the Gene Set Enrichment Analysis (GSEA) [##REF##16199517##8##], Significance Analysis of Function and Expression (SAFE) [##REF##15647293##9##], and the globaltest package [##REF##14693814##10##] These methods focus on groups of genes that share common biological functions such as cell cycle regulation; metabolic or signaling pathways defined by Gene Ontology (GO); online databases such as BioCarta, KEGG and signaling data base; or a literature-defined gene set subject to experimental perturbations such as a drug treatment or an oncogene-activation. In addition, Rhodes <italic>et al</italic>. [##REF##17534450##11##] introduced a Molecular Concepts Map (MCM) providing an expanded analytic framework to explore the network of relationships among biologically related gene sets.</p> <p>The motivation of this study came from a recent paper by Fan <italic>et al</italic>. [##REF##16899776##12##], which addresses the second question described above. The authors demonstrated a high degree of prediction concordance of five breast cancer gene-signatures despite minimal gene-wise overlap. In an independent data set of 295 tumors, the authors showed that the intrinsic subtypes [##REF##11553815##3##] of basal-like, HER2+/ER-, and luminal B were consistently classified as poor 70-gene profile [##REF##15598354##13##] prognosis, activated wound response [##UREF##0##4##] and high recurrence score [##REF##15591335##14##]. It raises the question that perhaps the gene-overlap is not the most relevant measure of robustness and reproducibility of the gene-signatures. There may be common themes shared across these signatures that are not apparent on the individual gene level. As an example, the cell cycle gene Cyclin E1 (CCNE1) was included in the 70-gene profile while Cyclin E2 (CCNE2) in the intrinsic subtype signature. The two signatures apparently share commonality in the activation of the Cyclin family genes. For another example, ERBB2 and EGFR are both receptor tyrosine kinase involved in estrogen pathway. Inclusion of one or the other in two different signatures apparently converges at the pathway level both indicating the activation of the estrogen-signaling pathway.</p> <p>In this study, we assess the potential functional convergence of these gene-signatures on the basis of activated oncogenic pathways. This involves first annotating each gene-signature to identify significantly enriched functional modules (e.g., cell growth, response to estrogen, myb-regulated pathways, etc.). Definition of the modules can be based on Gene Ontology (GO) terms, online pathway databases such as BioCarta and KEGG, or literature-defined concepts. In the next step, the overlapping functional modules are obtained by intersecting the annotated sets. We investigated six breast cancer signatures (four of which were compared in Fan <italic>et al</italic>. [##REF##16899776##12##]) that share high prediction concordance. We found eighteen common modules including estrogen-signaling, responses to tamoxifen treatment, and BRCA1 expression. The degree of the functional overlap across the six BR-signatures is highly significant (P = 0.0002) under a bootstrapped null distribution.</p>
<title>Methods</title> <p>Table ##TAB##0##1## lists the six BR-signatures that are compared in this study. Fan <italic>et al</italic>. [##REF##14693814##10##] showed high prediction concordance of signature 1–3 and signature 5. In addition, a 90-gene meta-signature [##REF##15598354##13##] is included. This signature was derived in a meta-analysis framework by integrating four microarray data sets, which included the van't Veer data set and the Sorlie data set. Another signature included here is the subset of 60-gene profile from Wang <italic>et al</italic>. [##UREF##1##18##] that was derived in tumors with estrogen receptor (ER) positive status. The recurrence-score signature [##REF##15591335##14##] is also an ER+ disease signature that has been shown in a clinical trial to be able to identify patients with very low risk of recurrence on hormone therapy using tamoxifen alone, and do not require adjuvant chemotherapy.</p> <title>Annotation</title> <p>The set of signature genes were annotated using two different annotation sources:</p> <title>Literature-defined module</title> <p>A collection of 661 literature-defined modules from the Molecular Concept (MCM) database MCM that focuses on human cancer studies. These include gene sets from peer-reviewed publications using microarrays to study gene expression changes subject to experimental perturbation such as drug treatment or candidate gene activation.</p> <title>Pathway module</title> <p>The functional subsets from the molecular signature database or MSigDB GSEA, including modules representing metabolic and signaling pathways imported from online pathway databases such as BioCarta [##UREF##2##19##], signalling pathway database [##UREF##3##20##] and the Kyoto Encyclopedia of Genes and Genomes (KEGG) [##UREF##4##21##].</p> <p>Enrichment analysis was performed using hypergeometric tests. In particular, the procedure tests the significance of the proportion of module genes (e.g., estrogen pathway) in the signature being greater than the "population"-proportion of the module genes in the experimental set from which the signature was selected. Multiple testing was adjusted by using the Benjamini-Hochberg procedure [##UREF##5##22##].</p> <title>Notation and methods</title> <p>For a set of <italic>K </italic>gene-signatures, let <italic>n<sub><italic>i </italic></sub></italic>be the number of genes in signature <italic>i </italic>and <italic>N<sub><italic>i </italic></sub></italic>be the total number of genes in the experimental set from which the signature genes were selected. Furthermore, let <italic>J </italic>= 661 or 552 denote the number of literature-defined concepts and the number of metabolic and signaling pathways in the two annotation database MCM and MsigDB respectively. For a gene signature, we first perform a module enrichment analysis using a hypergeometric test. As mentioned earlier, the basic idea is to test whether the proportion of the module genes in the signature of size <italic>n<sub><italic>i </italic></sub></italic>is significantly larger than the "population"-fraction of the module genes in the experimental set of size <italic>N<sub><italic>i</italic></sub></italic>. The <italic>j</italic>th module is enriched in the <italic>i</italic>th signature if the hypergeometric test p-value is less than 0.3. Across the <italic>K </italic>signatures under comparison, this threshold correspond to a p-value of less than 0.05 under a conventional meta-analysis of combining the hypergeometric p-values  across the four signatures based on a chi-square distribution with 2<italic>K </italic>degrees of freedom. Let <italic>X<sub><italic>ij </italic></sub></italic>be the indicator variable where <italic>X</italic><sub><italic>ij </italic></sub>= 1 if the <italic>j</italic>th module is enriched in the <italic>i</italic>th (<italic>i </italic>= 1, ..., <italic>K</italic>) signature and <italic>X</italic><sub><italic>ij </italic></sub>= 0 otherwise. As a result,</p> <p></p> <p>is the total number of enriched modules in signature <italic>i</italic>. Then for the set of <italic>K </italic>signatures, the amount of functional overlap is .</p> <p>The significance of overlap is defined as <italic>P </italic>(<italic>Y </italic>&gt; <italic>y</italic><sup><italic>obs</italic></sup>) under a bootstrapped null distribution. The bootstrap procedure is described elsewhere [see Additional File ##SUPPL##3##4##].</p> <p>We used <italic>B </italic>= 100,000 in the procedure. The bootstrapped null distribution of <italic>Y </italic>preserves 1) potential correlation of the signature size <italic>n<sub><italic>i </italic></sub></italic>and the number of enriched modules <italic>m<sub><italic>i</italic></sub></italic>, and 2) the module-module dependence due to the one-to-many mapping of a gene to the annotation data.</p>
<title>Results and Discussion</title> <title>Prediction concordance across five breast-cancer gene-signatures</title> <p>In a similar fashion as in [##REF##16899776##12##], we cross-tabulated the prediction results of the gene-signatures listed in Table ##TAB##0##1## in the 295 breast cancer patients in the van de Vijver study [##REF##12490681##15##]. In Table ##TAB##1##2##, all the signatures consistently classified the basal-like and HER2/ER- subtype tumors as having high risk of recurrence outcome. The 70-gene profile and wound-response signatures both classify luminal B subtype to be a low risk group, while the meta-signature classifies the luminal A and the normal-like subtypes as low risk groups. Overall, the signatures showed a certain degree of prediction concordance. The kappa coefficient measuring the classification agreement across the signatures is estimated to be 0.67.</p> <title>Common "oncogenic" sets underpinning breast cancer outcome prediction</title> <p>For pairs of the six signatures, there is a fair amount of overlapping literature concepts (MCMs). Many of the overlaps are highly significant (Figure ##FIG##0##1A##). For example, there is a set of 142 enriched MCM modules shared between the 70-gene profile and the wound-response signature (<italic>P </italic>&lt; 0.00001) while only two genes were identified by both. Furthermore, signatures 1–3 showed marginal significance in metabolic and signaling pathway overlaps (Figure ##FIG##0##1B##).</p> <p>We found 18 common MCMs (<italic>P </italic>= 0.0002) and 5 common metabolic and signaling pathways (<italic>P </italic>= 0.04) across signatures 1–4. Table ##TAB##2##3## and ##TAB##3##4## list these common sets ordered by the overall significance of enrichment (summarized hypergeometric test P-value adjusted for multiple testing). Among the top are deregulated genes in androgen-sensitive prostate cancer cell lines in response to MSA (MCM 258), Myb-regulated transcriptional changes in the estrogen-dependent human breast cancer cell line MCF7 (MCM 458), several MCMs comprising responsive genes upon antiestrogen hormonal treatment (MCM 691, 379, 375, 673). Clearly a dominant common characteristic underpinning the four breast-cancer signatures is closely related to the estrogen-receptor status of the tumor which is a main prognostic factor in breast cancer. Another common prognostic set of interest is response to BRCA1 expression (MCM513), which many studies have shown a characteristic of sporadic basal-like cancers.</p> <p>Table ##TAB##3##4## listed the five common metabolic and signaling pathways using the functional subset of the MsigDB annotation data. All of the signatures apparently enlisted genes customized on a commercial array platform that represent the breast cancer estrogen signaling pathway [see Additional File ##SUPPL##0##1##].</p> <p>For the gene signatures listed in Table ##TAB##0##1##, it should be pointed out that they were constructed using different types of endpoints, along with differing supervised learning algorithms. In attempting to combine results across the signatures, we make the assumption that there exists an underlying tumorigenetic mechanism that manifests itself in terms of the endpoints used by the different authors. One such mechanism might be tumor metastasis.</p> <title>ER-positive relapse signatures</title> <p>Both ER+ relapse-signatures showed evidence of E2F activation, response to Interleukin-6 (IL6), and activation of insulin-signaling pathways, some of which have been reported in the literature to be specific to ER+ disease [see Additional Files ##SUPPL##1##2## and ##SUPPL##2##3##]. For example, studies have shown in estrogen-sensitive breast cancer cell lines, the widely used antiestrogen tamoxifen treatment inhibits insulin-signaling. The degree of such inhibition can reflect the effectiveness of the tamoxifen treatment and thus correlate with a patient's risk of recurrence [##REF##9205064##16##,##REF##18245484##17##].</p>
<title>Results and Discussion</title> <title>Prediction concordance across five breast-cancer gene-signatures</title> <p>In a similar fashion as in [##REF##16899776##12##], we cross-tabulated the prediction results of the gene-signatures listed in Table ##TAB##0##1## in the 295 breast cancer patients in the van de Vijver study [##REF##12490681##15##]. In Table ##TAB##1##2##, all the signatures consistently classified the basal-like and HER2/ER- subtype tumors as having high risk of recurrence outcome. The 70-gene profile and wound-response signatures both classify luminal B subtype to be a low risk group, while the meta-signature classifies the luminal A and the normal-like subtypes as low risk groups. Overall, the signatures showed a certain degree of prediction concordance. The kappa coefficient measuring the classification agreement across the signatures is estimated to be 0.67.</p> <title>Common "oncogenic" sets underpinning breast cancer outcome prediction</title> <p>For pairs of the six signatures, there is a fair amount of overlapping literature concepts (MCMs). Many of the overlaps are highly significant (Figure ##FIG##0##1A##). For example, there is a set of 142 enriched MCM modules shared between the 70-gene profile and the wound-response signature (<italic>P </italic>&lt; 0.00001) while only two genes were identified by both. Furthermore, signatures 1–3 showed marginal significance in metabolic and signaling pathway overlaps (Figure ##FIG##0##1B##).</p> <p>We found 18 common MCMs (<italic>P </italic>= 0.0002) and 5 common metabolic and signaling pathways (<italic>P </italic>= 0.04) across signatures 1–4. Table ##TAB##2##3## and ##TAB##3##4## list these common sets ordered by the overall significance of enrichment (summarized hypergeometric test P-value adjusted for multiple testing). Among the top are deregulated genes in androgen-sensitive prostate cancer cell lines in response to MSA (MCM 258), Myb-regulated transcriptional changes in the estrogen-dependent human breast cancer cell line MCF7 (MCM 458), several MCMs comprising responsive genes upon antiestrogen hormonal treatment (MCM 691, 379, 375, 673). Clearly a dominant common characteristic underpinning the four breast-cancer signatures is closely related to the estrogen-receptor status of the tumor which is a main prognostic factor in breast cancer. Another common prognostic set of interest is response to BRCA1 expression (MCM513), which many studies have shown a characteristic of sporadic basal-like cancers.</p> <p>Table ##TAB##3##4## listed the five common metabolic and signaling pathways using the functional subset of the MsigDB annotation data. All of the signatures apparently enlisted genes customized on a commercial array platform that represent the breast cancer estrogen signaling pathway [see Additional File ##SUPPL##0##1##].</p> <p>For the gene signatures listed in Table ##TAB##0##1##, it should be pointed out that they were constructed using different types of endpoints, along with differing supervised learning algorithms. In attempting to combine results across the signatures, we make the assumption that there exists an underlying tumorigenetic mechanism that manifests itself in terms of the endpoints used by the different authors. One such mechanism might be tumor metastasis.</p> <title>ER-positive relapse signatures</title> <p>Both ER+ relapse-signatures showed evidence of E2F activation, response to Interleukin-6 (IL6), and activation of insulin-signaling pathways, some of which have been reported in the literature to be specific to ER+ disease [see Additional Files ##SUPPL##1##2## and ##SUPPL##2##3##]. For example, studies have shown in estrogen-sensitive breast cancer cell lines, the widely used antiestrogen tamoxifen treatment inhibits insulin-signaling. The degree of such inhibition can reflect the effectiveness of the tamoxifen treatment and thus correlate with a patient's risk of recurrence [##REF##9205064##16##,##REF##18245484##17##].</p>
<title>Conclusion</title> <p>Cancer gene-expression signatures derived from microarray experiments are beginning to be tested in clinical trials, while the exact biology that enables these gene-signatures to accurately predict tumor metastasis and patient survival is unclear. Microarray experiments are often limited in power by the small number of samples used to derive a panel of prognostic genes relative to the large number of features on the array. In addition, sets of biologically-related genes are often co-regulated while many feature selection procedures are univariate in nature. As a result, gene-signatures developed by different studies typically share very few common components. A recent study showed high prediction concordance of several breast cancer gene-signatures despite minimal overlap in gene identity. It gave main motivation to investigate common oncogenic themes that may not be apparent at the individual gene level. This study explored this hypothesis by evaluating the functional overlap of the signatures on the basis of annotated gene sets. When the gene signatures are mapped to the deregulated pathway space, two things become clear. First, there is a significant degree of functional overlap in oncogenic and prognostic pathways. Second, many of these common pathways provide plausible explanations of tumor biology through which these signatures predict patient outcome. There are several conclusions to be gleaned from this study. First, this work explains why independent signatures that appear to perform equally well at predicting patient prognosis show minimal overlap in gene membership. This is because such genes are different members of pathways and processes that are relevant to prognosis. Thus, the lack of gene overlap found between the various signatures listed in Table ##TAB##0##1## should not be considered problematic. The implication of our study is that most of these signatures will do well in clinical trials given that they seem to be picking up the same pathway signals. We can thus be assured that the gene lists found by different investigators are consistent, even if they do not contain the same genes.</p> <p>Second, the results have suggested that the interpretability and delineation of how diverse cancer gene expression signatures work are more likely attainable at the pathway level rather than the individual gene level. On the other hand, as many studies have already suggested so, feature selection methods need to be based on biologically related gene sets that are deregulated as a group [##REF##16199517##8##, ####REF##15647293##9##, ##REF##14693814##10##, ##REF##17534450##11####17534450##11##]. However, it is not a straightforward task to construct a prognostic signature based on pathways that are composed of overlapping sets of genes. New statistical methods need to be established in this area. This is beyond the scope of the study and is currently under investigation.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>A recent study has shown high concordance of several breast-cancer gene signatures in predicting disease recurrence despite minimal overlap of the gene lists. It raises the question if there are common themes underlying such prediction concordance that are not apparent on the individual gene-level. We therefore studied the similarity of these gene-signatures on the basis of their functional annotations.</p> <title>Results</title> <p>We found the signatures did not identify the same set of genes but converged on the activation of a similar set of oncogenic and clinically-relevant pathways. A clear and consistent pattern across the four breast cancer signatures is the activation of the estrogen-signaling pathway. Other common features include BRCA1-regulated pathway, reck pathways, and insulin signaling associated with the ER-positive disease signatures, all providing possible explanations for the prediction concordance.</p> <title>Conclusion</title> <p>This work explains why independent breast cancer signatures that appear to perform equally well at predicting patient prognosis show minimal overlap in gene membership.</p>
<title>Competing interests</title> <p>The authors declare that they have no competing interests.</p> <title>Authors' contributions</title> <p>DG and RS conceived the method and prepared the manuscript. RS performed the analyses. AC contributed to the discussion. All authors have read and approved the final manuscript.</p> <title>Pre-publication history</title> <p>The pre-publication history for this paper can be accessed here:</p> <p><ext-link ext-link-type="uri" xlink:href="http://www.biomedcentral.com/1755-8794/1/28/prepub"/></p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We would like to thank D. Rhodes and S. Kalyana-Sundaram for providing the MCM data set. RS, DG, and AMC participated in the conception and design of the study. RS performed the analysis and drafted the manuscript. DG and AMC reviewed the manuscript. RS is supported in part by NCI 2 P30 CA008748-43; DG is supported in part by NIH grant GM72007 and the Huck Institute for Life Sciences; AMC is supported by a Clinical Translational Science Award from the Burroughs Welcome Foundation.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Pair-wise functional overlap of the six breast cancer gene-signatures</bold>. 1. 70-gene profile 2. Wound response 3. Intrinsic subtype 4. Meta90 5. Recurrence score 6. Wang ER+ profile. A. The number of overlapping literature-defined oncogenic concepts (MCM) and the corresponding P-value heatmap indicating the significance of the overlap under bootstrapped null distribution. B. The number of overlapping pathway sets (MsigDB) and the corresponding P-value heatmap.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Breast cancer gene-signatures.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>Gene-signature</bold></td><td align="center"><bold>Number of genes</bold></td><td align="center"><bold>number of samples</bold></td><td align="center"><bold>Experiment summary</bold></td></tr></thead><tbody><tr><td align="left">1. 70-gene profile [##REF##11823860##1##]</td><td align="center">70</td><td align="center">78</td><td align="left">Inkjet oligonucleotide array on 25,000 genes</td></tr><tr><td align="left">2. Wound-response [##UREF##0##4##]</td><td align="center">512</td><td align="center">50</td><td align="left">cDNA microarrays profiled over 36,000 genes</td></tr><tr><td align="left">3. Intrinsic subtype [##REF##10963602##2##,##REF##11553815##3##]</td><td align="center">427</td><td align="center">78</td><td align="left">cDNA microarrays on a core set of 8,102 genes</td></tr><tr><td align="left">4. meta-90 [##REF##17534450##11##]</td><td align="center">90</td><td align="center">305</td><td align="left">Integrative analysis of 4 microarray studies on a set of 2,555 genes</td></tr><tr><td align="center" colspan="3">ER+ signature</td><td/></tr><tr><td align="left">5. Recurrence score [##REF##15591335##14##]</td><td align="center">21</td><td align="center">2892</td><td align="left">RT-PCR on 250 genes selected from the literature</td></tr><tr><td align="left">6. Wang ER+ profile [##UREF##1##18##]</td><td align="center">60</td><td align="center">80</td><td align="left">Affymetrix GeneChips on 22,000 transcripts</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Classification concordance of the breast cancer gene signatures (Kappa coefficient = 0.67)</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center"><bold>Intrinsic Subtype</bold></td><td align="center"><bold>No. of Patients</bold></td><td align="center" colspan="2"><bold>70-Gene Profile</bold></td><td align="center" colspan="2"><bold>Wound Response</bold></td><td align="center" colspan="2"><bold>Meta90</bold></td></tr></thead><tbody><tr><td/><td align="center">Patients</td><td align="center">Classification</td><td align="center">No. of Patients</td><td align="center">Classification</td><td align="center">No. of Patients</td><td align="center">Classification</td><td align="center">No. of Patients</td></tr><tr><td colspan="8"><hr/></td></tr><tr><td align="center">Basal-like</td><td align="center">36</td><td align="center">Good</td><td align="center">0</td><td align="center">Quiescent</td><td align="center">0</td><td align="center">Low</td><td align="center">0</td></tr><tr><td/><td/><td align="center">Poor</td><td align="center">36</td><td align="center">Activated</td><td align="center">36</td><td align="center">High</td><td align="center">36</td></tr><tr><td align="center">Luminal A</td><td align="center">91</td><td align="center">Good</td><td align="center">69</td><td align="center">Quiescent</td><td align="center">34</td><td align="center">Low</td><td align="center">89</td></tr><tr><td/><td/><td align="center">Poor</td><td align="center">22</td><td align="center">Activated</td><td align="center">57</td><td align="center">High</td><td align="center">2</td></tr><tr><td align="center">Luminal B</td><td align="center">41</td><td align="center">Good</td><td align="center">5</td><td align="center">Quiescent</td><td align="center">1</td><td align="center">Low</td><td align="center">16</td></tr><tr><td/><td/><td align="center">Poor</td><td align="center">36</td><td align="center">Activated</td><td align="center">40</td><td align="center">High</td><td align="center">25</td></tr><tr><td align="center">HER2+ and ER-</td><td align="center">28</td><td align="center">Good</td><td align="center">3</td><td align="center">Quiescent</td><td align="center">0</td><td align="center">Low</td><td align="center">8</td></tr><tr><td/><td/><td align="center">Poor</td><td align="center">25</td><td align="center">Activated</td><td align="center">28</td><td align="center">High</td><td align="center">20</td></tr><tr><td align="center">Normal-like</td><td align="center">23</td><td align="center">Good</td><td align="center">12</td><td align="center">Quiescent</td><td align="center">12</td><td align="center">Low</td><td align="center">22</td></tr><tr><td/><td/><td align="center">Poor</td><td align="center">11</td><td align="center">Activated</td><td align="center">11</td><td align="center">High</td><td align="center">1</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Eighteen common literature-defined oncogenic concepts (MCM) across the four breast cancer gene signatures (significance of overlap, P = 0.0002)</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="left"><bold>70-gene profile</bold></td><td align="left"><bold>Wound Response</bold></td><td align="left"><bold>Intrinsic Subtype</bold></td><td align="left"><bold>Meta-90</bold></td><td/><td/></tr></thead><tbody><tr><td align="left"><bold>Common GeneSet</bold></td><td align="center" colspan="4"><bold>No. of mapped genes (Enrichment p-value**)</bold></td><td align="center"><bold>MCM size</bold></td><td align="center"><bold>Description</bold></td></tr><tr><td colspan="7"><hr/></td></tr><tr><td align="left">MCM258</td><td align="left">7 (3e-04)</td><td align="left">25 (2e-04)</td><td align="left">19 (0.21)</td><td align="left">6 (0.09)</td><td align="left">350</td><td align="left">Downregulated genes in prostate cancer cells in response to MSA (full list)</td></tr><tr><td align="left">MCM458</td><td align="left">2 (0.16)</td><td align="left">24 (1e-04)</td><td align="left">34 (0.005)</td><td align="left">6 (0.24)</td><td align="left">322</td><td align="left">Differentially expressed genes in MCF7 cells expressing Myb</td></tr><tr><td align="left">MCM396</td><td align="left">10 (0.13)</td><td align="left">89 (2e-04)</td><td align="left">117 (0.005)</td><td align="left">20 (0.21)</td><td align="left">2265</td><td align="left">Upregulated genes in U937 cells expressing the PLZF/RAR fusion protein</td></tr><tr><td align="left">MCM691</td><td align="left">1 (0.1)</td><td align="left">6 (0.06)</td><td align="left">19 (6e-04)</td><td align="left">3 (0.16)</td><td align="left">101</td><td align="left">Up-regulated genes in untreated or permanently tamoxifen-treated MaCa 3366/TAM compared with MaCa 3366</td></tr><tr><td align="left">MCM513</td><td align="left">2 (0.22)</td><td align="left">24 (6e-04)</td><td align="left">29 (0.13)</td><td align="left">12 (0.04)</td><td align="left">375</td><td align="left">Differentially expressed genes in EcR-293 cells in response to BRCA1 expression</td></tr><tr><td align="left">MCM277</td><td align="left">1 (0.01)</td><td align="left">3 (0.02)</td><td align="left">5 (0.05)</td><td align="left">1 (0.16)</td><td align="left">22</td><td align="left">Upregulated genes in NCCIT cells in response to Wnt-3A</td></tr><tr><td align="left">MCM30</td><td align="left">1 (0.01)</td><td align="left">4 (0.004)</td><td align="left">3 (0.25)</td><td align="left">1 (0.15)</td><td align="left">24</td><td align="left">Upregulated genes in colorectal cancer cells</td></tr><tr><td align="left">MCM673</td><td align="left">2 (0.01)</td><td align="left">7 (0.008)</td><td align="left">7 (0.28)</td><td align="left">3 (0.15)</td><td align="left">79</td><td align="left">Androgen</td></tr><tr><td align="left">MCM6209872</td><td align="left">2 (0.003)</td><td align="left">2 (0.13)</td><td align="left">5 (0.09)</td><td align="left">1 (0.16)</td><td align="left">34</td><td align="left">Skin</td></tr><tr><td align="left">MCM349</td><td align="left">1 (0.01)</td><td align="left">2 (0.05)</td><td align="left">2 (0.25)</td><td align="left">2 (0.04)</td><td align="left">23</td><td align="left">Downregulated genes in hSNF5/INI1-deficient malignant rhabdoid tumor cell line upon hSNF5/INI1 expression</td></tr><tr><td align="left">MCM12</td><td align="left">2 (0.008)</td><td align="left">2 (0.25)</td><td align="left">7 (0.09)</td><td align="left">2 (0.14)</td><td align="left">56</td><td align="left">Aniogenic and Non-angiogenic tumours Signature</td></tr><tr><td align="left">MCM363</td><td align="left">5 (0.13)</td><td align="left">37 (0.007)</td><td align="left">46 (0.28)</td><td align="left">14 (0.15)</td><td align="left">808</td><td align="left">Upregulated genes in monocytes in response to IL-10 stimulation for 1 and 4 hours</td></tr><tr><td align="left">MCM574</td><td align="left">4 (0.04)</td><td align="left">22 (0.03)</td><td align="left">23 (0.27)</td><td align="left">6 (0.13)</td><td align="left">497</td><td align="left">Upregulated genes in advanced papillary serous tumor specimens</td></tr><tr><td align="left">MCM683</td><td align="left">1 (0.1)</td><td align="left">6 (0.06)</td><td align="left">11 (0.06)</td><td align="left">2 (0.17)</td><td align="left">111</td><td align="left">Downregulated genes wrt 3,5-diaryl-1,2,4-oxadiazole (MX-126374)</td></tr><tr><td align="left">MCM379</td><td align="left">1 (0.13)</td><td align="left">7 (0.07)</td><td align="left">12 (0.12)</td><td align="left">2 (0.29)</td><td align="left">129</td><td align="left">Unique genes regulated by tamoxifen, but not estradiol in osteosarcoma cells</td></tr><tr><td align="left">MCM1067</td><td align="left">1 (0.05)</td><td align="left">3 (0.15)</td><td align="left">5 (0.27)</td><td align="left">1 (0.21)</td><td align="left">64</td><td align="left">Upregulated genes in immmortilized epithelial cells in respense to Ad5-GFP infection</td></tr><tr><td align="left">MCM375</td><td align="left">1 (0.13)</td><td align="left">6 (0.12)</td><td align="left">10 (0.14)</td><td align="left">2 (0.27)</td><td align="left">127</td><td align="left">Unique genes regulated by estradiol, but not raloxifene in osteosarcoma cells</td></tr><tr><td align="left">MCM402</td><td align="left">1 (0.11)</td><td align="left">4 (0.25)</td><td align="left">8 (0.28)</td><td align="left">3 (0.14)</td><td align="left">116</td><td align="left">Downregulated genes in HepG2 T1 treated cells resulting from MIZ depletion</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Five common pathway sets (MsigDB) across the four breast cancer gene signatures (significance of overlap, <italic>P </italic>= 0.04).</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center"><bold>70 Gene</bold></td><td align="center"><bold>Wound Response</bold></td><td align="center"><bold>Intrinsic Subtype</bold></td><td align="center"><bold>Meta90</bold></td><td/></tr></thead><tbody><tr><td align="left"><bold>Common GeneSet</bold></td><td align="center" colspan="4"><bold>No. of mapped genes (Enrichment p-value**)</bold></td><td align="left"><bold>Description</bold></td></tr><tr><td colspan="6"><hr/></td></tr><tr><td align="left">breast cancer estrogen signaling</td><td align="center">1 (0.07)</td><td align="center">11 (0.001)</td><td align="center">11 (0.22)</td><td align="center">4 (0.11)</td><td align="left">GEArray</td></tr><tr><td align="left">EMT DOWN</td><td align="center">1 (0.02)</td><td align="center">2 (0.19)</td><td align="center">4 (0.29)</td><td align="center">1 (0.21)</td><td align="left">Jechlinger et al 2003</td></tr><tr><td align="left">CR DNA MET AND MOD</td><td align="center">1 (0.01)</td><td align="center">1 (0.24)</td><td align="center">3 (0.24)</td><td align="center">1 (0.12)</td><td align="left">PNAS 2007</td></tr><tr><td align="left">SA REG CASCADE OF CYCLIN EXPR</td><td align="center">1 (0.006)</td><td align="center">1 (0.12)</td><td align="center">2 (0.23)</td><td align="center">1 (0.07)</td><td align="left">SigmaAldrich</td></tr><tr><td align="left">reckPathway</td><td align="center">1 (0.005)</td><td align="center">1 (0.09)</td><td align="center">1 (0.28)</td><td align="center">1 (0.09)</td><td align="left">BioCarta</td></tr></tbody></table></table-wrap>
<inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1" name="1755-8794-1-28-i1" overflow="scroll"><mml:semantics><mml:mrow><mml:mo>−</mml:mo><mml:mn>2</mml:mn><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>K</mml:mi></mml:munderover><mml:mrow><mml:mi>log</mml:mi><mml:mo>⁡</mml:mo><mml:msub><mml:mi>P</mml:mi><mml:mi>i</mml:mi></mml:msub></mml:mrow></mml:mstyle></mml:mrow></mml:semantics></mml:math></inline-formula> <disp-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2" name="1755-8794-1-28-i2" overflow="scroll"><mml:semantics><mml:mrow><mml:msub><mml:mi>m</mml:mi><mml:mi>i</mml:mi></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:mrow><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>J</mml:mi></mml:munderover><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mrow></mml:mstyle></mml:mrow></mml:semantics></mml:math></disp-formula> <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M3" name="1755-8794-1-28-i3" overflow="scroll"><mml:semantics><mml:mrow><mml:mi>Y</mml:mi><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∑</mml:mo><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>J</mml:mi></mml:munderover><mml:mrow><mml:mstyle displaystyle="true"><mml:munderover><mml:mo>∏</mml:mo><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mi>K</mml:mi></mml:munderover><mml:mrow><mml:msub><mml:mi>X</mml:mi><mml:mrow><mml:mi>i</mml:mi><mml:mi>j</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:mstyle></mml:mrow></mml:mstyle></mml:mrow></mml:semantics></mml:math></inline-formula>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>List of modules genes involved in A. estrogen signaling and B. response to MSA in androgen-dependent prostate cell lines.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>List of the fifty-two common MCM sets shared between the two ER+ gene-signatures.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>List of the five common metabolic and signaling pathway sets (MsigDB) shared between the two ER+ gene-signatures.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Description of algorithm used to test for significance of overlap of datasets.</p></caption></supplementary-material>
<table-wrap-foot><p>**Hypergeometic test enrichment P-value adjusted for multiple testing.</p></table-wrap-foot> <table-wrap-foot><p>** Hypergeometic test enrichment P-value adjusted for multiple testing.</p></table-wrap-foot>
<graphic xlink:href="1755-8794-1-28-1"/>
<media xlink:href="1755-8794-1-28-S1.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1755-8794-1-28-S2.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1755-8794-1-28-S3.doc" mimetype="application" mime-subtype="msword"><caption><p>Click here for file</p></caption></media> <media xlink:href="1755-8794-1-28-S4.rtf" mimetype="application" mime-subtype="rtf"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Chang", "Sneddon", "Alizadeh", "Sood", "West", "Montgomery", "Chi", "Rijn", "Botstein", "Brown"], "given-names": ["HY", "JB", "AA", "R", "RB", "K", "JT", "M van de", "D", "PO"], "article-title": ["Gene expression signature of fibroblast serum response predicts human cancer progression: Similarities between tumors and wounds"], "source": ["Plos Biology"], "year": ["2004"], "volume": ["2"], "fpage": ["206"], "lpage": ["214"], "pub-id": ["10.1371/journal.pbio.0020007"]}, {"surname": ["Wang", "Atkins", "Zhang", "Yang", "Jatkoe", "Talantov", "Sieuwerts", "Timmermans", "Berns", "Klijn"], "given-names": ["Y", "D", "Y", "F", "T", "D", "A", "M", "E", "J"], "article-title": ["Gene expression profiles and molecular markers to predict distant metastasis of early stage breast cancers"], "source": ["Breast Cancer Research and Treatment"], "year": ["2003"], "volume": ["82"], "fpage": ["S120"], "lpage": ["S120"]}, {"article-title": ["BioCarta"]}, {"article-title": ["Signalling pathway database"]}, {"article-title": ["Kyoto Encyclopedia of Genes and Genomes"]}, {"surname": ["Benjamini", "Hochberg"], "given-names": ["Y", "Y"], "article-title": ["Controlling the False Discovery Rate \u2013 a Practical and Powerful Approach to Multiple Testing"], "source": ["Journal of the Royal Statistical Society Series B-Methodological"], "year": ["1995"], "volume": ["57"], "fpage": ["289"], "lpage": ["300"]}]
BMC Med Genomics. 2008 Jun 27; 1:28