Skip to main content

Investigating the prognostic and predictive value of the type II cystatin genes in gastric cancer

Abstract

Background

Accumulating evidence indicates that type II cystatin (CST) genes play a pivotal role in several tumor pathological processes, thereby affecting all stages of tumorigenesis and tumor development. However, the prognostic and predictive value of type II CST genes in GC has not yet been investigated.

Methods

The present study evaluated the expression and prognostic value of type II CST genes in GC by using The Cancer Genome Atlas (TCGA) database and the Kaplan–Meier plotter (KM plotter) online database. The type II CST genes related to the prognosis of GC were then screened out. We then validated the expression and prognostic value of these genes by immunohistochemistry. We also used Database for Annotation, Visualization, and Integrated Discovery (DAVID), Gene Multiple Association Network Integration Algorithm (GeneMANIA), Search Tool for the Retrieval of Interacting Genes/Proteins (STRING), nomogram, genome-wide co-expression analysis, and other bioinformatics tools to analyze the value of type II CST genes in GC and the underlying mechanism.

Results

The data from the TCGA database and the KM plotter online database showed that high expression of CST2 and CST4 was associated with the overall survival (OS) of patients with GC. The immunohistochemical expression analysis showed that patients with high expression of CST4 in GC tissues have a shorter OS than those with low expression of CST4 (HR = 1.85,95%CI: 1.13–3.03, P = 0.015). Multivariate Cox regression analysis confirmed that the high expression level of CST4 was an independent prognostic risk factor for OS.

Conclusions

Our findings suggest that CST4 could serve as a tumor marker that affects the prognosis of GC and could be considered as a potential therapeutic target for GC.

Peer Review reports

Introduction

Globally, gastric cancer (GC) is the fifth most common malignant tumor in new incidence and the fourth leading cause of cancer-related death worldwide [1]. The incidence of GC is highly region-specific, with the highest incidence in East Asia, including China [2]. The 2015 statistics from the National Cancer Center in China showed that there were 679,000 new cases of GC and 498,000 deaths from GC in China, and the new incidence of morbidity and mortality due to GC ranked the second among all malignant tumors in China during the same period [3]. According to statistics, the 5-year overall survival (OS) rate of patients with GC in China in 2015 was 35.1%, while the 5-year OS rate of patients with GC in the United States in 2014 was only 33.1% [4, 5]. The reason for this low survival rate is that only a small number of patients with GC are diagnosed at an early stage, and most patients with GC are diagnosed at an advanced stage; moreover, there are only few effective treatment approaches for patients in advanced stage [6]. Therefore, it is critical to explore new molecular biomarkers for the diagnosis, treatment, and prognosis of GC.

The cystatin (CST) superfamily are a kind of thiol proteinase inhibitors (TPIs) that are widely present in human tissues and body fluids and function as competitive and reversible inhibitors of cysteine proteases. The CST superfamily is divided into three subfamilies according to their three-dimensional molecular structure and biochemical function [7]: (1) type I CSTs are also known as stefins and include stefin A and stefin B [8]; (2) type II CSTs are extracellular proteins; their members include CST C, CST D, CST E/M, CST F, CST S, CST SN, and CST SA, and their coding genes are CST3, CST5, CST6, CST7, CST4, CST1, and CST2, respectively [9]; and (3) type III CSTs, also known as kininogens, are multidomain proteins. There are three types of kininogens: high-molecular-weight kininogens, low-molecular-weight kininogens, and T-type kininogens [10]. Several studies have shown that type II CSTs play a pivotal role in many pathological processes and affect various stages of tumorigenesis and tumor development, including proliferation, apoptosis, invasion, metastasis, and angiogenesis [11,12,13]. To date, no relevant research studies have been conducted on the expression and prognostic value of type II CST genes in GC.

To determine the prognostic value and potential function of type II CST genes in GC, we analyzed the expression and prognostic value of type II CST genes by using The Cancer Genome Atlas (TCGA; https://portal.gdc.cancer.gov) database and the Kaplan–Meier plotter (KM plotter) online database. We also used immunohistochemistry to validate the expression of type II CST genes in GC tissues and their prognostic value. Finally, various bioinformatics tools were used to determine the potential function and mechanism of type II CST genes.

Materials and methods

Functional and co‑expression analyses

The Database for Annotation, Visualization, and Integrated Discovery (DAVID, v.6.8; https://david.ncifcrf.gov; accessed September 1, 2020) was used to analyze Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment of type II CST genes [14, 15]. The Gene Multiple Association Network Integration Algorithm (GeneMANIA v.3.6.0; http://www.genemania.org/; accessed May 20, 2019) was used to predict the function of type II CST genes [16]. The Search Tool for the Retrieval of Interacting Genes (STRING v.11.0; https://string-db.org/; accessed May 20, 2019) was used to detect protein–protein interaction (PPI) networks [17].

TCGA

We downloaded mRNA expression data and the corresponding clinical information related to GC from the TCGA (https://portal.gdc.cancer.gov; accessed August 22, 2018). The TCGA database included 407 patients diagnosed with GC, including 375 GC tissues and 32 adjacent normal tissues. After excluding cases with missing follow-up data, 351 patients with GC were included in our present analysis. We used the TCGA database to perform paired t-test on the mRNA expression levels of type II CST genes in GC and adjacent normal tissues, and scatter plots and receiver operating characteristic (ROC) curves were then plotted. For survival analysis, patients with GC were divided into two groups on the basis of the median expression level of type II CST genes: high expression group and low expression group. The mRNA expression of type II CST genes and clinical information in the TCGA database were used for the univariate survival analysis and the multivariate Cox regression model survival analysis.

KM plotter online database

The KM plotter online database (http://kmplot.com/analysis/index.php?p=service&cancer=gastric; accessed September 27, 2020) was used to analyze the relationship between the mRNA expression of type II CST genes and the OS of patients with GC. The KM plotter online database was established using the gene expression data and survival information of 875 GC patients; these data were downloaded from the Gene Expression Omnibus (GEO) (GSE14210, GSE15459, GSE22377, GSE29272, GSE51105, and GSE62254) [18]. The KM plotter online database is used for the discovery and verification of biomarkers. The members of the type II CST family were entered into the KM plotter online database for analysis. All GC patients were divided into two groups on the basis of the median mRNA expression level of type II CST genes: high expression group and low expression group. Finally, parameters such as survival analysis graph, hazard ratio (HR), 95% confidence interval (CI), and log-rank P were analyzed.

Nomogram

A nomogram was constructed, and the prognostic risk between type II CST genes and the OS grade in patients with GC was assessed. The possible utility of type II CST genes to predict clinical rank was also assessed. According to prognostic clinical indicators and Cox regression model survival analysis, related factors and type II CST gene expression levels were included in the risk model. The scores for each factor were counted, and the 1-, 3-, 5-, and 10-year survival rates were calculated [19]. The nomogram was constructed using the rms package (https://CRAN.R-project.org/package=rms), along with its dependencies in the R platform, and visualized using gplot [20].

Stratified analyses

The stratified analysis method was used to evaluate the prognostic value of type II CST genes in different stratifications in the TCGA and the KM plotter online databases.

Immunohistochemistry

GC sections with complete clinicopathological and prognostic data were collected from the Affiliated Hospital of Guangxi Medical University. The sections were placed in 65◦C oven for 2 h, dewaxed with xylene, hydrated with a graded series of ethanol, and repaired with antigen repair buffer by EDTA method. Endogenous antigens were blocked with 3% hydrogen peroxide. The sections were incubated with CST2 antibody (Thermo Fisher Scientific, 1:150 dilution) and CST4 antibody (ABclonal, 1:200 dilution) overnight at 4◦C. The sections were heated for 30 min to room temperature and incubated with secondary antibody for 20 min. The staining result was visualized using DAB color developing solution. The sections were stained with Hematoxylin, differentiated in 1% hydrochloric acid in alcohol, and dehydrated. Thereafter, the sections were naturally dried in a fume hood, transparentized by xylene, and sealed with neutral resin. The negative control used PBS instead of primary antibody, and the positive section was served as positive control. Immunohistochemistry was performed to determine the expression of type II CST genes in GC and adjacent normal tissues. Paired t-test was used to analyze the difference in the expression level between GC and adjacent normal tissues, and scatter plots and ROC curves were drawn. The X-tile software was used to determine the optimal cutoff value. X-tile, a statistical software developed by Yale University with a single function, can determine the cut-off point of continuous variables and facilitate the drawing of the Kaplan–Meier curve [21]. Finally, the univariate survival analysis and multivariate Cox regression model survival analysis of immunohistochemical expression were performed to analyze the prognosis of patients with GC who were positive or negative for type II CST genes.

Genome-wide co-expression analysis and functional enrichment of type II CST prognostic genes

Pearson’s correlation analysis was performed between type II CST prognostic genes and TCGA genome expression profile data, and the co-expressed genes were screened out. Pearson’s correlation coefficient was calculated by the cor function on the R 3.4.4 platform. The difference was considered to be significant at |r|> 0.6, P < 0.05. These co-expressed genes were then used to perform GO and KEGG functional enrichment analyses by using DAVID v6.8 [22,23,24].

Statistical analyses

Statistical analysis was performed using SPSS v.25.0 software (IBM Corp., Armonk, NY, USA) and GraphPad Prism software (GraphPad Software, Inc., La Jolla, CA, USA) to generate vertical scatterplots and survival curves. Nomograms and correlations were also generated using R software (Vienna, Austria).

Results

Functional enrichment of type II CST genes and co-expression analysis

The biological function of type II CST genes was assessed by DAVID, and biological process (BP), cellular component (CC), and molecular functionality (MF) were obtained from GO and KEGG pathway analyses (Fig. 1A and B). GO analysis showed that the type II CST gene family was enriched in protease binding, negative regulation of proteolysis, negative regulation of cysteine-type endopeptidase activity, extracellular space, extracellular exosome, endopeptidase inhibitor activity, and cysteine-type endopeptidase inhibitor activity (Fig. 1A). KEGG pathway analysis showed that the type II CST gene family was enriched in salivary secretion (Fig. 1B). The gene–gene interaction network of type II CST genes is shown in Fig. 2A; the figure shows that the members of the type II CST gene family co-express and interact with other genes. According to the information in the STRING database, the PPI interaction network showed direct or indirect connections between members of the type II CST gene family, among which CST7 was directly related to other family members (Fig. 2B).

Fig. 1
figure 1

GO and KEGG analysis. (A)GO analysis of Type II CST genes. (B)KEGG analysis of Type II CST genes. CST, Cystatins; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes

Fig. 2
figure 2

Gene and protein interaction networks. A Type II CST genes multiple association network integration algorithm. The size of the circle represents the strength of the co‑expression. B Protein‑protein interaction networks. CST, Cystatins

Survival analysis in the TCGA

From the vertical scatter plot of type II CST gene expression levels in the TCGA, it was observed that the expression levels of CST1, CST2, CST3, and CST4 were significantly different in GC tissues and adjacent normal tissues (P < 0.05, Fig. 3A). The ROC curve of the expression level of type II CST genes in the TCGA showed that the expression levels of the genes CST1, CST2, CST3, CST4, and CST5 were significantly different (P < 0.001, Fig. 3B). The vertical scatter plot of the expression levels of type II CST genes in the TCGA revealed a significant difference between the high and low expression of all type II CST genes (P < 0.001, Fig. 3C). Figure 3D shows the results of the univariate survival analysis of the expression levels of type II CST genes in the TCGA. The high expression level of CST2, CST4, CST5, and CST6 in patients with GC was significantly correlated with poor prognosis (P = 0.018, HR = 1.48, 95% CI = 1.07–2.06; P = 0.018, HR = 1.48, 95% CI = 1.07–2.06; P = 0.010, HR = 1.53, 95% CI = 1.10–2.13; P = 0.030, HR = 1.43, 95% CI = 1.03–1.99, respectively).

Fig. 3
figure 3

Survival analysis in The Cancer Genome Atlas. A Scatterplots for Type II CST genes expression levels in The Cancer Genome Atlas. B ROC curve of the expression level of type II CST genes in The Cancer Genome Atlas. C Scatterplots for Type II CST genes expression levels in The Cancer Genome Atlas. Red, high expression; blue, low expression. D Prognostic graphs illustrating the impact of Type II CST genes expression on overall survival in The Cancer Genome Atlas. Kaplan‑Meier survival curves for patients with gastric cancer according to median expression of CST1-7 (n = 351). CST, Cystatins; HR, hazard ratio (95% CI)

Univariate survival analysis was performed to explore the prognostic value of clinical factors in patients with GC, including gender, age, tumor location, Helicobacter pylori infection, histological type, histologic grade, microsatellite stability (MMS), tumor stage, cancer status, residual tumor, chemotherapy, and targeted therapy. The results showed that age, tumor stage, cancer status, residual tumor, chemotherapy, and targeted therapy were associated with OS (all P < 0.05, Supplemental Table S1). After including age, tumor stage, cancer status, residual tumor, chemotherapy, and targeted therapy as adjustment factors in the Cox proportional hazards regression models, the results of the multivariate survival analysis showed that high expression of the genes CST2, CST4, and CST5 significantly decreased OS (adjusted P = 0.030, HR = 1.57, 95% CI = 1.05–2.35; adjusted P = 0.010, HR = 1.69, 95% CI = 1.13–2.51; adjusted P = 0.001, HR = 1.94, 95% CI = 1.29–2.91, respectively; Table 1). This result was consistent with that of the univariate survival analysis (P = 0.018, HR = 1.48, 95% CI = 1.07–2.06; P = 0.018, HR = 1.48, 95% CI = 1.07–2.06; P = 0.010, HR = 1.53, 95% CI = 1.10–2.13, respectively; Table 1; Fig. 3D). However, the expression level of the gene CST6 did not show a significant difference in the multivariate survival analysis (P = 0.515, HR = 1.15, 95% CI = 0.76–1.72; Table 1), which was inconsistent with the results of the univariate survival analysis (P = 0.030, HR = 1.43, 95% CI = 1.03–1.99; Table 1; Fig. 3D).

Table 1 Analysis of the association between CST genes and the risk of death in The Cancer Genome Atlas gastric cancer cohort (n = 351)

Validation of the CST cohort using the KM plotter online database

The expression and prognostic value of type II CST genes were validated using the KM plotter online database. The KM curves of type II CST genes are shown in Fig. 4A. The analysis results showed that GC patients with high mRNA levels of CST1, CST2, CST3, CST4, CST5, and CST6 had significantly worse OS (P = 5.8e-06, HR = 1.48, 95% CI = 1.25–1.75; P = 3.6e-09, HR = 1.67, 95% CI = 1.40–1.98; P = 9.6e-05, HR = 1.40, 95% CI = 1.18–1.66; P = 3e-06, HR = 1.50, 95% CI = 1.26–1.78; P = 7.9e-06, HR = 1.47, 95% CI = 1.24–1.74; P < 0.001, HR = 1.40, 95% CI = 1.18–1.66, respectively). Consistent with the TCGA analysis, the expression of CST2, CST4, and CST5 was significantly associated with worse prognosis (adjusted P = 0.030, HR = 1.57, 95% CI = 1.05–2.35; adjusted P = 0.010, HR = 1.69, 95% CI = 1.13–2.51; adjusted P = 0.001, HR = 1.94, 95% CI = 1.29–2.91, respectively; Table 1). Therefore, we screened out CST2 and CST4 for further analyses.

Fig. 4
figure 4

Survival analysis in the Kaplan-Meier plotter online database and Nomogram analysis. A Prognostic graphs of Type II CST median expression for overall survival generated using the Kaplan-Meier plotter online database. B Nomogram for the association between clinicopathological data and risk score. CST, Cystatins; HR, hazard ratio (95% CI)

Nomogram and stratified analysis

Nomograms were drawn based on CST2 and CST4 expression levels, age, tumor stage, cancer status, residual tumor, and targeted therapy. Next, the 1-, 3-, 5-, and 10-year survival rates were calculated (Fig. 4B). The results of nomogram analysis revealed that CST2 and CST4 had certain predictive effects on the prognosis of patients with GC. In particular, the predictive effect of CST4 was higher than that of age and targeted therapy.

Stratified analysis was used to analyze the prognostic value of the CST2 and CST4 genes in GC, and the relationship between CST2 and CST4 and OS was analyzed using the TCGA database (Tables 2 and 3) and the KM plotter online database (Tables 4 and 5). As shown in Table 2, high levels of CST2 mRNA were significantly associated with poor prognosis in the histological signet ring type, tumor stage IV, residual tumor R0, radiotherapy, and target therapy in the TCGA. As noted in Table 3, high levels of CST4 mRNA were significantly associated with poor prognosis in female, age over 65 years, histologic grade G2 and G3, cancer status as with tumor, residual tumor R0 and R1, radiotherapy, and target therapy in the TCGA. As shown in Table 4, high CST2 mRNA levels showed significantly worse prognosis in male and female, tumor stage III and IV, Lauren classification as intestinal and diffuse, poorly differentiated tumor, surgery alone as treatment, and positive and negative HER2 status in the KM plotter online database. As noted in Table 5, high CST4 mRNA levels showed significantly worse prognosis in male and female, tumor stage III and IV, Lauren classification as intestinal and diffuse, well-differentiated tumor, surgery alone as treatment, and positive and negative HER2 status in the KM plotter online database.

Table 2 Stratified analysis of association between CST2 and overall survival in TCGA GC cohort
Table 3 Stratified analysis of association between CST4 and overall survival in TCGA GC cohort
Table 4 Stratified analysis of association between CST2 and overall survival in KM plotter GC cohort
Table 5 Stratified analysis of association between CST4 and overall survival in KM plotter GC cohort

Immunohistochemistry

A total of 115 cases of GC tissues and 20 cases of adjacent tissues with complete clinical data were collected from the Affiliated Hospital of Guangxi Medical University. The demographic and clinical data are shown in Table 6. The results showed that among the clinicopathological parameters of GC in our hospital, tumor stage is associated with patient survival.

Table 6 Demographic and clinical data for 115 GC patients

Figure 5A shows the immunohistochemical semi-quantitative score scatter plot of 20 cases of GC and the corresponding adjacent normal tissues. The expression of CST2 and CST4 was higher in GC tissues than in adjacent normal tissues, and the difference of CST4 was statistically significant (P < 0.05). The ROC curve of the expression levels of CST2 and CST4 is shown in Fig. 5B, and the difference was statistically significant (P < 0.05). Figure 5C shows the immunohistochemical staining images of CST2 and CST4. In the univariate survival analysis, patients with high CST4 immunohistochemical expression in GC tissues had a shorter OS than those with low expression (P = 0.010, HR = 1.78, 95% CI = 1.13–2.81, Fig. 5D). The results in Table 6 showed that among the clinicopathological parameters of GC in our immunohistochemistry experiments, tumor stage is associated with OS. Therefore, after adjusting for the tumor stage in the multivariate Cox model survival analysis, the survival time of GC patients with high CST4 immunohistochemical expression was significantly different from that of GC patients with low CST4 immunohistochemical expression (P = 0.015, HR = 1.76, 95% CI = 1.76–2.78) (Table 7). In both univariate survival analysis and multivariate Cox model survival analysis, no significant difference was observed in survival time of GC patients with high CST2 immunohistochemical expression in our hospital (P > 0.05, Table 7; Fig. 5D).

Fig. 5
figure 5

Immunohistochemical analysis. A The scatter plot of CST2, CST4 immunohistochemical staining semi-quantitative score of gastric cancer and adjacent normal tissues in affiliated hospital of guangxi medical university. B ROC curve of the expression level of CST2, CST4 in affiliated hospital of guangxi medical university. C The IHC staining for CST2 and CST4 in adjacent normal tissues and gastric cancer (IHC, ×200). D Survival analysis of CST2,CST4 immunohistochemical expression in gastric cancer in affiliated hospital of guangxi medical university. CST, Cystatins

Table 7 Analysis of the association between CST2,CST4 and the risk of death in affiliated hospital of guangxi medical university (n = 115)

Genome-wide co-expression analysis and functional enrichment of type II CST prognostic genes

To determine the underlying mechanism by which CST4 plays a role in GC, we performed a genome-wide co-expression analysis. The regulatory network of CST4 and its co-expressed related genes in GC tissues from the TCGA cohort are shown in Supplemental Fig. S1, Supplemental Data S1. The results showed that 561 genes had positive correlations with CST4 in GC.

The results of GO analysis of CST4 and its co-expressed genes of GC in the TCGA showed that it enriched in focal adhesion, negative regulation of cell proliferation, fascia adherens, cell-matrix adhesion, positive regulation of cell-substrate adhesion, negative regulation of cell migration, etc. (Table 8, Supplemental Data S2). Similarly, KEGG functional enrichment analysis showed that CST4 was enriched in the cGMP-PKG signaling pathway, Calcium signaling pathway, cAMP signaling pathway, and MAPK signaling pathway (Table 9).

Table 8 GO term enrichments of co-expressed genes of cystatin 4 in gastric cancer
Table 9 Kyoto Encyclopedia of Genes and Genomes term enrichments of co-expressed genes of cystatin 4 in gastric cancer

Discussion

Cancer development is a multifactorial long-term interaction process, and normal cells undergo multiple steps to transform into malignant cells [25]. Many stages of cancer growth and progression are associated with dysregulation of protease activity [26, 27]. Cysteine proteases are a class of intracellular proteins with a protein-degrading activity that affect a variety of biological processes, including inflammation, immune responses, and cancer [28, 29]. Changes in the expression of cysteine proteases are regulated by their endogenous inhibitor CST [30]. CST can effectively inhibit the activity of cysteine proteases by binding to their active site to form tight equimolar complexes.

The maintenance of the balance between CST and cysteine proteases is crucial to the normal function of biological systems. Once the balance is disrupted, it will cause cell death, leading to tumorigenesis of malignancies [31]. CSTs were originally considered as TPIs; however, this view has changed considerably over the past few decades with various studies targeting the value of CST in tumors. Recently, researchers have revealed that CST affects all stages of cancer progression, including tumor growth, apoptosis, invasion, metastasis, angiogenesis, and antitumor immune responses [32]. In particular, type II CSTs are involved in the metastasis of several types of malignant tumors [12, 33]. However, data regarding the relationship between the expression of type II CST genes and the prognosis of patients with GC are unavailable.

The combination of serum CST4 and DR-70 was reported to further enhance the diagnosis of early colorectal cancer [34]. Recent studies also indicated CST4 was a novel and improved diagnostic marker for colorectal cancer [35]. In the present study, the TCGA database was used to analyze the expression of type II CST genes in GC tissues and adjacent normal tissues. The results showed that the mRNA expression levels of CST1, CST2, and CST4 were higher in GC tissues than in adjacent normal tissues, while the mRNA expression level of CST3 was lower in GC tissues than in adjacent normal tissues. Our immunohistochemical analysis based on the tissues collected from the Affiliated Hospital of Guangxi Medical University also showed that the expression of the CST4 protein was higher in GC tissues than in the adjacent tissues. Dou et al. showed that the mRNA and protein expression of CST4 was significantly upregulated in GC tissues and cell lines, and its sensitivity and specificity for diagnosing GC were 69.0% and 85.6%, respectively [36]. The study of Zhang et al. [37] also supports this conclusion. These studies suggest that CST4 may be a potential novel tumor marker for GC.

Wang et al. [38] indicated that CST4 may be useful in predicting the prognosis of ovarian cancer. In our present study, both univariate survival analysis and multivariate Cox regression survival analysis in the TCGA database showed that the OS of GC patients with high expression of CST2, CST4, and CST5 mRNA was shorter than that of GC patients with low expression of these mRNAs. In the KM plotter online database, the results of survival analysis showed that GC patients with high expression of CST1, CST2, CST3, CST4, CST5, and CST6 had a shorter OS than those with low expression of these genes. The combined results of the TCGA database and the KM plotter online database indicated that CST2 and CST4 may play a role in the prognosis of GC. Accordingly, we screened out CST2 and CST4 for further immunohistochemical experimental validation. In the univariate survival analysis and multivariate Cox model in our study, GC patients with high CST4 immunohistochemical expression had a shorter OS than those with low CST4 expression; this finding was consistent with the results obtained for the TCGA database and the KM plotter online database. In both univariate survival analysis and multivariate Cox model survival analysis, no significant difference was observed in the survival time of patients with high CST2 immunohistochemical expression. According to the results of the nomogram study, the combination of CST2, CST4, and other markers may be considered as a new method to predict the prognosis of GC, and the predictive effect of CST4 was greater than that of age and targeted therapy. Stratified analysis also suggested that CST2 and CST4 have value as an independent prognostic indicator for GC. Taken together, we concluded that GC patients with high expression of CST4 show poor prognosis.

The mechanism by which CST4 plays a role in the occurrence and development of GC has not yet been elucidated. Zhang et al. [37] found that CST4 overexpression promoted the proliferation, migration, and invasion of GC cells in vitro and significantly promoted the tumorigenicity of GC cells in vivo, while silencing CST4 yielded the opposite results. Moreover, silencing CST4 significantly inhibited the migration and invasion ability of GC cell lines and lung metastasis in vivo, while overexpression of CST4 had the opposite effect. It is speculated that the mechanism may be related to the enhancement of the invasiveness of GC by CST4 by regulating the downstream target fibronectin type III domain containing 2 (ELFN2) signaling pathway. To further study the mechanism of CST4 in the development of GC, we used bioinformatics tools, namely genome-wide co-expression analysis, GO term analysis, and KEGG functional enrichment analysis. The genome-wide co-expression analysis of CST4 in GC indicated that CST4 had a co-expression relationship with several genes. The results of GO term analysis showed that CST4 was involved in focal adhesion, negative regulation of cell proliferation, fascia adherens, cell-matrix adhesion, positive regulation of cell-substrate adhesion, and negative regulation of cell migration. Abnormal expression of cell adhesion molecules is often associated with general carcinogenesis and may lead to abnormal proliferation of normal cells [39]. In the present study, the KEGG functional enrichment analysis revealed that CST4 was involved in the regulation of the cGMP-PKG signaling pathway, calcium signaling pathway, cAMP signaling pathway, and MAPK signaling pathway. Activation of the cGMP/PKG signaling pathway promotes the development and progression of GC [40, 41]. The calcium signaling pathway is related to the process of angiogenesis and proliferation of cancer cells [42, 43]. The spatial regulation of cAMP is critical for accurate signal coding [44]. The MAPK signaling pathway plays an important role in controlling cellular processes, including proliferation, differentiation, and apoptosis [45].

The present study has one potential limitation. We used clinical specimens to validate the results obtained by bioinformatics, which largely confirmed the differential expression of CST4 and its significance in the survival and prognosis of patients with GC. However, we were unable to conduct in vivo and in vitro experiments to further confirm the mechanism by which CST4 affects GC occurrence and development. This aspect will be the focus of our future studies.

Conclusion

In conclusion, this study systematically investigated the relationship between the expression of type II CST genes and OS in patients with GC. As the findings suggest, high expression of CST4 in GC is correlated with a poor prognosis. Therefore, CST4 may be a novel prognostic and predictive indicator for GC.

Availability of data and materials

The datasets generated and analysed during the current study are available in the KM plotter online database (http://kmplot.com/analysis/index.php?p=service&cancer=gastric) and TCGA (https://portal.gdc.cancer.gov). The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

References

  1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209-49.

  2. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  PubMed  Google Scholar 

  3. Chen W, Zheng R, Baade PD, Zhang S, Zeng H, Bray F, Jemal A, Yu XQ, He J. Cancer statistics in China, 2015. CA Cancer J Clin. 2016;66(2):115–32.

    Article  PubMed  Google Scholar 

  4. Zeng H, Chen W, Zheng R, Zhang S, Ji JS, Zou X, Xia C, Sun K, Yang Z, Li H, et al. Changing cancer survival in China during 2003–15: a pooled analysis of 17 population-based cancer registries. Lancet Glob Health. 2018;6(5):e555–67.

    Article  PubMed  Google Scholar 

  5. Allemani C, Matsuda T, Di Carlo V, Harewood R, Matz M, Nikšić M, Bonaventure A, Valkov M, Johnson CJ, Estève J, et al. Global surveillance of trends in cancer survival 2000–14 (CONCORD-3): analysis of individual records for 37 513 025 patients diagnosed with one of 18 cancers from 322 population-based registries in 71 countries. Lancet. 2018;391(10125):1023–75.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Lordick F, Kang Y-K, Chung H-C, Salman P, Oh SC, Bodoky G, Kurteva G, Volovat C, Moiseyenko VM, Gorbunova V, et al. Capecitabine and cisplatin with or without cetuximab for patients with previously untreated advanced gastric cancer (EXPAND): a randomised, open-label phase 3 trial. Lancet Oncol. 2013;14(6):490–9.

    Article  CAS  PubMed  Google Scholar 

  7. Dickinson DP. Salivary (SD-type) cystatins: over one billion years in the making–but to what purpose? Crit Rev Oral Biol Med. 2002;13(6):485–508.

    Article  CAS  PubMed  Google Scholar 

  8. Cox JL. Cystatins and cancer. Front Biosci (Landmark Ed). 2009;14(2):463–74.

    Article  CAS  PubMed  Google Scholar 

  9. Manconi B, Liori B, Cabras T, Vincenzoni F, Iavarone F, Castagnola M, Messana I, Olianas A. Salivary cystatins: exploring new post-translational modifications and polymorphisms by top-down high-resolution mass spectrometry. J Proteome Res. 2017;16(11):4196–207.

    Article  CAS  PubMed  Google Scholar 

  10. Kashuba E, Bailey J, Allsup D, Cawkwell L. The kinin-kallikrein system: physiological roles, pathophysiology and its relationship to cancer biomarkers. Biomarkers. 2013;18(4):279–96.

    Article  CAS  PubMed  Google Scholar 

  11. Shamsi A, Bano B. Journey of cystatins from being mere thiol protease inhibitors to at heart of many pathological conditions. Int J Biol Macromol. 2017;102:674–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Dai DN, Li Y, Chen B, Du Y, Li SB, Lu SX, Zhao ZP, Zhou AJ, Xue N, Xia TL, et al. Elevated expression of CST1 promotes breast cancer progression and predicts a poor prognosis. J Mol Med (Berl). 2017;95(8):873–86.

    Article  CAS  PubMed  Google Scholar 

  13. Oh BM, Lee SJ, Cho HJ, Park YS, Kim JT, Yoon SR, Lee SC, Lim JS, Kim BY, Choe YK, et al. Cystatin SN inhibits auranofin-induced cell death by autophagic induction and ROS regulation via glutathione reductase activity in colorectal cancer. Cell Death Dis. 2017;8(3):e2682.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):P3.

    Article  PubMed  Google Scholar 

  15. da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  PubMed  Google Scholar 

  16. Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(Web Server issue):W214-220.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, Santos A, Doncheva NT, Roth A, Bork P, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):D362–8.

    Article  CAS  PubMed  Google Scholar 

  18. Szász AM, Lánczky A, Nagy Á. Cross-validation of survival associated biomarkers in gastric cancer using transcriptomic data of 1,065 patients. Oncotarget. 2016;7(31):49322–33.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Balachandran VP, Gonen M, Smith JJ, DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol. 2015;16(4):e173–80.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Liao X, Yang C, Huang R, Han C, Yu T, Huang K, Liu X, Yu L, Zhu G, Su H, et al. Identification of potential prognostic long non-coding RNA biomarkers for predicting survival in patients with hepatocellular carcinoma. Cell Physiol Biochem. 2018;48(5):1854–69.

    Article  CAS  PubMed  Google Scholar 

  21. Camp RL, Dolled-Filhart M, Rimm DL. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. 2004;10(21):7252–9.

    Article  CAS  PubMed  Google Scholar 

  22. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–92.

    Article  CAS  PubMed  Google Scholar 

  25. Hanahan D, Weinberg RA. Hallmarks of cancer: the next generation. Cell. 2011;144(5):646–74.

    Article  CAS  PubMed  Google Scholar 

  26. Breznik B, Motaln H, Lah Turnsek T. Proteases and cytokines as mediators of interactions between cancer and stromal cells in tumours. Biol Chem. 2017;398(7):709–19.

    Article  CAS  PubMed  Google Scholar 

  27. Olson OC, Joyce JA. Cysteine cathepsin proteases: regulators of cancer progression and therapeutic response. Nat Rev Cancer. 2015;15(12):712–29.

    Article  CAS  PubMed  Google Scholar 

  28. Travis J, Potempa J. Bacterial proteinases as targets for the development of second-generation antibiotics. Biochim Biophys Acta - Protein Struct Mol Enzymol. 2000;1477:35–50.

    Article  CAS  Google Scholar 

  29. Koblinski JE, Ahram M, Sloane BF. Unraveling the role of proteases in cancer. Clin Chim Acta. 2000;291(2):113–35.

    Article  CAS  PubMed  Google Scholar 

  30. Turk V, Stoka V, Turk D. Cystatins, biochemical and structural properties, and medical relevance. Bioscience. 2008;13:5406–20.

    CAS  Google Scholar 

  31. Gole B, Huszthy PC, Popovic M, Jeruc J, Ardebili YS, Bjerkvig R, Lah TT. The regulation of cysteine cathepsins and cystatins in human gliomas. Int J Cancer. 2012;131(8):1779–89.

    Article  CAS  PubMed  Google Scholar 

  32. Breznik B, Mitrovic A, Lah TT, Kos J. Cystatins in cancer progression: more than just cathepsin inhibitors. Biochimie. 2019;166:233–50.

    Article  CAS  PubMed  Google Scholar 

  33. Shi D, Zhou Z, Zhang S. miRNA-6715-5p inhibits cellular proliferation and invasion in colorectal cancer by directly targeting CST4. J Oncol. 2021;2021:7615712.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Cai L, Tu M, Yin X, Zhang S, Zhuang W, Xia Y, Zhang Y, Zhang L, Yu L, Chi L, et al. Combination of serum CST4 and DR-70 contributes to early diagnosis of colorectal cancer. Clin Chim Acta. 2022;531:318–24.

    Article  CAS  PubMed  Google Scholar 

  35. Huang M, Yang Z, Ren J, Wang T, Chen D, Zhan Y, Le W, Deng S. The diagnosis significance of serum cysteine protease inhibitors (CST4) in colorectal cancer. Technol Cancer Res Treat. 2023;22:15330338231164232.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Dou Y, Lv Y, Zhou X, He L, Liu L, Li P, Sun Y, Wang M, Gao M, Wang C. Antibody-sandwich ELISA analysis of a novel blood biomarker of CST4 in gastrointestinal cancers. Onco Targets Ther. 2018;11:1743–56.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Zhang YQ, Zhang JJ, Song HJ, Li DW. Overexpression of CST4 promotes gastric cancer aggressiveness by activating the ELFN2 signaling pathway. Am J Cancer Res. 2017;7(11):2290–304.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Wang S, Wang C, Liu O, Hu Y, Li X, Lin B. Prognostic value of immune-related cells and genes in the tumor microenvironment of ovarian cancer, especially CST4. Life Sci. 2021;277:119461.

    Article  CAS  PubMed  Google Scholar 

  39. Janiszewska M, Primi MC, Izard T. Cell adhesion in cancer: beyond the migration of single cells. J Biol Chem. 2020;295(8):2495–505.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Xiang T, Yuan C, Guo X, Wang H, Cai Q, Xiang Y, Luo W, Liu G. The novel ZEB1-upregulated protein PRTG induced by Helicobacter pylori infection promotes gastric carcinogenesis through the cGMP/PKG signaling pathway. Cell Death Dis. 2021;12(2):150.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Piazza GA, Ward A, Chen X, Maxuitenko Y, Coley A, Aboelella NS, Buchsbaum DJ, Boyd MR, Keeton AB, Zhou G. PDE5 and PDE10 inhibition activates cGMP/PKG signaling to block Wnt/beta-catenin transcription, cancer cell growth, and tumor immunity. Drug Discov Today. 2020;25(8):1521–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Cui C, Merritt R, Fu L, Pan Z. Targeting calcium signaling in cancer therapy. Acta Pharm Sin B. 2017;7(1):3–17.

    Article  PubMed  Google Scholar 

  43. Li Q, Xu D, Gu Z, Li T, Huang P, Ren L. Rutin restrains the growth and metastasis of mouse breast cancer cells by regulating the microRNA-129-1-3p-mediated calcium signaling pathway. J Biochem Mol Toxicol. 2021;35(7):e22794.

    Article  CAS  PubMed  Google Scholar 

  44. Zaccolo M, Zerio A, Lobo MJ. Subcellular organization of the cAMP signaling pathway. Pharmacol Rev. 2021;73(1):278–309.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Kim EK, Choi EJ. Compromised MAPK signaling in human diseases: an update. Arch Toxicol. 2015;89(6):867–82.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors thank the contribution of TCGA (https://cancergenome.nih.gov/) database and the KM plotter (http://kmplot.com/analysis/index.php?p=service&cancer=gastric) for sharing the GC data on open access.

Funding

The study was supported by the Natural Science Foundation of Guangxi Zhuang Autonomous Region (NO. 2017AB45153); the Scientific Research and Technology Development Program of Guangxi (NO.1598011-4); the Natural Science Foundation of Guangxi (NO.2016GXNSFAA380180); the Guangxi Zhuang Autonomous Region Health and Family Planning Commission (NO. Z2015526 and NO. Z-K20221813), the National Natural Science Foundation of China (NO.82260465) and the Youth Science Foundation of Guangxi Medical University (NO. GXMUYSF201502).

Author information

Authors and Affiliations

Authors

Contributions

Jun-Qiang Chen, Ye-yang Chen and Bo-pei Li conceived, designed the research and drafted the manuscript; Jun-fu Wang were responsible for the acquisition of data; Ye Wang analyzed and interpreted the data; Shan-shan Luo performed the statistical analysis; Ru-jing Lin and Xi-wen Liao revised the manuscript for important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Jun-Qiang Chen.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of the First Affiliated Hospital of Guangxi Medical University (Approval Number: 2020 (KY-E-079)). Written informed consent was obtained from the patients to participating in this study. The study was performed under the principles of the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Supplemental Table 1.

Demographic and clinical data for 351 GC patients.

Additional file 2: Supplemental Figure S1.

Pearson correlation analysis of correlation between TCGA genome-wide expression profile data and CST4.

Additional file 3: Supplemental Data S1.

CST4 and its co-expressed related genes in GC tissues from the TCGA cohort.

Additional file 4: Supplemental Data S2.

GO term enrichments of co-expressed genes of CST4 in gastric cancer.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, Yy., Li, Bp., Wang, Jf. et al. Investigating the prognostic and predictive value of the type II cystatin genes in gastric cancer. BMC Cancer 23, 1122 (2023). https://doi.org/10.1186/s12885-023-11550-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-11550-6

Keywords