Identification and validation of novel prognostic signatures based on m5C methylation patterns and tumor EMT profiles in head and neck squamous cell carcinoma

Abstract The role of 5-methylcytosine (m5C) in tumor initiation and progression has been increasingly recognized. However, the precise association between the regulation of m5C and the progression, metastasis, and prognosis of head and neck squamous cell carcinoma (HNSCC) has not yet been fully explored. Data from 545 HNSCC patients obtained from The Cancer Genome Atlas (TCGA) database were analyzed. Unsupervised cluster analysis was conducted using the expression levels of m5C regulatory genes. Additionally, gene set variation analysis (GSVA), single-sample gene set enrichment analysis (ssGSEA), and Cox regression analysis were utilized. Quantitative reverse transcription polymerase chain reaction (RT-qPCR), colony formation assay, transwell experiments and western blots were performed in the HNSCC cell line UM-SCC-17B to assess the expression and functional role of one of the novel signatures, CNFN. Significant expression differences were found in m5C regulatory genes between tumor and normal tissues in HNSCC. Two distinct m5C modification patterns, characterized by substantial prognostic differences, were identified. Cluster-2, which exhibited a strong association with epithelial-mesenchymal transition (EMT), was found to be associated with a poorer prognosis. Based on the m5C clusters and EMT status, differentially expressed genes (DEGs) were identified. Using DEGs, an 8-gene signature (CAMK2N1, WNT7A, F2RL1, AREG, DEFB1, CNFN, TGFBI, and CAV1) was established to develop a prognostic model. The performance of this signature was validated in both the training and external validation datasets, demonstrating its promising efficacy. Furthermore, additional investigations using RT-qPCR on clinical specimens and experimental assays in cell lines provided compelling evidence suggesting that CNFN, one of the genes in the signature, could play a role in HNSCC progression and metastasis through the EMT pathway. This study highlighted the role of m5C in HNSCC progression and metastasis. The relationship between m5C and EMT has been elucidated for the first time. A robust prognostic model was developed for accurately predicting HNSCC patients’ survival outcomes. Potential molecular mechanisms underlying these associations have been illuminated through this research.