Enhancing the efficacy of proteasome inhibitors (PI) is a central goal in myeloma therapy. We proposed that signaling-level responses after PI may reveal new mechanisms of action that can be therapeutically exploited. Unbiased phosphoproteomics after treatment with the PI carfilzomib surprisingly demonstrates the most prominent phosphorylation changes on splicing related proteins. Spliceosome modulation is invisible to RNA or protein abundance alone. Transcriptome analysis after PI demonstrates broad-scale intron retention, suggestive of spliceosome interference, as well as specific alternative splicing of protein homeostasis machinery components. These findings lead us to evaluate direct spliceosome inhibition in myeloma, which synergizes with carfilzomib and shows potent anti-tumor activity. Functional genomics and exome sequencing further support the spliceosome as a specific vulnerability in myeloma. Our results propose splicing interference as an unrecognized modality of PI mechanism, reveal additional modes of spliceosome modulation, and suggest spliceosome targeting as a promising therapeutic strategy in myeloma.

While advancements in genome sequencing have identified millions of somatic mutations in cancer, their functional impact is poorly understood. We previously developed the expression-based variant impact phenotyping (eVIP) method to use gene expression data to characterize the function of gene variants. The eVIP method uses a decision tree-based algorithm to predict the functional impact of somatic variants by comparing gene expression signatures induced by introduction of wild-type (WT) versus mutant cDNAs in cell lines. The method distinguishes between variants that are gain-of-function, loss-of-function, change-of-function, or neutral. We present eVIP2, software that allows for pathway analysis (eVIP Pathways) and usage with RNA-seq data. To demonstrate the eVIP2 software and approach, we characterized two recurrent frameshift variants in RNF43, a negative regulator of Wnt signaling, frequently mutated in colorectal, gastric, and endometrial cancer. RNF43 WT, RNF43 R117fs, RNF43 G659fs, or GFP control cDNA were overexpressed in HEK293T cells. Analysis with eVIP2 predicted that the frameshift at position 117 was a loss-of-function mutation, as expected. The second frameshift at position 659 has been previously described as a passenger mutation that maintains the RNF43 WT function as a negative regulator of Wnt. Surprisingly, eVIP2 predicted G659fs to be a change-of-function mutation. Additional eVIP Pathways analysis of RNF43 G659fs predicted 10 pathways to be significantly altered, including TNF-α via NFκB signaling, KRAS signaling, and hypoxia, highlighting the benefit of a more comprehensive approach when determining the impact of gene variant function. To validate these predictions, we performed reporter assays and found that each pathway activated by expression of RNF43 G659fs, but not expression of RNF43 WT, was identified as impacted by eVIP2, supporting that RNF43 G659fs is a change-of-function mutation and its effect on the identified pathways. Pathway activation was further validated by Western blot analysis. Lastly, we show primary colon adenocarcinoma patient samples with R117fs and G659fs variants have transcriptional profiles similar to BRAF missense mutations with activated RAS/MAPK signaling, consistent with KRAS signaling pathways being GOF in both variants. The eVIP2 method is an important step towards overcoming the current challenge of variant interpretation in the implementation of precision medicine. eVIP2 is available at https://github.com/BrooksLabUCSC/eVIP2.

While advancements in genome sequencing have identified millions of somatic mutations in cancer, their functional impact is poorly understood. We previously developed the expression-based variant impact phenotyping (eVIP) method to use gene expression data to characterize the function of gene variants. The eVIP method uses a decision tree-based algorithm to predict the functional impact of somatic variants by comparing gene expression signatures induced by introduction of wild-type versus mutant cDNAs in cell lines. The method distinguishes between variants that are gain-of-function, loss-of-function, change-of-function, or neutral. We present eVIP2, software that allows for pathway analysis (eVIP Pathways) and usage with RNA-seq data. To demonstrate the eVIP2 software and approach, we characterized two recurrent frameshift variants in RNF43, a negative regulator of Wnt signaling, frequently mutated in colorectal, gastric and endometrial cancer. RNF43 WT, RNF43 R117fs, RNF43 G659fs, or GFP control cDNA were overexpressed in HEK293T cells. Analysis with eVIP2 predicted that the frameshift at position 117 was a loss-of-function mutation, as expected. The second frameshift at position 659, was, surprisingly, predicted to be a gain-of-function mutation.Additional eVIP Pathways analysis of RNF43 G659fs predicted 10 pathways to be significantly altered, including TNF alpha via NFKB signaling, KRAS signaling, and hypoxia. To validate these predictions, we performed reporter assays and found that all eVIP2 impactful pathways tested in the assay were activated by expression of RNF43 G659fs, but not by expression of RNF43 WT, supporting that RNF43 G659fs is a gain-offunction mutation and its effect on the identified pathways. The eVIP2 method is an important step towards overcoming the current challenge of variant interpretation in the 3 implementation of precision medicine. eVIP2 is available at https://github.com/BrooksLabUCSC/eVIP2.

No abstract