Identification of an Autophagy-Related Signature Based on Whole Bone Marrow Sequencing for the Prognosis and Immune Microenvironment Characterization of Multiple Myeloma

Myeloma (MM) is a malignant plasma cell disorder, which is incurable owing to its drug resistance. Autophagy performs an integral function in homeostasis, survival, and drug resistance in multiple myeloma (MM). Therefore, the purpose of the present research was to identify potential autophagy-related genes (ARGs) in patients with MM. We downloaded the transcriptomic data (GSE136400) of patients with MM, as well as the corresponding clinical data from the Gene Expression Omnibus (GEO); the patients were classified at random into two groups in a ratio of 6: 4, with 212 samples in the training dataset and 142 samples in the test dataset. Both multivariate and univariate Cox regression analyses were performed to identify autophagy-related genes. The univariate Cox regression analysis demonstrated that 26 ARGs had a significant correlation with overall survival (OS). We constructed an autophagy-related risk prognostic model based on six ARGs: EIF2AK2 (ENSG00000055332), KIF5B (ENSG00000170759), MYC (ENSG00000136997), NRG2 (ENSG00000158458), PINK1 (ENSG00000158828), and VEGFA (ENSG00000112715) using LASSO-Cox regression analysis to predict risk outcomes, which revealed substantially shortened OS duration in the high-risk cohort in contrast with that in the low-risk cohort. Therefore, the ARG-based model significantly predicted the MM patients' prognoses and was verified in an internal test set. Differentially expressed genes were found to be predominantly enriched in pathways associated with inflammation and immune regulation. Immune infiltration of tumor cells resulted in the formation of a strong immunosuppressive microenvironment in high-risk patients. The potential therapeutic targets of ARGs were subsequently analyzed via protein–drug network analysis. Therefore, a prognostic model for MM was established via a comprehensive analysis of ARGs, through using the clinical models; we have further revealed the molecular landscape features of multiple myeloma.