Integrative analysis of cuproptosis‐associated genes for predicting immunotherapy response in single‐cell and multi‐cohort studies

Abstract Background The role of genes associated with the cuproptosis cell signaling pathway in prognosis and immunotherapy in ovarian cancer (OC) has been extensively investigated. In this study, we aimed to explore these mechanisms and establish a prognostic model for patients with OC using bioinformatics techniques. Methods We obtained the single cell sequencing data of ovarian cancer from the Gene Expression Omnibus (GEO) database and preprocessed the data. We analyzed a variety of factors including cuproptosis cell signal score, transcription factors, tumorigenesis and progression signals, gene set variation analysis (GSVA) and intercellular communication. Differential gene analysis was performed between groups with high and low cuproptosis cell signal scores, as well as Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses. Using bulk RNA sequencing data from The Cancer Genome Atlas, we used the least absolute shrinkage and selection operator (LASSO)‐Cox algorithm to develop cuproptosis cell signaling pathword‐related gene signatures and validated them with GEO ovarian cancer datasets. In addition, we analyzed the inherent rules of the genes involved in building the model using a variety of bioinformatics methods, including immune‐related analyses and single nucleotide polymorphisms. Molecular docking is used to screen potential therapeutic drugs. To confirm the analysis results, we performed various wet experiments such as western blot, cell counting kit 8 (CCK8) and clonogenesis tests to verify the role of the Von Willebrand Factor (VWF) gene in two ovarian cancer cell lines. Results Based on single‐cell data analysis, we found that endothelial cells and fibroblasts showed active substance synthesis and signaling pathway activation in OC, which further promoted immune cell suppression, cancer cell proliferation and metastasis. Ovarian cancer has a high tendency to metastasize, and cancer cells cooperate with other cells to promote disease progression. We developed a signature consisting of eight cuproptosis‐related genes (CRGs) (MAGEF1, DNPH1, RARRES1, NBL1, IFI27, VWF, OLFML3 and IGFBP4) that predicted overall survival in patients with ovarian cancer. The validity of this model is verified in an external GEO validation set. We observed active infiltrating states of immune cells in both the high‐ and low‐risk groups, although the specific cells, genes and pathways of activation differed. Gene mutation analysis revealed that TP53 is the most frequently mutated gene in ovarian cancer. We also predict small molecule drugs associated with CRGs and identify several potential candidates. VWF was identified as an oncogene in ovarian cancer, and the protein was expressed at significantly higher levels in tumor samples than in normal samples. The high‐score model of the cuproptosis cell signaling pathway was associated with the sensitivity of OC patients to immunotherapy. Conclusions Our study provides greater insight into the mechanisms of action of genes associated with the cuproptosis cell signaling pathway in ovarian cancer, highlighting potential targets for future therapeutic interventions.