In-depth analysis of alternative splicing landscape in multiple myeloma and potential role of dysregulated splicing factors

Abstract

Splicing changes are common in cancer and are associated with dysregulated splicing factors. Here, we analyzed RNA-seq data from 323 newly diagnosed multiple myeloma (MM) patients and described the alternative splicing (AS) landscape. We observed a large number of splicing pattern changes in MM cells compared to normal plasma cells (NPC). The most common events were alterations of mutually exclusive exons and exon skipping. Most of these events were observed in the absence of overall changes in gene expression and often impacted the coding potential of the alternatively spliced genes. To understand the molecular mechanisms driving frequent aberrant AS, we investigated 115 splicing factors (SFs) and associated them with the AS events in MM. We observed that ~40% of SFs were dysregulated in MM cells compared to NPC and found a significant enrichment of SRSF1, SRSF9, and PCB1 binding motifs around AS events. Importantly, SRSF1 overexpression was linked with shorter survival in two independent MM datasets and was correlated with the number of AS events, impacting tumor cell proliferation. Together with the observation that MM cells are vulnerable to splicing inhibition, our results may lay the foundation for developing new therapeutic strategies for MM. We have developed a web portal that allows custom alternative splicing event queries by using gene symbols and visualizes AS events in MM and subgroups. Our portals can be accessed at http://rconnect.dfci.harvard.edu/mmsplicing/ and https://rconnect.dfci.harvard.edu/mmleafcutter/.