A transcriptomic map of EGFR-induced epithelial-to-mesenchymal transition identifies prognostic and therapeutic targets for head and neck cancer

High glutamine increases stroke risk by inducing the endothelial-to-mesenchymal transition in moyamoya disease

At present, there is limited research on the mechanisms underlying moyamoya disease (MMD). Herein, we aimed to determine the role of glutamine in MMD pathogenesis, and 360 adult patients were prospectively enrolled. Human brain microvascular endothelial cells (HBMECs) were subjected to Integrin Subunit Beta 4 (ITGB4) overexpression or knockdown and atorvastatin. We assessed factors associated with various signaling pathways in the context of the endothelial-to-mesenchymal transition (EndMT), and the expression level of related proteins was validated in the superficial temporal arteries of patients. We found glutamine levels were positively associated with a greater risk of stroke (OR = 1.599, p = 0.022). After treatment with glutamine, HBMECs exhibited enhanced proliferation, migration, and EndMT, all reversed by ITGB4 knockdown. In ITGB4-transfected HBMECs, the MAPK-ERK-TGF-β/BMP pathway was activated, with Smad4 knockdown reversing the EndMT. Furthermore, atorvastatin suppressed the EndMT by inhibiting Smad1/5 phosphorylation and promoting Smad4 ubiquitination in ITGB4-transfected HBMECs. We also found the protein level of ITGB4 was upregulated in the superficial temporal arteries of patients with MMD. In conclusion, our study suggests that glutamine may be an independent risk factor for hemorrhage or infarction in patients with MMD and targeting ITGB4 could potentially be therapeutic approaches for MMD.

Tertiary Lymphoid Structures Gene Signature Predicts Prognosis and Immune Infiltration Analysis in Head and Neck Squamous Cell Carcinoma

Objectives

This study aims to assess the prognostic implications of gene signature of the tertiary lymphoid structures (TLSs) in head and neck squamous cell carcinoma (HNSCC) and scrutinize the influence of TLS on immune infiltration.

Methods

Patients with HNSCC from the Cancer Genome Atlas were categorized into high/low TLS signature groups based on the predetermined TLS signature threshold. The association of the TLS signature with the immune microenvironment, driver gene mutation status, and tumor mutational load was systematically analyzed. Validation was conducted using independent datasets (GSE41613 and GSE102349).

Results

Patients with a high TLS signature score exhibited better prognosis compared to those with a low TLS signature score. The group with a high TLS signature score had significantly higher immune cell subpopulations compared to the group with a low TLS signature score. Moreover, the major immune cell subpopulations and immune circulation characteristics in the tumor immune microenvironment were positively correlated with the TLS signature. Mutational differences in driver genes were observed between the TLS signature high/low groups, primarily in the cell cycle and NRF2 signaling pathways. Patients with TP53 mutations and high TLS signature scores demonstrated a better prognosis compared to those with TP53 wild-type. In the independent cohort, the relationship between TLS signatures and patient prognosis and immune infiltration was also confirmed. Additionally, immune-related biological processes and signaling pathways were activated with elevated TLS signature.

Conclusion

High TLS signature is a promising independent prognostic factor for HNSCC patients. Immunological analysis indicated a correlation between TLS and immune cell infiltration in HNSCC. These findings provide a theoretical basis for future applications of TLS signature in HNSCC prognosis and immunotherapy.

Evaluation of 5'-tRF-His-GTG As a Molecular Biomarker in Breast Cancer Diagnoses and Prognosis

Background: Breast cancer (BC) is the most prevalent cancer among women worldwide. Although advances have been made in the identification of predictive biomarkers, current options for early diagnosis and prognostic analysis are still suboptimal. Recently, transfer-RNA-derived RNA fragments (tRFs) have emerged as a class of small noncoding RNAs that play a role in the cancer progression. The authors aimed to identify a specific class of tRFs as a molecular marker for BC diagnosis and prognosis in clinical management. Methods: The levels of 5'-tRF-His-GTG were quantified in BC tissue (n = 101) and inflammatory normal breast tissue (n = 22) using in situ hybridization. Clinicopathological parameters were obtained, including age, tumor node metastasis stage, hormone receptor status, histopathological grade, lymphovascular invasion, and recurrence. The correlation between the expression of 5'-tRF-His-GTG and these parameters in different BC subtypes was analyzed. Patient death and cancer progression were regarded as clinical endpoints in the survival analysis. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were also performed to predict the involvement in pivotal biological process. Results: The expression of 5'-tRF-His-GTG was significantly downregulated in BC tissues and was in connection with T stage in human epidermal growth factor 2-positive and basal-like BC, as well as N stage and histopathological grade in luminal BC. Patients with low expression of 5'-tRF-His-GTG had a poor overall survival rate. Statistics of GO and KEGG pathway revealed that cation channel activity, protein catabolic process, response to temperature stimulus, cell cycle, focal adhesion, and glycerophospholipid metabolism were significantly enriched. Conclusions: This study suggests that the assessment of 5'-tRF-His-GTG expression could serve as a novel biomarker for individual diagnosis and prognosis in BC.

Exercise potentially prevents colorectal cancer liver metastases by suppressing tumor epithelial cell stemness via RPS4X downregulation

Background

Colorectal cancer (CRC) is the third most prevalent tumor globally. The liver is the most common site for CRC metastasis, and the involvement of the liver is a common cause of death in patients with late-stage CRC. Consequently, mitigating CRC liver metastasis (CRLM) is key to improving CRC prognosis and increasing survival. Exercise has been shown to be an effective method of improving the prognosis of many tumor types. However, the ability of exercise to inhibit CRLM is yet to be thoroughly investigated.

Methods

The GSE157600 and GSE97084 datasets were used for analysis. A pan-cancer dataset which was uniformly normalized was downloaded and analyzed from the UCSC database: TCGA, TARGET, GTEx (PANCAN, n = 19,131, G = 60,499). Several advanced bioinformatics analyses were conducted, including single-cell sequencing analysis, correlation algorithm, and prognostic screen. CRC tumor microarray (TMA) as well as cell/animal experiments are used to further validate the results of the analysis.

Results

The greatest variability was found in epithelial cells from the tumor group. RPS4X was generally upregulated in all types of CRC, while exercise downregulated RPS4X expression. A lowered expression of RPS4X may prolong tumor survival and reduce CRC metastasis. RPS4X and tumor stemness marker-CD44 were highly positively correlated and knockdown of RPS4X expression reduced tumor stemness both in vitro and in vivo.

Conclusion

RPS4X upregulation may enhance CRC stemness and increase the odds of metastasis. Exercise may reduce CRC metastasis through the regulation of RPS4X.

Monoclonal antibody cetuximab impairs matrix metalloproteinases 2 and 9, epithelial-mesenchymal transition and cell migration in feline oral squamous cell carcinoma in vitro

Feline oral squamous cell carcinoma (FOSCC) is characterised by invasive and metastatic behaviour and is poorly responsive to current treatments, hence the need for new therapeutic strategies. FOSCC shares molecular targets with human head and neck squamous cell carcinoma (HNSCC), among these the epidermal growth factor receptor. Cetuximab is an anti-epidermal growth factor receptor monoclonal antibody employed in the therapy of HNSCC and, interestingly, previous work in vitro suggested that it displays cytostatic and cytotoxic properties also against FOSCC. With the present study, we aimed at further investigating the effects of cetuximab on invasion and metastasis pathways proven to be relevant in human patients. To this purpose, FOSCC cell lines SCCF1, SCCF2 and SCCF3 were treated with cetuximab for 48/72 h and subjected to Western blot for matrix metalloproteinases-2/9 (MMP-2/9) and epithelial-mesenchymal transition markers vimentin, E-, P- and N-cadherin. Treatment with cetuximab resulted in downregulation of MMP-2/-9 in all of the three cell lines in a dose-dependent manner. Moreover, cetuximab downregulated vimentin and P-cadherin in SCCF1, upregulated E-cadherin whilst downregulating P-/N-cadherins in SCCF2, and impaired P-/N-cadherins in SCCF3. An in vitro scratch test also demonstrated that cetuximab delayed cell migration in SCCF3. These data suggest that cetuximab mitigates invasion and metastasis processes by impairing MMPs and epithelial-mesenchymal transition pathways in FOSCC, indicating that this monoclonal antibody may help to counteract malignant progression and improve the management of locally invasive disease.

Development of a promising PPAR signaling pathway-related prognostic prediction model for hepatocellular carcinoma

The peroxisome proliferator-activated receptor (PPAR) signaling pathway plays a crucial role in systemic cell metabolism, energy homeostasis and immune response inhibition. However, its significance in hepatocellular carcinoma (HCC) has not been well documented. In our study, based on the RNA sequencing data of HCC, consensus clustering analyses were performed to identify PPAR signaling pathway-related molecular subtypes, each of which displaying varying survival probabilities and immune infiltration status. Following, a prognostic prediction model of HCC was developed by using the random survival forest method and Cox regression analysis. Significant difference in survival outcome, immune landscape, drug sensitivity and pathological features were observed between patients with different prognosis. Additionally, decision tree and nomogram models were adopted to optimize the prognostic prediction model. Furthermore, the robustness of the model was verified through single-cell RNA-sequencing data. Collectively, this study systematically elucidated that the PPAR signaling pathway-related prognostic model has good predictive efficacy for patients with HCC. These findings provide valuable insights for further research on personalized treatment approaches for HCC.

Resistance of HNSCC cell models to pan-FGFR inhibition depends on the EMT phenotype associating with clinical outcome

BACKGROUND

Focal adhesion signaling involving receptor tyrosine kinases (RTK) and integrins co-controls cancer cell survival and therapy resistance. However, co-dependencies between these receptors and therapeutically exploitable vulnerabilities remain largely elusive in HPV-negative head and neck squamous cell carcinoma (HNSCC).

METHODS

The cytotoxic and radiochemosensitizing potential of targeting 10 RTK and β1 integrin was determined in up to 20 3D matrix-grown HNSCC cell models followed by drug screening and patient-derived organoid validation. RNA sequencing and protein-based biochemical assays were performed for molecular characterization. Bioinformatically identified transcriptomic signatures were applied to patient cohorts.

RESULTS

Fibroblast growth factor receptor (FGFR 1-4) targeting exhibited the strongest cytotoxic and radiosensitizing effects as monotherapy and combined with β1 integrin inhibition, exceeding the efficacy of the other RTK studied. Pharmacological pan-FGFR inhibition elicited responses ranging from cytotoxicity/radiochemosensitization to resistance/radiation protection. RNA sequence analysis revealed a mesenchymal-to-epithelial transition (MET) in sensitive cell models, whereas resistant cell models exhibited a partial epithelial-to-mesenchymal transition (EMT). Accordingly, inhibition of EMT-associated kinases such as EGFR caused reduced adaptive resistance and enhanced (radio)sensitization to FGFR inhibition cell model- and organoid-dependently. Transferring the EMT-associated transcriptomic profiles to HNSCC patient cohorts not only demonstrated their prognostic value but also provided a conclusive validation of the presence of EGFR-related vulnerabilities that can be strategically exploited for therapeutic interventions.

CONCLUSIONS

This study demonstrates that pan-FGFR inhibition elicits a beneficial radiochemosensitizing and a detrimental radioprotective potential in HNSCC cell models. Adaptive EMT-associated resistance appears to be of clinical importance, and we provide effective molecular approaches to exploit this therapeutically.

Radiomics Features on Enhanced Computed Tomography Predict FOXP3 Expression and Clinical Prognosis in Patients with Head and Neck Squamous Cell Carcinoma

Forkhead box P3 (FOXP3) has been identified as a novel molecular marker in various types of cancer. The present study assessed the expression of FOXP3 in patients with head and neck squamous cell carcinoma (HNSCC) and its potential as a clinical prognostic indicator, and developed a radiomics model based on enhanced computed tomography (CT) imaging. Data from 483 patients with HNSCC were downloaded from the Cancer Genome Atlas for FOXP3 prognostic analysis and enhanced CT images from 139 patients included in the Cancer Imaging Archives, which were subjected to the maximum relevance and minimum redundancy and recursive feature elimination algorithms for radiomics feature extraction and processing. Logistic regression was used to build a model for predicting FOXP3 expression. A prognostic scoring system for radiomics score (RS), FOXP3, and patient clinicopathological factors was established to predict patient survival. The area under the receiver operating characteristic (ROC) curve (AUC) and calibration curve and decision curve analysis (DCA) were used to evaluate model performance. Furthermore, the relationship between FOXP3 and the immune microenvironment, as well as the association between RS and immune checkpoint-related genes, was analyzed. Results of analysis revealed that patients with HNSCC and high FOXP3 mRNA expression exhibited better overall survival. Immune infiltration analysis revealed that FOXP3 had a positive correlation with CD4 + and CD8 + T cells and other immune cells. The 8 best radiomics features were selected to construct the radiomics model. In the FOXP3 expression prediction model, the AUC values were 0.707 and 0.702 for the training and validation sets, respectively. Additionally, the calibration curve and DCA demonstrated the positive diagnostic utility of the model. RS was correlated with immune checkpoint-related genes such as ICOS, CTLA4, and PDCD1. A predictive nomogram was established, the AUCs were 0.87, 0.787, and 0.801 at 12, 24, and 36 months, respectively, and DCA demonstrated the high clinical applicability of the nomogram. The enhanced CT radiomics model can predict expression of FOXP3 and prognosis in patients with HNSCC. As such, FOXP3 may be used as a novel prognostic marker to improve individualized clinical diagnosis and treatment decisions.

Homeobox-D 1 and FTO form a transcriptional-epigenetic feedback loop to promote head and neck cancer proliferation

Transcription factors (TFs) and N6-methyladenosine (m6A) modifiers are involved in tumor progression through transcriptional regulation and posttranscriptional regulation of genes, respectively. However, the crosstalk and role of these two types of gene expression regulators in head and neck squamous cell carcinoma (HNSC) remains poorly understood. In this study, we demonstrate that the TF homeobox-D1 (HOXD1) and the m6A demethylase fat mass and obesity-associated protein (FTO) form a positive feedback loop to promote cell proliferation and survival in HNSC. Clinically, HOXD1 expression is dysregulated in multiple cancer types and is associated with worse prognosis in patients with HNSC, stomach adenocarcinoma, uterine corpus endometrial carcinoma, and pheochromocytoma and paraganglioma. Mechanistically, FTO is overexpressed in HNSC tumor samples and positively regulates HOXD1 expression in an m6A-dependent manner. Functionally, deficiency of HOXD1 relieved the resistance of HNSC cells to apoptosis and arrested tumor cells at the G0/G1 phase, thereby inhibiting cell growth, whereas overexpression of HOXD1 caused the opposite effect. Furthermore, HOXD1 activates the transcription of the oncogenic factor FTO by directly targeting its promoter. Downregulation of FTO mimicked the biological effect of HOXD1 knockdown on HNSC. Importantly, overexpression of HOXD1 significantly rescued the proliferation inhibition and apoptosis promotion of HNSC cells induced by deficiency of FTO. Together, our findings reveal HOXD1 as a novel prognostic predictor and a potential target for HNSC, providing mechanistic insights into the role of the HOXD1-FTO circuit in this cancer.

Clinical implication of genetic intratumor heterogeneity for targeted therapy in head and neck cancer

BACKGROUND

Genomic profiling is increasingly used both in therapeutic decision-making and as inclusion criteria for trials testing targeted therapies. However, the mutational landscape may vary across different areas of a tumor and intratumor heterogeneity will challenge treatments or clinical decisions based on single tumor biopsies. The purpose of this study was to assess the clinical relevance of genetic intratumor heterogeneity in head and neck squamous cell carcinomas (HNSCC) using the ESMO Scale for Clinical Actionability of Molecular Targets (ESCAT).

MATERIALS AND METHODS

This prospective study included 33 whole tumor specimens from 28 patients with primary or recurrent HNSCC referred for surgery. Three tumor blocks were selected from central, semi-peripheral, and peripheral positions, mimicking biopsies in three different locations. Genetic analysis of somatic copy number alterations (SCNAs) was performed on the three biopsies using Oncoscan, focusing on 45 preselected HNSCC genes of interest. Clinical relevance was assessed using the ESCAT score to investigate whether and how treatment decisions would change based on the three biopsies from the same tumor.

RESULTS

The SCNAs identified among 45 preselected genes within the three tumor biopsies derived from the same tumor revealed distinct variations. The detected discrepancies could potentially influence treatment approaches or clinical decisions in 36% of the patients if only one tumor biopsy was used. Recurrent tumors exhibited significantly higher variation in SCNAs than primary tumors (p = .024). No significant correlation between tumor size and heterogeneity (p = .7) was observed.

CONCLUSION

In 36% of patients diagnosed with HNSCC, clinically significant intratumor heterogeneity was observed which may have implications for patient management. This finding substantiates the need for future studies that specifically investigate the clinical implications associated with intratumor heterogeneity.

A comprehensive prognostic score for head and neck squamous cancer driver genes and phenotype traits

BACKGROUND

Head and neck squamous cancer (HNSCC) presents variable phenotype and progression features. Clinically applicable, high-accuracy multifactorial prognostic models for HNSCC survival outcomes are warranted and an active area of research. This study aimed to construct a comprehensive prognostic tool for HNSCC overall survival by integrating cancer driver genes with tumor clinical and phenotype information.

METHODS

Key overall survival-related cancer driver genes were screened from among main effector and reciprocal gene pairs using TCGA data using univariate Cox proportional hazard regression analysis. Independent validation was performed using the GSE41613 dataset. The main effector genes among these were selected using LASSO regression and transcriptome score modeling was performed using multivariate Cox regression followed by validation analysis of the prognostic score. Next, multivariate Cox regression analysis was performed using the transcriptome score combined with age, grade, gender, and stage. An 'Accurate Prediction Model of HNSCC Overall Survival Score' (APMHO) was computed and validated. Enriched functional pathways, gene mutational landscape, immune cell infiltration, and immunotherapy sensitivity markers associated with high and low APMHO scores were analyzed.

RESULTS

Screening 107 overall survival-related cancer genes and 402 interacting gene pairs, 6 genes: CRLF2, HSP90AA1, MAP2K1, PAFAH1B2, MYCL and SET genes, were identified and a transcriptional score was obtained. Age, stage and transcriptional score were found to be significant predictors in Cox regression analysis and used to construct a final APMHO model showing an AUC > 0.65 and validated. Transcriptional score, age, pathologic_N, pathologic_T, stage, and TCGA_subtype were significantly different in distribution between high and low APMHO groups. High APMHO samples showed significantly higher mutation rate, enriched tumor-related pathways including Hypoxia, unfold_protein_response, Glycolysis, and mTORC1 signaling, along with differences in immune cell infiltration and immune checkpoint, interferon-γ pathway and m6A regulator expression patterns.

CONCLUSION

The APMHO score combining transcriptional and clinical variables showed good prognostic ability for HNSCC overall survival outcomes and was associated with different patterns of phenotypical features, immune and mutational landscape, and immunotherapy sensitivity marker expression. Future studies should validate this score in independent clinical cohorts.

Enhancing head and neck tumor management with artificial intelligence: Integration and perspectives

Head and neck tumors (HNTs) constitute a multifaceted ensemble of pathologies that primarily involve regions such as the oral cavity, pharynx, and nasal cavity. The intricate anatomical structure of these regions poses considerable challenges to efficacious treatment strategies. Despite the availability of myriad treatment modalities, the overall therapeutic efficacy for HNTs continues to remain subdued. In recent years, the deployment of artificial intelligence (AI) in healthcare practices has garnered noteworthy attention. AI modalities, inclusive of machine learning (ML), neural networks (NNs), and deep learning (DL), when amalgamated into the holistic management of HNTs, promise to augment the precision, safety, and efficacy of treatment regimens. The integration of AI within HNT management is intricately intertwined with domains such as medical imaging, bioinformatics, and medical robotics. This article intends to scrutinize the cutting-edge advancements and prospective applications of AI in the realm of HNTs, elucidating AI's indispensable role in prevention, diagnosis, treatment, prognostication, research, and inter-sectoral integration. The overarching objective is to stimulate scholarly discourse and invigorate insights among medical practitioners and researchers to propel further exploration, thereby facilitating superior therapeutic alternatives for patients.

Constructing heterogeneous single-cell landscape and identifying microenvironment molecular characteristics of primary and lymphatic metastatic head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) accounts for more than half of head and neck tumors, roughly 90%. This study focused on constructing the heterogeneous landscape using single-cell and bulk transcriptomic data to identify molecular characteristics of the microenvironment in primary and lymphatic metastatic head and neck squamous carcinomas.

METHOD

The study enrolled 23 HNSCC samples with scRNA-seq data and 546 HNSCC samples from TCGA. We used Monocle to sort the cells and used CellPhoneDB to explore the cell-cell interactions. Infercnv, which was used to infer cells with apparent copy number variation based on single-cell sequencing transcriptome data. We re-evaluated HNSCC bulk RNA transcriptome data to characterize the functions of different cell types in shaping the immune microenvironment of HNSCC.

RESULTS

We combined genealogical reconstruction, CNV inference, and cellular interactions to uncover the characteristics of distinct cell populations in different disease states, differences in cancer and immune cell lineages of differentiation trajectories, and interactions between non-immune and immune cell. PD-1 and PD-L1/PD-L2 expressed extremely rare in T cells, the immune checkpoint molecule KLRB1-CLEC2D achieved a high-level expression. We identified three microenvironment-based HNSCC subtypes associated with the prognosis of HNSCC patients.

CONCLUSIONS

In summary, the present study dissected the intratumoral heterogeneity and immune microenvironment of primary and metastatic HNSCC, which is crucial to reveal the mechanisms of resistance to immunotherapy in HNSCC in different disease states and is expected to assist in the further investigation of the mechanism of HNSCC cell metastasis and guide the treatment of clinical patients.

Predictive value of tumor budding in head and neck squamous cell carcinoma: an update

Head and neck squamous cell carcinoma forms an anatomically and functionally complex group of malignancies. The significant local aggressiveness and frequent regional relapses motivate ongoing research to identify more reliable and sensitive prognostic and predictive biomarkers. One emerging area of cancer biology is the evaluation of tumor budding at the advancing invasive front of various types of epithelial cancers. Recent studies suggest that tumor budding is a relatively common phenomenon in cancer progression and that it may have important prognostic implications for patients due to its potential to provide valuable insights into the biology and clinical behavior of head and neck cancer. In this review, we aim to provide information about tumor budding in head and neck squamous cell carcinoma. Thus, we hope to shed light on the complex biology of these malignancies, as well as aiding diagnostic, classification, and better characterization and thereby, looking for new avenues for improving patient outcomes.

Identification of CKS2 as a novel prognostic biomarker and potential therapeutic target for oral squamous cell carcinoma

Background

The cyclin-dependent kinase subunit 2 (CKS2) is recognized to have a substantial impact on the pathogenesis and advancement of several malignant neoplasms. Nevertheless, its biological function and prognostic significance in oral squamous cell carcinoma (OSCC) have yet to be thoroughly investigated. Our primary objective was to clarify the contribution of CKS2 in the progression and prognosis of OSCC.

Methods

We first conducted a thorough examination of online databases to investigate the expression of CKS2, and subsequently corroborated our discoveries by analyzing clinical specimens that we collected. According to the clinicopathological data, we then explored the prognostic significance of CKS2. Furthermore, we predicted the role of CKS2 in OSCC progression by employing weighted gene co-expression network analysis (WGCNA) in conjunction with functional enrichment analysis. We conducted functional experiments in vitro to confirm our speculations. Additionally, we explored other potential functions of CKS2 in immune infiltration, tumor mutation burden (TMB), and drug sensitivity. Finally, we established and validated a nomogram that effectively integrated CKS2-related genes and other relevant clinical factors.

Results

Our findings indicated a significant upregulation of CKS2 expression in OSCC tissues compared to normal groups, which was positively associated with poor clinical outcomes. We also predicted and validated the role of CKS2 in promoting proliferation by regulating the cell cycle. Additionally, its upregulation was significantly correlated to enhanced immune cell infiltration, high TMB, and increased sensitivity of anti-tumor agents. Following verification, the nomogram was conducted to quantify an individual's survival probability.

Conclusions

In general, our study indicates that CKS2 is a novel prognostic biomarker and potential therapeutic target in OSCC.

An integrated genomic approach identifies follistatin as a target of the p63-epidermal growth factor receptor oncogenic network in head and neck squamous cell carcinoma

Although numerous putative oncogenes have been associated with the etiology of head and neck squamous cell carcinoma (HNSCC), the mechanisms by which these oncogenes and their downstream targets mediate tumor progression have not been fully elucidated. We performed an integrative analysis to identify a crucial set of targets of the oncogenic transcription factor p63 that are common across multiple transcriptomic datasets obtained from HNSCC patients, and representative cell line models. Notably, our analysis revealed FST which encodes follistatin, a secreted glycoprotein that inhibits the transforming growth factor TGFβ/activin signaling pathways, to be a direct transcriptional target of p63. In addition, we found that FST expression is also driven by epidermal growth factor receptor EGFR signaling, thus mediating a functional link between the TGF-β and EGFR pathways. We show through loss- and gain-of-function studies that FST predominantly imparts a tumor-growth and migratory phenotype in HNSCC cells. Furthermore, analysis of single-cell RNA sequencing data from HNSCC patients unveiled cancer cells as the dominant source of FST within the tumor microenvironment and exposed a correlation between the expression of FST and its regulators with immune infiltrates. We propose FST as a prognostic biomarker for patient survival and a compelling candidate mediating the broad effects of p63 on the tumor and its associated microenvironment.

5'-Ectonucleotidase CD73/NT5E supports EGFR-mediated invasion of HPV-negative head and neck carcinoma cells

BACKGROUND

Epithelial-to-mesenchymal transition (EMT) of malignant cells is a driving force of disease progression in human papillomavirus-negative (HPV-negative) head and neck squamous cell carcinomas (HNSCC). Sustained hyper-activation of epidermal growth factor receptor (EGFR) induces an invasion-promoting subtype of EMT (EGFR-EMT) characterized by a gene signature ("'EGFR-EMT_Signature'") comprising 5´-ectonucleotidase CD73. Generally, CD73 promotes immune evasion via adenosine (ADO) formation and associates with EMT and metastases. However, CD73 regulation through EGFR signaling remains under-explored and targeting options are amiss.

METHODS

CD73 functions in EGFR-mediated tumor cell dissemination were addressed in 2D and 3D cellular models of migration and invasion. The novel antagonizing antibody 22E6 and therapeutic antibody Cetuximab served as inhibitors of CD73 and EGFR, respectively, in combinatorial treatment. Specificity for CD73 and its role as effector or regulator of EGFR-EMT were assessed upon CD73 knock-down and over-expression. CD73 correlation to tumor budding was studied in an in-house primary HNSCC cohort. Expression correlations, and prognostic and predictive values were analyzed using machine learning-based algorithms and Kaplan-Meier survival curves in single cell and bulk RNA sequencing datasets.

RESULTS

CD73/NT5E is induced by the EGF/EGFR-EMT-axis and blocked by Cetuximab and MEK inhibitor. Inhibition of CD73 with the novel antagonizing antibody 22E6 specifically repressed EGFR-dependent migration and invasion of HNSCC cells in 2D. Cetuximab and 22E6 alone reduced local invasion in a 3D-model. Interestingly, combining inefficient low-dose concentrations of Cetuximab and 22E6 revealed highly potent in invasion inhibition, substantially reducing the functional IC50 of Cetuximab regarding local invasion. A role for CD73 as an effector of EGFR-EMT in local invasion was further supported by knock-down and over-expression experiments in vitro and by high expression in malignant cells budding from primary tumors. CD73 expression correlated with EGFR pathway activity, EMT, and partial EMT (p-EMT) in malignant single HNSCC cells and in large patient cohorts. Contrary to published data, CD73 was not a prognostic marker of overall survival (OS) in the TCGA-HNSCC cohort when patients were stratified for HPV-status. However, CD73 prognosticated OS of oral cavity carcinomas. Furthermore, CD73 expression levels correlated with response to Cetuximab in HPV-negative advanced, metastasized HNSCC patients.

CONCLUSIONS

In sum, CD73 is an effector of EGF/EGFR-mediated local invasion and a potential therapeutic target and candidate predictive marker for advanced HPV-negative HNSCC.

Tanshinone IIA inhibits cell viability and promotes PUMA-mediated apoptosis of oral squamous cell carcinoma

Apoptosis alteration is responsible for tumorigenesis and tumor resistance to therapies. The natural product Tanshinone IIA (Tan IIA) exhibits potent inhibitory effects against various tumors. However, the effect of Tan IIA on apoptosis and its underlying mechanism remains elusive in oral squamous cell carcinoma (OSCC). Here, we demonstrated that Tan IIA dose-dependently suppressed cell viability and colony formation in CAL27, SCC4, and SCC25 cells. Moreover, Tan IIA inhibited Akt activation from inducing Foxo3a dephosphorylation and PUMA-mediated apoptosis. PUMA or Foxo3a knockdown compromised the inhibitory effect of Tan IIA on OSCC cells. Tan IIA administration inhibited CAL27-deprived xenograft tumor growth and increased PUMA expression in vivo. Tan IIA synergistically intensified the efficacy of CDDP/5-FU-based chemotherapy on OSCC cells. Overall, our results revealed that Tan IIA exerted potent antitumor effects via promoting PUMA-mediated apoptosis in OSCC cells.

Cellular hierarchy framework based on single-cell/multi-patient sample sequencing reveals metabolic biomarker PYGL as a therapeutic target for HNSCC

BACKGROUND

A growing body of research has revealed the connection of metabolism reprogramming and tumor progression, yet how metabolism reprogramming affects inter-patient heterogeneity and prognosis in head and neck squamous cell carcinoma (HNSCC) still requires further explorations.

METHODS

A cellular hierarchy framework based on metabolic properties discrepancy, METArisk, was introduced to re-analyze the cellular composition from bulk transcriptomes of 486 patients through deconvolution utilizing single-cell reference profiles from 25 primary and 8 metastatic HNSCC sample integration of previous studies. Machine learning methods were used to identify the correlations between metabolism-related biomarkers and prognosis. The functions of the genes screened out in tumor progression, metastasis and chemotherapy resistance were validated in vitro by cellular functional experiments and in vivo by xenograft tumor mouse model.

RESULTS

Incorporating the cellular hierarchy composition and clinical properties, the METArisk phenotype divided multi-patient cohort into two classes, wherein poor prognosis of METArisk-high subgroup was associated with a particular cluster of malignant cells with significant activity of metabolism reprogramming enriched in metastatic single-cell samples. Subsequent analysis targeted for phenotype differences between the METArisk subgroups identified PYGL as a key metabolism-related biomarker that enhances malignancy and chemotherapy resistance by GSH/ROS/p53 pathway, leading to poor prognosis of HNSCC.

CONCLUSION

PYGL was identified as a metabolism-related oncogenic biomarker that promotes HNSCC progression, metastasis and chemotherapy resistance though GSH/ROS/p53 pathway. Our study revealed the cellular hierarchy composition of HNSCC from the cell metabolism reprogramming perspective and may provide new inspirations and therapeutic targets for HNSCC in the future.

A ferroptosis-related prognostic model with excellent clinical performance based on the exploration of the mechanism of oral squamous cell carcinoma progression

As a hot topic today, ferroptosis is closely involved in the progression and treatment of cancer. Accordingly, we built a prognostic model around ferroptosis to predict the overall survival of OSCC patients. We used up to 6 datasets from 3 different databases to ensure the credibility of the model. Then, through differentially expressed, Univariate Cox, and Lasso regression analyses, a model composed of nine prognostic-related differently expressed ferroptosis-related genes (CISD2, DDIT4, CA9, ALOX15, ATG5, BECN1, BNIP3, PRDX5 and MAP1LC3A) were constructed. Moreover, Kaplan-Meier curves, Receiver Operating Characteristic curves and principal component analysis used to verify the model's predictive ability showed the model's superiority. To deeply understand the mechanism of ferroptosis affecting the occurrence, development and prognosis of OSCC, we performed enrichment analysis in different risk groups identified by the model. The results showed that numerous TP53-related, immune-related and ferroptosis-related functions and pathways were enriched. Further immune microenvironment analysis and mutation analysis have once again revealed the correlation between risk score and immunity and TP53 mutation. Finally, the correlation between risk score and OSCC clinical treatment, as well as Nomogram show the brilliant clinical application prospects of the prognostic model.

Transcriptomic and Proteomic Profiles for Elucidating Cisplatin Resistance in Head-and-Neck Squamous Cell Carcinoma

To identify the novel genes involved in chemoresistance in head and neck squamous cell carcinoma (HNSCC), we explored the expression profiles of the following cisplatin (CDDP) resistant (R) versus parental (sensitive) cell lines by RNA-sequencing (RNA-seq): JHU029, HTB-43 and CCL-138. Using the parental condition as a control, 30 upregulated and 85 downregulated genes were identified for JHU029-R cells; 263 upregulated and 392 downregulated genes for HTB-43-R cells, and 154 upregulated and 68 downregulated genes for CCL-138-R cells. Moreover, we crossed-checked the RNA-seq results with the proteomic profiles of HTB-43-R (versus HTB-43) and CCL-138-R (versus CCL-138) cell lines. For the HTB-43-R cells, 21 upregulated and 72 downregulated targets overlapped between the proteomic and transcriptomic data; whereas in CCL-138-R cells, four upregulated and three downregulated targets matched. Following an extensive literature search, six genes from the RNA-seq (CLDN1, MAGEB2, CD24, CEACAM6, IL1B and ISG15) and six genes from the RNA-seq and proteomics crossover (AKR1C3, TNFAIP2, RAB7A, LGALS3BP, PSCA and SSRP1) were selected to be studied by qRT-PCR in 11 HNSCC patients: six resistant and five sensitive to conventional therapy. Interestingly, the high MAGEB2 expression was associated with resistant tumours and is revealed as a novel target to sensitise resistant cells to therapy in HNSCC patients.

Effect of the m6ARNA gene on the prognosis of thyroid cancer, immune infiltration, and promising immunotherapy

BACKGROUND

Accumulating evidence suggests that N6-methyladenosine (m6A) RNA methylation plays an important role in tumor proliferation and growth. However, its effect on the clinical prognosis, immune infiltration, and immunotherapy response of thyroid cancer patients has not been investigated in detail.

METHODS

Clinical data and RNA expression profiles of thyroid cancer were extracted from the Cancer Genome Atlas-thyroid carcinoma (TCGA-THCA) and preprocessed for consensus clustering. The risk model was constructed based on differentially expressed genes (DEGs) using Least Absolute Shrinkage and Selection Operator (LASSO) and Cox regression analyses. The associations between risk score and clinical traits, immune infiltration, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Set Enrichment Analysis (GSEA), immune infiltration, and immunotherapy were assessed. Immunohistochemistry was used to substantiate the clinical traits of our samples.

RESULTS

Gene expression analysis showed that 17 genes, except YHTDF2, had significant differences (vs healthy control, P<0.001). Consensus clustering yielded 2 clusters according to their clinical features and estimated a poorer prognosis for Cluster 1 (P=0.03). The heatmap between the 2 clusters showed differences in T (P<0.01), N (P<0.001) and stage (P<0.01). Based on univariate Cox and LASSO regression, a risk model consisting of three high-risk genes (KIAA1429, RBM15, FTO) was established, and the expression difference between normal and tumor tissues of three genes was confirmed by immunohistochemical results of our clinical tissues. KEGG and GSEA analyses showed that the risk DEGs were related mainly to proteolysis, immune response, and cancer pathways. The levels of immune infiltration in the high- and low-risk groups were different mainly in iDCs (P<0.05), NK cells (P<0.05), and type-INF-II (P<0.001). Immunotherapy analysis yielded 30 drugs associated with the expression of each gene and 20 drugs associated with the risk score.

CONCLUSIONS

Our risk model can act as an independent marker for thyroid cancer and provides promising immunotherapy targets for its treatment.