Single-cell transcriptomics reveal the intratumoral landscape of infiltrated T-cell subpopulations in oral squamous cell carcinoma

T cell proliferation-related subtypes, prognosis model and characterization of tumor microenvironment in head and neck squamous cell carcinoma

Background

Thirty-three synthetic driver genes of T-cell proliferation have recently been identified through genome-scale screening. However, the tumor microenvironment (TME) cell infiltration, prognosis, and response to immunotherapy mediated by multiple T cell proliferation-related genes (TRGs) in patients with head and neck squamous cell carcinoma (HNSC) remain unclear.

Methods

This study examined the genetic and transcriptional changes in 771 patients with HNSC by analyzing the TRGs from two independent datasets. Two different subtypes were analyzed to investigate their relationship with immune infiltrating cells in the TME and patient prognosis. The study also developed and validated a risk score to predict overall survival (OS). Furthermore, to enhance the clinical utility of the risk score, an accurate nomogram was constructed by combining the characteristics of this study.

Results

The low-risk score observed in this study was associated with high levels of immune checkpoint expression and TME immune activation, indicating a favorable OS outcome. Additionally, various factors related to risk scores were depicted.

Conclusion

Through comprehensive analysis of TRGs in HNSC, our study has revealed the characteristics of the TME and prognosis, providing a basis for further investigation into TRGs and the development of more effective immunotherapy and targeted therapy strategies.

Unravelling immune microenvironment features underlying tumor progression in the single-cell era

The relationship between the immune cell and tumor occurrence and progression remains unclear. Profiling alterations in the tumor immune microenvironment (TIME) at high resolution is crucial to identify factors influencing cancer progression and enhance the effectiveness of immunotherapy. However, traditional sequencing methods, including bulk RNA sequencing, exhibit varying degrees of masking the cellular heterogeneity and immunophenotypic changes observed in early and late-stage tumors. Single-cell RNA sequencing (scRNA-seq) has provided significant and precise TIME landscapes. Consequently, this review has highlighted TIME cellular and molecular changes in tumorigenesis and progression elucidated through recent scRNA-seq studies. Specifically, we have summarized the cellular heterogeneity of TIME at different stages, including early, late, and metastatic stages. Moreover, we have outlined the related variations that may promote tumor occurrence and metastasis in the single-cell era. The widespread applications of scRNA-seq in TIME will comprehensively redefine the understanding of tumor biology and furnish more effective immunotherapy strategies.

Spatiotemporal characteristics of tissue derived small extracellular vesicles is associated with tumor relapse and anti-PD-1 response

Small extracellular vesicles (sEVs) residing at tumor tissues are valuable specimens for biopsy. Tumor heterogeneity is common across all cancer types, but the heterogeneity of tumor tissue-derived sEVs (Ti-sEVs) is undefined. This study aims to discover the spatial distributions of Ti-sEVs in oral squamous cell carcinoma (OSCC) tissues and explore how these vesicle distributions affect the patients' prognosis. Multi-regional sampling enabled us to uncover that Ti-sEVs' accumulation at peritumoral sites correlates with a higher disease-free survival rate, and conversely, sparse peritumoral Ti-sEVs tend to forecast a higher risk of relapse. Of those relapsed patients, Ti-sEVs strongly bind to extracellular matrix and subsequently degrade it for allowing themselves enter the bloodstream rather than staying in situ. In advanced OSCC patients, the quantity and spatial distribution of Ti-sEVs prior to anti-PD-1 treatment, as well as the temporal variance of Ti-sEVs before and after immunotherapy, strongly map the clinical response and can help to distinguish the patients with shrinking tumors from those with growing tumors. Our work elucidates the correlation of spatiotemporal features of Ti-sEVs with patients' therapeutic outcomes and exhibit the potential for using Ti-sEVs as a predictor to forecast prognosis and screen the responders to anti-PD-1 therapy.

Tumor-infiltrating immune cells and survival in head and neck squamous cell carcinoma: a retrospective computational study

The immune infiltration profiles of the tumor microenvironment have effects on the prognosis of head and neck squamous cell carcinoma (HNSCC). Whereas, HNSCC is a heterogeneous group of tumors, but past work has not taken this into consideration. Herein, we investigate the associations between survival and the function of immune cells in different tumorigenic sites of HNSCC. 1149 samples of HNSCC were collected from publicly accessible databases. Based on gene expression data, CIBERSORTx was applied to determine the proportion of 22 immune cell subpopulations. In the Cox regression model, the associations between overall survival, disease-free survival, and immune cells were examined, modeling gene expression and immune cell proportion as quartiles. Consensus cluster analysis was utilized to uncover immune infiltration profiles. Regardless of tumor sites, CD8+ T cells and activated CD4 memory T cells were associated with favorable survival, while eosinophils were the opposite. The survival of the hypopharynx, oral cavity, and larynx subsites was somewhat affected by immune cells, while the survival of the oropharynx subsite potentially was the most impacted. High expression of TIGIT, CIITA, and CXCR6 was linked to better survival, mainly in the oropharynx subsite. Immune cell clusters with four distinct survival profiles were discovered, of which the cluster with a high CD8+ T cell content had a better prognosis. The immune-infiltration pattern is related to the survival of HNSCC to varying degrees depending on the tumor sites; forthcoming studies into immune-mediated infiltration profiles will lay the groundwork for treating HNSCC with precision therapy.

An AI-assisted integrated, scalable, single-cell phenomic-transcriptomic platform to elucidate intratumor heterogeneity against immune response

We present a novel framework combining single-cell phenotypic data with single-cell transcriptomic analysis to identify factors underpinning heterogeneity in antitumor immune response. We developed a pairwise, tumor-immune discretized interaction assay between natural killer (NK-92MI) cells and patient-derived head and neck squamous cell carcinoma (HNSCC) cell lines on a microfluidic cell-trapping platform. Furthermore we generated a deep-learning computer vision algorithm that is capable of automating the acquisition and analysis of a large, live-cell imaging data set (>1 million) of paired tumor-immune interactions spanning a time course of 24 h across multiple HNSCC lines (n = 10). Finally, we combined the response data measured by Kaplan-Meier survival analysis against NK-mediated killing with downstream single-cell transcriptomic analysis to interrogate molecular signatures associated with NK-effector response. As proof-of-concept for the proposed framework, we efficiently identified MHC class I-driven cytotoxic resistance as a key mechanism for immune evasion in nonresponders, while enhanced expression of cell adhesion molecules was found to be correlated with sensitivity against NK-mediated cytotoxicity. We conclude that this integrated, data-driven phenotypic approach holds tremendous promise in advancing the rapid identification of new mechanisms and therapeutic targets related to immune evasion and response.

CCDC71L as a novel prognostic marker and immunotherapy target via lipid metabolism in head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is the sixth most widespread cancer globally with high rate and poor prognosis. Coiled-coil domain containing 71 like (CCDC71L) exerts an important role in cellular lipid metabolic process. However, its function in HNSCC remains unclear. To this end, we examined the CCDC71L implications for prognosis and tumor microenvironment in HNSCC.

MATERIALS AND METHODS

First, CCDC71L expression was explored through The Cancer Genome Atlas (TCGA), Human Protein Atlas (HPA), the Gene Expression Profiling Interactive Analysis 2 (GEPIA2), Clinical Proteomic Tumor Analysis Consortium (CPTAC) and the Gene Expression Omnibus (GEO) databases. The clinicpathological information were obtained from the dataset of TCGA. The Kaplan-Meier Plotter databases and Cox model were performed for the determination of prognostic values of CCDC71L, including the overall survival (OS), progress free interval (PFI), recurrence-free survival (RFS) and disease specific survival (DSS). Then, the potential mechanism of CCDC71L in HNSCC development was elucidated by means of Gene set enrichment analysis (GSEA) and Metascape databases. Furthermore, the relevance of CCDC71L to immune cells infiltration and immune checkpoints was assessed. The correlations among CCDC71L expression, mutational landscape and genome heterogeneity [mutant-allele tumor heterogeneity (MATH) and tumor purity] were detected by the data in TCGA.

RESULTS

CCDC71L expression was significantly upregulated in HNSCC, and positively associated with age, gender and N stage. Higher CCDC71L expression resulted in poor OS, RFS, DSS and PFI. Multivariate Cox regression analysis showed CCDC71L would be an independent prognostic marker in patients with HNSCC. Moreover, CCDC71L and the level of macrophage and neutrophil cells infiltration were significantly correlated in HNSCC. High expression of CCDC71L was related to immune checkpoint genes, oncogene mutations and genome heterogeneity markers.

CONCLUSION

These results implied that CCDC71L plays vital roles in HNSCC progression, which could be used as a underlying biomarker for the diagnosis and prognosis of HNSCC. Meanwhile, CCDC71L participates in immune regulation, which has a potential value for the immunotherapy of HNSCC patients.

TLR9 activation induces immunosuppression and tumorigenesis via PARP1/PD-L1 signaling pathway in oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is the most common type of oral cancer, and metastasis and immunosuppression are responsible for the poor prognosis of OSCC. Previous studies have shown that poly(ADP-ribose) polymerase (PARP)1 plays a key role in the pathogenesis of OSCC. Therefore, PARP1 may serve as an important research target for the potential treatment of OSCC. Here, we aimed to investigate the role of PARP1 in the tumorigenesis of OSCC and elucidate the key molecular mechanisms of its upstream and downstream regulation in vivo and in vitro. In human OSCC tissues and cells, Toll-like receptor (TLR)9 and PD-L1 were highly expressed and PARP1 was lowly expressed. Suppression of TLR9 remarkably repressed CAL27 and SCC9 cell proliferation, migration, and invasion. After coculture, we found that low expression of TLR9 inhibited PD-L1 expression and immune escape. In addition, TLR9 regulated PD-L1 expression through the PARP1/STAT3 pathway. PARP1 mediated the effects of TLR9 on OSCC cell proliferation, migration, and invasion and immune escape. Additionally, in vivo experiments further verified that TLR9 promoted tumor growth and immune escape by inhibiting PARP1. Collectively, TLR9 activation induced immunosuppression and tumorigenesis via PARP1/PD-L1 signaling pathway in OSCC, providing important insights for subsequent in-depth exploration of the mechanism of OSCC.NEW & NOTEWORTHY In this research, we took PARP1 as the key target to explore its regulatory effect on oral squamous cell carcinoma (OSCC). The key molecular mechanisms involved in its upstream and downstream regulation were elucidated in OSCC cell lines in vitro and tumor-bearing mice in vivo, combined with clinical OSCC tissues.

New perspectives on biology, disease progression, and therapy response of head and neck cancer gained from single cell RNA sequencing and spatial transcriptomics

Head and neck squamous cell carcinoma (HNSCC) is one of the most frequent cancers worldwide. The main risk factors are consumption of tobacco products and alcohol, as well as infection with human papilloma virus. Approved therapeutic options comprise surgery, radiation, chemotherapy, targeted therapy through epidermal growth factor receptor inhibition, and immunotherapy, but outcome has remained unsatisfactory due to recurrence rates of ~50% and the frequent occurrence of second primaries. The availability of the human genome sequence at the beginning of the millennium heralded the omics era, in which rapid technological progress has advanced our knowledge of the molecular biology of malignant diseases, including HNSCC, at an unprecedented pace. Initially, microarray-based methods, followed by approaches based on next-generation sequencing, were applied to study the genetics, epigenetics, and gene expression patterns of bulk tumors. More recently, the advent of single-cell RNA sequencing (scRNAseq) and spatial transcriptomics methods has facilitated the investigation of the heterogeneity between and within different cell populations in the tumor microenvironment (e.g., cancer cells, fibroblasts, immune cells, endothelial cells), led to the discovery of novel cell types, and advanced the discovery of cell-cell communication within tumors. This review provides an overview of scRNAseq, spatial transcriptomics, and the associated bioinformatics methods, and summarizes how their application has promoted our understanding of the emergence, composition, progression, and therapy responsiveness of, and intercellular signaling within, HNSCC.

Application of single-cell RNA sequencing in head and neck squamous cell carcinoma

Single-cell RNA sequencing has been broadly applied to head and neck squamous cell carcinoma (HNSCC) for characterizing the heterogeneity and genomic mutations of HNSCC benefiting from the advantage of single-cell resolution. We summarized most of the current studies and aimed to explore their research methods and ideas, as well as how to transform them into clinical applications. Through single-cell RNA sequencing, we found the differences in tumor cells' expression programs and differentiation tracks. The studies of immune microenvironment allowed us to distinguish immune cell subpopulations, the extensive expression of immune checkpoints, and the complex crosstalk network between immune cells and non-immune cells. For cancer-associated fibroblasts (CAFs), single-cell RNA sequencing had made an irreplaceable contribution to the exploration of their differentiation status, specific CAFs markers, and the interaction with tumor cells and immune cells. In addition, we demonstrated in detail how single-cell RNA sequencing explored the HNSCC epithelial-to-mesenchymal transition (EMT) model and the mechanism of drug resistance, as well as its clinical value.

Prognostic model based on telomere-related genes predicts the risk of oral squamous cell carcinoma

BACKGROUND

This study investigated a potential prognostic model based on telomere-related genes (TRGs) for the clinical prediction of oral squamous cell carcinoma (OSCC).

METHODS

Gene expression data and associated clinical phenotypes were obtained from online databases. Differentially expressed (DE)-TRGs were identified between OSCC and normal samples, followed by protein-protein interaction and enrichment analyses. Subsequently, the prognostic genes explored based on the DE-TRGs and survival data were applied in the establishment of the current prognostic model, and an integrated analysis was performed between high- and low-risk groups using a prognostic model. The expression of certain prognostic genes identified in the present study was validated using qPCR analysis and/or western blot in OSCC cell lines and clinical samples.

RESULTS

169 DE-TRGs were identified between the OSCC samples and controls. DE-TRGs are mainly involved in functions such as hypoxia response and pathways such as the cell cycle. Eight TRGs (CCNB1, PDK4, PLOD2, RACGAP1, MET, PLK1, KPNA2, and CCNA2) associated with OSCC survival and prognosis were used to construct a prognostic model. qPCR analysis and western blot showed that most of the eight prognostic genes were consistent with the current bioinformatics results. Analysis of the high- and low-risk groups for OSCC determined by the prognostic model showed that the current prognostic model was reliable.

CONCLUSIONS

A novel prognostic model for OSCC was constructed by TRGs. PLOD2 and APLK1 may participate in the progression of OSCC via responses to hypoxia and cell cycle pathways, respectively. TRGs, including KPNA2 and CCNA2, may serve as novel prognostic biomarkers for OSCC.

TOX regulates T lymphocytes differentiation and its function in tumor

Thymocyte selection-associated high mobility group box protein (TOX) is expressed differently at all T lymphocytes development stages. Owing to more advanced scientific and technological means, including single-cell sequencing technology, heterogeneity of T lymphocytes and TOX has gradually been revealed. Further exploration of such heterogeneity will help us comprehend the developmental stage and functional characteristics of T lymphocytes in greater detail. Emerging evidence supports its regulation not only in exhausting, but also in activating T lymphocytes, thereby verifying TOX heterogeneity. TOX can be used not only as a latent intervention target for tumor diseases and chronic infections, and a therapeutic strategy for autoimmune diseases, but also as a critical factor predicting the drug response and overall survival of patients with malignant tumors.

T cell effects and mechanisms in immunotherapy of head and neck tumors

Head and neck tumors (HNCs) are a common tumor in otorhinolaryngology head and neck surgery, accounting for 5% of all malignant tumors in the body and are the sixth most common malignant tumor worldwide. In the body, immune cells can recognize, kill, and remove HNCs. T cell-mediated antitumor immune activity is the most important antitumor response in the body. T cells have different effects on tumor cells, among which cytotoxic T cells and helper T cells play a major killing and regulating role. T cells recognize tumor cells, activate themselves, differentiate into effector cells, and activate other mechanisms to induce antitumor effects. In this review, the immune effects and antitumor mechanisms mediated by T cells are systematically described from the perspective of immunology, and the application of new immunotherapy methods related to T cells are discussed, with the objective of providing a theoretical basis for exploring and forming new antitumor treatment strategies. Video Abstract.

Review of biomarkers for response to immunotherapy in HNSCC microenvironment

Head and neck squamous cell carcinoma are one of the most common types of cancer worldwide. Although a variety of treatment methods such as surgery, radiotherapy, chemotherapy, and targeted therapy are widely used in diagnosing and treating HNSCC, the survival prognosis of patients has not been significantly improved in the past decades. As an emerging treatment approach, immunotherapy has shown exciting therapeutic effects in R/M HNSCC. However, the current screening methods are still insufficient, and there is a significant need for reliable predictive biomarkers for personalized clinical management and new therapeutic strategies. This review summarized the application of immunotherapy in HNSCC, comprehensively analyzed the existing bioinformatic studies on immunotherapy in HNSCC, evaluated the current methods of tumor immune heterogeneity and immunotherapy, and aimed to screen molecular markers with potential predictive significance. Among them, PD-1 has obvious predictive relevance as the target of existing immune drugs. Clonal TMB is a potential biomarker for HNSCC immunotherapy. The other molecules, including IFN-γ, CXCL, CTLA-4, MTAP, SFR4/CPXM1/COL5A1, TILs, CAFs, exosomes, and peripheral blood indicators, may have suggestive significance for tumor immune microenvironment and prognosis of immunotherapy.

Immunoexpression of Tumor Infiltrating Lymphocytes(TILS) CD4 + and CD8 + in Oral Squamous Cell Carcinoma (OSCC) in Correlations with Clinicopathological Characteristics and Prognosis

OBJECTIVE

This study aimed to analyze CD4 +and CD8 + TILs in oral squamous cell carcinoma (OSCC) and to correlate it with histologic grade of malignancy and clinicopathologic data.

METHODS

The sample was composed of 43 archived specimens. Clinical features and histological grade of malignancy were obtained. The infiltrating intensity of CD4 +, CD8 positive cells were assessed by immunohistochemistry. One-way ANOVA was used to study the association between CD4 +, CD8 + and the grade of OSCC. The cut-off values of the proposed diagnostic indices were received from calculating the coordinates of the receiver operating characteristic (ROC) curve. For clinicopathologic data Independent-Samples T test, Pearson Correlation Coefficient, Correlation Coefficient were used clinicopathologic characteristics.

RESULTS

CD4 +and CD8 + were observed in all specimens. CD4 + were more frequent in poorly differentiated specimens (74.14) (P= 0.021<0.05). CD8 + were more frequent in well- differentiated  specimens (51.18). None of these correlations were significant (P=0.454>0.05). CD4 +/ CD8 ratio was higher in low-grade specimens (180.28) (P=0.017<0.05). No differences between CD4 +, CD8 +and CD4 +/ CD8 ratio between poorly- differentiated   and moderately- differentiated   groups ROC P value (0.370, 0.248, 0.126) respectively. there is a difference between CD4 +, CD4 +/ CD8ratio between poorly- differentiated   and well- differentiated   groups ROC P value (0.022, 0.341, 0.012) Sensitivity (0.857, 0.882), specificity (0.706, 0.857) respectively. and no differences between CD8 + poorly- differentiated and well- differentiated groups ROC P value (0.341). there is a difference between CD4 + between moderately - differentiated   and well- differentiated groups ROC P value (0.038) Sensitivity (0.368), specificity (0,765). No significant correlation was obtained with clinicopathologic findings of OSCC.

CONCLUSION

CD4 + and CD4 +/ CD8 + ratio are independent prognostic factor in OSCC.

Evaluation of CD4+ tumor-infiltrating lymphocyte association with some clinicopathological indices of oral squamous cell carcinoma

BACKGROUND

The delayed diagnosis of oral squamous cell carcinoma (OSCC) affects therapeutic and prognostic strategies, and provides regional recurrence or distant metastasis. The tumor-infiltrating lymphocytes (TILs) are known as a critical diagnostic biomarker in antitumor immune response. We evaluated the association between CD4⁺ T-lymphocyte marker, some clinicopathological indices, and the impact of TILs on the stage and grade of OSCC.

MATERIALS AND METHODS

In this cross-sectional study, 37 OSCC specimens including 16 early and 21 advanced stages (categorized base-on recent clinical oncology references) and their related healthy surgical margin (as internal control group) were collected. Obtained histochemical data were analyzed by SPSS V.23 software. The expression of CD4⁺ marker in tumor microenvironment (TME) was compared by nonparametric Mann-Whitney and Kruskal-Wallis as well as Fisher's exact tests. P < 0.05 was remarked statistically significant.

RESULTS

The low-grade patients represented more CD4⁺ TIL that was statistically significant (P = 0.011). However, there was no statistically significant difference in CD4⁺ TIL between various stages (P = 0.404), tumor size, and lymph node involvement (P > 0.05). Moreover, there was no significant relation between TIL infiltration, age, and tumor localization (P > 0.05), however CD4⁺ expression in women was more than men (P = 0.008). The CD4⁺ T-lymphocyte infiltration in TME was more significant than healthy surgical margin (P < 0.001). There was a statistically significant difference between healthy surgical margin and different grades and stages of OSCCs that lower grades demonstrated more CD4⁺ TIL infiltration (P < 0.001).

CONCLUSION

The CD4⁺ T-lymphocytes may play important role in differentiation and maturity of epithelial cell, tumorigenesis, and progression of OSCC.

Inhibition of Autophagy Promotes the Anti-Tumor Effect of Metformin in Oral Squamous Cell Carcinoma

Oral Squamous Cell Carcinoma (OSCC) is the most common malignant tumor in the head and neck. Due to its high malignancy and easy recurrence, the five-year survival rate is only 50–60%. Currently, commonly used chemotherapy drugs for OSCC include cisplatin, paclitaxel, and fluorouracil, which are highly cytotoxic and cause drug resistance in patients. Therefore, a safe and effective treatment strategy for OSCC is urgent. To address this issue, our study investigated the anti-tumor activity of metformin (the first-line diabetes drug) in OSCC. We found that metformin could inhibit OSCC cell proliferation by promoting apoptosis and blocking the cell cycle in G1 phase. Additionally, we also found that metformin could induce protective autophagy of OSCC cells. After inhibiting autophagy with hydroxychloroquine (HCQ), the metformin-induced apoptosis was enhanced. In vitro, metformin inhibited the growth of subcutaneous xenograft tumor in nude mice and HCQ enhanced this effect of metformin. Therefore, metformin combined with HCQ may become a safe and effective treatment strategy for OSCC.

Prognostic Value of Molecular Intratumor Heterogeneity in Primary Oral Cancer and Its Lymph Node Metastases Assessed by Mass Spectrometry Imaging

Different aspects of intra-tumor heterogeneity (ITH), which are associated with the development of cancer and its response to treatment, have postulated prognostic value. Here we searched for potential association between phenotypic ITH analyzed by mass spectrometry imaging (MSI) and prognosis of head and neck cancer. The study involved tissue specimens resected from 77 patients with locally advanced oral squamous cell carcinoma, including 37 patients where matched samples of primary tumor and synchronous lymph node metastases were analyzed. A 3-year follow-up was available for all patients which enabled their separation into two groups: with no evidence of disease (NED, n = 41) and with progressive disease (PD, n = 36). After on-tissue trypsin digestion, peptide maps of all cancer regions were segmented using an unsupervised approach to reveal their intrinsic heterogeneity. We found that intra-tumor similarity of spectra was higher in the PD group and diversity of clusters identified during image segmentation was higher in the NED group, which indicated a higher level of ITH in patients with more favorable outcomes. Signature of molecular components that correlated with long-term outcomes could be associated with proteins involved in the immune functions. Furthermore, a positive correlation between ITH and histopathological lymphocytic host response was observed. Hence, we proposed that a higher level of ITH revealed by MSI in cancers with a better prognosis could reflect the presence of heterotypic components of tumor microenvironment such as infiltrating immune cells enhancing the response to the treatment.

Single‐cell RNA sequencing in oral science: Current awareness and perspectives

The emergence of single‐cell RNA sequencing enables simultaneous sequencing of thousands of cells, making the analysis of cell population heterogeneity more efficient. In recent years, single‐cell RNA sequencing has been used in the investigation of heterogeneous cell populations, cellular developmental trajectories, stochastic gene transcriptional kinetics, and gene regulatory networks, providing strong support in life science research. However, the application of single‐cell RNA sequencing in the field of oral science has not been reviewed comprehensively yet. Therefore, this paper reviews the development and application of single‐cell RNA sequencing in oral science, including fields of tissue development, teeth and jaws diseases, maxillofacial tumors, infections, etc., providing reference and prospects for using single‐cell RNA sequencing in studying the oral diseases, tissue development, and regeneration.

Development of a novel signature derived from single cell RNA-sequencing for preoperative prediction of lymph node metastasis in head and neck squamous cell carcinoma

BACKGROUND

Lymph node metastasis (LNM) is considered as an adverse prognostic indicator for cancer patients. Preoperative knowledge of LNM is valuable for pretreatment decision making. Here, we sought to develop and validate an LNM signature for preoperative prediction of LNM in patients with head and neck squamous cell carcinoma (HNSCC).

METHODS

By studying single cell RNA-sequencing data (scRNA-seq), differentially expressed mRNA were selected and analyzed through univariate logistic regression and least absolute shrinkage and selection operator (LASSO) to identify an LNM signature. Multivariate logistic regression was utilized to establish an LNM nomogram incorporating LNM signature and T-classification.

RESULTS

The LNM signature was significantly associated with lymph node status and prognosis. The LNM signature and LNM nomogram displayed a robust predictive effect.

CONCLUSION

Our study reveals that LNM signature is a powerful biomarker for preoperative prediction of LNM in patients with HNSCC, which may be effective to realize individualized outcome prediction.

Dynamic changes of exhaustion features in T cells during oral carcinogenesis

OBJECTIVES

This study aimed to clarify the dynamic changes of exhaustion features in T cells during oral carcinogenesis.

MATERIALS AND METHODS

Mice were randomly divided into 4NQO group and control group. The exhaustion features of CD4+ and CD8+ T cells of both groups were detected by flow cytometry. Furthermore, multiplex immunohistochemistry was used to evaluate the expression of inhibitory receptors in human normal, dysplastic, and carcinogenesis tissues. Finally, anti-PD-1 antibody treatment was performed at the early premalignant phase of oral carcinogenesis.

RESULTS

The proportion of naive T cells in 4NQO group was lower than those in control group, while the proportion of effector memory T cells was higher in 4NQO group. The expression of inhibitory receptors on CD4+ and CD8+ T cells increased gradually during carcinogenesis. In contrast, the secretion of cytokines by CD4+ and CD8+ T cells decreased gradually with the progression stage. Strikingly, those changes occurred before the onset of oral carcinogenesis. The expression of inhibitory receptors on T cells increased gradually as the human tissues progressed from normal, dysplasia to carcinoma. Interestingly, PD-1 blockade at the early premalignant phase could reverse carcinogenesis progression by restoring T cell function.

CONCLUSIONS

T-cell dysfunction was established at the early premalignant phase of oral carcinogenesis; PD-1 blockade at the early premalignant phase can effectively reverse T-cell exhaustion features and then prevent carcinogenesis progression.

Single-Cell Profiling Reveals Heterogeneity of Primary and Lymph Node Metastatic Tumors and Immune Cell Populations and Discovers Important Prognostic Significance of CCDC43 in Oral Squamous Cell Carcinoma

Although substantial progress has been made in biological research and clinical treatment in recent years, the clinical prognosis of oral squamous cell carcinoma (OSCC) is still not satisfactory. Tumor immune microenvironment (TIME) is a potential target, which plays an essential role in the response of anti-tumor immunity and immunotherapy. In this study, we used scRNA-seq data, revealing the heterogeneity of TIME between metastatic and primary site. We found that in the metastatic site, the content of cytotoxic T cells and classical activated macrophages (M1 macrophages) increases significantly, while alternately activated macrophages (M2 macrophages) and inflammatory cancer-associated fibroblasts (iCAFs) decrease, which may be due to the increased immunogenicity of OSCC cells in the metastatic site and the changes in some signal pathways. We also found that iCAFs may recruit alternately activated macrophages (M2 macrophages) by secreting CXCL12. Then, we described a regulatory network for communication between various TIME cells centered on OSCC cells, which can help to clarify the possible mechanism of lymph node metastasis in OSCC cells. By performing pseudotime trajectory analysis, we found that the expression CCDC43 is upregulated in more advanced OSCC cells and is an independent prognostic factor for poor living conditions. Other than this, the high expression of CCDC43 may impair the antitumor immunity of the human body and promote the metastasis of OSCC cells. Our research provides a profound insight into the immunological study of OSCC and an essential resource for future drug discovery.

MEK1/2 inhibition transiently alters the tumor immune microenvironment to enhance immunotherapy efficacy against head and neck cancer

Background

Although the mitogen-activated protein kinases (MAPK) pathway is hyperactive in head and neck cancer (HNC), inhibition of MEK1/2 in HNC patients has not shown clinically meaningful activity. Therefore, we aimed to characterize the effect of MEK1/2 inhibition on the tumor microenvironment (TME) of MAPK-driven HNC, elucidate tumor-host interaction mechanisms facilitating immune escape on treatment, and apply rationale-based therapy combination immunotherapy and MEK1/2 inhibitor to induce tumor clearance.

Methods

Mouse syngeneic tumors and xenografts experiments were used to analyze tumor growth in vivo. Single-cell cytometry by time of flight, flow cytometry, and tissue stainings were used to profile the TME in response to trametinib (MEK1/2 inhibitor). Co-culture of myeloid-derived suppressor cells (MDSC) with CD8⁺ T cells was used to measure immune suppression. Overexpression of colony-stimulating factor-1 (CSF-1) in tumor cells was used to show the effect of tumor-derived CSF-1 on sensitivity to trametinib and anti-programmed death- 1 (αPD-1) in mice. In HNC patients, the ratio between CSF-1 and CD8A was measured to test the association with clinical benefit to αPD-1 and αPD-L1 treatment.

Results

Using preclinical HNC models, we demonstrated that treatment with trametinib delays HNC initiation and progression by reducing tumor cell proliferation and enhancing the antitumor immunity of CD8⁺ T cells. Activation of CD8⁺ T cells by supplementation with αPD-1 antibody eliminated tumors and induced an immune memory in the cured mice. Mechanistically, an early response to trametinib treatment sensitized tumors to αPD-1-supplementation by attenuating the expression of tumor-derived CSF-1, which reduced the abundance of two CSF-1R⁺CD11c⁺ MDSC populations in the TME. In contrast, prolonged treatment with trametinib abolished the antitumor activity of αPD-1, because tumor cells undergoing the epithelial to mesenchymal transition in response to trametinib restored CSF-1 expression and recreated an immune-suppressive TME.

Conclusion

Our findings provide the rationale for testing the trametinib/αPD-1 combination in HNC and highlight the importance of sensitizing tumors to αPD-1 by using MEK1/2 to interfere with the tumor–host interaction. Moreover, we describe the concept that treatment of cancer with a targeted therapy transiently induces an immune-active microenvironment, and supplementation of immunotherapy during this time further activates the antitumor machinery to cause tumor elimination.

Pyroptosis patterns and immune infiltrates characterization in head and neck squamous cell carcinoma

Background

Pyroptosis plays an essential role in tumor immune responses and inflammation related to chemotherapy. Herein, we studied the characteristic patterns of pyroptosis in head and neck squamous cell carcinoma (HNSCC) to determine their prognostic and therapeutic effects.

Methods

Consensus clustering analysis was performed to classify patients into pyroptosis or gene clusters. A novel pyroptosis score was constructed by principal component analysis. Kaplan–Meier survival curves were used to show the prognostic value. We also assessed the functional enrichment, tumor mutation burden, immune cell infiltration, and the sensitivity to chemotherapy and immunotherapy between high and low pyroptosis score group.

Results

Two distinct pyroptosis clusters were defined based on the mRNA expression profiles of PRGs, which were related to immune activation in HNSCC. Notably, a pyroptosis score was constructed according to different expression gene signatures, and then, each HNSCC patient was classified into a low or high pyroptosis score group. Patients with low pyroptosis scores had better immunotherapeutic responses and higher sensitivities to chemotherapeutic agents (paclitaxel, docetaxel, and gemcitabine). Kaplan–Meier survival curves showed that the pyroptosis patterns were independent prognostic indicators regardless of the level of tumor mutation burden.

Conclusions

Pyroptosis plays an essential role in immune infiltration in HNSCC. Quantifying the pyroptosis score of individual patients might suggest prognostic, immunotherapeutic, and chemotherapeutic strategies for HNSCC.

An Immune Feature-Based, Three-Gene Scoring System for Prognostic Prediction of Head-and-Neck Squamous Cell Carcinoma

Head-and-neck squamous cell carcinoma (HNSCC) is characterized by a high frequency of neck lymph node metastasis (LNM), a key prognostic factor. Therefore, identifying the biological processes during LNM of HNSCC has significant clinical implications for risk stratification. This study performed Gene Ontology enrichment analysis of differentially expressed genes between tumors with LNM and those without LNM and identified the involvement of immune response in the lymphatic metastasis of HNSCC. We further identified greater infiltrations of CD8⁺ T cells in tumors than in adjacent normal tissues through immunochemistry in the patient cohort (n = 62), indicating the involvement of CD8⁺ T cells in the antitumor immunity. Hierarchical clustering analysis was conducted to initially identify the candidate genes relevant to lymphocyte-mediated antitumor response. The candidate genes were applied to construct a LASSO Cox regression analysis model. Three genes were eventually screened out as progression-related differentially expressed candidates in HNSCC and a risk scoring system was established based on LASSO Cox regression model to predict the outcome in patients with HNSCC. The score was calculated using the formula: 0.0636 × CXCL11 - 0.4619 × CXCR3 + 0.2398 × CCR5. Patients with high scores had significantly worse overall survival than those with low scores (p < 0.001). The risk score showed good performance in characterizing tumor-infiltrating lymphocytes and provided a theoretical basis for stratifying patients receiving immune therapies. Additionally, a nomogram including the risk score, age, and TNM stage was constructed. The prediction model displayed marginally better discrimination ability and higher agreement in predicting the survival of patients with HNSCC compared with the TNM stage.

The RNA Methylation Modification 5-Methylcytosine Impacts Immunity Characteristics, Prognosis and Progression of Oral Squamous Cell Carcinoma by Bioinformatics Analysis

Disorders pertaining to 5-methylcytosine (m5C) modifications are involved in the pathological process of many diseases. However, the effect of m5C on the tumorigenesis and progression of oral squamous cell carcinoma (OSCC) remains unclear. In this study, we integrated the genomic and clinical data of 558 OSCC samples to comprehensively evaluate m5C modification patterns. Based on 16 m5C methylation regulators, two m5C modification clusters were identified with distinct tumor immune microenvironment (TIME) characteristics and prognosis in OSCC. We then performed weighted gene co-expression network analysis (WGCNA) to identify m5C modification cluster-related modules. Genes in the selected module were chosen to construct the m5Cscore scoring system for evaluating m5C modification pattern in individual OSCC patients. Patients with a high m5Cscore had higher immune, stromal, and ESTIMATE scores; lower tumor purity score; lower immune activity; and higher tumor mutational burden. The overall survival rate and progression-free survival rate were markedly worse and the tumor recurrence rate was higher in OSCC patients with a high m5Cscore. Furthermore, patients with oral leukoplakia who also had a high m5Cscore had a higher risk of deterioration to OSCC. This study demonstrated that m5C modification patterns might affect the TIME in OSCC. m5Cscore may provide a new approach for predicting the prognosis and progression of OSCC.

The Impact of m1A Methylation Modification Patterns on Tumor Immune Microenvironment and Prognosis in Oral Squamous Cell Carcinoma

N1-methyladenosine (m1A) modification widely participates in the occurrence and progression of numerous diseases. Nevertheless, the potential roles of m1A in the tumor immune microenvironment (TIME) are still not fully understood. Based on 10 m1A methylation regulators, we comprehensively explored the m1A modification patterns in 502 patients with oral squamous cell carcinoma (OSCC). The m1A modification patterns were correlated with TIME characteristics and the m1A score was established to evaluate the effect of the m1A modification patterns on individual OSCC patients. The TIME characteristics and survival outcomes under the three m1A modification patterns were significantly distinct. OSCC patients in the high m1A score group were characterized by poorer prognosis, lower immune infiltration, lower ssGSEA score, lower expression levels of immune checkpoint molecules, and higher tumor mutation loads. The present study revealed that m1A modification might be associated with the TIME in OSCC, and has potential predictive ability for the prognosis of OSCC.