KCNN4 is a potential prognostic marker and critical factor affecting the immune status of the tumor microenvironment in kidney renal clear cell carcinoma.

BACKGROUND: The tumor microenvironment (TME) has emerged as a crucial factor in cancer development and progression. Recent findings have indicated that tumor-infiltrating immune cells (TICs) in the TME may predict cancer prognosis and response to treatment. Herein, we sought to identify critical modulators of the kidney renal clear cell carcinoma (KIRC) TME.
METHODS: KIRC datasets from The Cancer Genome Atlas (TCGA) were analyzed using the ESTIMATE algorithm to determine the ImmuneScore and StromalScore. By profiling the differentially expressed genes (DEGs) in the ImmuneScore and StromalScore, we finally identified the immune- and stromal-related DEGs of the cases, through which we then performed intersection analysis to determine the immune-related genes (IRGs). Cox regression analysis and least absolute shrinkage and selection operator (LASSO) regression analysis were used to identify critical IRGs and construct a prognostic model. The CIBERSORT algorithm was used to calculate the relative content of 22 immune cell types. Finally, the datasets from the Gene Expression Omnibus (GEO) database were analyzed to validate results from the above analyses. Experimental validation was used on KIRC tissues by quantitative polymerase chain reaction (qPCR) and western blot.
RESULTS: We found that the ImmuneScore was negatively correlated with patients' prognosis. Intersection analysis of the ImmuneScore and StromalScore identified 118 IRGs that were enriched in immune-related functions. Following IRGs screening by Cox and LASSO regression analyses, six genes were identified and used to construct a KIRC prognostic model. Intersection analysis of these six genes and protein-protein interaction (PPI) were performed and obtained the most critical gene: Potassium Calcium-Activated Channel Subfamily N Member 4 (KCNN4). Further analysis showed that KCNN4 expression was higher in tumor samples relative to normal controls, and was negatively correlated with prognosis. CIBERSORT analysis revealed significant correlation between KCNN4 expression and multiple types of TICs, demonstrating that KCNN4 may affect KIRC prognosis by influencing the TME immune status. Ultimately, the GEO datasets and validation experiments confirmed that KCNN4 was highly expressed in tumor tissues compared to the corresponding normal tissues.
CONCLUSIONS: Our study demonstrated that KCNN4 might be a potential prognostic marker in KIRC, offering a novel therapeutic avenue. 2021 Translational Andrology and Urology. All rights reserved.

Introduction

Renal cell carcinoma (RCC) originates from renal epithelium and accounts for 2–3% of all adult malignancies. According to global cancer statistics, there were approximately 403,262 new RCC cases and 175,098 deaths in 2018 (1). Kidney renal clear cell carcinoma (KIRC) is the most common subtype of RCC, accounting for 80–90% of RCCs (2). In recent years, the incidence of KIRC has gradually climbed, while age at diagnosis continues to fall (3). KIRC prognosis is relatively poor and approximately 30% of patients have distant metastases at the time of diagnosis (4,5). Although multiple targeted therapies are used as first-line treatments for advanced and recurrent KIRC, improvements on disease free survival (DFS) and overall survival (OS) have been minimal due to variations in drug efficacy (6). Thus, it is of crucial to better characterize KIRC pathogenesis for improved clinical outcomes.

A growing body of studies has indicated that stromal cells are major components of tumors as they constitute the tumor microenvironment (TME) (7). Non-hematopoietic stromal cells in the TME include endothelial cells, mesenchymal stem cells (MSCs), fibroblasts, and pericytes (8). Moreover, various innate and adaptive immune cells, referred to as tumor-infiltrating immune cells (TICs), are also present in the TME (9). Stromal and tumor cells closely interact and influence each other, with the former secreting extracellular matrix (ECM) and inflammatory factors that affect tumor cell proliferation and metastasis (10), as well as the response and access to therapy (11,12). The mutual influence between tumor and stromal elements begins in the early stage of TME formation and changes dynamically over time (8). Intriguingly, TICs in the TME have been suggested as predictors of cancer prognosis and response to therapy. For instance, in colorectal cancer, increased M1 macrophage levels indicated poor prognosis (13), while low lymphocyte number may reflect poor prognosis in non-small cell lung cancer (14). Likewise, in high-grade serous ovarian cancer, neutrophils were associated with poor prognosis of patients (15), signifying that TME immune cell components were closely correlated with malignancy and may influence the response to treatment and clinical outcomes. Currently, high throughput sequencing has made it feasible to carry out genome-wide examination of biological processes and disease mechanisms, and this technique has been used to study the TMEs of various tumors (16,17).

In this article, we downloaded the data of 539 KIRC patients from The Cancer Genome Atlas (TCGA) database. The ESTIMATE and CIBERSORT algorithms were performed to calculate the ImmuneScore and StromalScore, as well as the relative proportions of various TICs in the respective cases, combined with the least absolute shrinkage and selection operator (LASSO) algorithm to screen immune-related genes (IRGs). Finally, the most critical gene, Potassium Calcium-Activated Channel Subfamily N Member 4 (KCNN4), was identified, which can predict the prognosis of KIRC patients and the immune status in the TME. The protein encoded by KCNN4 is part of a potentially heterotetrameric voltage-independent potassium channel that is activated by intracellular calcium and may be part of the predominant calcium-activated potassium channel in T-lymphocytes. KCNN4 was highly expressed in various cancers, including thyroid cancer (18), colorectal cancer (19), and hepatocellular carcinoma (20). Moreover, KCNN4 was strongly implicated in tumor invasion, metastasis, and epithelial-mesenchymal transition (EMT) (20,21). In this study, we found that KCNN4 may affect stromal and immune components of the TME in KIRC. These findings were validated using independent datasets from the Gene Expression Omnibus (GEO) database, and experimentally using quantitative polymerase chain reaction (qPCR) and western blot analyses of four KIRC patient samples collected at the First Affiliated Hospital of Guangxi Medical University. We present the following article in accordance with the REMARK reporting checklist (available at http://dx.doi.org/10.21037/tau-21-332).

Methods

Data extraction

Transcriptomic data on the KIRC cohort (539 tumor samples; 72 normal samples) and associated clinical data were obtained from TCGA (https://portal.gdc.cancer.gov/).

Analysis of ImmuneScore, StromalScore, and ESTIMATEScore

The ESTIMATE algorithm was used to calculate the proportion of immune and stromal cells in the TME using R (https://www.r-project.org/). Kaplan-Meier survival analysis was used to assess prognosis. Log-rank P<0.05 was considered significant.

Screening for IRGs

Tumor samples were classified into high- and low-score group based on the median ImmuneScore and StromalScore. Differentially expressed genes (DEGs) were identified with limma package using a cut-off threshold of |log 2 FC| >1 and a false discovery rate (FDR) <0.05. Following intersection analyses to identify IRGs, univariate Cox regression analysis was used to identify genes that affect the OS of patients. LASSO and multivariate Cox regression analyses were used to ascertain the critical IRGs.

Enrichment function analysis

Enrichment functional analyses of the IRGs were performed on R using the packages clusterProfiler, enrichplot, and ggplot2. P<0.05 and q<0.05 were considered statistically significant.

Construction of heatmaps for DEGs and IRGs

The package pheatmap was used for heatmap visualization of DEGs and IRGs.

Correlation between clinical features and scores

After excluding cases with incomplete clinical data, 530 KIRC cases remained, and the correlation between clinical features and three kinds of scores were analyzed on R. Wilcoxon rank sum or Kruskal-Wallis rank sum tests were used for significance tests.

Construction of protein-protein interaction (PPI) networks

PPI networks were constructed using the Search Tool for the Retrieval of Interacting Genes (STRING, https://string-db.org/) online database and reconstructed using Cytoscape version 3.8.1. Nodes with interaction confidence >0.4 were selected in building the network.

TIC abundance profile

TIC abundance in KIRC samples was determined using the CIBERSORT algorithm. 403 cancer samples (P<0.05) were applied to the following analysis.

Verification of results via the GEO database

The GSE66270, GSE57357, and GSE29609 datasets were downloaded from the GEO database and used as validation KIRC datasets. The Wilcoxon rank sum test was used to compare gene expression differences between tumors and corresponding normal tissues in the GSE66270 and GSE57357 datasets. Kaplan-Meier analysis was used to assess KIRC prognosis in the GSE29609 dataset.

Specimen collection

KIRC specimens and corresponding adjacent non-cancer specimens from four patients were acquired from the First Affiliated Hospital of Guangxi Medical University, and the tumor pathologic types were verified by immunohistochemistry. The patients had not received preoperative anti-cancer therapy. Written informed consent was obtained from the patients. All procedures performed in this study involving human participants were in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved by the Ethics Committee of the First Affiliated Hospital of Guangxi Medical University.

qPCR analysis

Total RNA was isolated from tissues using Trizol reagent (Thermo Fisher Scientific, MA, USA). Reverse transcription was carried out on 2 µg RNA using a cDNA synthesis kit (Thermo Fisher Scientific, MA, USA). qPCR was performed using SYBR Green Master (Roche, Switzerland) on a Mx-3000P quantitative PCR system (Stratagene, USA). KCNN4 primer sequences were 5'-ctgggtgctgtccgtgg-3' (forward) and 5'-agccgatggtcaggaatgtg-3' (reverse). β-actin was used as a housekeeping reference gene. Relative mRNA levels were assessed using the 2-ΔΔCt method.

Western blotting

Total protein was extracted from tissues using the RIPA lysis buffer (Solarbio, China) supplemented with PMSF (Solarbio, China). Equal amounts of denatured protein were separated by 10% SDS-PAGE and transferred onto PVDF membranes (Millipore, MA, USA). Membranes were then incubated overnight with the following primary antibodies: anti-β-actin [1:1,000, Cell Signaling Technology (CST), MA, USA]; and KCNN4 (1:1,000, Affinity Biosciences, OH, USA). They were then incubated with the secondary antibody (CST, MA, USA) at room temperature. The signal was developed by enhanced chemiluminescence (Millipore, MA, USA) and analyzed on ImageJ.

Statistical analysis

Statistical analysis was performed using SPSS 26.0 (IBM Corporation, USA) and GraphPad Prism 8 (GraphPad Software, USA). Student’s t-test was used to analyze the expression levels of KCNN4 in the patient samples of tumor tissue and control tissue. A P value <0.05 was considered to be statistically significant.

Results

Correlation between three kinds of scores and prognosis of KIRC

To assess the relationship between immune and stromal cell proportions and the patients’ prognosis, the ESTIMATE algorithm, combined with Kaplan-Meier analysis, were used to calculate the ImmuneScore, StromalScore, and ESTIMATEScore of the KIRC cases. Here, a higher ImmuneScore manifests a higher proportion of immune cells in the TME, while higher StromalScore reflects a higher proportion of stromal cells. The ESTIMATEScore is the sum of ImmuneScore and StromalScore. This analysis showed that the ImmuneScore was negatively correlated with KIRC OS (, P=0.033). However, the correlation between the StromalScore and ESTIMATEScore and prognosis did not differ significantly ().

Figure 1

Association between scores and survival time of KIRC patients. (A) Kaplan-Meier survival analysis for high- and low-score KIRC patients, grouped based on the median ImmuneScore. (B) Kaplan-Meier survival analysis of the StromalScore. (C) Kaplan-Meier survival analysis of the ESTIMATEScore. KIRC, kidney renal clear cell carcinoma.

Correlation of the scores with clinical features in the KIRC patients

Next, we analyzed correlation between the proportion of immune and stromal cells and patients’ clinical features like gender, tumor stage, and grade. This analysis revealed that females had a lower ImmuneScore relative to males (, P=0.032), and the ImmuneScore was positively correlated with tumor grade and stage (). Further analysis demonstrated that ImmuneScore increased significantly from T1 to T2 and T1 to T3 in T classification of tumor stage (). M classification also showed a positive correlation (, P=0.011). Nonetheless, no statistically significant correlation between the StromalScore and most of the clinical features (). The ESTIMATEScore showed a gender bias similar to that of the Immune score (, P=0.038). Meanwhile, the ESTIMATEScore was positively correlated with tumor grade and tumor stage (), and ESTIMATScore increased significantly from T1 to T3 in T classification of tumor stage (, P=0.0065). No statistically significant correlation between the ESTIMATScore and M classification (). Taken together, these data suggested that the proportions of immune and stromal cells in the TME of KIRC were correlated with tumor progression and metastasis.

Figure 2

Association between scores and clinical features of KIRC patients. (A,B,C,D,E) Association between gender, clinical stage, tumor grade, T classification, and M classification of KIRC patients and the ImmuneScore. (F,G,H,I,J) Association between different clinical features of KIRC patients and the StromalScore. (K,L,M,N,O) Association between different clinical features of KIRC patients and the ESTIMATEScore. KIRC, kidney renal clear cell carcinoma.

IRGs shared by the ImmuneScore and StromalScore and their enrichment analysis

To outline gene expression differences between immune and stromal cells of the TME, the ImmuneScore and StromalScore were divided into high- and low-score groups based on the median. Values higher than the median were classified into the high-score group and vice versa, and the DEGs were visualized in heatmaps (). This analysis identified 1,267 DEGs from the ImmuneScore and 119 DEGs from the StromalScore, which included up-regulated and down-regulated genes. We then intersected the ImmuneScore and StromalScore DEGs and obtained 118 overlapping IRGs that may affect both the stromal and immune components of the TME (). Subsequently, gene ontology (GO) enrichment analysis of the biological roles of the IRGs revealed that they were enriched in immune-related activities such as immune response and antigen binding (), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found the IRGs to be enriched in the chemokine signaling pathway and Fc gamma R-mediated phagocytosis (). Taken together, the results above clarified that IRGs were mainly enriched in immune-related functions, and also exert influence on stromal components, highlighting their significance in the TME.

Figure 3

DEGs shared by the ImmuneScore and StromalScore and functional enrichment analysis. (A) The heat map of DEGs generated based on the ImmuneScore. (B) The heat map of DEGs generated based on the StromalScore. (C) Venn plot of all 118 IRGs shared by the ImmuneScore and StromalScore. (D,E) GO and KEGG analyses for the 118 IRGs. DEG, differentially expressed gene; IRG, immune-related gene; GO, gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes.

IRGs screening and prediction model construction

Next, we conducted a downstream analysis of the 118 overlapping IRGs. Univariate Cox regression analysis was used to screen IRGs that affected KIRC prognosis and identified 26 statistically significant IRGs (). Further analysis of these 26 IRGs using LASSO regression analysis uncovered 11 key IRGs (). To further screen the IRGs, multivariate Cox regression analysis of these 11 IRGs was carried out, thereby obtaining KCNN4, AC022079.1, LINC02611, TNFSF14, FREM1, and PLG (). Finally, we constructed a KIRC prognostic model based on these six IRGs, as well as a nomogram for determining the 1-, 3-, and 5-year KIRC survival rate ().

Figure 4

Univariate Cox regression analysis of 118 IRGs. The analysis was based on the 26 key genes with a threshold P value <0.05. IRG, immune-related gene.

Figure 5

Screening of key IRGs using LASSO regression and multivariate Cox regression. (A) LASSO coefficient profiles of the IRGs associated with the overall survival of KIRC. (B) 10-time cross-validation for tuning parameter selection in the LASSO model for overall survival. (C) Forest plots of six genes in prediction model obtained by multivariate Cox regression analysis. IRG, immune-related gene; LASSO, least absolute shrinkage and selection operator; KIRC, kidney renal clear cell carcinoma.

Figure 6

Prediction model construction. (A) A nomogram for predicting the 1-, 3- and 5-year OS rate of KIRC patients based on six key genes. (B) Calibration curves for predicting patient survival rates at 1-, 3- and 5-years. (C) Heat map of the six genes in the prediction model. (D) The distribution of patients’ risk scores. (E) Kaplan-Meier analysis of KIRC patients with high and low risk scores. (F) Time-dependent ROC curves at 1-, 3- and 5-years. OS, overall survival; KIRC, kidney renal clear cell carcinoma; ROC, receiver operating characteristic.

According to median risk score (RS), all cases were divided into high- and low-RS groups, the expression level of each gene of cases was visualized by heat map (), survival status of cases were shown in . Kaplan-Meier analysis demonstrated that the prognosis of the high-RS group was significantly poorer than that of the low-RS group (). The time-dependent receiver operating characteristic (ROC) curves analysis confirmed that the model had strong clinical predictive ability (), with 1-, 3-, and 5-year area under the curves (AUCs) of 0.716, 0.686, and 0.713, respectively.

Intersection analysis of the PPI network with IRGs in the prediction model

To identify key genes affecting the TME status, we used the STRING database to construct a PPI network of the 118 overlapping IRGs, and displayed gene names and node numbers on bar plots of the top 42 genes (). Intersecting the 42 genes with the six IRGs in the prediction model, we determined only one gene from the above analysis, KCNN4, reflecting it as a potentially key factor affecting the TME ().

Figure 7

Construction of PPI networks. (A) Construction of interaction network with nodes of interaction confidence >0.4. (B) Top 42 genes sorted by the number of connected nodes. PPI, protein-protein interaction.

Figure 8

Venn plot showing the intersection of the top 42 genes sorted by number of connected nodes and the six genes in the prediction model. PPI, protein-protein interaction.

Comprehensive analysis of KCNN4 expression and its prognostic value in the KIRC

Further KCNN4 characterization showed that the expression of KCNN4 was remarkably higher in tumor tissues relative to normal controls (). Similar results were obtained from pairwise comparative analysis (). To validate the correlation between KCNN4 and KIRC prognosis, we divided all of the tumor samples into high- and low-KCNN4 groups based on the median KCNN4 expression. Survival analysis revealed that KIRC patients with high KCNN4 levels had a worse prognosis than those with low KCNN4 expression (). Further analysis of the relationship between KCNN4 levels and clinical features showed that both tumor stage and grade were positively correlated with KCNN4 expression (). Similar results were obtained based on T and M classifications of tumor stage separately (). These results demonstrated that KCNN4 may be a reliable indicator of KIRC prognosis.

Figure 9

Differential expression of KCNN4 in different tissues and its relationship with the prognosis and clinical features of KIRC patients. (A) Differential expression of KCNN4 between normal and tumor tissues. (B) Paired comparative analysis of KCNN4 expression in normal and tumor samples obtained from the same patient. (C) Survival analysis for KIRC patients with high and low expression of KCNN4. (D-H) The relationship between KCNN4 expression and gender, tumor stage, tumor grade, T classification, and M classification of KIRC patients. KIRC, kidney renal clear cell carcinoma.

Relationship between TIC distribution and KIRC prognosis

To expound the TIC distribution pattern in the TME, we used the CIBERSORT algorithm to evaluate the relative levels of 22 immune cell subtypes in the TME from KIRC samples. The relative proportions of each TIC subtype within tumor samples and the correlations between them were shown in . To identify key TICs that may affect prognosis, we divided KIRC patients into high and low groups based on the median level of each TIC subtype, and found that relative abundance of regulatory T cells and follicular helper T cells were negatively correlated with patient survival, suggesting that they promote progression and metastasis of KIRC (). In contrast, the relative abundance of resting mast cells, resting dendritic cells and resting memory CD4 T cells showed a positive correlation with patient survival, indicating that they sustained anti-tumor immunity ().

Figure 10

Profiling of TICs in the TME. (A) Histogram showing the relative proportions of 22 TICs in tumor tissue. (B) Heat map showing the correlation between the 22 types of TICs. The number and color of the boxes indicates the correlation coefficient between two different immune cells. (C,D,E,F,G) Five key immune cells that influenced patients’ survival. TIC, tumor-infiltrating immune cell; TME, tumor microenvironment.

KCNN4 was a critical factor influencing the immune status of the TME

Next, we analyzed correlation between KCNN4 expression and various TIC subtypes, the differences in the proportion of 22 types of TICs in tumor tissues for the high- and low-KCNN4 expression groups were shown in violin plots (). Following analysis found that resting memory CD4 T cells, activated dendritic cells, M1 and M2 macrophages, resting mast cells, monocytes, and resting NK cells were negatively correlated with KCNN4 expression. In contrast, activated memory CD4 T cells, CD8 T cells, regulatory T cells, follicular helper T cells, memory B cells, and plasma cells were positively correlated with KCNN4 expression, suggesting that KCNN4 affecting KIRC prognosis by influencing the immune status of the TME ().

Figure 11

Correlation between the relative proportion of TICs and KCNN4 expression. (A) Violin plot showing the difference in the proportion of 22 immune cells in KIRC tumor tissues for the high and low KCNN4 expression groups. (B) Scatter plot showing the correlation between the proportions of 13 TICs and the expression of KCNN4. TIC, tumor-infiltrating immune cell; KIRC, kidney renal clear cell carcinoma.

GEO datasets confirmed KCNN4 expression in tumors and its impact on patient prognosis

To validate the reliability and authenticity of the results in TCGA, we analyzed the GSE57357 and GSE66270 datasets. These GEO datasets are comprised of sequencing data from 72 tumor tissues and 14 corresponding normal samples. This analysis revealed that KCNN4 expression was higher in tumor samples relative to controls (). Analysis of the GSE29609 dataset to determine the relationship between KCNN4 and patients’ survival revealed that high KCNN4 levels were correlated with poorer survival (). However, this relationship was not statistically significant (P=0.088), and may be due to the small sample size in this dataset (n=39).

Figure 12

Verification of results via GEO database. (A) KCNN4 was up-regulated in tumor tissues based on analysis of the GSE53757 dataset. (B) KCNN4 was differentially expressed between tumors and corresponding normal tissues in the GSE66270 dataset. (C) KCNN4 expression was associated with the survival of KIRC patients in the GSE29609 dataset. GEO, Gene Expression Omnibus.

Validation experiments verified the difference in KCNN4 expression

We examined KCNN4 expression in four KIRC samples and normal adjacent non-tumor controls by qPCR and western blot analyses, and found KCNN4 levels to be significantly higher in tumor tissues than in control tissues (P<0.05, ). The clinical features of the patients were shown in .

Figure 13

Expression of KCNN4 increased in KIRC samples. (A) Expression of KCNN4 in KIRC samples and normal adjacent non-tumor controls were assessed by qPCR (*P<0.05). (B,C) Representative blots and quantification results of western blot showed that expression of the KCNN4 protein in KIRC was significantly higher than controls (*P<0.05). KIRC, kidney renal clear cell carcinoma.

Table 1

The clinical features of the four KIRC patients

Patients	Gender	Age	Tumor size/cm	Location of tumor	T stage	N stage	M stage	Pathological type
Patient 1	Male	64	4.7×4.3×4.5	Left	T1b	N0	M0	KIRC
Patient 2	Female	59	7.9×6.2×7.0	Left	T2a	N0	M0	KIRC
Patient 3	Male	34	4.8×4.0×3.5	Left	T1b	N0	M0	KIRC
Patient 4	Female	69	4.7×3.5×4.3	Left	T1b	N0	M0	KIRC

KIRC, kidney renal clear cell carcinoma.

Discussion

In recent years, the role of the TME in tumorigenesis and cancer progression has become increasingly clear (22). The TME is infiltrated by a variety of TICs, which may have pro- or anti-tumor effects. Thus, identification of the factors affecting the immune status of the TME may uncover novel anti-tumor avenues. In this study, we analyzed KIRC patient datasets from TCGA using ESTIMATE and found that the proportion of immune components in the TME was correlated with the prognosis of KIRC. These findings have potential therapeutic implications.

Recent studies have confirmed that the anti-tumor properties of immune cells are compromised within the TME of KIRC, as tumor-secreted factors could influence dendritic cell differentiation and induce anergy-related genes in T cells, thereby inhibiting the anti-tumor reactivity in the tumor milieu (23). PD-1, an immunosuppressive receptor belonging to the CD28/CTLA-4 family and can exert inhibitory activity during the effective stage of T cell activation in the TME by regulating inhibitory signals (24,25). Hence, monoclonal antibodies against PD-1 or its ligands are being used to activate T cells. Better understanding of tumor immune escape mechanisms and the TME may expand immunotherapy strategies. KIRC immunotherapy has developed considerably, and many immune checkpoint inhibitors (ICIs) have entered clinical trials and shown encouraging results (26). Notwithstanding, there is still a proportion of patients who respond poorly to treatment. Additionally, side effects and adverse events have been identified in clinical trials, which may limit the use of such treatments (27). Through analyzing KIRC transcriptomic data from TCGA, we found that KCNN4 was positively correlated with clinical stage and tumor grade, illustrating its association with poor patient prognosis. Our data highlighted KCNN4 as a potential prognostic marker and therapeutic target in KIRC.

KCNN4 is a member of the calcium-dependent potassium channel family, which are crucial regulators of membrane potential, hormone secretion, epithelial function, cell proliferation, and apoptosis (28,29). Recent studies discovered that KCNN4 promoted the proliferation and invasion of lung cancer (30), endometrial cancer (31), and glioblastoma (32). Notably, high KCNN4 expression was correlated with poor prognosis in these cancers. Our data demonstrated that KCNN4 expression in KIRC was significantly higher than in normal tissues, and its expression level was positively correlated with clinical stage and tumor grade. These findings shed the light on the fact that KCNN4 may be associated with tumor occurrence, development, and metastasis, which is consistent with findings of Rabjerg et al. (33). To study the relationship between KCNN4 expression and the TME, we analyzed the correlation between KCNN4 and various TIC subtypes. The results showed that the expression of KCNN4 was correlated with a variety of immune cells, indicating that KCNN4 may affect the immune status of the TME and KIRC prognosis. Further evidence was derived from various studies, which suggested that KCNN4 altered the antigen presentation and functions of various immune cells (34,35). Chimote et al. reported that KCNN4 upregulation improved cancer immune monitoring and response to immunotherapy (35), which may be one of the mechanisms through which KCNN4 affects immune activity in the TME.

We subsequently validated our results using GEO datasets and clinical samples, confirming that KCNN4 mRNA and protein levels were higher in tumors than in normal tissues. Nevertheless, this study has limitations that should be noted. Firstly, we did not experimentally test whether high KCNN4 levels promote tumor cell invasion and metastasis. Secondly, we did not elucidate the mechanism by which KCNN4 affects the immune status of the TME.

Conclusions

In this study, using bioinformatics to evaluate TCGA datasets, we concluded that KCNN4 may influence KIRC prognosis by affecting the immune status of the TME. Our study highlighted KCNN4 as a potential prognostic factor and therapeutic target in KIRC.

The article’s supplementary files as