Functional Impact of Genomic Complexity on the Transcriptome of Multiple Myeloma

Multiple myeloma is a biologically heterogenous plasma-cell disorder. In this study, we aimed at dissecting the functional impact on transcriptome of gene mutations, copy-number abnormalities (CNA), and chromosomal rearrangements (CR). Moreover, we applied a geno-transcriptomic approach to identify specific biomarkers for personalized treatments. We analyzed 514 newly diagnosed patients from the IA12 release of the CoMMpass study, accounting for mutations in multiple myeloma driver genes, structural variants, copy-number segments, and raw-transcript counts. We performed an drug sensitivity screen (DSS), interrogating the Cancer Dependency Map (DepMap) dataset after anchoring cell lines to primary tumor samples using the algorithm. Immunoglobulin translocations, hyperdiploidy and chr(1q)gain/amps were associated with the highest number of deregulated genes. Other CNAs and specific gene mutations had a lower but very distinct impact affecting specific pathways. Many recurrent genes showed a hotspot (HS)-specific effect. The clinical relevance of double-hit multiple myeloma found strong biological bases in our analysis. Biallelic deletions of tumor suppressors and chr(1q)-amplifications showed the greatest impact on gene expression, deregulating pathways related to cell cycle, proliferation, and expression of immunotherapy targets. Moreover, our DSS showed that not only t(11;14) but also chr(1q)gain/amps and inactivation predicted differential expression of transcripts of the BCL2 axis and response to venetoclax. The multiple myeloma genomic architecture and transcriptome have a strict connection, led by CNAs and CRs. Gene mutations impacted especially with HS-mutations of oncogenes and biallelic tumor suppressor gene inactivation. Finally, a comprehensive geno-transcriptomic analysis allows the identification of specific deregulated pathways and candidate biomarkers for personalized treatments in multiple myeloma.