- Research article
- Open Access
Construction of an 11-microRNA-based signature and a prognostic nomogram to predict the overall survival of head and neck squamous cell carcinoma patients
BMC Genomics volume 21, Article number: 691 (2020)
Head and neck squamous cell carcinoma (HNSCC) is a fatal malignancy owing to the lack of effective tools to predict overall survival (OS). MicroRNAs (miRNAs) play an important role in HNSCC occurrence, development, invasion and metastasis, significantly affecting the OS of patients. Thus, the construction of miRNA-based risk signatures and nomograms is desirable to predict the OS of patients with HNSCC. Accordingly, in the present study, miRNA sequencing data of 71 HNSCC and 13 normal samples downloaded from The Cancer Genome Atlas (TCGA) were screened to identify differentially expressed miRNAs (DEMs) between HNSCC patients and normal controls. Based on the exclusion criteria, the clinical information and miRNA sequencing data of 67 HNSCC samples were selected and used to establish a miRNA-based signature and a prognostic nomogram. Forty-three HNSCC samples were assigned to an internal validation cohort for verifying the credibility and accuracy of the primary cohort. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to explore the functions of 11 miRNA target genes.
In total, 11 DEMs were successfully identified. An 11-miRNA risk signature and a prognostic nomogram were constructed based on the expression levels of these 11 DEMs and clinical information. The signature and nomogram were further validated by calculating the C-index, area under the curve (AUC) in receiver-operating characteristic curve analysis, and calibration curves, which revealed their promising performance. The results of the internal validation cohort shown the reliable predictive accuracy both of the miRNA-based signature and the prognostic nomogram. GO and KEGG analyses revealed that a mass of signal pathways participated in HNSCC proliferation and metastasis.
Overall, we constructed an 11-miRNA-based signature and a prognostic nomogram with excellent accuracy for predicting the OS of patients with HNSCC.
Head and neck squamous cell carcinoma (HNSCC), the sixth most common and eighth most fatal malignancy worldwide  is an epithelial tumor arising from the oral cavity, nasal cavity, larynx, hypopharynx, and pharynx. Excessive consumption of tobacco and alcohol is considered a major risk factor for the occurrence and development of HNSCC . In addition, human papillomavirus infection was recently confirmed as an important factor underlying HNSCC progression [2, 3]. Despite the rapid development in examination techniques, HNSCC is generally detected at advanced stages owing to the lack of awareness of regular inspections and no or mild symptoms at early stages. Hence, HNSCC is associated with high mortality . Many patients with HNSCC develop distant metastases within 5 years of receiving comprehensive and systematic chemotherapy . This serves as a significant contributor to death. Thus, improvement in the screening rate of early tumors may be useful as an effective measure to reduce HNSCC-related mortality.
MicroRNAs (miRNAs) are short nonprotein-coding RNAs involved in post-transcriptional regulation of protein-coding gene expression via binding to the 3′-untranslated regions of target mRNAs . miRNAs participate in various physiological and pathological activities in the human body, including cell development, differentiation, cycle regulation, and apoptosis [6, 7]. Several studies have reported the potential diagnostic or prognostic roles of miRNAs in HNSCC, including those of miR-let-7a-5p, miR-3928, miR-936, miR-383, miR-615, miR-877, miR-9-5p, and miR-9-3p [8,9,10,11]. The suppression of miR-30a and miR-379 expression could facilitate the oncogenic activity via upregulating the DNMT3B expression and activating the hypermethylation of ADHFE1 and ALDH1A genes in oral squamous cell carcinoma . In addition, circ-0000495 has been shown to sponge miR-488-3p expression and epigenetically silence TROP2 expression, resulting in the weakening of the proliferative capacity of HNSCC . Thus, the functions of miRNAs affect HNSCC generation, development, and metastasis and are highly associated with the overall survival (OS) of patients with HNSCC.
In the present study, we investigated the miRNAs that were closely bound up with the OS of patients with HNSCC. A miRNA-based signature based on differentially expressed miRNAs (DEMs) as well as a novel miRNA-based prognostic model were constructed to reliably predict the OS of HNSCC patients and provide an important tool for clinicians to improve treatment regimens.
Identification of DEMs associated with HNSCC patients
Raw HNSCC datasets, consisting of 71 HNSCC samples and 13 normal samples, were downloaded from The Cancer Genome Atlas (TCGA) database. In total, 797 miRNAs were acquired after eliminating those with expression levels < 1. In the heatmap (Fig. 1), the expression levels of 50 miRNAs were visually displayed. The differential expression of 797 miRNAs was visually observed using a volcano plot (Fig. 1). Of these, 90 miRNAs with |log2FC| ≥ 2 and an adjusted P-value < 0.05, including 54 upregulated and 36 downregulated miRNAs, showed significant differential expression. After eliminating the miRNAs detected in 13 normal samples and 4 patients, 90 DEMs were subjected to a univariate Cox proportional hazard regression (CPHR) analysis to determine the independent prognostic impact of individual genes. The results of the univariate CPHR analysis showed that 16 DEMs had the capacity to influence prognosis. Next, these 16 DEMs were subjected to LASSO Cox analysis, and a LASSO Cox regression model with a 10-fold cross validation result was proposed (Fig. 1c and d). In total, 11 DEMs were identified the close correlation with the prognosis of patients with HNSCC (Table 1).
Construction of a risk signature
The 11 DEMs verified from the LASSO regression analysis were used to generate a risk signature as per the following formula: Risk score = (0.236 × expression miR-204-5p) + (0.059 × expression miR-499a-5p) + (0.212 × expression miR-498-5p) − (0.062 × expression miR-155-3p) + (0.434 × expression miR-4714-3p) − (0.141 × expression miR-365a-5p) + (0.321 × expression miR-30a-5p) + (0.123 × expression miR-1-5p) + (0.240 × expression miR-548f-3p) + (0.196 × expression miR-518a-3p) − (0.140 × expression miR-196b-5p). Patients with HNSCC were distributed into high-risk and low-risk groups according to the median of risk score value. The new heatmap generated (Fig. 2a) clearly revealed the differences in the expression levels of the 13 DEMs between high-risk and low-risk groups. Eight DEMs (miR-204-5p, miR-499a-5p, miR-498-5p, miR-4714-3p, miR-30a-5p, miR-1-5p, miR-548f-3p, and miR-518a-3p) in the primary and internal validation cohorts showed higher expression in the high-risk group than that in the low-risk group. Contrarily, miR-155-3p, miR-365a-5p, and miR-196b-5p were overexpressed in the low-risk group, suggesting that they might might function as tumor suppressors. The survival status and risk score distribution analyses further demonstrated the high risk in the high-risk group (Fig. 2b and c). We established a prognostic nomogram associated with 11 DEMs (Fig. 2d) and found that miR-4714-3p, miR-30a-5p, and miR-548f-3p strongly affected the OS of patients.
Estimation of the reliability of the risk signature
To estimate the reliability of the risk signature established herein, a Kaplan–Meier survival analysis (Fig. 3a) was performed. The result of this analysis revealed the shorter OS for patients from the high-risk group than for those from the low-risk group both in the primary (P = 5.393e− 06) and internal validation cohorts (P = 5.176e− 04). In addition, the area under the curve (AUC) value of the risk signature for 5-year OS had reliable predictive accuracy (Fig. 3b). In the primary cohort, the AUC values of the receiver-operating characteristic (ROC) curve analysis for the risk signature for 1-, 3-, and 5-year OS were 0.802, 0.804, and 0.825, respectively. These values were reported to be 0.724, 0.811, and 0.829 for 1-, 3-, and 5-year OS in the internal validation cohort (Fig. 3c). The calibration curves of the risk signature in the two cohorts revealed excellent agreement between the expected and actual outcomes for 3- and 5-year OS (Fig. 3d). Furthermore, the C-index value for both the primary and internal cohorts was 0.77, indicating considerable accuracy.
Establishment and evaluation of a nomogram
The clinical information, including age, sex, TNM stage and grade, and hypoxia score, was remarkably associated with the OS of patients with HNSCC (Table 2). Univariate and multivariate CPHR analyses were carried out to obtain information and risk scores for the primary cohort (Fig. 4a and b). In the primary cohort, we proved those factors were independent prognostic variables of OS, including TNM stage, hypoxia score, and risk score. Further, a prognostic nomogram was established using three independent prognostic variables (Fig. 4c). The miRNA signature was more effective in predicting the OS of HNSCC patients, followed by TNM stage and hypoxia score. Further, the AUC values of the ROC curves of the two independent prognostic variables demonstrated that each variable had credible predictive accuracy, especially the miRNA signature (Fig. 5a). The AUC values of the ROC analysis for nomogram were 0.705, 0.729, and 0.827 at 1-, 3-, and 5-year OS, respectively, for the primary cohort, and 0.723, 0.748, and 0.837 at 1-, 3-, and 5-year OS, respectively, for the internal validation cohort (Fig. 5b). To assess the calibration capability of this prognostic model, we established calibration curves and found that the predicted and actual survival in the two cohorts corresponded using this prognostic model (Fig. 5c). The C-index values of the nomogram were 0.776 and 0.744 in the primary and internal validation cohorts, respectively.
Target genes functional enrichment analysis
We predicted the corresponding target genes using three independent databases to confirm the potential biological functions of the 11 DEMs. In total, 38,191 target genes were detected, of which 305 genes were overlapping. Thus, these overlapping genes potentially modulated by the 11 DEMs were subjected to GO and KEGG enrichment analyses. Based on the criterion of a P-value < 0.05, 42 categories, including nucleus, cytoplasm, and membrane, showed a strong influence on the occurrence and development of cancer (Fig. 6a). In addition, the KEGG enrichment analysis revealed the multiple pathways that play key roles in HNSCC progression, especially the neurotrophin signaling pathway and protein processing in endoplasmic reticulum (Fig. 6b).
Many reports have suggested the participation of miRNAs in different pathological processes related to HNSCC [8, 13]. However, there is no single representative miRNA that predicts the OS of patients with HNSCC. In the present study, we collected information on patients with HNSCC and evaluated their miRNA expression levels to systematically analyze the miRNAs and clinical characteristics associated with their OS.
We performed the CPHR and LASSO regression analyses and identified 11 DEMs, including miR-204-5p, miR-499a-5p, miR-498-5p, miR-155-3p, miR-4714-3p, miR-365a-5p, miR-30a-5p, miR-1-5p, miR-548f-3p, miR-518a-3p, and miR-196b-5p, that were confirmed to influence the OS of patients with HNSCC. Kimura et al. conducted a microarray analysis and found that the expression of miR-30a-5p was upregulated both in HNSCC and esophagus squamous cell carcinoma cell lines compared to that in normal squamous epithelial cell lines . Another study reported the significantly upregulated expression of miR-30a-5p and miR-769-5p in the plasma of patients with oral squamous cell carcinoma (OSCC) and their effectiveness as minimally invasive biomarkers for OSCC diagnosis . However, the function of miR-30a-5p in patients with HNSCC is yet unclear. A previous study found that miR-769-5p promotes the growth of glioma cells through the suppression of KMT2A  and could be useful as a therapeutic target for glioma treatment. The expression of miR-518a-3p has been found to be upregulated in hepatocellular carcinoma, as demonstrated by Taqman low-density miRNA array and quantitative real-time polymerase chain reaction [17, 18], indicative of its tumor-inducing properties. In addition, miR-518a-3p has been shown to promote the proliferation, migration, and invasion of choriocarcinoma cells; knockdown of miR-518a-3p expression has been found to induce S phase arrest in choriocarcinoma cells . miR-155-3p expression downregulation in breast cancer is related to resistance to tumor invasion and metastasis and reduction in paclitaxel resistance . miR-365a-5p inhibits the viability, colony formation, migration, and invasion of non-small cell lung cancer cells by negatively regulating Pellino E3 ubiquitin protein ligase family member 3 (PELI3) , and PELI3 silencing promotes the miR-365a-5p-mediated suppression of pro-tumoral effects. However, the expression of miR-499a-5p was down-regulated in OSCC tissues . During glioma progression, miR-499a-5p exerts tumor-suppressive effects by regulating the formation of vasculogenic mimicry . In addition, the overexpression of miR-499a-5p can markedly inhibit tumor cell proliferation, invasion, and migration . The expression of miR-204-5p is suppressed in HNSCC, wherein it serves as a tumor suppressor and inhibits tumor growth and metastasis . miR-196b-5p expression upregulation in HNSCC has also been verified using TCGA database . Although these authors described three different miRNAs, their study had fewer samples for verification or used different calculation methods to select miRNAs. Therefore, it was imperative to further verify the functions of these miRNAs. The miRNAs miR-498-5p, miR-1-5p, miR-548f-3p, and miR-4714-3p have not been reported in previous HNSCC-related studies; thus, additional research is warranted to determine their functions in HNSCC.
GO analysis revealed 11 DEMs that were primarily enriched at the nucleus, cytoplasm, and membrane and were related to the positive regulation of transcription from the RNA polymerase II promoter. These locations and pathways were associated with different physiological and pathological activities of the human body. Uncontrolled disease activity accelerates the occurrence and development of HNSCC, highlighting the importance of pathological changes related to these locations and pathways. Furthermore, the KEGG enrichment analysis revealed that the neurotrophin signaling pathway, protein processing in the endoplasmic reticulum, and microRNAs in cancer significantly affected HNSCC occurrence and development. We found that several targeted genes exert their functions through the neurotrophin signaling pathway, which performs key functions in the occurrence and development of HNSCC by regulating cell survival, angiogenesis, tissue invasion, DNA damage resistance, and epithelial-to-mesenchymal transition [27,28,29,30]. Gkouveris et al. demonstrated that interferon-gamma inhibits OSCC cell viability and migration and increases cellular apoptosis through the endoplasmic reticulum stress pathway . Further, dentin sialophosphoprotein induces OSCC cell viability  and HIV protease inhibitors with radiosensitizing activities mediate HNSCC cell apoptosis  through the activation of the endoplasmic reticulum stress mechanisms. Overall, these pathways predicted herein may play an important role in HNSCC development and progression.
In the present study, we analyzed the miRNA expression levels and clinical information of 71 HNSCC patients and 13 normal subjects. A mess of samples would effectively reduce bias. The univariate CPHR and LASSO regression analyses yielded 11 DEMs. Thus, we constructed and systematically validated the risk signature and nomogram, which showed high predictive accuracy. In comparison with invasive tests or imageological examination, the risk signature and predicted nomogram were convenient and economic but must be verified in clinical studies and cell or animal experiments.
We constructed an 11-miRNA risk signature and a nomogram that showed high correlation with HNSCC prognosis. The 3- and 5-year OS of patients with HNSCC could be accurately predicted through the risk signature and nomogram, thus serving as an effective tool for clinicians to facilitate decision-making.
Data resources and preprocessing
The miRNA sequencing data and related clinical information of patients with HNSCC were obtained from TCGA database (https://portal.gdc.cancer.gov/). The “edgeR” package in R software was used to identify DEMs from 797 miRNAs in 71 HNSCC samples and 13 normal samples based on |log2FC| ≥ 2 and adjusted P-value < 0.05. We then integrated the miRNA expression data and clinical information based on sample ID, and excluded four patients because they had no significant clinical date or their OS was less than 30 days to avoid introduction of mixed factors. Finally, 67 patients with HNSCC were enrolled as the primary cohort. In order to unbiasedly appraise the final model, 44 patients randomly selected from the primary cohort were designed as the internal validation cohort by using the “caret” package in R software.
Construction of risk score formula
DEMs were detected from the univariate CPHR analysis of the primary cohort. These DEMs with a P-value < 0.05 were subjected to LASSO CPHR analysis to produce a miRNA signature using the “glmnet” package in R software. The risk score of each patient relied on the miRNA expression level and the regression coefficient obtained from the LASSO regression analysis. The risk score was calculated as follows: Risk score (miRNA-based classifier) = sum of coefficients × expression level of miRNAs, refer to previous description . According to the median risk score value, 67 patients with HNSCC were divided into high-risk and low-risk groups.
Assessment of miRNA signature
The survival rate between the high- and low-risk groups were compared via the Kaplan–Meier survival analysis. A time-dependent ROC curve analysis performed using “survival ROC” packages estimated the sensitivity and specificity of each miRNA signature on the basis of the AUC value.
The differences in the signature between predicted survival and actual survival were assessed and calibration curves were constructed. The C-index values were calculated to measure the predictive performance of the risk signature. All analyses were performed to compare the consistency between the primary and internal validated cohorts, refer to previous description .
Construction and evaluation of a novel nomogram
Due to the T stage, M stage, N stage exist serious information loss, those clinical informations were excluded for guaranteeing the reality of this study. Risk signature and other clinical factors, such as age, sex, neoplasm histologic grade (grade), TNM stage, and Ragnum hypoxia score (hypoxia score), were included in the univariate and multivariate CPHR analyses of the primary cohort. Factors with a P-value < 0.05 both in the univariate and multivariate CPHR analyses were confirmed as independent prognostic variables and then used to establish a novel nomogram by using the package “rms” in R software. To evaluate the predictive ability of the novel nomogram, a calibration plot was constructed and AUC in ROC analysis and C-index were calculated. All analyses were implemented for the primary and internal validation cohorts.
Functional enrichment analysis
Three online databases, including TargetScan, miRTarBase, and miRDB, were used to predict the potential target genes of prognostic miRNAs. The overlapping miRNA target genes predicted by these three databases were undergo to GO analysis and KEGG pathway enrichment analysis using DAVID (http://david.ncifcrf.gov/) refer to previous description . P < 0.05 was considered the cut-off value.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Head and neck squamous cell carcinoma
Kyoto Encyclopedia of Genes and Genomes
Differentially expressed miRNA
The Cancer Genome Atlas
Cox proportional hazard regression
Oral squamous cell carcinoma
Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136(5):E359–86.
Leemans C, Braakhuis B, Brakenhoff R. The molecular biology of head and neck cancer. Nat Rev Cancer. 2011;11:9–22.
Castellsagué X, Alemany L, Quer M, Halec G, Quirós B, Tous S, et al. HPV involvement in head and neck cancers: comprehensive assessment of biomarkers in 3680 patients. J Natl Cancer Inst. 2016;108:1–12.
Kim L, King T, Agulnik M. Head and neck cancer: changing epidemiology and public health implications. Head Neck. 2010;24(10):915–9.
Gregory RI, Shiekhattar R. MicroRNA biogenesis and cancer. Cancer Res. 2005;65(9):3509–12.
Gambari R, Brognara E, Spandidos DA, Fabbri E, et al. Targeting oncomiRNAs and mimicking tumor suppressor miRNAs: new trends in the development of miRNA therapeutic strategies in oncology (review). Int J Oncol. 2016;49(1):5–32.
Saliminejad K, Khorram Khorshid HR, Soleymani Fard S, Ghaffari SH, et al. An overview of microRNAs: biology, functions, therapeutics, and analysis methods. J Cell Physiol. 2019;234(5):5451–65.
Fadhil RS, Wei MQ, Nikolarakos D, Good D, Nair RG. Salivary microRNA miR-let-7a-5p and miR-3928 could be used as potential diagnostic bio-markers for head and neck squamous cell carcinoma. PLoS One. 2020;15(3):e221779.
Lin XJ, Liu H, Li P, Wang HF, Yang AK, Di JM, et al. miR-936 suppresses cell proliferation, invasion, and drug resistance of laryngeal squamous cell carcinoma and targets GPR78. Front. Oncol. 2020;10:60.
Liu C, Yu Z, Huang S, Zhao Q, Sun Z, Fletcher C, et al. Combined identification of three miRNAs in serum as effective diagnostic biomarkers for HNSCC. EBioMedicine. 2019;50:135–43.
Kovaříková J, Baranová I, Laco J, Rozkošová K, Vošmíková H, Vošmík M, et al. Deregulation of Selected MicroRNAs in Sinonasal Squamous Cell Carcinoma: Searching for Potential Prognostic Biomarkers. Folia Biol (Praha). 2019;65(3):142–51.
Shiah SG, Hsiao JR, Chang HJ, Hsu YM, Wu GH, Peng HY, et al. MiR-30a and miR-379 modulate retinoic acid pathway by targeting DNA methyltransferase 3B in oral cancer. J Biomed Sci. 2020;27(1):46.
Hao Y, Zhang D, Guo Y, Fu Z, Yu D, Guan G, et al. miR-488-3p sponged by circ-0000495 and mediated upregulation of TROP2 in head and neck squamous cell carcinoma development. J Cancer. 2020;11(11):3375–86.
Kimura S, Naganuma S, Susuki D, Hirono Y, Yamaguchi A, Fujieda S, et al. Expression of microRNAs in squamous cell carcinoma of human head and neck and the esophagus: miR-205 and miR-21 are specific markers for HNSCC and ESCC. Oncol. Rep. 2010;23(6):1625–33.
Pedersen NJ, Jensen DH, Lelkaitis G, Kiss K, Charabi BW, Ullum H, et al. MicroRNA-based classifiers for diagnosis of oral cavity squamous cell carcinoma in tissue and plasma. Oral Oncol. 2018;83:46–52.
Chang M, Yan P, Zhang B, Zhang G, Wang J, Ge H, et al. MicroRNA-769-5p promotes the growth of Glioma cells by targeting lysine methyltransferase 2A. Onco Targets Ther. 2019;12:9177–87.
Wang W, Zhao LJ, Tan YX, Ren H, Qi ZT. Identification of deregulated miRNAs and their targets in hepatitis B virus-associated hepatocellular carcinoma. World J. Gastroenterol. 2012;18(38):5442–53.
Wang W, Zhao LJ, Tan YX, Ren H, Qi ZT. MiR-138 induces cell cycle arrest by targeting cyclin D3 in hepatocellular carcinoma. Carcinogenesis. 2012;33(5):1113–20.
Zhao JR, Cheng WW, Wang YX, Cai M, Wu WB, Zhang HJ. Identification of microRNA signature in the progression of gestational trophoblastic disease. Cell Death Dis. 2018;9(2):94.
Zhang L, Chen T, Yan L, Xu H, Wang Y, Li Y, et al. MiR-155-3p acts as a tumor suppressor and reverses paclitaxel resistance via negative regulation of MYD88 in human breast cancer. Gene. 2019;700:85–95.
He Y, Shi Y, Liu R, Wang Z, Wang B, Li S, et al. PELI3 mediates pro-tumor actions of down-regulated miR-365a-5p in non-small cell lung cancer. Biol Res. 2019;52(1):24.
Hou YY, Lee JH, Chen HC, Yang CM, Huang SJ, Liou HH, et al. The association between miR-499a polymorphism and oral squamous cell carcinoma progression. Oral Dis. 2015;21(2):195–206.
Wang M, Yang C, Liu X, Zheng J, Xue Y, Ruan X, et al. An upstream open reading frame regulates vasculogenic mimicry of glioma via ZNRD1-AS1/miR-499a-5p/ELF1/EMI1 pathway. J Cell Mol Med. 2020;24(11):6120–36.
Gu X, Dong M, Liu Z, Yang J, Shi Y. MiR-499a-5p inhibits proliferation, invasion, migration, and epithelial-Mesenchymal transition, and enhances Radiosensitivity of cervical Cancer cells via targeting eIF4E. Onco Targets Ther. 2020;13:2913–24.
Zhuang Z, Yu P, Xie N, Wu Y, Liu H, Zhang M, et al. MicroRNA-204-5p is a tumor suppressor and potential therapeutic target in head and neck squamous cell carcinoma. Theranostics. 2020;10(3):1433–53.
Nunez Lopez YO, Victoria B, Golusinski P, Golusinski W, Masternak M. Characteristic miRNA expression signature and random forest survival analysis identify potential cancer-driving miRNAs in a broad range of head and neck squamous cell carcinoma subtypes. Rep Pract Oncol Radiother. 2018;23(1):6–20.
Chen-Tsai CP, Colome-Grimmer M, Wagner RF. Correlations among neural cell adhesion molecule, nerve growth factor, and its receptors, TrkA, TrkB, TrkC, and p75, in perineural invasion by basal cell and cutaneous squamous cell carcinomas. Dermatol Surg. 2004;30(7):1009–16.
Kupferman ME, Jiffar T, El-Naggar A, Yilmaz T, Zhou G, Xie T, et al. TrkB induces EMT and has a key role in invasion of head and neck squamous cell carcinoma. Oncogene. 2010;29(14):2047–59.
Sasahira T, Ueda N, Yamamoto K, Bhawal UK, Kurihara M, Kirita T, et al. Trks are novel oncogenes involved in the induction of neovascularization, tumor progression, and nodal metastasis in oral squamous cell carcinoma. Clin. Exp. Metastasis. 2013;30(2):165–76.
Yilmaz T, Jiffar T, de la Garza G, Lin H, Milas Z, Takahashi Y, et al. Theraputic targeting of Trk supresses tumor proliferation and enhances cisplatin activity in HNSCC. Cancer biol. Ther. 2010;10(6):644–53.
Gkouveris I, Nikitakis NG, Aseervatham J, Ogbureke KUE. Interferon γ suppresses dentin sialophosphoprotein in oral squamous cell carcinoma cells resulting in antitumor effects, via modulation of the endoplasmic reticulum response. Int. J. Oncol. 2018;53(6):2423–32.
Gkouveris I, Nikitakis NG, Aseervatham J, Ogbureke K. The tumorigenic role of DSPP and its potential regulation of the unfolded protein response and ER stress in oral cancer cells. Int J Oncol. 2018;53(4):1743–51.
Liu R, Zhang L, Yang J, Zhang X, Mikkelsen R, Son S, et al. HIV Protease Inhibitors Sensitize Human Head and Neck Squamous Carcinoma Cells to Radiation by Activating Endoplasmic Reticulum Stress. PLoS One. 2015;10(5):e125928.
Zheng Z, Wen Y, Nie K, Tang S, Chen X, Lan S, et al. Construction of a 13-microRNA-based signature and prognostic nomogram for predicting overall survival in patients with hepatocellular carcinoma. Hepatol Res. 2020. https://0-doi-org.brum.beds.ac.uk/10.1111/hepr.13538.
This work was supported by the National Natural Science Foundation of China (No. 81700512).
Ethics approval and consent to participate
Consent for publication
The authors declared no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Huang, Y., Liu, Z., Zhong, L. et al. Construction of an 11-microRNA-based signature and a prognostic nomogram to predict the overall survival of head and neck squamous cell carcinoma patients. BMC Genomics 21, 691 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12864-020-07104-w
- Head and neck squamous cell carcinoma
- Overall survival
- Risk signature