Integrating Phosphoproteomics and Transcriptional Classifiers Reveals "Hidden Signaling" in Multiple Myeloma Including Differential KRAS and NRAS Mutant Effects
Authors: Lin Yu-Hsiu T. , Way Gergory P. , Mariano Margarette C. , Marcoulis Makeba , Ferguson Ian , Greene Casey S. , Wiita Arun P 
Published: October 2019

Abstract
Introduction: Multiple myeloma (MM) is a complex disease that requires a sophisticated treatment strategy. Currently, no kinase inhibitors have been approved for MM despite their potential for supplementing current combination therapies. Previous functional studies have explored kinase dependency in MM by either a small molecule inhibitor library (Dhimolea et al. 2014 ASH) or RNA interference (Tiedemann et al. 2010, Blood 115:1594). However, owing to their off-target effects, these approaches are imprecise at dissecting signaling networks driving MM growth and survival. Here, we aim to improve diagnostic and prognostic measures and recommend small molecule-based treatments for MM patients by identifying vulnerable signaling patterns in disease using integrated transcriptome- and phosphoproteome-based predictive models.
Methods: We employed two methods for measuring cellular signaling activity within a tumor sample. The first involves an unbiased phosphoproteome profiling of eight human myeloma cell lines (HMCL) in biological triplicates. Cells were lysed and digested with trypsin prior to enrichment for phosphorylated peptides by Fe3+-IMAC. Samples were analyzed on a Thermo Q-Exactive Plus mass spectrometer and data processed in MaxQuant to identify and quantify phosphorylation sites. To infer relative kinase activities, we applied kinase set enrichment analysis (KSEA). Drug screening was performed in 384-well plates with CellTiter-Glo as viability readout as previously described (Lam et al. 2018, Haematologica 103:1218). For transcriptome analysis, we implemented the gene expression-based signaling pathway prediction model called PROGENy (Schubert et al. 2018, Nat Comm. 9:20). We applied PROGENy to RNAseq data on 64 MM cell lines (www.keatslab.org) as well as plasma cells from >1000 patients in the MMRF CoMMpass study (research.themmrf.org). Furthermore, using our recently described approach (Way et al. 2018, Cell Rep. 23:172), we built a machine learning classifier to predict RAS genotype from the transcriptomic profiles of MM patients. A tenth of the data set was withheld for testing while the rest was used to train the multiclass logistic regression classifier with sparse penalty.
Results: We measured the KSEA-predicted activities of 297 kinases across eight tested MM cell lines. Initially, we were surprised to find high predicted activity in the Ras signaling pathway for KRAS-codon 12 mutant cell lines but low predicted activity in NRAS-mutant cell lines (Fig. A: AMO1 harbors a non-canonical KRAS mutation at codon 146). We further explored this finding with our machine learning-based Ras classifier built on transcriptional data in the CoMMpass study. We identified 311, 405, and 390 genes whose expressions are characteristic of the WT RAS, KRAS mutant, and NRAS mutant genotype, respectively, with surprisingly limited overlap between KRAS and NRAS transcriptional signatures. Building on our KSEA analysis, we next performed a kinase inhibitor screen to evaluate the predictive value of the inferred kinase activities for drug sensitivity. Of 12 screened compounds, mTOR inhibitor INK128 displayed the strongest correlation between drug response and predicted kinase activity. Furthermore, we probed the potential of using pathway activity signatures as prognostic and therapeutic markers. To this end, we applied PROGENy to RNAseq data derived from CoMMpass patients and found that the MAPK signature stratifies patient survival with statistical significance, while the presence and absence of RAS mutations carry no prognostic value (Fig. B). Finally, by integrating RNAseq and drug screen data from the Cancer Dependency Map, we identified three compounds whose inhibitory effects strongly correlate with MAPK activity scores while no significant difference in drug sensitivity was detected between RAS WT and mutants.
Conclusion: Both phosphoproteomics and a machine learning-based transcriptional classifier highlight a striking difference in the pattern of signaling between NRAS and KRAS mutants. In addition, we have demonstrated that PROGENy scores possess clinical value for prognostic and therapeutic use based on patient transcriptome data. Taken together, uncovering the cellular signaling networks dysregulated in MM may lead to improved precision medicine, particularly in stratifying patients who may benefit most from kinase inhibitor therapy.
Figure. Figure.

Disclosures
Wiita: TeneoBio: Research Funding; Sutro Biopharma: Research Funding.