Integrative single-cell and bulk transcriptomics reveal cuproptosis subtypes, prognostic models, and COL5A1 as a key gene in glioblastoma

  • Cell Signal. 2026 Oct:146:112691. doi: 10.1016/j.cellsig.2026.112691.
Ying Zhan  1 Yige Zhang  2 Jixing Miao  3 Xiaochuan Qu  3 Yuheng Li  2 Xiangyu Li  4 Yuhan Dong  1 Baoqin Liu  2 Yichen Song  5
Affiliations
  • 1. Department of Ultrasound, Shengjing Hospital of China Medical University, Shenyang 110004, Liaoning, China.
  • 2. Department of Biochemistry and Molecular Biology, School of Life Sciences, China Medical University, Shenyang 110122, Liaoning, China.
  • 3. Department of Biochemistry and Molecular Biology, School of Life Sciences, China Medical University, Shenyang 110122, Liaoning, China; The First Clinical Medical School, China Medical University, Shenyang 110122, Liaoning, China.
  • 4. Department of Neurosurgery, Shengjing Hospital of China Medical University, Shenyang 110004, Liaoning, China.
  • 5. Department of Neurosurgery, Shengjing Hospital of China Medical University, Shenyang 110004, Liaoning, China. Electronic address: [email protected].
Abstract

Glioblastoma (GBM) remains one of the most aggressive primary brain tumors with limited therapeutic options. Cuproptosis, a recently identified copper-dependent form of regulated cell death, has emerged as a potential therapeutic vulnerability in Cancer. Here, we systematically investigated the molecular and clinical significance of cuproptosis-related genes (CRGs) in GBM using integrated bulk and single-cell transcriptomic analyses. Based on the expression profiles of 31 differentially expressed CRGs, GBM patients were classified into two molecular subtypes with distinct prognostic and immune characteristics. Cluster 2 exhibited poorer survival and increased immune and stromal infiltration. To further assess the prognostic value of CRGs, we constructed a risk model using 110 subtype-associated genes and systematically compared 101 machine-learning algorithms. The SuperPC model achieved the best predictive performance and was successfully validated in two independent CGGA cohorts. The resulting risk score was closely associated with tumor microenvironment infiltration. Integrated single-cell analysis and experimental validation identified COL5A1 as a key GBM-associated gene. COL5A1 was specifically enriched in tumor cells across independent scRNA-seq datasets and positively correlated with the cuproptosis-related genes LOXL2 and SPARC. Functional assays demonstrated that COL5A1 knockdown suppressed GBM cell migration and invasion, reduced LOXL2 and SPARC expression, and enhanced sensitivity to Cuproptosis induction. Collectively, our study establishes a robust cuproptosis-associated prognostic framework for GBM and identifies COL5A1 as a potential regulator linking extracellular matrix remodeling, malignant progression, and cuproptosis-related vulnerability, providing new insights into GBM biology and potential therapeutic targets.

Keywords
COL5A1; Cuproptosis; Glioblastoma; Prognostic risk model; Single-cell RNA sequencing.