Upregulation of STC2 in colorectal cancer and its clinicopathological significance

Abstract

Background

Stanniocalcin 2 (STC2) is a glycoprotein hormone involved in many biological processes and a secretory protein that regulates malignant tumor progression. The aim of the present study was to further explore the clinicopathological significance and prognostic role of STC2 in colorectal cancer (CRC).

Methods

In this study, STC2 expression was first investigated in Gene Expression Omnibus and The Cancer Genome Atlas, and then validated with the data from our medical center. Univariate and multivariate analyses were performed to assess the association between prognostic factors and survival outcome.

Results

In Gene Expression Omnibus and The Cancer Genome Atlas databases, bioinformatics analysis confirmed that STC2 was significantly increased in CRC compared with that in normal tissues (P<0.01), and CRC patients with high STC2 expression had a shorter overall survival. By analyzing data from our medical center, the results also showed that STC2 expression of CRC tissues was higher than that in normal tissues, whether the transcriptional or protein levels. In the CRC tissues, high STC2 expression was significantly correlated with lymph node metastasis (P=0.047), distant metastasis (P=0.040), and advanced clinical stage (P=0.047). Moreover, Kaplan–Meier analyses indicated that high STC2 expression predicted a worse prognosis, and multivariate Cox regression analysis revealed that STC2 was an independent prognostic factor for overall survival (HR =1.976, 95% CI: 1.092–3.576, P=0.024) in patients with CRC.

Conclusion

Our results suggested that STC2 played an important role in CRC progression and prognosis, and could be a useful biomarker for survival prediction.

Keywords:

• STC2
• prognosis
• bioinformatics analysis
• colorectal cancer

Introduction

Colorectal cancer (CRC) has been the fourth leading cause of cancer death worldwide for several decades.¹ The 5-year survival rate for localized stage CRC is 90%, whereas as CRC spreads to the regional lymph nodes or distant parts of the body, the 5-year survival rate plunges from 71% to 14%.² Although improvements in treatment protocols, including the discovery of targeted therapies, have substantially improved the clinical outcomes, most CRC patients fail treatment due to local recurrence and distant metastasis.¹ Therefore, predicting the risk of distant metastasis and poor prognosis is essential to improve clinical management of CRC. Unfortunately, most of the current medical standard tumor markers lack sensitivity and specificity for the detection of CRC. For the advances of the diagnostic methods, molecular-based diagnostics have been studied and proposed by many groups. Proteomic or genomic approach, based on profiling differentially expressed proteins or genes, has been used for novel biomarker discovery.^3–5

We had analyzed differently expressed genes (DEGs) between CRC and paired adjacent normal tissues (PANT) by gene expression microarray and found high expression of the stanniocalcin 2 (STC2) gene in CRC.⁶ STC2 is a glycoprotein hormone first known to be involved in calcium and phosphate homeostasis.^7–9 In addition, STC2 was involved in glutamine or glucose deprivation, and it also had been found to be upregulated under hypoxia, endoplasmic reticulum stress, and radiation.^10–13 Microenvironment characterized by hypoxia and low glucose and glutamine supply facilitated the process of solid tumors. STC2 overexpression contributes to tumor cells’ adaptation to such stress conditions, thus promoting tumor progression.¹⁰ Recent studies demonstrated that STC2 expression was upregulated in various tumors, including lung cancer,¹⁴ breast cancer,¹⁵ hepatocellular carcinoma,¹⁶ gastric cancer,¹⁷ esophageal carcinoma,¹⁸ and nasopharyngeal carcinoma.¹⁹ Clinical and pathological studies reveal that STC2 overexpression correlates with advanced tumor grade, tumor invasiveness, metastasis, and poor prognosis.^16–19 On the other hand, STC2 expression was also upregulated in breast cancer patients; however, the prognosis for these cases was good.^20,21 The clinical significance of STC2 in cancer is disputed and depends on cancer type. Furthermore, clinicopathological significance of STC2 gene expression in CRC had been reported in 2009.²² However, its clinical significance and molecular mechanism in carcinogenesis are still not completely understood. Therefore, to further explore the precise role of STC2 for CRC diagnosis and prognosis is essential.

In order to draw a solid conclusion, we first studied STC2 expression in the publicly available Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases, and then validated it with the data from our medical center. We evaluated the expression of STC2 in several CRC specimens and investigated its association with clinicopathological parameters and overall survival (OS).

Materials and methods

Microarray data processing

The microarray datasets of CRC, including GSE21510, GSE32323, GSE39582, and GSE41328, were downloaded from the public GEO databases (http://www.ncbi.nlm.nih. gov/geo/). Robust Multichip Average, an algorithm used to create an expression matrix from Affymetrix data, was used to adjust the raw files of background. Median levels of transcript expressions were calculated. Gene-level data were then filtered to include only those probe sets with annotations. Probe set annotation mainly reference the new version annotation files that can be download on Affymetrix official website (http://www.affymetrix.com/support/technical/byproduct.affx).

Differentially expressed genes mining

The limma package (http://bioconductor.org/packages/ release/bioc/html/limma.html) was used for selection of DEGs.²³ The limma package was a correct and popular method for gene selection through differential expression analyses of microarray. Genes that met the cutoff criteria of adjusted P<0.01 and |log₂fold change (FC)| ≥2.0 were screened out as DEGs. Among DEGs, STC2 gene was highly expressed in CRC tissues, similar to the result of our microarray data.⁶ Based on the previous research foundation, the STC2 gene would be further explored.

Validation and survival analysis

To further validate our results, we employed the TCGA database. Expression of the STC2 mRNA in CRC tissues and normal tissues was compared. Kaplan–Meier plotting data of OS of patients with CRC were generated. The patients were classified into high and low STC2 expression groups according to the median STC2 expression level.

Tissue collection

All tissues were acquired from patients who had been diagnosed with CRC by pathological assessment and undergone surgeries at the Zhongnan Hospital of Wuhan University. The study was approved by the Research Ethics Committee of Zhongnan Hospital of Wuhan University. All patients had been adequately informed and written consent was obtained. This was conducted in accordance with the Declaration of Helsinki.

Six CRC tissues and PANT (distance to cancer >5 cm), obtained from patients after surgery between June 2012 and November 2012, made into tissue microarray.⁶ None of the patients received preoperative radiotherapy or chemotherapy. One-hundred fifteen CRC paraffin-embedded tissues were obtained from the Department of Pathology. The patients received surgical treatment between January 2014 and May 2015 with complete clinicopathological and prognostic information. STC2 protein expression was assessed in these tissues using immunohistochemical staining. Follow-ups were terminated by August 2018.

Immunohistochemical staining

Immunohistochemistry (IHC) of the tissues was performed as described previously.²⁴ Briefly, deparaffinized sections were pretreated with 10 mM sodium citrate buffer for antigen unmasking (pH 6.0, boiling temperature, 30 minutes), blocked in normal serum, incubated with primary antibody of STC2 (1:100, ab63057; Abcam, Cambridge, UK) at 4°C overnight, rinsed, and incubated with secondary antibody (MaxVision™ HRP-Polymer IHC Kit, Maixin-Bio, Fuzhou, China). Signals were amplified using the kit per the manufacturer’s instructions. The positive brown staining was visualized and then photographed using a light microscope at 200× magnification (BX51; Olympus Corporation, Tokyo, Japan).

The IHC evaluation was performed by two independent observers (Xiang Ma and Xiang Shu, the co-authors) who were blinded to the clinical and pathological characteristics associated with the specimens. The extent of STC2 staining was scored by a semiquantitative method in which staining of >10% of the tumor cells was considered positive. The staining intensity was scored as 0 (negative,−), 1 (weak positive, +), 2 (moderated positive, ++), and 3 (strong positive, +++). For analysis, the STC2 protein expression levels were divided into two groups: low expression level (score value <2) and high expression level (score value ≥2).

Western blot

Western blot was completed as described formerly.²⁵ Total protein was extracted with cell lysis buffer (Beyotime, Shanghai, China), and the protein concentration was quantified using an Enhanced BCA Protein Assay Kit (Beyotime). Protein was separated by 10% SDS-PAGE and transferred by semidry electroblotting to polyvinylidene difluoride membranes (Millipore, Billerica, MA, USA). The membrane was blocked for 1 hour with 5% BSA in TBS-T, and probed with corresponding primary antibodies overnight at 4°C, followed by incubation with rabbit or mouse horseradish peroxidase-coupled secondary antibodies for 1 hour. Specific bands were detected using Clarity™ Western ECL Substrate (Bio-Rad, Richmond, CA, USA) on an autoradiographic film. The primary antibodies used were as follows: anti-STC2 (Abcam) and anti-GAPDH (Proteintech, Wuhan, China).

Quantitative real-time polymerase chain reaction

RNA extraction and quantitative real-time polymerase chain reaction (qRT-PCR) were carried out as previously described.²⁵ Total RNA was isolated using RNA plus reagent (Vazyme, Nanjing, China). Complementary DNA was prepared using oligodT primers according to the protocol supplied with the HiScript II Q RT SuperMix for qPCR (Vazyme). Expression of STC2 was determined by qRT-PCR using AceQqPCR SYBR Green Master Mix (Vazyme). Oligonucleotide sequences of the primer sets used were as follows: STC2 forward (5′-TGAAATGTAAGGCCCACGCT-3′) and reverse (5′-CGAGGTGCAGAAGCTCAAGA-3′) and ACTIN forward (5′-CACCCAGCACAATGAAGATCAAGAT-3′) and reverse (5′-CCAGTTTTTAAATCCTGAGTCAAGC-3′).

Statistical analyses

R 3.5.0 was used for microarray analysis. Statistical analyses were performed with the IBM SPSS 22.0 statistical software package (IBM Inc.). The relationship between STC2 expression and clinicopathological features was analyzed using the Chi-square tests. Kaplan–Meier survival curves were constructed, and the log-rank test was carried out using univariate analysis. Multivariate analysis was performed using Cox’s proportional hazards model. P<0.05 was considered statistically significant.

Results

Microarray datasets

The mRNA expression profile datasets, GSE21510, GSE32323, GSE39582, and GSE41328, were generated using the GPL570 microarray platform. One hundred thirty-two samples, consisting of 69 normal tissues and 63 CRC tissues, were employed (Table S1).

Differentially expressed genes

Using R 3.5.0, we identified 1,627, 1,236, 1,607, and 762 DEGs between CRC and normal tissues from the GSE21510, GSE32323, GSE39582, and GSE41328 datasets, respectively (P<0.01, Figures 1A and S1). We found that 381 genes were differentially expressed in all four datasets (Figure 1B). Among the DEGs, STC2 gene was highly expressed in CRC tissue, similar to the result of our microarray data (Figure 2A). Based on the previous research foundation, the STC2 gene underwent further evaluation.⁶

Figure 1 Increase of STC2 expression in GEO and TCGA databases in CRC.

Notes: (A) DEGs in CRC vs normal tissues across each independent dataset. Each column represents a sample and each row represents the expression level of a given mRNA. The color scale represents the raw Z score ranging from blue (low expression) to red (high expression). (B) Venn diagrams of overlapping DEGs between GSE21510, GSE32323, GSE39582, and GSE41328. (C) Analysis of STC2 expression in unpaired CRC (N=349) and normal tissues (N=41) in TCGA (P<0.001). (D) Analysis of STC2 expression in CRC and PANT (N=41) in TCGA (P<0.001). (E) High STC2 expression was associated with poor OS in CRC patients in TCGA.

Abbreviations: CRC, colorectal cancer; DEGs, differently expressed genes; GEO, Gene Expression Omnibus; OS, overall survival; PANT, paired adjacent normal tissues; STC2, stanniocalcin 2; TCGA, The Cancer Genome Atlas.

Figure 2 High expression of STC2 validated by our own database in CRC.

Notes: (A) DEGs in six pairs of CRC and PANT. (B) STC2 protein and mRNA levels were measured in CRC and PANT by Western blotting and qRT-PCR, respectively.

Abbreviations: CRC, colorectal cancer; DEGs, differently expressed genes; PANT, paired adjacent normal tissues; qRT-PCR, quantitative real-time polymerase chain reaction; STC2, stanniocalcin 2.

Validation and survival analyses

To evaluate the expression of STC2 in CRC, we further analyzed STC2 expression in TCGA datasets. The mRNA level of STC2 in the CRC tissues (N=349) and normal tissues (N=41) from TCGA showed that STC2 mRNA expression was significantly upregulated in cancer tissues (P<0.001, Figure 1C). Meanwhile, in 41 paired tissues, the expression of STC2 mRNA was dramatically higher in cancer tissues than that in the PANT (P<0.001, Figure 1D). To further investigate whether the deregulated STC2 correlates with the survivals of the CRC patients, we performed Kaplan–Meier and Cox’s proportional hazards regression model analysis, and found that high STC2 was significantly correlated with poor OS (Figure 1E; Table 1).

Table 1 Univariate and multivariate analyses of clinicopathological characteristics and STC2 with overall survival in TCGA

	Univariate analysis		Multivariate analysis
	HR (95% CI)	P-value	HR (95% CI)	P-value
Gender (male vs female)	0.520 (0.264–1.022)	0.058
Age, years (≥median vs <median)	1.756 (0.903–3.413)	0.097
Vascular invasion (present vs absent)	2.459 (1.280–4.725)	0.007	1.275 (0.472–3.448)	0.632
Lymphatic invasion (present vs absent)	2.361 (1.228–4.537)	0.010	1.327 (0.486–3.627)	0.581
Tumor infiltration (T3 and T4 vs T1 and T2)	2.060 (0.489–8.673)	0.324
Lymph node metastasis (present vs absent)	1.743 (0.911–3.333)	0.093
Node counts (≥12 vs <12)	0.486 (0.227–1.042)	0.064
Distant metastasis (present vs absent)	3.562 (1.804–7.031)	<0.001	2.444 (1.004–5.950)	0.049
Clinical stage (III and IV vs I and II)	2.237 (1.147–4.361)	0.018	1.174 (0.482–2.861)	0.725
STC2 (≥median vs <median)	2.338 (1.189–4.599)	0.014	2.064 (1.035–4.113)	0.040

Note: Bold indicates statistically significant values (P<0.05).

Abbreviations: STC2, stanniocalcin 2; TCGA, The Cancer Genome Atlas.

STC2 is overexpressed in colorectal cancer

As STC2 showed a consistent pattern of overexpression in GEO and TCGA datasets, our next research direction was to explore the expression pattern of STC2 in the CRC patients of our medical center. In six pairs of specimens, microarray data showed that STC2 gene was significantly upregulated in CRC tissues (Figure 2A). STC2 protein expression was also obviously increased in CRC tissues as determined by Western blot (Figure 2B). Similarly, according to the results of qRT-PCR, STC2 mRNA expression was significantly upregulated in CRC tissues, compared with that in the PANT (Figure 2B).

Associations of STC2 expression with clinicopathological parameters in CRC

The expression of STC2 protein in archived CRC tissue samples was detected by IHC assay. We firstly evaluated STC2 protein expression in 15 pairs of CRC tissues and PANT. STC2 protein was mainly distributed in the cytoplasm and partially in the nucleus. The protein expression of STC2 was significantly increased in CRC compared with that in PANT tissues (Figure 3A).

Figure 3 Detection of STC2 protein expression in tissues by IHC staining, and high STC2 expression in tumors predicted poor prognosis in CRC patients.

Notes: (A) Representative images of STC2 protein expression in CRC and PANT. STC2 protein expression was significantly increased in tumor specimens compared with PANT. (B) IHC staining of STC2 protein in CRC tissues. The staining intensity was scored as 0 (negative), 1 (weak positive), 2 (moderated positive), and 3 (strong positive). (C) High expression of STC2 had a worse OS for CRC patients.

Abbreviations: CRC, colorectal cancer; IHC, immunohistochemistry; OS, overall survival; PANT, paired adjacent normal tissues; STC2, stanniocalcin 2.

The association of STC2 expression and clinicopathological features including outcome was analyzed in 115 cases of CRC patients in our medical center. Of the 115 CRC patients, 28 patients (24.35%) had negative STC2 staining (IHC score, 0), 35 patients (30.43%) had weak staining (IHC score, 1), 39 patients (33.91%) had moderated staining (IHC score, 2), and 13 patients (10.43%) had strong staining (IHC score, 3; Figure 3B). Negative and weak stainings were defined as low STC2 expression, whereas moderate and strong stainings were defined as high STC2 expression. The correlations between STC2 expression level and clinicopathological features of CRC were calculated by the chi-squared test and are summarized in Table 2. As shown in Table 2, high STC2 expression in CRC tissues was significantly correlated with lymph node metastasis (P=0.047), distant metastasis (P=0.040), and clinical stage (P=0.047). However, there were no statistically significant relationships between STC2 expression and other clinicopathological variables, such as gender (P=0.595), age (P=0.932), tumor location (P=0.136), tumor size (P=0.105), serum CEA (P=0.904), and tumor infiltration (P=0.205).

Table 2 Correlations between STC2 expression level and clinicopathological variables of 115 CRC patients

	Number	STC2 expression level		P-value
		Low (n=63)	High (n=52)
Gender				0.595
Male	61	32	29
Female	54	31	23
Age, years				0.932
<60	57	31	26
≥60	58	32	26
Tumor location				0.136
Colon	62	30	32
Rectum	53	33	20
Tumor size, cm				0.105
<5	47	30	17
≥5	68	33	35
Serum CEA, ng/mL				0.904
<5	56	31	25
≥5	59	32	27
Vascular invasion				0.655
Present	33	17	16
Absent	82	46	36
Lymphatic invasion				0.137
Present	34	15	19
Absent	81	48	33
Tumor infiltration				0.205
T1 and T2	38	24	14
T3 and T4	77	39	38
Lymph node metastasis				0.047
Present	46	20	26
Absent	69	43	26
Node counts				0.055
<12	38	16	22
≥12	77	47	30
Distant metastasis				0.040
Present	17	5	12
Absent	98	58	40
Clinical stage				0.047
I	19	14	5
II	45	27	18
III	34	17	17
IV	17	5	12

Note: Bold indicates statistically significant values (P<0.05).

Abbreviations: CEA, carcinoembryonic antigen; CRC, colorectal cancer; STC2, stanniocalcin 2.

High STC2 expression indicated a worse prognosis in colorectal cancer

For all the study subjects, the follow-up period ranged from 1 to 52 months, with a median survival time of 42.8 months. The OS rate was 53.0%, and in patients with high STC2 expression level and those with low level, it was 38.5% and 65.1%, respectively. Additionally, Kaplan–Meier and Cox’s proportional hazards regression model survival analysis revealed that CRC patients with high STC2 expression had a shorter OS (Figure 3C).

To determine whether STC2 expression was an independent prognostic predictor for CRC patients, univariate and multivariate analyses were performed to compare the impact of STC2 expression and other clinicopathological factors on the prognosis. Univariate analysis revealed that clinical variables, including vascular invasion (HR =2.311, 95% CI: 1.302–4.100, P=0.004), lymphatic invasion (HR =1.876, 95% CI: 1.054–3.341, P=0.033), lymph node metastasis (HR =2.391, 95% CI: 1.337–4.275, P=0.003), distant metastasis (HR =5.057, 95% CI: 2.647–9.662, P<0.001), clinical stage (HR =3.015, 95% CI: 1.650–5.508, P<0.001), and STC2 (HR =2.273, 95% CI: 1.278–4.042, P=0.005), were significantly associated with OS of CRC patients (Table 3). Furthermore, in multivariate Cox regression analysis, only distant metastasis (HR =3.742, 95% CI: 1.633–8.575, P=0.002) and STC2 expression (HR =1.976, 95% CI: 1.092–3.576, P=0.024) were independent prognostic factors for CRC (Table 3). The high expression of STC2 had been consistently observed in the CRC patients with poor OS, indicating that STC2 may be functionally important in CRC pathogenesis.

Table 3 Univariate and multivariate analyses of clinicopathological characteristics and STC2 with overall survival in 115 CRC patients

	Univariate analysis		Multivariate analysis
	HR (95% CI)	P-value	HR (95% CI)	P-value
Gender (male vs female)	1.155 (0.659–2.022)	0.615
Age, years (≥median vs <median)	0.749 (0.427–1.313)	0.313
Tumor location (rectum vs colon)	0.840 (0.478–1.475)	0.682
Tumor size, cm (≥5 vs <5)	1.730 (0.949–3.153)	0.074
Serum CEA, ng/mL (≥5 vs <5)	0.842 (0.481–1.474)	0.546
Vascular invasion (present vs absent)	2.311 (1.302–4.100)	0.004	1.875 (0.992–3.545)	0.053
Lymphatic invasion (present vs absent)	1.876 (1.054–3.341)	0.033	1.740 (0.940–3.221)	0.078
Tumor infiltration (T3 and T4 vs T1 and T2)	1.626 (0.861–3.072)	0.112
Lymph node metastasis (present vs absent)	2.391 (1.337–4.275)	0.003	1.303 (0.391–4.339)	0.666
Node counts (≥12 vs <12)	1.364 (0.733–2.537)	0.327
Distant metastasis (present vs absent)	5.057 (2.647–9.662)	<0.001	3.742 (1.633–8.575)	0.002
Clinical stage (III and IV vs I and II)	3.015 (1.650–5.508)	<0.001	1.335 (0.344–5.176)	0.676
STC2 (≥median vs <median)	2.273 (1.278–4.042)	0.005	1.976 (1.092–3.576)	0.024

Note: Bold indicates statistically significant values (P<0.05).

Abbreviations: CEA, carcinoembryonic antigen; CRC, colorectal cancer; STC2, stanniocalcin 2.

Discussion

CRC local recurrence and distant metastasis leads to the death of 50,260 patients annually in the United States in 2017.¹ Therefore, it is critical to distinguish the high risk of recurrence and metastasis to improve the survival of patients with CRC. Based on our previous microarray data, we identified that STC2 could be a potential new biomarker to diagnose and evaluate prognosis of CRC patients in the present study, which represented the association of STC2 expression with clinicopathological features and prognosis in CRC. By analyzing the data from the GEO and TCGA databases, we confirmed that STC2 mRNA level was significantly upregulated in CRC tissues compared with that in normal tissues. The data from our medical center also demonstrated a higher STC2 level in CRC tissue. Similar to the cellular location of STC2 in other types of carcinoma, STC2 protein was located mainly in the cytoplasm and partially in the nucleus. Additionally, STC2 expression was significantly correlated with lymph node metastasis, distant metastasis, and clinical stage. Moreover, high STC2 expression was associated with worse prognosis in patients undergoing surgery for CRC. These findings suggested that STC2 may be involved in CRC progression, and serve as a useful predictor of prognosis.

STC2 belongs to stanniocalcin gene family, which is a family of secreted glycoprotein hormones that was initially discovered in the corpuscles of stannius, an endocrine gland of bony fish.^26,27 STC2 and its homologues STC1, are reported to be involved in calcium and phosphate homeostasis, metabolism, reproduction, stress response, and development.^27,28 The human STC2 gene encoding 302 amino acid long protein is widely expressed in the kidney, heart, pancreas, and spleen.²⁹ Emerging evidence has revealed abnormal expression of STC2 in tumor tissues, which might be used as a new biomarker to evaluate prognosis of patients with malignant tumors.^{16–19,30,31} Previous studies have reported that elevated expression of STC2 in gastric cancer tissues was significantly correlated with depth of tumor invasion, lymph node metastasis, venous invasion, and clinical stage, and predicted a worse outcome.^17,32 Similarly, aberrant STC2 expression was also correlated with advanced tumor progression in esophageal squamous-cell cancer, hepatocellular carcinoma, nasopharyngeal cancer, lung cancer, and ovarian cancer.^{14,18,19,33,34} Ieta et al reported that high STC2 expression in CRC tissue was significantly associated with lymph node metastasis, lymphatic invasion, tumor depth, clinical stage, and worse OS.²² Compared with this study, our study employed 596 CRC samples, 132 from GEO, 349 from TCGA, and 115 from our medical center. We accessed more cases from the public database. Furthermore, we first identified STC2 gene was highly expressed in CRC than in normal tissue and high STC2 in cancer tissue was significantly correlated with poor OS by bioinformatics analysis. Subsequently, data from our medical center further confirmed this result. More cases and secondary validation ensured the reliability of our results. Hence, we proved that STC2 may serve as a useful prognosis biomarker in CRC patients.

The impact of STC2 on malignant biological properties has been preliminarily studied. It is reported that knockdown of STC2 in H460 lung cancer cells suppressed growth, migration, invasion, and G0/G1 cell cycle progression, and attenuated the H₂O₂-induced oxidative stress on cell viability with a subsequent increase in intracellular ROS levels.¹⁴ In ovarian cancer, overexpressed STC2 promoted epithelial– mesenchymal transition, as revealed by the increase of N-cadherin/vimentin but a decrease of E-cadherin levels, and enhanced degree of invasiveness under hypoxic condition, mediated by increasing MMP2 and MMP9 expressions and activating ERK1/2 levels.¹³ Mechanistically, STC2 could trigger the activation of MMP2, Mus81, PI3K/AKT/Snail, and MAPK signaling pathway in neuroblastoma, hepatocellular carcinoma, head and neck squamous cell carcinoma, and cervical cancer.^16,35–37 Interestingly, STC2 acts as an oncogene in many kinds of cancers, but it has been reported to be a polite gene in breast cancer. Studies had reported that STC2 overexpression could suppresses breast cancer cell migration and invasion via the PKC/Claudin-1-mediated signaling, and be associated with favorable prognosis, particularly ER-positive breast cancers.^21,38 ER-positive breast cancers are usually low-grade malignancies and can be effectively treated with hormonotherapy, which may explain the good prognosis.^39,40 Moreover, STC2 could promote growth and migration under hypoxic conditions, and stimulate P-glycoprotein via the PI3K/Akt signaling pathway to influence oxaliplatin resistance in CRC.^41,42 There is little research paper for the molecular mechanism of STC2 on CRC behaviors. Therefore, future studies need to elucidate the molecular mechanisms of STC2 in CRC.

Conclusion

Our results revealed that STC2 is an independent prognostic factor for OS in CRC patients. High STC2 expression is strongly correlated with lymph node metastasis, distant metastasis, advanced clinical stage, and worse clinical outcome. Our findings suggested that STC2 could be a potential new biomarker to diagnose and evaluate prognosis of patients with CRC.

Author contributions

CZ and SC contributed equally to this work. CZ, SC, WX, MF, and BX were responsible for conception and design of the study. CZ, SC, QY, and FS did the experiments, data collection, statistical analyses, and writing of the manuscript. XM and XS participated in data collection and provided statistical expertise. All authors contributed to data analysis, drafting and revising the article, gave final approval of the version to be published, and agree to be accountable for all aspects of the work.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (81072152, 81770283), Natural Science Foundation of Hubei Province (2015CFA027), Research Foundation of Health and Family Planning Commission of Hubei Province (WJ2015MA010, WJ2017M249), and Clinical Medical Research Center of Peritoneal Cancer of Wuhan (2015060911020462).

Supplementary materials

Figure S1 Differentially expressed genes in CRC vs normal tissues across each independent dataset. Each column represents a sample and each row represents the expression level of a given mRNA. The color scale represents the raw Z score ranging from blue (low expression) to red (high expression).

Abbreviation: CRC, colorectal cancer.

Table S1 Microarray datasets of mRNA expression profiles

Dataset	Year	Country	Organization	Platform	Sample		PMID
					Normal tissue	CRC tissue
GSE21510^Citation1	2011	Japan	Juntendo University	GPL570	25	19	21270110
GSE32323^Citation2	2012	Japan	Juntendo University	GPL570	17	17	22399497
GSE39582^Citation3	2013	France	Ligue Nationale Contre le Cancer	GPL570	17	17	23700391
GSE41328^Citation4	2012	USA	Illinois University	GPL570	10	10	17160039

Abbreviation: CRC, colorectal cancer.

References

• TsukamotoSIshikawaTIidaSClinical significance of osteoprotegerin expression in human colorectal cancerClin Cancer Res20111782444245021270110
   PubMed Web of Science ®Google Scholar
• KhamasAIshikawaTShimokawaKScreening for epigenetically masked genes in colorectal cancer Using 5-Aza-2′-deoxycytidine, microarray and gene expression profileCancer Genom Proteom2012926775
   PubMed Web of Science ®Google Scholar
• MarisaLReynièsADDuvalAGene expression classification of colon cancer into molecular subtypes: characterization, validation, and prognostic valuePLoS Med2013105e100145323700391
   PubMed Web of Science ®Google Scholar
• LinGHeXJiHReproducibility probability score – incorporating measurement variability across laboratories for gene selectionNat Biotechnol200624121476147717160039
   PubMed Web of Science ®Google Scholar

Disclosure

The authors report no conflicts of interest in this work.