Head and Neck Squamous Cell Carcinoma Subtypes Based on Immunologic and Hallmark Gene Sets in Tumor and Non-tumor Tissues

Characterization of macrophages in head and neck squamous cell carcinoma and development of MRG-based risk signature

Macrophages are immune cells in the TME that can not only inhibit angiogenesis, extracellular matrix remodeling, cancer cell proliferation, and metastasis but also mediate the phagocytosis and killing of cancer cells after activation, making them key targets in anti-tumor immunotherapy. However, there is little research on macrophages and their relation to disease prognosis in HNSCC. Initially, we collected scRNA-seq, bulk RNA-seq, and clinical data. Subsequently, we identified macrophages and distinguished MRGs. Using the K-means algorithm, we performed consensus unsupervised clustering. Next, we used ssGSEA analysis to assess immune cell infiltration in MRG clusters. A risk model was established using multivariate Cox analysis. Then, Kaplan-Meier, ROC curves, univariate and multivariate COX analyses, and C-index was used to validate the predictive power of the signature. The TIDE method was applied to assess the response to immunotherapy in patients diagnosed with HNSCC. In addition, drug susceptibility predictions were made for the GDSC database using the calcPhenotype function. We found that 8 MRGs had prognostic potential. Patients in the MRG group A had a higher probability of survival, and MRG clusters A and B had different characteristics. Cluster A had a higher degree of expression and infiltration in MRG, indicating a closer relationship with MRG. The accuracy of the signature was validated using univariate and multivariate Cox analysis, C-index, and nomogram. Immune landscape analysis found that various immune functions were highly expressed in the low-risk group, indicating an improved response to immunotherapy. Finally, drugs with high sensitivity to HNSCC (such as 5-Fluorouracil, Temozolomide, Carmustine, and EPZ5676) were explored and analyze the malignant characteristics of HNSCC. We constructed a prognostic model using multivariate Cox analysis, consisting of 8 MRGs (TGM2, STC1, SH2D3C, PIK3R3, MAP3K8, ITGA5, ARHGAP4, and AQP1). Patients in the low-risk group may have a higher response to immunotherapy. The more prominent drugs for drug selection are 5-fluorouracil, temozolomide and so on. Malignant features associated with HNSCC include angiogenesis, EMT, and the cell cycle. This study has opened up new prospects for the prognosis, prediction, and clinical treatment strategy of HNSCC.

Multiple programmed cell death patterns and immune landscapes in bladder cancer: Evidence based on machine learning and multi-cohorts

BACKGROUND

Bladder cancer (BLCA) is the most prevalent malignant neoplasm of the urinary tract, and ranks seventh as the most frequent systemic neoplasm in males. Dysregulation of programmed cell death (PCD) has been implicated in various stages of cancer progression, including tumorigenesis, invasion, and metastasis. However, the correlation between multiple PCD modes and BLCA is lacking. Thus, a risk prediction model was built based on 12 models of PCD to predict prognosis and immunotherapy response in patients with BLCA.

METHODS

The RNA sequencing transcriptome data of BLCA were collected from the Cancer Genome Atlas Program (TCGA) and GEO datasets. Univariate Cox and LASSO regression analyzes were performed to identify PCD-related genes (PCDRGs) significant for prognosis. Multivariate Cox regression analysis was used to develop a prognostic model for PCD. Survival analysis and chi-squared test were employed to analyze the survival variations between different risk groups. Univariate and multivariate Cox analyses were performed to evaluate the model as an independent prognostic predictor. A nomogram was formulated using both clinical data and the model to predict the survival rates of BLCA patients. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes were performed to analyze and elucidate the molecular mechanisms and pathways operating within different risk score groups. Furthermore, the immune landscape was investigated and the efficacy of various anti-tumor drugs was evaluated for BLCA. Finally, consensus clustering analysis was adopted to explore the association between different PCD clusters and clinical characteristics.

RESULTS

Assessment of the public datasets and multivariate Cox analysis yielded 1254 PCDRGs, of which 10 PCDRGs for BLCA were identified. Based on the PCDRGs, a prognostic model was built for BLCA patient prognosis. Compared with the low-risk group, the high-risk group had a poorer prognosis. The model predicted area under the curve (AUC) values of 0.751, 0.753, and 0.763, respectively, for 1-, 3-, and 5-year survival of BLCA patients. The nomogram further demonstrated the credibility of the prognosis model. The low-risk group patients exhibited lower TIDE scores and higher TMB scores, implying better response of the low-risk group to immunotherapy. The consensus clustering analysis indicated that compared with PCD cluster A, PCD cluster B was significantly more expressed in PCDRGs, suggesting a closer relation of PCD cluster B to PCDRGs. Patients in PCD cluster B had lower risk scores.

CONCLUSION

To summarize, the effects of 12 PCD patterns on BLCA were synthesized and the correlation between PCD and BLCA was explored. These findings provide new and convincing evidence for individualized treatment of BLCA, and help guide the treatment strategy and improve the prognosis of BLCA patients.

Identification of molecular classification and gene signature for predicting prognosis and immunotherapy response in HNSCC using cell differentiation trajectories

Head and neck squamous cell carcinoma (HNSCC) is a highly heterogeneous malignancy with poor prognosis. This article aims to explore the clinical significance of cell differentiation trajectory in HNSCC, identify different molecular subtypes by consensus clustering analysis, and develop a prognostic risk model on the basis of differentiation-related genes (DRGs) for predicting the prognosis of HNSCC patients. Firstly, cell trajectory analysis was performed on single-cell RNA sequencing (scRNA-seq) data, four molecular subtypes were identified from bulk RNA-seq data, and the molecular subtypes were predictive of patient survival, clinical features, immune infiltration status, and expression of immune checkpoint genes (ICGs)s. Secondly, we developed a 10-DRG signature for predicting the prognosis of HNSCC patients by using weighted correlation network analysis (WGCNA), differential expression analysis, univariate Cox regression analysis, and multivariate Cox regression analysis. Then, a nomogram integrating the risk assessment model and clinical features can successfully predict prognosis with favorable predictive performance and superior accuracy. We projected the response to immunotherapy and the sensitivity of commonly used antitumor drugs between the different groups. Finally, we used the quantitative Reverse Transcription-Polymerase Chain Reaction (qRT-PCR) analysis and western blot to verify the signature. In conclusion, we identified distinct molecular subtypes by cell differentiation trajectory and constructed a novel signature based on differentially expressed prognostic DRGs, which could predict the prognosis and response to immunotherapy for patients and may provide valuable clinical applications in the treatment of HNSCC.

Gastric Cancer Subtypes in Tumour and Nontumour Tissues by Immunologic and Hallmark Gene Sets

A previous research study on differentiating gastric cancer (GC) into distinct subtypes or prognostic models was mostly based on GC tissues, which neglected the role of nontumour tissues in GC subtypes. The purpose of the research was to identify GC subtypes on the basis of tumour and adjacent nontumour tissues to assess the prognosis of GC patients. We characterized three GC subtypes on the basis of the immunologic and hallmark gene sets in GC and adjacent nontumour tissues: among them, the GC patients with subtype I had the longest survival time compared to patients with other subtypes. The classification was closely associated with T stage and pathological stage of GC patients. A prognostic model containing two gene sets was constructed by LASSO analysis. Kaplan–Meier analysis showed that patients in the high-risk group survived longer than those in the low-risk group and the two prognostic genes sets in the model were strongly correlated with survival status. Then, GO and KEGG analyses and PPI network show that nontumour and tumour tissues are influencing the prognosis of GC patients in separate manners. In summary, we emphasized the prognostic value of nontumour tissue in GC patients and proposed a novel insight that both changes in tumour and nontumour tissues should be taken into account when selecting a treatment strategy for GC.