Abstract
The fatty acid metabolism (FAM) is known to impact tumorigenesis, tumor progression and treatment resistance via enhancing lipid synthesis, storage and catabolism. However, the role of FAM in head and neck squamous cell carcinoma (HNSCC) has remained elusive. In the present study, we obtained a total of 69 differentially expressed FAM-related genes between 502 HNSCC samples and 44 normal samples from The Cancer Genome Atlas (TCGA) database. The HNSCC samples were divided into 2 clusters according to 69 differentially expressed genes (DEGs) via cluster analysis. Then DEGs in the two clusters were found, and 137 prognostic DEGs were identified by univariate analysis. Subsequently, combined with the clinical information of 546 HNSCC patients from TCGA database, a 12-gene prognostic risk model was established (FEPHX3, SPINK7, FCRLA, MASP1, ZNF541, CD5, BEST2 and ZAP70 were down-regulation, ADPRHL1, DYNC1I1, KCNG1 and LINC00460 were up-regulation) using multivariate Cox regression and LASSO regression analysis. The risk scores of 546 HNSCC samples were calculated. According to the median risk score, 546 HNSCC patients were divided into the high- and low-risk (high- and low score) groups. The Kaplan-Meier survival analysis showed that the survival time of HNSCC patients was significantly shorter in the high-risk group than that in the low-risk group (p < 0.001). The same conclusion was obtained in the Gene Expression Omnibus (GEO) dataset. After that, the multivariate Cox regression analysis indicated that the risk score was an independent factor for patients with HNSCC in the TCGA cohort. In addition, single-sample gene set enrichment analysis (ssGSEA) indicated that the level of infiltrating immune cells was relatively low in the high-risk group compared with the low-risk group. In summary, FAM-related gene expression-based risk signature could predict the prognosis of HNSCC independently.