KLF4 loss in hepatocellular carcinoma: Improving prognostic prediction and correlating immune infiltrates

Introduction: Although the molecular mechanisms of Krüpple-like factor 4 (KLF4) as a tumor suppressor in HCC tumorigenesis have been thoroughly examined, its clinical application in terms of precise prognostication and its influence on tumor immune microenvironment in patients with HCC require further investigation. Methods: Bioinformatics and immunohistochemistry (IHC) were used to validate KLF4 expressions in a tissue microarray (TMA) containing HCC samples. Using Cox regression models, independent prognostic factors were identified and employed in the development of nomograms. Decision curve analysis (DCA) demonstrated the superiority of the nomograms. GO and KEGG pathway analyses were applied to the functional study of KLF4. The GSVA program explored the link between KLF4 expression and tumor-infiltrating immune cells, and CAMOIP was used to construct KLF4 expression immune scores. Changes in immune-related gene markers were also investigated in relation to KLF4 expression. The association between immune cell infiltration and KLF4 expression was validated by IHC in TMA. Results: HCC was reported to have a notable depletion of KLF4. The absence of KLF4 was associated with advanced clinicopathological characteristics of HCC and predicted a bad prognosis for patients. Nomograms constructed using KLF4 expression, tumor differentiation, and TNM stage provided a more accurate prognostic assessment of HCC patients than TNM stage alone. KLF4 expression was associated with immunological-related functions, infiltration of macrophages, CD8+ T cells, and other immune cells, and elevation of immune checkpoints. Higher levels of CD8+ T cells and macrophage infiltration are associated with increased KLF4 expression in HCC TMA. Conclusion: KLF4 loss in HCC is a prognostic biomarker that influences the tumor immune microenvironment (TIME).


Introduction
Primary liver cancer is estimated to account for over one million deaths by 2030, with hepatocellular carcinoma (HCC) accounting for 85% of these cases (Sung et al., 2021). Surgical treatment continues to be the primary curative treatment for HCC patients. Nevertheless, even when HCC is treated with curative intent, tumor recurrence and metastasis are common, and almost 70% of HCC patients will relapse within 5 years after surgery, thus negatively impacting patient survival (Chen M. Y. et al., 2020;Chen et al., 2021). For cost-effective follow-up and personalized adjuvant therapy, prompt and accurate prediction of HCC recurrence and patient prognosis after resection are crucial. However, traditional pathological factors (e.g., tumor number, tumor size, and vascular invasion) cannot fully capture the biology of HCC and are insufficient for accurate prediction of patient prognosis (Mazzoccoli et al., 2016;Chen M. et al., 2020). Therefore, the development of a method that can more accurately predict the prognosis of HCC patients is urgently required in clinical practice.
Critical genetic and/or molecular events underlying hepatocarcinogenesis have been intensively studied over the past decade. KLF4 is a transcription factor that plays significant tumorsuppressive roles in the pathogenesis and progression of HCC, such as promoting tumor differentiation, suppressing epithelialmesenchymal transition, and restraining tumor growth and metastasis (Sun et al., 2017). While KLF4 has been extensively studied in cancer, the potential clinical application of KLF4 for cancer diagnosis and prognosis has not been adequately studied, and it is still a still far from clinical transition. Notably, the role of KLF4 is significantly greater than was previously believed and is highly dependent on the microenvironment in which KLF4 directs its cadre of transcriptional targets. Previous studies have found KLF4 may be a potential target for manipulating immune cells. KLF4 deficiency inhibits the functional activation of neutrophils, which are essential components of innate immunity and key regulators of adaptive immunity (Shen et al., 2017). In contrast, KLF4 can mediate the hedgehog signaling pathway-dependent polarization of M2 macrophages and disrupt the cellular NF-κB signaling pathway in macrophages to reduce Th1-type chemokine production, thereby promoting intratumoral immune suppression. Due to variations in the immune microenvironment, it appears that KLF4 has distinct transcriptional mechanisms in each cell type (Petty et al., 2019). KLF4 modulates the recruitment of myeloidderived suppressor cells via the CXCL5/CXCR2 axis in mammary tumor cells (Yu et al., 2013). Collectively, KLF4 plays significant multifaceted roles in tumor cells, immune cells infiltrating tumors, and their crosstalk. Therefore, it is reasonable to speculate that KLF4 may serve as a prognostic biomarker associated with immune infiltrations in HCC.
In general, the HCC tumor microenvironment exhibits substantial immunosuppressive activity, which correlates closely with HCC progression and immune-based treatment resistance. Biomarkers that correlate with tumor immune infiltrates should aid in the accurate prediction of prognosis and the development of effective immunotherapies for HCC. The objective of this study is to assess the role of KLF4 as a potential biomarker for the prognostic prediction of HCC after resection and its association with tumor immune infiltrates.

TMA construction and IHC
In this retrospective study, we included 141 samples from patients who received curative treatment at Shanghai General Hospital. All patients met the following inclusion criteria: 1) available baseline data and histologic material, 2) HCC diagnosed by a physician based on unambiguous histologic characteristics (Sun et al., 2011). Age, gender, hepatitis B virus surface antigens (HBsAg), liver cirrhosis, α-fetoprotein (AFP), Child-Pugh grade, tumor size, tumor number, tumor differentiation grade, tumor distribution, and pathological Tumor Node Metastasis stage (TNM stage) were all systematically recorded. All patients provided informed consent, and the academic study was conducted in accordance with ethical standards.
The primary antibody was anti-KLF4 (1:500, Santa Cruz Biotechnology). The IHC procedures were performed carried out precisely as described previously (Sun et al., 2016;Sun et al., 2017). KLF4 staining was classified into three categories based on the intensity of the staining and the percentage of positive cells: weak/negative, moderate, and strong (Wei et al., 2005). In addition, in order to differentiate between low and high KLF4 expression, negative/weak expression was deemed low, whereas moderate/strong expression was deemed high.

Tumor cohort
The TCGA database (https://portal.gdc.cancer.gov/) was used to obtain genomic data from 10,363 tumors and 730 adjacent normal tissues, as well as corresponding clinical information from 33 tumor types (Cancer Genome Atlas Research Network et al., 2013). We utilized level 3 RNA-seq data from LIHC for HCC, which included normal (n = 50) and tumor (n = 374) samples. All fragments per kilobase per million (FPKM) data were converted to transcripts per million (TPM) format and log2 (TPM+1) conversion for further research. KLF4 expression levels were assessed for each cancer type using log2 (TPM+1) of that gene, and mRNA expression data were characterized by median ± upper and lower quartiles. The ggplot2 package was used to create the visualization.

Clinical cohort analysis
The prognostic ability of KLF4 was also determined by survival analyses between the low-and high-expression groups, which were separated by median expression thresholding. Subsequently, we performed various subgroup analyses based on clinical characteristics to investigate the overall survival (OS) and disease-free survival (RFS) of KLF4.
The clinical characteristics and demographics of the HCC cohort were presented using descriptive statistics, with comparative statistical methods described in the table legend. Survival analyses were performed using the survminer (version 0.4.9) and survival package (version 3.2-10). To test the independent prognostic value of KLF4, we used Cox regression via the survival package's Cox regression function (coxph).   Significant variables associated with survival on univariate analysis were included in the multivariate models for OS and RFS. Nomograms for prognosis prediction were also created using the RMS package (version 6.2-0) and survival package (version 3.2-10), and all independent predictors from Cox models were integrated. Finally, calibration curves and DCA were plotted to evaluate the calibration ability and the clinical significance.

Identification of KLF4 roles in HCC associated with differentially expressed genes
The differentially expressed genes (DEGs) between two groups of patients (low-/high-KLF4 expression groups) were identified using the DESeq2 package (version 1.26.0) based on the median counts of KLF4 expression. Genes with an absolute value of log2 of fold change >1 and adjusted p-value < 0.05 were considered significantly differentially expressed, and their potential functions were investigated by gene ontology (GO) and KEGG enrichment using the clusterProfiler package (version 3.14.3) (Yu et al., 2012).

Immune cell infiltration analysis based on KLF4 roles in HCC
Immune cell infiltration levels in the HCC cohort were estimated using ssGSEA with the GSVA package (Hänzelmann et al., 2013). In our study, the deconvolution algorithm allowed us to distinguish 24 distinct subtypes of immune cells. The relationship between the level of expressed KLF4 and the infiltration levels of each distinct immune cell population was then characterized using Spearman's correlation analysis. In addition, the levels of infiltration of various immune cells associated with antitumor immunity were compared. Spearman's correlation was also used to examine the association between the KLF4 gene, the classical immune cell gene marker, and immune checkpoints. Concurrently, the interaction between KLF4 and immune scores was investigated in the CAMOIP database (https://www.camoip.net/) (Lin et al., 2022).

Validation of the relationship between KLF4 and immune-cell infiltration in HCC
The validation cohort TMA consisted of 268 cases of HCC with matched non-tumor adjacent tissue obtained from Shanghai General Hospital. Primary antibodies used included anti-CD8 (Invitrogen, 14-0081-82), anti-CD68 (Abcam, ab955), and anti-KLF4 (Abcam, ab215036). CD8 and CD68 were used to identify CD8 + T cells and macrophages, respectively. TMA was processed for IHC analysis of immune cell infiltration and KLF4 expression. For the score of stained markers, we used a histologic score (H-score) based on the intensity and percentage score (Paschalis et al., 2019). H-score for each scanned section was determined automatically by Servicebio (Ding et al., 2019). In addition, patients were grouped

Statistical analysis
All statistical analyses and visualizations were carried out using R version 3.6.3 (R Foundation) and the R package ggplot2. The Mann-Whitney U-test or Wilcoxon signed-rank test was applied to compare KLF4 expression between tumor and nontumor tissues. The Wilcoxon signed-rank test, the Chi-square test, and the Fisher exact test were used to assess the association between KLF4 expression and clinicopathologic features. Logrank test was utilized to compare survival curves. Cox regression was employed in both univariate and multivariate analyses.
Three methods were used to estimate the property of a nomogram: DCA, calibration curve, and Concordance index (C-index). Infiltration levels of macrophages and CD8 + T cells were estimated using the Wilcoxon signed-rank test. All correlation analyses incorporated either Pearson's or Spearman's correlation tests. The difference was considered statistically significant when the p-value was less than 0.05 (*p < 0.05, **p < 0.01, ***p < 0.001).

Low expression of KLF4 in HCC
Compared to normal tissues, the pan-cancer data analysis found that KLF4 expression was low in tumor tissues. 11 tumors exhibited low KLF4 expression: bladder urothelial carcinoma (BLCA), breast  Figure 1A). Notably, we found that KLF4 was significantly underexpressed in HCC ( Figure 1C). For additional validation, a larger sample size study of unpaired samples was done. Several tumors, including BLCA, BRCA, COAD, HNSC, KICH, LIHC, LUAD, LUSC, pheochromocytoma and paraganglioma (PCPG), rectum adenocarcinoma (READ), STAD, THCA, and UCEC, exhibited a lower expression level than normal tissue, although cholangiocarcinoma (CHOL) had a greater expression level ( Figure 1B). Notably, KLF4 expression appears to be considerably decreased in HCC ( Figure 1D). 141 HCC patients were classified according to four levels of KLF4 expression, as illustrated by representative IHC images in Figure 1E. KLF4 expression was used to categorize samples into negative/weak, moderate, and strong groups. Tumor tissues exhibited 46.1%, 47.52%, and 6.38%, whereas normal tissues exhibited 31.21% in the moderate group and 68.79% in the strong group, indicating that KLF4 expression was lower in tumors than in normal tissues ( Figure 1F). These findings suggested that KLF4 could be a useful biomarker for detecting HCC.

KLF4 impact on HCC prognosis
Based on their IHC scores, 141 patients with HCC were divided into two groups. We categorized weak/negative as the group with low expression, and the remainder (moderate/strong) as the group with high expression. According to Table 1, low KLF4 expression was related to a high AFP level, a large tumor, and an advanced TNM stage. Other clinical characteristics in our HCC sample did not correlate significantly with KLF4 expression. Reduced KLF4 expression was significantly related to rapid HCC progression to advanced stages. Figure 2A demonstrated that patients with a higher KLF4 expression had a better OS (HR = 0.24). This difference was also observed in the subgroup of patients with HCC ( Figures  2B-I). In addition, high KLF4 expression was associated with an improved RFS prognosis in HCC (HR = 0.27) ( Figure 2J). Analysis of subgroups revealed that highly expressed KLF4 had a longer RFS ( Figures 2K-R). According to the results, KLF4 may have farreaching effects on the prognosis prediction of HCC.

KLF4 clinically unique value assessment
The results of univariate and multivariate analyses of OS and RFS were used to evaluate the clinical significance of KLF4 expression in HCC patients (Table 2). Age, liver cirrhosis, tumor number, tumor size, AFP level, tumor differentiation grade, TNM stage, and KLF4 expression were all significantly associated with OS and RFS in the univariate analysis. On the contrary hand, gender, HBsAg, tumor distribution, and Child-Pugh grade were not associated with OS or RFS. In addition, in multivariate analysis, only three of 8 prognostic factors remained independent, including KLF4 expression, tumor differentiation, and TNM stage, remained independent prognostic factors. In HCC patients, high KLF4 expression was consistently found to be an independent prognostic predictor of OS (HR = 0.274, 95% CI = 0.152-0.495) and RFS (HR = 0.325, 95% CI = 0.187-0.563).

Establishment of KLF4-based prognostic nomograms for HCC
We developed two nomograms that incorporate KLF4 expression, tumor differentiation, and TNM stage to establish more accurate prognostic models for OS and RFS in HCC patients. In OS and RFS nomograms, a worse prognosis was associated with a higher total sum point of the assigned number of points for each factor (Figures 3A, B). For instance, HCC with low KLF4 expression, tumor stage III, and tumor differentiation grade I/II would have a total of 191 points (91 for KLF4 expression, 0 for tumor differentiation grade I/II, and 100 for tumor stage III), corresponding to a predicted 1-year OS of 42.0%. To report the clinical net benefit of our models, we also developed a conventional prognostic staging system (TNM status). C-indices of OS and RFS in TNM status were 0.710 (0.685-0.735) and 0.706, respectively, demonstrating the prognostic capability of our model (0.682-0.730). The new models based on KLF4 expression had a greater net benefit for predicting OS (Figures 3C, D) and RFS ( Figures 3E, F) than the TNM models for specific threshold probabilities (0.20-0.87) exhibiting a wider range in DCA at 1 and 3 years. The 5-year cut-off probabilities ranged between 0.25 and 0.60 in OS ( Figure 3G) and between 0.25 and 0.75 in RFS ( Figure 3I). Notably, RFS prediction models exacerbated these differences. The accuracy of our nomograms was then evaluated using calibration plots. Calibration plots for the 1-, 3-, and 5-year OS ( Figure 3H) and RFS ( Figure 3J) run very close to the diagonal, indicating that predicted and actual observed outcomes are consistent. Overall, the models we developed could aid clinicians in promoting optimal clinical outcomes for patients with HCC. Frontiers in Genetics frontiersin.org

Signaling pathways related to KLF4
According to Figure 4A, 1,536 genes had a significant upregulated expression in high-KLF4 expression groups (red dots), while 357 genes had a significant downregulated expression (blue dots). Enrichment analysis of all DEGs subsequently revealed that KLF4 was involved in numerous pathways ( Figure 4B). The relationship between KLF4 and immune response-activating cell  Frontiers in Genetics frontiersin.org surface receptor signaling pathway, regulation of lymphocyte immunity, phagocytosis, T cell activation, external side of the plasma membrane, antigen binding, and so on, was revealed by GO and KEGG annotation ( Figures 4C-F). These findings suggested that the KLF4 expression network influenced the TIME in HCC.

KLF4 role to increase immune infiltration
In our study, we found a correlation between KLF4 expression and immune population levels. HCC exhibited significant variations in macrophages, active DC, T cells, T helper cells, and other immune cells ( Figure 5A). Overall, KLF4 expression increased antitumor immune cell infiltration, suggesting KLF4 may play a significant role in antitumor immune response ( Figure 5B). KLF4 was also found to be significantly associated with CD8+T cells, Th1 cells, DC, B cells, NK cells, and macrophages gene markers (Table 3). When measured using immune-related scores, all outcomes remained constant ( Figure 5C). Patients with elevated KLF4 expression had elevated scores for IFN-y response, macrophage regulation, Th1 cells, and lymphocyte infiltration. In the meantime, it was discovered that KLF4 is strongly associated with the genes involved in immune checkpoints, particularly the genes responsible for T-cell exhaustion ( Figures 5D, E). It is widely believed that HCC patients with elevated KLF4 levels are more susceptible to immune checkpoint blockades (ICBs) and that targeting KLF4 could affect exhausted T cell reprogramming to improve immune responses to cancer cells.

KLF4 expression related to macrophages and CD8 + T cells infiltration in HCC
As previously stated, the TIME appeared to have some distinct characteristics in HCC patients with high KLF4 expression. The remarkable correlation between immune cell infiltration levels and KLF4 expression was confirmed by our clinical cohorts. The high KLF4 expression group in HCC tissues contained significantly more CD8 + T cells and macrophages than the low KLF4 expression group. Additionally, CD8 + T cells and macrophage infiltration were positively correlated with KLF4 expression levels ( Figures  6A, B). Notably, only macrophages infiltrated adjacent normal tissues similarly. KLF4 expression was associated with increased macrophage infiltration, but not CD8 + T cell infiltration (Figures 6C,  D). Again, the significant relationship between KLF4 expression and immune infiltration in HCC, specifically macrophages, was evident, suggesting that KLF4 may play crucial roles in the TIME of HCC.

Discussion
Despite the availability of adequate therapeutic options for HCC, the 5-year OS remains dismal at less than 5% (Yang et al., 2019), which is in part due to a lack of biomarkers to precisely predict the prognosis and guide treatment. Although AFP is commonly used for patient risk assessment and surveillance, its sensitivity and specificity are not optimal. For example, several HCC patients have normal AFP levels. Thus, a comprehensive investigation of potential prognostic biomarkers is necessary.
The HCC tumor microenvironment (TME) is a complex and dynamic network composed of cancer cells, stromal cells, chemokines, and extracellular matrix, and the interactions between these components fuel the progression of HCC and influence treatment efficacy. As a tumor-suppressor molecule, KLF4 has been identified in gastrointestinal cancers (Katz et al., 2005). Nevertheless, tumor suppression is insufficient to recapitulate its function, and the precise molecular basis for KLF4 function in the tumor remains obscure. We confirmed that KLF4 expression levels in HCC samples were significantly lower than in non-tumor tissues, and increased KLF4 expression was inversely correlated with advanced tumor grade and tumor stage, indicating that decreased KLF4 expression may promote the progression of HCC. Subgroup The prognostic value of KLF4 suggested that it plays a crucial role in the TME of HCC. Several immune pathways, such as phagocytosis, T cell activation, and other immune responses,

FIGURE 6
Validation of immune cell infiltration in tumors with KLF4 high expression (A) Representative images of KLF4 expression, CD8 + T cells, and macrophage infiltration in HCC tumor tissues with low-/high-KLF4 expression, respectively (B) The CD8 + T cells, and macrophage infiltration level in HCC tumor tissues between low-/high-KLF4 expression group as well as correlation analysis (C) Representative images of KLF4 expression, CD8 + T cells, and macrophage infiltration in HCC adjacent normal tissues with low-/high-KLF4 expression, respectively (D) The CD8 + T cells, and macrophages infiltration level in HCC adjacent normal tissues between low-/high-KLF4 expression group as well as correlation analysis.
Frontiers in Genetics frontiersin.org were found to be associated with TIME in HCC in functional analyses of KLF4, accumulating evidence indicates that KLF4 regulates immune responses. KLF4 expression in DCs is necessary for Th2 cell responses, whereas IFN-γ produced by activated immune cells could increase KLF4 expression (Nagai et al., 2015;Tussiwand et al., 2015). However, it remains unclear how the differential expression of KLF4 affects the infiltration of immune cells. Our study revealed that KLF4 had a significant positive association with macrophages, CD8+T cells, DC cells, and Th1 cells. These findings supported the previous GO and KEGG analyses' findings regarding enriched pathways. According to the results, KLF4 and the TIME were closely related and constantly interacted.
Notably, the TIME of HCC patients with high KLF4 expression was characterized by increased infiltration of CD8 + T cells and macrophages, both of which exhibited a significant positive correlation with increased KLF4 expression. Intratumoral CD8 + T cell infiltration improves patient survival because it is the most important executor of adaptive immunity against HCC (Prieto et al., 2015). Our findings indicated that increased KLF4 expression may be associated with CD8 + T cell infiltration and that this machinery may serve as antitumor immunity in HCC. In addition, we discovered that only a higher degree of macrophage infiltration was associated with increased expression of KLF4 in liver cancer tissues as well as adjacent non-tumor tissues, indicating that KLF4 was likely more important for macrophage infiltration in HCC TIME. Although multiple mechanisms may be responsible for the antitumor effect of macrophages, phagocytosis appears to be the most important one. Consistent with this perception, functional analysis identified phagocytosis as the pathway that was significantly enriched. KLF4 may collectively modulate the antitumor activity of macrophages during phagocytosis, thereby halting the progression of HCC. In this study, we hypothesized that the mechanism underlying KLF4's antitumor effect was more complex than promoting the infiltration of CD8 + T cells or macrophages and that the relationships between KLF4 and immune cells could play a crucial role in TIME in HCC.
Nearly 5 years ago, the landscape of treatments for HCC changed dramatically due to rapid progress in the field of ICBs, including Atezolizumab in 2020 (Lee et al., 2020), Nivolumab in 2020 (Yau et al., 2020), and Pembrolizumab in 2018 (Zhu et al., 2018). The rationale for employing ICBs in HCC is based on the immunosuppressive nature of the HCC TME, which has proved effective treatment for advanced HCC (Huang et al., 2020). In the latest clinical trials RATIONAL 301 study (tislelizumab versus sorafenib for HCC), tislelizumab monotherapy demonstrated positive results with a 14.3% overall response rate (Qin et al., 2022). Yet the precise biomarkers for predicting tumor immunotherapy efficiency are still lacking and the development of prognostic markers is necessary (Galun et al., 2022). Recent research indicates that T cells in a state of exhaustion can serve as medication guidance for tumor immunotherapy (Binnewies et al., 2018). Our findings highlighted that patterns of exhausted T cell hallmarks expression was distinctly tied to KLF4 in HCC, implying that high KLF4 expression may indicate increased sensitivity of HCC to ICBs treatment. The significant tumor suppressive effect of KLF4 and its prognostic value highlighted future application of KLF4 as a therapeutic target or response biomarker in HCC. A small molecule named APTO-253, an inducer of KLF4, has been approved to have anti-tumor activity against cancer cells and in animal models through KLF4 upregulation (Cercek et al., 2015). Phase1 clinical trial has verified the safety, pharmacokinetics and antitumor activity of APTO-253 in adults with advanced solid tumors (e.g., colon cancer, and non-small cell lung cancer) (Cercek et al., 2015). Furthermore, KLF4 could sensitize the therapeutic response to many known drugs, e.g., cetuximab, cisplatin and mesalazine, in the treatment of cancer (Taracha-Wisniewska et al., 2020), suggesting the combined use of APTO-253 with other anticancer drugs in future clinical trials. The present study reflected the positive correlation between KLF4 expression and immune infiltrates as well as immune exhaustion markers, which were indicators for ICBs sensitivity in HCC. Therefore, KLF4 expression might represent a candidate biomarker in predicting the response to ICBs. Besides, the synergistic therapeutic effect of APTO-253 in combination with ICBs also deserves expectation, although further mechanistic studies and clinical validation are needed.
Although our study provided important evidence regarding the value of KLF4 loss in prognostic prediction and correlation with immune infiltrates in TIME, it did have some limitations. First, the hepatitis B virus was the main cause of HCC in this cohort, so whether our results still hold true in patients dominated by other etiologies needs further investigation. Second, our conclusions were derived from a single-center cohort, and future multi-center validation was needed to substantiate our findings. Third, while this study was the first to report the relationship between KLF4 and immune cell infiltration in HCC, the underlying molecular mechanisms remained unknown. With this notion in mind, more detailed analysis was needed to explore how KLF4 affected tumor immunity in HCC.
In conclusion, KLF4 loss significantly correlated with advanced clinical features and independently predicted unfavorable OS and RFS of HCC patients. Our prognostic model incorporating KLF4 expression, tumor differentiation, and TNM stage showed better accuracy than TNM stage alone, which may enable clinicians to make more cost-effective follow-up program and precise treatment decisions for the patients. Our research uncovered a correlation between KLF4 expression, immune cells, and immune checkpoints for the first time. Specifically, CD8+T and macrophage levels were significantly elevated in HCC patients with high KLF4 expression. These results highlighted the unique significance of KLF4 in the prognosis and TIME of HCC. Therefore, stratifying patients further based on their KLF4 expression level would provide greater insight into prognosis, the likelihood of immunotherapy response, and individualized treatment strategies.

Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Frontiers in Genetics frontiersin.org