Comprehensive Analysis of the Role of <i>SLC2A3</i> on Prognosis and Immune Infiltration in Head and Neck Squamous Cell Carcinoma

Chu, Manru; Zheng, Ke; Li, Xiaojie; Luo, Zhiqiang; Yang, Xin; Wei, Changbo

doi:https://doi.org/10.1155/2022/2371057

Analytical Cellular Pathology

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 2371057 | https://doi.org/10.1155/2022/2371057

Comprehensive Analysis of the Role of SLC2A3 on Prognosis and Immune Infiltration in Head and Neck Squamous Cell Carcinoma

Manru Chu,¹Ke Zheng,²Xiaojie Li,¹Zhiqiang Luo,¹Xin Yang,³and Changbo Wei¹

Academic Editor: Elena Andreucci

Received30 May 2022

Revised01 Aug 2022

Accepted26 Aug 2022

Published19 Sept 2022

Abstract

Background. SLC2A3 is upregulated in various cancer types and promotes proliferation, invasion, and metabolism. However, its role in the prognosis and immune regulation of head and neck squamous cell carcinoma (HNSCC) is still obscure. This study is aimed at exploring the prognostic and immunotherapeutic potential of SLC2A3 in HNSCC. Methods. All data were downloaded from TCGA database and integrated via R software. SLC2A3 expression was evaluated using R software, TIMER, CPTAC, and HPA databases. The association between SLC2A3 expression and clinicopathologic characteristics was assessed by R software. The effect of SLC2A3 on survival was analyzed by R software and Kaplan-Meier Plotter. Genomic alterations in SLC2A3 were investigated using the cBioPortal database. Coexpression of SLC2A3 was studied using LinkedOmics and STRING, and enrichment analyses were performed with R software. The relationship between SLC2A3 expression and immune infiltration was determined using TIMER and TISIDB databases. Immune checkpoints and ESTIMATE score were analyzed via the SangerBox database. Results. SLC2A3 expression was upregulated in HNSCC tissues compared to normal tissues. It was significantly related to TNM stage, histological grade, and alcohol history. High SLC2A3 expression was associated with poor prognosis in HNSCC. Coexpression analysis indicated that SLC2A3 mostly participated in the HIF-1 signaling pathway and glycolysis. Furthermore, SLC2A3 expression strongly correlated with tumor-infiltrating lymphocytes in HNSCC. Conclusion. SLC2A3 could serve as a potential prognostic biomarker for tumor immune infiltration in HNSCC.

1. Introduction

Head and neck squamous cell carcinoma (HNSCC) is one of the top 10 cancer types worldwide and has a high mortality rate [1]. Despite advances in comprehensive cancer therapies, including traditional surgery, radiotherapy, and chemotherapy, the 5-year survival rate for HNSCC remains around 50%; this substantially increases the burden on family and society because most patients are middle-aged adults [2, 3]. The high mortality due to HNSCC accounts for the aggressive proliferation resulting in high recurrence and metastasis. The lack of early screening and diagnostic biomarkers is responsible for the low survival rate, while most patients with metastatic HNSCC die within 1 year after diagnosis [4].

Recently, many endeavors have been made to explore the biomarkers for more effective treatments for HNSCC [4–6]. Therapeutic advances in immunotherapy have revolutionized as promising methods for cancer treatment, and a small number of patients have benefited from them. The tumor immune microenvironment (TIME) is a dynamic structure comprising complex components, in which the immune system can influence the growth and mutation of tumor cells by yielding antitumor immune responses [7]. At the same time, tumor cells also damage immune cells differently. HNSCC is one of the most immunosuppressive human cancer types, so immune-based therapy could effectively prolong patients’ survival by preventing tumor immune escape and stimulating antitumor immune responses. Thus, it is imperative to explore novel noninvasive biomarkers for improving the prognosis and immunotherapy in HNSCC.

SLC2A3, also annotated as GLUT3, is a member of the solute carrier 2A (SLC2A) family participating in glucose transport across the hydrophobic cell membrane in mammalian cells. Cancer cells mainly rely on glycolysis for nutrient intake and energy production [8]. There is a growing body of literature that focusing on the upregulating of SLC2A3 to research cancer diagnosis and therapy [9]. In colorectal cancer, SLC2A3 plays a crucial role in perineural invasion via the regulation of epithelial-mesenchymal transition [10]. SLC2A3 could predict poor prognosis in gastric cancer, associated with macrophage infiltration in the tumor microenvironment (TME) [11]. Similarly, SLC2A3 overexpression in triple-negative breast cancer could activate inflammatory M1 tumor-associated macrophages resulting in the aggressiveness of tumor cells [12]. With the deep interrogation of glycolysis-related gene SLC2A3 in diverse cancer types, there is a pressing need to research the correlation between SLC2A3 and immune infiltration. However, the relationship between SLC2A3 and HNSCC remains unclear.

In the present study, we identified the differential expression levels of SLC2A3 gene and protein between tumor tissues and normal tissues based on The Cancer Genome Atlas (TCGA) and several online databases. The interaction networks of SLC2A3 and enrichment analyses were performed. Prognostic value and clinicopathologic features of SLC2A3 in HNSCC were evaluated. We conducted immune infiltration analyses to explore the potential immune mechanisms of SLC2A3 in HNSCC. Our study suggests that SLC2A3 could act as a valuable prognostic biomarker. The study also sheds light on the association of SLC2A3 with immune infiltration in HNSCC.

2. Materials and Methods

2.1. Data Extraction

The transcriptome profiling of SLC2A3 and related clinical information was extracted from TCGA database (https://portal.gdc.cancer.gov/). The RNA-seq gene expression with the format of FPKM was transformed into TPM format, and log2 conversion was conducted for further analyses. Gene expression data of SLC2A3 in HNSCC were acquired from TCGA and involved 502 tumor tissue and 44 adjacent normal tissue data. According to the median expression level of SLC2A3, gene expression data were divided into two subgroups: the high expression group and the low expression group. The samples contained several clinical characteristics, including TNM stage, clinical stage, histologic grade, gender, age, smoke condition, alcohol history, radiation therapy, primary therapy outcome, and lymphovascular invasion (Table 1). The RNA expression data were performed by using R software (version 3.6.3). All the data were downloaded from public databases; therefore, no ethical approval document was required.

2.2. Gene Expression of SLC2A3 in Human Cancers and HNSCC

The gene expression of SLC2A3 in human cancers/HNSCC and paired normal tissues were analyzed using the “ggplot2” package. Tumor Immune Estimation Resource (TIMER) is a user-friendly online web resource for the systematic analysis of immune infiltrates in different cancers (https://cistrome.Shinyapps.io/timer/). We used the TIMER database to verify the differential expression of SLC2A3 in pan-cancer [13].

2.3. Protein Expression of SLC2A3

UALCAN is a comprehensive web resource referring to publicly available cancer omics data (http://ualcan.path.uab.edu). Clinical Proteomic Tumor Analysis Consortium (CPTAC) is focused on molecular mechanisms of cancer via the application of proteome and genome analysis or proteogenomics (https://proteomics.cancer.gov/programs/cptac) [14]. The Human Protein Atlas (HPA) is aimed at exploring all the human proteins in cells, tissues, and organs, and all the data were freely available to the public for research on the human proteome (https://proteinatlas.org/) [15]. In the present study, UALCAN was used to analyze the protein expression of SLC2A3 from CPTAC, and HPA was performed to analyze the protein expression of SLC2A3 in normal tissue and tumor tissue of HNSCC.

2.4. Correlation between Clinical Pathological Characteristics and SLC2A3

To explore the association between SLC2A3 expression and clinical pathological characteristics including TNM stage, clinical stage, gender, age, histologic grade, smoker condition, and alcohol history, we applied the “ggplot2” package based on the TCGA database. The statistical significance was set at .

2.5. The Prognostic Value of SLC2A3

The ROC curve was drawn to examine the cutoff value of SLC2A3 using the “pROC” and “ggplot2” packages. The overall survival (OS) and disease-specific survival (DSS) between the high expression group and the low expression group were evaluated using the “survival” package and visualized using the “survminer” package. was regarded as statistically significant. The Kaplan-Meier Plotter collected discovery and validation of survival biomarkers to analyze the correlations between gene expression and survival in multiple tumor types (http://kmplot.com/analysis/) [16]. According to the Kaplan-Meier Plotter database, we explored the relapse-free survival (RFS) associated with SLC2A3 expression levels in HNSCC. The prognostic value was statistically significant when the .

2.6. Construction and Evaluation of Nomogram

To predict 1-, 3-, and 5-year survival probability associated with SLC2A3 in HNSCC, we applied the “RMS” package to construct a predictive model of nomogram including several clinical characteristics. Calibration and discrimination analysis is a very effective predictive model evaluation method. In this study, variables of radiation therapy, primary therapy outcome, lymphovascular invasion, and SLC2A3 expression were brought in to map the nomogram-predicted probabilities. The 45-degree line in calibration curves represented the best predictive value. The C-index and ROC curves were used to examine the predictive accuracies of the nomogram and individual prognostic factors. The statistical significance was set at .

2.7. Genomic Alterations of SLC2A3 in HNSCC

cBio Cancer Genomics Portal (cBioPortal) is an online website for analyzing multidimensional cancer genomics data (http://cbioportal.org). Thus, cBioPortal was used to investigate the genomic alterations in HNSCC [17].

2.8. Coexpression and Enrichment Analysis of SLC2A3

The LinkedOmics database provides a public platform to analyze cancer multiomics data from all 32 TCGA cancer types and 10 CTPCA cancer cohorts (http://www.linkedomics.org/login.php). Thus, we applied LinkedOmics to examine the coexpression of SLC2A3 with Pearson’s correlation coefficient [18]. Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) is a database for computational predicted physical and functional protein-protein interactions (http://string-db.org, version 11.5). In this study, STRING analysis was conducted to explore the coexpression genes associated with SLC2A3 [19]. The network was built, and the minimum required interaction score was set at 0.7 (high confidence). Functional enrichment analyses were performed using the “ggplot2” and “ClusterProfiler” packages. Briefly, Gene Oncology (GO) was annotated, and the Kyoto Encyclopedia of Genes and Genomes (KEGG) was analyzed to illustrate the potential mechanisms.

2.9. Correlation between SLC2A3 Expression and Immune Infiltration Cells

We applied TIMER to investigate the association between immune infiltrates (B cell, CD8⁺ T cell, CD4⁺ T cell, macrophage, neutrophil, and dendritic cell) and SLC2A3 expression or clinical outcome in HNSCC. Tumor Immune System Interaction Database (TISIDB) integrates diverse data types for exploring the interaction between tumors and the immune system (https://cis.hku.hk/TISIDB) [20]. In this study, we performed TISIDB to examine the correlation between SLC2A3 expression and the abundance of 28 tumor-infiltrating lymphocytes (TILs) acquired by gene set variation analysis. The association between immune or molecular subtypes and SLC2A3 expression in HNSCC was also explored via TISIDB. The difference at a was regarded as statistically significant.

2.10. Relationship between SLC2A3 Expression and Immune Checkpoint (ICP) Genes, and ESTIMATE Score of TME in HNSCC

SangerBox website (http://sangerbox.com/Tool) was used to investigate the correlation between SLC2A3 expression and immune checkpoint (ICP) genes, and ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data) score based on TCGA and GTEx databases. 60 ICP genes were analyzed in this study. ESTIMATE is an algorithm to predict levels of infiltrating stroma and immune cells by calculating stroma and immune scores, further deducing tumor purity [21]. Differences with a were regarded to be statistically significant.

2.11. Statistical Analyses

All statistical analyses were performed by R software (version 3.6.3). Mann–Whitney test and Wilcoxon signed rank test were used to compare SLC2A3 expression levels in normal tissues and tumor tissues. The Wilcoxon rank sum test was performed to evaluate the relationships between clinical-pathologic characteristics. The “pROC” package was used to examine the cutoff value of SLC2A3 in the ROC curve. The log-rank test and Cox regression analyses were used to evaluate the OS.

3. Results

3.1. The Differential Expression of SLC2A3 in Pan-Cancer

Analysis of the mRNA expression of SLC2A3 in 33 different types of cancers showed that SLC2A3 was upregulated in six cancer types, including cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), esophageal carcinoma (ESCA), head and neck carcinoma (HNSC), kidney renal clear cell carcinoma (KIRC), and stomach adenocarcinoma (STAD) compared to normal tissues (Figure 1(a)). However, SLC2A3 expression levels were significantly lower in nine types of cancer, including bladder urothelial carcinoma (BLCA), breast invasive carcinoma (BRCA), cervical and endocervical cancer (CESC), glioblastoma multiforme (GBM), kidney chromophobe (KICH), kidney renal papillary carcinoma (KIRP), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC) and uterine corpus endometrial carcinoma (UCEC). As displayed in Figure 1(b), the results of the TIMER database were similar to those obtained by R software. The SLC2A3 expression level was significantly higher in CHOL, COAD, ESCA, HNSC, KIRC, and STAD than in normal tissue. However, the SLC2A3 expression level was lower in BCLA, BRCA, KICH, KIRP, liver hepatocellular carcinoma (LIHC), LUAD, LUSC, and UCEC.

(a)

(b)

3.2. mRNA and Protein Overexpression of SLC2A3 in HNSCC

Unpaired data analyses showed that the mRNA expression levels of SLC2A3 in 502 HNSCC tissues were significantly higher than those in 44 normal tissues (; Figure 2(a)). In the paired samples, the results were consistent with the above analyses (Figure 2(b)). SLC2A3 was significantly upregulated in 43 tumor tissues than in paired normal tissues (). To explore the protein expression of SLC2A3, we performed analysis from HPA and CPTAC with UALCAN. The results revealed that the protein of SLC2A3 was overexpressed in HNSCC compared to normal tissues (Figures 2(c) and 2(d)).

(a)

(b)

(c)

(d)

3.3. Expression of SLC2A3 and Different Clinicopathologic Variables

The Wilcoxon rank sum test showed that the expression level of SLC2A3 was significantly related to T stage (), N stage (), M stage (), histological grade (), and alcohol history (). However, no significant correlations were detected between SLC2A3 expression and other clinicopathologic variables, such as clinical stage (), gender (), age (), and smoker condition () (Figure 3).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

3.4. High Expression of SLC2A3 Related to Poor Prognosis of HNSCC

Several methods were applied to analyze the prognosis associated with SLC2A3 expression in HNSCC. Firstly, the ROC curve showed promising predictive power with an AUG value of 0.727 (95% CI: 0.635-0.819; Figure 4(a)). Secondly, the results from “survival” and “survminer” packages showed that higher expression level of SLC2A3 was significantly related to poorer OS (, ; Figure 4(b)) and DSS (, ; Figure 4(c)). Meanwhile, analyses from the Kaplan-Meier Plotter database corresponded with the above results. Overexpression of SLC2A3 was correlated with RFS (, ; Figure 4(d)). Lastly, a nomogram was constructed to analyze the survival independently associated with several clinical characteristics, including radiation therapy, primary therapy outcome, lymphovascular invasion, and SLC2A3 expression (Figure 4(e)). The points scale was used to distribute every single variable, and the total points scale represented the sum of a point from every single variable. The survival probability at 1, 3, and 5 years corresponded to the total points. Thus, SLC2A3 was a poor prognostic factor for OS in HNSCC. The C-index of the nomogram was 0.714 (95% CI: 0.687-0.741), and the time-dependent ROC curve to 1-, 3-, and 5-year OS showed the AUG value of 0.548, 0.578, and 0.535, respectively (Figure 4(f)). In addition, the bias-corrected line in the calibration was close to the 45-degree line (Figure 4(g)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 4

Effect of SLC2A3 on HNSCC prognosis. ROC curve showed promising predictive power (a). High expression of SLC2A3 related to poor OS and DSS, respectively, analyzed using R software (b, c). High expression of SLC2A3 related to RFS, analyzed using Kaplan-Meier Plotter database (d). A nomogram for predicting 1-, 3-, and 5-year survival probability (e). ROC curve for 1-, 3-, and 5-year OS showed the AUG values of 0.548, 0.578, and 0.535, respectively (f). Calibration plot of the nomogram-predicted survival (g).

3.5. Genomic Alterations of SLC2A3 in Six Data Sets

The exploration of genomic alterations of SLC2A3 in six HNSCC data sets was conducted with the cBioPortal website. Genomic alteration of SLC2A3 occurred in 1.8% of patients with HNSCC (Figure 5(a)). The gene alteration type was diverse associated with SLC2A3 mRNA expression (Figure 5(b)). Copy number variation (CNV) mainly occurred in HNSCC. While narrowing the cancer field, no CNV was detected in oral squamous cell carcinoma (OSCC) (Figure 5(c)).

(a)

(b)

(c)

3.6. SLC2A3 Coexpression Networks and Enrichment Analyses

With the LinkedOmics database, coexpression genes of SLC2A3 were displayed as volcano plots in which 7089 genes were significantly positively correlated to SLC2A3, and 4057 genes were negatively related to SLC2A3 (Figure 6(a)). The heat maps showed the top 50 genes positively and negatively associated with SLC2A3, respectively (Figures 6(b) and 6(c)). Figure 6(d) shows a PPI network of SLC2A3 and its eight coexpression genes. We conducted “clusterProfiler” and “ggplot2” with these eight genes to explore the GO and KEGG enrichment. Coexpression genes were mainly involved in several key functions, including biological processes of response to oxygen levels, response to decreased oxygen levels, and response to hypoxia; cellular components of nuclear transcription factor complex, RNA polymerase II transcription factor complex, and parallel fiber to Purkinje cell synapse; and molecular functions of hexose transmembrane transporter activity, glucose transmembrane transporter activity, and glucose binding (Figure 6(a)). At the same time, KEGG results showed that coexpression genes enriched in central carbon metabolism in cancer, the HIF-1 signaling pathway, and glycolysis/gluconeogenesis. In addition, the correlation analyses between coexpression genes and SLC2A3 in HNSCC are displayed in Figure 6(f)–6(m).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

(j)

(k)

(l)

(m)

3.7. Correlation between SLC2A3 Expression and Immune Cell Infiltration

The correlation between SLC2A3 expression and six immune cells was analyzed using the TIMER database (Figure 7(a)). In HNSCC, SLC2A3 expression level was significantly related to B cell (, ), CD8⁺ T cell (, ), CD4⁺ T cell (, ), macrophage (, ), neutrophil (, ), and dendritic cell (, ). The cumulative survival was significantly associated with the abundance of B cells according to Figure 7(b). While the correlations between the expression levels of SLC2A3 and 28 TIL types across human cancers were explored in the TISIDB database (Figure 7(c)). Twenty-two types of TILs significantly related to SLC2A3 expression levels in HNSCC, such as central memory CD8 T (Tcm CD8) cells, type 1 helper (Th1) cells, regulatory T (Treg) cells, memory B (Mem B) cells, natural killer (NK) cells, natural killer T (NKT) cells, plasmacytoid dendritic cells (pDCs), macrophages, eosinophils, and mast cells (mast) (Figure 7(d)). Additionally, the role of SLC2A3 on immune and molecular subtypes in HNSCC was analyzed (Figures 7(e) and 7(f)). Immune subtypes were divided into six types: C1 (wound healing), C2 (IFN-gamma dominant), C3 (inflammatory), C4 (lymphocyte depleted), C5 (immunologically quiet), and C6 (TGF-b dominant). SLC2A3 showed high expression in the C6 type and low expression in C1 and C2 types. The molecular subtypes of cancers were classified into four types: atypical, basal, classical, and mesenchymal. SLC2A3 was upregulated in classical and mesenchymal types and downregulated in atypical and basal types. For different immune and molecular subtypes, SLC2A3 expression was significantly different in HNSCC.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 7

Correlation between SLC2A3 expression and immune cell infiltration. SLC2A3 expression was significantly related to the abundances of B cells, CD8⁺ T cells, CD4⁺ T cells, macrophages, neutrophils, and dendritic cells (a). Cumulative survival was poor significantly associated with the abundance of B cells (b). Correlations between expression of SLC2A3 and 28 TIL types across human cancers, investigated using the TISIDB database (c). SLC2A3 was significantly correlated with abundances of Tcm CD8 cells, Th1 cells, Treg cells, Mem B cells, NK cells, NKT cells, pDCs, macrophages, eosinophils, and mast cells (d). Relationship between SLC2A3 expression and immune subtypes (e). Relationship between SLC2A3 expression and molecular subtypes (f).

3.8. SLC2A3 Expression Is Related to Immune Checkpoint (ICP) Genes and ESTIMATE Score in HNSCC

We found that 53 of 60 ICP genes significantly correlated to SLC2A3 expression in HNSCC, in which 21 genes were inhibitors and 32 were stimulators (Figure 8(a)). The associations between SLC2A3 expression and three kinds of ESTIMATE score (Figure 8(b)–8(d)). The results indicated significant positive relationships with stromal score (), immune score (), and ESTIMATE score () in HNSCC.

(a)

(b)

(c)

(d)

4. Discussion

The SLC2A family comprises 14 members from SLC2A1-SLC2A14, among which those genes exhibit a tissue-specific expression. For instance, SCL2A1 is widely expressed in most tissues, while SLC2A2 is predominantly expressed in liver tissues [22]. SLC2A3 was mainly expressed in the central nervous system, which relies on steady support for ATP generation [23]. However, overexpression of SLC2A3 was found in several cancers such as colorectal cancer, gastric cancer and breast cancer [10–12]. In acute myeloid leukemia (AML), a SLC2A3-knockdown cell line accounted for the less effective demethylation and vitamin C in the AML and lymphoma cell lines suggesting SLC2A3 could act as a predictive biomarker in vitamin C treatment in AML [24]. In GBM, SLC2A3 played a vital oncogenic role in attenuating cell proliferation and glucose uptake via miR-106a [25]. Thus, According to the mediation in glucose transport, SLC2A3 was regarded as a cancer-promoting gene. The prognostic role of SLC2A3 in cancers has recently gained attention. High SLC2A3 expression resulted in shorter OS and disease-free survival in colorectal carcinoma (CRC) [26]. We performed bioinformatics to study the correlation between SLC2A3 and prognosis in HNSCC based on transcriptome profiling using the TCGA database and several online websites. It was found that SLC2A3 is differentially expressed in many cancer types, which was consistent with the above results. Moreover, upregulation of SLC2A3 in HNSCC was demonstrated, and high expression of SLC2A3 was significantly related to TNM stage, histologic grade, and alcohol history, suggesting SLC2A3 functioned as an oncogene in HNSCC. A previous study revealed that siRNA-mediated silencing of SLC2A3 could suppress the growth and promote apoptosis in AML cells [27]. Further, high expression of SLC2A3 was associated with poor OS, DSS, and RFS compared to low expression. Similarly, the immunohistochemical analysis of SLC2A3 was conducted in OSCC tissues, and the results showed that positive expression of SLC2A3 was correlated to unfavorable prognosis outcomes [28]. We concluded that SLC2A3 could act as a biomarker for prognosis in HNSCC.

Many studies have demonstrated that SLC2A3 could serve as an oncogene because of the synergistic action with upstream and downstream molecules. In GBM, knockdown of HDAC2 could induce GBM cell death by miR-3189-targeted SLC2A3 repression through reprogramming of glucose metabolism [29]. YAP played a crucial role in promoting the proliferation and migration of CRC cells by activating the SLC2A3/AMPK signaling pathway and epithelial-to-mesenchymal transition (EMT) [30]. Other studies have shown that SLC2A3 was directly related to SOX9, TRIM66, and HMGA1, which could promote cancer cell proliferation, migration, and invasion [31–33]. Given the condition, we constructed the coexpression network and performed the enrichment analysis for SLC2A3 in HNSCC. Interestingly, it was found that SLC2A3 was significantly correlated to glycolysis-related gene HK2 and hypoxia-related gene HIF1A, which indicated that the expression of SLC2A3 in HNSCC was regulated via glycolysis and HIF-1 signaling pathway according to enrichment analysis. In prostate cancer, SLC2A3 was a direct target of miR-29c to inhibit cell growth and glucose metabolism [34]. Similar results showed that miR-129-5p directly targets SLC2A3 to suppress glucose consumption, lactate generation, cellular ATP levels, and glucose uptake in gastric cancer cells [35]. Furthermore, immune cells use the ATPs generated by increasing glycolysis to modulate functional activities, growth, volume, and ion integrity, and antiglycolytics limit immunity [36]. However, research based on bioinformatics showed that high glycolysis status significantly correlated with less immune function, further related to poor prognosis in patients with osteosarcoma [37]. Hypoxia is a common feature in the TME because the uncontrolled proliferation of tumor cells mismatches the relatively insufficient supply of blood as well as oxygen, leading to stress-tolerant tumor cell survival [38]. In hypoxic-cultured prostate cell lines and hypoxic regions of xenograft tumors, the expression of SLC2A3 was significantly increased [39]. Combining previous results with our present study, we speculate that SLC2A3 predicting unfavorable prognosis in HNSCC could associate with glucose metabolism and HIF-1 signaling pathway.

It has long been acknowledged that the comprehensive immunobiology of the TME is modulated by several factors, including TILs, location of the tumor, and tumor stromal and vasculature [40–44]. The relationship between SLC2A3 and immune infiltration in human cancers is not widely studied. A recent study observed that SLC2A3 was one of an-eight gene prognostic risk score model to predict the survival significantly related to immune cell infiltration [45]. In this study, we found that SLC2A3 strongly correlated to TILs in HNSCC by TIMER, such as B cell, CD8⁺ T cell, CD4⁺ T cell, macrophage, neutrophil, and dendritic cell. For example, macrophages have been proven to account for immune evasion and suppression, and CD8⁺ T cells could mediate the recognition of the antigens expressed by most tumor cells [46]. Interestingly, B cell was significantly associated with the cumulative survival in HNSCC between high expression and low expression of SLC2A3 suggesting that SLC2A3 could serve as an immune-related prognostic biomarker. We also demonstrated a positive correlation between SLC2A3 expression and Tcm CD8, Th1, Treg, Mem B, NK, NKT, pDC, eosinophil, and mast cells. In prostate cancer, SLC2A3 was one of the high-risk genes, and the risk score was significantly related to NK cell infiltration [47]. While functional Treg was generated from naive T cells lacking Glut1, despite expressing a low level of SLC2A3 which mediates glucose uptake [48]. Then, the expression of SLC2A3 was explored, and the results showed significantly differential expression of SLC2A3 in different immune and molecular subtypes. A previous study demonstrated that glucose metabolism could be applied in conjunction with several inflammatory diseases like rheumatoid arthritis, in which rheumatoid arthritis fibroblast-like synoviocytes were cultured with CD4 T cells or T cell conditioned medium. Then it was observed that the secretion of inflammatory genes including TNF-α, IFN-γ, IL-17A, IL-6, and IL-8 was augmented, and glucose transporters like GLUT1 and SLC2A3 were increased. [49]. According to these results, we hypothesized that SLC2A3 could play an essential role in HNSCC immune microenvironment. Moreover, ICP genes are important immune regulators in maintaining self-tolerance and modulating the immune response by both stimulatory and inhibitory pathways [50]. T cell-associated inhibitory molecules like lymphocyte activation gene 3 (LAG3) could suppress the antitumor immune responses, while upregulation of inducible costimulator could contribute to favorable clinical outcomes [51]. When it comes to SLC2A3, a recent study observed that SLC2A3 knockdown significantly inhibited the expression of PD-L1 indicating SLC2A3 could participate in the immune response by regulating the PD-L1 immune checkpoint in CRC [10]. Similarly, upregulation of SLC2A3 could enhance CD28 costimulation reprogrammed renal cell carcinoma CD8⁺ TIL metabolism, which is antagonized by the inhibitory and checkpoint immunotherapy receptors CTLA4 and PD-1 [52]. We found a strong relationship between SLC2A3 and ICP genes which provided a novel insight for targeted immunotherapy in HNSCC. In addition, we found that SLC2A3 was positively related to stromal, immune, and ESTIMATE score associated with tumor invasion.

However, we used diverse online databases and R software to cross-validate the results systematically in our study, but several issues remain to be addressed. First, it would be of value to perform in vitro/in vivo studies to verify these results. Second, an in-depth analysis of the molecular mechanism of SLC2A3 in HNSCC is required. Moreover, further research is deserved for the lack of direct evidence of the association between SLC2A3 and immune regulation in HNSCC.

5. Conclusion

In summary, our findings indicated an intimate connection between SLC2A3 expression and HNSCC. We uncovered a role for SLC2A3 as a potential prognostic marker and immune therapeutic target in HNSCC via bioinformatics methods for the first time, and it suggested avenues for the precise treatment of patients with HNSCC.

Data Availability

The data analyzed in this study are available in the public online website have been mentioned in this paper.

Conflicts of Interest

There are no conflicts of interest among all the authors.

Authors’ Contributions

Manru Chu and Ke Zheng contributed to the conception and manuscript writing of this article and were responsible for data preparation, analyses, and interpretation. Xiaojie Li and Zhiqiang Luo revised the manuscript at different times critically. Xin Yang provided important guidance for research. Changbo Wei was responsible for funding acquisition and project administration. All authors have read and approve to publish the final version of the manuscript. Manru Chu and Ke Zheng contributed equally to this work.

Acknowledgments

This study was supported by Suzhou Science and Technology Projects (SYSD2020071), Science Foundation of Suzhou Stomatological Hospital (SZKQYY20A001), Science and Technology Development Fund of Nanjing Medical University (NMUB20210303) and Fund of Wuxi Health Commission (Q202137).

References

R. L. Siegel, K. D. Miller, and A. Jemal, “Cancer statistics, 2020,” CA: A Cancer Journal for Clinicians, vol. 70, no. 1, pp. 7–30, 2020.
View at: Publisher Site | Google Scholar
L. H. Chen, K.-L. Tsai, Y.-W. Chen et al., “MicroRNA as a novel modulator in head and neck squamous carcinoma,” Journal of Oncology, vol. 2010, Article ID 135632, pp. 1–15, 2010.
View at: Publisher Site | Google Scholar
A. Tandon, K. Sandhya, N. N. Singh, and A. Kumar, “Prognostic relevance of lymphatic vessel density in squamous cell carcinomas of the oral cavity: a systematic review and meta-analysis,” A Systematic Review and Meta-Analysis. Head Neck Pathology, vol. 16, 2022.
View at: Publisher Site | Google Scholar
J. Hu, D. Qiu, A. Yu et al., “YTHDF1 is a potential pan-cancer biomarker for prognosis and immunotherapy,” Frontiers in Oncology, vol. 11, article 607224, 2021.
View at: Publisher Site | Google Scholar
L. Liang, Q. Huang, M. Gan et al., “HighSEC61Gexpression predicts poor prognosis in patients with head and neck squamous cell carcinomas,” Journal of Cancer, vol. 12, no. 13, pp. 3887–3899, 2021.
View at: Publisher Site | Google Scholar
G. Lu, Y. Zhou, C. Zhang, and Y. Zhang, “Upregulation of LIMK1 is correlated with poor prognosis and immune infiltrates in lung adenocarcinoma,” Frontiers in Genetics, vol. 12, article 671585, 2021.
View at: Publisher Site | Google Scholar
B. Lv, Y. Wang, D. Ma et al., “Immunotherapy: reshape the tumor immune microenvironment,” Frontiers in Immunology, vol. 13, article 844142, 2022.
View at: Publisher Site | Google Scholar
D. Jia, M. Lu, K. H. Jung et al., “Elucidating cancer metabolic plasticity by coupling gene regulation with metabolic pathways,” Proceedings of the National Academy of Sciences of the United States of America, vol. 116, no. 9, pp. 3909–3918, 2019.
View at: Google Scholar
D. Deng, P. Sun, C. Yan et al., “Molecular basis of ligand recognition and transport by glucose transporters,” Nature, vol. 526, no. 7573, pp. 391–396, 2015.
View at: Publisher Site | Google Scholar
H. Gao, J. Liang, J. Duan et al., “A prognosis marker SLC2A3 correlates with EMT and immune signature in colorectal cancer,” Frontiers in Oncology, vol. 11, article 638099, 2021.
View at: Publisher Site | Google Scholar
X. Yao, Z. He, C. Qin et al., “SLC2A3 promotes macrophage infiltration by glycolysis reprogramming in gastric cancer,” Cancer Cell International, vol. 20, no. 1, p. 503, 2020.
View at: Publisher Site | Google Scholar
T. H. Tsai, C. C. Yang, T. C. Kou et al., “Overexpression of GLUT3 promotes metastasis of triple-negative breast cancer by modulating the inflammatory tumor microenvironment,” Journal of Cellular Physiology, vol. 236, no. 6, pp. 4669–4680, 2021.
View at: Publisher Site | Google Scholar
T. Li, J. Fan, B. Wang et al., “TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells,” Cancer Research, vol. 77, no. 21, pp. e108–e110, 2017.
View at: Publisher Site | Google Scholar
D. S. Chandrashekar, B. Bashel, S. A. H. Balasubramanya et al., “UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses,” Neoplasia, vol. 19, no. 8, pp. 649–658, 2017.
View at: Publisher Site | Google Scholar
E. Sjostedt, W. Zhong, L. Fagerberg et al., “An atlas of the protein-coding genes in the human, pig, and mouse brain,” Science, vol. 367, no. 6482, 2020.
View at: Publisher Site | Google Scholar
A. Lanczky and B. Gyorffy, “Web-based survival analysis tool tailored for medical research (KMplot): development and implementation,” Journal of Medical Internet Research, vol. 23, no. 7, article e27633, 2021.
View at: Publisher Site | Google Scholar
J. Gao, B. A. Aksoy, U. Dogrusoz et al., “Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal,” Science Signaling, vol. 6, no. 269, p. pl1, 2013.
View at: Publisher Site | Google Scholar
S. V. Vasaikar, P. Straub, J. Wang, and B. Zhang, “LinkedOmics: analyzing multi-omics data within and across 32 cancer types,” Nucleic Acids Research, vol. 46, no. D1, pp. D956–D963, 2018.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, K. C. Nastou et al., “The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets,” Nucleic Acids Research, vol. 49, no. D1, pp. D605–D612, 2021.
View at: Publisher Site | Google Scholar
B. Ru, C. N. Wong, Y. Tong et al., “TISIDB: an integrated repository portal for tumor-immune system interactions,” Bioinformatics, vol. 35, no. 20, pp. 4200–4202, 2019.
View at: Publisher Site | Google Scholar
K. Yoshihara, M. Shahmoradgoli, E. Martínez et al., “Inferring tumour purity and stromal and immune cell admixture from expression data,” Nature Communications, vol. 4, no. 1, p. 2612, 2013.
View at: Publisher Site | Google Scholar
R. A. Medina and G. I. Owen, “Glucose transporters: expression, regulation and cancer,” Biological Research, vol. 35, no. 1, pp. 9–26, 2002.
View at: Publisher Site | Google Scholar
G. C. Ziegler, P. Almos, R. V. McNeill, C. Jansch, and K. P. Lesch, “Cellular effects and clinical implications of SLC2A3 copy number variation,” Journal of Cellular Physiology, vol. 235, no. 12, pp. 9021–9036, 2020.
View at: Publisher Site | Google Scholar
J. Liu, J. Hong, H. Han et al., “Decreased vitamin C uptake mediated by _SLC2A3_ promotes leukaemia progression and impedes TET2 restoration,” British Journal of Cancer, vol. 122, no. 10, pp. 1445–1452, 2020.
View at: Publisher Site | Google Scholar
D. W. Dai, Q. Lu, L. X. Wang et al., “Decreased miR-106a inhibits glioma cell glucose uptake and proliferation by targeting SLC2A3 in GBM,” BMC Cancer, vol. 13, no. 1, p. 478, 2013.
View at: Publisher Site | Google Scholar
E. Kim, S. Jung, W. S. Park et al., “Upregulation of SLC2A3 gene and prognosis in colorectal carcinoma: analysis of TCGA data,” BMC Cancer, vol. 19, no. 1, p. 302, 2019.
View at: Publisher Site | Google Scholar
Y. Zhuang, J. Zhao, X. Xu, and L. Bi, “Downregulation of GLUT3 promotes apoptosis and chemosensitivity of acute myeloid leukemia cells via EGFR signaling,” Archives of Iranian Medicine, vol. 21, no. 2, pp. 73–78, 2018.
View at: Google Scholar
F. R. Ayala, R. M. Rocha, K. C. Carvalho et al., “GLUT1 and GLUT3 as potential prognostic markers for oral squamous cell carcinoma,” Molecules, vol. 15, no. 4, pp. 2374–2387, 2010.
View at: Publisher Site | Google Scholar
S. Kwak, S. H. Park, S. H. Kim et al., “miR-3189-targeted GLUT3 repression by HDAC2 knockdown inhibits glioblastoma tumorigenesis through regulating glucose metabolism and proliferation,” Journal of Experimental & Clinical Cancer Research, vol. 41, no. 1, p. 87, 2022.
View at: Publisher Site | Google Scholar
L. Jiang, J. Zhang, Q. Xu et al., “YAP promotes the proliferation and migration of colorectal cancer cells through the Glut3/AMPK signaling pathway,” Oncology Letters, vol. 21, no. 4, p. 312, 2021.
View at: Publisher Site | Google Scholar
S. Chen, L. Yang, Z. Li et al., “EGFR/EGFRvIII partly regulates the tumourigenesis of glioblastoma through the SOX9-GLUT3 axis,” American Journal of Translational Research, vol. 13, no. 6, pp. 6055–6065, 2021.
View at: Google Scholar
Y. Song, L. Meng, J. Yu, Z. Cao, J. Sun, and H. Zhao, “TRIM66 overexpression promotes glioma progression and regulates glucose uptake through cMyc/GLUT3 signaling,” Cancer Management and Research, vol. 13, pp. 5187–5201, 2021.
View at: Publisher Site | Google Scholar
M. Yang, Y. Guo, X. Liu, and N. Liu, “HMGA1 promotes hepatic metastasis of colorectal cancer by inducing expression of glucose transporter 3 (GLUT3),” Medical Science Monitor, vol. 26, article e924975, 2020.
View at: Publisher Site | Google Scholar
J. Li, F. Fu, X. Wan, S. Huang, D. Wu, and Y. Li, “Up-regulated miR-29c inhibits cell proliferation and glycolysis by inhibiting SLC2A3 expression in prostate cancer,” Gene, vol. 665, pp. 26–34, 2018.
View at: Publisher Site | Google Scholar
D. Chen, H. Wang, J. Chen et al., “MicroRNA-129-5p regulates glycolysis and cell proliferation by targeting the glucose transporter SLC2A3 in gastric cancer cells,” Frontiers in Pharmacology, vol. 9, p. 502, 2018.
View at: Publisher Site | Google Scholar
K. S. Gill, P. Fernandes, T. R. O’Donovan et al., “Glycolysis inhibition as a cancer treatment and its role in an anti-tumour__immune response,” Biochimica et Biophysica Acta-Bioenergetics, vol. 1866, pp. 87–105, 2016.
View at: Publisher Site | Google Scholar
K. Tian, W. Qi, Q. Yan, M. Lv, and D. Song, “Signature constructed by glycolysis-immune-related genes can predict the prognosis of osteosarcoma patients,” Investigational New Drugs, vol. 40, no. 4, pp. 818–830, 2022.
View at: Publisher Site | Google Scholar
E. E. Wicks and G. L. Semenza, “Hypoxia-inducible factors: cancer progression and clinical translation,” The Journal of Clinical Investigation, vol. 132, no. 11, 2022.
View at: Publisher Site | Google Scholar
J. M. Ryniawec, M. R. Coope, E. Loertscher et al., “GLUT3/SLC2A3 is an endogenous marker of hypoxia in prostate cancer cell lines and patient-derived xenograft tumors,” Diagnostics (Basel), vol. 12, no. 3, p. 676, 2022.
View at: Publisher Site | Google Scholar
F. Azimi, R. A. Scolyer, P. Rumcheva et al., “Tumor-infiltrating lymphocyte grade is an independent predictor of sentinel lymph node status and survival in patients with cutaneous melanoma,” Journal of Clinical Oncology, vol. 30, no. 21, pp. 2678–2683, 2012.
View at: Publisher Site | Google Scholar
J. Li, K. T. Byrne, F. Yan et al., “Tumor cell-intrinsic factors underlie heterogeneity of immune cell infiltration and response to immunotherapy,” Immunity, vol. 49, no. 1, pp. 178–193.e7, 2018.
View at: Publisher Site | Google Scholar
B. L. Horton, T. B. Fessenden, and S. Spranger, “Tissue site and the cancer immunity cycle,” Trends Cancer, vol. 5, no. 10, pp. 593–603, 2019.
View at: Publisher Site | Google Scholar
M. Yamauchi, T. H. Barker, D. L. Gibbons, and J. M. Kurie, “The fibrotic tumor stroma,” The Journal of Clinical Investigation, vol. 128, no. 1, pp. 16–25, 2018.
View at: Publisher Site | Google Scholar
M. B. Schaaf, A. D. Garg, and P. Agostinis, “Defining the role of the tumor vasculature in antitumor immunity and immunotherapy,” Cell Death & Disease, vol. 9, no. 2, p. 115, 2018.
View at: Publisher Site | Google Scholar
B. Wu, L. Tao, D. Yang, W. Li, H. Xu, and Q. He, “Development of an immune infiltration-related eight-gene prognostic signature in colorectal cancer microenvironment,” BioMed Research International, vol. 2020, Article ID 2719739, pp. 1–43, 2020.
View at: Publisher Site | Google Scholar
T. F. Gajewski, H. Schreiber, and Y. X. Fu, “Innate and adaptive immune cells in the tumor microenvironment,” Nature Immunology, vol. 14, no. 10, pp. 1014–1022, 2013.
View at: Publisher Site | Google Scholar
L. Zhang, B. Wang, Z.-S. Wang, Y.-L. Guo, and H. Shen, “Construction of glycolytic regulator gene signature to predict the prognosis and tumor immune cell infiltration levels for prostate cancer,” Computational and Mathematical Methods in Medicine, vol. 2022, Article ID 9273559, pp. 1–13, 2022.
View at: Publisher Site | Google Scholar
A. N. Macintyre, V. A. Gerriets, A. G. Nichols et al., “The glucose transporter Glut1 is selectively essential for CD4 T cell activation and effector function,” Cell Metabolism, vol. 20, no. 1, pp. 61–72, 2014.
View at: Publisher Site | Google Scholar
A. Petrasca, J. J. Phelan, S. Ansboro, D. J. Veale, U. Fearon, and J. M. Fletcher, “Targeting bioenergetics prevents CD4 T cell-mediated activation of synovial fibroblasts in rheumatoid arthritis,” Rheumatology (Oxford), vol. 59, no. 10, pp. 2816–2828, 2020.
View at: Publisher Site | Google Scholar
J. A. Marin-Acevedo, B. Dholaria, A. E. Soyano, K. L. Knutson, S. Chumsri, and Y. Lou, “Next generation of immune checkpoint therapy in cancer: new developments and challenges,” Journal of Hematology & Oncology, vol. 11, no. 1, p. 39, 2018.
View at: Publisher Site | Google Scholar
A. Hanson, K. Elpek, E. Duong et al., “ICOS agonism by JTX-2011 (vopratelimab) requires initial T cell priming and fc cross-linking for optimal T cell activation and anti-tumor immunity in preclinical models,” PLoS One, vol. 15, no. 9, article e0239595, 2020.
View at: Publisher Site | Google Scholar
K. E. Beckermann, R. Hongo, X. Ye et al., “CD28 costimulation drives tumor-infiltrating T cell glycolysis to promote inflammation,” Insight, vol. 5, no. 16, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Manru Chu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

391

Downloads

565

Citations