Skip to main content

Expression and gene regulatory network of S100A16 protein in cervical cancer cells based on data mining

Abstract

S100A16 protein belongs to the S100 family of calcium-binding proteins, which is widely distributed in human tissues and highly conserved. S100 calcium-binding proteins possess broad biological functions, such as cancer cell proliferation, apoptosis, tumor metastasis, and inflammation (Nat Rev Cancer 15:96–109, 2015). The S100A16 protein was initially isolated from a cell line derived from astrocytoma. The S100A16 protein, consisting of 103 amino acids, is a small acidic protein with a molecular weight of 11,801.4 Da and an isoelectric point (pI) of 6.28 (Biochem Biophys Res Commun 313:237–244, 2004). This protein exhibits high conservation among mammals and is widely expressed in various human tissues (Biochem Biophys Res Commun 322:1111–1122, 2004). Like other S100 proteins, S100A16 contains two EF-hand motifs that form a helix-loop-helix structural domain. The N-terminal domain and the C-terminal domain of S100A16 are connected by a "hinge" linker.S100A16 protein exhibits distinct characteristics that distinguish it from other S100 proteins. A notable feature is the presence of a single functional Ca2 + binding site located in the C-terminal EF-hand, consisting of 12 amino acids per protein monomer (J Biol Chem 281:38905–38917, 2006). In contrast, the N-terminal EF-hand of S100A16 comprises 15 amino acids instead of the typical 14, and it lacks the conserved glutamate residue at the final position. This unique attribute may contribute to the impaired Ca2 + binding capability in the N-terminal region (J Biol Chem 281:38905–38917, 2006). Studies have shown an integral role of S100 calcium-binding proteins in the diagnosis, treatment, and prognosis of certain diseases (Cancers 12:2037, 2020). Abnormal expression of S100A16 protein is implicated in the progression of breast and prostate cancer, but an inhibitor of oral cancer and acute lymphoblastic leukemia tumor cell proliferation (BMC Cancer 15:53, 2015; BMC Cancer 15:631, 2015). Tu et al. (Front Cell Dev Biol 9:645641, 2021) indicate that the overexpression of S100A16 mRNA in cervical cancer(CC) such as cervical squamous cell carcinoma and endocervical adenocarcinoma as compared to the control specimens. Tomiyama N. and co-workers (Oncol Lett 15:9929–9933, 2018) (Tomiyama, N) investigated the role of S100A16 in cancer stem cells using Yumoto cells (a CC cell line),The authors found upregulation of S100A16 in Yumoto cells following sphere formation as compared to monolayer culture.Despite a certain degree of understanding, the exact biological function of S100A16 in CC is still unclear. This article explores the role of S100A16 in CC through a bioinformatics analysis. Referencing the mRNA expression and SNP data of cervical cancer available through The Cancer Genome Atlas (TCGA) database, we analyzed S100A16 and its associated regulatory gene expression network in cervical cancer. We further screened genes co-expressed with S100A16 to hypothesize their function and relationship to the S100A16 cervical cancer phenotype.

Our results showed that data mining can effectively elucidate the expression and gene regulatory network of S100A16 in cervical cancer, laying the foundation for further investigations into S100A16 cervical tumorigenesis.

Peer Review reports

Introduction

Cervical cancer is one of the most common malignant tumors in gynecology and a serious threat to women’s health worldwide; it ranks fourth among cancers in women and is especially prevalent in the developing world [1]. Improvement in cervical cancer screening and early detection methods and the availability and increased administration of the HPV vaccine in recent years, has reduced the cervical cancer burden on public health, but its pathogenic mechanisms remain unclear. The main treatment modalities for cervical cancer include surgery, chemotherapy, radiation therapy, and targeted therapy. Clinical research on the etiology of cervical cancer is still ongoing.

Recent advancements in gene microarray, high-throughput sequencing techniques, and bioinformatics continue to provide novel approaches for disease research [1]. Our study screened the differential gene expressions in cervical cancer and healthy tissues based on the TCGA database. Through analysis of the copy number variation of hub genes, gene set variations, immune cell infiltration, weighted gene co-expression network analysis (WGCNA), and gene module function enrichment, we explored the related molecular mechanisms and provided a new research approach to the diagnosis and treatment of cervical cancer.

Materials

TCGA data acquisition and differential analysis

The TCGA database (https://portal.gdc.cancer.gov/.) is the largest cancer gene information resource for information pertaining to gene expression data, copy number variation, and single nucleotide polymorphisms (SNP). We downloaded the original mRNA expression data and known SNPs in cervical cancer for subsequent analysis. A total of 309 specimens were collected (normal group, n = 3; tumor group, n = 306) to analyze the expression differences of S100A16. We downloaded the series matrix files of GSE44001 from the NCBI GEO public database. The annotation platform was GPL14951. The data for 300 cervical cancer patients complete with expression profiles and survival information were retrieved. This study was approved by the Ethics Committee.

Expression of S100A16 protein in cervical cancer tissues

A total of 63 pairs of pathological slides were collected from cervical cancer surgical patients at the Second Hospital of Lanzhou University from January 2020 to July 2022. The pathology confirmed the presence of cervical cancer. A microscopic imaging system was used to capture images of the slides. Initially, the entire tissue was observed at 100 × magnification, followed by image acquisition at 400 × magnification. The results were photographed and recorded. After completing immunohistochemical staining on tissue microarrays (primary antibody at a ratio of 200:1), the cell nuclei were stained blue with hematoxylin, and positive expression shown by DAB staining appeared as brownish-yellow. Each specimen was photographed at 200 × magnification. The staining results were blindly assigned by three pathologists with independent diagnostic qualifications. The scoring criteria included: (1) staining intensity of the protein, and (2) percentage of positive cells. The product of (1) and (2) was used to determine the score, which was then averaged and grouped. The presence of a yellow or brown area in S100A16 protein immunohistochemical staining was considered positive. The scoring criteria for protein staining intensity were as follows: no staining (0 points), light staining (1 point), moderate staining (2 points), and strong staining (3 points). The scoring criteria for the percentage of positive cells were as follows: > 10% (1 point), 11%-50% (2 points), 51%-75% (3 points), and > 75% (4 points). Additionally, in order to understand the relationship between the expression level of S100A16 protein and clinical indicators, the expression of S100A16 protein was divided into a high-expression group (IHC score ≥ 12) and a low-expression group (IHC score < 12).

Co-expression analysis

The co-expression status of S100A16 in cervical cancer was analyzed with a pre-established filter and selection threshold of a 0.3 correlation coefficient and p-value of 0.05. After selecting the most significant genes with S100A16 expression, the “corrplot” and “circlize” packages were used to develop the heatmap and circular plot for S100A16 correlative analysis.

Immune cell infiltration analysis

The CIBERSORT algorithm was used to analyze the RNA-seq data of cervical cancer patients and identify relative proportions of 22 immune cells to uncover any correlation in gene expression level and immune response.

Gene set variation analysis (GSVA)

GSVA is a non-parametric unsupervised method to assess the enrichment of transcriptomic gene sets. Through comprehensively scoring the gene sets of interest, GSVA converts gene-level changes to pathway-level changes and then determines the biological function of the sample. We downloaded gene sets from the Molecular Signatures Database (MSigDB), developed a comprehensive score for each gene set using the GSVA algorithm, and evaluated the potential biological function changes across different samples.

Gene set enrichment analysis (GSEA)

GSEA is commonly used for disease classification and is closely related to biological significance. In this study, patients were divided into high and low S100A16 expression groups. The differences in the signaling pathways between the high and low expression groups were further analyzed using GSEA. The background gene set and pathway subtypes were the version 7.0 annotated gene sets downloaded from the MSigDB database. Differential expression analysis of signaling pathways between subtypes was performed, and the significantly enriched gene sets (adjusted p value < 0.05) were sorted based on consensus scoring.

Drug sensitivity analysis

Based on the largest pharmacogenomics database (https://www.cancerrxgene.org/) we used the R software package “pRRophetic” to predict the chemotherapy sensitivity for each tumor specimen. The estimated IC50 with each specific chemotherapy drug was obtained using the regression method, and 10-fold cross-validation was performed to test the regression and prediction accuracy using the GDSC training sets. All parameters had the default and “combat” values selected removing the batch effect and taking the average of duplicate gene expression.

Nomogram model construction

A nomogram was built based on regression analysis of gene expression level and clinical symptoms. Line segments were used with tick marks to draw on the same plane according to a predetermined scale, thereby conveying the interrelationships between variables in a predictive model. The multivariate regression model assigned scores to individual value levels of each influencing factor based on its degree of contribution in the model to the outcome variable (i.e., magnitude of the regression coefficient). The individual scores were added together to obtain the total score, ultimately allowing for calculation of the predicted value.

WGCNA

By constructing a weighted gene co-expression network, we investigated S100A16 and any associated genes. The co-expression networks for all genes in the cervical cancer datasets were constructed using the WGCNA-R package, and the top 10,000 genes with variance were screened using this algorithm with a soft threshold set to 5. The weighted adjacency matrix was converted to a topological overlap matrix (TOM) that estimated network connectivity, and hierarchical clustering constructed a tree of the TOM matrix. Different branches of the clustering tree represented different gene modules that were further differentiated by color. Based on their weighted correlation coefficients, the genes were classified by expression patterns and combined with similarly expressing genes into one module, grouping all genes into multiple modules.

Functional enrichment analysis of gene modules

R package “ClusterProfiler” was used to annotate the biological function of the modular genes selected from the WGCNA key (the magenta module was most correlated with S100A16). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) was used to evaluate the associated functional categories. The GO and KEGG enriched pathways with both p and Q values below 0.05 were considered significant categories.

Statistical analysis

Statistical analysis was performed using R programming language (version 4.2.1). For immunohistochemical scores, non-parametric test for paired data (Wilcoxon signed rank test) was performed using SPSS (version 25.0). P < 0.05 was considered statistically significant.

Results

The expression pattern and prognosis of S100A16 in cervical cancer

S100A16 expression was significantly upregulated in cervical cancer specimens (Fig. 1A). Our survival analysis, ranked by degree of S100A16 gene expression, indicated that in the GSE44001 dataset, the overall survival (OS) within the high S100A16 expression group was significantly shorter than groups with low S100A16 expression (Fig. 1B). Using clinical information and S100A16 gene expression level, we established univariate and multivariate Cox regression models and constructed forest plots that showed N in cervical cancer patients was associated with the risk and had a statistical difference between the groups (Fig. 1C).

Fig. 1
figure 1

A Expression of S100A16 mRNA in tumor and normal tissues. Blue represents control tissues, and red represents tumor tissues. B Increased S100A16 expression in cervical cancer is associated with a poor prognosis. C Univariate and multivariate forest plots. Green indicates protective factors, and red indicates risk factors

Immunohistochemical analysis of S100A16 in cervical cancer

Immunohistochemical staining of cervical cancer cells and the corresponding adjacent normal healthy cervical tissue specimens (ANT) was performed to investigate whether the S100A16 protein expression was altered in the tissues of cervical cancer patients. Our research revealed a statistically significant (p < 0.05) increase in the levels of S100A16 protein in 63 tumor specimens compared to the corresponding cervical paracancerous specimens (Fig. 2). According to the scoring criteria of immunohistochemistry, the 63 cervical cancer patients were divided into two groups: a high-expression group with 43 cases (68%) and a low-expression group with 20 cases (32%). The results indicated a significant correlation between the expression level of S100A16 protein and the pathological grading and depth of muscular infiltration in the cervix (P < 0.05). However, there was no significant difference in age, lymph node metastasis, tumor size, and tissue type among the patients (P > 0.05) (Table 1).

Fig. 2
figure 2

Expression level of S100A16 protein in human cervical cancer tumor specimens (tumor) and corresponding adjacent normal tissue specimens (ANT) detected by immunohistochemical staining

Table 1 The relationship between the expression of S100A16 and the clinical pathological indicators in cervical cancer patients

Co-expression gene analysis of S100A16 in cervical cancer

We further explored the co-expression network of S100A16 using the expression profiles of cervical cancer patients annotated in the TCGA database. A total of 2,934 genes that were significantly correlated (correlation coefficient filter criterion of 0.3 and p value of 0.05) with S100A16 expression were screened and resultant heatmaps of the top 5 genes with positive or negative correlation coefficients (Fig. 3A) along with co-expression correlative circular charts are shown in Fig. 3B.

Fig. 3
figure 3

A Gene heat map. The heatmap shows the top 10 genes with the most significant correlation in CESC patients with high or low expression of S100A16. B Circular plot showing the network of important S100A16-related genes in the high or low expression group of S100A16 in CESC patients

The relationship between S100A16 and immune cell infiltration in cervical cancer

The tumor microenvironment consists primarily of cancer cells and tumor-associated fibroblasts, immune cells, extracellular matrix, growth factors, and inflammatory factors. The composition and unique physiochemical properties of the tumor microenvironment significantly affects cancer diagnosis, survival outcome, and clinical sensitivity to treatment. By analyzing the relationship between S100A16 and the subsequent immune response within the TCGA datasets we probed the potential cancer progression molecular mechanisms influenced by S100A16. Specifically, S100A16 had significant positive correlations to population of resting mast cells and activated dendritic cells and had significant negative correlations with naïve B cells and resting CD4 memory T cells (Fig. 4).

Fig. 4
figure 4

Correlation between S100A16 gene expression and immune cell content

The signaling pathways related to S100A16 in cervical cancer

We examined the specific signaling pathways associated with S100A16 to explore the underlying molecular mechanisms by which S100A16 influenced tumor progression. The GSVA results showed that patients with high S100A16 expression had enrichment in several signaling pathways, such as the p53 regulatory pathway and increased levels of reactive oxygen species (Fig. 5A). Moreover, GSEA indicated S100A16 enrichment in such signaling pathways as oxidative phosphorylation, proteasomes, and ribosomes (Fig. 5B, C). These findings suggest S100A16 involvement in cervical cancer development through such pathways.

Fig. 5
figure 5

A Important biological pathways related to S100A16 were obtained by GSVA in CESC; B S100A16 expression was positively correlated with the oxidative phosphorylation, proteasome, and ribosome signaling pathways. C Correlation of the genes of the oxidative phosphorylation, proteasome, and ribosome

Mutation analysis

We downloaded the processed SNP-related data of cervical cancer, selected the top 30 genes with relatively high mutation frequency, compared the differences of the mutated genes between the two groups, and constructed the mutational landscape map using the R package ComplexHeatmap (Fig. 6). Our results showed differences in gene mutation frequency between high and low S100A16 expression levels.

Fig. 6
figure 6

Waterfall plot of S100A16-related mutation genes obtained through mutation analysis in CESC

The relationship between chemotherapy drugs and S100A16 expression in cervical cancer

The therapeutic benefit of surgery combined with chemotherapy in early-stage cervical cancer is clear. Based on drug sensitivity data in the GDSC database, our study predicted the chemotherapy sensitivity of each tumor specimen using the R package pRRophetic to further investigate any relation between S100A16 expression levels and sensitivity to common chemotherapeutic drugs. Our results showed a relation in S100A16 and sensitivity to gemcitabine, ABT.263, ABT.888, AP.24534, AS601245, and axitinib (Fig. 7).

Fig. 7
figure 7

Chemotherapy drug sensitivity analysis based on S100A16 expression

Nomogram model for cervical cancer patients

Based on the S100A16 expression levels, we presented the results of our regression analysis in the form of a nomogram. Regression analysis of patients in our study showed that the values of the different clinical indicators of cervical cancer and the distribution of S100A16 expression had varying degrees of contribution throughout the scoring process (Fig. 8A). Predictive analysis for 3-year and 5-year OS (Fig. 8B) rates indicated a trend in our nomogram model-predicted OS that was similar to the actual observed OS.

Fig. 8
figure 8

A Nomogram model for predicting 3-year and 5-year survival rates of CESC patients (B) Calibration curve for 3-year OS based on the nomogram model

Weighted gene co-expression network analysis and functional enrichment analysis

To determine the co-expression network of S100A16, we performed WGCNA. The soft threshold β was determined by the function “sft$powerEstimate,” with the soft threshold set to 5. Subsequently, the gene modules were inspected based on the TOM matrix. During analysis, a total of 16 gene modules were detected, and we found that the magenta module had the highest correlation with S100A16 (Fig. 9A, B). Subsequent pathway analysis using magenta modular genes was performed, and the GO results showed enriched genes primarily in pathways such as chromosome segregation, nuclear division, and organelle fission (Fig. 9C). Additionally, KEGG results showed enrichment mainly in pathways such as cell cycle, DNA replication, and oocyte meiosis (Fig. 9D).

Fig. 9
figure 9

A Dendrogram of sample clustering in CESC; B Gene set enrichment analysis (GSEA) of S100A16-related biological pathways in CESC; C GO analysis results of differentially expressed genes; D KEGG analysis results of differentially expressed genes

Discussion

Cervical cancer is one of the most common malignant tumors in women. At present, persistent HPV infection is the most important risk factor leading to cervical cancer. However, other factors such as age and genetic factors also influence the progression of cervical cancer [2, 3]. Current standard of care for cervical cancer is surgery and chemotherapy; however, the prognosis of cervical cancer patients remains uncertain. Therefore, it is particularly imperative to search for new diagnostic methods, markers of tumor progression, and novel therapeutic targets.

As a member of S100 proteins, S100A16 was initially isolated from astrocytomas [4, 5]. It can promote chromosomal rearrangement and instability and in this respect, it is oncogenic [6, 7]. In addition, S100A16 showed aberrant expressions in different neoplastic tissues further suggesting its involvement in tumorigenesis [8,9,10,11,12,13,14,15]. Previous studies have found a correlation between S100A14 and the occurrence and progression of cervical cancer. Overexpression of S100A14 is closely associated with the staging and lymph node metastasis of cervical cancer. S100A14 has been shown to promote cell cycle progression, cell growth, migration, and invasion of cervical cancer cells [16]. In this study, a close correlation was observed between the high expression of S100A16 and GJB3. GRAEBER et al. [17] conducted in vitro experiments using HeLa cells transfected with GJB3 and found that the expression of the junction protein supports the invasion of tumor cells into normal tissues, and this process does not require heterotypic gap junction coupling. As a member of the S100 protein family,Studies have shown that S100A16 overexpression was associated with varying tumor occurrence and progression rates [10, 11, 18]. High S100A16 expression in lung adenocarcinoma, ovarian cancer, cervical cancer and breast cancer were associated with a poor prognosis, whereas high expression in colorectal cancer and oral squamous cell carcinoma was associated with a favorable prognosis [8,9,10,11, 19, 20]. By analyzing the differences in S100A16 expression between tumor tissues and normal tissues in the TCGA datasets, we found that S100A16 expression was significantly increased in tumor tissue samples, and cervical cancer patients with high S100A16 expression had more rapid disease progression and shorter survival time. These results implicate S100A16 as potential therapeutic target in cervical cancer. However, the exact mechanisms by which it affects tumor progression remain unclear and a topic for further investigation.

The tumor microenvironment has a significant impact on cancer diagnosis, survival outcome, and sensitivity to clinical treatment. Among them, the immune system establishes the microenvironment adaptive to the cancer cells through modulating various cytokines and chemokines, thereby promoting or repressing tumor growth [21]. Tumor cells leverage the autogenous regulation of the immune system to build an immunosuppressive network in the tumor microenvironment. Cells regulating the immunosuppressive network include regulatory T cells (Tregs), dendritic cells (DCs), and natural killer cells (NKs) [22]. Because of the interaction between tumor cells and immune cells, large quantities of regulatory cells and inhibitory factors that are detrimental to immune response and favorable to tumor growth, migrate to the vicinity of tumor cells [23]. These immune effectors and suppressors are closely related to cancer development. On the other hand, the tumor microenvironment may also influence the metabolic reprogramming of tumors, maintaining immune cell function and modulating the immune system through metabolic competition and symbiosis [24]. Through the relationship between hub genes and tumor immune infiltration from the TCGA datasets, this study investigated the potential molecular mechanisms of hub genes and cervical cancer progression. We discovered that S100A16 was significantly positively correlated with resting mast cells, dendritic cells, and activated t cells, and significantly negatively correlated with naïve B cells and resting CD4 memory T cells. Inside tumors [25], MCs interact with the infiltrated immune cells, cancer cells, and extracellular matrix through direct cell-to-cell interactions or release of mediators capable of reshaping the tumor microenvironment. By releasing the classical or non-classical proangiogenic factors, MCs actively promote angiogenesis and induce neovascularization. In addition, MCs support tumor invasion by releasing extensive matrix metalloproteinases (MMPs) [18]. Furthermore, tumor cells induced the differentiation of neighboring dendritic precursors to the Gr-1( +) conventional dendritic cell subpopulation and binding with cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4), suppressing their proliferation and promoting the immune escape of cancer cells. Thus, inhibiting Gr-1( +) and CTLA-4 can improve tumor immune response [26]. In the present study, the positive correlation of dendritic cells and S100A16 confirmed previous findings. T cells are produced by activated monocytes/macrophages in lymph nodes. Through the expression of T cell receptor (TCR) α/β, CD4+ or CD8+ T cells recognize tumor antigens and autoantigens and act on specific cancer cells, thereby exerting anti-tumor immunity. For instance, miR-18a acts by inhibiting proliferation and inducing cell death of CD4+ T cells [27]. Moreover, Tregs often accumulate in tumor tissues maintaining the immunosuppressive environment of the tumor and facilitating invasion and metastasis [28, 29].

Genomic instability is one of the hallmarks of cancer, characterized by an increased rate of changes in the cell genome, which promotes cancer progression and resistance to treatment. The most common form of genomic instability in cancer is chromosomal instability, driving uncontrolled cell proliferation and tumor development [30]. CC cells also exhibit abnormal chromosomes, which contain various gene rearrangements, including translocations, deletions, and gene amplifications. Given the importance of spindle-kinetochore interactions during cell division, any defects in mitosis are associated with chromosomal instability. The instability may be attributed to chromosomal segregation defects [31].

During tumor development, the enhancement of glycolysis in cells is mainly due to irreversible damage to oxidative phosphorylation function. The mechanisms leading to oxidative phosphorylation pathway damage and the primary energy metabolic pathways (glycolysis and oxidative phosphorylation) relied upon differ among different types of tumor cells. However, almost all tumor cells exhibit varying degrees of impairment in mitochondrial oxidative phosphorylation function. The specific mechanisms still require further in-depth research.

High expression of the S100A16 protein was correlated to disease progression and prognosis of cervical cancer, and thus we investigated the hub genes and sensitivities to common chemotherapeutic drugs [32]. Our study results suggest that cervical cancer patients with high S100A16 expression were highly sensitive to drugs, such as gemcitabine and axitinib. Treating SiHa cells using an agonistic anti-CD40 monoclonal antibody or gemcitabine alone could not inhibit the proliferation of SiHa cells in vitro, whereas the activation of CD40 on SiHa cells enhanced their sensitivity to gemcitabine [32]. Preliminary research showed that when used concomitantly with cisplatin, gemcitabine may be helpful in the treatment of recurrent or advanced cervical cancer [33].

The biological functions of S100A16 in tumors remain obscure. Studies have shown that S100A16 activates the AKT signaling pathway in prostate cancer to promote cell invasion, metastasis, and proliferation [18]. Additionally, the activation of the AKT cellular pathway promotes cell survival and inhibits apoptosis [34, 35]. S100A16 plays a crucial role in cisplatin resistance during chemotherapy for the treatment of lung cancer. Zhou et al. reported that S100A16 promotes epithelial to mesenchymal transition (EMT) through the Notch pathway in breast cancer [36], while EMT enhanced the invasion and metastasis of epithelial cells and was associated with resistance to chemotherapy in a variety of tumors [37, 38]. On the contrary, in oral squamous carcinoma, the prognosis was better for patients with high S100A16 expression [10] suggesting that S100A16 had different buffering or inhibitory functions in adenocarcinoma and squamous carcinoma in an expression-dependent manner. A study by Li et al. showed that S100A16 enhanced the expression of TWIST1 by activating the STAT3 signaling pathway, which subsequently promoted the EMT and invasiveness of pancreatic cancer cells [39]. In addition, inhibition of S100A16 expression could slow the metastasis of pancreatic cancer cells. Ou et al. showed that S100A16 could repress the proliferation, migration, and invasion of colorectal cancer cells through the JNK/p38 MAPK pathway [40]. A study by Zhang et al. showed that S100A16 [41] promoted the proliferation, migration, and tumor angiogenesis of cervical cancer HeLa cells via regulating the PI3K/PKB signal transduction pathway, while the PI3K/PKB signaling pathway was closely related to cancer proliferation, invasion, differentiation, and drug resistance [42, 43]. Zhu et al. showed that S100A16 could promote the proliferation of prostate cancer cells through the AKT and ERK signaling pathways [18]. The results of this study showed that high S100A16 expression significantly enriched the p53 and apoptosis pathways. In addition, research has shown that CIZAR induced apoptosis not only by reinstating the p53/Rb-dependent pathway in HPV-positive cells, but also by activating the p53/Rb-independent pathway and mitochondrial death signaling pathway in cervical carcinoma cells, which was unrelated to HPV infection [44]. SOX14 overexpression in cervical carcinoma triggered the accumulation of p53, indicating that potential interactions exist between the SOX14 and p53 signaling pathways [45]. 

Conclusions

In summary, S100A16 mRNA and protein were abnormally upregulated in cervical cancer, and their overexpression indicated a poor prognosis and malignant tumor progression in cervical cancer patients. These results preliminarily present an opportunity for improvement in the diagnosis and treatment of cervical cancer, by revealing S100A16 as a new target for cervical cancer therapy and uncovering its apparent role in cellular susceptibility to an array of current chemotherapeutic drugs. Global results of this study also implicated S100A16 as an up-regulator of multiple tumor growth progression processes, highlighting the potential importance of continued investigation and further studies into the pathogenic mechanisms of S100A16 in cervical cancer cells and its relationship to immune response and key cellular pathways. 

Availability of data and materials

The initial data used to support the findings of this study are available from the corresponding author upon request.

Publicly available datasets were analyzed in this study. This data can be found here: TCGA-CESC (https://portal.gdc.cancer.gov/).

GSE44001 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE44001).

References

  1. 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–424r.

    Article  PubMed  Google Scholar 

  2. Ahn WS, Bae SM, Lee JM, Namkoong SE, Han SJ, Cho YL, Nam GH, Seo JS, Kim CK, Kim YW. Searching for pathogenic gene functions to cervical cancer. Gynecol Oncol. 2004;93(1):41–8.

    Article  CAS  PubMed  Google Scholar 

  3. Panjkovic M, Ivkovic-Kapicl T. Etiology and pathogenesis of precancerous lesions and invasive cervical carcinoma. Med Pregl. 2008;61(7–8):364–8.

    Article  PubMed  Google Scholar 

  4. Marenholz I, Heizmann CW, Fritz G. S100 proteins in mouse and man: from evolution to function and pathology (including an update of the nomenclature). Biochem Biophys Res Commun. 2004;322(4):1111–22.

    Article  CAS  PubMed  Google Scholar 

  5. 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.

  6. Chen H, Xu C, Jin Q, Liu Z. S100 protein family in human cancer. Am J Cancer Res. 2014;4(2):89–115.

    PubMed  PubMed Central  Google Scholar 

  7. 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 

  8. Tanaka M, Ichikawa-Tomikawa N, Shishito N, Nishiura K, Miura T, Hozumi A, Chiba H, Yoshida S, Ohtake T, Sugino T. Coexpression 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  Google Scholar 

  9. Bai Y, Li LD, Li J, Lu X. Prognostic values of S100 family members in ovarian cancer patients. BMC Cancer. 2018;18(1):1256.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. 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  Google Scholar 

  11. Saito K, Kobayashi M, Nagashio R, Ryuge S, Katono K, Nakashima H, Tsuchiya B, Jiang SX, Saegusa M, Satoh Y, et al. S100A16 is a prognostic marker for lung adenocarcinomas. Asian Pac J Cancer Prev. 2015;16(16):7039–44.

    Article  PubMed  Google Scholar 

  12. 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 

  13. 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 

  14. Allgower C, Kretz AL, von Karstedt S, Wittau M, Henne-Bruns D, Lemke J. Friend or Foe: S100 Proteins in Cancer. Cancers. 2020;12(8):2037.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Tu G, Gao W, Li Y, Dian Y, Xue B, Niu L, Yu X, Zhu H. Expressional and prognostic value of S100A16 in pancreatic cancer via integrated bioinformatics analyses. Front Cell Dev Biol. 2021;9.

  16. Wang X, Yang J, Qian J, Liu Z, Chen H, Cui Z. S100A14, a mediator of epithelial-mesenchymal transition, regulates proliferation, migration and invasion of human cervical cancer cells. Am J Cancer Res. 2015;5(4):1484–95.

    PubMed  PubMed Central  Google Scholar 

  17. Graeber SH, Hülser DF. Connexin transfection induces invasive properties in HeLa cells. Exp Cell Res. 1998;243(1):142–9.

    Article  CAS  PubMed  Google Scholar 

  18. Zhu W, Xue Y, Liang C, Zhang R, Zhang Z, Li H, Su D, Liang X, Zhang Y, Huang Q, et al. S100A16 promotes cell proliferation and metastasis via AKT and ERK cell signaling pathways in human prostate cancer. Tumour Biol. 2016;37(9):12241–50.

    Article  CAS  PubMed  Google Scholar 

  19. 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. 2018;117(2):275–83.

    Article  CAS  PubMed  Google Scholar 

  20. Tomiyama N, Ikeda R, Nishizawa Y, Masuda S, Tajitsu Y, Takeda Y. S100A16 up-regulates Oct4 and Nanog expression in cancer stem-like cells of Yumoto human cervical carcinoma cells. Oncol Lett. 2018;15(6):9929–33.

    PubMed  PubMed Central  Google Scholar 

  21. Fan JQ, Wang MF, Chen HL, Shang D, Das JK, Song J. Current advances and outlooks in immunotherapy for pancreatic ductal adenocarcinoma. Mol Cancer. 2020;19(1):32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Rabinovich GA, Gabrilovich D, Sotomayor EM. Immunosuppressive strategies that are mediated by tumor cells. Annu Rev Immunol. 2007;25:267–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Zou W. Immunosuppressive networks in the tumour environment and their therapeutic relevance. Nat Rev Cancer. 2005;5(4):263–74.

    Article  CAS  PubMed  Google Scholar 

  24. Wang T, Liu G, Wang R. The intercellular metabolic interplay between tumor and immune cells. Front Immunol. 2014;5:358.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Komi DEA, Redegeld FA. Role of Mast Cells in Shaping the Tumor Microenvironment. Clin Rev Allergy Immunol. 2020;58(3):313–25.

    Article  CAS  PubMed  Google Scholar 

  26. Diao J, Zhao J, Winter E, Cattral MS. Tumors suppress in situ proliferation of cytotoxic T cells by promoting differentiation of Gr-1(+) conventional dendritic cells through IL-6. J Immunol. 2011;186(9):5058–67.

    Article  CAS  PubMed  Google Scholar 

  27. Egeland NG, Jonsdottir K, Aure MR, Sahlberg K, Kristensen VN, Cronin-Fenton D, Skaland I, Gudlaugsson E, Baak JPA, Janssen EAM. MiR-18a and miR-18b are expressed in the stroma of oestrogen receptor alpha negative breast cancers. BMC Cancer. 2020;20(1):377.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Michalek RD, Gerriets VA, Jacobs SR, Macintyre AN, MacIver NJ, Mason EF, Sullivan SA, Nichols AG, Rathmell JC. Cutting edge: distinct glycolytic and lipid oxidative metabolic programs are essential for effector and regulatory CD4+ T cell subsets. J Immunol. 2011;186(6):3299–303.

    Article  CAS  PubMed  Google Scholar 

  29. Sharma P, Hu-Lieskovan S, Wargo JA, Ribas A. Primary, adaptive, and acquired resistance to cancer immunotherapy. Cell. 2017;168(4):707–23.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Shah P, Wolf K, Lammerding J. Bursting the bubble - nuclear envelope rupture as a path to genomic instability? Trends Cell Biol. 2017;27(8):546–55.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Hämälistö S, Stahl JL, Favaro E, Yang Q, Liu B, Christoffersen L, Loos B, Guasch Boldú C, Joyce JA, Reinheckel T, et al. Spatially and temporally defined lysosomal leakage facilitates mitotic chromosome segregation. Nat Commun. 2020;11(1):229.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Huang Q, Qu QX, Xie F, Hu JM, Chen YG, Zhang XG. Sensitization of SiHa cell to gemcitabine by CD40 activation and its overexpression in cervical carcinoma. Med Oncol. 2011;28(3):781–8.

    Article  CAS  PubMed  Google Scholar 

  33. Mutch DG, Bloss JD. Gemcitabine in cervical cancer. Gynecol Oncol. 2003;90(2 Pt 2):S8–15.

  34. Papadimitrakopoulou V. Development of PI3K/AKT/mTOR pathway inhibitors and their application in personalized therapy for non-small-cell lung cancer. J Thorac Oncol. 2012;7(8):1315–26.

    Article  CAS  PubMed  Google Scholar 

  35. Lin Y, Wang Z, Liu L, Chen L. Akt is the downstream target of GRP78 in mediating cisplatin resistance in ER stress-tolerant human lung cancer cells. Lung Cancer (Amsterdam, Netherlands). 2011;71(3):291–7.

  36. 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(1):97.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Du B, Shim JS. Targeting epithelial-mesenchymal transition (EMT) to overcome drug resistance in cancer. Molecules. 2016;21(7):965r.

  38. Sung WJ, Park KS, Kwak SG, Hyun DS, Jang JS, Park KK. Epithelial-mesenchymal transition in patients of pulmonary adenocarcinoma: correlation with cancer stem cell markers and prognosis. Int J Clin Exp Pathol. 2015;8(8):8997–9009.

    PubMed  PubMed Central  Google Scholar 

  39. Li T, Ren T, Huang C, Li Y, Yang P, Che G, Luo L, Chen Y, Peng S, Lin Y, et al. S100A16 induces epithelial-mesenchymal transition in human PDAC cells and is a new therapeutic target for pancreatic cancer treatment that synergizes with gemcitabine. Biochem Pharmacol. 2021;189:114396.

  40. Ou S, Liao Y, Shi J, Tang J, Ye Y, Wu F, Wang W, Fei J, Xie F, Bai L. S100A16 suppresses the proliferation, migration and invasion of colorectal cancer cells in part via the JNK/p38 MAPK pathway. Mol Med Rep. 2021;23(2):164.

    Article  CAS  PubMed  Google Scholar 

  41. Zhang H, Yang Y, Ma X, Xin W, Fan X. S100A16 regulates HeLa cell through the phosphatidylinositol 3 Kinase (PI3K)/AKT signaling pathway. Med Sci Monit. 2020;26:e919757.

  42. Chen Y, Sun Z, Qi M, Wang X, Zhang W, Chen C, Liu J, Zhao W. INPP4B restrains cell proliferation and metastasis via regulation of the PI3K/AKT/SGK pathway. J Cell Mol Med. 2018;22(5):2935–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Shu XR, Wu J, Sun H, Chi LQ, Wang JH. PAK4 confers the malignance of cervical cancers and contributes to the cisplatinresistance in cervical cancer cells via PI3K/AKT pathway. Diagn Pathol. 2015;10:177.

  44. Bae SN, Lee KH, Kim JH, Lee SJ, Park LO. Zinc induces apoptosis on cervical carcinoma cells by p53-dependent and -independent pathway. Biochem Biophys Res Commun. 2017;484(1):218–23.

    Article  CAS  PubMed  Google Scholar 

  45. Stanisavljevic D, Petrovic I, Vukovic V, Schwirtlich M, Gredic M, Stevanovic M, Popovic J. SOX14 activates the p53 signaling pathway and induces apoptosis in a cervical carcinoma cell line. PLoS ONE. 2017;12(9):e0184686.

Download references

Acknowledgements

Not applicable.

Funding

1.Natural Science Foundation of Gansu Province(22JR11RA073).

2.Lanzhou University Education Development Foundation donation(2021620005003643).

Author information

Authors and Affiliations

Authors

Contributions

Haibin zhang was the lead PI and refined the manuscript. Yongxiu Yang,Wenhu Xing and Yufeng drafted the manuscript. Shan zhang and Haibin zhang produced all the figures. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yongxiu Yang.

Ethics declarations

Ethics approval and consent to participate

The study was conducted in accordance with the Declaration of Helsinki and was approved by the Institutional Review Board of the Ethics Review Committee of the Second Hospital of the lanzhou University [approval number:2022A-262]. Informed consent was obtained from all study participants or their legal guardian.

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.

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

Zhang, H., Yang, Y., Xing, W. et al. Expression and gene regulatory network of S100A16 protein in cervical cancer cells based on data mining. BMC Cancer 23, 1124 (2023). https://doi.org/10.1186/s12885-023-11574-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-023-11574-y

Keywords