P-205: Identification of therapy induced evolution in MRD clones and resistance pathways in multiple myeloma through single-cell sequencing

Introduction: In multiple myeloma (MM), the leading cause of disease progression, exhibits intratumoral heterogeneity that enables adaptability, limits therapeutic success, and remains incompletely understood. Previous studies have confirmed that clonal evolution frequently occurs at disease follow-up, and the patterns of clonal evolution between diagnosis and relapse correlates with patient outcomes. However, it remains unclear whether clonal evolution occurs as early as after induction therapy.

Methods: In this study, we analyzed single-cell RNA sequencing from 19 patients with newly diagnosed MM (NDMM), 17 paired bone marrow (BM) samples after 2-4 cycles of proteasome inhibitor (PI) plus immunomodulatory drugs (IMiDS)-based induction therapy were collected and underwent scRNA seq. BM samples from one healthy donor was also collected. By tracing transcriptional and cytogenetic PCs clones over time and performing differential expression analysis, we defined different patterns of clonal evolution and identified potential resistant pathways.

Results: In our cohort, 12 patients were identified with more than 30 PCs in the BM samples sequenced after induction therapy. Among them, our analysis of the transcriptional and cytogenetic clonal dynamics in MM patients identified three main trajectories: patients with sensitive PC clones or clone that responds to treatment with >90% of the malignant PCs replaced with seemingly health PCs (3/12 patients, 25%); in patients with a resistant clone that did not respond or only very partly responded to treatment, > 50% of the malignant PCs were found in the BM post-treatment (5/12 patients, 42%); and patients with clonal selection, indicating that the major clone has been replaced by a small or undetectable clone at baseline (4/12 patients, 33%). For these four patients identified with clonal selection, one was observed with branching evolution, and the other three were observed with differential evolution.

Transcriptional differences among sensitive clones, resistant clones and selective clones were detected based on a pairwise comparison of the gene expressions. A large number of differentially expressed genes with reported MM resistant-related functions were observed in the resistant clones, including previously reported 1q-related genes such as CKS1B, HNRNPU and H3F3A; unfolded protein-stress response-related genes such as NOP56 and HSP90B1; cell cycle- and cell proliferation-related genes such as TUBA1B, STMN1 and HMGB2. For selective clones, an evident activation of NF-kB signaling pathway was observed. Furthermore, an increase in CD74-MIF cellular interaction was found between selective clones and the tumor microenvironment.

Conclusions: Together, our study confirms that clonal dynamics of the evolving PC clones may occur early after upfront therapy, and reveals that the acquisition of therapeutic resistant pathways is associated with early adaptation to treatment.