Comprehensive analysis of KLF2 as a prognostic biomarker associated with fibrosis and immune infiltration in advanced hepatocellular carcinoma
BMC Bioinformatics volume 24, Article number: 270 (2023)
Most Hepatocellular carcinoma (HCC) patients are in advanced or metastatic stage at the time of diagnosis. Prognosis for advanced HCC patients is dismal. This study was based on our previous microarray results, and aimed to explore the promising diagnostic and prognostic markers for advanced HCC by focusing on the important function of KLF2.
The Cancer Genome Atlas (TCGA), Cancer Genome Consortium database (ICGC), and the Gene Expression Comprehensive Database (GEO) provided the raw data of this study research. The cBioPortal platform, CeDR Atlas platform, and the Human Protein Atlas (HPA) website were applied to analyze the mutational landscape and single-cell sequencing data of KLF2. Basing on the results of single-cell sequencing analyses, we further explored the molecular mechanism of KLF2 regulation in the fibrosis and immune infiltration of HCC.
Decreased KLF2 expression was discovered to be mainly regulated by hypermethylation, and indicated a poor prognosis of HCC. Single-cell level expression analyses revealed KLF2 was highly expressed in immune cells and fibroblasts. The function enrichment analysis of KLF2 targets indicated the crucial association between KLF2 and tumor matrix. 33-genes related with cancer associated fibroblasts (CAFs) were collected to identify the significant association of KLF2 with fibrosis. And SPP1 was validated as a promising prognostic and diagnostic marker for advanced HCC patients. CXCR6 CD8+ T cells were noted as a predominant proportion in the immune microenvironment, and T cell receptor CD3D was discovered to be a potential therapeutic biomarker for HCC immunotherapy.
This study identified that KLF2 is an important factor promoting HCC progression by affecting the fibrosis and immune infiltration, highlighting its great potential as a novel prognostic biomarker for advanced HCC.
Hepatocellular carcinoma (HCC), one of the most common and invasive solid malignancies, accounts for the most proportion of liver cancers. The incidence rate of primary liver cancer ranks fifth and the mortality is the third highest, causing approximately 906,000 new cases and over 830,000 deaths per year . Most HCC patients are in advanced or metastatic stage at the time of diagnosis. However, more desperately, the current clinical treatments for advanced HCC patients does not yield good results. Once the tumor has been advanced and metastasized, the HCC patient's prognosis is very bleak . Lymph node is one the most frequent metastasis site for primary HCC. Accordingly, lymph node metastasis occurs in about half of HCC patients with extrahepatic metastasis, and these patients only have a median survival of less than 1 year . Thus, understanding the mechanism of the advancement and metastasis of HCC is essential to propose new therapeutic strategies for patients. And it is imperative to exploit new diagnostic approaches and treatment strategies for advanced HCC patients.
Based on our previous microarray results , we identified several lncRNAs as important factors in HCC progress. Lnc-EPS15L1-2:1 is strongly associated with metastasis and advancement in HCC . However, the molecular regulation of lnc-EPS15L1-2:1 in advanced HCC is still unclear. Accordingly, lncRNAs regulate both coding and noncoding genes via cis- and trans- regulatory signals [6, 7]. Surprisingly, we found that KLF2 was a predictive target gene of lnc-EPS15L1-2:1 from both cis- and trans- analysis, which implied a potential interaction between KLF2 and the progress of HCC.
KLF2, once termed lung Krüppel-like factor (LKLF), is a transcription factor from Krüppel-like factor family. The family genes, characterized by a zinc-finger-containing DNA binding domain, regulate downstream gene transcription via binding to GC-rich DNA sequences . KLF2 is involved in many major biological processes, including proinflammatory activation, cell proliferation, apoptosis, and metabolism (such as glucose metabolism, fatty acid and cholesterol metabolism, amino acid and protein metabolism and so on) [9,10,11,12,13]. In addition, researches have reported that KLF2 is significantly dysregulated in many solid malignancies, including gastric cancer , non-small-cell lung cancer , pancreatic ductal adenocarcinoma and prognostic cancer [16, 17]. Moreover, KLF2 has been demonstrated to be a significantly terminal factor for tumor progress and metastasis in multiple cancers [18,19,20]. Therefore, based on previous results and existing studies, we speculated that KLF2 exert an important biological effect in the lnc-EPS15L1-2:1-related pathway axis to promote HCC progression.
Here, we focus on exploring how KLF2 influences biological changes to promote HCC progression. Interestingly, we found that KLF2 is highly expressed in fibroblasts and immune cells, but low in hepatocytes. Therefore, we further studied the association between cancer associated fibroblasts (CAFs) and KLF2. Results showed that the CAFs-related target gene, SPP1, is a key regulatory gene of KLF2 affecting cancer fibrosis, and is also a promising diagnostic factor for HCC. On the other hand, additionally, the important relationship between immune cells, especially T cells, and KLF2 was also analyzed. We found that an important T-cell receptor (TCR) molecule, CD3D, plays a significant role in KLF2 affecting the immune microenvironment. The results also showed that CD3D has a strong predictive role in HCC immunotherapy. In conclusion, our study provides a novel clue that KLF2 is a considerable contributor for advanced HCC by affecting the fibrosis and immune infiltration, providing new perspectives on exploring the molecular mechanism for HCC advancement, and emphasizing the potential of KLF2 for improving the prognosis of advanced HCC patients in clinical practice.
Materials and methods
Data source and processing
The gene expression profile was derived from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov), the International Cancer Genome Consortium database (ICGC, https://daco.icgc.org/), and the Gene Expression Comprehensive Database (GEO, http://www.ncbi.nlm.nih.gov/geo). TCGA_LIHC and ICGC_LIRI transcriptome data were used for expression correlation analysis and prognostic analysis. GSE25097 and GSE6764 from GEO database were used to explore genes expression distribution. The sequencing data and clinical prognosis data of the HCC multi-cohorts were aggregated by the BEST tool (https://rookieutopia.com/) with batch effects removed to analyze the correlation of gene expression with drug susceptibility, immune infiltration, and prognosis of ICI treatment.
The HCC cohhort, Firehose Legacy cohort, from the cBioPortal platform (https://www.cbioportal.org/) was applied to analyze the mutational landscape and relevant correlation. Single-cell sequencing data were analyzed in CeDR Atlas platform (https://ngdc.cncb.ac.cn/cedr/) and the Human Protein Atlas (HPA, https://www.proteinatlas.org/). The identified KLF2 transcription factor target gene sets were collected from Harmonizome platform (https://maayanlab.cloud/Harmonizome/) to go further analyzing.
For unnormalized RNA-seq data, the raw expression values of genes were log2 transformed. If a gene is traced to multiple probes, the average expression value was taken to representative gene expression levels when using microarray data from the GEO database. Probe entries mapping to unrecorded gene IDs or multiple gene IDs was eliminated. Finally, according to the annotation of the corresponding microarray platform, probe IDs were converted to gene symbols. For the overall expression level of the gene signature, the single sample Gene Set Enrichment Analysis (ssGSEA) algorithm is applied to evaluate the gene enrichment fraction in each sample, thus differentiating high and low expression groups of the gene signature.
Analysis of KLF2 methylation and m6A correlation
N6-methyladenosine (m6A) is an RNA modification that involves the addition of a methyl group to the nitrogen atom at the sixth position (N6) of adenosine. This modification plays a critical role in various biological processes, including mRNA splicing, translation, stability, and degradation. Moreover, it has been implicated in numerous physiological and pathological processes, such as cancer, immune responses, and viral infection [21, 22]. To explore comprehensively the methylation of KLF2 promoter, we used the deoxyribonucleic acid (DNA) methylation data from the online MethSurv tool (https://biit.cs.ut.ee/methsurv/) to analyze different methylation sites of KLF2 and the survival data in TCGA cohort of HCC. The m6A-related genes were derived from the research by Juan Xu et al. on the molecular characterization and clinical significance of m6A modulators across 33 cancer types .
Analysis of differential expression genes and functional enrichment
Limma package in the R software was used to study the differentially expressed mRNA. The threshold was defined as “adjusted P < 0.05 and log2 (Fold Change) > 1 or log2 (Fold Change) < − 1” for the differential expression of mRNAs. For functional enrichment, Gene Ontology (GO), a widely-used tool, is utilized to annotate genes with functions, especially molecular function, biological pathways, and cellular components, and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis , a practical resource, is used to study gene functions and associated high-level genome functional information. ClusterProfiler package (version: 3.18.0) in R was employed to analyze the GO function of potential targets and enrich the KEGG pathway. The R software ggplot2 package was used to draw boxplot and the pheatmap package was used to draw heatmap.
Analysis of immune infiltration and Kaplan–Meier survival
To assess the immune score evaluation, we used immuneeconv, an R software package integrating six latest algorithms, including TIMER, xCell, MCP-counter, CIBERSORT, EPIC and quanTIseq. The immune scoring results were displayed visually through R software package “ggplot2” and “pheatmap”. Kaplan–Meier survival curves were based on RNA sequencing data and corresponding clinical information, which were analyzed and visualized by the “survival” and “surviviner” R packages.
Based on the RNA sequencing data and corresponding clinical information of 371 HCC samples in the TCGA dataset, consistency clustering was performed using the R software package “ConsensesclusterPlus” (v1.54.0), and the parameters were set as follows: the maximum number of clusters was 6, 80% of the total samples were extracted 100 times, clusterAlg = “HC”, innerlinkage = “ward”, D2'. The cluster heatmap was visualized using the “pheatmap” (v1.0.12) R package. Gene expression heatmap was drawn by the “survival” and “surviviner” R packages, and genes with a variance greater than 0.1 were retained.
Screening of prognostic factors and establishment of nomogram
The “forestplot” R package was used to construct the forest map after performed multivariate cox regression analysis. A nomogram was developed based on the results of multivariate cox proportional hazards analysis to predict the 1-, 3-, and 5-year overall recurrence. And the nomogram was established through “rms” R package.
Decreased KLF2 expression is associated with poor survival outcome
Firstly, we analyzed KLF2 expression distribution between different human tumor and normal tissues in TCGA database. The result presented that KLF2 expression is significantly down-regulated in multiple solid cancers (Fig. 1A), such as breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), kidney renal papillary cell carcinoma (KIRP), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), rectum adenocarcinoma (READ) and so on. Additionally, Fig. 1B showed a significant down-regulation of KLF2 expression in the ICGC_LIRI cohort.
In addition, as a supplement, we also explore the expression distribution and prognostic value of several KLF family members (including zinc finger transcription factors) and regulators of NOS enzymes, including KLF2, KLF4, KLF5, KLF6, KLF8, KLF9, KLF10, KLF11, KLF12, NOS2 and NOS3 [25,26,27,28]. The expression distribution of these genes in TCGA database and LIRI databases was initially examined. The results indicated that, as depicted in the Additional file 1: Fig. S1, KLF5 and KLF12 did not exhibit significant differential expression in HCC in either database.
Then survival prognosis of the analyzed genes was assessed using Kaplan–Meier curves in TCGA_LIHC and ICGC- LIRI. Results showed that KLF2 was the only gene found to have a significant correlation with overall survival (OS), progression free survival (PFS), disease free survival (DSS), and disease specific survival (DSS). This means that HCC patients with higher level of KLF2 expression had a significantly better OS, PFS, DFS, and DSS (Fig. 1C–G, Additional file 1: Fig. S2).
Analysis of gene expression regulation of KLF2 in HCC from genomic alteration landscape and methylation modifications
To analyze the expression regulation of KLF2 in HCC comprehensively, we used cBioPortal platform to investigate the genetic mutation status of KLF2. As shown in Fig. 2A, B, KLF2 genomics presented a hypo-mutation condition in most of HCC patients in Firehose Legacy cohort. Only a small percent (1.4%) patients were accompanied by the KLF2 amplification of copy number variation (CNV). The sequence features of KLF2 presented that there were highly conserved classical Cys2/His2 zinc fingers (Fig. 2C). Furthermore, Fig. 2D showed KLF2 expression was negatively correlated with its methylation (Spearman = − 0.57, P < 0.001).
Then, we used Methsurv platform to further analyze methylation level of KLF2 in HCC. The heatmap result showed 10 CpG islands of KLF2 and the corresponding methylation levels. We observed that cg22247553 had the highest level of DNA hypermethylation (Fig. 2E). Besides, we noticed that the methylation level of five CpG sites, cg22247553, cg10819847, cg03725130, cg02668248, and cg18473733, were associated with prognosis (Fig. 2F). And patients with KLF2 hypermethylation of cg22247553, had a worse survival prognosis. Afterwards, we assessed the variance of the expression levels of the m6A-related genes between HCC cancer and normal tissues. As expected, a significantly higher percentage of m6A genes expression was detected in cancer tissues (Fig. 2G). Furthermore, Fig. 2H demonstrated there was a strong linkage between KLF2 exposure and m6A-related genes. These foregoing results indicate that the low expression level of KLF2 in HCC mainly is regulated by hypermethylation, instead of genetic mutation.
Single-cell level expression analysis reveals a high KLF2 expression in fibroblasts and immune cells, while low expression in hepatocytes
To analyze KLF2 expression distribution in different cell types in liver cancer and normal tissue, we used Uniform Manifold Approximation and Projection (UMAP) and t-distributed Stochastic Neighbor Embedding (tSNE) algorithm for the single cell expression analysis in multiple platforms. Figure 3A and B showed the single-cell expression level of KLF2 in different cell lines of the SCP542 analysis  in the Cellular Drug Response (CeDR) Atlas and Single cell portal (SCP) database. Using the CeDR Atlas, we further analyzed the single-cell sequencing data of hepatic cells derived from GSE115469  and GSE130073 . The UMAP and the Cell Fraction plots showed the epithelial cells and fibroblasts make up a substantial proportion among all liver cells, which together accounted for nearly 90% (Fig. 3C, D). Then we used HPA platform to further analyze the single-cell expression level of KLF2 in liver tissues. The bar graph results exhibited that KLF2 was mostly expressed in fibroblasts, epithelial cells and immune cells (T-cells, B-cells, plasma cell, NK-cells and so on) (Fig. 3E, F). Surprisingly, we noticed a low KLF2 expression in hepatocytes. These results suggest that KLF2 expression levels are significantly different in different cell types in the liver tissue, mainly focusing on fibroblasts, epithelial cells and immune cells, which are the main components of liver tissue.
Functional enrichment analysis of KLFTs
KLF2 is an important transcription factor from the Krüppel-like family of proteins, which exerts regulation in various cell types in the processes of activation, differentiation and migration by targeting different downstream genes. We collected 98 target genes of the KLF2 from the CHEA Transcription Factor Targets dataset in Harmonizome platform, namely KLFTs (Additional file 2: Table S1). The single sample Gene Set Enrichment Analysis (ssGSEA) was used to calculate absolute enrichment scores, and Gene Set Enrichment Analysis (GSEA) was used to demonstrate the biological processes. We observed that the top several enrichment entries were concentrated in the biological characteristics related to tumor matrix and immunity. For GO analysis shown in Fig. 4A, the major items were “External encapsulating structure organization”, “Collagen fibril organization”, “Lymph vessel development”. For KEGG analysis in Fig. 4B, enriched items mainly were “Focal adhesion”, “ECM receptor adhesion”, “Leishmania infection”. Furthermore, we examined the correlation between KLF2 and common signaling pathways associated with tumor matrix formation. As we expected, KLF2 was in a high positive correlation with these pathways, such as EMT_markers, ECM_related genes, collagen_ formation, degradation_of_ECM and TGF β (Fig. 4C). EMT is primarily executed by EMT-activating transcription factors that belong to families such as SNAIL, TWIST, and ZEB and so on [32, 33]. Subsequently, a correlation analysis was conducted between KLF2 and EMT-markers, such as SNAI1, SNAI2, TWIST1, TWIST2, ZEB1, ZEB2, VIM and ID1. As depicted in Additional file 1: Fig. S3A–H, the study findings revealed a positive association between KLF2 expression levels and these EMT markers.
Fibroblasts account for the major stromal cell type in the microenvironment of liver diseases, including liver cirrhosis and liver cancers. Cirrhosis is a predominant contributor to the development of HCC. Therefore, we further analyzed the distribution of KLF2 expression in two datasets related with cirrhosis development, GSE 25097 and GSE 6764, to investigate the different expression of KLF2 in the tumor tissue compared to cirrhosis tissue. In Fig. 4D, the results analyzed from GSE 25097 showed KLF2 expression level continued to decreasing, as liver tissue became progressively cirrhotic and then progressed to HCC (G1: normal liver tissue; G2: liver cirrhosis tissue; G3: HCC tissue). The identical results were also observed in GSE 6764 in Fig. 4E (G1: liver cirrhosis tissue; G2: early HCC tissue; G3: advanced HCC tissue). In addition, we examined the expression distribution of SNAI1, ZEB2, and VIM, which are the top three markers exhibiting strong correlation with KLF2, in the GSE25097 dataset related to cirrhosis development. The findings revealed that normal tissues exhibited higher expression levels of SNAI1, ZEB2, and VIM compared to cancer tissues, and cirrhotic tissues showed higher expression levels than cancerous tissues as well (Additional file 1: Fig. S3I–K). The above results suggest that KLF2 is involved in the regulation of biological processes associated with tumor matrix. These results further demonstrate an essential role that KLF2 performed during development of HCC associated with liver fibrosis/cirrhosis.
Identification a prognostic signature based on the KLFTs
In order to identify the key target genes of KLF2 in HCC, we successively performed the co-expression analysis and univariate Cox regression analysis of KLFTs (Additional file 3: Table S2, Additional file 1: Fig. S4). Then we obtained 17 prognostic significant genes, named KLFTs-17. Figure 4F showed that these 17 targets were over-expressed in HCC tumor tissues compared to the normal tissues.
Fibrosis of the tumor stroma is an important biological process contributing to the advancement of solid tumors. And cancer associated fibroblasts (CAFs) are the primary cells involved in this process. An increasing number of studies have reported that CAFs play a crucial role in promoting solid tumorigenesis. To further investigate the association between KLFTs and CAFs, we first collected 33 common CAFs-related marker genes by consulting relevant references [34,35,36,37] (Additional file 4: Table S3). Figure 4G showed the strong association of KLF2 with CAFs-related marker genes. Furthermore, we selected three KLF2 target genes (PPP1R12A, SP1 and SMARCAD1) with the strongest association with CAFs-related genes among KLFTs for expression level analysis in the GSE25097 and GSE6764 (Fig. 4H, I). Obviously, the expressions of PPP1R12A, SP1 and SMARCAD1 in HCC tissues were higher than those in the liver cirrhosis tissues, which were opposite to the expression trend of KLF2 in HCC. The above results suggest that KLF2 negatively regulates the expression of those downstream genes associated with fibrosis.
Multidimensional analysis of subgroups classified by KLFTs-17
We used cumulative distribution function (CDF) based on consensus clustering and classfied TCGA_LIHC samples into two subtype groups with k = 2 selected as the optimum (Fig. 5A–C). The cluster heatmap in Fig. 5D showed that the expression level of genes of the subgroup C1 were higher than that in the subgroup C2. Then we explored the expression distribution of m6A-related marker genes between subgroups. The results showed that the gene expression level differed distinctly between subgroups, and the expression level was higher in C1 subtype (Fig. 5E).
In addition, Kaplan–Meier survival analysis of subgroups presented those patients with higher KLFTs-17 expression achieved the worse OS (Fig. 5F). The clinical characteristics comparison analysis showed that the subgroup C1 was correlated to the more advanced T stage and higher histologic grade (Fig. 5G). The heatmap and volcano maps showed the differential gene expression (DEGs) between subgroups (Fig. 6H, I). Specifically, the functional enrichment analysis in Fig. 5J indicated that enriched terms were related to infection, cell proliferation, oncogenesis and so on. For GO analysis, the DEGs-up were highlighted in the entries “organelle fission”, “nuclear division”, “extracellular structure organization”, “extracellular matrix organization” and so on. However, the DEGs-down were mostly enriched in biological process related with metabolism. Taken together, the results above listed imply that highly expressed genes of subgroup classified by KLFTs-17 is more closely related to tumorigenesis and progression.
Screening and identification of CAFs-related prognostic marker in HCC
To further derive promising markers associated with CAFs in HCC progression, we selected the C1 subgroup with the most striking features of KLFTs-17 to analyze the important role of CAF. COX regression analysis was applied, and SPP1 was identified as an independent prognostic factor associated with HCC fibrosis (Additional file 5: Table S4, Fig. 6A). Additionally, we depicted a nomogram to help predict the 1-, 3-, and 5-year survival of HCC patients (P < 0.001, C-index = 0.632, 95% CI 0.563–1) (Fig. 6B). The calibration curve in Fig. 6C showed a fine prognostic prediction performance of the nomogram established.
To further demonstrate the role of SPP1 in HCC, we explored the distribution of SPP1 expression across different types of tumor and normal tissues in TCGA. The result showed that SPP1 was increased not only in LIHC tumor tissue, but also in multiple cancers (Fig. 6D), such as BRCA, COAD, KIRP, LUAD, LUSC, READ and so on. To explore the molecular biological function of SPP1 in HCC, we used GSEA analysis to enrich the SPP1-related genes. Intriguingly, Fig. 6E revealed that SPP1 was involved not only in tumor stroma-related pathways, but also in the immune-related biological processes, such as “Allograft rejection” and “Inflammatory response”.
Furthermore, we analyzed the association between SPP1 expression and clinical characteristics. The results in Fig. 6F and G indicated a significant interaction of SPP1 expression with the tumor Grade and T stage. Additionally, we assessed the prognostic value of SPP1 in TCGA_LIHC and ICGC_LIRI. The KM curves presented HCC patients with higher level of SPP1 expression had a significantly worse OS, DSS, PFS, and RFS (Fig. 6H–K). What’s more, we noticed that SPP1, also known as osteopontin (OPN), had been reported to be a promising tumor marker for detecting metastatic disease in many tumors [38, 39]. Therefore, we assessed the diagnosis value of SPP1 in TCGA_LIHC, and the result of ROC curve showed SPP1 had a good sensitivity and high specificity for the diagnosis of HCC (AUC 0.732; CI 0.684–0.779) (Fig. 6L,M). The foregoing results support that patients with increased SPP1 expression have a poor prognosis. And SPP1 shows a favorable ability for HCC diagnosis.
Relationship of KLF2 expression with immune infiltration in HCC
The contribution of the tumor microenvironment to tumor prognosis cannot be negligible. ESTIMATE method was utilized to estimate non-tumor cell infiltration level involved in tumor microenvironment. The findings indicated that KLF2 expression was positively associated with Immune Score, Stromal Score and Estimate Score (Fig. 7A). Then we further investigated correlation of KLF2 expression with different types of immune cells. Scatter plots results showed a significantly positive correlation between KLF2 expression and immune cells infiltration in HCC based on the TIMER algorithm (Fig. 7B).
Immune checkpoint molecules expressed on the immune cells play a critical role in the immune responses and immunotherapy. Figure 7C showed the significant association between KLF2 and immune checkpoint molecules. Moreover, KLF2 expression was negatively correlated with tumor mutational burden (TMB) (ρ = − 0.29, P < 0.001) (Fig. 7D). Immune checkpoint inhibits (ICIs) therapy have dramatically improved outcomes of cancer patients in clinical practice, especially for PD-1 and PD-L1 (CD274). The scatter plots revealed the positive relationship of KLF2 with PD-1 and PD-L1 (CD274) (Fig. 7E, F). Then, we assessed the responsive of KLF2 in the immunotherapy cohorts, IMvigor210 cohort 2018 and Kim cohort 2019. Boxplots showed that patients with low KLF2-expressing level had a stronger response to these ICIs (Fig. 7G, H). Prognostic KM curves showed that patients with low KLF2 expression achieved better prognostic survival after ICIs treatments (Fig. 7I, J). Taken together, the above results suggest that KLF2 is highly correlated with the immune microenvironment in HCC, and advanced HCC patients with lower KLF2 expression derive greater benefit from ICIs treatments.
Identification and analysis of a therapeutic biomarker CD3D
To better understand the association between KLF2 and heterogeneity of immune infiltration cells in liver tissue, we further analyzed the single-cell sequencing data of liver-resident immune cells derived from GSE125188  using the CeDR Atlas database. As shown in Fig. 8A, the UMAP and the Cell Fraction plots showed the clustering of cell types and the proportion of cell types. The results showed that T cells accounted for the major proportion among all immune infiltration cells in the liver tissue. In addition, as Fig. 3E and F shown, KLF2 expression was higher in the T cells than in other immunocytes. In the anti-tumor immune response, CD8+ cytotoxic T lymphocytes (CTLs) played a main cellular effector role. Accordingly, we selected CXCR6 CD8+ T cells for further analysis. The heatmap in Fig. 8C indicated a strong association of KLF2 with the main marker genes expressed in CXCR6 CD8+ T cells . As shown in Fig. 8B, a network of cell types and drug response was constructed, and for CXCR6 CD8+ T cells, the drug Isoflupredone showed the most significant statistical significance (Spearman cor: − 0.25: P value: 0.042). Then marker genes of CXCR6 CD8+ T cells and differential genes expressed by the drug Isoflupredone-induced in the CXCR6 CD8+ T cells were analyzed applying GSEA method (Fig. 8C–E). As Fig. 8F presented, the matrix plot visualized the main signature genes involved in associations between Isoflupredone and CXCR6 CD8+ T cells, with CD3D highlighted.
CD3D, an essential part of the T-cell receptor/CD3 complex (TCR/CD3 complex), is involved in development and signal transduction of T-cell . Therefore, we further analyzed the association of CD3D with other TCRs. Scatter plots showed CD3D was statistically correlated with ICs, PDCD1 and CD274 (PDLD1) (Fig. 8G, H). Moreover, GSEA analysis was applied to enrich CD3D-related genes in HCC to better understand the molecular biological function of CD3D in HCC. Figure 8I showed the enriched entries were mainly immune-related biological processes. However, we noticed that the pathway “epithelial mesenchymal transition” was also included in the top ten items. Surprisingly, when comparing the enrichment results of SPP1-related genes and CD3D-related genes, we found that nearly half of the two pathways overlapped in the top ten enriched entries (Figs. 6E, 8I).
What’s more, Fig. 8J showed an essential distinction of CD3D expression in different Grades of HCC patients. And as tumor Grade advanced, CD3D expression levels were increasingly up-regulated. We noticed that the higher expression of CD3D was more responsive to ICIs (Fig. 8K, L). KM prognostic curves showed that higher CD3D expression predicted a better survival in HCC patients after ICIs treatments (Fig. 8M, N). Furthermore, area under curve (AUC) of ROC curves revealed CD3D had a moderate capacity to predicate HCC patients’ response to ICB therapy (Fig. 8O, P).
In our study, we focus on the exploration of KLF2, a predictive target gene of lnc-EPS15L1-2:1, which has been identified to have a strong association with HCC advancement in our previous research . Surprisingly, we discovered that methylation, rather than genetic mutations is responsible for the down regulation of KLF2 expression in HCC. In addition, it has been reported that lncRNA ANRIL can regulate cell growth in vitro and in vivo through epigenetic silencing of KLF2 in HCC . KM survival analysis reveals that HCC patients with decreased KLF2 expression tend to achieve a much worse OS, DSS, PFS, and RFS. This highlights the strong influence of KLF2 on the prognosis of HCC patients and its potential use as a dependable marker for prognostication.
KLF2, an important transcription factor, participates in many biological processes. The function enrichment results of the KLF2 targets, namely KLFTs, indicate that highly enriched entries focus on biological features associated with tumor stroma and immune responses. And we also assessed the differential levels of KLF2 expressed in cirrhosis tissue and HCC tissue. The findings further demonstrate that the expression level of KLF2 decreased continuously during the process of liver fibrosis to cirrhosis and then to hepatocarcinogenesis. Several studies have reported that KLF2 plays an important role in maintaining hepatic endothelial cell homeostasis and vascular integrity, and protects the liver from fibrosis or cirrhosis [43, 44]. TGF-β, a strong inducer of EMT, has been implicated as a key cytokine mediating liver fibrosis . TGF-β is able to promote both fibrosis and carcinogenesis and shows increased levels in cirrhosis and advanced HCC . KLF2 has been shown to function as a tumor suppressor through TGF-β/Smad signaling in HCC cells . What’s more, other members of KLF family have been reported to significantly regulate the fibrotic process by transcriptionally regulating TGF-β expression, such as KLF4, KLF5, KLF6, and KLF15 [47,48,49,50], which indicates the great potential of KLF2 as a key gene involved in HCC fibrosis.
For most solid tumors, tumor matrix is an indispensable factor to promote tumor progression. Cancer associated fibroblasts (CAFs) are the predominant cells to improve tumor stromal microenvironment to facilitate tumor outgrowth. A growing number of studies have reported the essential contribution of CAFs to tumor progression [51, 52]. Most importantly, our study identifies a CAFs-related marker, SPP1, which is strongly associated with clinical characteristics and prognostic survival of HCC patients. Osteopontin (OPN), encoded by SPP1, has been implicated in multiple human diseases and has been shown to play an important regulatory role in HCC progression [53,54,55]. Therefore, we believe that the KLF2-SPP1 pathway is an important signaling axis that promotes liver tissue fibrosis, thus leading to HCC progression. However, whether SPP1 is a potential target of KLF2 involved in fibrosis deserves to be investigated experimentally.
OPN has been reported to be a promising tumor marker for detection in the diagnosis of many tumor progressions [38, 39]. And the ROC result reveals the high sensitivity and specificity of SPP1 diagnosis value in HCC, indicting the strong capacity of SPP1 as a prospective biomarker for future diagnosis and prognosis prediction for HCC patients.
Multidimensional analyses of the correlation between KLF2 and immune infiltration showed an essential involvement of KLF2 in regulating tumor immune microenvironment. We found that patients with low KLF2 expression levels were more responsive to ICIs therapy and achieved a better prognosis and survival, indicating that advanced HCC patients with lower KLF2 expression levels are more suitable for ICIs therapy.
The results of single-cell level expression analysis from multiple platforms reveal the expression distribution of KLF2 is highly distinct in different cells of liver tissue. High levels of expression in immune cells (especially T cells), endothelial cells, and fibroblasts, but low levels in hepatocytes. Kuo et al. have demonstrated that KLF2 is expressed in both CD4+ and CD8+ T cells. The expression of KLF2 mRNA and protein is significantly down-regulated upon activation of resting T cells via T cell receptor (TCR) , which suggests an essential role of KLF2 on T cell function. And KLF2 has also been shown to inhibit the proliferation and growth of Jurkat T leukemia cells [57, 58]. These existed researches have suggested that KLF2 plays a critical role in maintaining the function of T cells. CD8+ T cells performs important functions in the immune response, and cluster analysis shows CXCR6 CD8+ T cell is the predominant T cell type in liver tissue. And KLF2 is significantly associated with biomarkers of CXCR6 CD8+ T cell. In addition, we found that CD3D, an important membrane protein for CXCR6 CD8+ T cells to exert immune responsive, was associated with HCC progression and immunotherapy. Furthermore, many studies have reported CD3D as a promising prognostic and therapeutic biomarker [59,60,61]. These studies further confirmed our findings that CD3D is prospective to be an emerging marker for HCC immunotherapy. Therefore, we speculate CD3D is a key mediator of KLF2 involvement in HCC immune response.
Our study identifies the important function of KLF2 for advanced HCC by affecting the fibrosis and immune infiltration, and provides new perspectives on exploring the molecular mechanism for HCC advancement, emphasizing the potential of KLF2 as a new biomarker for improving the prognosis of advanced HCC patients in clinical practice.
Availability of data and materials
The data of this study are available in the The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov), the International Cancer Genome Consortium database (ICGC, https://daco.icgc.org/), the Gene Expression Comprehensive Database (GEO, http://www.ncbi.nlm.nih.gov/geo), cBioPortal platform (https://www.cbioportal.org/), CeDR Atlas platform (https://ngdc.cncb.ac.cn/cedr/), the Human Protein Atlas (HPA, https://www.proteinatlas.org/), and Harmonizome platform (https://maayanlab.cloud/Harmonizome/).
Differentially expressed genes
The Cancer Genome Atlas
The gene expression comprehensive database
International cancer genome consortium
Human Protein Atlas
Relapse free survival
Progression free survival
Disease free survival
Breast invasive carcinoma
Kidney renal papillary cell carcinoma
Lung squamous cell carcinoma
Cancer associated fibroblasts
Receiver operating characteristic
Area under curve
Gene set enrichment analysis
Kyoto Encyclopedia of Genes and Genomes
Secreted phosphoprotein 1
T cell receptor
Tumor mutational burden
Immune checkpoint inhibits
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: globocan estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49. https://doi.org/10.3322/caac.21660.
Aino H, Sumie S, Niizeki T, Kuromatsu R, Tajiri N, Nakano M, et al. Clinical characteristics and prognostic factors for advanced hepatocellular carcinoma with extrahepatic metastasis. Mol Clin Oncol. 2014;2(3):393–8. https://doi.org/10.3892/mco.2014.259.
Xia F, Wu L, Lau WY, Li G, Huan H, Qian C, et al. Positive lymph node metastasis has a marked impact on the long-term survival of patients with hepatocellular carcinoma with extrahepatic metastasis. PLoS ONE. 2014;9(4):e95889. https://doi.org/10.1371/journal.pone.0095889.
Ma J, Zhang L, Yang P, Zeng ZC, Xiang ZL. Integrated analysis of long noncoding RNA expression profiles in lymph node metastasis of hepatocellular carcinoma. Gene. 2018;676:47–55. https://doi.org/10.1016/j.gene.2018.07.001.
Ma J, Zhang L, Bian HR, Lu ZG, Zhu L, Yang P, et al. A noninvasive prediction nomogram for lymph node metastasis of hepatocellular carcinoma based on serum long noncoding RNAs. Biomed Res Int. 2019;2019:1710670. https://doi.org/10.1155/2019/1710670.
Engreitz JM, Haines JE, Perez EM, Munson G, Chen J, Kane M, et al. Local regulation of gene expression by lncRNA promoters, transcription and splicing. Nature. 2016;539(7629):452–5. https://doi.org/10.1038/nature20149.
Kopp F, Mendell JT. Functional classification and experimental dissection of long noncoding RNAs. Cell. 2018;172(3):393–407. https://doi.org/10.1016/j.cell.2018.01.011.
Huang Z, He H, Qiu F, Qian H. Expression and prognosis value of the KLF family members in colorectal cancer. J Oncol. 2022;2022:6571272. https://doi.org/10.1155/2022/6571272.
Black AR, Black JD, Azizkhan-Clifford J. Sp1 and Krüppel-like factor family of transcription factors in cell growth regulation and cancer. J Cell Physiol. 2001;188(2):143–60. https://doi.org/10.1002/jcp.1111.
Manavski Y, Abel T, Hu J, Kleinlützum D, Buchholz CJ, Belz C, et al. Endothelial transcription factor KLF2 negatively regulates liver regeneration via induction of activin A. Proc Natl Acad Sci U S A. 2017;114(15):3993–8. https://doi.org/10.1073/pnas.1613392114.
SenBanerjee S, Lin Z, Atkins GB, Greif DM, Rao RM, Kumar A, et al. KLF2 is a novel transcriptional regulator of endothelial proinflammatory activation. J Exp Med. 2004;199:1305–15.
Yerra VG, Drosatos K. Specificity proteins (SP) and Krüppel-like factors (KLF) in liver physiology and pathology. Int J Mol Sci. 2023;24:4682. https://doi.org/10.3390/ijms24054682.
Lin J, Tan H, Nie Y, Wu D, Zheng W, Lin W, et al. Krüppel-like factor 2 inhibits hepatocarcinogenesis through negative regulation of the hedgehog pathway. Cancer Sci. 2019;110:1220–31. https://doi.org/10.1111/cas.13961.
Xu TP, Liu XX, Xia R, Yin L, Kong R, Chen WM, et al. SP1-induced upregulation of the long noncoding RNA TINCR regulates cell proliferation and apoptosis by affecting KLF2 mRNA stability in gastric cancer. Oncogene. 2015;34(45):5648–61. https://doi.org/10.1038/onc.2015.18.
Yin L, Wang JP, Xu TP, Chen WM, Huang MD, Xia R, et al. Downregulation of Kruppel-like factor 2 is associated with poor prognosis for nonsmall-cell lung cancer. Tumour Biol J Int Soc Oncodev Biol Med. 2015;36(4):3075–84. https://doi.org/10.1007/s13277-014-2943-4.
Zhang D, Dai Y, Cai Y, Suo T, Liu H, Wang Y, et al. KLF2 is downregulated in pancreatic ductal adenocarcinoma and inhibits the growth and migration of cancer cells. Tumour Biol J Int Soc Oncodev Biol Med. 2016;37(3):3425–31. https://doi.org/10.1007/s13277-015-4053-3.
Wang B, Liu M, Song Y, Li C, Zhang S, Ma L. KLF2 inhibits the migration and invasion of prostate cancer cells by downregulating mmp2. Am J Men’s Health. 2019;13(1):1817094795. https://doi.org/10.1177/1557988318816907.
Ma Z, Peng P, Zhou J, Hui B, Ji H, Wang J, et al. Long non-coding RNA SH3PXD2A-AS1 promotes cell progression partly through epigenetic silencing P57 and KLF2 in colorectal cancer. Cell Physiol Biochem Int J Exp Cell Physiol Biochem Pharmacol. 2018;46(6):2197–214. https://doi.org/10.1159/000489589.
Xia E, Bhandari A, Shen Y, Zhou X, Wang O. lncRNA LINC00673 induces proliferation, metastasis and epithelial-mesenchymal transition in thyroid carcinoma via Kruppel-like factor 2. Int J Oncol. 2018;53(5):1927–38. https://doi.org/10.3892/ijo.2018.4524.
Lian Y, Yan C, Xu H, Yang J, Yu Y, Zhou J, et al. A novel lncRNA, LINC00460, affects cell proliferation and apoptosis by regulating KLF2 and CUL4A expression in colorectal cancer. Mol Ther Nucleic Acids. 2018;12:684–97. https://doi.org/10.1016/j.omtn.2018.06.012.
He PC, He C. m6A RNA methylation: from mechanisms to therapeutic potential. EMBO J. 2021;40:e105977. https://doi.org/10.15252/embj.2020105977.
He L, Li H, Wu A, Peng Y, Shu G, Yin G. Functions of n6-methyladenosine and its role in cancer. Mol Cancer. 2019;18:1–15. https://doi.org/10.1186/s12943-019-1109-9.
Li Y, Xiao J, Bai J, Tian Y, Qu Y, Chen X, et al. Molecular characterization and clinical relevance of m6A regulators across 33 cancer types. Mol Cancer. 2019;18(1):137. https://doi.org/10.1186/s12943-019-1066-3.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28:1947–51. https://doi.org/10.1002/pro.3715.
McConnell BB, Yang VW. Mammalian Krüppel-like factors in health and diseases. Physiol Rev. 2010;90:1337–81. https://doi.org/10.1152/physrev.00058.2009.
Kim C, He P, Bialkowska AB, Yang VW. SP and KLF transcription factors in digestive physiology and diseases. Gastroenterology. 2017;152:1845–75. https://doi.org/10.1053/j.gastro.2017.03.035.
Wang R, Geller DA, Wink DA, Cheng B, Billiar TR. No and hepatocellular cancer. Br J Pharmacol. 2020;177:5459–66. https://doi.org/10.1111/bph.14838.
Chen T, Shi Z, Zhao Y, Meng X, Zhao S, Zheng L, et al. LncRNA Airn maintains LSEC differentiation to alleviate liver fibrosis via the KLF2-eNOS-sGC pathway. BMC Med. 2022;20:335. https://doi.org/10.1186/s12916-022-02523-w.
Madissoon E, Wilbrey-Clark A, Miragaia RJ, Saeb-Parsy K, Mahbubani KT, Georgakopoulos N, et al. scRNA-seq assessment of the human lung, spleen, and esophagus tissue stability after cold preservation. Genome Biol. 2019;21(1):1. https://doi.org/10.1186/s13059-019-1906-x.
MacParland SA, Liu JC, Ma XZ, Innes BT, Bartczak AM, Gage BK, et al. Single cell RNA sequencing of human liver reveals distinct intrahepatic macrophage populations. Nat Commun. 2018;9(1):4383. https://doi.org/10.1038/s41467-018-06318-7.
Ouchi R, Togo S, Kimura M, Shinozawa T, Koido M, Koike H, et al. Modeling steatohepatitis in humans with pluripotent stem cell-derived organoids. Cell Metab. 2019;30(2):374–84. https://doi.org/10.1016/j.cmet.2019.05.007.
Brabletz T, Kalluri R, Nieto MA, Weinberg RA. EMT in cancer. Nat Rev Cancer. 2018;18:128–34. https://doi.org/10.1038/nrc.2017.118.
Giannelli G, Koudelkova P, Dituri F, Mikulits W. Role of epithelial to mesenchymal transition in hepatocellular carcinoma. J Hepatol. 2016;65:798–808. https://doi.org/10.1016/j.jhep.2016.05.007.
Kalluri R. The biology and function of fibroblasts in cancer. Nat Rev Cancer. 2016;16(9):582–98. https://doi.org/10.1038/nrc.2016.73.
Affo S, Yu LX, Schwabe RF. The role of cancer-associated fibroblasts and fibrosis in liver cancer. Annu Rev Pathol. 2017;12:153–86. https://doi.org/10.1146/annurev-pathol-052016-100322.
Zou B, Liu X, Gong Y, Cai C, Li P, Xing S, et al. A novel 12-marker panel of cancer-associated fibroblasts involved in progression of hepatocellular carcinoma. Cancer Manag Res. 2018;10:5303–11. https://doi.org/10.2147/CMAR.S176152.
Baglieri J, Brenner DA, Kisseleva T. The role of fibrosis and liver-associated fibroblasts in the pathogenesis of hepatocellular carcinoma. Int J Mol Sci. 2019. https://doi.org/10.3390/ijms20071723.
Fedarko NS, Jain A, Karadag A, Van Eman MR, Fisher LW. Elevated serum bone sialoprotein and osteopontin in colon, breast, prostate, and lung cancer. Clin Cancer Res Off J Am Assoc Cancer Res. 2001;7(12):4060–6.
Reiniger IW, Wolf A, Welge-Lüssen U, Mueller AJ, Kampik A, Schaller UC. Osteopontin as a serologic marker for metastatic uveal melanoma: results of a pilot study. Am J Ophthalmol. 2007;143(4):705–7. https://doi.org/10.1016/j.ajo.2006.11.040.
Zhao J, Zhang S, Liu Y, He X, Qu M, Xu G, et al. Single-cell RNA sequencing reveals the heterogeneity of liver-resident immune cells in human. Cell Discov. 2020;6:22. https://doi.org/10.1038/s41421-020-0157-z.
Smith-Garvin JE, Koretzky GA, Jordan MS. T cell activation. Annu Rev Immunol. 2009;27:591–619. https://doi.org/10.1146/annurev.immunol.021908.132706.
Huang M, Chen W, Qi F, Xia R, Sun M, Xu T, et al. Erratum to: long non-coding RNA ANRIL is upregulated in hepatocellular carcinoma and regulates cell proliferation by epigenetic silencing of KLF2. J Hematol Oncol. 2017;10:1–1. https://doi.org/10.1186/s13045-017-0513-0.
Gracia-Sancho J, Russo L, García-Calderó H, García-Pagán JC, García-Cardeña G, Bosch J. Endothelial expression of transcription factor Kruppel-like factor 2 and its vasoprotective target genes in the normal and cirrhotic rat liver. Gut. 2011;60:517–24. https://doi.org/10.1136/gut.2010.220913.
Marrone G, Russo L, Rosado E, Hide D, García-Cardeña G, García-Pagán JC, et al. The transcription factor KLF2 mediates hepatic endothelial protection and paracrine endothelial-stellate cell deactivation induced by statins. J Hepatol. 2013;58:98–103. https://doi.org/10.1016/j.jhep.2012.08.026.
Bataller R, Brenner DA. Liver fibrosis. J Clin Investig. 2005;115:209–18.
Li Y, Tu S, Zeng Y, Zhang C, Deng T, Luo W, et al. KLF2 inhibits TGF-β-mediated cancer cell motility in hepatocellular carcinoma. Acta Biochim Biophys Sin. 2020;52:485–94. https://doi.org/10.1093/abbs/gmaa024.
Chen WC, Lin HH, Tang MJ. Matrix-stiffness-regulated inverse expression of Krüppel-like factor 5 and Krüppel-like factor 4 in the pathogenesis of renal fibrosis. Am J Pathol. 2015;185(9):2468–81. https://doi.org/10.1016/j.ajpath.2015.05.019.
King KE, Iyemere VP, Weissberg PL, Shanahan CM. Krüppel-like factor 4 (KLF4/GKLF) is a target of bone morphogenetic proteins and transforming growth factor beta 1 in the regulation of vascular smooth muscle cell phenotype. J Biol Chem. 2003;278(13):11661–9. https://doi.org/10.1074/jbc.M211337200.
Holian J, Qi W, Kelly DJ, Zhang Y, Mreich E, Pollock CA, et al. Role of Kruppel-like factor 6 in transforming growth factor-beta1-induced epithelial-mesenchymal transition of proximal tubule cells. Am J Physiol Renal Physiol. 2008;295(5):F1388–96. https://doi.org/10.1152/ajprenal.00055.2008.
Gao X, Wu G, Gu X, Fu L, Mei C. Kruppel-like factor 15 modulates renal interstitial fibrosis by ERK/MAPK and JNK/MAPK pathways regulation. Kidney Blood Press Res. 2013;37(6):631–40. https://doi.org/10.1159/000355743.
Yang F, Tuxhorn JA, Ressler SJ, McAlhany SJ, Dang TD, Rowley DR. Stromal expression of connective tissue growth factor promotes angiogenesis and prostate cancer tumorigenesis. Cancer Res. 2005;65(19):8887–95. https://doi.org/10.1158/0008-5472.CAN-05-1702.
Orimo A, Gupta PB, Sgroi DC, Arenzana-Seisdedos F, Delaunay T, Naeem R, et al. Stromal fibroblasts present in invasive human breast carcinomas promote tumor growth and angiogenesis through elevated sdf-1/cxcl12 secretion. Cell. 2005;121(3):335–48. https://doi.org/10.1016/j.cell.2005.02.034.
Icer MA, Gezmen-Karadag M. The multiple functions and mechanisms of osteopontin. Clin Biochem. 2018;59:17–24. https://doi.org/10.1016/j.clinbiochem.2018.07.003.
Ye QH, Qin LX, Forgues M, He P, Kim JW, Peng AC, et al. Predicting hepatitis B virus-positive metastatic hepatocellular carcinomas using gene expression profiling and supervised machine learning. Nat Med. 2003;9(4):416–23. https://doi.org/10.1038/nm843.
Wang J, Hao F, Fei X, Chen Y. SPP1 functions as an enhancer of cell growth in hepatocellular carcinoma targeted by mir-181c. Am J Transl Res. 2019;11(11):6924–37.
Kuo CT, Veselits ML, Leiden JM. LKLF: a transcriptional regulator of single-positive t cell quiescence and survival. Science (New York, NY). 1997;277(5334):1986–90. https://doi.org/10.1126/science.277.5334.1986.
Wu J, Lingrel JB. KLF2 inhibits Jurkat T leukemia cell growth via upregulation of cyclin-dependent kinase inhibitor p21WAF1/CIP1. Oncogene. 2004;23(49):8088–96. https://doi.org/10.1038/sj.onc.1207996.
Buckley AF, Kuo CT, Leiden JM. Transcription factor LKLF is sufficient to program T cell quiescence via a c-Myc-dependent pathway. Nat Immunol. 2001;2(8):698–704. https://doi.org/10.1038/90633.
Zhu Z, Ye W, Wu X, Lin S, Xu J, Li L, et al. Comprehensive analysis reveals a prognostic and therapeutic biomarker CD3D in the breast carcinoma microenvironment. 2021. Biosci Rep. https://doi.org/10.1042/BSR20202898.
Yang Y, Zang Y, Zheng C, Li Z, Gu X, Zhou M, et al. CD3D is associated with immune checkpoints and predicts favorable clinical outcome in colon cancer. Immunotherapy. 2020;12(1):25–35. https://doi.org/10.2217/imt-2019-0145.
Hu Y, Lu Y, Xing F, Hsu W. FGFR1/MAPK-directed brachyury activation drives PD-L1-mediated immune evasion to promote lung cancer progression. Cancer Lett. 2022;547:215867. https://doi.org/10.1016/j.canlet.2022.215867.
This work was supported by the National Natural Science Foundation of China (Grant No. 81960525 and 82160591), Science and the Technology Commission of Shanghai Municipality (Grant No. 19DZ1930900), The Top-level Clinical Discipline Project of Shanghai Pudong (PWYgf2021-07), Key Project of Clinical Research of Shanghai East Hospital, Tongji University (Grant No. DFLC2022012) and Key Specialty Construction Project of Shanghai Pudong New Area Health Commission (Grant No. PWZzk2022-02).
Ethical approval and consent to participate
The public database used in this study is available for public access and re-analysis without requiring ethical approval from the local ethics committees. Therefore, this study does not require ethical approval.
Consent for publication
The authors declare no conflicts of interest in this work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The expression distribution of KLF family members and regulators of NOS enzymes in tumor tissues and normal tissues of TCGA-LIHC and ICGC-LIRI. Fig. S2. The prognostic survival value of KLF family members and regulators of NOS enzymes from TCGA-LIHC and ICGC-LIRI. Fig. S3. Analysis of the relationship between EMT-markers and KLF2 and the expression distribution of EMT-markers in the GSE25097. Fig. S4. Forest plots presentation of univariate Cox regression of KLFTs. Table S1: Target genes of the KLF2 from the CHEA Transcription Factor Targets dataset in Harmonizome platform, namely KLFTs.
Target genes of the KLF2 from the CHEA Transcription Factor Targets dataset in Harmonizome platform, namely KLFTs.
The co-expression analysis of KLF2 and KLFTs.
The common CAFs-related marker genes by consulting relevant references.
Univariate COX regression analysis of CAFs in the C1 subgroup.
About this article
Cite this article
Chen, XQ., Ma, J., Xu, D. et al. Comprehensive analysis of KLF2 as a prognostic biomarker associated with fibrosis and immune infiltration in advanced hepatocellular carcinoma. BMC Bioinformatics 24, 270 (2023). https://doi.org/10.1186/s12859-023-05391-0