Abstract
Multiple myeloma (MM) is the second most frequent hematological cancer and is characterized by the clonal proliferation of malignant plasma cells. Genome-wide expression profiling (GEP) analysis with DNA microarrays has emerged as a powerful tool for biomedical research, generating a huge amount of data. Microarray analyses have improved our understanding of MM disease and have led to important clinical applications. In MM, GEP has been used to stratify patients, define risk, identify therapeutic targets, predict treatment response, and understand drug resistance. In this study, we built a gene risk score for 267 genes using RNA-seq data that demonstrated a prognostic value in two independent cohorts (n = 674 and n = 76) of newly diagnosed MM patients treated with high-dose Melphalan and autologous stem cell transplantation. High-risk patients were associated with the expression of genes involved in several major pathways implicated in MM pathophysiology, including interferon response, cell proliferation, hypoxia, IL-6 signaling pathway, stem cell genes, MYC, and epigenetic deregulation. The RNA-seq-based risk score was correlated with specific MM somatic mutation profiles and responses to targeted treatment including EZH2, MELK, TOPK/PBK, and Aurora kinase inhibitors, outlining potential utility for precision medicine strategies in MM.