Identification of TRPM2 as a Marker Associated With Prognosis and Immune Infiltration in Kidney Renal Clear Cell Carcinoma

TRPM2 (transient receptor potential melastatin-2), a Ca2+ permeable, non-selective cation channel, is highly expressed in cancers and regulates tumor cell migration, invasion, and proliferation. However, no study has yet demonstrated the association of TRPM2 with the prognosis of cancer patients or tumor immune infiltration, and the possibility and the clinical basis of TRPM2 as a prognostic marker in cancers are yet unknown. In the current study, we first explored the correlation between the mRNA level of TRPM2 and the prognosis of patients with different cancers across public databases. Subsequently, the Tumor Immune Estimation Resource (TIMER) platform and the TISIDB website were used to assess the correlation between TRPM2 and tumor immune cell infiltration level. We found that 1) the level of TRPM2 was significantly elevated in most tumor tissues relative to normal tissues; 2) TRPM2 upregulation was significantly associated with adverse clinical characteristics and poor survival of kidney renal clear cell carcinoma (KIRC) patients; 3) the level of TRPM2 was positively related to immune cell infiltration. Moreover, TRPM2 was closely correlated to the gene markers of diverse immune cells; 4) a high TRPM2 expression predicted worse prognosis in KIRC based on different enriched immune cell cohorts; and 5) TRPM2 was mainly implemented in the T-cell activation process indicated by Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. In conclusion, TRPM2 can serve as a marker to predict the prognosis and immune infiltration in KIRC through the regulation of T-cell activation. The current data may provide additional information for further studies surrounding the function of TRPM2 in KIRC.


INTRODUCTION
Transient receptor potential (TRP) ion channels are a family of membrane proteins that play diverse physiological and pathological roles. Previous studies have indicated that TRP channels have a strong diagnostic potential for various cancer types, especially in head and neck, kidney, and lung cancers, in which clinically useful diagnostic markers are not available (Park et al., 2016). Aberrant regulation of TRP channels results in various diseases, including numerous types of cancer.
TRPM2 is a Ca 2+ -permeable, non-selective cation channel activated by ADP-ribose (ADPR), temperature, oxidative stress, and Ca 2+ (Sumoza-Toledo and Penner, 2011b). The activation of TRPM2 results in the transport of Ca 2+ across the plasma membrane and the release of cytosolic Ca 2+ from organelles of the endoplasmic reticulum store and lysosome (Lange et al., 2009;Sun et al., 2012;Li F. et al., 2016). TRPM2 is ubiquitously distributed in the body, especially highly expressed in many energy-demanding tissues, including the brain, heart, and vasculature. The physiological functions of TRPM2 include insulin secretion in the pancreas, warm sensitivity in neurons (Tan and McNaughton, 2016;Paricio-Montesinos et al., 2020), and the induction of dendritic cell maturation and chemotaxis (Sumoza-Toledo et al., 2011a). TRPM2 expressed in immune cells, such as macrophages and polymorphonuclear neutrophils (PMNs), is crucial for regulating the bactericidal activity of phagocytic cells and PMN migration in tissues (Mittal et al., 2017).
Recent studies have demonstrated that TRPM2 was highly expressed in melanoma, breast cancer, prostate cancer, tongue cancer, neuroblastoma, and kidney cancer (Hopkins et al., 2015;Park et al., 2016). Knockdown/inhibition of TRPM2 impaired mitochondrial function and autophagy, reduced cellular bioenergetics, and increased the levels of reactive oxygen species (ROS), resulting in decreased tumor proliferation and/ or viability in many malignancies, suggesting a role of TRPM2 in cancer cell propagation and growth (Hopkins et al., 2015;Miller, 2019). Next, we questioned whether TRPM2 could be used as a prognostic marker in cancers. We screened the mRNA level of TRPM2 in the Gene Expression Profiling and Interactive Analysis (GEPIA) and TISIDB databases and found that the expression of TRPM2 was significantly enhanced in several cancers; of these, kidney renal clear cell carcinoma (KIRC) was the most TRPM2relevant cancer with respect to tumor subtypes, cancer stages, nodal metastasis, and tumor grades.
KIRC is the most common histological subtype that accounts for almost 90% of all kidney tumors Arora and Limaiem, 2021). Currently, surgical resection is the first-line treatment for KIRC. However, local recurrence or distant metastasis persists in 25% of patients with "local" disease following curative nephrectomy (Wang et al., 2020). Therefore, KIRC is one of the most aggressive kidney cancers mainly due to the high risk of tumor recurrence, metastasis, chemotherapy, and radiotherapy resistance (Linehan, 2012). Although multidisciplinary synthetic therapy has been used to treat KIRC, the prognosis and clinical outcomes have been unsatisfactory. Moreover, there are no credible predictive markers for the prognosis and treatment of individual sensitivity or resistance in KIRC, although some prognostic factors associated with the survival of KIRC patients have been described (Park et al., 2016;Marquardt et al., 2021). To the best of our knowledge, there have been no studies in the literature that investigated the role of TRPM2 in kidney cancers. Only one article has mentioned that the overexpression of TRPM2 might be used as a diagnostic marker for kidney cancer in terms of sensitivity and specificity, indicating strong diagnostic potential of TRPM2 (Park et al., 2016). In addition, a high expression of TRPM2 was closely associated with poor prognosis in bladder, head and neck, liver, and lung cancers (adenocarcinoma) , implying that TRPM2 significantly affects cancer progression (Park et al., 2016).
Features of the tumor immune microenvironment (TIME) are highly heterogeneous and have emerged as critical biomarkers in predicting the efficacy of and the response to systemic therapy. KIRC is one of the most immune-infiltrated tumors in pancancer comparisons (Vuong et al., 2019). Recent studies have considered KIRC as an immunogenic tumor with immune dysfunction partly as a result of the infiltration of immuneinhibitory cells, such as regulatory T cells (Tregs) and myeloid-derived suppressor cells, into the tumor microenvironment (TME) (Diaz-Montero et al., 2020). The association of TRPM2 with immune landscapes in KIRC may have clinical implications in defining distinct classes and subclasses of TIME.
In the present study, we described the performance of TRPM2 as a predictive marker of KIRC on its association with immune infiltration and then further explored the implicated mechanism based on relevant clinical background. Firstly, we used the Tumor Immune Estimation Resource (TIMER) and Gene Expression Profiling and Interactive Analysis (GEPIA) databases to systematically analyze the mRNA level of TRPM2 in different types of tumors and assessed its prognostic value using the Kaplan-Meier plotter and the GEPIA database. Furthermore, we used the TIMER and TISIDB databases to assess the correlation between TRPM2 and the abundance of tumorinfiltrating immune cells in the TME. Finally, TRPM2correlated genes were enriched using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), indicating that TRPM2 may be largely involved in T-cell activation in KIRC.

TIMER Database Analysis
TIMER (http://timer.comp-genomics.org/) is a comprehensive resource for systematic analysis of immune infiltrates across diverse cancer types. This platform provides four modules for investigating the association between the immune infiltrates and genetic or clinical features and four modules for exploring cancerrelated associations in The Cancer Genome Atlas (TCGA) cohorts (Li B. et al., 2016;Li et al., 2017;Li et al., 2020).
We used the Gene_DE module to investigate the mRNA expression of TRPM2 in different cancer types and matched normal tissues across all TCGA cancer types. The level of TRPM2 was expressed as log2 TPM, and the significance was determined by differential gene expression analysis using the edgeR software package.
We also used this platform to evaluate the correlation of the expression level of TRPM2 with tumor purity and tumor infiltration of B cells, CD4 + T cells, CD8 + T cells, neutrophils, macrophages, and dendritic cells. A heatmap of Spearman's correlations between the level of TRPM2 and immune cell infiltration across diverse cancer types was also generated using the TIMER database. Furthermore, the correlation between TRPM2 expression and the gene markers of immune cells identified the potential subtypes of infiltrating immune cells.
GEPIA GEPIA (http://gepia.cancer-pku.cn/) is an interactive web server for analyzing the RNA sequencing expression data of 9,736 tumors and 8,587 normal samples from TCGA and the Genotype-Tissue Expression (GTEx) project, respectively (Tang et al., 2017). The GEPIA database was utilized to evaluate the tumor/normal differential expression levels of TRPM2 and to conduct patient survival analysis based on the expression levels of TRPM2 across various cancer types. The GEPIA database utilizes log2(TPM + 1) for log scale with a | Log2FC| cutoff of 1 and a p-value cutoff of 0.01 in determining the differential expression of TRPM2, shown in Figure 1B.

Kaplan-Meier Plotter Database Analysis
The Kaplan-Meier plotter (http://kmplot.com/analysis/) can assess the effect of 54 k genes (mRNA, miRNA, and protein) on survival in 21 cancer types and discover and validate survival biomarkers (Gyorffy, 2021). The databases Gene Expression Omnibus (GEO), European Genome-Phenome Archive (EGA), and TCGA were used to assess the correlation between clinical outcomes and TRPM2 expression in different cancers and between TRPM2 expression and immune cell infiltration. We also calculated the hazard ratios (HRs) of 95% confidence intervals (CIs) and the log-rank p-value.

TISIDB
The TISIDB database (http://cis.hku.hk/TISIDB/index.php) is a portal for analyzing tumor and immune cell interactions that integrates multiple heterogeneous data types (Ru et al., 2019). Next, we analyzed the correlation between TRPM2 expression and clinical features, such as overall survival (OS), cancer stage, and cancer grade. Heatmap and dot plots showing the Spearman's correlations between the expression of TRPM2 and the abundance of tumor immune-infiltrating cells across diverse cancer types were also generated using the TISIDB database.

UALCAN Database Analysis
The UALCAN database (http://ualcan.path.uab.edu/index.html) is available for online analysis of cancer omics data (TCGA, MET500, and CPTAC). This database allows users to identify biomarkers, perform pan-cancer gene expression analysis, obtain patient survival information, and analyze epigenetic regulation of gene expression (Chandrashekar et al., 2017). We used this database to validate the results of the correlation between TRPM2 expression and clinical features in the TISIDB database. A p < 0.05 indicated statistically significant differences.

Statistical Analysis
The expression of TRPM2 was analyzed via the TIMER and GEPIA databases. The correlation between TRPM2 expression and survival prognosis, including OS, disease-free survival (DFS), and relapse-free survival (RFS) in pan-cancer, were analyzed in the GEPIA database and the Kaplan-Meier plotter platform. To compare the survival curves, we used the log-rank test to calculate the HRs, 95% CIs, and p-values in the Kaplan-Meier plotter and GEPIA. We used Spearman's correlation analysis to evaluate the correlation between gene expression and immune cell infiltration in the TIMER and TISIDB databases. A p < 0.05 was considered statistically significant.

Prognostic Significance of TRPM2 Expression in Various Cancers
Next, we investigated the prognostic value of TRPM2 in different cancers using GEPIA and the Kaplan-Meier plotter platform. In GEPIA, the Kaplan-Meier plots revealed that a high TRPM2 level was associated with a short OS and/or DFS in cancers, including KIRC ( Figure S1). These findings indicated that TRPM2 was a prognostic factor for short OS in KIRC. In contrast, the expression level of TRPM2 was positively correlated with the prognosis of BLCA, READ, UCES, THCA, and STAD, as predicted by GEPIA and the Kaplan-Meier plotter (Supplementary Figure S2).

Clinical Characteristics of TRPM2 in Cancers
In order to examine the association between TRPM2 and clinical features, we assessed the multiple clinical prognostic values of TRPM2 in various cancers using the TIDISB database. Among the human cancers, KIRC was the most correlated cancer with regard to OS ( Figure 3A), stage ( Figure 3B), and tumor grade ( Figure 3C). The association plots derived from the 533 KIRC cases in TCGA showed that the upregulation of TRPM2 expression was significantly associated with shorter OS, poor pathological stage, and tumor grade in KIRC ( Figure 3D-F, respectively). Using the UALCAN database, we further validated the finding that the expression level of TRPM2 was negatively correlated with the prognosis of KIRC patients, in consideration of the tumor subtype ( Figure 3G), cancer stage ( Figure 3H), nodal metastasis ( Figure 3I), and tumor grade ( Figure 3J). The results demonstrated a significant correlation between high TRPM2 levels and these adverse clinicopathological parameters, which was consistent with those described above ( Figure 2A). Therefore, we mainly focused on the function of TRPM2 in KIRC to examine its prognostic performance, clinical significance, and mechanism. To elucidate whether TRPM2 is an independent risk factor for clinical outcomes of KIRC patients, univariate and multivariate Cox analyses were performed. In the univariate Cox analysis, the pathological stage, age, histologic grade, serum calcium, and TRPM2 expression were significantly correlated with OS (p < 0.001, p < 0.001, p < 0.001, p < 0.001, and p 0.005, respectively) ( Figure 4A). To exclude the confounder effect, the significant prognostic factors detected in the univariate analysis ( Figure 4A) were evaluated further in the multivariate analysis ( Figure 4B). Strikingly, multivariate Cox analysis revealed that the pathological stage, age, serum calcium, and TRPM2 expression were independent risk factors for OS (p 0.039, p 0.034, p 0.023, and p 0.001, respectively) ( Figure 4B) of KIRC patients. Since calcium influx plays a critical role in tumor occurrence and development through TRPM2 channels, we considered serum calcium for the above Cox analysis. Interestingly, both univariate and multivariate analyses verified that serum calcium was an independent risk factor for the OS of KIRC patients (HR 2.329, 95% CI 1.125-4.824, p 0.023) ( Figure 4B). Based on the clinicopathological features and TRPM2, we also constructed a nomogram to predict the 1-, 5-, and 10-year OS rates using the Cox regression algorithm (Supplementary Figure S3).

Frontiers in
Then, the TISIDB database was used to further explore the correlations between the level of TRPM2 and the 28 tumor immune-infiltrating cell subtypes. The results showed that the level of TRPM2 was associated with 25 immune cell subtypes in KIRC ( Figure 5B). Notably, myeloid-derived suppressor cells (MDSCs) (r 0.843, p < 2.2e−16), T follicular helper (Tfh) cells (r 0.74, p < 2.2e−16), and macrophages (r 0.689, p < 2.2e−16) displayed relative strong correlation with TRPM2 expression in KIRC ( Figure 5C) Figure 5B). In addition, TRPM2 showed a weak correlation with immune cells of activated CD8 + cells (r 0.198, p 5.68e−05), CD56dim (r 0.161, p 0.00112), and CD56bright (r 0.133, p 0.00737) in BLCA ( Figure 5D). In addition, the correlation between the level of TRPM2 and the degree of immune infiltration in various cancers were analyzed using the TIMER database and displayed as a heatmap (Supplementary Figure S4). These findings strongly indicated that TRPM2 serves as a major tumor immune infiltration regulator in KIRC.

TRPM2 Expression Was Correlated With Immune Cell Type Markers
We assessed the correlation between the expression of TRPM2 and the level of tumor-infiltrating immune cell gene markers in KIRC and BLCA using the TIMER database. The results showed that the level of TRPM2 in KIRC tissues was strongly associated with immune markers of B cells, CD8 + T cells, dendritic cells, M1/ M2 macrophages, monocytes, neutrophils, general T cells, T-cell  Purity and none are the Spearman's rho between TRPM2 and immune infiltration level after or without purity adjustment, respectively. TAMs, tumor-associated macrophages; cor, correlation. *p < 0.01, **p < 0.001, ***p < 0.0001. exhaustion, tumor-associated macrophages (TAMs), Th1, Th2, and Tregs, but not with NK, Tfh, and Th17 cells (Table 1). A major factor determining tumor progression over time is the overall proportion and property of T cells within the TIME (Binnewies et al., 2018). Notably, the level of TRPM2 was significantly correlated with various subtypes of T cells based on cell markers, including CD8 + T-cell markers (CD8A and CD8B), T-cell (general) markers (CD3D, CD3E, and CD2), exhausted T-cell markers (CTLA4, GZMB, LAG-3, and PDCD 1), Th1 markers (TBX21, STAT4, STAT1, and IFNG), Th2 markers (GATA3 and STAT5A), Treg markers (FOXP3 and CCR8), and with neutrophil markers (ITGAM and CCR7), dendritic cell markers (ITGAX, HLA-DPA1, HLA-DRA, and HLA-DPB1), and B-cell markers (CD79A and CD19) in KIRC (Table 1). Also, a significant correlation was established between the level of TRPM2 and the expressions of marker genes in different subsets of macrophages, including M1 macrophage markers (IRF5), M2 macrophage markers (MS4A4A, VSIG4, and CD163), and TAMs (IL10 and CD68), in KIRC (Table 1). However, only nine immune cell marker genes were significantly associated with the levels of TRPM2 in BLCA. Furthermore, the expression of TRPM2 was not markedly related to most marker genes of CD8 + T cells, NK cells, Th2, and Th17 cells in BLCA. These findings revealed that TRPM2 is involved in the regulation of tumor immune infiltration in KIRC.

Role of TRPM2 Level in Prognosis Prediction of KIRC Based on Different Enriched Immune Cell Cohorts
TRPM2 is involved in the regulation of tumor immune infiltration in KIRC ( Figure 5 and Table 1). Moreover, upregulated TRPM2 indicated poor prognosis in KIRC/ PAAD patients, but better prognosis in BLCA patients. Thus, we hypothesized that TRPM2 may affect the prognosis of patients partially through the regulation of immune infiltration. Analysis of the Kaplan-Meier plotter revealed that the level of TRPM2 was associated with OS in KIRC, PAAD, and BLCA patients based on the enrichment of different immune cells. In KIRC, a high TRPM2 level in enriched CD4 + memory T cells, CD8 + T cells, macrophages, and Tregs indicated poor prognosis of patients with KIRC ( Figure 6A). However, no significant correlation was established between the expression of TRPM2 and OS in enriched B cells, mesenchymal stem cells, Th1 cells, and Th2 cells ( Figure 6A). On the other hand, in PAAD, a high expression of TRPM2 was related to poor prognosis only in CD8 + T cells and enriched B cells, but showed no statistically significant difference ( Figure 6B). The clinical samples for Th1 cell enrichment were too few to analyze in PAAD. In BLCA, low-TRPM2 expression groups showed short OS in enriched CD4 + memory T cells, CD8 + T cells, Tregs, mesenchymal stem cells, NK T cells, Th2 cells, Th1 cells, B cells, and macrophages ( Figure 6C). The above results suggested that TRPM2 may affect the prognosis of patients by potentially regulating the infiltration of immune cells in KIRC.

Enrichment of TRPM2-Correlated Genes in KIRC Patients by Gene Ontology/Kyoto Encyclopedia of Genes and Genomes Analysis
To gain in-depth insights into the molecular mechanism of TRPM2 as a prognostic marker in KIRC, the gene expression correlation was calculated with Pearson's correlation analysis to reveal the genes associated with TRPM2 in KIRC patients. Subsequently, we conducted GO and KEGG analysis of the TRPM2-related genes in KIRC. Pearson's correlation coefficients between TRPM2 and the genes of interest were calculated, and only genes with |r| > 0.5 and p < 0.05 were subjected to GO classification and KEGG pathway enrichment. The top 5 GO categories are shown in Figure 7. In KIRC, the above genes were enriched in several biological processes, such as T-cell activation, regulation of T-cell activation, lymphocyte activation, leukocyte cell-cell adhesion, and lymphocyte differentiation. For molecular function (MF), these genes were mainly associated with cytokine receptor activity, MHC protein binding, MHC protein complex binding, IgG binding, and cytokine binding. Figure 7 also shows the most significant KEGG pathways, namely, hematopoietic cell lineage, Staphylococcus aureus infection, leishmaniasis, osteoclast differentiation, and cell adhesion molecules. In addition, GO cell component (CC) analysis indicated that TRPM2 plays a key role in the regulation of T-cell activation on the external side of the plasma membrane, secretory granule membrane, tertiary granule membrane, plasma membrane receptor complex, and tertiary granule membrane in KIRC. Conversely, in BLCA, TRPM2-correlated genes (|r| > 0.5 and p < 0.05) failed to enrich in the GO and KEGG analyses. We expanded the range of TRPM2-correlated genes (|r| > 0.3 and p < 0.05) and found that only the CC was enriched, indicating that TRPM2 may function mainly in the preribosome, lamellipodium, large subunit precursor of preribosome, and tertiary granule membrane (Supplementary Figure S5).

DISCUSSION
TRPM2, a highly Ca 2+ -permeable cation channel of the TRPM family, regulates cancer cell growth and survival. Although, in a small number of malignancies, the activation rather than the inhibition of TRPM2 has been reported to reduce tumor cell survival (Di et al., 2012), most of the studies considered that inhibition of the expression or function of TRPM2 resulted in decreased tumor proliferation and/or viability in several malignancies (Klumpp et al., 2016;Belrose and Jackson, 2018). In this study, for the first time, we analyzed the correlation between TRPM2 expression and its prognostic value in various cancers and the tumor immune cell infiltration in KIRC.
The expression of TRPM2 was elevated in various cancers compared to that in normal tissues (Figure 1). Survival analysis showed that a high TRPM2 expression was associated with poor prognosis in KIRC, LGG, OV, THYM, UVM, and LIHC patients ( Figure 2) and better prognosis in BLCA, READ, UCEC, THCA, and STAD patients (Supplementary Figure S2). Among these cancer types, the expression level of TRPM2 was most relevant with the tumor subtype, stage, grade, and lymph node metastasis in KIRC (Figure 3). Moreover, TRPM2 was associated with TNM stage, especially with metastasis to lymph nodes and distant metastasis (Supplementary Table S1). These results suggested that TRPM2 plays a critical role in the progression and metastasis of KIRC, which was consistent with previous research studies showing that TRPM2 affects tumor growth and invasion and is correlated with poor prognosis in patients with breast, gastric, pancreatic, prostate, and head and neck cancers, melanoma, and neuroblastoma (Belrose and Jackson, 2018). Li et al. showed that TRPM2 mediated cancer cell migration through Ca 2+ and Zn 2+ (Li F. et al., 2016). In addition, TRPM2 downregulation inhibited the tumor growth of lung cancer in a xenograft NOD/SCID mouse model (Almasi et al., 2019). Furthermore, the strong association of TRPM2 with the clinical characteristics of KIRC and the results of the multivariate analysis indicated that TRPM2 is an independent risk factor for OS in KIRC patients (Figures 3 and 4). According to research findings, when TRPM2-mediated calcium influx is inhibited, the mitochondria are dysfunctional, cellular bioenergetics is reduced, and the production of ROS is increased, thereby decreasing tumor growth and increasing chemotherapy sensitivity (Miller, 2019). Therefore, serum calcium was taken into account for evaluation. Intriguingly, we found that serum calcium was also an independent risk factor in KIRC patients (Figure 4). These data implicated aberrant intracellular Ca 2+ signaling in the progression of KIRC.
Previously, TRPM2 was shown to regulate TME. It is expressed in monocytes and regulates CXCL2 production and neutrophil migration (Ambale-Venkatesh et al., 2017). Herein, we reported that the expression of TRPM2 was associated with several immune-infiltrating cells (B cells, CD8 + T cells, CD4 + T cells, macrophages, neutrophils, and dendritic cells) in KIRC through the correlation analysis in the TIMER database ( Figure 5 and Table 1). Notably, an increased TRPM2 level was positively associated with Tregs and Th cells ( Table 1). FOXP3 is a target for identifying Tregs in the TME, contributes to Treg differentiation, and mediates tumor immune escape (Sadlon et al., 2010). Consistently, our data revealed a positive association between the level of TRPM2 and FOXP3 (r 0.516, after cell purity correction). Moreover, the positive correlation between TRPM2 and several Th cell (Th1, Th2, Tfh, and Th17) markers implied a significant regulatory role of TRPM2 in T-cell function in KIRC. These observations suggested that TRPM2 affects the prognosis of KIRC patients by regulating the recruitment of various T-cell subtypes.
Further investigation of the correlation of TRPM2 with T-cell inhibitory receptors indicated that TRPM2 was also positively correlated with T-cell exhaustion signature markers after cell purity correction (PDCD1: r 0.583; CTLA4: r 0.525; HAVCR2: r 0.101; GZMB: r 0.277; LAG3: r 0.545) ( Table 1). Collectively, these data indicated that TRPM2 might be the key factor that modulates T-cell exhaustion and inhibition of antitumor immune responses. This finding was consistent with that of a previous study showing that TRPM2 critically influences T-cell proliferation and pro-inflammatory cytokine secretion following polyclonal T-cell receptor stimulation (Melzer et al., 2012).
In order to investigate the molecular mechanism underlying TRPM2 regulating immune infiltration and tumor progression, we analyzed the network and focused on the specific functions and pathways of TRPM2-related genes in KIRC using GO and KEGG analyses. Based on TCGA-KIRC data, the TRPM2associated genes with |r| > 0.5 and p < 0.05 were selected for enrichment. Interestingly, T-cell activation was ranked on top among biological processes (Figure 7). These results suggested that TRPM2 influences the progression of KIRC via T-cell activation.
In this study, we highlighted the role of TRPM2 as a potential prognostic marker by regulating immune cell infiltration in KIRC. Together with the results of the GO and KEGG enrichment analyses, the putative mechanism involved, but not limited to the regulation of T-cell activation and exhaustion, could provide clinical relevance of the function of TRPM2 in KIRC. However, the mechanism and clinical application of TRPM2 need to be substantiated further using in vitro and in vivo studies.

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. FIGURE 7 | Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of the TRPM2-related genes in kidney renal clear cell carcinoma (KIRC). Gene expression was calculated using Pearson's correlation analysis based on KIRC RNA-seq data from The Cancer Genome Atlas (TCGA), and then the TRPM2-related genes (|r| > 0.5 and p < 0.05) in KIRC were subjected to GO classification and KEGG pathway enrichment analysis (https://www.xiantao.love/).