Comprehensive analysis reveals COPB2 and RYK associated with tumor stages of larynx squamous cell carcinoma

Background Laryngeal squamous cell carcinoma (LSCC) is one of the highly aggressive malignancy types of head and neck squamous cell carcinomas; genes involved in the development of LSCC still need exploration. Methods We downloaded expression profiles of 96 (85 in advanced stage and 11 in early stage) LSCC patients from TCGA-HNSC. Function enrichment and protein-protein interactions of genes in significant modules were conducted. Univariate and multivariate Cox regression analyses were performed to explore potential prognostic biomarkers for LSCC. The expression levels of genes at different stages were compared and visualized via boxplots. Immune infiltration was examined by the CIBERSORTx web-based tool and depicted with ggplot2. Gene set enrichment analysis (GSEA) was utilized to analyze functional enrichment terms and pathways. Immunohistochemical staining (IHC) was used to verify the expression of genes in the LSCC samples. Results We identified 25 modules, including 3 modules significantly related to tumor stages of LSCC via weighted gene co-expression network analysis (WGCNA). UIMC1, NPM1, and DCTN4 in the module 'cyan', TARS in the module 'darkorange', and COPB2 and RYK in the module 'lightyellow' showed statistically significant relation to overall survival. The expression of COPB2, DCTN4, RYK, TARS, and UIMC1 indicated association with the change of fraction of immune cells in LSCC patients; two genes, COPB2 and RYK, indicated different expression in various tumor stages of LSCC. Finally, COPB2 and RYK showed high-expression in tumor tissues of advanced LSCC patients. Conclusions Our study provided a potential perceptive in analyzing progression of LSCC cells and exploring prognostic genes.