- Research article
- Open Access
Identification of S100A8-correlated genes for prediction of disease progression in non-muscle invasive bladder cancer
BMC Cancer volume 10, Article number: 21 (2010)
S100 calcium binding protein A8 (S100A8) has been implicated as a prognostic indicator in several types of cancer. However, previous studies are limited in their ability to predict the clinical behavior of the cancer. Here, we sought to identify a molecular signature based on S100A8 expression and to assess its usefulness as a prognostic indicator of disease progression in non-muscle invasive bladder cancer (NMIBC).
We used 103 primary NMIBC specimens for microarray gene expression profiling. The median follow-up period for all patients was 57.6 months (range: 3.2 to 137.0 months). Various statistical methods, including the leave-one-out cross validation method, were applied to identify a gene expression signature able to predict the likelihood of progression. The prognostic value of the gene expression signature was validated in an independent cohort (n = 302).
Kaplan-Meier estimates revealed significant differences in disease progression associated with the expression signature of S100A8-correlated genes (log-rank test, P < 0.001). Multivariate Cox regression analysis revealed that the expression signature of S100A8-correlated genes was a strong predictor of disease progression (hazard ratio = 15.225, 95% confidence interval = 1.746 to 133.52, P = 0.014). We validated our results in an independent cohort and confirmed that this signature produced consistent prediction patterns. Finally, gene network analyses of the signature revealed that S100A8, IL1B, and S100A9 could be important mediators of the progression of NMIBC.
The prognostic molecular signature defined by S100A8-correlated genes represents a promising diagnostic tool for the identification of NMIBC patients that have a high risk of progression to muscle invasive bladder cancer.
Non-muscle invasive bladder cancer (NMIBC) is the most common histological subtype of bladder cancer, accounting for approximately 80% of all cases. Approximately 20% of these patients experience disease progression to muscle invasive bladder cancer (MIBC) after treatment, a development that is associated with a very poor prognosis for survival. Conventional histopathological parameters, such as tumor stage or grade, are generally considered to be prognostic factors, and numerous biomarkers have been investigated as prognostic indicators of the likelihood that NMIBC will develop into MIBC [1–6].
Members of the S100 family of calcium-binding proteins play essential roles in epithelial tissues and participate in a wide range of cellular processes, including transcription, proliferation, and differentiation [7–11]. At least 16 genes that encode members of the S100 family, including the gene for S100A8, are clustered on human chromosome 1q21 [12, 13], in a region that frequently experiences chromosomal rearrangement during tumor development [14, 15]. S100A8 is reportedly up-regulated in many cancers, including bladder cancer [16–23], and has been implicated in the regulation of tumor cell proliferation and metastasis [16, 24–26].
Although numerous diagnostic markers have been investigated as indicators of the risk of disease progression [1–6], none are able to sufficiently predict the behavior of NMIBC [1–6]. S100A8 has been suggested to be a predictive biomarker of bladder cancer outcome in several studies [21–23], however, the regulation of S100A8 gene expression and whether the genes associated with its expression provides additional insight into the mechanisms of disease progression or tumor invasion have not been studied. Therefore, we analyzed the expression pattern of S100A8 and its correlated genes to assess whether their molecular signature could identify patients with a higher likelihood of disease progression.
Patients and tissue samples
Primary NMIBC tissue samples from 103 consecutive cases of patients with histologically diagnosed transitional cell carcinoma were obtained from Chungbuk National University Hospital. To reduce confounding factors for affecting the analyses, any patients diagnosed with a concomitant carcinoma in situ (CIS) lesion or only CIS lesion were excluded. All tumors were macro-dissected, typically within 15 minutes of surgical resection. Each bladder cancer specimen was confirmed as representative by analysis of adjacent tissue in fresh frozen sections from transurethral resection (TUR) specimens, and then frozen in liquid nitrogen and stored at -80°C until use. The collection and analysis of all samples was approved by the Institutional Review Board of Chungbuk National University, and informed consent was obtained from each subject.
Tumors were staged and graded according to the 2002 TNM classification and the 2004 WHO grading system, respectively . A second TUR was performed 2-4 week after the initial resection when it was incomplete or when a high-grade or T1 tumor was detected . Patients with intermediate- or high-risk NMIBC received one cycle of intravesical BCG [27, 28]. All patients were followed and managed according to the standard recommendation for treatment of NMIBC [27–29]. In this study, we defined progression of the disease as an increase in stage from either Ta or T1 to T2 or higher after disease relapse .
RNA extraction, microarray experiments, and data processing
Total RNA was isolated by TRIzol reagent (Life Technologies, NY), according to the manufacturer's protocol. The quality and integrity of the RNA were confirmed by agarose gel electrophoresis and ethidium bromide staining, followed by visual examination under ultraviolet light. Five-hundred nanograms of total RNA were used for labeling hybridization according to the manufacturer's protocols (Illumina HumanWG-6 BeadChip, version 2). Arrays were scanned with an Illumina Bead Array Reader confocal scanner (BeadStation 500GXDW; Illumina, Inc., San Diego, CA) according to the manufacturer's instructions. After scanning, the microarray data were normalized using quantile normalization in the R language environment (version 2.8.1, available at http://www.r-project.org/). Measured gene expression values were log2 transformed and median centered across genes and samples. The full microarray data set is available in the NCBI Gene Expression Omnibus public database under the data series accession number GSE13507.
To classify patients into two groups, we used the median gene expression value of S100A8 as the cut-off. Pearson correlation coefficients were calculated to evaluate the association between S100A8 and its correlated genes. A hierarchical clustering algorithm, using the uncentered correlation coefficient as the measure of similarity and average linkage clustering, was applied as described in Eisen et al . The Kaplan-Meier method was used to calculate the time to progression, and differences between the times was assessed using log rank statistics. The prognostic value of the S100A8-correlated gene signature was determined with multivariate Cox proportional hazard regression models.
To validate the expression signature of S100A8 and its correlated genes identified in our original cohort, we applied it to independent microarray data from 302 patients with NMIBC, reported by Dyrskjot et al . For validation of the prognostic value of the molecular signature, we developed prediction models using the compound covariate predictor , Bayesian compound covariate predictor , linear discriminator analysis , nearest centroid classification , and support vector machines . The models incorporated genes that were differentially expressed between the two classes using a two-sample t-test. Genes were considered to have statistically significant differences in expression if the P-value was less than 0.001. We estimated the prediction error of each model using leave-one-out cross-validation (LOOCV), as described by Simon et al . For each LOOCV training set, the entire model-building procedure was repeated, including the gene selection process. Validation procedure was performed in BRB ArrayTools (version 3.7.1).
To explore the relationships between S100A8-correlated genes, we examined functional associations among the genes and generated gene networks based on whether they had more interconnected genes than would be expected to occur by chance. The significance of each network was estimated using the scoring system provided by the Ingenuity Pathway Analysis Tool (version 7.5). The scores were determined by the number of differentially expressed genes within each of the networks and the strength of the associations among the network members.
Table 1 details the baseline characteristics of the 103 primary NMIBC patients. The median age was 66.0 years (range: 24 to 88 years), and the median follow-up period after surgery was 57.6 months (range: 3.2 to 137.0 months). During the follow-up period, 11 of the 103 (10.7%) experienced disease progression.
Prognostic utility of S100A8-correlated genes
We analyzed 103 NMIBC samples and divided them into two groups according to the expression level of S100A8. The frequency of progression was significantly higher in the group with S100A8 expression levels in the upper 50th percentile than in the group with S100A8 expression levels in the lower 50th percentile (log-rank test, P = 0.003; Figure 1-A and 1-B).
We next sought to identify a gene expression signature that directly correlated with S100A8 expression levels during disease progression and then use that signature to predict the likelihood of tumor progression. We identified 1,015 genes whose change in expression correlated with S100A8 expression (Pearson correlation test, P < 0.001, r < -0.3 or r > 0.3). Based on hierarchical clustering analysis of the expression patterns of these genes, we divided the NMIBC samples into two groups: high S100A8 cluster (HSC) and low S100A8 cluster (LSC). The progression rate of HSC patients was significantly higher than that of the LSC patients (log-rank test, P < 0.001; Figure 2).
To evaluate the prognostic efficacy of the newly identified signature, we applied multivariate Cox regression analysis to the signature and known clinical and pathologic prognostic factors for NMIBC (Table 2). This analysis revealed that the molecular signature of S100A8-correlated genes (hazard ratio 15.225, 95% confidence interval = 1.746 to 133.52, P = 0.014) was the only strong predicator of bladder cancer progression.
Validation of the signature in an independent cohort
We next sought to validate our findings by using gene expression data from an independent cohort of European patients with bladder cancer . Patients with NMIBC from this cohort (n = 302) were divided into two groups according to the expression level of S100A8. Consistent with our results, progression of NMIBC was significantly higher in patients with S100A8 expression levels in the upper 50th percentile than in those with S100A8 expression levels in the lower 50th percentile (log-rank test, P < 0.001; Figure 1-C).
We also validated the association of the S100A8-correlated gene signature with disease progression in the European cohort. To overcome the peculiarities of any one particular prediction algorithm, we applied five different statistical methods to test the accuracy of our signature-based prediction of disease progression (Figure 3). We identified the genes with the greatest difference in expression level between the HSC and LSC subgroups in our original cohort (the training set). These genes were pooled to form a series of classifiers able to estimate the probability that a particular bladder cancer sample belonged to the HSC or LSC subgroup. The number of genes in the classifier set was optimized to minimize misclassification during LOOCV of the tumors in the training set. The performance of each prediction model is illustrated in Table 3. When applied to the European cohort (the test set), all five models produced consistent prediction patterns. Kaplan-Meier estimations in the test set revealed significant differences in the risk of progression between patients in subgroups HSC and LSC (Figure 3-B).
Biological insights into the signature for disease progression
To identify the predominant signalling networks active in the advancement from NMIBC to MIBC, gene network analysis of the 1,015 genes featured in the progression signature (Figure 2) was carried out using Ingenuity™ Pathways Analysis software. Of the 1,015 genes, 768 were mapped to gene networks defined by this tool. This analysis revealed a series of putative networks and associated functional categories. The 10 putative networks with highest scores are listed in Additional file 1 and their associated functions are illustrated in Figure 4.
As expected, genes involved in cellular growth and proliferation, cell cycle, cancer, and cell death were enriched, providing confidence in our results. We also found that genes involved in cell-mediated immune responses, immune cell trafficking, humoral immune responses, inflammatory responses, and immunological disease were significantly enriched (Figure 4).
Interestingly, among the top 10 putative networks, we observed functional connectivity between S100A8 and IL1B in the third network, in which S100A8 is a downstream effector of IL1B  (Figure 5). The expression level of IL1B was significantly higher in the HSC than in the LSC subgroup (two-sample t-test, P < 0.001; Figure 6-B), indicating its potential role in the progression of NMIBC to MIBC. We also found functional relationship between S100A8 and S100A9, in which S100A9 is a direct binding partner of S100A8 (Figure 5). S100A9 was more highly expressed in the HSC than in the LSC subgroup (two-sample t-test, P < 0.001; Figure 6-C).
The data presented in this study demonstrates that the expression signature of S100A8-correlated genes is able to predict the likelihood of bladder cancer progression. The validity of this signature as a prognostic indicator was confirmed by analysis of 302 cancers from an independent European cohort. In addition, based on the results of gene network analysis, we have identified a putative mechanism that may be responsible for disease progression.
Although considerable effort has been devoted to the establishment of a prognostic model of NMIBC that can provide information concerning survival and treatment options at diagnosis [1–6], the ability to predict the course of disease progression for patients with bladder cancer remains a major clinical challenge. A number of the molecular markers that have been identified to date have been explored as potential predictors of disease progression. Some of these, such as p53, have been suggested to be independent markers, while others do not appear to be effective as prognostic indicators [4–6]. Thus, there is a crucial need for methods capable of identifying patients with NMIBC that is likely to develop into MIBC. In the present study, we developed a method to predict the progression of primary NMIBC based on a gene expression signature. We showed that our method has strong predictive value through multivariate regression analysis and a validation study in an independent cohort. These data underscore the effectiveness of this molecular signature as a prognostic indicator in NMIBC, and suggest that this signature could be clinically useful.
The accumulation of recent reports that document the deregulated expression of calcium-binding proteins of the S100 family in a variety of human cancers implies that S100 proteins are required for neoplastic or metastatic transformation of tumor cells . Previous observations that S100A8 expression is enhanced in bladder cancer suggest a role for this protein in the progression of bladder cancer [21–23]. The data obtained from the current study demonstrated that upregulation of S100A8 expression was strongly associated with disease progression. This result substantiates the involvement of S100A8 during the progression of NMIBC.
Expression of S100A8 is not the only indicator of S100A8 activity, because it is regulated by many different mechanisms. The identification of stable and reliable human gene-to-gene relationships is an essential step towards unraveling the interactions and functional correlations between human genes . Therefore, we performed gene network analysis to identify the association of S100A8-correlated genes with progression of NMIBC. Strikingly, expression of both S100A8 and its correlated genes were strong predictors of the progression of bladder cancer (Figures 1 and 2). This finding was further supported by multivariate analysis, which showed that the molecular signature of S100A8-correlated genes was a strong predictor of cancer progression-free survival (Table 2). This result suggests that the S100A8-correlated gene signature retains its prognostic relevance even after additional pathological prognostic features have been taken into account.
The standard strategy for evaluating the accuracy of classification methods is application of a training-validation approach, in which a training set is utilized to identify the molecular signature and a validation set is used to estimate the degree of reliability. For the validation assay, a large sample size is needed before expression profiling can be utilized in a clinical setting . In the current study, the accuracy of a predictive S100A8-correlated gene signature was tested using five independent algorithms with gene expression data from a large independent European cohort (Figure 3). This result demonstrated not only a strong association between gene expression patterns and progression, but also provided strong evidence of the reliability of the prediction.
Based on an analysis of the S100A8-correlated genes in the context of gene networks, we identified a putative role for S100A8 in disease progression of bladder cancer. Because S100A8 is the downstream target of IL1B , our results strongly support the potential involvement of IL1B in tumor progression. Recent studies indicate that IL1B is associated with tumor invasiveness and metastasis [40, 41]. Many of the satellite genes associated with IL1B (i.e., CXCL2, PLA2G7, CCL11, S100A9, CD163, RNASE7, and OSM) (Figure 5) participate in inflammatory and immune responses, which are the best known activities of IL1B. Previous studies have demonstrated that IL1B activates inflammation that promotes tumor invasiveness . It was reported previously that S100A9 is up-regulated in conjunction with S100A8 in many cancers, including gastric cancer , prostate cancer [16, 17], and colorectal cancer [18, 19]. Both S100A8 and S100A9 have been implicated in the regulation of cell proliferation [16, 24] and metastatic processes . Thus, S100A9 may prove a valuable target for prevention of the migration of tumor cells to pre-metastatic sites .
Our results demonstrate that an expression signature consisting of S100A8 and its correlated genes can be a reliable prognostic indicator of progression in NMIBC, independent of traditional pathologic prognostic parameters. The use of this signature as a predictive indicator could potentially enable the prognosis of heterogeneous NMIBC patients to be determined at diagnosis, which would allow for individualized treatment and evaluation.
We describe a prognostic molecular signature based on gene expression of S100A8-correlated genes that can identify patients with high-risk NMIBC that is likely to develop into MIBC. Identification of patients with high-risk NMIBC may improve the effectiveness of currently available treatments and provide opportunities for the development of new treatment modalities.
(non-muscle invasive bladder cancer)
(muscle invasive bladder cancer)
(leave-one-out cross validation)
(high S100A8 cluster)
(low S100A8 cluster).
Pow-Sang JM, Seigne JD: Contemporary management of superficial bladder cancer. Cancer Control. 2000, 7 (4): 335-339.
Holmang S, Johansson SL: Stage Ta-T1 bladder cancer: the relationship between findings at first followup cystoscopy and subsequent recurrence and progression. J Urol. 2002, 167 (4): 1634-1637. 10.1016/S0022-5347(05)65168-3.
Herr HW: Tumor progression and survival of patients with high grade, noninvasive papillary (TaG3) bladder tumors: 15-year outcome. J Urol. 2000, 163 (1): 60-61. 10.1016/S0022-5347(05)67972-4.
Liou LS: Urothelial cancer biomarkers for detection and surveillance. Urology. 2006, 67 (3 Suppl 1): 25-33. 10.1016/j.urology.2006.01.034.
Cordon-Cardo C, Zhang ZF, Dalbagni G, Drobnjak M, Charytonowicz E, Hu SX, Xu HJ, Reuter VE, Benedict WF: Cooperative effects of p53 and pRB alterations in primary superficial bladder tumors. Cancer Res. 1997, 57 (7): 1217-1221.
Habuchi T, Marberger M, Droller MJ, Hemstreet GP, Grossman HB, Schalken JA, Schmitz-Drager BJ, Murphy WM, Bono AV, Goebell P, et al: Prognostic markers for bladder cancer: International Consensus Panel on bladder tumor markers. Urology. 2005, 66 (6 Suppl 1): 64-74. 10.1016/j.urology.2005.08.065.
Berridge MJ, Bootman MD, Roderick HL: Calcium signalling: dynamics, homeostasis and remodelling. Nat Rev Mol Cell Biol. 2003, 4 (7): 517-529. 10.1038/nrm1155.
Schafer BW, Heizmann CW: The S100 family of EF-hand calcium-binding proteins: functions and pathology. Trends Biochem Sci. 1996, 21 (4): 134-140.
Donato R: S100: a multigenic family of calcium-modulated proteins of the EF-hand type with intracellular and extracellular functional roles. Int J Biochem Cell Biol. 2001, 33 (7): 637-668. 10.1016/S1357-2725(01)00046-2.
Donato R: Intracellular and extracellular roles of S100 proteins. Microsc Res Tech. 2003, 60 (6): 540-551. 10.1002/jemt.10296.
Heizmann CW, Fritz G, Schafer BW: S100 proteins: structure, functions and pathology. Front Biosci. 2002, 7: d1356-1368. 10.2741/heizmann.
Marenholz I, Heizmann CW, Fritz G: S100 proteins in mouse and man: from evolution to function and pathology (including an update of the nomenclature). Biochem Biophys Res Commun. 2004, 322 (4): 1111-1122. 10.1016/j.bbrc.2004.07.096.
Ravasi T, Hsu K, Goyette J, Schroder K, Yang Z, Rahimi F, Miranda LP, Alewood PF, Hume DA, Geczy C: Probing the S100 protein family through genomic and functional analysis. Genomics. 2004, 84 (1): 10-22. 10.1016/j.ygeno.2004.02.002.
Schafer BW, Wicki R, Engelkamp D, Mattei MG, Heizmann CW: Isolation of a YAC clone covering a cluster of nine S100 genes on human chromosome 1q21: rationale for a new nomenclature of the S100 calcium-binding protein family. Genomics. 1995, 25 (3): 638-643. 10.1016/0888-7543(95)80005-7.
Schutte BC, Carpten JD, Forus A, Gregory SG, Horii A, White PS: Report and abstracts of the sixth international workshop on human chromosome 1 mapping 2000. Iowa City, Iowa, USA. 30 September-3 October 2000. Cytogenet Cell Genet. 2001, 92 (1-2): 23-41. 10.1159/000056867.
Hermani A, De Servi B, Medunjanin S, Tessier PA, Mayer D: S100A8 and S100A9 activate MAP kinase and NF-kappaB signaling pathways and trigger translocation of RAGE in human prostate cancer cells. Exp Cell Res. 2006, 312 (2): 184-197. 10.1016/j.yexcr.2005.10.013.
Hermani A, Hess J, De Servi B, Medunjanin S, Grobholz R, Trojan L, Angel P, Mayer D: Calcium-binding proteins S100A8 and S100A9 as novel diagnostic markers in human prostate cancer. Clin Cancer Res. 2005, 11 (14): 5146-5152. 10.1158/1078-0432.CCR-05-0352.
Stulik J, Koupilova K, Osterreicher J, Knizek J, Macela A, Bures J, Jandik P, Langr F, Dedic K, Jungblut PR: Protein abundance alterations in matched sets of macroscopically normal colon mucosa and colorectal carcinoma. Electrophoresis. 1999, 20 (18): 3638-3646. 10.1002/(SICI)1522-2683(19991201)20:18<3638::AID-ELPS3638>3.0.CO;2-W.
Stulik J, Osterreicher J, Koupilova K, Knizek , Macela A, Bures J, Jandik P, Langr F, Dedic K, Jungblut PR: The analysis of S100A9 and S100A8 expression in matched sets of macroscopically normal colon mucosa and colorectal carcinoma: the S100A9 and S100A8 positive cells underlie and invade tumor mass. Electrophoresis. 1999, 20 (4-5): 1047-1054. 10.1002/(SICI)1522-2683(19990101)20:4/5<1047::AID-ELPS1047>3.0.CO;2-E.
Yong HY, Moon A: Roles of calcium-binding proteins, S100A8 and S100A9, in invasive phenotype of human gastric cancer cells. Arch Pharm Res. 2007, 30 (1): 75-81. 10.1007/BF02977781.
Yao R, Lopez-Beltran A, Maclennan GT, Montironi R, Eble JN, Cheng L: Expression of S100 protein family members in the pathogenesis of bladder tumors. Anticancer Res. 2007, 27 (5A): 3051-3058.
Tolson JP, Flad T, Gnau V, Dihazi H, Hennenlotter J, Beck A, Mueller GA, Kuczyk M, Mueller CA: Differential detection of S100A8 in transitional cell carcinoma of the bladder by pair wise tissue proteomic and immunohistochemical analysis. Proteomics. 2006, 6 (2): 697-708. 10.1002/pmic.200500033.
Yao R, Davidson DD, Lopez-Beltran A, MacLennan GT, Montironi R, Cheng L: The S100 proteins for screening and prognostic grading of bladder cancer. Histol Histopathol. 2007, 22 (9): 1025-1032.
Cross SS, Hamdy FC, Deloulme JC, Rehman I: Expression of S100 proteins in normal human tissues and common cancers using tissue microarrays: S100A6, S100A8, S100A9 and S100A11 are all overexpressed in common cancers. Histopathology. 2005, 46 (3): 256-269. 10.1111/j.1365-2559.2005.02097.x.
Rafii S, Lyden D: S100 chemokines mediate bookmarking of premetastatic niches. Nat Cell Biol. 2006, 8 (12): 1321-1323. 10.1038/ncb1206-1321.
Salama I, Malone PS, Mihaimeed F, Jones JL: A review of the S100 proteins in cancer. Eur J Surg Oncol. 2008, 34 (4): 357-364.
Babjuk M, Oosterlinck W, Sylvester R, Kaasinen E, Bohle A, Palou-Redorta J: EAU guidelines on non-muscle-invasive urothelial carcinoma of the bladder. Eur Urol. 2008, 54 (2): 303-314. 10.1016/j.eururo.2008.04.051.
Hall MC, Chang SS, Dalbagni G, Pruthi RS, Seigne JD, Skinner EC, Wolf JS, Schellhammer PF: Guideline for the management of nonmuscle invasive bladder cancer (stages Ta, T1, and Tis): 2007 update. J Urol. 2007, 178 (6): 2314-2330. 10.1016/j.juro.2007.09.003.
Stenzl A, Cowan NC, De Santis M, Jakse G, Kuczyk MA, Merseburger AS, Ribal MJ, Sherif A, Witjes JA: The Updated EAU Guidelines on Muscle-Invasive and Metastatic Bladder Cancer. Eur Urol. 2009, 55 (4): 815-25. 10.1016/j.eururo.2009.01.002.
Malmstrom PU, Sylvester RJ, Crawford DE, Friedrich M, Krege S, Rintala E, Solsona E, Di Stasi SM, Witjes JA: An individual patient data meta-analysis of the long-term outcome of randomised studies comparing intravesical mitomycin C versus bacillus Calmette-Guerin for non-muscle-invasive bladder cancer. Eur Urol. 2009, 56 (2): 247-256. 10.1016/j.eururo.2009.04.038.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.
Dyrskjot L, Zieger K, Real FX, Malats N, Carrato A, Hurst C, Kotwal S, Knowles M, Malmstrom PU, de la Torre M, et al: Gene expression signatures predict outcome in non-muscle-invasive bladder carcinoma: a multicenter validation study. Clin Cancer Res. 2007, 13 (12): 3545-3551. 10.1158/1078-0432.CCR-06-2940.
Radmacher MD, McShane LM, Simon R: A paradigm for class prediction using gene expression profiles. J Comput Biol. 2002, 9 (3): 505-511. 10.1089/106652702760138592.
Dudoit S, Fridlyand J, Speed TP: Comparison of Discrimination Methods for the Classification of Tumors Using Gene Expression Data. J Am Stat Assoc. 2002, 97 (457): 10-10.1198/016214502753479248.
Ramaswamy S, Tamayo P, Rifkin R, Mukherjee S, Yeang CH, Angelo M, Ladd C, Reich M, Latulippe E, Mesirov JP, et al: Multiclass cancer diagnosis using tumor gene expression signatures. Proc Natl Acad Sci USA. 2001, 98 (26): 15149-15154. 10.1073/pnas.211566398.
Simon R, Radmacher MD, Dobbin K, McShane LM: Pitfalls in the use of DNA microarray data for diagnostic and prognostic classification. J Natl Cancer Inst. 2003, 95 (1): 14-18.
Rahimi F, Hsu K, Endoh Y, Geczy CL: FGF-2, IL-1beta and TGF-beta regulate fibroblast expression of S100A8. FEBS J. 2005, 272 (11): 2811-2827. 10.1111/j.1742-4658.2005.04703.x.
Prieto C, Risueno A, Fontanillo C, De las Rivas J: Human gene coexpression landscape: confident network derived from tissue transcriptomic profiles. PLoS One. 2008, 3 (12): e3911-10.1371/journal.pone.0003911.
Thorgeirsson SS, Lee JS, Grisham JW: Molecular prognostication of liver cancer: end of the beginning. J Hepatol. 2006, 44 (4): 798-805. 10.1016/j.jhep.2006.01.008.
Apte RN, Dotan S, Elkabets M, White MR, Reich E, Carmi Y, Song X, Dvozkin T, Krelin Y, Voronov E: The involvement of IL-1 in tumorigenesis, tumor invasiveness, metastasis and tumor-host interactions. Cancer Metastasis Rev. 2006, 25 (3): 387-408. 10.1007/s10555-006-9004-4.
Yano S, Nokihara H, Yamamoto A, Goto H, Ogawa H, Kanematsu T, Miki T, Uehara H, Saijo Y, Nukiwa T, et al: Multifunctional interleukin-1beta promotes metastasis of human lung cancer cells in SCID mice via enhanced expression of adhesion-, invasion- and angiogenesis-related molecules. Cancer Sci. 2003, 94 (3): 244-252. 10.1111/j.1349-7006.2003.tb01428.x.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/10/21/prepub
This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2009-0063255 and 2009-0068819).
The authors declare that they have no competing interests.
SKK conceived of the study, carried out the statistical analyses, and drafted the manuscript. EJK and SHL performed microarray preparation, sample selection, and RNA isolation. YSH investigated the clinical records of the considered patients and contributed to the interpretation of the results. YJK performed the statistical analyses, contributed to the interpretation of the results, and helped to draft the manuscript. WJK designed the study concept, interpreted the results and approved the final manuscript. All authors have read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
About this article
Cite this article
Kim, SK., Kim, EJ., Leem, SH. et al. Identification of S100A8-correlated genes for prediction of disease progression in non-muscle invasive bladder cancer. BMC Cancer 10, 21 (2010). https://doi.org/10.1186/1471-2407-10-21
- Bladder Cancer
- Gene Expression Signature
- Invasive Bladder Cancer
- Muscle Invasive Bladder Cancer
- S100 Family