Identification of mutation gene prognostic biomarker in multiple myeloma through gene panel exome sequencing and transcriptome analysis in Chinese population.

Background: The 5-year survival rate of multiple myeloma (MM) in China is less than 40%, with considerable individual heterogeneity. Gene mutations are important predictive biomarkers that influence MM treatment decision. The aim of our study was to uncover the clinical significance of mutated genes in MM in the Chinese population.

Methods: Targeted exon panel sequencing was performed of 400 genes to detect the gene mutation status in plasma cells from 50 patients with MM. DAVID was used to explore the functions and pathways of mutated genes. Detection of mutant gene expression, prognosis and immune cell infiltration with GSE6477. GEO2R was utilized to identify differentially expressed genes (DEGs). Kaplan-Meier and CIBERSORT were applied to compare survival distributions and evaluate the gene expression associated with immune cell infiltration, respectively.

Results: Mutations of 337 genes were identified in MM. The mutation types included SNP, INS, and DEL, but the dominant mutation type was SNP. Function and pathway analysis of mutant genes were performed to elucidate DNA modifications. We identified a total number of 660 downregulated and 587 upregulated genes from the GSE6477 dataset. Thirty-three common genes were present in both the mutant genes and DEGs. The functions and pathways of the mutated genes were enriched in myeloid cell differentiation, regulation of hemopoiesis, etc. Moreover, we found that the low expression of BCL6, BIRC3, HLA-DQA1, and VCAN was correlated with poor prognosis in MM.

Conclusions: The mutations and low expression of BCL6, BIRC3, HLA-DQA1, and VCAN were correlated with poor prognosis and immune cell infiltration in MM. This study is the first to reveal the spectrum of mutations in the Chinese population by the use of an NGS panel.