Integrative single cell transcriptomic analysis reveals 3p deletion associated tumor microenvironment and chemoresistance in head and neck squamous cell carcinoma.

Head and neck squamous cell carcinoma (HNSCC) remains a prevalent and lethal malignancy, with a five-year survival rate of just 50% for cases of locally advanced disease. Chromosomal aberrations, particularly the deletion of the short arm of chromosome 3 (3p), have been strongly associated with poor prognosis and more aggressive tumor phenotypes. The tumor microenvironment (TME) plays a pivotal role in tumor progression and resistance to therapy. This study aims to elucidate the impact of 3p deletion on the TME, immune cell infiltration, and treatment response in HNSCC, to identify novel therapeutic targets to improve patient outcomes. We analyzed single-cell RNA sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) and bulk transcriptome data from The Cancer Genome Atlas (TCGA). Pseudo-time trajectory and cell-cell communication analyses were performed with the Monocle and CellChat packages. The Wilcoxon test was used to evaluate the differential gene expression between wild-type (wt) and mutant (mut) groups. Prognostic models were developed using the Least Absolute Shrinkage and Selection Operator (LASSO) algorithm and Cox regression analyses to find the genes related to survival, with survival analysis conducted via Kaplan-Meier curves. Gene set enrichment analysis (GSEA) was employed to investigate pathway dysregulation, and immune cell infiltration was assessed using various immune scoring methodologies to explore the differences immune environment. The Tumor Immune Dysfunction and Exclusion (TIDE) database was utilized to predict the potential efficacy of immune checkpoint inhibitors. mRNA and protein expression levels of SPP1 were examined by RT-qPCR and Western blotting, while cell proliferation was assessed using the CCK8 assay. The mut group demonstrated significant alterations in cellular composition, characterized by increased endothelial cells and macrophages and decreased fibroblasts and CD8 + T cells, indicative of an immunosuppressive TME. Differential expression analysis revealed downregulation of immune pathways, including antigen processing and presentation, T cell receptor signaling, and B cell receptor signaling pathways in the mut group, along with enhanced metabolic activity in glycolysis and lipid metabolism. The prognostic model identified nine key genes associated with poor survival in HNSCC. The mut group exhibited poorer overall survival and a more immunosuppressive microenvironment compared to the wt group, which correlated with the outcomes observed in high-risk versus low-risk groups. High-risk patients also showed a diminished response to immunotherapy compared to low-risk patients. Additionally, SPP1 emerged as a critical gene associated with chemotherapy resistance and macrophage M2 polarization. This study demonstrates that 3p deletion significantly reshapes the TME, contributing to poor prognosis in HNSCC by fostering an immunosuppressive environment and enhancing chemoresistance. These findings highlight the potential for developing targeted therapies that address the genetic and immunological landscape of HNSCC.

Chen X ，Xu S ，Pan J ，Xu W ，Yang H ，Chen X ，Chen R ，Wang Y ，Qiu S ... -  《Scientific Reports》

Characterizing tumor biology and immune microenvironment in high-grade serous ovarian cancer via single-cell RNA sequencing: insights for targeted and personalized immunotherapy strategies.

High-grade serous ovarian cancer (HGSOC), the predominant subtype of epithelial ovarian cancer, is frequently diagnosed at an advanced stage due to its nonspecific early symptoms. Despite standard treatments, including cytoreductive surgery and platinum-based chemotherapy, significant improvements in survival have been limited. Understanding the molecular mechanisms, immune landscape, and drug sensitivity of HGSOC is crucial for developing more effective and personalized therapies. This study integrates insights from cancer immunology, molecular profiling, and drug sensitivity analysis to identify novel therapeutic targets and improve treatment outcomes. Utilizing single-cell RNA sequencing (scRNA-seq), the study systematically examines tumor heterogeneity and immune microenvironment, focusing on biomarkers influencing drug response and immune activity, aiming to enhance patient outcomes and quality of life. scRNA-seq data was obtained from the GEO database in this study. Differential gene expression was analyzed using gene ontology and gene set enrichment methods. InferCNV identified malignant epithelial cells, while Monocle, Cytotrace, and Slingshot software inferred subtype differentiation trajectories. The CellChat software package predicted cellular communication between malignant cell subtypes and other cells, while pySCENIC analysis was utilized to identify transcription factor regulatory networks within malignant cell subtypes. Finally, the analysis results were validated through functional experiments, and a prognostic model was developed to assess prognosis, immune infiltration, and drug sensitivity across various risk groups. This study investigated the cellular heterogeneity of HGSOC using scRNA-seq, focusing on tumor cell subtypes and their interactions within the tumor microenvironment. We confirmed the key role of the C2 IGF2+ tumor cell subtype in HGSOC, which was significantly associated with poor prognosis and high levels of chromosomal copy number variations. This subtype was located at the terminal differentiation of the tumor, displaying a higher degree of malignancy and close association with stage IIIC tissue types. The C2 subtype was also associated with various metabolic pathways, such as glycolysis and riboflavin metabolism, as well as programmed cell death processes. The study highlighted the complex interactions between the C2 subtype and fibroblasts through the MK signaling pathway, which may be closely related to tumor-associated fibroblasts and tumor progression. Elevated expression of PRRX1 was significantly connected to the C2 subtype and may impact disease progression by modulating gene transcription. A prognostic model based on the C2 subtype demonstrated its association with adverse prognosis outcomes, emphasizing the importance of immune infiltration and drug sensitivity analysis in clinical intervention strategies. This study integrates molecular oncology, immunotherapy, and drug sensitivity analysis to reveal the mechanisms driving HGSOC progression and treatment resistance. The C2 IGF2+ tumor subtype, linked to poor prognosis, offers a promising target for future therapies. Emphasizing immune infiltration and drug sensitivity, the research highlights personalized strategies to improve survival and quality of life for HGSOC patients.

Zhao F ，Jiang X ，Li Y ，Huang T ，Xiahou Z ，Nie W ，Li Q ... -  《-》

SLC2A3 promotes head and neck squamous cancer developing through negatively regulating CD8(+) T cell in tumor microenvironment.

Recent studies have identified SLC2A3 as being abnormally upregulated in multiple tumor types, correlating with poor survival and disrupted microenvironments. However, its prognostic significance in head and neck squamous cell carcinoma (HNSC) remains underexplored. In this study, SLC2A3 was screened as a potential risk gene influencing both immune and tumor components within the tumor microenvironment (TME) of 504 HNSC patients from the TCGA database. Immune infiltration analyses and clinical significance on SLC2A3 were conducted using ESTIMATE, CIBERSORT, ssGSEA, TIMER and clinical prognosis parameters. Additionally, the single-cell dataset is used to analyze the expression of SLC2A3 in various subpopulations. The magnetic activated cell sorting (MACS) is used to isolate CD8+ T cells from PBMCs or tumor tissues. Flow cytometry is used to identify purified and activated CD8+ T cells. GSEA and WB were used to investigate the molecular mechanism of SLC2A3 in CD8+ T cells. The co-culture system of CD8+ T cells and TU686 was used to investigate the effects of SLC2A3 on immune cells and tumor development. In this study, SLC2A3 was identified as a potential risk gene affecting both immune cells and tumor components within the TME of 504 HNSC patients derived from the TCGA database. We conducted immune infiltration analyses and assessed the clinical significance of SLC2A3 using various bioinformatics tools, including ESTIMATE, CIBERSORT, ssGSEA, and TIMER, along with clinical prognosis parameters. The single-cell RNA sequencing dataset was utilized to examine SLC2A3 expression across different cellular subpopulations. Magnetic activated cell sorting (MACS) was employed to isolate CD8+ T cells from peripheral blood mononuclear cells (PBMCs) or tumor tissues. Flow cytometry was implemented to confirm the purity and activation state of the isolated CD8+ T cells. GSEA and Western blot were applied to explore the molecular mechanisms underlying SLC2A3's role in CD8+ T cells. Lastly, a co-culture system involving CD8+ T cells and TU686 tumor cells was established to study the impact of SLC2A3 on immune cell function and tumor progression. SLC2A3 emerges as an actively variable gene within the immune and stromal components of the TME, linked to aggravated immune infiltration and poor clinical outcomes. The upregulated expression of SLC2A3 is predominantly enriched in immune-related biological processes and linked to the suppression of CD8+ T cells, which are crucial for the survival of HNSC patients. Furthermore, SLC2A3 exhibits specific overexpression in CD8+ T cells and may potentially trigger ferroptosis. Knockdown of SLC2A3 led to a significant increase in the proliferation of CD8+ T cells compared to those without knockdown. In co-culture systems, CD8+ T cells with SLC2A3 knockdown demonstrated an enhanced ability to eliminate tumor cells compared to those without the knockdown. SLC2A3 is associated with changes in the TME and prognostic indicators. Moreover, high SLC2A3 expression in CD8+ T cells may drive cell death through ferroptosis, fostering tumor progression.

Jiang W ，Xu S ，Zhao M ，Li C ... -  《Scientific Reports》

CD8 + T-Cell-Related Genes: Deciphering Their Role in the Pancreatic Adenocarcinoma TME and Their Effect on Prognosis.

Because of the unique tumor microenvironment (TME), immunotherapy and targeted therapies have shown limited efficacy in treating pancreatic adenocarcinoma (PAAD). CD8 + T cells play crucial roles in regulating the TME in PAAD; therefore, exploring the function of CD8 + T-cell-related genes (CD8RGs) in PAAD has high potential clinical value and could provide a comprehensive understanding of the microenvironment of PAAD. We employed the weighted gene coexpression network analysis and CIBERSORT algorithms to assess PAAD transcriptome data from The Cancer Genome Atlas (TCGA) dataset and identify modules strongly associated with CD8 + T cell infiltration. Using least absolute shrinkage and selection operator regression analysis and Kaplan-Meier curves, we developed a prognostic risk score model for patients with PAAD. We validated this model using single-cell and transcriptome datasets obtained from the Gene Expression Omnibus (GEO). We also examined the correlations between the risk score and factors such as the TME, clinical characteristics, and tumor mutation burden (TMB). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis were performed on differentially expressed genes between the high- and low-risk groups. In addition, the Tumor Immune Dysfunction and Exclusion website and "pRRophetic" R package were used to predict response to immunotherapy and chemotherapy in the high- and low-risk groups, respectively. Finally, we analyzed the expressions of hub genes at the cellular level with quantitative real-time PCR. A risk model based on five CD8RGs was established and validated using TCGA and GEO datasets. The low-risk group exhibited significantly longer overall and progression-free survival. A positive correlation between the TMB and the risk score was observed. The TME analysis revealed a significant correlation between the risk score and immune function, as well as immune checkpoints. The expression of hub genes was significantly correlated with the infiltration level of CD8 + T cells. The high-risk group responded better to immunotherapy, paclitaxel, cisplatin, mitomycin C, afatinib (BIBW2992), and gefitinib. In contrast, the low-risk group showed higher sensitivity to sunitinib, MK.2206, palbociclib (PD.0332991), and axitinib. Compared with that in normal pancreatic epithelial cells, the expression levels of BCL11A, PHOSPHO1, and GNG7 were significantly decreased, while those of KLK11 and VCAM1 were significantly increased in pancreatic tumor cells. CD8RGs play an important role in regulating the TME of PAAD. Five hub genes-BCL11A, KLK11, GNG7, PHOSPHO1, and VCAM1-are closely associated with the prognosis of PAAD patients, providing new references for the exploration of biomarkers. Furthermore, our findings offer novel insights for clinical decision-making.

Zhang Y ，Hou H ，Zhang X ，Lan H ，Huo X ，Duan X ，Li Y ，Zhang X ，Zhou N ... -  《-》

[A novel glycolysis-related prognostic risk model for colorectal cancer patients based on single-cell and bulk transcriptomic data].

Yao K ，Xia J ，Zhang S ，Sun Y ，Ma J ，Zhu B ，Ren L ，Zhang C ... -  《-》