LINC00958 and HOXC13-AS as key candidate biomarkers in head and neck squamous cell carcinoma by integrated bioinformatics analysis

Background: Head and neck squamous cell carcinoma (HNSCC) is a malignant tumor with a strong tendency for metastasis and recurrence. Finding effective biomarkers for the early diagnosis of HNSCC is critical for the early treatment and prognosis of patients. Methods: RNA sequencing data including long non-coding RNAs (lncRNAs), messenger RNA (mRNAs) and microRNAs (miRNAs) of 141 HNSCC and 44 adjacent normal tissues were obtained from the TCGA. Differentially expressed genes were analyzed using the R package DESeq. GO terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted. A competing endogenous RNAs (ceRNA) network was constructed. The most differentially expressed genes in the main ceRNA network were chosen for nasopharyngeal carcinoma (NPC) cell lines and NPEC2 Bmi-1 cell line verification. A receiver operating characteristic (ROC) curve was constructed for 141 specimens of HNSCC tissues from 44 control samples. Results: In our study, 79 HNSCC-associated abnormally expressed lncRNAs , 86 abnormally expressed miRNAs and 324 abnormally expressed mRNAs were identified. The public microarray results showed that LINC00958 and HOXC13-AS expression levels were upregulated in HNSCC tissues compared with the adjacent normal tissues in this study (p < 0.0001). LINC00958 and HOXC13-AS expression levels in NPC cell lines were higher than those in the NPEC2 Bmi-1 cell line (p < 0.05). The results showed that the area under the ROC curve (AUC) of LINC00958 reached up to 0.906 at a cutoff value of 7.96, with a sensitivity and specificity of 80.85% and 90.91%, respectively. The AUC of HOXC13-AS reached up to 0.898 at a cutoff value of 0.695, with sensitivity and specificity values of 86.23% and 83.78%, respectively. Conclusion: The current study indicates that LINC00958 and HOXC13-AS are new candidate diagnostic biomarkers for HNSCC patients.