Skip to main content
  • Research article
  • Open access
  • Published:

Prognostic values of S100 family members in ovarian cancer patients

Abstract

Objective

Exhibiting high consistence in sequence and structure, S100 family members are interchangeable in function and they show a wide spectrum of biological processes, including proliferation, apoptosis, migration, inflammation and differentiation and the like. While the prognostic value of each individual S100 in ovarian cancer is still elusive. In current study, we investigated the prognostic value of S100 family members in the ovarian cancer.

Methods

We used the Kaplan Meier plotter (KM plotter) database, in which updated gene expression data and survival information are from 1657 ovarian cancer patients, to assess the relevance of individual S100 family mRNA expression to overall survival in various ovarian cancer subtypes and different clinicopathological features.

Results

It was found that high expression of S100A2 (HR = 1.18, 95%CI: 1.04–1.34, P = 0.012), S100A7A (HR = 1.3, 95%CI: 1.04–1.63, P = 0.02),S100A10 (HR = 1.2, 95%CI: 1.05–1.38, P = 0.0087),and S100A16 (HR = 1.23, 95%CI: 1–1.51, P = 0.052) were significantly correlated with worse OS in all ovarian cancer patients, while the expression of S100A1 (HR = 0.87, 95%CI: 0.77–0.99, P = 0.039), S100A3 (HR = 0.83, 95%CI: 0.71–0.96, P = 0.0011), S100A5 (HR = 0.84, 95%CI: 0.73–0.97, P = 0.017), S100A6 (HR = 0.84, 95%CI: 0.72–0.98, P = 0.024), S100A13 (HR = 0.85, 95%CI:0.75–0.97, P = 0.014) and S100G (HR = 0.86, 95%CI: 0.74–0.99, P = 0.041) were associated with better prognosis. Furthermore, we assessed the prognostic value of S100 expression in different subtypes and the clinicopathological features, including pathological grades, clinical stages and TP53 mutation status, of ovarian cancer patients.

Conclusion

Comprehensive understanding of the S100 family members may have guiding significance for the diagnosis and outcome of ovarian cancer patients.

Peer Review reports

Background

As the leading cause of gynecologic cancer-related deaths worldwide, almost 75% ovarian cancer patients are diagnosed in late stage, losing the best time for operation. Thus, it is urgent to find new prognosis or therapeutic targets for ovarian cancer.

In humans, S100 family is consisted of 21 acidic-Ca2+ binding proteins, which are highly conserved and functioned as both intracellular Ca2+ sensors and extracellular factors [1]. Among 21 members of S100 family, four are found dispersed throughout the genome (including S100B-21q22, S100G-Xp22, S100P-4p16 and S100Z-5q14); and the rest ones (from S100A1 to A100A14, S100A7 and S100A16) are encoded in two series clusters within a 2 Mb region on chromosome locus 1q21 [1]. Exhibiting high consistence in sequence and structure, S100 family members are interchangeable in function and they show a wide spectrum of biological processes, including proliferation, apoptosis,migration,inflammation and differentiation and the like [2].

Dysregulated expression of multiple S100 proteins have been found to be involved in the progression of numerous tumors. In addition, the protein expression of S100 family members has been used to facilitate diagnosis, predict prognosis, determine treatment options and monitor therapeutic efficacy [2, 3].

The expression and function of S100 protein vary across different malignant tumors or tumor subtypes. For instance, in oral carcinoma, S100A2 acts as a tumor-suppressor, while tumor-promoter in lung carcinoma [4, 5]. In ovarian cancer, many S100 family members have been reported, such as S100A1, S100A4, S100A6, S100B and S100P. DeRycke MS’s report indicated that S100A1 mRNA and protein expression elevated with increasing Silverberg grade but not stage in serous ovarian cancer, while in endometrial subtype of ovarian cancer, there was no correlation between S100A1 expression and clinical stage or grade, but it was related with decreased relapse-free survival (RFS) [6]. Serum S100A6 concentration has been found that it predicts peritoneal tumor burden and is associated with advanced stage in ovarian cancer [7]. S100A4, S100B, S100A14, S100P expressions are associated with poor prognosis [8,9,10,11,12], moreover,high cytoplasmic S100A10 staining is significantly associated with reduced overall survival (OS) and progression-free survival (PFS) in serous ovarian cancer [13, 14]. Whereas, some other S100 family members, including S100A3, S100A5, S100A7, S100A7A, S100A8, S100A16 and S100G are less reported in ovarian carcinoma. The prognosis value of each S100 member, particularly at the transcriptional level in ovarian cancer patients needs further to be explored. In this report, we evaluated the prognosis value of each S100 member mRNA level in ovarian cancer patients by Kaplan Meier plotter platform.

Methods

We used an online database to assess the relevance of individual S100 mRNA expression to OS. The database was developed by use of gene-chip or RNA sequencing data with clinical information of the 1647 EOC (epithelial ovarian cancer) patients from GEO (Gene Expression Omnibus; www.ncbi.nlm.nih.gov/geo/) and TCGA (The Cancer Genome Atlas; http://cancergenome.nih.gov), which integrates the gene expression with clinical information. In order to examine the prognosis effect of an individual gene, the patients were divided into two groups in accordance with the expression level of this gene.

Briefly, the 21 S100 family members were individually filled into the Kaplan Meier plotter database to obtain the survival plots. HR (hazard ratio), 95 %CI (confidence intervals) and log rank P value were provided.

Differential expression data between normal and ovarian cancer tissues was from the GEPIA (Gene Expression Profiling Interactive Analysis) database (http://gepia.cancer-pku.cn/index.html), which is a newly developed interactive web server for analyzing the RNA sequencing expression data of 9736 tumors and 8587 normal samples from the TCGA and the GTEx (Genotype-Tissue Expression) projects [15].

Results

Prognostic values of S100 members in all ovarian cancer patients

Respectively, we evaluated the prognosis role of each S100 family member in http://kmplot.com/analysis/. Among all the 21 family members, 10 members were correlated with significant OS in all EOC patients. Survival curves were plotted for all patients, and were shown in Fig. 1. It was found that S100A2 (HR = 1.18, 95%CI: 1.04–1.34, P = 0.012), S100A7A (HR = 1.3, 95%CI: 1.04–1.63, P = 0.02), S100A10 (HR = 1.2, 95%CI: 1.05–1.38, P = 0.0087) and S100A16 (HR = 1.23, 95%CI: 1–1.51, P = 0.052) mRNA high expression was significantly correlated with worse OS in all ovarian cancer patients (Fig. 1). We also observed that the high expression of S100A1 (HR = 0.87, 95%CI: 0.77–0.99, P = 0.039), S100A3 (HR = 0.83, 95%CI: 0.71–0.96, P = 0.0011), S100A5 (HR = 0.84, 95%CI: 0.73–0.97, P = 0.017), S100A6 (HR = 0.84, 95%CI: 0.72–0.98, P = 0.024), S100A13 (HR = 0.85, 95%CI:0.75–0.97, P = 0.014), S100G (HR = 0.86, 95%CI: 0.74–0.99, P = 0.041) were associated with longer survival time (Fig. 1). In addition, the other members of the S100 family were found no correlations with the patient’s survival (Additional file 1: Table S1).

Fig. 1
figure 1

The prognostic value of S100A mRNA expression in all ovarian cancer patients in www.kmplot.com. Notes: a Prognostic HRs of individual S100 members in all ovarian cancer patients. b-k Survival curves of S100A2(the desired Affymetrix IDs is valid: 204268_at), S100A7A (Affymetrix IDs: 232170_at), S100A10(Affymetrix IDs: 200872_at), S100A3(Affymetrix IDs: 206027_at), S100A1(Affymetrix IDs: 205334_at), S100A5(Affymetrix IDs: 207763_at), S100A6(217728_at), S100A13(Affymetrix IDs: 202598_at) and S100G(207885_at) are plotted for all ovarian cancer patients. Abbreviation: HR: hazard ratio

Prognostic values of S100 members in different ovarian cancer subtypes

Next, we also examined the prognosis effects of S100 family members in two subtypes of ovarian cancer patients, including serous and endometrioid ovarian cancer. As shown in Fig. 2, for S100A1 (HR = 0.86, 95%CI: 0.74–1.0, P = 0.055), S100A5 (HR = 0.83, 95%CI: 0.7–0.98, P = 0.025), S100A6 (HR = 0.8, 95%CI: 0.67–0.95, P = 0.013), S100A8 (HR = 0.84, 95%CI: 0.71–0.99, P = 0.038) and S100A13 (HR = 0.8, 95%CI: 0.68–0.94, P = 0.0055), high mRNA expression predicted favorable OS in serous ovarian cancers respectively. However, S100A10 (HR = 1.25, 95%CI: 1.05–1.49, P = 0.011) predicted worse survival. The rest members of S100 family were not related with prognosis effect in serous ovarian cancer (Additional file 1: Table S2).

Fig. 2
figure 2

Prognostic values of S100 members in different ovarian cancer subtypes. Notes: a-f Survival curves of S100A1 (the desired A ffymetrix IDs is valid: 205334_at), S100A5(Affymetrix IDs: 207763_at), S100A6(Affymetrix IDs: 217728_at), S100A8 (Affymetrix IDs: 202917_s_at), S100A13(Affymetrix IDs: 202598_at) and S100A10 (Affymetrix IDs: 200872_at) are plotted for serous ovarian cancer patients. g-k Survival curves of S100A12 (the desired Affymetrix IDs is valid: 205863_at), S100A13(Affymetrix IDs: 202598_at), S100G (Affymetrix IDs: 207885_at), S100A2(Affymetrix IDs: 204268_at) and S100A10 (Affymetrix IDs: 200872_at) are plotted for endometrioid ovarian cancer patients. Abbreviation: HR: hazard ratio

In endometrioid ovarian cancer, higher mRNA expression of S100A12 (HR = 0.1, 95%CI: 0.01–0.89, P = 0.011), S100A13 (HR = 0.16, 95%CI: 0.03–0.97, P = 0.022) and S100G (HR = 0.1, 95%CI: 0.01–0.86, P = 0.0092) were associated with better OS, while S100A2 (HR = 4.9, 95%CI: 0.82–29.32, P = 0.054) and S100A10(HR =1,691,454,124.4, 95%CI: (0-Inf), P = 0.022) predicted worse survival (Fig. 2). We found no correlation between the rest S100 family members with endometrioid ovarian cancer patients’ survival. (Additional file 1: Table S3).

Prognostic values of S100 members in ovarian cancer patients with different clinicopathological features

Furthermore, we examined the relationship between S100 expression and the clinicopathological features in the ovarian cancer patients, including pathologic grades, clinical stages and TP53 mutation status (Additional file 2: Table S7). With less number of cases in grade I, we focused on the grade II and III EOC patients. As was shown in Table 1, higher S100A4 mRNA expression (HR = 1.46, 95%CI: 1.06–2.03, P = 0.021), S100A6 (HR = 1.4, 95%CI: 1.02–1.93, P = 0.038), S100A16 (HR = 1.75, 95%CI: 1.02–3.01, P = 0.04) and S100B (HR = 1.41, 95%CI: 1.02–1.94, P = 0.036) was found to be associated with significantly worse OS in grade II ovarian cancer patients respectively. While, the S100A12 (HR = 0.56, 95%CI: 0.41–0.78, P = 0.0043), S100A13 (HR = 0.73, 95%CI: 0.53–1.01, P = 0.058) and S100Z (HR = 0.6, 95%CI: 0.37–0.96, P = 0.032) expression predicted better prognosis in grade II ovarian cancer patients (Table 1). For grade III patients, the S100A2 (HR = 1.27, 95%CI: 1.07–1.5, P = 0.0051), S100A10 (HR = 1.38, 95%CI: 1.15–1.65, P = 0.00058), S100A16 (HR = 1.29, 95%CI: 1–1.67 P = 0.05) and S100P (HR = 1.25, 95%CI: 1.06–1.48, P = 0.0096) expression were found to be correlated with shorter overall survival time. While, the S100A3 (HR = 0.78, 95%CI: 0.65–0.93, P = 0.0056), S100A5 (HR = 0.79, 95%CI: 0.66–0.95, P = 0.011), S100A6 (HR = 0.83, 95%CI: 0.7–1, P = 0.051), S100A8 (HR = 0.75, 95%CI: 0.63–0.89, P = 0.0012) and S100B (HR = 0.83, 95%CI: 0.7–0.99, P = 0.037) expression were related to better prognosis respectively.

Table 1 Correlation of S100 with different pathological grade status of ovarian cancer patients

Then, we further investigated the relationship between the prognosis and the clinical stages. In order to facilitate the research,we combined the stage of phase I and phase II, and, we also combined phase III and phase IV. As was shown in Table 2, S100B (HR = 2.24, 95%CI: 1.02–4.92, P = 0.039) and S100P (HR = 2.38, 95%CI: 1.1–5.18, P = 0.024) were associated with worse survival, while, S100A7 (HR = 0.34, 95%CI: 0.13–0.85, P = 0.016), S100A13 (HR = 0.23, 95%CI: 0.08–0.66, P = 0.0028) and S100G (HR = 0.27, 95%CI: 0.12–0.61, P = 0.00076) were associated with better prognosis in stage I + II patients. In stage III + IV ovarian cancer patients, the expression of S100A10 (HR = 1.3, 95%CI: 1.11–1.51, P = 0.0011), S100A16 (HR = 1.37, 95%CI: 1.04–1.81, P = 0.026) and S100P (HR = 1.24, 95%CI: 1.06–1.46, P = 0.006) were correlated with worse OS and S100A11 (HR = 1.16, 95%CI: 1.0–1.35, P = 0.057) with modestly worse OS; while, for S100A1 (HR = 0.85, 95%CI:0.74–0.99, P = 0.036), S100A3 (HR = 0.8, 95%CI: 0.68–0.94, P = 0.008), S100A5 (HR = 0.81, 95%CI: 0.68–0.97, P = 0.022), S100A8 (HR = 0.84, 95%CI: 0.71–0.99, P = 0.036), S100A12 (HR = 0.82, 95%CI: 0.7–0.96, P = 0.012) and S100B (HR = 0.85, 95%CI: 0.74–0.99, P = 0.036), the mRNA expression levels of them predicted favorable overall survival. S100G (HR = 0.85, 95%CI: 0.73–1, P = 0.052) was only slightly correlated with favorable OS but with no statistic difference. Figure 3 has shown that mRNA expression of S100A8 (HR = 1.84, 95%CI: 1.05–3.23, P = 0.031), S100A11 (HR = 2.02, 95%CI: 1.06–3.85, P = 0.029) and S100B (HR = 2.04, 95%CI: 1.16–3.58, P = 0.012) were related with worse survival in wild-TP53-type ovarian cancer patients, while, S100A2 (HR = 0.45, 95%CI: 0.22–0.93, P = 0.027), S100A3 (HR =0.41, 95%CI: 0.23–0.74, P = 0.0023) and S100A5 (HR = 0.39, 95%CI: 0.18–0.86, P = 0.016) were related with better OS in wild-p53-type EOC patients. We also found that the level of S100A7A (HR = 1.71, 95%CI: 1.15–2.56, P = 0.0079), S100A12 (HR = 1.31, 95%CI: 1.04–1,64, P = 0.019), S100A14 (HR = 1.3, 95%CI: 1.02–1.65, P = 0.035), S100P (HR = 1.28, 95%CI: 1.01–1.62, P = 0.038) and S100G (HR = 1.27, 95%CI: 1.01–1.61, P = 0.043) predicted worse OS in mutant-p53-type EOC patients; however, S100A1 (HR = 0.66, 95%CI: 0.51–0.85, P = 0.0011), S100A13 (HR = 0.75, 95%CI: 0.6–0.94, P = 0.013) and S100A16 (HR = 0.64, 95%CI: 0.42–0.98, P = 0.08) were associated with better prognosis in ovarian cancer patients (Fig. 4).

Table 2 Correlation of S100 with different clinical stages of ovarian cancer patients
Fig. 3
figure 3

Prognostic values of S100 members in wild-TP53-type ovarian cancer patients. Notes: a-f Survival curves of S100A2 (the desired Affymetrix IDs is valid: 204268_at), S100A3(Affymetrix IDs: 206027_at), S100A5(Affymetrix IDs: 207763_at), S100A8 (Affymetrix IDs: 202917_s_at), S100A11(Affymetrix IDs: 200660_at) and S100B (Affymetrix IDs: 209686_at) are plotted for wild-TP53-type ovarian cancer patients. Abbreviation: HR: hazard ratio

Fig. 4
figure 4

Prognostic values of S100 members in mutant-TP53-type ovarian cancer patients. Notes: a-i Survival curves of S100A1 (the desired Affymetrix IDs is valid: 205334_at), S100A7A (Affymetrix IDs: 232170_at), S100A12(Affymetrix IDs: 205863_at), S100A13(Affymetrix IDs: 202598_at), S100A14(Affymetrix IDs: 218677_at), S100A16(Affymetrix IDs: 227998_at), S100B(Affymetrix IDs: 209686_at),S100G(Affymetrix IDs: 207885_at) and S100P (Affymetrix IDs: 204351_at) are plotted for mutant-TP53-type ovarian cancer patients. Abbreviation: HR: hazard ratio

Prognostic values of S100 members in different surgery and chemotherapy

Finally, we investigated the relationships between the S100 members’ prognostic values and different surgery, chemotherapy, which were shown in Additional file 1: Table S4 and S5.

Compared with suboptimal debulk, the high expression of S100A3 (HR = 0.61, 95%CI: 0.49–0.76, P = 1.10E-05), S100A7 (HR = 0.76, 95%CI: 0.61–0.94, P = 0.012), S100A9 (HR = 0.78, 95%CI: 0.63–0.96, P = 0.018), S100A12 (HR = 0.71, 95%CI: 0.57–0.9, P = 0.0034) and S100G (HR = 0.69, 95%CI: 0.55–0.86, P = 0.0011) were associated with better OS for patients with surgery of optimal debulk.

As for chemotherapy strategies containing platin, S100A5 (HR = 0.85, 95%CI: 0.72–0.99, P = 0.042), S100B (HR = 0.86, 95%CI: 0.74–0.99, P = 0.041) and S100G (HR = 0.84, 95%CI: 0.72–0.99, P = 0.035) high expression were associated with better OS. As for chemotherapy containing taxol. S100G (HR = 0.80, 95%CI: 0.65–0.97, P = 0.027) mRNA predicted better OS. As for chemotherapy containing platin and taxol, S100A12 (HR = 0.81, 95%CI: 0.66–1.00, P = 0.047) and S100G (HR = 0.81, 95%CI: 0.66–0.99, P = 0.039) were associated with better OS for EOC patients. With regard to chemotherapy containing docetaxel, S100A3 (HR = 0.52, 95%CI: 0.28–0.94, P = 0.027), S100A4 (HR = 0.55, 95%CI: 0.31–0.97, P = 0.036), S100A7 (HR = 0.40, 95%CI: 0.21–0.73, P = 0.0022), S100A8 (HR = 0.54, 95%CI: 0.31–0.92, P = 0.022), S100A12 (HR = 0.56, 95%CI: 0.32–1, P = 0.047) and S100G (HR = 0.49, 95%CI: 0.28–0.85, P = 0.0098) were associated with better OS.

In order to better understand the relationship between the S100 expression and prognosis, we learned the differential expression information of each S100 family member between normal and ovarian cancer tissues (Additional file 1: Table S6). And we found that there were obvious expression differences of S100A1, S100A2, S100A4, S100A9, S100A11, S100A13, S100A14 and S100P (P < 0.05). Interestingly, we also found that, compared to the other normal organs across the whole body, each S100 family member has extremely low expression in the normal ovarian tissue. However, in the tumor tissues, most members have a relative higher expression level in the ovarian cancer versus other organ tumor tissues (Additional file 3: Figure S2A and Additional file 4: Table S8). To some extent, the result above indicates the specificity of the S100 family members in ovarian cancer.

Discussion

Until now, a vast majority of previous researches has focused on overall survival of ovarian cancer patients, and transcriptome analysis has been used to hunt for molecular signatures correlated with prognosis of ovarian cancer patients. However, the small sample sizes in some researches raised doubts about creditability and generalizability of the results. In our work, we identified the relationship between the prognosis and transcriptional level of the S100 family members in a large scale of patients by use of Kaplan Meier plotter database.

Variation in S100A1 expression have been reported in various diseases, including malignant tumor, heart failure, diabetes, ischemia, and chronic pulmonary hypertension [16, 17]. It has been reported that S100A1 expression was significantly related with favorable OS in breast cancer patients [3], which is in accordance with all ovarian cancer patients, regardless of specific subtypes. However, DeRycke M. S et al. have reported that S100A1 expression was a prognostic indicator of relapse-free survival in endometrioid subtypes of ovarian cancer [6]. In addition, although we found no significant OS prognostic value of S100A1 expression, it signifies poor PFS (progressive free survival) for endometrioid ovarian cancer patients (Additional file 5: Figure S1). Thus, we suppose that S100A1 may play a more important role in disease relapse or progression in endometrioid subtypes of ovarian cancer patients. The specific mechanisms behind this situation need further investigation.

S100A2 has been identified as a markedly downregulated gene in tumor-derived mammary epithelial cell lines and was assumed as a tumor suppressor gene [18]. Meanwhile,a down-regulation of S100A2 can be verified in many human cancers and certified as a prognostic marker, such as melanoma, prostate and breast cancer [19,20,21,22]. On the other hand, some other type of tumors, including pancreas, colorectal, non-small-cell lung carcinoma (NSCLC) exhibit an upregulation of S100A2, which are associated with poor patient survival [5, 23, 24]. The role of S100A2 in ovarian cancer has rarely been explored. Hough, C. D., et al. have shown the expression of S100A2 is elevated in ovarian cancer [25]. In accordance, in our study, for all ovarian cancer patients, S100A2 mRNA expression was correlated with poor prognosis. However, for wild-TP53-type ovarian cancer patients, the higher level of S100A2 predicted better OS, but showed no prognostic value for mutant-p53-type ovarian cancer patients. Buckley NE has reported that S100A2 modulated binding of mutant TP53 to HSP90, and loss of S100A2 leaded to an HSP90-dependent stabilization of mutant TP53 [26]. It was also reported S100A2 protein co-localized with the tumor suppressor wild type TP53, and increased TP53 transcription activity through a Ca2+-dependent interaction [27]. We guess that the interaction between S100A2 and TP53 may mediate the tumor suppressor effect of S100A2 in wild-TP53-type ovarian cancer patients. As to its poor prognosis for all patients, Buckley NE also demonstrated that, as downstream of the BRCA1/ΔNp63 signaling axis, S100A2 played an important role in regulating transcriptional responses and enhancing growth control mechanisms through destabilization of mutant TP53 [26]. Thus, we speculated that, as the downstream of the BRCA1/ΔNp63 signaling axis, S100A2 might quite possibly be upregulated in modulating transcriptional responses in those patients with genomic instability, including BRCA1 or TP53 mutation, but not limited to those, who had worse survival than those without genomic instability.

S100A3 is a matricellular protein, which has been identified in various tissues, such as diaphragm, heart, liver, placenta, skeletal muscle and so on [28, 29]. S100A3 promotes tumor progression in a variety of tumors. In colorectal cancer, the high level of S100A3 were related to the tumor occurrence and progression [30]. In castration-resistant prostate cancer cells, S100A3 inhibition reduced invasion and tumor growth [31]. It was also reported that S100A3 was relatively highly expressed in poorly differentiated and advanced gastric cancer [28]. While, in breast cancer, the loss of S100A3 expression was associated with malignant development [29], which was in accordance with its favorable prognosis in ovarian cancer in our study. Up to now, the prognostic role of S100A3 in ovarian cancer has not been reported.

S100A5 is a novel member of S100 protein family which is characterized by its interaction with Ca2+, Zn2+, and Cu[2+.[32] It was detected that S100A5 expressed in very restricted region of adult brain [32, 33]. Hancq, S et al. showed an inverse relationship between the S100A5 expression and the risk of recurrence of totally resected WHO grade I meningiomas [34]. In consistent with this result, this study indicated that S100A5 predicted favorable OS in all EOC patients.

As a member of S100 protein family, S100A6 is highly expressed in epithelial cells, fibroblasts and in several types of tumor cells, which can regulate the cellular function of proliferation, apoptosis, cytoskeleton dynamics, tumorigenesis [35, 36]. In several human carcinomas, the expression of S100A6 appears to be associated with aggressive disease and high levels of S100A6 was an independent marker for poor prognostic, such as colorectal, gastric, thyroid, hepatocellular carcinoma, and so on [37,38,39,40]. Serum S100A6 concentration has been found that it predicts peritoneal tumor burden and is associated with advanced stage in ovarian cancer [7]. Nevertheless, in our study, high expression of S100A6 predicted favorable OS in ovarian cancer patients, especially in serous type ovarian cancer. Interestingly, we observed a positive relationship between S100A6 mRNA level and OS in stage II patients, but negative in stage IV patients. Thus, we suggest that S100A6 may play different roles in early and late stage ovarian cancer patients. It also reported that S100A6 could interact with p53 and affect its biological activity [41]. While, S100A6 expression did not show any correlation with overall survival in patients with wild or mutant TP53 status.

S100A8 and S100A9 originally existed as homodimers similar with many other S100 family members, and preferentially expressed in myeloid cells, such as granulocytes, monocytes, osteoclasts, and so on [42,43,44]. It has been longer recognized the relationship between the S100A8/S100A9 expression and inflammatory disorders, including rheumatoid arthritis, inflammatory bowel disease, multiple sclerosis [42, 45, 46]. As is well known, acute and chronic inflammation has been confirmed to increase the risk of tumorigenesis. Strong up-regulation of S100A8/S100A9 has also been detected in various human cancers and promotes tumors progression, including gastric, colon, breast, and live cancer [46,47,48,49,50]. While, several other studies have indicated that extracellular S100A8/S100A9 exhibits powerful anti-tumor responses by promoting cytotoxicity and apoptosis [51]. Until now, little was known about the relationship between the S100A8/S100A9 and prognosis in ovarian cancer. In our research, the higher expression of S100A8 predicted favorable OS in serous, as well as grade III and stage III + IV ovarian cancer patients. However, higher expression of S100A9 was not correlated with prognostic value in ovarian cancer patients.

S100A10 is known to form (annexinA2)2-(S100A10)2 heterotetramer (AIIt), which results in the translocation of S100A10 to the plasma membrane [52]. Extracellularly, the S100A10/annexin II complex functions as a plasminogen receptor, and regulates the tissue plasminogen activator (tPA)-dependent plasminogen activation and the plasmin formation on cancer cells. Intracellularly, S100A10 implicates in various cell processes, including cell growth, cell cycle, transcription, and cell differentiation [53].

Furthermore, these processes mentioned above were thought to be associated with tumor invasion, metastasis and angiogenesis [53, 54]. In addition, S100 plays a role in the regulation of cell proliferation by negatively affecting BAD-induced apoptosis. In colorectal cancer cells, S100A10 expression are associated with resistance to oxaliplatin (L-OHP) [55]. It has been also reported that S100A10 expression contributed to the aggressive characteristics of thyroid anaplastic carcinoma [56]. Recently, Lokman, N. A., et al. reported high S100A10 mRNA levels powerfully predicted poor outcome in serous ovarian cancer [13]. Here, our research also supported that the high level of S100A10 may indicated reduced OS in serous, grade III, or stage III + IV ovarian cancer patients.

S100A13 is a novel member of S100 family that characterized by its specificity in various form of cancer. Previously, S100A13 has been indicated to be a powerfully angiogenic biomarker for melanoma and astrocytic gliomas [57,58,59]. Recently, it has been reported that S100A13 expression contributes to more aggressive invasive phenotype in lung cancer cells [60]. S100A13 also influences chemoresistance by regulating the secretion of FGF1 and IL1A, or functions in several key signaling pathways such as cytokine and NF-kB signaling. Contrary to these functions, our research showed that upregulated S100A13 was associated better OS in all ovarian cancer, grade II, stage I + II, and mutant-p53-type ovarian cancer patients.

S100A16 was a recent addition to the S100 family that isolated from astrocytoma [61]. The expression of S100A16 was upregulated in various tumors, such as bladder, lung, thyroid gland, pancreas, and so on [62, 63]. In lung adenocarcinoma, S100A16 expression was associated with vessel invasion and poor outcome and confirmed as an independent prognostic marker [62]. Similarity,higher expression of S100A16 predicted worse OS in breast cancer by promoting epithelial-mesenchymal transition [64, 65]. However, high membrane S100A16 predicted better OS in colorectal cancer as well as oral squamous cell carcinoma [66, 67]. Our results indicated that higher level of S100A16 predicted worse OS in ovarian cancer, especially in grade II, grade III, and stage III EOC patients.

S100P was first discovered in placental tissue. Accumulating evidences have suggested a link between S100P expression and progression of various types of cancers [68,69,70,71,72]. In previous studies, statistical analysis has shown S100P positively correlated with serum level of CA125, worse survival and progressive disease in EOC patients [8, 73]. While, Umezaki, Y., et al. reported that S100P low expression was associated with bad prognosis in clear cell adenocarcinoma of the ovary [74]. In our study, we find that high expression of S100P predicted worse OS in grade III, stage I + II and mutant-p53-type of ovarian cancer patients.

As intracellular Ca2+ sensors and as extracellular factors, this family of proteins modulates a wide spectrum of other cellular functions in the normal or non-cancerous conditions. We had reviewed some recently published literatures, and found that S100A12 has pro-inflammatory properties that are likely to be stable in an oxidative environment. Conversely, S100A8 and S100A9 may have important protective mechanisms in inflammation [75]. In endometrial stromal cells, S100A6 over-expression promoted beta-catenin expression at the RNA and protein levels [76]. Additionally, the expression of RAGE and EN-RAGE, the downstream effectors of S100A12, was significantly increased, as evidenced by the significantly greater mRNA and protein expression in the cells of the endometriosis patients [77]. As to hepatic or renal illnesses and diabetes, S100A4 promotes liver fibrosis by activating HSCs (hepatic stellate cell) [78]; serum concentrations of S100B are higher in patients with cirrhosis than in healthy volunteers; silencing S100A9 clearly alleviated G-MDSCs (Granuloid medullary sources inhibit cells) expansion and hepatic steatosis caused by ATF3 deficiency or GC (glucocorticoid) treatment [79]; the sunitinib may protect against renal damage from diabetes mellitus through regulating the levels of vimentin, E-cadherin and S100 [80]; there is coexistence of increased expression of the S100B and the type 2 diabetes mellitus gene in patients [81]; patients with type 2 diabetes have reticulated thrombocytosis that correlates with glycated hemoglobin as well as increased plasma S100A8/A9 levels [82]; S100A8/A9 regulates renal damage and fibrosis, probably through loss of tubular epithelial cell contacts and irreversible damage [83]; higher urine S100 levels are associated with increased lupus nephritis activity in childhood-onset systemic lupus erythematosus, whereas serum S100 (S100A4, S100A6, S100A8/9, and S100A12) levels do not correlate with disease activity [84]; high S100A12 levels are associated with the presence and severity of coronary artery disease in patients with T2DM (type 2 diabetes mellitus) [85].

Our research has some limitations. Firstly, all clinical data were collected by Kaplan-Meier database where the Kaplan-Meier survival curves and the log-rank P value were fully analyzed, and the multivariate analysis by Logistical or COX regression could not be achieved for the reason of data limitation. Thus, the survival according to expression level of S100 family members might be cross- influenced by variable factors such as grade, stage, surgery and so on. Secondly, as for grade, secondary classification system about ovarian serous cancer, namely the low-grade serous carcinoma (LGSC) and high-grade serous carcinoma (HGSC), were recommended by 2014 WHO for the superiority that can be simple and easy to use, high repeatability, better able to guide the clinical treatment. While, the Kaplan-Meier database still used 2003 edition of the WHO classification method, namely, Silverberg grading system, which were divided into tumor Grade 1 (high differentiation), Grade 2 (moderate differentiation), Grade 3 (poorly differentiation). Thirdly, we can only obtain the differential expression information of each S100 family member between normal and ovarian cancer tissues, which provides clues to enhance the points indirectly. (Additional file 1: Table S6).

Conclusion

Comprehensive understanding of the S100 family members may have guiding significance for the diagnosis and outcome in ovarian cancer patients. Based on our study, further discovery of the systematic molecular mechanisms that how S100 interacts with different signaling and other molecules or leads to different prognosis of ovarian cancer patients can pave a way for more effective tumor diagnosis and serve as a genetic treatment target.

Abbreviations

CI:

Confidence interval

EOC:

Epithelial ovarian cancer

GEO:

Gene expression omnibus

GEPIA:

Gene expression profiling interactive analysis

GTEx:

Genotype-tissue expression

HR:

Hazard ratio

KM plotter:

Kaplan Meier plotter

NSCLC:

Non-small-cell lung carcinoma

OS:

Overall survival

PFS:

Progression-free survival

RFS:

Relapse free survival

TCGA:

The cancer genome atlas

References

  1. Zimmer DB, Eubanks JO, Ramakrishnan D, Criscitiello MF. Evolution of the S100 family of calcium sensor proteins. Cell Calcium. 2013;53(3):170–9.

    Article  CAS  PubMed  Google Scholar 

  2. Bresnick AR, Weber DJ, Zimmer DB. S100 proteins in cancer. Nat Rev Cancer. 2015;15(2):96–109.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Zhang S, Wang Z, Liu W, Lei R, Shan J, Li L, Wang X. Distinct prognostic values of S100 mRNA expression in breast cancer. Sci Rep. 2017;7:39786.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Tsai WC, Tsai ST, Jin YT, Wu LW. Cyclooxygenase-2 is involved in S100A2-mediated tumor suppression in squamous cell carcinoma. Mol Cancer Res. 2006;4(8):539–47.

    Article  CAS  PubMed  Google Scholar 

  5. Bulk E, Sargin B, Krug U, Hascher A, Jun Y, Knop M, Kerkhoff C, Gerke V, Liersch R, Mesters RM, et al. S100A2 induces metastasis in non-small cell lung cancer. Clin Cancer Res. 2009;15(1):22–9.

    Article  CAS  PubMed  Google Scholar 

  6. DeRycke MS, Andersen JD, Harrington KM, Pambuccian SE, Kalloger SE, Boylan KL, Argenta PA, Skubitz AP. S100A1 expression in ovarian and endometrial endometrioid carcinomas is a prognostic indicator of relapse-free survival. Am J Clin Pathol. 2009;132(6):846–56.

    Article  CAS  PubMed  Google Scholar 

  7. Wei BR, Hoover SB, Ross MM, Zhou W, Meani F, Edwards JB, Spehalski EI, Risinger JI, Alvord WG, Quinones OA, et al. Serum S100A6 concentration predicts peritoneal tumor burden in mice with epithelial ovarian cancer and is associated with advanced stage in patients. PLoS One. 2009;4(10):e7670.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  8. Wang X, Tian T, Li X, Zhao M, Lou Y, Qian J, Liu Z, Chen H, Cui Z. High expression of S100P is associated with unfavorable prognosis and tumor progression in patients with epithelial ovarian cancer. Am J Cancer Res. 2015;5(8):2409–21.

    PubMed  PubMed Central  Google Scholar 

  9. Maelandsmo GM, Florenes VA, Nguyen MT, Flatmark K, Davidson B. Different expression and clinical role of S100A4 in serous ovarian carcinoma at different anatomic sites. Tumour Biol. 2009;30(1):15–25.

    Article  CAS  PubMed  Google Scholar 

  10. Yang T, Cheng J, Yang Y, Qi W, Zhao Y, Long H, Xie R, Zhu B. S100B mediates Stemness of ovarian Cancer stem-like cells through inhibiting p53. Stem Cells. 2017;35(2):325–36.

    Article  CAS  PubMed  Google Scholar 

  11. Zhao H, Guo E, Hu T, Sun Q, Wu J, Lin X, Luo D, Sun C, Wang C, Zhou B, et al. KCNN4 and S100A14 act as predictors of recurrence in optimally debulked patients with serous ovarian cancer. Oncotarget. 2016;7(28):43924–38.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Cho H, Shin HY, Kim S, Kim JS, Chung JY, Chung EJ, Chun KH, Hewitt SM, Kim JH. The role of S100A14 in epithelial ovarian tumors. Oncotarget. 2014;5(11):3482–96.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Lokman NA, Pyragius CE, Ruszkiewicz A, Oehler MK, Ricciardelli C. Annexin A2 and S100A10 are independent predictors of serous ovarian cancer outcome. Transl Res. 2016;171:83–95.e1–2.

    Article  CAS  PubMed  Google Scholar 

  14. Nymoen DA, Hetland Falkenthal TE, Holth A, Ow GS, Ivshina AV, Trope CG, Kuznetsov VA, Staff AC, Davidson B. Expression and clinical role of chemoresponse-associated genes in ovarian serous carcinoma. Gynecol Oncol. 2015;139(1):30–9.

    Article  CAS  PubMed  Google Scholar 

  15. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98–W102.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wright NT, Cannon BR, Zimmer DB, Weber DJ. S100A1: structure, function, and Therapeutic Potential. Curr Chem Biol. 2009;3(2):138–45.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Li G, Barthelemy A, Feng G, Gentil-Perret A, Peoc'h M, Genin C, Tostain J. S100A1: a powerful marker to differentiate chromophobe renal cell carcinoma from renal oncocytoma. Histopathology. 2007;50(5):642–7.

    Article  CAS  PubMed  Google Scholar 

  18. Lee SWTC, Sager R. Positive selection of candidate tumor-suppressor genes by subtractive hybridization. Proc Natl Acad Sci U S A. 1991;88:2825–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Wolf S, Haase-Kohn C, Pietzsch J. S100A2 in cancerogenesis: a friend or a foe? Amino Acids. 2011;41(4):849–61.

    Article  CAS  PubMed  Google Scholar 

  20. Rehman I, Cross SS, Catto JW, Leiblich A, Mukherjee A, Azzouzi AR, Leung HY, Hamdy FC. Promoter hyper-methylation of calcium binding proteins S100A6 and S100A2 in human prostate cancer. Prostate. 2005;65(4):322–30.

    Article  CAS  PubMed  Google Scholar 

  21. Klopper JP, Sharma V, Bissonnette R, Haugen BR. Combination PPARgamma and RXR agonist treatment in melanoma cells: functional importance of S100A2. PPAR Res. 2010;2010:729876.

    Article  PubMed  CAS  Google Scholar 

  22. Liu D, Rudland PS, Sibson DR, Platt-Higgins A, Barraclough R. Expression of calcium-binding protein S100A2 in breast lesions. Br J Cancer. 2000;83(11):1473–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Ohuchida K, Mizumoto K, Miyasaka Y, Yu J, Cui L, Yamaguchi H, Toma H, Takahata S, Sato N, Nagai E, et al. Over-expression of S100A2 in pancreatic cancer correlates with progression and poor prognosis. J Pathol. 2007;213(3):275–82.

    Article  CAS  PubMed  Google Scholar 

  24. Masuda T, Ishikawa T, Mogushi K, Okazaki S, Ishiguro M, Iida S, Mizushima H, Tanaka H, Uetake H, Sugihara K. Overexpression of the S100A2 protein as a prognostic marker for patients with stage II and III colorectal cancer. Int J Oncol. 2016;48(3):975–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Hough CD, Cho KR, Zonderman AB, Schwartz DR, Morin PJ. Coordinately up-regulated genes in ovarian cancer. Cancer Res. 2001;61(10):3869–76.

    CAS  PubMed  Google Scholar 

  26. Buckley NE, D'Costa Z, Kaminska M, Mullan PB. S100A2 is a BRCA1/p63 coregulated tumour suppressor gene with roles in the regulation of mutant p53 stability. Cell Death Dis. 2014;5:e1070.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Mueller A, Schafer BW, Ferrari S, Weibel M, Makek M, Hochli M, Heizmann CW. The calcium-binding protein S100A2 interacts with p53 and modulates its transcriptional activity. J Biol Chem. 2005;280(32):29186–93.

    Article  CAS  PubMed  Google Scholar 

  28. Liu J, Li X, Dong GL, Zhang HW, Chen DL, Du JJ, Zheng JY, Li JP, Wang WZ. In silico analysis and verification of S100 gene expression in gastric cancer. BMC Cancer. 2008;8:261.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  29. Lloyd BH, Ruddell C, Rudland PS, Barraclough R. S100A3 mRNA expression displays an inverse correlation to breast cancer progression. Biochem Soc Trans. 1996;24(3):340S.

    Article  CAS  PubMed  Google Scholar 

  30. Liu B, Sun WY, Zhi CY, Lu TC, Gao HM, Zhou JH, Yan WQ, Gao HC. Role of S100A3 in human colorectal cancer and the anticancer effect of cantharidinate. Exp Ther Med. 2013;6(6):1499–503.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Kang M, Lee HS, Lee YJ, Choi WS, Park YH, Jeong CW, Ku JH, Kim HH, Kwak C. S100A3 suppression inhibits in vitro and in vivo tumor growth and invasion of human castration-resistant prostate cancer cells. Urology. 2015;85(1):273.e9–15.

    Article  PubMed  Google Scholar 

  32. Schafer BW, Fritschy JM, Murmann P, Troxler H, Durussel I, Heizmann CW, Cox JA. Brain S100A5 is a novel calcium-, zinc-, and copper ion-binding protein of the EF-hand superfamily. J Biol Chem. 2000;275(39):30623–30.

    Article  CAS  PubMed  Google Scholar 

  33. Cho CC, Chou RH, Yu C. Pentamidine blocks the interaction between mutant S100A5 and RAGE V domain and inhibits the RAGE signaling pathway. Biochem Biophys Res Commun. 2016;477(2):188–94.

    Article  CAS  PubMed  Google Scholar 

  34. Hancq S, Salmon I, Brotchi J, De Witte O, Gabius HJ, Heizmann CW, Kiss R, Decaestecker C. S100A5: a marker of recurrence in WHO grade I meningiomas. Neuropathol Appl Neurobiol. 2004;30(2):178–87.

    Article  CAS  PubMed  Google Scholar 

  35. Donato R, Sorci G, Giambanco I. S100A6 protein: functional roles. Cell Mol Life Sci. 2017.

  36. Lesniak W, Slomnicki LP, Filipek A. S100A6 - new facts and features. Biochem Biophys Res Commun. 2009;390(4):1087–92.

    Article  CAS  PubMed  Google Scholar 

  37. Nedjadi T, Kitteringham N, Campbell F, Jenkins RE, Park BK, Navarro P, Ashcroft F, Tepikin A, Neoptolemos JP, Costello E. S100A6 binds to annexin 2 in pancreatic cancer cells and promotes pancreatic cancer cell motility. Br J Cancer. 2009;101(7):1145–54.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Wang XH, Zhang LH, Zhong XY, Xing XF, Liu YQ, Niu ZJ, Peng Y, Du H, Zhang GG, Hu Y, et al. S100A6 overexpression is associated with poor prognosis and is epigenetically up-regulated in gastric cancer. Am J Pathol. 2010;177(2):586–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Hua Z, Chen J, Sun B, Zhao G, Zhang Y, Fong Y, Jia Z, Yao L. Specific expression of osteopontin and S100A6 in hepatocellular carcinoma. Surgery. 2011;149(6):783–91.

    Article  PubMed  Google Scholar 

  40. Brown LM, Helmke SM, Hunsucker SW, Netea-Maier RT, Chiang SA, Heinz DE, Shroyer KR, Duncan MW, Haugen BR. Quantitative and qualitative differences in protein expression between papillary thyroid carcinoma and normal thyroid tissue. Mol Carcinog. 2006;45(8):613–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Slomnicki LP, Nawrot B, Lesniak W. S100A6 binds p53 and affects its activity. Int J Biochem Cell Biol. 2009;41(4):784–90.

    Article  CAS  PubMed  Google Scholar 

  42. Roth J, Vogl T, Sorg C, Sunderkotter C. Phagocyte-specific S100 proteins: a novel group of proinflammatory molecules. Trends Immunol. 2003;24(4):155–8.

    Article  CAS  PubMed  Google Scholar 

  43. Cheng P, Corzo CA, Luetteke N, Yu B, Nagaraj S, Bui MM, Ortiz M, Nacken W, Sorg C, Vogl T, et al. Inhibition of dendritic cell differentiation and accumulation of myeloid-derived suppressor cells in cancer is regulated by S100A9 protein. J Exp Med. 2008;205(10):2235–49.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Srikrishna G. S100A8 and S100A9: new insights into their roles in malignancy. J Innate Immun. 2012;4(1):31–40.

    Article  CAS  PubMed  Google Scholar 

  45. Foell D, Wittkowski H, Vogl T, Roth J. S100 proteins expressed in phagocytes: a novel group of damage-associated molecular pattern molecules. J Leukoc Biol. 2007;81(1):28–37.

    Article  CAS  PubMed  Google Scholar 

  46. Gebhardt C, Nemeth J, Angel P, Hess J. S100A8 and S100A9 in inflammation and cancer. Biochem Pharmacol. 2006;72(11):1622–31.

    Article  CAS  PubMed  Google Scholar 

  47. Arai K, Takano S, Teratani T, Ito Y, Yamada T, Nozawa R. S100A8 and S100A9 overexpression is associated with poor pathological parameters in invasive ductal carcinoma of the breast. Curr Cancer Drug Targets. 2008;8(4):243–52.

    Article  CAS  PubMed  Google Scholar 

  48. Kim HJ, Kang HJ, Lee H, Lee ST, Yu MH, Kim H, Lee C. Identification of S100A8 and S100A9 as serological markers for colorectal cancer. J Proteome Res. 2009;8(3):1368–79.

    Article  CAS  PubMed  Google Scholar 

  49. Nemeth J, Stein I, Haag D, Riehl A, Longerich T, Horwitz E, Breuhahn K, Gebhardt C, Schirmacher P, Hahn M, et al. S100A8 and S100A9 are novel nuclear factor kappa B target genes during malignant progression of murine and human liver carcinogenesis. Hepatology. 2009;50(4):1251–62.

    Article  CAS  PubMed  Google Scholar 

  50. Kwon CH, Moon HJ, Park HJ, Choi JH, Park DY. S100A8 and S100A9 promotes invasion and migration through p38 mitogen-activated protein kinase-dependent NF-kappaB activation in gastric cancer cells. Mol Cells. 2013;35(3):226–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Ghavami S, Chitayat S, Hashemi M, Eshraghi M, Chazin WJ, Halayko AJ, Kerkhoff C. S100A8/A9: a Janus-faced molecule in cancer therapy and tumorgenesis. Eur J Pharmacol. 2009;625(1–3):73–83.

    Article  CAS  PubMed  Google Scholar 

  52. Kwon M, Yoon CS, Jeong W, Rhee SG, Waisman DM. Annexin A2-S100A10 heterotetramer, a novel substrate of thioredoxin. J Biol Chem. 2005;280(25):23584–92.

    Article  CAS  PubMed  Google Scholar 

  53. Kwon M, MacLeod TJ, Zhang Y, Waisman DM. S100A10, annexin A2, and annexin a2 heterotetramer as candidate plasminogen receptors. Front Biosci. 2005;10:300–25.

    Article  CAS  PubMed  Google Scholar 

  54. Zhang L, Fogg DK, Waisman DM. RNA interference-mediated silencing of the S100A10 gene attenuates plasmin generation and invasiveness of Colo 222 colorectal cancer cells. J Biol Chem. 2004;279(3):2053–62.

    Article  CAS  PubMed  Google Scholar 

  55. Suzuki S, Tanigawara Y. Forced expression of S100A10 reduces sensitivity to oxaliplatin in colorectal cancer cells. Proteome Sci. 2014;12:26.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  56. Ito Y, Arai K, Nozawa R, Yoshida H, Higashiyama T, Takamura Y, Miya A, Kobayashi K, Kuma K, Miyauchi A. S100A10 expression in thyroid neoplasms originating from the follicular epithelium: contribution to the aggressive characteristic of anaplastic carcinoma. Anticancer Res. 2007;27(4C):2679–83.

    CAS  PubMed  Google Scholar 

  57. Azimi A, Pernemalm M, Frostvik Stolt M, Hansson J, Lehtio J, Egyhazi Brage S, Hertzman Johansson C. Proteomics analysis of melanoma metastases: association between S100A13 expression and chemotherapy resistance. Br J Cancer. 2014;110(10):2489–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Massi D, Landriscina M, Piscazzi A, Cosci E, Kirov A, Paglierani M, Di Serio C, Mourmouras V, Fumagalli S, Biagioli M, et al. S100A13 is a new angiogenic marker in human melanoma. Mod Pathol. 2010;23(6):804–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Landriscina M, Schinzari G, Di Leonardo G, Quirino M, Cassano A, D'Argento E, Lauriola L, Scerrati M, Prudovsky I, Barone C. S100A13, a new marker of angiogenesis in human astrocytic gliomas. J Neuro-Oncol. 2006;80(3):251–9.

    Article  CAS  Google Scholar 

  60. Pierce A, Barron N, Linehan R, Ryan E, O'Driscoll L, Daly C, Clynes M. Identification of a novel, functional role for S100A13 in invasive lung cancer cell lines. Eur J Cancer. 2008;44(1):151–9.

    Article  CAS  PubMed  Google Scholar 

  61. Sturchler E, Cox JA, Durussel I, Weibel M, Heizmann CW. S100A16, a novel calcium-binding protein of the EF-hand superfamily. J Biol Chem. 2006;281(50):38905–17.

    Article  CAS  PubMed  Google Scholar 

  62. Katono K, Sato Y, Kobayashi M, Nagashio R, Ryuge S, Igawa S, Ichinoe M, Murakumo Y, Saegusa M, Masuda N. S100A16, a promising candidate as a prognostic marker for platinum-based adjuvant chemotherapy in resected lung adenocarcinoma. Onco Targets Ther. 2017;10:5273–9.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Marenholz I, Heizmann CW. S100A16, a ubiquitously expressed EF-hand protein which is up-regulated in tumors. Biochem Biophys Res Commun. 2004;313(2):237–44.

    Article  CAS  PubMed  Google Scholar 

  64. Tanaka M, Ichikawa-Tomikawa N, Shishito N, Nishiura K, Miura T, Hozumi A, Chiba H, Yoshida S, Ohtake T, Sugino T. Co-expression of S100A14 and S100A16 correlates with a poor prognosis in human breast cancer and promotes cancer cell invasion. BMC Cancer. 2015;15:53.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  65. Zhou W, Pan H, Xia T, Xue J, Cheng L, Fan P, Zhang Y, Zhu W, Xue Y, Liu X, et al. Up-regulation of S100A16 expression promotes epithelial-mesenchymal transition via Notch1 pathway in breast cancer. J Biomed Sci. 2014;21:97.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  66. Sapkota D, Bruland O, Parajuli H, Osman TA, Teh MT, Johannessen AC, Costea DE. S100A16 promotes differentiation and contributes to a less aggressive tumor phenotype in oral squamous cell carcinoma. BMC Cancer. 2015;15:631.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  67. Sun X, Wang T, Zhang C, Ning K, Guan ZR, Chen SX, Hong TT, Hua D. S100A16 is a prognostic marker for colorectal cancer. J Surg Oncol. 2017.

  68. Peng C, Chen H, Wallwiener M, Modugno C, Cuk K, Madhavan D, Trumpp A, Heil J, Marme F, Nees J, et al. Plasma S100P level as a novel prognostic marker of metastatic breast cancer. Breast Cancer Res Treat. 2016;157(2):329–38.

    Article  CAS  PubMed  Google Scholar 

  69. Ezzat NE, Tahoun NS, Ismail YM. The role of S100P and IMP3 in the cytologic diagnosis of pancreatic adenocarcinoma. J Egypt Natl Canc Inst. 2016;28(4):229–34.

    Article  PubMed  Google Scholar 

  70. Chung L, Phillips L, Lin MZ, Moore K, Marsh DJ, Boyle FM, Baxter RC. A novel truncated form of S100P predicts disease-free survival in patients with lymph node positive breast cancer. Cancer Lett. 2015;368(1):64–70.

    Article  CAS  PubMed  Google Scholar 

  71. Jiang H, Hu H, Lin F, Lim YP, Hua Y, Tong X, Zhang S. S100P is overexpressed in squamous cell and Adenosquamous carcinoma subtypes of endometrial Cancer and promotes Cancer cell proliferation and invasion. Cancer Investig. 2016;34(10):477–88.

    Article  CAS  Google Scholar 

  72. Shen ZY, Fang Y, Zhen L, Zhu XJ, Chen H, Liu H, Jiang B, Li GX, Deng HJ. Analysis of the predictive efficiency of S100P on adverse prognosis and the pathogenesis of S100P-mediated invasion and metastasis of colon adenocarcinoma. Cancer Genet. 2016;209(4):143–53.

    Article  CAS  PubMed  Google Scholar 

  73. Surowiak P, Maciejczyk A, Materna V, Drag-Zalesinska M, Wojnar A, Pudelko M, Kedzia W, Spaczynski M, Dietel M, Zabel M, et al. Unfavourable prognostic significance of S100P expression in ovarian cancers. Histopathology. 2007;51(1):125–8.

    Article  CAS  PubMed  Google Scholar 

  74. Umezaki Y, Ito M, Nakashima M, Mihara Y, Kurohama H, Yatsunami N, Yasuhi I. Low expression of S100P is associated with poor prognosis in patients with clear cell adenocarcinoma of the ovary. Int J Gynecol Cancer. 2015;25(9):1582–6.

    Article  PubMed  Google Scholar 

  75. Perera C, McNeil HP, Geczy CL. S100 Calgranulins in inflammatory arthritis. Immunol Cell Biol. 2010;88(1):41–9.

    Article  CAS  PubMed  Google Scholar 

  76. Liu Z, Zhang X, Chen M, Cao Q, Huang D. Effect of S100A6 over-expression on beta-catenin in endometriosis. J Obstet Gynaecol Res. 2015;41(9):1457–62.

    Article  CAS  PubMed  Google Scholar 

  77. Sharma I, Dhawan V, Saha SC, Rashmi B, Dhaliwal LK. Implication of the RAGE-EN-RAGE axis in endometriosis. Int J Gynaecol Obstet. 2010;110(3):199–202.

    Article  CAS  PubMed  Google Scholar 

  78. Chen L, Li J, Zhang J, Dai C, Liu X, Wang J, Gao Z, Guo H, Wang R, Lu S, et al. S100A4 promotes liver fibrosis via activation of hepatic stellate cells. J Hepatol. 2015;62(1):156–64.

    Article  CAS  PubMed  Google Scholar 

  79. Liu YF, Wei JY, Shi MH, Jiang H, Zhou J. Glucocorticoid induces hepatic steatosis by inhibiting activating transcription factor 3 (ATF3)/S100A9 protein signaling in granulocytic myeloid-derived suppressor cells. J Biol Chem. 2016;291(41):21771–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Akarca-Dizakar SO, Aktug H, Oltulu F, Oktem G, Yavasoglu A, Acikgoz E, Yigitturk G, Demir K, Uysal A. Effects of sunitinib on immunoreactivity of vimentin, E-cadherin and S100 in kidneys of streptozotocin induced diabetic mice. Biotech Histochem. 2018;93(5):328–35.

    Article  CAS  PubMed  Google Scholar 

  81. Hamasaki MY, Hirata MH, Hirata RD, Himelfarb ST, Campos LM, Nogueira MI. Analysis of the mRNA expression of the S100beta protein in adipocytes of patients with diabetes mellitus, type 2. Arq Bras Endocrinol Metabol. 2012;56(7):435–40.

    Article  PubMed  Google Scholar 

  82. Kraakman MJ, Lee MK, Al-Sharea A, Dragoljevic D, Barrett TJ, Montenont E, Basu D, Heywood S, Kammoun HL, Flynn M, et al. Neutrophil-derived S100 calcium-binding proteins A8/A9 promote reticulated thrombocytosis and atherogenesis in diabetes. J Clin Invest. 2017;127(6):2133–47.

    Article  PubMed  PubMed Central  Google Scholar 

  83. Tammaro A, Florquin S, Brok M, Claessen N, Butter LM, Teske GJD, de Boer OJ, Vogl T, Leemans JC, Dessing MC. S100A8/A9 promotes parenchymal damage and renal fibrosis in obstructive nephropathy. Clin Exp Immunol. 2018;193(3):361–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Turnier JL, Fall N, Thornton S, Witte D, Bennett MR, Appenzeller S, Klein-Gitelman MS, Grom AA, Brunner HI. Urine S100 proteins as potential biomarkers of lupus nephritis activity. Arthritis Res Ther. 2017;19(1):242.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  85. Zhao P, Wu M, Yu H, Huang Y, Wang Y, Wang W, Yin W. Serum S100A12 levels are correlated with the presence and severity of coronary artery disease in patients with type 2 diabetes mellitus. J Investig Med. 2013;61(5):861–6.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank all of the patients for participating in our study.

Funding

This work was supported by the Natural Science Foundation of China (NSFC) grants awarded to Xin Lu (17ZR1403500) and by the Shanghai Sailing Program grants awarded to Jun Li (16YF1401100).

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

YB and LDL participated in the design of the study, collected and analyzed the data. JL and XL participated in its design and coordination, and supervised the study. XL drafted the manuscript. YB and LDL contributed equally to this work. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xin Lu.

Ethics declarations

Ethics approval and consent to participate

It is not applicable for this section, as our manuscript does not contain any individual personal information, and data from “Kaplan Meier Plotter (http://kmplot.com/analysis/)” and “GEPIA (http://gepia.cancer-pku.cn/index.html)” are reproduced with their permissions.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. The mRNA expression of some S100 family members is of no apparent prognostic significance in all ovarian cancer patients. Abbreviations: exp.: expression; OS: overall survival; HR: hazard ratio; CI: confidence interval; Table S2. The mRNA expression of some S100 family members is of no apparent prognostic significance in all serous ovarian cancer patients. Abbreviations: exp.: expression; OS: overall survival; HR: hazard ratio; CI: confidence interval; Table S3. The mRNA expression of some S100 family members is of no apparent prognostic significance in all endometrioid ovarian cancer patients. Abbreviations: exp.: expression; OS: overall survival; HR: hazard ratio; CI: confidence interval; Table S4. Prognostic value of mRNA expression of S100 family members in ovarian cancer patients with optimal or suboptimal surgery. Abbreviations: exp. -expression, HR-hazard ratio, CI-confidence intervals; Table S5. Prognostic value of mRNA expression of S100 family members in ovarian cancer patients with chemotherapy. Abbreviations: exp. -expression, HR-hazard ratio, CI-confidence intervals; Table S6. The mRNA expression levels of S100 family members in ovarian normal vs. tumor patients. Abbreviations: Min-minimal expression level, Q1-one quarter, Q3-three quarter, Max-max expression level, n-number (XLSX 29 kb)

Additional file 2:

Table S7. Patient numbers of low or high mRNA expression of S100 family members in different clinical and pathological features. Abbreviations: exp. -expression, n-number. (XLSX 41 kb)

Additional file 3:

Figure S2. Differential expression of S100 family members in normal and tumor tissues across 31 tumor types from the TCGA and GTEx data. Abbreviation: N-normal tissues; T-tumor tissues; ACC-Adrenocortical carcinoma; BLCA-Bladder Urothelial Carcinoma; BRCA-Breast invasive carcinoma; CESC-Cervical squamous cell carcinoma and endocervical adenocarcinoma; CHOL-Cholangio carcinoma; COAD-Colon adenocarcinoma; DLBC-Lymphoid Neoplasm Diffuse Large B-cell Lymphoma; ESCA-Esophageal carcinoma; GBM-Glioblastoma multiforme; HNSC-Head and Neck squamous cell carcinoma; KICH-Kidney Chromophobe; KIRC-Kidney renal clear cell carcinoma; KIRP-Kidney renal papillary cell carcinoma; LAML-Acute Myeloid Leukemia; LGG-Brain Lower Grade Glioma; LIHC-Liver hepatocellular carcinoma; LUAD-Lung adenocarcinoma; LUSC-Lung squamous cell carcinoma; MESO-Mesothelioma; OV-Ovarian serous cystadenocarcinoma; PAAD-Pancreatic adenocarcinoma; PCPG-Pheochromocytoma and Paraganglioma; PRAD-Prostate adenocarcinoma; READ-Rectum adenocarcinoma; SARC-Sarcoma; SKCM-Skin Cutaneous Melanoma; STAD-Stomach adenocarcinoma; TGCT-Testicular Germ Cell Tumors; THCA-Thyroid carcinoma; THYM-Thymoma; UCEC-Uterine Corpus Endometrial Carcinoma; UCS-Uterine Carcinosarcoma; UVM-Uveal Melanoma. (TIF 849 kb)

Additional file 4:

Table S8. Differential expression of S100 family members in normal and tumor tissues across 31 tumor types from the TCGA and GTEx data. Abbreviations: exp. -expression, n-number. (XLSX 16 kb)

Additional file 5:

Figure S1. Progressive free survival of S100 members in different ovarian cancer subtypes. Notes: (A-C) Survival curves of S100A1 (the desired Affymetrix IDs is valid: 205334_at) are plotted for all/endometrioid/serous ovarian cancer patients. Abbreviation: HR: hazard ratio (TIF 131 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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

Bai, Y., Li, LD., Li, J. et al. Prognostic values of S100 family members in ovarian cancer patients. BMC Cancer 18, 1256 (2018). https://doi.org/10.1186/s12885-018-5170-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-018-5170-3

Keywords