Background Epidermal growth factor receptor (EGFR) is both a driver oncogene and a therapeutic target in advanced head and neck squamous cell carcinoma (HNSCC). However, response to EGFR treatment is inconsistent and lacks markers for treatment prediction. This study investigated EGFR-induced epithelial-to-mesenchymal transition (EMT) as a central parameter in tumor progression and identified novel prognostic and therapeutic targets, and a candidate predictive marker for EGFR therapy response. Methods Transcriptomic profiles were analyzed by RNA sequencing (RNA-seq) following EGFR-mediated EMT in responsive human HNSCC cell lines. Exclusive genes were extracted via differentially expressed genes (DEGs) and a risk score was determined through forward feature selection and Cox regression models in HNSCC cohorts. Functional characterization of selected prognostic genes was conducted in 2D and 3D cellular models, and findings were validated by immunohistochemistry in primary HNSCC. Results An EGFR-mediated EMT gene signature composed of n = 171 genes was identified in responsive cell lines and transferred to the TCGA-HNSCC cohort. A 5-gene risk score comprising DDIT4, FADD, ITGB4, NCEH1, and TIMP1 prognosticated overall survival (OS) in TCGA and was confirmed in independent HNSCC cohorts. The EGFR-mediated EMT signature was distinct from EMT hallmark and partial EMT (pEMT) meta-programs with a differing enrichment pattern in single malignant cells. Molecular characterization showed that ITGB4 was upregulated in primary tumors and metastases compared to normal mucosa and correlated with EGFR/MAPK activity in tumor bulk and single malignant cells. Preferential localization of ITGB4 together with its ligand laminin 5 at tumor-stroma interfaces correlated with increased tumor budding in primary HNSCC tissue sections. In vitro, ITGB4 knock-down reduced EGFR-mediated migration and invasion and ITGB4-antagonizing antibody ASC8 impaired 2D and 3D invasion. Furthermore, a logistic regression model defined ITGB4 as a predictive marker of progression-free survival in response to Cetuximab in recurrent metastatic HNSCC patients. Conclusions EGFR-mediated EMT conveyed through MAPK activation contributes to HNSCC progression upon induction of migration and invasion. A 5-gene risk score based on a novel EGFR-mediated EMT signature prognosticated survival of HNSCC patients and determined ITGB4 as potential therapeutic and predictive target in patients with strong EGFR-mediated EMT.

Partial epithelial-to-mesenchymal transition (pEMT) contributes to cellular heterogeneity that is associated with nodal metastases and unfavorable clinical parameters in head and neck squamous cell carcinomas (HNSCCs). We developed a single-cell RNA sequencing signature-based pEMT quantification through cell type-dependent deconvolution of bulk RNA sequencing and microarray data combined with single-sample scoring of molecular phenotypes (Singscoring). Clinical pEMT-Singscores served as molecular classifiers in multivariable Cox proportional hazard models and high scores prognosticated poor overall survival and reduced response to irradiation as independent parameters in large HNSCC cohorts [The Cancer Genome Atlas (TCGA), MD Anderson Cancer Centre (MDACC), Fred Hutchinson Cancer Research Center (FHCRC)]. Differentially expressed genes confirmed enhanced cell motility and reduced oxidative phosphorylation and epithelial differentiation in pEMT high patients. In patients and cell lines, the EMT transcription factor SLUG correlated most strongly with pEMT-Singscores and promoted pEMT, enhanced invasion, and resistance to irradiation in vitro. SLUG protein levels in HNSCC predicted disease-free survival, and its peripheral expression at the interphase to the tumor microenvironment was significantly increased in relapsing patients. Hence, pEMT-Singscores represent a novel risk predictor for HNSCC stratification regarding clinical outcome and therapy response that is partly controlled by SLUG.

Ulceration and immune status are independent prognostic factors for survival in melanoma patients. Herein univariate Cox regression analysis revealed 53 ulcer-immunity-related DEGs. We performed consensus clustering to divide The Cancer Genome Atlas (TCGA) cohort (n = 467) into three subtypes with different prognosis and biological functions, followed by validation in three merged Gene Expression Omnibus (GEO) cohorts (n = 399). Multiomics approach was used to assess differences among the subtypes. Cluster 3 showed relatively lesser amplification and expression of immune checkpoint genes. Moreover, Cluster 3 lacked immune-related pathways and immune cell infiltration, and had higher proportion of non-responders to immunotherapy. We also constructed a prognostic model based on ulceration and immune related genes in melanoma. EIF3B was a hub gene in the intersection between genes specific to Cluster 3 and those pivotal for melanoma growth (DepMap, https://depmap.org/portal/download/). High EIF3B expression in TCGA and GEO datasets was related to worst prognosis. In vitro models revealed that EIF3B knockdown inhibited melanoma cell migration and invasion, and decreased TGF-β1 level in supernatant compared with si-NC cells. EIF3B expression was negatively correlated with immune-related signaling pathways, immune cell gene signatures, and immune checkpoint gene expression. Moreover, its low expression could predict partial response to anti-PD-1 immunotherapy. To summarize, we established a prognostic model for melanoma and identified the role of EIF3B in melanoma progression and immunotherapy resistance development.

Rationale The M2-like tumor-associated macrophages (TAMs) are independent prognostic factors in melanoma. Methods We performed weighted gene co-expression network analysis (WGCNA) to identify the module most correlated with M2-like TAMs. The Cancer Genome Atlas (TCGA) patients were classified into two clusters that differed based on prognosis and biological function, with consensus clustering. A prognostic model was established based on the differentially expressed genes (DEGs) of the two clusters. We investigated the difference in immune cell infiltration and immune response-related gene expression between the high and low risk score groups. Results The risk score was defined as an independent prognostic value in melanoma. VARS1 was a hub gene in the M2-like macrophage-associated WGCNA module that the DepMap portal demonstrated was necessary for melanoma growth. Overexpressing VARS1 in vitro increased melanoma cell migration and invasion, while downregulating VARS1 had the opposite result. VARS1 overexpression promoted M2 macrophage polarization and increased TGF-β1 concentrations in tumor cell supernatant in vitro. VARS1 expression was inversely correlated with immune-related signaling pathways and the expression of several immune checkpoint genes. In addition, the VARS1 expression level helped predict the response to anti-PD-1 immunotherapy. Pan-cancer analysis demonstrated that VARS1 expression negatively correlated with CD8 T cell infiltration and the immune response-related pathways in most cancers. Conclusion We established an M2-like TAM-related prognostic model for melanoma and explored the role of VARS1 in melanoma progression, M2 macrophage polarization, and the development of immunotherapy resistance.

No abstract