Abstract
Introduction
S100 calcium-binding protein A2 (S100A2) is associated with various tumors. However, its expression profile, clinical relevance, and prognostic value in hepatocellular carcinoma (HCC) remain unclear; therefore, this study assessed S100A2 expression levels in HCC and adjacent normal tissues.
Methods
To investigate the role of S100A2 in HCC, RNA sequencing and DNA methylation data were obtained from The Cancer Genome Atlas (TCGA)-Liver Hepatocellular Carcinoma (LIHC) cohort comprising 374 tumor and 50 normal liver tissues. A retrospective cohort of 216 HCC patients was also evaluated for correlations between S100A2 expression and clinicopathological characteristics. In a subset of 62 paired tumor and adjacent normal tissues, S100A2 protein and mRNA levels were assessed by immunohistochemistry (IHC) and quantitative RT-PCR. Finally, the relationship between S100A2 overexpression and clinicopathological variables was examined using the Cox proportional hazards regression model.
Results
Analysis of the TCGA-LIHC dataset revealed a marked elevation in S100A2 expression in tumor tissues compared to normal liver tissues. Consistently, DNA methylation analysis showed hypomethylation of several S100A2-associated CpG sites in liver hepatocellular carcinoma, suggesting a potential epigenetic mechanism for its upregulation. Correlation analysis demonstrated that increased S100A2 expression was associated with advanced histological grade, lymph node metastasis, serum alpha-fetoprotein level, microvascular invasion, tyrosine kinase inhibitor level, concurrent treatment, and higher Tumor, Node, Metastasis stage. Univariate analysis showed that elevated S100A2 levels were associated with significantly poorer recurrence-free survival (RFS) and overall survival (OS). Moreover, multivariate analysis identified S100A2 as an independent prognostic indicator for both RFS and OS. Kaplan–Meier survival curves also confirmed that patients with high S100A2 protein levels had significantly worse 5-year OS and RFS rates.
Conclusion
These findings indicate that S100A2 overexpression is associated with poor prognosis in patients with HCC, highlighting its potential utility as a diagnostic biomarker.
Plain Language Summary
Keywords
1. Introduction
Hepatocellular carcinoma (HCC) is a common malignant solid tumor worldwide, accounting for 80–90% of all primary liver cancer cases. 1 Treatment modalities for HCC include liver transplantation, hepatic resection, radiotherapy, percutaneous ablation, and intra-arterial and systemic therapies. 2 However, owing to its insidious onset, less than 30% of patients are eligible for curative treatment at initial diagnosis. 3 Although new systemic therapies have emerged, survival rates for advanced HCC remain low. 4 Therefore, identifying novel therapeutic targets and reliable prognostic biomarkers for HCC is essential.
S100A2, a calcium-binding protein member of the S100 family encoded on chromosome 1q21,5,6 significantly influences multiple oncogenic processes in HCC. 7 It functions as a key regulator of intracellular calcium homeostasis and is involved in cell proliferation, 8 migration, 6 invasion, 9 and apoptosis. 10 S100 proteins, including S100A2, have been implicated in key aspects of HCC progression, including tumor cell proliferation, 11 metastasis, 12 and apoptosis. 13 For instance, S100A2 silencing abrogates the suppressive effect of NFIC1 on cancer cell invasion and migration, thereby promoting metastasis by activating the MEK/ERK signaling pathway. 14 The role of S100A2 in tumorigenesis is context-dependent, exhibiting both tumor-promoting and -suppressive functions in different cancer types. For instance, S100A2 is upregulated in endometrial 15 and nasopharyngeal carcinomas, 16 where it facilitates tumor progression and correlates with poor prognosis, whereas in gastric cancer, its loss is associated with lymph node metastasis and unfavorable outcomes. 17 In HCC, S100A2 upregulation in HCC is driven by the hypoxia/hypoxia-inducible factor-1 alpha pathway and epigenetic modulation, 18 and high S100A2 expression is linked to shorter overall survival (OS) in patients with HCC. 19 Nevertheless, compared to well-characterized S100 family members, such as S100A4 20 and S100A6, 21 the expression profile, clinical relevance, and prognostic significance of S100A2 in HCC remain unclear.
This study aimed to assess S100A2 expression levels in HCC and adjacent normal tissues using data from The Cancer Genome Atlas (TCGA), immunohistochemistry (IHC), and quantitative reverse transcription polymerase chain reaction (RT-qPCR), and to further analyze its associations with clinicopathological features and patient prognosis. The primary objective of this study was to determine its prognostic significance and potential as a therapeutic target for HCC.
2. Materials and Methods
2.1. Patients
This retrospective study included 216 patients with HCC treated at Longyan Second People’s Hospital between January 2012 and December 2019. All the participants signed a written informed consent form allowing the use of their clinical data and tissue samples for research purposes. The medical data of all patients were carefully evaluated and approved for use by the Longyan Second People’s Hospital Medical Ethics Committee (IRB approval number: LYEYEC 2023-026, IRB approval date: May 21, 2023). The inclusion criteria were as follows: (i) histopathologically confirmed HCC; (ii) primary liver resection; (iii) no concurrent malignancies; and (iv) complete clinical data. Clinicopathological parameters, including sex, age, cirrhosis, HBsAg status, tumor size, pathological stage, histological grade, serum alpha-fetoprotein (AFP) level, microvascular invasion, tyrosine kinase inhibitor (TKI) level, concurrent treatment, and lymph node metastasis (LNM) status, were extracted from the medical records. All cases were pathologically confirmed and staged according to the Union for International Cancer Control (UICC) TNM, 8th edition. All patient information was de-identified, and individual identities could not be recognized. This study adhered to REMARK guidelines. 22
Postoperative follow-up was conducted via telephone interviews and clinical visits. The median follow-up duration was 66 months (range: 5–95 months), spanning from the date of surgery to the last follow-up (October 26, 2024) or death. OS is defined as the time elapsed from the commencement of treatment or surgery until death, the patient’s final follow-up (October 26, 2024), or death for any reason. According to this definition, the recurrence-free survival (RFS) period is the interval from the start of treatment or surgery to the recurrence or reappearance of the tumor.
2.2. TCGA Data Analysis
RNA sequencing datasets and DNA methylation profiles for liver hepatocellular carcinoma (LIHC), including 374 tumor and 50 normal liver tissues, were obtained from TCGA. Normalized S100A2 expression levels were analyzed for differential expression between tumor and normal groups. For DNA methylation, CpG sites associated with S100A2 were identified, and their beta values were visualized using a heatmap.
2.3. Immunohistochemistry
Following dewaxing and rehydration, IHC for S100A2 was performed on 62 paired 4-μm paraffin-embedded HCC and adjacent normal tissue sections according to the manufacturer’s instructions. Subsequently, sections were subjected to antigen retrieval and peroxidase blocking. They were then incubated for 1 h with rabbit monoclonal anti-S100A2 antibody (dilution 1:200, ab109494; Abcam, Cambridge, UK). Positive staining, characterized by brown chromogenic deposition in the nucleus and cytoplasm of the target cells, was confirmed by counterstaining with Mayer’s hematoxylin. Finally, all slides were examined under an optical microscope (Olympus BX53; Olympus Corporation, Tokyo, Japan) at 200× magnification, and staining specificity was verified using positive and negative controls.
2.4. Evaluation of the IHC Results
IHC results were evaluated by recording the percentage of S100A2-stained cells, followed by an assessment of intracellular staining intensity as described previously. 23 Two independent observers evaluated the fraction of immunopositive cells within the tumor matrix using the following criteria: 0 (0–10%), 1 (>10–30%), 2 (>30–50%), 3 (>50–70%), and 4 (>70–100%). Subsequently, staining intensity was scored semi-quantitatively as 0 (negative), 1 (mild), 2 (moderate), or 3 (intense). The overall score (range, 0–7) was calculated by summing the scores for positive cell percentage and staining intensity. The final score for each section was established by consensus among the observers. To assess the interrater reliability of the IHC scoring system, 30 cases were randomly selected for double-blind evaluation by two independent pathologists. Interobserver agreement was assessed using Cohen’s kappa (κ) coefficient.
2.5. RT-qPCR Assay
After RNA extraction with TRIzol reagent (Invitrogen, Waltham, MA, USA) and cDNA synthesis using RevertAid™ reverse transcriptase (Takara, Shiga, Japan), qPCR was conducted on an ABI PRISM 7500 system (Applied Biosystems, Waltham, MA, USA). The thermal cycling profile consisted of a 2-min hold at 50°C, a 10-min denaturation at 95°C, and 50 cycles of 95°C for 15 s and 60°C for 1 min. S100A2 was amplified using the primers 5′-GAAGGAACTTCTGCACAAGG-3′ (forward) and 5′-GTGCCAGGAAAACAGCATAC-3′ (reverse). Experiments were performed in triplicate, and relative mRNA expression was quantified using the 2−ΔΔCt method, with data presented as the mean of three independent replicates. A standard curve was generated using 10-fold serial dilutions of the cDNA. Ct values were plotted against log10 of template concentration, and the slope was used to calculate primer efficiency using formula E = (10−1/slope−1) × 100%E = (10^{-1/slope}-1)\times 100\%E = (10−1/slope−1) × 100%. The efficiency was XX% and R2 = 0.99.
2.6. Statistical Analyses
SPSS (version 26.0) was used for all data analyses, and P < 0.05 was considered statistically significant. The relationship between S100A2 IHC staining in adjacent normal (paired adjacent non-cancerous liver tissues) and tumor tissue was determined by conducting either a chi-square test or a t-test. The association between S100A2 expression and clinicopathological variables was explored using Fisher’s exact test. OS was defined as tumor-related mortality, whereas the onset of any local or distant recurrence defined RFS, whichever occurred initially. The discrepancies in RFS and OS between the low and high S100A2 groups were assessed using univariate and multivariate analyses.
3. Results
3.1. S100A2 Protein Expression is Upregulated in Patients with HCC
Analysis of TCGA-LIHC dataset revealed that S100A2 mRNA expression was significantly higher in tumor tissues than in normal liver tissues (374 vs. 50, P < 0.001; Figure 1A). DNA methylation analysis revealed distinct patterns at the CpG sites of S100A2 in tumor and normal tissues, with multiple sites exhibiting differential methylation levels (Figure 1B). Representative IHC images of S100A2 are shown in Figure 1C. The average IHC staining intensity of S100A2 was significantly elevated in HCC tissues compared with normal tissues (4.113 ± 0.24 vs. 2.806 ± 0.22; Figure 1D). Furthermore, S100A2 protein expression was significantly higher in patients with HCC than in normal tissues (34/62 vs. 22/62, P < 0.001). Specifically, 34 cases showed higher S100A2 expression in tumor tissue, 22 cases showed lower expression, and six pairs showed no significant difference. The inter-rater agreement analysis demonstrated substantial consistency, with κ values of 0.800 for normal tissues and 0.726 for tumor tissues (Supplementary Material). S100A2 mRNA expression was markedly increased in tumor tissues compared to that in normal tissues (P < 0.001; Figure 1E). Transcriptional expression and methylation analysis of S100 calcium-binding protein A2 (S100A2) in liver hepatocellular carcinoma (LIHC). (A) S100A2 mRNA expression in tumor and normal liver tissues from The Cancer Genome Atlas (TCGA)-Liver Hepatocellular Carcinoma (LIHC) cohort (n = 374 tumors, n = 50 normal samples). S100A2 was significantly upregulated in tumor tissues (**P < 0.01). (B) DNA methylation patterns of CpG sites within S100A2 in tumor and normal tissues. Multiple sites showed differential methylation levels. (C) Representative immunohistochemical staining of S100A2 in 62 paired tumor and adjacent normal tissues. Positive staining was primarily cytoplasmic. (D) Scatter plot comparing S100A2 protein staining scores (mean ± standard error of the mean) between normal and tumor tissues. Statistical significance was evaluated using a paired t-test (***P < 0.001). (E) Quantitative RT-PCR analysis of S100A2 mRNA expression level in paired liver hepatocellular carcinoma (LIHC) and normal tissues. Tumor tissues exhibited significantly higher S100A2 levels (***P < 0.001)
3.2. S100A2 Expression and Clinical Criteria Correlation in HCC
Relationship Between High and Low S100A2 Expression and Different Clinical Indicators in 216 Patients With Hepatocellular Carcinoma (HCC)
S100A2, S100 calcium-binding protein A2; HBsAg, hepatitis B surface antigen; AFP, alpha-fetoprotein; TKI, tyrosine kinase Inhibitor; LNM, lymph node metastasis. Bold values indicate P < 0.05.
3.3. Survival Analysis
Univariate and Multivariable Analysis of Factors Potentially Predictive of Overall Survival in 216 Patients With HCC
HCC, hepatocellular carcinoma; S100A2, S100 calcium-binding protein A2; HBsAg: hepatitis B surface antigen; AFP, alpha-fetoprotein; TKI, tyrosine kinase inhibitor; LNM, lymph node metastasis. Bold values indicate P < 0.05.
Univariate and Multivariate Analysis of Factors Potentially Predictive of Recurrence-free Survival in 216 Patients With HCC
HCC, hepatocellular carcinoma; S100A2, S100 calcium-binding protein A2; HBsAg: Hepatitis B surface antigen; AFP, alpha-fetoprotein; TKI, tyrosine kinase inhibitor; LNM, lymph node metastasis. Bold values indicate P < 0.05.
Kaplan–Meier analysis indicated that increased S100A2 expression was associated with unfavorable outcomes in patients with HCC (all P < 0.05) (Figure 2). In TCGA-LIHC cohort (n = 374), Kaplan–Meier survival analysis showed that patients with high S100A2 expression had significantly poorer OS than those with low S100A2 expression (P = 0.009, Figure 2A). Spearman’s correlation analysis showed that S100A2 expression was negatively correlated with OS (r = –0.16; Figure 2B). Moreover, the estimated 5-year OS and RFS rates in patients with HCC with different S100A2 expression levels (OS: 48.4% for low vs. 33.6% for high; RFS: 45.4% for low vs. 26.8% for high) corroborate this observation (Figure 2C–D). Correlation and survival analysis of S100A2 expression in hepatocellular carcinoma (HCC). (A) Kaplan–Meier overall survival curve of S100A2 expression in The Cancer Genome Atlas (TCGA)-Liver Hepatocellular Carcinoma (LIHC) cohort. (B) Spearman’s correlation heatmap between S100A2 expression and overall survival (OS) in TCGA-LIHC. Survival analysis based on S100A2 expression in patients with HCC. (C) Kaplan–Meier curves showing 5-year OS in patients with HCC, as well as OS according to S100A2 levels in 216 patients with HCC. (D) Kaplan–Meier curves showing the 5-year recurrence-free survival (RFS) rate in patients with HCC, as well as RFS according to S100A2 levels in 216 patients with HCC
4. Discussion
In this study, S100A2 expression was notably increased in HCC tissues compared with adjacent normal tissues. Increased S100A2 expression was significantly associated with TNM stage, histological grade, serum AFP levels, microvascular invasion, TKI levels, concurrent treatment, and LNM status, but not with age, sex, tumor size, HBsAg, cirrhosis, or recurrence. Furthermore, high S100A2 expression was associated with shorter OS and RFS. These findings suggest that S100A2 is a potential prognostic biomarker for HCC, highlighting its value in facilitating adjuvant therapy, particularly in high-risk subpopulations. Although previous in vitro studies (Jia et al 18 ) have suggested potential regulatory mechanisms, to the best of our knowledge, this study is the first to demonstrate a direct correlation between hypomethylation of the S100A2 promoter and its protein overexpression in a clinical HCC cohort, offering novel mechanistic insight into its epigenetic regulation in human tumor tissues.
The prognostic significance of S100A2 in HCC is likely attributable to its involvement in multiple oncogenic signaling pathways. S100A2 acts as a potent activator of core oncogenic cascades, such as PI3K/AKT and TGF-β/Smad, which collectively drive key tumorigenic phenotypes, including sustained proliferation, enhanced migration, tissue invasion, and metastatic dissemination.24,25 Furthermore, S100A2 has also been implicated in broader oncogenic mechanisms, such as metabolic reprogramming, 26 epithelial–mesenchymal transition, 27 and modulation of anti-tumor immune responses, 28 all of which are highly relevant to HCC pathogenesis and progression. In HCC specifically, PI3K/AKT and TGF-β/Smad pathway hyperactivation are major drivers of disease progression and directly promote proliferation and metastatic spread.29,30 Given the characteristic glycolytic phenotype and high metabolic demand of HCC, S100A2 overexpression may further enhance tumor aggressiveness by reinforcing PI3K/AKT-mediated survival signaling and facilitating metabolic reprogramming, which supports rapid proliferation under nutrient stress.31,32 Moreover, HCC develops within a unique microenvironment characterized by chronic inflammation, immune tolerance, and frequent exposure to viral hepatitis- or cirrhosis-related injuries. 33 In this liver-specific context, S100A2 may contribute to an immunosuppressive tumor microenvironment by altering immune cell infiltration, thereby promoting immune evasion and potentially diminishing the efficacy of both conventional therapies and immunotherapies. 34 Collectively, these mechanistic insights provide a biologically plausible foundation for the observed association between elevated S100A2 expression and poor prognosis in HCC.
This study confirmed that S100A2 is an independent prognostic marker in HCC, consistent with previous reports (Wan et al. 19 ) linking its expression to poor survival outcomes. Unlike earlier studies that primarily emphasized correlations with survival, this study demonstrated that high S100A2 expression was systematically associated with multiple adverse clinicopathological features, providing direct tissue-level evidence of its role in tumor aggressiveness. These clinical observations are consistent with previous experimental findings indicating an oncogenic function of S100A2; its overexpression promotes tumor progression, and its silencing disrupts TGF-β–mediated invasion and migration, ultimately contributing to poor survival outcomes. 35 By integrating quantitative IHC scoring with survival analysis and methylation data, this study advances S100A2 from molecular observations to a clinically actionable biomarker framework by bridging mechanistic understanding with potential diagnostic and prognostic applications.
Numerous studies have demonstrated the roles of S100A2 in cancer therapy and drug resistance. For instance, targeting S100A2 during drug design may provide a novel therapeutic strategy for pancreatic cancer by disrupting its interaction with p53, thereby suppressing cancer cell proliferation and survival. 36 In head and neck squamous cell carcinoma, elevated S100A2 levels are a key factor in cisplatin-specific chemoresistance, where its overexpression sustains resistance even after cisplatin withdrawal. 37 Similarly, in colorectal cancer, the KRT6A/S100A2 axis promotes tumor progression, and inhibiting this pathway significantly enhances the anti-tumor efficacy of sinigrin, suggesting its potential as a therapeutic target. 38 Although these findings suggest a broader role of S100A2 in therapy resistance, direct evidence supporting the causal role of S100A2 in modulating chemotherapy sensitivity in HCC remains limited. Therefore, although S100A2 dysregulation may represent a biologically relevant feature of HCC, further mechanistic investigations and in vivo validation, including patient-derived xenograft models, are required to determine its therapeutic relevance.
This study has several limitations: First, its retrospective, single-center design may introduce potential confounding factors and biases owing to the lack of standardized follow-up protocols. Second, the limited sample size may compromise statistical power and limit the generalizability of the findings. Third, this study was primarily correlative, relying on clinical association and bioinformatics analyses without experimental validation of the mechanistic role of S100A2 in HCC.
5. Conclusion
This study identified S100A2 as an important biomarker for diagnosing and predicting HCC prognosis, and it may function as a gene that promotes tumor progression. Elucidating the molecular mechanisms underlying HCC progression could provide critical insights that guide the development of personalized targeted therapies for this disease. However, future multicenter studies with larger cohorts and more comprehensive experimental approaches are required to confirm the prognostic significance of S100A2 and further elucidate its molecular mechanisms in HCC.
Supplemental Material
Supplemental Material - Overexpression of S100 Calcium-Binding Protein A2 is Associated With Poor Prognosis in Hepatocellular Carcinoma
Supplemental Material for Overexpression of S100 Calcium-Binding Protein A2 is Associated With Poor Prognosis in Hepatocellular Carcinoma by Xiaopeng Chen, Ma Shaoqing, Zeng Wenlong, Huang Chuiguo, Jianyang Guo in Cancer Control.
Footnotes
Ethical Considerations
The medical data of all patients were carefully evaluated and sanctioned for use approved by the Longyan Second People’s Hospital’s medical ethics committee (IRB approval number: LYEYEC 2023-026, IRB approval date: May 21, 2023). This study was conducted in accordance with the Declaration of Helsinki.
Consent to Participate
Informed consent was obtained from all subjects involved in the study.
Authors’ Contributions
Conceptualization, Xiaopeng Chen and Shaoqing Ma; Methodology, Wenlong Zeng; Software, Jianyang Guo; Validation, Xiaopeng Chen, Chuiguo Huang and Shaoqing Ma; Formal analysis, Wenlong Zeng; Investigation, Shaoqing Ma; Resources, Jianyang Guo; Data curation, Chuiguo Huang; Writing—original draft preparation, Xiaopeng Chen; Writing—review and editing, Shaoqing Ma; Project administration, Xiaopeng Chen; All authors have read and agreed to the published version of the manuscript.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was supported by Natural Science Foundation of Fujian Province (2022LYF17045).
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Data Availability Statement
The data used to support the findings of this study are included within the article and the supplementary information. The data and materials in the current study are available from the corresponding author on reasonable request.
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
