Oral submucous fibrosis (OSF) as one of the premalignant disorders endures a series of histopathological stages to invasive oral squamous cell carcinoma (OSCC) eventually. However, the role of long non-coding RNA (lncRNA) expression in OSF malignant progression still remains poorly understood. Through RNA-sequencing normal mucous, OSF and OSCC tissues, we found 687 lncRNA transcripts significantly and differentially expressed during OSF progression, including 231 upregulated lncRNAs and 456 downregulated lncRNAs, indicating that lncRNAs are involved in the regulation of different stages of OSF development. Further functional enrichment analysis showed these differentially expressed lncRNAs participated in inflammation signaling, Wnt signaling, angiogenesis, CCKR signaling, integrin signaling, PDGF signaling, p53 signaling, and EGF receptor (EGFR) signaling pathways, which contribute to inflammatory and fibroelastic pathogenetic changes of OSF and further malignant progression. Five novel lncRNAs were differentially expressed during OSF progression with varied expression levels, indicating the importance of these lncRNAs in OSF malignant development. Moreover, some lncRNAs have been previously identified to be associated with OSCC pathogenesis, including HCG22, RP11-397A16.1, LINC00271, CTD-3179P9.1, and ZNF667-AS1. Thus, our study firstly comprehensively elucidated lncRNAs expression profile of malignant procession from OSF premalignant lesion to OSCC, which will enlighten our understanding of the importance of lncRNA involved in OSF malignant development.

Background Head and neck squamous cell carcinoma (HNSCC) is a common tumor worldwide with poor prognosis. The pathogenesis of human papillomavirus (HPV)-positive and HPV-negative HNSCCs differs. However, few studies have considered the HPV status when identifying biomarkers for HNSCC. Thus, the identification of biomarkers for HPV-positive and HPV-negative HNSCCs is urgently needed. Methods Three microarray datasets from Gene Expression Omnibus (GEO) were analyzed, and the differentially expressed genes (DEGs) were obtained. Then, functional enrichment pathway analysis was performed and protein–protein interaction (PPI) networks were constructed. The expression of hub genes at both the mRNA and protein level was determined in Oncomine, The Cancer Genome Atlas (TCGA) and the Human Protein Atlas (HPA). In addition, survival analysis of the patient stratified by HPV status and the expression levels of key genes were performed based on TCGA data. The role of AREG, STAG3, CAV1 and C19orf57 in cancer were analyzed through Gene set enrichment analysis (GSEA). The top ten small molecule drugs were identified and the therapeutic value of zonisamide, NVP-AUY922, PP-2 and fostamatinib was further evaluated in six HPV-negative HNSCC cell lines. Finally, the therapeutic value of NVP-AUY922 was tested in vivo based on three HPV-negative HNSCC models, and statistical analysis was performed. Results In total, 47 DEGs were obtained, 11 of which were identified as hub genes. Biological process analysis indicated that the hub genes were associated with the G1/S transition of the mitotic cell cycle. Survival analysis uncovered that the prognostic value of AREG, STAG3, C19orf57 and CAV1 differed between HPV-positive and HPV-negative patients. Gene set enrichment analysis (GSEA) showed the role of AREG, STAG3 and CAV1 in dysregulated pathways of tumor. Ten small molecules were identified as potential drugs specifically for HPV-positive or HPV-negative patients; three—NVP-AUY922, fostamatinib and PP-2—greatly inhibited the proliferation of six HPV-negative HNSCC cell lines in vitro, and NVP-AUY922 inhibited three HPV-negative HNSCC xenografts in vivo. Conclusions In conclusion, AREG, STAG3, C19orf57 and CAV1 are key prognostic factors and potential therapeutic targets in HPV-negative HNSCC. NVP-AUY922, fostamatinib and PP-2 may be effective drugs for HPV-negative HNSCC.