Integrative analysis of the prognostic value and immune microenvironment of mitophagy-related signature for multiple myeloma

Abstract

Abstract Background Multiple myeloma (MM) is a fatal malignant tumor in hematology. Mitophagy plays vital roles in the pathogenesis and drug sensitivity of MM. Methods We acquired transcriptomic expression data and clinical index of MM patients from NCI public database, and 36 genes involved in mitophagy from the gene set enrichment analysis (GSEA) database. Least absolute shrinkage and selection operator (LASSO) Cox regression analysis was conducted to construct a risk score prognostic model. Kaplan–Meier survival analysis and receiver operation characteristic curves (ROC) were conducted to identify the efficiency of prognosis and diagnosis. ESTIMATE algorithm and immune-related single-sample gene set enrichment analysis (ssGSEA) was performed to uncover the level of immune infiltration. QRT-PCR was performed to verify gene expression in clinical samples of MM patients. The sensitivity to chemotherapy drugs was evaluated upon the database of the genomics of drug sensitivity in cancer (GDSC). Results Fifty mitophagy-related genes were differently expressed in two independent cohorts. Ten out of these genes were identified to be related to MM overall survival (OS) rate. A prognostic risk signature model was built upon on these genes: VDAC1, PINK1, VPS13C, ATG13, and HUWE1, which predicted the survival of MM accurately and stably both in training and validation cohorts. MM patients suffered more adverse prognosis showed more higher risk core. In addition, the risk score was considered as an independent prognostic element for OS of MM patients by multivariate cox regression analysis. Functional pathway enrichment analysis of differentially expressed genes (DEGs) based on risk score showed terms of cell cycle, immune response, mTOR pathway, and MYC targets were obviously enriched. Furthermore, MM patients with higher risk score were observed lower immune scores and lower immune infiltration levels. The results of qRT-PCR verified VDAC1, PINK1, and HUWE1 were dysregulated in new diagnosed MM patients. Finally, further analysis indicated MM patients showed more susceptive to bortezomib, lenalidomide and rapamycin in high-risk group. Conclusion Our research provided a neoteric prognostic model of MM based on mitophagy genes. The immune infiltration level based on risk score paved a better understanding of the participation of mitophagy in MM.