Background
Immunoregulatory drugs regulate the ubiquitin-proteasome system, which is the main treatment for multiple myeloma (MM) at present. In this study, bioinformatics analysis was used to construct the risk model and evaluate the prognostic value of ubiquitination-related genes in MM.
Methods and results
The data on ubiquitination-related genes and MM samples were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The consistent cluster analysis and ESTIMATE algorithm were used to create distinct clusters. The MM prognostic risk model was constructed through single-factor and multiple-factor analysis. The ROC curve was plotted to compare the survival difference between high- and low-risk groups. The nomogram was used to validate the predictive capability of the risk model. A total of 87 ubiquitination-related genes were obtained, with 47 genes showing high expression in the MM group. According to the consistent cluster analysis, 4 clusters were determined. The immune infiltration, survival, and prognosis differed significantly among the 4 clusters. The tumor purity was higher in clusters 1 and 3 than in clusters 2 and 4, while the immune score and stromal score were lower in clusters 1 and 3. The proportion of B cells memory, plasma cells, and T cells CD4 naïve was the lowest in cluster 4. The model genes KLHL24, HERC6, USP3, TNIP1, and CISH were highly expressed in the high-risk group. AICAr and BMS.754,807 exhibited higher drug sensitivity in the low-risk group, whereas Bleomycin showed higher drug sensitivity in the high-risk group. The nomogram of the risk model demonstrated good efficacy in predicting the survival of MM patients using TCGA and GEO datasets.
Conclusions
The risk model constructed by ubiquitination-related genes can be effectively used to predict the prognosis of MM patients. KLHL24, HERC6, USP3, TNIP1, and CISH genes in MM warrant further investigation as therapeutic targets and to combat drug resistance.