Single Cell Profiling Unravels Plasma Cell Gene Expression Signatures Predictive of Response to Daratumumab Based Therapy, and Prognosis in Multiple Myeloma

Background: Multiple Myeloma (MM) is the second most common hematologic malignancy and remains incurable. Daratumumab (dara) is a potent anti-CD38 monoclonal antibody used for MM treatment, but responses are heterogeneous, and resistance is inevitable. There is hence an urgent need for biomarkers to predict the response to dara. We hypothesized that gene expression profiles (GEP) of tumor subpopulations could be used to predict the response to dara-based therapy. Methodology: We used the Enhanced Single-Cell Analysis with Protein Expression (ESCAPE) RNA-Seq platform (Singleron Biotechnologies) to simultaneously profile the gene and protein expression of plasma cells from the bone marrow of 15 MM patients prior to treatment with dara-based regimens. These data were combined with publicly available single-cell data from Cohen et al. 2021, resulting in a total of 32 samples. Response was defined based on the international myeloma working group (IMWG) criteria, considering patients achieving a very good partial response (VGPR) or better as “responders,” while those achieving a partial response (PR) or worse were deemed “non-responders”. The data were split into 23 training samples (15 responders and 8 non-responders) and 9 validation samples (3 responders and 6 non-responders), each set containing a combination of cases from the two sources of single-cell data. An interpretable machine-learning model was then built to predict the response of each individual cell. The median response score for each sample was used to predict the patient-level response. The Shapley Additive Explanations (SHAP) values from the dara response prediction model were then used as a starting point for creating a general prognostic model. The 200 most influential genes from the dara response prediction model were used to create a linear response model based on the TCGA (MMRF-Compass) bulk RNAseq dataset. The bulk RNAseq dataset was split 80:20 into training and test datasets. Results from the test dataset are presented here. Results: We analyzed single-cell data from a total of 13,094 plasma cells from 32 BM samples. These included 17 cases of relapsed MM described by Cohen et al., who were treated with dara, carfilzomib, lenalidomide, and dexamethasone. The patients treated at our centre comprised 7 cases of relapsed MM treated with dara, thalidomide, and dexamethasone, and 8 samples from patients with newly diagnosed MM treated with dara, bortezomib, and dexamethasone on clinical tria