Next Article in Journal
Glial Cultures Differentiated from iPSCs of Patients with PARK2-Associated Parkinson’s Disease Demonstrate a Pro-Inflammatory Shift and Reduced Response to TNFα Stimulation
Next Article in Special Issue
Gain of Aggressive Histological and Molecular Patterns after Acquired Resistance to Novel Anti-EGFR Therapies in Non-Small Cell Lung Cancer
Previous Article in Journal
Coordinated Loss and Acquisition of NK Cell Surface Markers Accompanied by Generalized Cytokine Dysregulation in COVID-19
Previous Article in Special Issue
Breast Cancer with Increased Drug Resistance, Invasion Ability, and Cancer Stem Cell Properties through Metabolism Reprogramming
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Glycolysis-Related Gene Analyses Indicate That DEPDC1 Promotes the Malignant Progression of Oral Squamous Cell Carcinoma via the WNT/β-Catenin Signaling Pathway

1
State Key Laboratory of Oral Diseases, National Clinical Research Center for Oral Diseases, West China Hospital of Stomatology, Sichuan University, Chengdu 610044, China
2
Department of Head and Neck Oncology, West China Hospital of Stomatology, Sichuan University, Chengdu 610044, China
3
Division of Oral Ecology and Biochemistry, Tohoku University Graduate School of Dentistry, Sendai 980-8575, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(3), 1992; https://doi.org/10.3390/ijms24031992
Submission received: 16 December 2022 / Revised: 13 January 2023 / Accepted: 17 January 2023 / Published: 19 January 2023
(This article belongs to the Special Issue Advance in Targeted Cancer Therapy and Mechanisms of Resistance)

Abstract

:
Increasing evidence suggests that aerobic glycolysis is related to the progression of oral squamous cell carcinoma (OSCC). Hence, we focused on glycolysis-related gene sets to screen for potential therapeutic targets for OSCC. The expression profiles of OSCC samples and normal controls were obtained from The Cancer Genome Atlas (TCGA). Then, the differentially expressed gene sets were selected from the official GSEA website following extraction of the differentially expressed core genes (DECGs). Subsequently, we tried to build a risk model on the basis of DECGs to predict the prognosis of OSCC patients via Cox regression analysis. Furthermore, crucial glycolysis-related genes were selected to explore their biological roles in OSCC. Two active glycolysis-related pathways were acquired and 66 DECGs were identified. Univariate Cox regression analysis showed that six genes, including HMMR, STC2, DDIT4, DEPDC1, SLC16A3, and AURKA, might be potential prognostic factors. Subsequently, a risk formula consisting of DEPDC1, DDIT4, and SLC16A3 was established on basis of the six molecules. Furthermore, DEPDC1 was proven to be related to advanced stage cancer and lymph node metastasis. Moreover, functional experiments suggested that DEPDC1 promoted the aerobic glycolysis, migration, and invasion of OSCC via the WNT/β-catenin pathway. The risk score according to glycolysis-related gene expression might be an independent prognostic factor in OSCC. In addition, DEPDC1 was identified as playing a carcinogenic role in OSCC progression, suggesting that DEPDC1 might be a novel biomarker and therapeutic target for OSCC.

1. Introduction

OSCC is one of the prevalent types of head and neck squamous cell carcinomas. Recently, OSCC has attracted increasing attention owing to its severe impact on patient quality of life. With advances in diagnosis and operative treatment, the prognosis and quality of life have significantly improved. However, for some advanced malignancies, the 5-year overall survival is just under 50% [1]. In addition, lymph node metastasis and local recurrence contribute to the poor prognosis. Studies have indicated that the overall survival of patients with early-stage OSCC could be improved to more than 85% [2]. Nevertheless, there are few early effective biomarkers and therapeutic targets for OSCC. Therefore, it is vital to explore the mechanism of OSCC progression and identify effective biomarkers to improve the overall survival for advanced-stage OSCC patients.
Cancer cells are characterized by infinite proliferation. Metabolism includes a variety of interconnected cellular chemical transformations that provide energy to sustain life [3]. To meet their abundant energy demands, cancer cells show metabolic plasticity in the reprogramming of metabolic pathways such as glycolysis [4]. In 1956, Warburg indicated that cancer cells performed glycolysis to acquire energy even under aerobic conditions, which is called the Warburg effect or aerobic glycolysis [5]. Aerobic glycolysis leads to an increase in glucose uptake and accumulation of lactate in cancer cells. Meanwhile, it also speeds up the production of adenosine triphosphate (ATP) [6]. Studies have demonstrated that a faster energy supply through aerobic glycolysis contributed to the proliferation of cancer cells, and the acidic microenvironment formed by accumulation of lactate also facilitated cancer invasion and metastasis [7]. At present, many studies have shown that aerobic glycolysis plays a crucial role in a variety of cancers. For instance, the activation of aerobic glycolysis in pancreatic cancer cells was related to their growth [8]. In addition, accelerated glycolysis also promoted the metastasis of breast cancer [9]. In OSCC, previous studies have indicated that lactate dehydrogenase A, which catalyzes the conversion of L-lactate to pyruvate, could promote malignant progression and the epithelial–mesenchymal transition (EMT) by facilitating glycolysis [10]. Meanwhile, activation of aerobic glycolysis induced by protein disulfide isomerase family 6 promoted the progression of OSCC [11]. In addition, aerobic glycolysis was also relevant to the treatment of OSCC with compounds such as Tanshinone IIA [12] and Nobiletin [13]. Therefore, aerobic glycolysis may be significantly associated with tumorigenesis, progression, and drug therapy in OSCC. However, there are few studies which focus on glycolysis-related pathways and genes to explore their prognostic values.
In our study, we focused on glycolysis-related pathways and carried out Cox regression analysis to explore their prognostic value. As a result, the glycolysis-related biomarker DEP domain containing 1 (DEPDC1), which is a part of a transcription repressor complex and a newly discovered tumor-associated gene [14], was identified as being associated with the prognosis of OSCC. Moreover, the knockdown of DEPDC1 significantly inhibited migration, invasion, glucose uptake and accumulation of lactate. In addition, we demonstrated that DEPDC1 promoted OSCC progression via the WNT/β-catenin signaling pathway. These findings suggest that DEPDC1 may be a viable glycolysis-related biomarker and therapeutic target for OSCC.

2. Results

2.1. Glycolysis-Related Gene Sets in OSCC

The gene expression profiles and corresponding clinical data of 319 OSCC samples and 32 normal controls were obtained from TCGA. In addition, four glycolysis-related pathways including the biocarta-glycolysis-pathway, kegg-glycolysis-gluconeogenesis, hallmark-glycolysis, and reactome-glycolysis were gained from the official GSEA website. Bioinformatics analyses indicated that the gene sets hallmark-glycolysis (Figure 1A), and reactome-glycolysis (Figure 1B) were significantly activated in OSCC. To explore the prognostic value of the two gene sets, the core genes were screened out (Table 1). Furthermore, differential analyses indicated that there were 66 DEGs including 19 downregulated and 49 upregulated genes in OSCC (Figure 1C,D, Table 2 and Table 3).

2.2. Identification of Prognostic Cox Model

To investigate the prognostic value of glycolysis-related gene sets, 66 DEGs were enrolled in univariate Cox regression, which indicated that six molecules might be candidate biomarkers (Table S1). Furthermore, multivariate Cox regression and stepwise Cox analyses were used to construct the prognostic model based on the expression of six molecules. The results showed that three genes, including DNA damage inducible transcript 4 (DDIT4), DEPDC1 and solute carrier family 16 member 3 (SLC16A3), were all verified to be upregulated in TCGA OSCC samples (Figure 2A–C), and the risk-score= (DDIT4× 0.0022) + (DEPDC1× 0.0788) + (SLC16A3× 0.0316) (Table S2). Subsequently, to explore the prognostic value of the Cox model, we divided OSCC samples into high- and low-risk groups according to the median value of the risk score. Survival analysis showed that low-risk group patients had better prognosis than high-risk group patients (Figure 2D). In addition, similar results were also observed in stratification analyses (Figure S1). Moreover, the risk score might be visualized as an independent prognostic factor, in line with the univariate and multivariate Cox regression analyses (Figure 2E,F). Meanwhile, a high risk score meant poor differentiation in OSCC (Figure 2G). To intuitively understand the difference between high- and low-risk groups, the expression of risk genes combined with survival time of each patient was visualized in R software (Version R4.2.1, Chengdu, China) (Figure S2A–C).

2.3. DEPDC1 Was Identified as a Potential Biomarker in OSCC

To further understand the prognostic value of the three risk genes, survival analyses were carried out according to the expression of each gene. The results indicated that patients with high DEPDC1 expression had worse overall survival (Figure 3A). Meanwhile, high expression of DEPDC1 was related to lymph node metastasis (Figure 3B) and advanced stage (Figure 3C). Therefore, DEPDC1 might be an advanced biomarker in OSCC. To further explore the potential mechanism of DEPDC1, we performed GSEA and showed that high expression of DEPDC1 was significantly relevant to the activation of the WNT/β-catenin signaling pathway (Figure 3D,E). In addition, the carcinogenic roles of DEPDC1 in OSCC were further explored via GSEA (Figure S3).

2.4. Knockdown of DEPDC1 Inhibited the Migration and Invasion of OSCC

To determine the tumorigenic role of DEPDC1 in vitro, siRNAs targeting DEPDC1 and their corresponding negative controls were transfected into SCC9 and SCC15. The efficiency of knockdown was identified by qRT-PCR and Western blot (Figure S4A,B). As shown in the results, a decrease in DEPDC1 expression inhibited OSCC migration and invasion (Figure 4A,B). Meanwhile, knockdown of DEPDC1 decreased the protein expression level of β-catenin and N-cadherin and enhanced the expression of E-cadherin (Figure 4C).

2.5. DEPDC1 Promoted the Progression of OSCC via the WNT/β-Catenin Pathway

GSEA showed that high DEPDC1 expression was related to activation of the WNT/β-catenin pathway (Figure 3D,E). To further determine the effect of DEPDC1 in OSCC, the overexpression vectors of DEPDC1 and their controls were transfected into OSCC cells (Figure S4C,D), followed by treatment with XAV939 (Selleck Chemicals, Chengdu, China), an inhibitor of the WNT/β-catenin pathway, at a concentration of 10 μM. The results indicated that overexpression of DEPDC1 promoted the migration and invasion of OSCC. However, the migration and invasion abilities were reversed by XAV939 (Figure 5A,B). Meanwhile, upregulation of DEPDC1 increased the protein levels of β-catenin and N-cadherin and decreased E-cadherin expression (Figure 5C). The protein levels were also reversed by XAV939.

2.6. DEPDC1 Is Related to Glucose Uptake and Accumulation of Lactate in OSCC

To explore whether DEPDC1 influences aerobic glycolysis, we detected the glucose uptake and lactate production in OSCC cells after knockdown or overexpression of DEPDC1. As shown in the figure, knockdown of DEPDC1 decreased the glucose uptake and accumulation of lactate (Figure 6A). This phenomenon was reversed after overexpression of DEPDC1. Nevertheless, the activation of aerobic glycolysis by upregulated DEPDC1 was abolished following treatment with XAV939 (Figure 6B). In addition, knockdown of DEPDC1 decreased the expression of glucose transporter 1 (Glut1), a major glucose transporter, while overexpression of DEPDC1 enhanced the Glut1 protein level, which was also abolished by XAV939 (Figure 6C,D).

3. Discussion

Generally, cells depend on mitochondrial oxidative phosphorylation to produce energy for cellular processes. However, cancer cells always demand large amounts of energy for their fast proliferation and metastasis [15]. Therefore, cancer cells reprogram metabolic pathways including aerobic glycolysis, glutamine metabolism, and fatty acid metabolism [4]. Among them, aerobic glycolysis, also called the “Warburg effect”, is a common method of metabolic reprogramming [7]. This process leads to an increase in glucose uptake, accumulation of lactate and faster production of ATP, which plays a crucial role in tumorigenesis and tumor progression [6]. In addition, glycolysis helps cancer cells survive under cellular stress and accelerated aerobic glycolysis is related to the treatment resistance of cancers [16]. Hence, these mechanisms need to be further explored in OSCC, which might be helpful for improving prognosis.
Owing to the lack of effective early biomarkers, the 5-year overall survival of OSCC is still poor. Recently, a number of studies have indicated that aerobic glycolysis is associated with the proliferation [17], invasion [18], metastasis [19] and EMT [20] of OSCC. Therefore, it may be appropriate to explore glycolysis-related pathways to screen for potential biomarkers of OSCC. In our study, two active glycolysis-related pathways were filtered following extraction of 255 core genes and 66 DEGs. The hallmark-glycolysis gene set includes genes encoding proteins involved in glycolysis and gluconeogenesis, while reactome-glycolysis is related to glycolysis and canonical pathways. According to the DEGs from two gene sets, a Cox risk model including DDIT4, DEPDC1, and SLC16A3 was constructed, which predicted the prognosis of OSCC patients accurately. In addition, the risk score might be an independent prognostic factor, indicating that the three genes might play a crucial carcinogenic role in the progression of OSCC. DDIT4, located in the cytosol, is related to the response to hypoxia, which plays a vital role in aerobic glycolysis in the tumor microenvironment. Recently, Han et al. also showed that DDIT4 could predict the prognosis of OSCC patients [21]. In addition, DDIT4 is known to participate in a variety of cancers. In colorectal cancer, DDIT4 has been identified as an advanced stage and metastasis biomarker [22]. Overexpression of DDIT4 could accelerate the proliferation and tumorigenesis of gastric cancer. Generally, inhibition of DDIT4 inactivated P53 signaling pathways and mitogen-activated protein kinase (MAPK) activity, which was associated with drug resistance to 5-fluorouracil [23]. In addition, DDIT4 expression could be induced by nicotinamide nucleotide transhydrogenase antisense RNA 1 (NNT-AS1)/miR-496 axis in prostate cancer and played a carcinogenic role in proliferation and migration [24]. Similarly, SLC16A3, also known as monocarboxylate transporter 3 (MCT3) or MCT4, was associated with lactate and pyruvate transport across plasma membranes, which played a crucial role in diverse malignant tumors. For instance, SLC16A3, mediated by LINC00035/CCAAT enhancer binding protein beta (CEBPB), enhanced the glycolysis and promoted the development of ovarian cancer [25]. Levels of aberrant N6−methylation of adenosine (m6A) at SLC16A3 were also relevant to the efficacy of immunotherapy in melanoma [26]. In the present study, DEPDC1 was identified as a glycolysis-related biomarker via integrated bioinformatics. DEPDC1 is a newly identified cancer-associated gene that is part of a transcription repressor complex and related to cellular processes such as the cell cycle, transcription, mitosis and apoptosis [27]. Recently, DEPDC1 was identified as participating in tumorigenesis and cancer progression. Guo et al. indicated that DEPDC1 was overexpressed in OSCC tissues and related to overall survival [28]. They also indicated that DEPDC1 promoted the proliferation of OSCC via inhibiting cytochrome P450 family 27 subfamily B member 1 (CYP27B1) expression [28]. In addition, Qiu et al. showed that DEPDC1 could facilitate development and metastasis of OSCC [29]. In our study, we demonstrated that DEPDC1 might target the WNT/β−catenin signaling pathway, and then promoted the migration, invasion, and aerobic glycolysis of OSCC. Geng et al. also demonstrated that DEPDC1 could promote nephroblastoma progression through the WNT/β-catenin signaling pathway [30]. Moreover, high DEPDC1 expression was also shown in advanced stage cancer and lymph nodes metastasis, suggesting that DEPDC1 might be an advanced biomarker in OSCC. Unfortunately, we failed to validate these data in real patient samples. Generally, there may be a small difference between the database and real patients. Therefore, a further study should be done to clarify their role in OSCC. In addition, DEPDC1 is known to be associated with various biological processes of other cancers. For instance, DEPDC1 promoted cancer metastasis and differentiation via accelerating the cell cycle from G1 to S phase in gastric cancer [31]. In hepatocellular carcinoma, DEPDC1, as a metabolic gene related to glycolysis, has been validated as a contributor to the reconstruction of the tumor microenvironment [32]. In addition, high DEPDC1 expression mediated by LincRNA regulator of reprogramming (Linc-ROR) accelerated the development of hepatocellular carcinoma and angiogenesis [33]. Meanwhile, DEPDC1 also promoted chemotherapy resistance [34]. These results suggest that DEPDC1 might be a potential therapeutic target.
In our study, GSEA showed that DEPDC1 was relevant to various signaling pathways in OSCC. Recently, Zhao et al. demonstrated that DEPDC1 overexpression led to a significant promotion of proliferation by regulating the cell cycle in breast cancer, and that DEPDC1 was associated with a high activation of the PI3K-AKT-mTOR signaling pathway [35]. In addition, DEPDC1 could interact with E2F1 and promote its transcriptional activity in prostate cancer, resulting in the activation of the E2F signaling pathway to regulate the G1/S phase cell cycle transition and then increase cell proliferation [36]. In nasopharyngeal carcinoma, DEPDC1 knockdown downregulated various downstream targets, such as c−Myc, BCL2, MMP2 and MMP9, which participated in proliferation, tumorigenesis, and metastasis [37], suggesting that DEPDC1 might be a vital factor in these biological processes. In OSCC, higher activation of the PI3K-AKT-mTOR pathway was reported to accelerate progression by increasing aerobic glycolysis [38]. Similarly, the E2F signaling pathway was also highly active in OSCC tissues, and this high activation drove radioresistance in OSCC patients, resulting in worse recurrence-free survival of radiotherapy patients [39]. Moreover, many studies suggest that adipogenesis [40], fatty acid metabolism [41], the MYC pathway [42], and the MTORC1 pathway [43] play a crucial role in the initiation and progression of OSCC. In conclusion, all of these pathways were remarkably relevant to OSCC initiation, progression, and therapy resistance. DEPDC1 was identified as activating these pathways through GSEA. Therefore, DEPDC1 is expected to be a novel biomarker and therapeutic target in OSCC.

4. Materials and Methods

4.1. Data Acquisition from TCGA

TCGA is a landmark cancer genomics program which contains gene expression profiles and clinical data for 33 cancers [44]. All of the OSCC data were downloaded from the TCGA database. A total of 351 samples including 319 OSCC samples and 32 normal controls were obtained. Two OSCC samples were excluded because of insufficient clinical data.

4.2. Gene Set Enrichment Analysis (GSEA)

A total of four glycolysis-related pathways (biocarta-glycolysis-pathway, kegg-glycolysis-gluconeogenesis, hallmark-glycolysis, reactome-glycolysis) were obtained from the GSEA website (http://www.gseamsigdb.org/gsea/, accessed on 30 September 2021). Subsequently, GSEA was used to screen the activated glycolysis-related gene sets in OSCC. Furthermore, the glycolysis-related genes were extracted from the activated glycolysis-related gene sets and then we carried out differentially expressed analysis. In addition, GSEA was also used to investigate the potential pathways which were activated by DEPDC1 in OSCC.

4.3. Identification of Potential Prognostic Biomarker

Based on the differentially expressed genes (DEGs) from the glycolysis−related gene sets, univariate and multivariate Cox regression analyses were used to construct the Cox risk model and screen potential biomarkers to predict the prognosis of OSCC patients. Furthermore, the risk score of each patient was calculated on the basis of gene expression profiles, and independent prognosis analyses were carried out based on clinical parameters and risk score.

4.4. Survival Analysis

According to the risk score of each patient, OSCC samples were divided into high and low risk groups. Combined with the survival data in TCGA, survival analyses with the log rank p test were applied to verify the accuracy of the Cox model in R software with the survival package (Version R-4.2.1, Chengdu, China).

4.5. Cell Culture

The human OSCC cell lines SCC9 and SCC15 were purchased from ATCC. Cells were cultured in DMEM (Gibco, Cat#11995500TB, Chengdu, China) with 10% fetal bovine serum (ExCell Bio, Inc., Chengdu, China) and incubated at 37 °C with 5% CO2. The medium of all cells was exchanged about every 1–2 days according to the cell density.

4.6. RNA Isolation and qRT-PCR

The protocol for RNA isolation and quantitative real-time PCR (qRT-PCR) has been described previously [45]. Briefly, total RNA was extracted from OSCC cells with TRIZOL reagent (Vazyme, Chengdu, China) and then reversed to cDNA. Subsequently, relative expression profiles of target genes were detected by qRT-PCR with ChamQ Universal SYBR qPCR Master Mix kit (Vazyme Biotech Co., Ltd., Chengdu, China). The primer sequences were as follows: DEPDC1 forward, 5′-ATGCGTATGATTTCCCGAATGAG-3′; DEPDC1 reverse, 5′-CACAGCATAACACACATCGAGAA-3′; GAPDH forward, 5′-CGCTGAGTACGTCGTGGAGTC-3′; and GAPDH reverse, 5′-GCTGATGATCTTGAGGCTGTTGTC-3′.

4.7. Cell Transfection

The small interfering RNAs (siRNAs) and DEPDC1 expression vectors were designed and synthesized in TINGKE (Cheng, China). According to the lipofectamine 3000 protocol (Invitrogen, Cat# L3000-015, Chengdu, China), cells were seeded in 6-well or 12-well plates for 12–24 h at a density of nearly 70–90%, and approximately 2500 ng siRNAs or DEPDC1 overexpression vectors were transfected into OSCC cells. The transfection system would last for 8–12 h, after which the medium was exchanged with fresh medium. After incubation for 2–4 days, RNA and protein were collected for analysis. The siRNA sequence of DEPDC1 was as follows: sense 5′-3′: CGCAGCCCUCUUGCUAUUTT, antisense 5′-3′: AAUAGCAUAGAGGGCUGCGTT.

4.8. Transwell Assay

The protocol for the transwell assay has been described previously [46]. Briefly, 600 μL DMEM with 10% FBS was placed into the bottom wells, and about 5 × 104 cells resuspended in 200 μL serum-free medium were put in the upper chambers. In addition, the invasion assay was performed with 100 μL 1:8 diluted matrigel in the upper chamber. The device was cultured at 37 °C with 5% CO2 for 24–72 h. Finally, the cells in the upper chamber were removed and the migration and invasion cells were stained with crystal violet.

4.9. Glucose Uptake and Lactate Detection

The cells of glucose uptake and lactate were detected using a detection kit (Cat#ADS-W-TDX002, ADS-W-T009-96, Chengdu, China) following the manufacturer’s instructions. Briefly, the same amounts of cells, about 5 × 106 cells, were collected to obtain intracellular glucose and lactate. Cells were broken up by cell ultrasound equipment at 200 W for 3 s, with a 10 s pause between each cycle, for a total of 30 times, and then centrifuged at 12,000 rpm and 4 °C for 10 min to collect the supernatant. the absorbance of the supernatant at 450 nm (lactic acid production) and 520 nm (glucose uptake) were determined by a microplate reader. Higher absorbance meant a higher level of glucose uptake and accumulation of lactate.

4.10. Western Blotting

The Western blot methods have been elaborated upon [46]. In short, the proteins of the OSCC cell lines were extracted, separated by polyacrylamide gel electrophoresis (PAGE), and transferred to PVDF membranes which were then sealed in 5% skim milk for 1 h at room temperature. Subsequently, the membranes were incubated overnight with primary antibodies at 4 °C and second antibodies for 1 h at room temperature. Finally, ECL was used to analyze the protein level. The antibody information was as follows: GAPDH (Ray ANTIBODU, 1:2000, Cat#RM2002, Chengdu, China); DEPDC1 (Bioss, 1:1000, Cat#bs-6525R, Chengdu, China); N-cadherin (Proteintech, 1:2000, Cat# 22018-1-AP, Chengdu, China); E-cadherin (Proteintech, 1:2000, Cat#20874-1-AP, Chengdu, China); Glut1 (Wanleibio, 1:1000, Cat#WL03141, Chengdu, China); and β-catenin (Wanleibio, 1:1000, Cat#WL0962a, Chengdu, China).

4.11. Statiscal Analysis

SPSS software (IBM, Version 23.0, Armonk, NY, USA) was used to perform statistical analysis. The difference between two groups was investigated with a Student’s t-test, and the difference in one factor among multiple groups was explored with a one-way ANOVA. Results were expressed as the mean ± SD (standard deviation). Significant differences were considered at p < 0.05 *; p < 0.01 **; p < 0.001 ***; and p < 0.0001 ****.

5. Conclusions

A total of two glycolysis-related gene sets with aberrant activity in OSCC were identified following extraction of the core genes. Subsequently, we tried to construct a Cox model which could predict the prognosis of OSCC patients more accurately, and found that the risk score according DEPDC1, DDIT4 and SLC16A3 expression profiles might be an independent prognostic factor, indicating that these three genes might play a crucial role in the prognosis of OSCC. Moreover, DEPDC1 was identified as accelerating aerobic glycolysis, migration, and invasion via the WNT/β-catenin signaling pathway in vitro. Meanwhile, DEPDC1 was reported to be relevant to various cancer associated pathways. Our findings, along with previous reports, suggest that DEPDC1 might be a promising novel biomarker and therapeutic target in OSCC.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24031992/s1.

Author Contributions

Conceptualization: G.H., C.L. and N.T.; Methodology: G.H.; Formal analysis: S.C.; Software: G.P.L.; Writing—original draft preparation: G.H.; Writing—review and editing: G.H., C.L., J.W. and N.T.; Supervision: C.L., J.W. and N.T.; Funding acquisition: C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (81870782), the Science and Technology Department of Sichuan Province (2021YJ0129) the Outstanding Young Scholars foundation of Sichuan University (2082604194311), and the Exploration and Research Project of West China Hospital of Stomatology, Sichuan University (LCYJ2023-DL-1).

Data Availability Statement

All the data presented in this study are available on request from the corresponding author. Publicly available data were obtained from official TCGA (https://portal.gdc.cancer.gov/, accessed on 30 September 2021) and GSEA websites (https://www.gsea-msigdb.org/gsea/index.jsp, accessed on 30 September 2021).

Conflicts of Interest

All the authors declare that they have no competing interests.

References

  1. Kaminagakura, E.; Tango, R.N.; Cruz-Perez, D.; Bonan, R.; de Almeida, L.Y.; Lança, M.L.A.; Bonan, P.; Martins, H.; Takahama, A., Jr.; Ito, F.A.; et al. Oral squamous cell carcinoma outcome in adolescent/young adult: Systematic review and m eta-analysis. Head Neck 2022, 44, 548–561. [Google Scholar] [CrossRef] [PubMed]
  2. Takahashi, H.; Yanamoto, S.; Yamada, S.; Umeda, M.; Shigeta, T.; Minamikawa, T.; Shibuya, Y.; Komori, T.; Shiraishi, T.; Asahina, I.; et al. Effects of postoperative chemotherapy and radiotherapy on patients with squamous cell carcinoma of the oral cavity and multiple regional lymph node metastases. Int. J. Oral Maxillofac. Surg. 2014, 43, 680–685. [Google Scholar] [CrossRef] [PubMed]
  3. Judge, A.; Dodd, M.S. Metabolism. Essays Biochem. 2020, 64, 607–647. [Google Scholar] [CrossRef] [PubMed]
  4. Sun, L.; Suo, C.; Li, S.-T.; Zhang, H.; Gao, P. Metabolic reprogramming for cancer cells and their microenvironment: Beyond the Warburg Effect. Biochim. Biophys. Acta (BBA) Rev. Cancer 2018, 1870, 51–66. [Google Scholar] [CrossRef]
  5. Warburg, O. On the Origin of Cancer Cells. Science 1956, 123, 309–314. [Google Scholar] [CrossRef]
  6. Icard, P.; Shulman, S.; Farhat, D.; Steyaert, J.-M.; Alifano, M.; Lincet, H. How the Warburg effect supports aggressiveness and drug resistance of cancer cells? Drug Resist. Updates 2018, 38, 1–11. [Google Scholar] [CrossRef]
  7. Huang, P.; Zhu, S.; Liang, X.; Zhang, Q.; Luo, X.; Liu, C.; Song, L. Regulatory Mechanisms of LncRNAs in Cancer Glycolysis: Facts and Perspectives. Cancer Manag. Res. 2021, 13, 5317–5336. [Google Scholar] [CrossRef]
  8. Chen, L.; Yuan, R.; Wen, C.; Liu, T.; Feng, Q.; Deng, X.; Du, Y.; Peng, X. E3 ubiquitin ligase UBR5 promotes pancreatic cancer growth and aerobic glycolysis by downregulating FBP1 via destabilization of C/EBPα. Oncogene 2021, 40, 262–276. [Google Scholar] [CrossRef]
  9. Park, M.K.; Zhang, L.; Min, K.-W.; Cho, J.-H.; Yeh, C.-C.; Moon, H.; Hormaechea-Agulla, D.; Mun, H.; Ko, S.; Lee, J.W.; et al. NEAT1 is essential for metabolic changes that promote breast cancer growth and metastasis. Cell Metab. 2021, 33, 2380–2397.e9. [Google Scholar] [CrossRef]
  10. Cai, H.; Li, J.; Zhang, Y.; Liao, Y.; Zhu, Y.; Wang, C.; Hou, J. LDHA Promotes Oral Squamous Cell Carcinoma Progression through Facilitating Glycolysis and Epithelial–Mesenchymal Transition. Front. Oncol. 2019, 9, 1446. [Google Scholar] [CrossRef]
  11. Mao, L.; Wu, X.; Gong, Z.; Yu, M.; Huang, Z. PDIA6 contributes to aerobic glycolysis and cancer progression in oral squamous cell carcinoma. World J. Surg. Oncol. 2021, 19, 88. [Google Scholar] [CrossRef] [PubMed]
  12. Li, M.; Gao, F.; Zhao, Q.; Zuo, H.; Liu, W.; Li, W. Tanshinone IIA inhibits oral squamous cell carcinoma via reducing Akt-c-Myc signaling-mediated aerobic glycolysis. Cell Death Dis. 2020, 11, 381. [Google Scholar] [CrossRef]
  13. Lin, C.; Tu, C.; Ma, Y.; Ye, P.; Shao, X.; Yang, Z.; Fang, Y. Nobiletin inhibits cell growth through restraining aerobic glycolysis via PKA-CREB pathway in oral squamous cell carcinoma. Food Sci. Nutr. 2020, 8, 3515–3524. [Google Scholar] [CrossRef] [PubMed]
  14. Shen, L.; Li, H.; Liu, R.; Zhou, C.; Bretches, M.; Gong, X.; Lu, L.; Zhang, Y.; Zhao, K.; Ning, B.; et al. DEPDC1 as a crucial factor in the progression of human osteosarcoma. Cancer Med. 2022. [Google Scholar] [CrossRef]
  15. Vander Heiden, M.G.; Cantley, L.C.; Thompson, C.B. Understanding the Warburg Effect: The Metabolic Requirements of Cell Proliferation. Science 2009, 324, 1029–1033. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Lin, J.; Xia, L.; Liang, J.; Han, Y.; Wang, H.; Oyang, L.; Tan, S.; Tian, Y.; Rao, S.; Chen, X.; et al. The roles of glucose metabolic reprogramming in chemo- and radio-resistance. J. Exp. Clin. Cancer Res. 2019, 38, 218. [Google Scholar] [CrossRef] [Green Version]
  17. Zhang, X.; Dong, Y.; Zhao, M.; Ding, L.; Yang, X.; Jing, Y.; Song, Y.; Chen, S.; Hu, Q.; Ni, Y. ITGB2-mediated metabolic switch in CAFs promotes OSCC proliferation by oxidation of NADH in mitochondrial oxidative phosphorylation system. Theranostics 2020, 10, 12044–12059. [Google Scholar] [CrossRef]
  18. Zheng, M.; Cao, M.; Luo, X.; Li, L.; Wang, K.; Wang, S.; Wang, H.; Tang, Y.; Liang, X. EZH2 promotes invasion and tumour glycolysis by regulating STAT3 and FoxO1 signalling in human OSCC cells. J. Cell. Mol. Med. 2019, 23, 6942–6954. [Google Scholar] [CrossRef] [Green Version]
  19. Guo, J.; Su, Y.; Zhang, M. Circ_0000140 restrains the proliferation, metastasis and glycolysis metabolism of oral squamous cell carcinoma through upregulating CDC73 via sponging miR-182-5p. Cancer Cell Int. 2020, 20, 1–12. [Google Scholar] [CrossRef]
  20. Zhang, Y.; Cai, H.; Liao, Y.; Zhu, Y.; Wang, F.; Hou, J. Activation of PGK1 under hypoxic conditions promotes glycolysis and increases stem cell-like properties and the epithelial-mesenchymal transition in oral squamous cell carcinoma cells via the AKT signalling pathway. Int. J. Oncol. 2020, 57, 743–755. [Google Scholar] [CrossRef]
  21. Han, Y.; Wang, X.; Xia, K.; Su, T. A novel defined hypoxia-related gene signature to predict the prognosis of oral squamous cell carcinoma. Ann. Transl. Med. 2021, 9, 1565. [Google Scholar] [CrossRef] [PubMed]
  22. Fattahi, F.; Kiani, J.; Alemrajabi, M.; Soroush, A.; Naseri, M.; Najafi, M.; Madjd, Z. Overexpression of DDIT4 and TPTEP1 are associated with metastasis and advanced stages in colorectal cancer patients: A study utilizing bioinformatics prediction and experimental validation. Cancer Cell Int. 2021, 21, 303. [Google Scholar] [CrossRef] [PubMed]
  23. Du, F.; Sun, L.; Chu, Y.; Li, T.; Lei, C.; Wang, X.; Jiang, M.; Min, Y.; Lu, Y.; Zhao, X.; et al. DDIT4 promotes gastric cancer proliferation and tumorigenesis through the p53 and MAPK pathways. Cancer Commun. 2018, 38, 45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Yao, C.; Cheng, X.; Guo, X.; Lu, X.; Bu, F.; Xu, Y. NNT-AS1 modulates prostate cancer cell proliferation, apoptosis and migration through miR-496/DDIT4 axis. Cancer Cell Int. 2020, 20, 463. [Google Scholar] [CrossRef]
  25. Yang, S.; Wang, J.; Cheng, R.; Pang, B.; Sun, P. LINC00035 Transcriptional Regulation of SLC16A3 via CEBPB Affects Glycolysis and Cell Apoptosis in Ovarian Cancer. Evid. Based Complement. Altern. Med. 2021, 2021, 5802082. [Google Scholar] [CrossRef]
  26. Li, N.; Kang, Y.; Wang, L.; Huff, S.; Tang, R.; Hui, H.; Agrawal, K.; Gonzalez, G.M.; Wang, Y.; Patel, S.P.; et al. ALKBH5 regulates anti–PD-1 therapy response by modulating lactate and suppressive immune cell accumulation in tumor microenvironment. Proc. Natl. Acad. Sci. USA 2020, 117, 20159–20170. [Google Scholar] [CrossRef]
  27. Harada, Y.; Kanehira, M.; Fujisawa, Y.; Takata, R.; Shuin, T.; Miki, T.; Fujioka, T.; Nakamura, Y.; Katagiri, T. Cell-Permeable Peptide DEPDC1-ZNF224 Interferes with Transcriptional Repression and Oncogenicity in Bladder Cancer Cells. Cancer Res. 2010, 70, 5829–5839. [Google Scholar] [CrossRef] [Green Version]
  28. Guo, J.; Zhou, S.; Huang, P.; Xu, S.; Zhang, G.; He, H.; Zeng, Y.; Xu, C.-X.; Kim, H.; Tan, Y. NNK-mediated upregulation of DEPDC1 stimulates the progression of oral squamous cell carcinoma by inhibiting CYP27B1 expression. Am. J. Cancer Res. 2020, 10, 1745–1760. [Google Scholar]
  29. Qiu, J.; Tang, Y.; Liu, L.; Yu, J.; Chen, Z.; Chen, H.; Yuan, R. FOXM1 is regulated by DEPDC1 to facilitate development and metastasis of oral squamous cell carcinoma. Front. Oncol. 2022, 12, 815998. [Google Scholar] [CrossRef]
  30. Geng, G.; Li, Q.; Guo, X.; Ni, Q.; Xu, Y.; Ma, Z.; Wang, Y.; Ming, M. FOXO3a-modulated DEPDC1 promotes malignant progression of nephroblastoma via the Wnt/β-catenin signaling pathway. Mol. Med. Rep. 2022, 26, 272. [Google Scholar] [CrossRef]
  31. Gong, Z.; Chu, H.; Chen, J.; Jiang, L.; Gong, B.; Zhu, P.; Zhang, C.; Wang, Z.; Zhang, W.; Wang, J.; et al. DEPDC1 upregulation promotes cell proliferation and predicts poor prognosis in patients with gastric cancer. Cancer Biomark. 2021, 30, 299–307. [Google Scholar] [CrossRef] [PubMed]
  32. Zhang, L.; Li, Y.; Dai, Y.; Wang, D.; Wang, X.; Cao, Y.; Liu, W.; Tao, Z. Glycolysis-related gene expression profiling serves as a novel prognosis risk predictor for human hepatocellular carcinoma. Sci. Rep. 2021, 11, 1–21. [Google Scholar] [CrossRef] [PubMed]
  33. Tian, C.; Abudoureyimu, M.; Lin, X.; Chu, X.; Wang, R. Linc-ROR facilitates progression and angiogenesis of hepatocellular carcinoma by modulating DEPDC1 expression. Cell Death Dis. 2021, 12, 18875. [Google Scholar] [CrossRef] [PubMed]
  34. Zhou, C.; Wang, P.; Tu, M.; Huang, Y.; Xiong, F.; Wu, Y. DEPDC1 promotes cell proliferation and suppresses sensitivity to chemotherapy in human hepatocellular carcinoma. Biosci. Rep. 2019, 39, 1047. [Google Scholar] [CrossRef] [Green Version]
  35. Zhao, H.; Yu, M.; Sui, L.; Gong, B.; Zhou, B.; Chen, J.; Gong, Z.; Hao, C. High Expression of DEPDC1 Promotes Malignant Phenotypes of Breast Cancer Cells and Predicts Poor Prognosis in Patients with Breast Cancer. Front. Oncol. 2019, 9, 262. [Google Scholar] [CrossRef] [Green Version]
  36. Huang, L.; Chen, K.; Cai, Z.-P.; Chen, F.-C.; Shen, H.-Y.; Zhao, W.-H.; Yang, S.-J.; Chen, X.-B.; Tang, G.-X.; Lin, X. DEPDC1 promotes cell proliferation and tumor growth via activation of E2F signaling in prostate cancer. Biochem. Biophys. Res. Commun. 2017, 490, 707–712. [Google Scholar] [CrossRef]
  37. Feng, X.; Zhang, C.; Zhu, L.; Zhang, L.; Li, H.; He, L.; Mi, Y.; Wang, Y.; Zhu, J.; Bu, Y. DEPDC1 is required for cell cycle progression and motility in nasopharyngeal carcinoma. Oncotarget 2017, 8, 63605–63619. [Google Scholar] [CrossRef] [Green Version]
  38. Li, Z.; Liu, J.; Que, L.; Tang, X. The immunoregulatory protein B7-H3 promotes aerobic glycolysis in oral squamous carcinoma via PI3K/Akt/mTOR pathway. J. Cancer 2019, 10, 5770–5784. [Google Scholar] [CrossRef]
  39. Lin, C.-H.; Lee, H.-H.; Kuei, C.-H.; Lu, L.-S.; Chang, J.; Wang, J.-Y.; Hsu, K.-C. Nicotinic Acetylcholine Receptor Subunit Alpha-5 Promotes Radioresistance via Recruiting E2F Activity in Oral Squamous Cell Carcinoma. J. Clin. Med. 2019, 8, 1454. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Patil, S.; Al-Brakati, A.; Abidi, N.H.; Almasri, M.A.; Almeslet, A.S.; Patil, V.R.; Raj, A.T.; Bhandi, S. CD44-positive cancer stem cells from oral squamous cell carcinoma exhibit reduced proliferation and stemness gene expression upon adipogenic induction. Med. Oncol. 2022, 39, 23. [Google Scholar] [CrossRef] [PubMed]
  41. Lai, Y.-H.; Liu, H.; Chiang, W.-F.; Chen, T.W.; Chu, L.J.; Yu, J.-S.; Chen, S.-J.; Chen, H.-C.; Tan, B.C.-M. MiR-31-5p-ACOX1 Axis Enhances Tumorigenic Fitness in Oral Squamous Cell Carcinoma Via the Promigratory Prostaglandin E2. Theranostics 2018, 8, 486. [Google Scholar] [CrossRef]
  42. Zhao, W.; Cui, Y.; Liu, L.; Ma, X.; Qi, X.; Wang, Y.; Liu, Z.; Ma, S.; Liu, J.; Wu, J. METTL3 Facilitates Oral Squamous Cell Carcinoma Tumorigenesis by Enhancing c-Myc Stability via YTHDF1-Mediated m6A Modification. Mol. Ther. Nucleic Acids 2020, 20, 1–12. [Google Scholar] [CrossRef] [PubMed]
  43. Liang, X.; Deng, M.; Zhang, C.; Ping, F.; Wang, H.; Wang, Y.; Fan, Z.; Ren, X.; Tao, X.; Wu, T.; et al. Combined class I histone deacetylase and mTORC1/C2 inhibition suppresses the initiation and recurrence of oral squamous cell carcinomas by repressing SOX2. Cancer Lett. 2019, 454, 108–119. [Google Scholar] [CrossRef] [PubMed]
  44. The Cancer Genome Atlas Research Network; Weinstein, J.N.; Collisson, E.A.; Mills, G.B.; Shaw, K.R.M.; Ozenberger, B.A.; Ellrott, K.; Shmulevich, I.; Sander, C.; Stuart, J.M. The Cancer Genome Atlas Pan-Cancer analysis project. Nat. Genet. 2013, 45, 1113–1120. [Google Scholar] [CrossRef] [PubMed]
  45. Huang, G.; Lu, Z.; Rao, Y.; Gao, H.; Lv, X. Screening and identification of autophagy-related biomarkers for oral squamous cell carcinoma (OSCC) via integrated bioinformatics analysis. J. Cell. Mol. Med. 2021, 25, 4444–4454. [Google Scholar] [CrossRef]
  46. Huang, G.-Z.; Wu, Q.-Q.; Zheng, Z.-N.; Shao, T.-R.; Chen, Y.-C.; Zeng, W.-S.; Lv, X.-Z. M6A-related bioinformatics analysis reveals that HNRNPC facilitates progression of OSCC via EMT. Aging 2020, 12, 11667–11684. [Google Scholar] [CrossRef]
Figure 1. Gene set enrichment analysis (GSEA) screening of glycolysis-related gene sets. (A,B); Hallmark−glycolysis and reactome−glycolysis were significantly activated in OSCC. Higher hallmark−glycolysis and reactome−glycolysis activation were positively related to tumor samples. (C,D), The heatmap and volcano plot of differentially expressed core genes (DECGs) in these 2 gene sets. Blue and red stand for downregulation and upregulation, respectively.
Figure 1. Gene set enrichment analysis (GSEA) screening of glycolysis-related gene sets. (A,B); Hallmark−glycolysis and reactome−glycolysis were significantly activated in OSCC. Higher hallmark−glycolysis and reactome−glycolysis activation were positively related to tumor samples. (C,D), The heatmap and volcano plot of differentially expressed core genes (DECGs) in these 2 gene sets. Blue and red stand for downregulation and upregulation, respectively.
Ijms 24 01992 g001
Figure 2. Identification of prognostic Cox model. The expression level of 3 genes, DEPDC1 ((A), p < 0.0001), DDIT4 ((B), p = 0.0008) and SLC16A3 ((C), p < 0.0001), included in the risk Cox formula in TCGA. Triangle and square stand for gene expression level in normal and tumor samples, respectively. (D); Survival analysis between high- and low-risk levels according to risk scores calculated from the DEPDC1, DDIT4 and SLC16A3 expression profiles of each patient. (E,F), Univariate and multivariate Cox regression analyses on the basis of risk score. (G), High risk score might mean poor grade level in OSCC. Significant differences were considered at p < 0.001 ***; and p < 0.0001 ****.
Figure 2. Identification of prognostic Cox model. The expression level of 3 genes, DEPDC1 ((A), p < 0.0001), DDIT4 ((B), p = 0.0008) and SLC16A3 ((C), p < 0.0001), included in the risk Cox formula in TCGA. Triangle and square stand for gene expression level in normal and tumor samples, respectively. (D); Survival analysis between high- and low-risk levels according to risk scores calculated from the DEPDC1, DDIT4 and SLC16A3 expression profiles of each patient. (E,F), Univariate and multivariate Cox regression analyses on the basis of risk score. (G), High risk score might mean poor grade level in OSCC. Significant differences were considered at p < 0.001 ***; and p < 0.0001 ****.
Ijms 24 01992 g002
Figure 3. DEPDC1 was identified as a potential biomarker. (A), Survival analysis showed that high DEPDC1 expression levels meant poor overall survival. (B,C), DEPDC1 expression level was related to advanced stage and lymph node metastasis in OSCC. (D,E), GSEA suggested that high DEPDC1 expression was associated with higher activation of the WNT/β−catenin signaling pathway.
Figure 3. DEPDC1 was identified as a potential biomarker. (A), Survival analysis showed that high DEPDC1 expression levels meant poor overall survival. (B,C), DEPDC1 expression level was related to advanced stage and lymph node metastasis in OSCC. (D,E), GSEA suggested that high DEPDC1 expression was associated with higher activation of the WNT/β−catenin signaling pathway.
Ijms 24 01992 g003
Figure 4. Knockdown of DEPDC1 inhibited the migration and invasion of OSCC. (A); (B), The migration (SCC9 p = 0.0012, SCC15 p = 0.0010) and invasion abilities (SCC9 p = 0.0002, SCC15 p = 0.0005) of SCC9 and SCC15 were inhibited after knockdown of DEPDC1. (C); Knockdown of DEPDC1 resulted in a decrease in N-cad (SCC9 p = 0.0024, SCC15 p = 0.0022) and β-catenin (SCC9 p = 0.0007, SCC15 p < 0.0001), and an increase in E-cad (SCC9 p = 0.0102, SCC15 p = 0.0040). Error bars in the graphs represent S.D. N-cad: N-cadherin, E-cad: E-cadherin, siNC: negative control of siRNA, si-DEPDC1: knockdown of DEPDC1. Significant differences were considered at p < 0.05 *; p < 0.01 **; p < 0.001 ***; and p < 0.0001 ****.
Figure 4. Knockdown of DEPDC1 inhibited the migration and invasion of OSCC. (A); (B), The migration (SCC9 p = 0.0012, SCC15 p = 0.0010) and invasion abilities (SCC9 p = 0.0002, SCC15 p = 0.0005) of SCC9 and SCC15 were inhibited after knockdown of DEPDC1. (C); Knockdown of DEPDC1 resulted in a decrease in N-cad (SCC9 p = 0.0024, SCC15 p = 0.0022) and β-catenin (SCC9 p = 0.0007, SCC15 p < 0.0001), and an increase in E-cad (SCC9 p = 0.0102, SCC15 p = 0.0040). Error bars in the graphs represent S.D. N-cad: N-cadherin, E-cad: E-cadherin, siNC: negative control of siRNA, si-DEPDC1: knockdown of DEPDC1. Significant differences were considered at p < 0.05 *; p < 0.01 **; p < 0.001 ***; and p < 0.0001 ****.
Ijms 24 01992 g004
Figure 5. DEPDC1 promoted OSCC progression via the WNT/β-catenin signaling pathway. (A,B), Overexpression of DEPDC1 enhanced the migration and invasion of SCC9 and SCC15, while XAV939, an inhibitor of the WNT/β-catenin signaling pathway, could weaken the migration (SCC9 p = 0.0039; p = 0.0081, SCC15 p = 0.0024; p = 0.0134) and invasion ability (SCC9 p = 0.0005; p = 0.0138, SCC15 p = 0.0005; p = 0.0011) after overexpression of DEPDC1. (C), Meanwhile, upregulation of DEPDC1 induced higher N-cad (SCC9 p = 0.0024; p = 0.0042, SCC15 p = 0.0002; p = 0.0012) and β-catenin expression (SCC9 p = 0.0021; p = 0.0114, SCC15 p = 0.0014; p = 0.0075) and lower E-cad expression (SCC9 p = 0.0047; p = 0.0131, SCC15 p = 0.0064; p = 0.0004), while the protein level could also be reversed by XAV939. N-cad: N-cadherin, E-cad: E-cadherin, Con: control, OE: overexpression of DEPDC1. Significant differences were considered at p < 0.05 *; p < 0.01 **; p < 0.001 ***.
Figure 5. DEPDC1 promoted OSCC progression via the WNT/β-catenin signaling pathway. (A,B), Overexpression of DEPDC1 enhanced the migration and invasion of SCC9 and SCC15, while XAV939, an inhibitor of the WNT/β-catenin signaling pathway, could weaken the migration (SCC9 p = 0.0039; p = 0.0081, SCC15 p = 0.0024; p = 0.0134) and invasion ability (SCC9 p = 0.0005; p = 0.0138, SCC15 p = 0.0005; p = 0.0011) after overexpression of DEPDC1. (C), Meanwhile, upregulation of DEPDC1 induced higher N-cad (SCC9 p = 0.0024; p = 0.0042, SCC15 p = 0.0002; p = 0.0012) and β-catenin expression (SCC9 p = 0.0021; p = 0.0114, SCC15 p = 0.0014; p = 0.0075) and lower E-cad expression (SCC9 p = 0.0047; p = 0.0131, SCC15 p = 0.0064; p = 0.0004), while the protein level could also be reversed by XAV939. N-cad: N-cadherin, E-cad: E-cadherin, Con: control, OE: overexpression of DEPDC1. Significant differences were considered at p < 0.05 *; p < 0.01 **; p < 0.001 ***.
Ijms 24 01992 g005
Figure 6. DEPDC1 promoted the aerobic glycolysis of OSCC. (A), Knockdown of DEPDC1 showed lower accumulation of lactate (SCC9 p = 0.0251, SCC15 p = 0.0023) and glucose uptake (SCC9 p = 0.0121, SCC15 p = 0.0035) in OSCC. (B), Lactic acid production (SCC9 p = 0.0051; p < 0.0001, SCC15 p = 0.0002; p < 0.0001) and glucose uptake (SCC9 p = 0.0067; p = 0.4904, SCC15 p = 0.0005; p = 0.0023) induced by DEPDC1 could be inhibited by XAV939. (C,D), Glut1 protein levels after knockdown (SCC9 p = 0.0004, SCC15 p = 0.0140) or overexpression (SCC9 p = 0.0013; p = 0.0072, SCC15 p = 0.0007; p = 0.1741) of DEPDC1 in SCC9 and SCC15. Con: control, OE: overexpression of DEPDC1. Significant differences were considered at p < 0.05 *; p < 0.01 **; p < 0.001 ***; and p < 0.0001 ****.
Figure 6. DEPDC1 promoted the aerobic glycolysis of OSCC. (A), Knockdown of DEPDC1 showed lower accumulation of lactate (SCC9 p = 0.0251, SCC15 p = 0.0023) and glucose uptake (SCC9 p = 0.0121, SCC15 p = 0.0035) in OSCC. (B), Lactic acid production (SCC9 p = 0.0051; p < 0.0001, SCC15 p = 0.0002; p < 0.0001) and glucose uptake (SCC9 p = 0.0067; p = 0.4904, SCC15 p = 0.0005; p = 0.0023) induced by DEPDC1 could be inhibited by XAV939. (C,D), Glut1 protein levels after knockdown (SCC9 p = 0.0004, SCC15 p = 0.0140) or overexpression (SCC9 p = 0.0013; p = 0.0072, SCC15 p = 0.0007; p = 0.1741) of DEPDC1 in SCC9 and SCC15. Con: control, OE: overexpression of DEPDC1. Significant differences were considered at p < 0.05 *; p < 0.01 **; p < 0.001 ***; and p < 0.0001 ****.
Ijms 24 01992 g006
Table 1. The core genes screened out from 2 glycolysis-related pathways.
Table 1. The core genes screened out from 2 glycolysis-related pathways.
AAASB4GALT2CTHGALEHAX1MDH2NUP210PFKPPPP2R5DSLC37A4
ABCB6B4GALT4CXCR4GALK1HDLBPME1NUP214PGAM1PRKACASOD1
ADORA2BB4GALT7CYB5AGALK2HK1ME2NUP35PGAM2PRKACBSOX9
ADPGKBIKDCNGAPDHHK2MED24NUP37PGK1PRKACGSPAG4
AGLBPGMDDIT4GAPDHSHK3MERTKNUP43PGK2PRPS1SRD5A3
AGRNBPNT1DEPDC1GCKHMMRMETNUP50PGLSPSMC4STC1
AK3CACNA1HDLDGCKRHOMER1MIFNUP54PGM2PYGBSTC2
AK4CAPN5DPYSL4GCLCHS2ST1MIOXNUP58PGM2L1PYGLSTMN1
AKR1A1CASP6DSC2GFPT1HS6ST2MPINUP62PGPQSOX1TALDO1
ALDH7A1CD44ECDGLCEHSPA5MXI1NUP85PHKA2RAE1TFF3
ALDH9A1CDK1EFNA3GLRXIDH1NANPNUP88PKLRRANBP2TGFA
ALDOACENPAEGFRGMPPAIDUANASPNUP93PKMRBCK1TGFBI
ALDOBCHPFEGLN3GMPPBIER3NDC1NUP98PKP2RPETKTL1
ALDOCCHPF2ELF3GNEIGFBP3NDST3P4HA1PLOD1RRAGDTPBG
ALG1CHST1ENO1GNPDA1IL13RA1NDUFV3P4HA2PLOD2SAP30TPI1
ANGCHST12ENO2GNPDA2IRS2NOL3PAMPMM2SDC1TPR
ANGPTL4CHST2ENO3GOT1ISG20NSDHLPAXIP1POLR3KSDC2TPST1
ANKZF1CHST4ERO1AGOT2KDELR3NT5EPCPOM121SDC3TSTA3
ARPP19CHST6EXT1GPC1KIF20ANUP107PDK3POM121CSDHCTXN
ARTNCITED2EXT2GPC3KIF2ANUP133PFKFB1PPFIA4SEC13UGP2
AURKACLDN3FAM162AGPC4LCTNUP153PFKFB2PPIASEH1LVCAN
B3GALT6CLDN9FBP2GPILDHANUP155PFKFB3PPP2CASLC16A3VEGFA
B3GAT1CLN6FKBP4GPR87LDHCNUP160PFKFB4PPP2CBSLC25A10VLDLR
B3GAT3COG2FUT8GUSBLHPPNUP188PFKLPPP2R1ASLC25A13XYLT2
B3GNT3COL5A1G6PDGYS1LHX9NUP205PFKMPPP2R1BSLC35A3ZNF292
B4GALT1COPB2GAL3ST1GYS2MDH1
Table 2. The downregulated glycolysis-related genes in OSCC.
Table 2. The downregulated glycolysis-related genes in OSCC.
GenelogFCp-ValueGenelogFCp-Value
FBP2−5.0625<0.0001CAPN5−1.8001<0.0001
ENO3−4.31650.0003AGL−1.40330.0002
GYS2−3.2913<0.0001GCKR−1.3610.0003
CHST4−2.77610.0031GCKR−1.3610.0003
PFKFB1−2.53520.0002ELF3−1.2759<0.0001
PFKFB1−2.53520.0002GOT1−1.1533<0.0001
GAPDHS−2.45990.0003GNE−1.11370.0011
RRAGD−2.414<0.0001ALDH9A1−1.0858<0.0001
TFF3−2.3824<0.0001ANG−1.08210.001
CLDN3−2.02540.0124
Table 3. The upregulated glycolysis-related genes in OSCC.
Table 3. The upregulated glycolysis-related genes in OSCC.
GenelogFCp ValueGenelogFCp Value
NUP1551.0005<0.0001CDK11.6558<0.0001
XYLT21.022<0.0001EXT11.6601<0.0001
BPGM1.0342<0.0001P4HA21.69<0.0001
SPAG41.045<0.0001DEPDC11.7072<0.0001
GPC11.048<0.0001P4HA11.7279<0.0001
TGFA1.0601<0.0001TPBG1.7836<0.0001
NUP1071.0995<0.0001HK31.8105<0.0001
PFKFB41.1089<0.0001VCAN1.8404<0.0001
KDELR31.1278<0.0001AGRN1.8455<0.0001
ISG201.13420.0001KIF20A1.8507<0.0001
EXT21.1395<0.0001PLOD21.8644<0.0001
DDIT41.162<0.0001IGFBP31.9107<0.0001
B4GALT41.1713<0.0001HMMR1.9118<0.0001
VEGFA1.1938<0.0001NT5E2.1584<0.0001
CHST11.21330.0001AURKA2.2105<0.0001
SLC16A31.2878<0.0001CENPA2.2121<0.0001
CHPF1.2915<0.0001CHST22.4394<0.0001
STC11.36860.0001ENO22.6407<0.0001
MET1.3857<0.0001STC22.9567<0.0001
STMN11.4038<0.0001COL5A13.5266<0.0001
EGLN31.4516<0.0001ARTN3.6158<0.0001
PYGL1.4539<0.0001TGFBI4.0522<0.0001
EGFR1.4653<0.0001MIOX4.2703<0.0001
TPST11.5324<0.0001LHX94.6398<0.0001
PLOD11.5884<0.0001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Huang, G.; Chen, S.; Washio, J.; Paka Lubamba, G.; Takahashi, N.; Li, C. Glycolysis-Related Gene Analyses Indicate That DEPDC1 Promotes the Malignant Progression of Oral Squamous Cell Carcinoma via the WNT/β-Catenin Signaling Pathway. Int. J. Mol. Sci. 2023, 24, 1992. https://doi.org/10.3390/ijms24031992

AMA Style

Huang G, Chen S, Washio J, Paka Lubamba G, Takahashi N, Li C. Glycolysis-Related Gene Analyses Indicate That DEPDC1 Promotes the Malignant Progression of Oral Squamous Cell Carcinoma via the WNT/β-Catenin Signaling Pathway. International Journal of Molecular Sciences. 2023; 24(3):1992. https://doi.org/10.3390/ijms24031992

Chicago/Turabian Style

Huang, Guangzhao, Su Chen, Jumpei Washio, Grace Paka Lubamba, Nobuhiro Takahashi, and Chunjie Li. 2023. "Glycolysis-Related Gene Analyses Indicate That DEPDC1 Promotes the Malignant Progression of Oral Squamous Cell Carcinoma via the WNT/β-Catenin Signaling Pathway" International Journal of Molecular Sciences 24, no. 3: 1992. https://doi.org/10.3390/ijms24031992

APA Style

Huang, G., Chen, S., Washio, J., Paka Lubamba, G., Takahashi, N., & Li, C. (2023). Glycolysis-Related Gene Analyses Indicate That DEPDC1 Promotes the Malignant Progression of Oral Squamous Cell Carcinoma via the WNT/β-Catenin Signaling Pathway. International Journal of Molecular Sciences, 24(3), 1992. https://doi.org/10.3390/ijms24031992

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop