Cuproptosis/OXPHOS tendency prediction of prognosis and immune microenvironment of esophageal squamous cell carcinoma: Bioinformatics analysis and experimental validation

Abstract

Background: Cuproptosis is a newly discovered cell death mechanism that relies on mitochondrial respiration, for which oxidative phosphorylation (OXPHOS) is an essential part. However, the detailed mechanisms of cuproptosis associated with OXPHOS in esophageal squamous cell carcinoma (ESCC) and how this correlation affects prognosis still remains unclear.

Methods: scRNA-seq data of ESCC were downloaded from SRA (Sequence Read Archive) database. "AUCell" algorithm was used to grouping epithelial cells according to cuproptosis and OXPHOS score. Cell-cell communication, Pseudo-time Trajectory and transcription factor enrichment analysis were repectively conducted by "CellChat", "monocle3" package and "pySCENIC" algorithm. Univariate and LASSO cox regression analysis were used to construct the prognostic cuproptosis-OXPHOS signature. Finally, CCK-8 assay and DCFH-DA staining assay were respectively validated the sensitive and ROS production of elesclomol.

Results: scRNA-seq data were analyzed to identify 10 core cell types. According to the median scores for cuproptosis and OXPHOS, malignant epithelial cells were divided into double high, double low, and mixed groups. The double high group distributed at the end of the pseudo-time trajectory and harbored HMGA1(+) as specific transcriptional regulons. Knockdown of HMGA1 partly reversed the inhibition of cell viability visualized by CCK-8 assay, while reactive oxygen species (ROS) production by elesclomol was enhanced after HMGA1 silencing. Furthermore, the immunosuppressive signal was significantly increased in the double high group detected by 'CellChat' in single-cell data and 'ssGSEA' in bulk data followed by 'CIBERSORTx' algorithm. Finally, a new cuproptosis-OXPHOS prognostic signature (CNN2, ATP6V0E1, PSMD6, CCDC25, IGFBP2, MT1E, and RPS4Y1) was constructed for the prediction of the prognosis, and a high-risk group corresponding to a more sensitive tendency to erlotinib, dasatinib, and bosutinib treatment was identified.

Conclusions: Our study revealed the relationship between OXPHOS and tendency of cuproptosis in ESCC, and malignant cells with this characteristic exerted immunosuppressive signals and indicated poor prognosis. Furthermore, we constructed the regulatory network in high cuproptosis-OXPHOS ESCC and identified HMGA1 as a potential regulator molecule of cuproptosis mediated by elesclomol.

Keywords: Cuproptosis; ESCC; Elesclomol; Immune microenvironment; OXPHOS; Prognosis; scRNA-seq.