A four-hypoxia-genes-based prognostic signature for oral squamous cell carcinoma
BMC Oral Health volume 21, Article number: 232 (2021)
Oral squamous cell carcinoma (OSCC) is one of the most common maligancies of the head and neck. The prognosis was is significantly different among OSCC patients. This study aims to identify new biomarkers to establish a prognostic model to predict the survival of OSCC patients.
The mRNA expression and corresponding clinical information of OSCC patients were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus. Additionally, a total of 26 hypoxia-related genes were also obtained from a previous study. Univariate Cox regression analysis and LASSO Cox regression analysis were performed to screen the optimal hypoxia-related genes which were associated with the prognosis of OSCC. to establish the predictive model (Risk Score) was established for estimating the patient's overall survival (OS). Multivariate Cox regression analysis was used to determine whether the Risk Score was an independent prognostic factor. Based on all the independent prognostic factors, nomogram was established to predict the OS probability of OSCC patients. The relative proportion of 22 immune cell types in each patient was evaluated by CIBERSORT software.
We determined that a total of four hypoxia-related genes including ALDOA, P4HA1, PGK1 and VEGFA were significantly associated with the prognosis of OSCC patients. The nomogram established based on all the independent factors could reliably predict the long-term OS of OSCC patients. In addition, our resluts indicated that the inferior prognosis of OSCC patients with high Risk Score might be related to the immunosuppressive microenvironments.
This study shows that high expression of hypoxia-related genes including ALDOA, P4HA1, PGK1 and VEGFA is associated with poor prognosis in OSCC patients, and they can be used as potential markers for predicting prognosis in OSCC patients.
Oral squamous cell carcinoma (OSCC) is one of the most universal tumors in the head and neck, with the global incidence rate is the eighth of all cancers . Statistically, more than 300,000 new cases are diagnosed worldwide annually . Its pathogenic factors include smoking, drinking, chewing betel nut, carcinogen exposure, viral infection, immunodeficiency, gene specificity, etc. [2, 3]. At present, there were certain advances in surgery, radiotherapy and chemotherapy for OSCC, but the fatality rate remains high . And the survival of OSCC patients is significantly different. Therefore, it is of great significance to screen biomarkers related to prognosis of OSCC.
Many original tumor areas have low concentrations of molecular oxygen, known as hypoxia . Hypoxia is the specific microenvironment of tumor, which is closely associated with the initiation and progression of cancer. The causes include disordered tumor blood vessels leading to insufficient oxygen supply and accelerated tumor metabolism leading to increased oxygen consumption. Hypoxia increases the migration of tumor cells. Hypoxia-inducible factors (HIFs) induce the expression of matrix metalloproteinases and other protease genes, leading to the degradation of the matrix around the tumor and providing channels for tumor metastasis . Hypoxia can change the abundance of gene expression, which will result in proteomics changes and affect cell physiological and biochemical functions . For example, the abundance of DNA double-stranded repair protein RAD51 is down-regulated in anoxic cells, resulting in inhibition of DNA repair. In addition, DNA mismatch repair genes MLH1 and MSH2 are down-regulated under hypoxia, leading to a large part of gene mutations [8, 9]. Furthermore, hypoxia affects the sensitivity of tumor cells to chemotherapy and radiotherapy, and its presence has a significant negative impact on the prognosis of patients after radiotherapy .
In recent years, the significance of hypoxia in the treatment of cancer has been identified. The earliest research indicated that hypoxia had effect on the proliferation of mouse embryonic fibroblasts, and the inactivation of HIF-1 enhances the inhibition of cell proliferation . Further, overexpression of HIF-1 was found in human tumor tissues commonly, including breast, stomach, skin, ovary, and pancreas . For instance, HIF expression level in pancreatic cancer is significantly increased and can be used as a clinical marker of pancreatic cancer survival rate to assess the prognosis of patients. Platinum insensitivity often occurs during platinum chemotherapy for non-small cell lung cancer (NSCLC), but there is evidence that A549 NSCLC cells are significantly more sensitive to platinum after reoxygenation . However, the prognostic value of hypoxia-related genes is very limited presently, especially in OSCC.
Therefore, in this study, the prognostic value of 26 hypoxia related genes in OSCC was explored, four optimal genes associated with the prognosis of OSCC were identified, and a reliable prognostic model was established, in order to provide some guidance for the clinical application of OSCC.
Material and methods
The information of 349 OSCC patients about genome-wide expressions and corresponding clinical data were downloaded from The Cancer Genome Atlas (TCGA, https://tcga-data.nci.nih.gov/tcga/), among which 348 patients had complete survival data, which were used for subsequent analysis. In addition, we also obtained mRNA expression and clinical information of 218 OSCC patients with the accession number of GSE65858 from Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/) . Clinical information for all patients in the TCGA and GEO datasets is shown in Table1.
Here the expression value of 26 hypoxia genes previously reported to associate with the efficacy of hypoxia therapy for laryngeal cancer  was used to clustered samples, and the 26 hypoxia related genes were shown in Additional file 1: Table S1.
Based on the mRNA expression of the 26 hypoxia-related genes, the function package factoextra (https://cran.cran.project.org/package=factoextra) in R language was used to perform sample clustering analysis.
LASSO Cox regression analysis
Basing on the mRNA expression values of 26 hypoxia-related genes, the univariate Cox regression analysis was conducted, and the genes significantly related to the prognosis of OSCC were selected with P < 0.05 as the threshold. Then, glmnet package of R language  was used for LASSO Cox regression analysis to further refine the hypoxia-related genes related to the prognosis of OSCC, and the Risk Score of each sample was calculated according to the refined genes:
In the formula, Coefi is the risk coefficient. Xi is the expression value of hypoxia-related genes. The samples were divided into Low Risk group and High Risk group based on the median of Risk Score.
The overall survival rate of different groups was estimated based on the Kaplan–Meier method using the survial package and the survminer package of R language, and log-rank test was used to determine the significance of the difference in survival rate between different groups with P < 0.05 as a significance threshold. Finally, multivariate Cox regression analysis was used to determine whether Risk Score was an independent factor to predict the survival of patients with OSCC.
Analysis of immune cell infiltration
The CIBERSORT software  was used to compute the relative ratios of 22 immune cell types in the samples. Based on the gene expression matrix, the deconvolution algorithm used the preset 547 barcode genes to characterize the proportion of the infiltrating immune cells. The sum of the estimated proportions of the 22 immune cells in each sample was 1.
Immune checkpoints among different groups
The corAndPvalue function in the WGCNA package  of R software was used to calculate the correlation between the immune checkpoint and Risk Score. The chordDiagram function in circlize package  was conducted for visualization of correlation.
Construction and evaluation of the nomogram
Based on the independent prognostic factors identified by multivariate Cox regression analysis with RMS (https://cran.cran.r-project.org/package=RMS) package, nomogram was established to predict overall survival probability for 1-year, 3-year and 5-year. A calibration curve was drawn to determine the divergency between Nomogram's predicted probability and actual incidence.
Wilcoxon signed-rank test was performed to compare the expression of immune checkpoints in different groups. P < 0.05 was used as the significance threshold. All statistical analyses were carried out by using R software V3.5.2.
Hypoxia-related genes could separate prognostically different OSCC patients
The factoextra function package in R language was used for clustering analysis for OSCC patients in TCGA based on the expression values of 26 hypoxia-related genes. According to the sum of the squared errors (SSE), the optimal number of clusters was determined as 3 (Fig. 1a). The clustering diagram (Fig. 1b) clearly showed that the OSCC patients from TCGA could classified be basically separated into three clusters, including 135, 159 and 54 cancer samples, respectively. The expression heatmap (Fig. 1c) indicated that the expression levels of Cluster 1 to Cluster 3 were gradually decreased, and the three clusters could be obviously identified. The above result was also supported by principal component analysis (PCA) (Fig. 1d), which revealed the differences of the three clusters of samples were obvious and the samples of the same Cluster were more closely similar. Kaplan–Meier survival analysis showed that there were obvious differences in the overall survival (OS) among the three OSCC groups (P = 0.008) (Fig. 1e), and patients in Cluster 1 had poorer survival. It indicated that hypoxia-related genes could discriminate the prognostically different OSCC samples effectively.
Hypoxia related genes-based prognostic model for OSCC
The expression values of 26 hypoxia-related genes in TCGA samples were used as continuous variables for univariate cox regression analysis to explore their associations with OSCC patients’ OS probability. There were 12 genes (ALDOA, ANLN, CA9, CDKN3, ENO1, HILPDA, MRGBP, P4HA1, PGAM1, PGK1, TPI1 and VEGFA) determined to be significantly associated with the OS of OSCC (Fig. 2a), which could be all risk genes for OSCC patients (HR > 1, P value < 0.05), indicating that high expression of these genes may cause poor prognosis.
Then the 12 hypoxia-related genes were used for LASSO Cox regression analysis. Based on the lambda value corresponding to the different number of genes, 4 optimal genes (ALDOA, P4HA1, PGK1, and VEGFA) were selected (the lambda value is the smallest) (Fig. 2b, c). Further, we established a predictive risk score model for OS of OSCC patients based on the regression coefficients weighted expressions of the four hypoxia-related genes (Risk Score = 0.2777 * Express Value of ALDOA + 0.1355 * Express Value of P4HA1 + 0.1237 * Express Value of PGK1 + 0.047 * Express Value of VEGFA). The OSCC patients in TCGA set and GEO set were all assigned risk scores based on the model and respectively divided into low- and high-risk groups according to the median of the Risk Score (Fig. 2d, g). It was found that the expression values of hypoxia-related genes in the model were significantly different between the high and low risk groups (Fig. 2e, h). Additionally, the OS of OSCC samples in the high-risk group was poorer than those in the low-risk group in both TCGA and GEO dataset (P < 0.0001 for TCGA, and P = 0.032 for GEO) (Fig. 2f, i). Overall, it showed that the prognostic model constructed using the four hypoxia-related genes, including ALDOA, P4HA1, PGK1, and VEGFA, could well predict the prognosis of OSCC patients.
Risk Score was significantly associated with OSCC patients’ immune microenvironment
In this study, CIBERSORT method was used to estimate the relative infiltrating proportions of 22 immune cells in OSCC patients from TCGA. The infiltrating landscape of the 22 immune cells exhibited clear discrepancies among OSCC patients (Fig. 3a). In addition, the correlations between the infiltrating ratios of most different types of immune cells is weak (Fig. 3b), indicating that there was great heterogeneity in the infiltration of different immune cells in OSCC patients. Additionally, it was found that there are 7 types of immune cells that have significant differences in the degree of infiltration between OSCC patients within high- and low-risk groups, including B cells naive, T cells CD8, T cells follicular helper, T cells regulatory (Tregs), Macrophages M0, Mast cells activated and Neutrophils. Moreover, compared to infiltration rates in low-risk group, only Macrophages M0 and Mast cells were significantly higher in the high-risk group, while the other five types of immune cells were markedly lower (Fig. 3c), and difference in the amount of immune cells might be the case for prognostic differences of OSSC patients.
As previously reported, the expression of immune checkpoints has become a biomarker for the selection of immunotherapy . Here we analyzed the correlation between the six immune checkpoints CTLA4, PDL1, PDL2, TIM3, LAG3, TIGIT and the risk score of TCGA samples, and found that the expression of PDL2 and TIM3 was significantly correlated with the Risk Score (Fig. 3d and Additional file 1: Table S2). We also studied the expression of these 6 immune checkpoints in OSCC patients in high and low risk groups, and found that the expressions of PDL2 and TIM3 were significantly different in OSCC patients in high- and low-risk groups, and the high-risk group has higher expression level (P value < 0.05) (Fig. 3e), indicating that the high-risk patients may more sensitive to the treatment of PDL2 and TIM3 immune checkpoint inhibitor.
Risk Score was an independent prognostic factor for OSCC
We conducted multivariate Cox regression analysis, using the characteristic related to OSCC progression, including Age, Gender, TNM Stage, Grade, Lymphovascular invasion, Perineural invasion and Risk Score as variables to analyze whether Risk Score was an independent prognostic factor. It was found that the risk score is still significantly correlated with the OS, and the greater the risk of death with a higher Risk Score, further indicates that Risk Score was a factor for poor prognosis (HR = 3.437, 95%CI 1.735–6.81, P < 0.001). In addition, Age and Perineural invasion were also independent prognostic factors (Fig. 4a).
We regrouped OSCC patients by Gender and Grade and performed Kaplan–Meier survival analysis. Among the female samples (Fig. 4b), male samples (Fig. 4c), early grade (G1 + G2) samples (Fig. 4d), and late grade (G2 + G3) samples (Fig. 4e), the OS of the high-risk group was significantly worse than that in the low-risk group. These results suggested that the Risk Score can be used as an independent factor to predict the prognosis of OSCC patients.
Nomogram model can predict the long-term OS of OSCC patients
The Nomogram model was constructed based on the three independent prognostic factors, including Age, Perineural invasion and Risk Score (Fig. 5a). For each OSCC sample of TCGA, draw vertical lines upwards to determine the points obtained from each factor in the Nomogram. And the sum of these points is located on the "Total Points" axis, and then a line is drawn down from the Total Points axis to obtain the predicted OS probability of OSCC patients for 1, 3, and 5 years. And the calibration curve was close to the ideal curve in the calibration chart (grey straight line), indicating that the predicted OS probability was consistent with the actual (Fig. 5b–d). The Nomogram could effectively predict the long-term survival probability of OSCC.
OSSC originates from the oral cavity and oropharyngeal mucosa, representing one of the most common malignancies in the head and neck. In addition to causing local hyperplasia and surrounding tissue erosion, it often metastasizes to lymph nodes, seriously affecting the health and lives of patients . we identified 4 hypoxia-related genes, which may related to the prognosis of OSCC patients.
Hypoxia is a key factor in tumor pathophysiology, causing changes in cell metabolism, triggering various molecular reactions, and promoting tumor spread. Hypoxia-inducible factors (HIFs) are the main molecules related to hypoxia research. It is known that HIFs regulate a large number of genes involved in cell proliferation, movement, metabolism, and angiogenesis through inducing the expression of its downstream target genes . The process of epithelial to mesenchymal transition (EMT) is related to cancer metastasis. According to reports, hypoxia can cause EMT in many cancers, including breast, prostate and oral cancer . Therefore, we speculated that hypoxia-related genes might have prognostic values in OSCC patients.
In this study, four hypoxia-related genes, ALDOA, P4HA1, PGK1 and VEGFA, were identified by LASSO Cox regression analysis, which were significantly associated with the overall survival rate of OSCC. Moreover, they were all risk genes for OSCC death event, i.e. their high expression would lead to poor prognosis. Human fructose diphosphate aldolase (ALDOA) and phosphoglycerate kinase 1 (PGK1) are the key rate-limiting enzymes in the process of aerobic glycolysis [23, 24]. Tumor cells tend to produce ATP through aerobic glycolysis rather than oxidative phosphorylation (Warburg Effect). In the hypoxic environment, Warburg Effect is induced by HIF-1 . Therefore, key enzymes of glycolysis in the hypoxic microenvironment can be used as signal molecules to participate in the regulation of tumor-related signal pathways and enhance the migration or proliferation of tumors. Previously, studies have shown that ALDOA was highly expressed in colon cancer, pancreatic cancer and osteosarcoma [26,27,28], and was associated with their poor prognosis. PGK1 is directly regulated by HIF-1 in many tumors. Early studies have shown that PGK1 was activated by HIF-1 in colorectal cancer and liver cancer under hypoxia stress . Due to HIF-1/PGK1-mediated epithelial-mesenchymal transformation (EMT), PGK1 is associated with the ability of cancer cells to metastasize , and the high expression of PGK1 is related to the inferior survival outcome of breast cancer patients . Therefore, we speculated that ALDOA and PGK1 as key enzymes in the glycolysis process may be regulated by HIF-1, and the high expression of ALDOA and PGK1 in the hypoxic environment could enhance the proliferation and migration of tumors. Indeed, several researches has reported that ALDOA and PGK1 also mediated in OSCC and head and neck squamous cells carcinoma (HNSCC) [30, 32].Ahluwalia also reported that HIF-1 under hypoxic conditions regulates the expression of vascular endothelial growth factor A (VEGFA) at the transcriptional level . The angiogenesis simulation (VM)  proposed by Maniotis in 1999 was a novel cancer marker, which has prognostic value in colorectal cancer, breast cancer, melanoma, etc. [35,36,37]. High VM expression was also considered to be a risk factor for poor prognosis, low survival rate, invasion and metastasis in patients with cancer, and HIF-1 expression was associated with VM . In addition, the angiogenesis caused by hypoxia is usually VEGFA-dependent. This suggests that HIF regulation of VEGFA under hypoxic conditions may lead to high VM expression and a high risk in OSSC patients. Prolyl 4-hydroxylase subunit α1 (P4HA1) was encoding the active catalytic component of prolyl 4-hydroxylase (P4H), the key enzyme for collagen production . It was reported that HIF-1 promotes P4HA1 expression to induce collagen deposition and a more aggressive phenotype in human breast cancer, leading to poor prognosis for breast cancer patients . Perhaps the high level of P4HA1 expression decreases the OS probability of OSSC patients through the same process. These studies proved that our results were reliable and the over-expression of 4 genes might involve in OSCC progression and prognosis.
Immune infiltration has always been an important indicator of tumor prognosis . Hypoxia was a prominent characteristic of chronically inflamed tissues, which could affect disease progression by inflammatory signaling pathways in immune and non-immune cells . Thus, this study analyzed the differences in immune infiltration of 22 types of immune cells between OSCC patients in high and low risk groups, and found 7 types of immune cells: B cells naive, T cells CD8, T cells follicular helper, T cells regulatory (Tregs), Macrophages M0, Mast Cell activated and Neutrophils have significant differences in the degree of infiltration between high and low risk groups, and the infiltration proportion of immune cell in the high-risk group is higher, which indicated that immune infiltration is related to the prognosis of OSCC. The immune checkpoint expression has become a biomarker for cancer patients in the selection of immunotherapy. The study also found that the expression of the two immune checkpoints PDL2 and TIM3 were significantly different in the high and low risk groups. The expression in the high-risk group was significantly higher, indicating that patients in the high-risk group may more sensitive to the treatment of PDL2 and TIM3 inhibitor. Recently, combined targeting of TIM-3 and PD-1/PDL-1 has been proposed in hematological malignancies. It can be synergistically administered to control tumor growth . It suggested that the poor prognosis of OSCC caused by hypoxia may be related to the immune microenvironment.
Overall, the prognostic model we constructed based on four optimal hypoxia-related genes that could effectively predict the survival probability of OSCC patients. However, the specific functions and potential mechanism of the above hypoxia-related genes in OSSC are required to be further explored.
Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM. Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Cancer. 2010;127(12):2893–917.
Argiris A, Karamouzis MV, Raben D, Ferris RL. Head and neck cancer. Lancet. 2008;371(9625):1695–709.
Hashibe M, Brennan P, Chuang SC, Boccia S, Castellsague X, Chen C, et al. Interaction between tobacco and alcohol use and the risk of head and neck cancer: pooled analysis in the International Head and Neck Cancer Epidemiology Consortium. Cancer Epidemiol Biomarkers Prev. 2009;18(2):541–50.
Gupta S, Kong W, Peng Y, Miao Q, Mackillop WJ. Temporal trends in the incidence and survival of cancers of the upper aerodigestive tract in Ontario and the United States. Int J Cancer. 2009;125(9):2159–65.
Bhandari V, Hoey C, Liu LY, Lalonde E, Ray J, Livingstone J, et al. Molecular landmarks of tumor hypoxia across cancer types. Nat Genet. 2019;51(2):308–18.
Sakamoto T, Seiki M. Integrated functions of membrane-type 1 matrix metalloproteinase in regulating cancer malignancy: Beyond a proteinase. Cancer Sci. 2017;108(6):1095–100.
Jing X, Yang F, Shao C, Wei K, Xie M, Shen H, et al. Role of hypoxia in cancer therapy by regulating the tumor microenvironment. Mol Cancer. 2019;18(1):157.
Nardinocchi L, Puca R, Sacchi A, D’Orazi G. Inhibition of HIF-1alpha activity by homeodomain-interacting protein kinase-2 correlates with sensitization of chemoresistant cells to undergo apoptosis. Mol Cancer. 2009;8:1.
Sullivan R, Pare GC, Frederiksen LJ, Semenza GL, Graham CH. Hypoxia-induced resistance to anticancer drugs is associated with decreased senescence and requires hypoxia-inducible factor-1 activity. Mol Cancer Ther. 2008;7(7):1961–73.
Horsman MR, Overgaard J. The impact of hypoxia and its modification of the outcome of radiotherapy. J Radiat Res. 2016;57(Suppl 1):i90–8.
Daponte A, Ioannou M, Mylonis I, Simos G, Minas M, Messinis IE, et al. Prognostic significance of Hypoxia-Inducible Factor 1 alpha(HIF-1 alpha) expression in serous ovarian cancer: an immunohistochemical study. BMC Cancer. 2008;8:335.
Deschoemaeker S, Di Conza G, Lilla S, Martin-Perez R, Mennerich D, Boon L, et al. PHD1 regulates p53-mediated colorectal cancer chemoresistance. EMBO Mol Med. 2015;7(10):1350–65.
Wichmann G, Rosolowski M, Krohn K, Kreuz M, Boehm A, Reiche A, et al. The role of HPV RNA transcription, immune response-related gene expression and disruptive TP53 mutations in diagnostic and prognostic profiling of head and neck cancer. Int J Cancer. 2015;137(12):2846–57.
Eustace A, Mani N, Span PN, Irlam JJ, Taylor J, Betts GN, et al. A 26-gene hypoxia signature predicts benefit from hypoxia-modifying therapy in laryngeal cancer but not bladder cancer. Clin Cancer Res. 2013;19(17):4879–88.
Friedman J, Hastie T, Tibshirani R. Regularization Paths for Generalized Linear Models via Coordinate Descent. J Stat Softw. 2010;33(1):1–22.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Zhang B, Horvath S. A general framework for weighted gene co-expression network analysis. Stat Appl Genet Mol Biol. 2005;4:Article17.
Gu Z, Gu L, Eils R, Schlesner M, Brors B. circlize Implements and enhances circular visualization in R. Bioinformatics. 2014;30(19):2811–2.
Sharma P, Wagner K, Wolchok JD, Allison JP. Novel cancer immunotherapy agents with survival benefit: recent successes and next steps. Nat Rev Cancer. 2011;11(11):805–12.
Liu Q, Zhang X, Zeng Y. Targeted and personalized therapy for cancer: theory and practice in China. Sci China Life Sci. 2011;54(12):1081–4.
Bertout JA, Patel SA, Simon MC. The impact of O2 availability on human cancer. Nat Rev Cancer. 2008;8(12):967–75.
Joseph JP, Harishankar MK, Pillai AA, Devi A. Hypoxia induced EMT: a review on the mechanism of tumor progression and metastasis in OSCC. Oral Oncol. 2018;80:23–32.
Telford DJ, Stewart BW. Characteristics of chromatin release during digestion of nuclei with micrococcal nuclease: preferential solubilization of nascent RNA at low enzyme concentration. Int J Biochem. 1989;21(11):1235–40.
McCarrey JR, Thomas K. Human testis-specific PGK gene lacks introns and possesses characteristics of a processed gene. Nature. 1987;326(6112):501–5.
Guney E, Menche J, Vidal M, Barabasi AL. Network-based in silico drug efficacy screening. Nat Commun. 2016;7:10331.
Cancer Genome Atlas Research Network. Electronic address wbe, Cancer Genome Atlas Research N. Comprehensive and integrative genomic characterization of hepatocellular carcinoma. Cell. 2017;169(7):1327–41 e23.
Bjornson E, Mukhopadhyay B, Asplund A, Pristovsek N, Cinar R, Romeo S, et al. Stratification of hepatocellular carcinoma patients based on acetate utilization. Cell Rep. 2015;13(9):2014–26.
Lee JS, Chu IS, Heo J, Calvisi DF, Sun Z, Roskams T, et al. Classification and prediction of survival in hepatocellular carcinoma by gene expression profiling. Hepatology. 2004;40(3):667–76.
Kress S, Stein A, Maurer P, Weber B, Reichert J, Buchmann A, et al. Expression of hypoxia-inducible genes in tumor cells. J Cancer Res Clin Oncol. 1998;124(6):315–20.
He Y, Luo Y, Zhang D, Wang X, Zhang P, Li H, et al. PGK1-mediated cancer progression and drug resistance. Am J Cancer Res. 2019;9(11):2280–302.
Fu D, He C, Wei J, Zhang Z, Luo Y, Tan H, et al. PGK1 is a potential survival biomarker and invasion promoter by regulating the HIF-1alpha-mediated epithelial-mesenchymal transition process in breast cancer. Cell Physiol Biochem. 2018;51(5):2434–44.
Lessa RC, Campos AH, Freitas CE, Silva FR, Kowalski LP, Carvalho AL, et al. Identification of upregulated genes in oral squamous cell carcinomas. Head Neck. 2013;35(10):1475–81.
Ahluwalia A, Tarnawski AS. Critical role of hypoxia sensor–HIF-1alpha in VEGF gene activation. Implications for angiogenesis and tissue injury healing. Curr Med Chem. 2012;19(1):90–7.
Maniotis AJ, Folberg R, Hess A, Seftor EA, Gardner LM, Pe’er J, et al. Vascular channel formation by human melanoma cells in vivo and in vitro: vasculogenic mimicry. Am J Pathol. 1999;155(3):739–52.
Zhang S, Fu Z, Wei J, Guo J, Liu M, Du K. Peroxiredoxin 2 is involved in vasculogenic mimicry formation by targeting VEGFR2 activation in colorectal cancer. Med Oncol. 2015;32(1):414.
Ruffini F, Graziani G, Levati L, Tentori L, D’Atri S, Lacal PM. Cilengitide downmodulates invasiveness and vasculogenic mimicry of neuropilin 1 expressing melanoma cells through the inhibition of alphavbeta5 integrin. Int J Cancer. 2015;136(6):E545–58.
Tu DG, Yu Y, Lee CH, Kuo YL, Lu YC, Tu CW, et al. Hinokitiol inhibits vasculogenic mimicry activity of breast cancer stem/progenitor cells through proteasome-mediated degradation of epidermal growth factor receptor. Oncol Lett. 2016;11(4):2934–40.
Li W, Zong S, Shi Q, Li H, Xu J, Hou F. Hypoxia-induced vasculogenic mimicry formation in human colorectal cancer cells: Involvement of HIF-1a, Claudin-4, and E-cadherin and Vimentin. Sci Rep. 2016;6:37534.
Myllyharju J, Kivirikko KI. Collagens, modifying enzymes and their mutations in humans, flies and worms. Trends Genet. 2004;20(1):33–43.
Gilkes DM, Chaturvedi P, Bajpai S, Wong CC, Wei H, Pitcairn S, et al. Collagen prolyl hydroxylases are essential for breast cancer metastasis. Cancer Res. 2013;73(11):3285–96.
Cummins EP, Keogh CE, Crean D, Taylor CT. The role of HIF in immunity and inflammation. Mol Aspects Med. 2016;47–48:24–34.
Koyama S, Akbay EA, Li YY, Herter-Sprie GS, Buczkowski KA, Richards WG, et al. Adaptive resistance to therapeutic PD-1 blockade is associated with upregulation of alternative immune checkpoints. Nat Commun. 2016;7:10501.
This study was supported by the grants from Tianjin Stomatology Hospital, School of Medicine, NanKai University. Tianjin Key Laboratory of Oral and Maxillofacial Function Reconstruction (Grant No. HWZX017).
Ethics approval and consent to participate
Consent for publication
Conflict of interest
The authors declare that there is no conflict of interest regarding the publication of this article.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Zhao, C., Zhou, Y., Ma, H. et al. A four-hypoxia-genes-based prognostic signature for oral squamous cell carcinoma. BMC Oral Health 21, 232 (2021). https://doi.org/10.1186/s12903-021-01587-z