Multi-omics characterization of a scoring system to quantify hypoxia patterns in patients with head and neck squamous cell carcinoma.

BACKGROUND: The 5-year survival rate of patients with head and neck squamous cell carcinoma (HNSCC) remains < 50%. Hypoxia patterns are a hallmark of HNSCC that are associated with its occurrence and progression. However, the precise role of hypoxia during HNSCC, such as the relationship between hypoxia, tumor immune landscape and cell communication orchestration remains largely unknown. The current study integrated data from bulk and single-cell RNA sequencing analyses to define the relationship between hypoxia and HNSCC. METHODS: A scoring system named the hypoxia score (HS) was constructed based on hypoxia-related genes (HRGs) expression. The predictive value of HS response for patient outcomes and different treatments was evaluated. Single-cell datasets and cell communication were utilized to rule out cell populations which hypoxia targeted on. RESULTS: The survival outcomes, immune/Estimate scores, responses to targeted inhibitors, and chemotherapeutic, and immunotherapy responses were distinct between a high HS group and a low HS group (all P < 0.05). Single-cell datasets showed different distributions of HS in immune cell populations (P < 0.05). Furthermore, HLA-DPA1/CD4 axis was identified as a unique interaction between CD4 + T Conv and pDC cells. CONCLUSIONS: Altogether, the quantification for hypoxia patterns is a potential biomarker for prognosis, individualized chemotherapeutic and immunotherapy strategies. The portrait of cell communication characteristics over the HNSCC ecosystem enhances the understanding of hypoxia patterns in HNSCC.

GAPDH: unveiling its impact as a key hypoxia-related player in head and neck squamous cell carcinoma tumor progression, prognosis, and therapeutic potential.

Head and neck squamous cell carcinoma (HNSCC), characterized by hypoxia patterns, ranks as the sixth most prevalent malignant tumor worldwide. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) plays a role in oncogenesis under hypoxic conditions in various cancers. However, its precise function in HNSCC, especially under varied hypoxic conditions, including at high altitudes, remains unclear. Elevated GAPDH mRNA and protein levels in HNSCC relative to normal tissues have been demonstrated through data from The Cancer Genome Atlas (TCGA), GSE29330, and the Human Protein Atlas (P<0.05). This elevation was further confirmed through in vitro experiments utilizing two HNSCC cell lines and a normal oral mucosal epithelial cell line. Additionally, data from TCGA and GSE41613 reveal a correlation between elevated GAPDH expression and diminished overall and progression-free survival in patients (P<0.05). Subsequent analysis identifies GAPDH as an independent risk factor for HNSCC (P<0.05). Using the ESTIMATE and single-sample gene set enrichment analysis (ssGSEA) algorithms, high GAPDH expression was found to be associated with reduced immune scores and diminished anti-tumor cell infiltration, such as CD8+ T cells, in TCGA and GSE41613 datasets (P<0.05). Analysis of single-cell RNA sequencing data from GSE139324 suggests that elevated GAPDH expression hinders communication between plasmacytoid dendritic cells and mast cells (P<0.05). Furthermore, in the TCGA and GSE41613 datasets, GAPDH's biological function is closely tied to hypoxia through Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Variation Analysis (GSVA) analyses. Moreover, its expression is linked to one cuproptosis-related gene, five N6-methyladenosine-related genes, six immune checkpoint genes, and pivotal pathways such as MYC and E2F (P<0.05). GAPDH showed excellent predictive value in estimating the efficacy of immunotherapy and 11 anti-tumor drugs (e.g., cisplatin) (P<0.05), using TIDE and pRRophetic algorithms on the TCGA and GSE41613 datasets. Under 1% O2 in vitro, HNSCC cells show elevated GAPDH expression, leading to decreased apoptosis and increased migration, clonogenicity, invasiveness, and resistance to cisplatin (P<0.05). At 5% O2, these effects persisted, albeit less pronouncedly. Inhibiting GAPDH reversed these effects under all oxygen concentrations (P<0.05). Overall, our findings reveal GAPDH as a key hypoxia-related player influencing tumor progression, prognosis, and therapeutic potential in HNSCC.