Identification of an Autophagy-Related Signature for Prognosis and Immunotherapy Response Prediction in Ovarian Cancer

Abstract

Background: Ovarian cancer (OC) is one of the most malignant tumors in the female reproductive system, with a poor prognosis. Various responses to treatments including chemotherapy and immunotherapy are observed among patients due to their individual characteristics. Applicable prognostic markers could make it easier to refine risk stratification for OC patients. Autophagy is closely implicated in the occurrence and development of tumors, including OC. Whether autophagy -related genes can be used as prognostic markers for OC patients remains unclear. Methods: The gene transcriptome data of 374 OC patients were downloaded from The Cancer Genome Atlas (TCGA) database. The correlation between the autophagy levels and outcomes of OC patients was identified through the single sample gene set enrichment analysis (ssGSEA). Recognized molecular markers of autophagy in different clinical specimens were detected by immunohistochemistry (IHC) assay. The gene set enrichment analysis (GSEA), ESTIMATE, and CIBERSORT analysis were applied to explore the correlation of autophagy with the tumor immune microenvironment (TIME). Single-cell RNA-sequencing (scRNA-seq) data from seven OC patients were included for characterizing cell-cell interaction patterns of autophagy-high or low tumor cells. Machine learning, Stepwise Cox regression and LASSO-Cox analysis were used to screen autophagy hub genes, which were used to establish an autophagy-related signature for prognosis evaluation. Four tumor immunotherapy cohorts were obtained from the GEO (Gene Expression Omnibus) database and the literature for autophagy risk score validation. Results: The autophagy levels were closely related to the prognosis of the OC patients. Additionally, the autophagy levels were correlated with TIME status including immune score, and immune-cell infiltration. The scRNA-seq analysis found that tumor cells with high or low autophagy levels had different interactions with immune cells, especially macrophages. Eight autophagy-hub genes (ZFYVE1, AMBRA1, LAMP2, TRAF6, PDPK1, ATG2B, DAPK1 and TP53INP2) were screened for an autophagy-related signature. According to this signature, higher risk score was correlated with poor prognosis and better immunotherapy response in the OC patients. Conclusions: The autophagy-related signature is applicable to predict the prognosis and immune checkpoint inhibitors (ICIs) therapy efficiency in OC patients. It is possible to identify OC patients who will respond to ICIs therapy and have a favorable prognosis, although more verification is needed.