Hub genes and associated drugs for multiple myeloma with 1q21+: identified by bioinformatic analysis

While 1q21+ was common genetic alteration and found to have adverse effect on prognosis, the underlying genes remain unclear. Identification of related genes may provide additional help for rational intervention. The microarray dataset GSE2658 associated with MM was downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were obtained, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to annotate their functions. The hub genes were derived from the combined results of up-regulated DEGs and weighted gene coexpression network analysis (WGCNA). The receiver operating characteristic (ROC) curves of hub genes were plotted to evaluate correlation with 1q21+. Survival analysis and drug-gene interaction of hub genes were performed separately to find the prognostic value and potential targeted drugs. A total of 55 DEGs were identified. GO and KEGG pathway analyses suggested that the DEGs were related to several pathways of cell proliferation. NVL, IL6R, DUSP23 were proven to be highly correlated with 1q21+ and have adverse effects on prognosis. IL6R, DUSP23 were matched to known interaction-drug. This study revealed potential roles of hub genes in the pathogenesis and progression of MM with 1q21+, further investigations are needed to elucidate the mechanisms.