Single-Cell RNA Sequencing of Circulating Tumor Cells in Precursor Myeloma Patients Reveals Mechanisms of Disease Dissemination
Authors: Lightbody Elizabeth D. , Sklavenitis-Pistofidis Romanos , Firer Danielle T. , Tsuji Junko , Agius Mike P. , Dutta Ankit K. , Wu Ting , Barr Hadley , Kim Sungjae , Coorens Tim , Alberge Jean-Baptiste , Aranha Michelle P. , Haradhvala Nicholas J , Kham Su Nang , Boehner Cody J. , Walsh Kelly , Rahmat Mahshid , Konishi Yoshinobu , Hevenor Laura , Towle Katherine , Horowitz Erica M. , Perry Jacqueline , Fleming Grace , El-Khoury Habib , Vinyard Michael , Cowan Annie , Auclair Daniel , Ready John E. , Marinac Catherine R. , Getz Gad , Ghobrial Irene M. 
Published: December 2023

Introduction
Multiple Myeloma (MM) is preceded by the precursors Monoclonal Gammopathy of Undetermined Significance (MGUS) and Smoldering Multiple Myeloma (SMM), where some patients with precursor disease progress to active MM faster than others. BM biopsies are useful to monitor disease progression, however, they cannot be repeated often for continuous monitoring of tumor burden especially in the precursor disease management setting. Numerous studies have identified circulating tumor cell (CTC) levels in peripheral blood (PB) are a powerful biomarker of disease aggressiveness; CTC levels have associated with disease stage, PFS and OS in various trials and even outperformed BM to assess tumor burden in the SMM setting. Genomic profiling also demonstrated CTCs harbor the same copy number abnormalities, translocations and mutations as BM tumor cells, and CTCs can be serially monitored to detect the emergence of high-risk subclones in the blood. Further investigations of the molecular profile of both normal plasma cells (NPCs) and CTCs are needed to improve our understanding of myeloma pathogenesis and mechanisms of CTC dissemination across the MM disease continuum.
Methods
A cohort of 309 samples including CD138-enriched paired PB and BM aspirates and matched CD138- BM aspirates were collected from 103 individuals, including patients with MGUS (n=11), SMM (n=46), NDMM (n=17) and healthy donors (n=29). Malignant PCs from BM and PB underwent 5' single-cell RNA sequencing (scRNA-seq) and single-cell B-cell receptor sequencing (scBCR-seq) (10x Genomics), and CD138- BM immune cells underwent 5' scRNA-seq and single-cell T-cell receptor sequencing (scTCR-seq) to study immune alterations related the tumors' ability to circulate. To differentiate malignant from normal PCs, we used clonal V(D)J rearrangements. Differential expression (DE) and composition analyses were conducted using Wilcoxon's rank-sum tests.
Results
We successfully captured and profiled 774,647 BM tumor cells and 93,878 CTCs from the PB. The percent malignant plasma cells within the PB increased with disease stage, where NDMM patients had significantly more CTCs compared to MGUS patients (p=0.017) and low-risk SMM patients (p=0.0052). A median of 13.51%, 15.53%, and 18.34% CTCs were present from low (n=20), intermediate (n=10), and high-risk (n=16) SMM patients as defined by the International Myeloma Working Group's 2/20/20 criteria, suggesting sequencing-based CTC enumeration captures prognostically relevant differences in tumor burden. High expression of driver genes upregulated in patients with translocations, including CCND1, NSD2, and MAF, were detected in both BM tumor cells and CTCs in patients with t(11;14), t(4;14), and t(14;16) as identified by fluorescence in situ hybridization (FISH). In 5 patients with normal or inconclusive FISH results, we observed high levels of CCND2 and MAF in both BM and CTCs, indicating scRNA-seq can detect missed prognostically relevant cytogenetic abnormalities. A cytogenetic classifier model was developed to determine the accuracy of translocation calling using minimal numbers of CTCs, where the classifier's performance was able to consistently identify translocations downsampled to 50 CTCs. DE analysis of CTCs vs. BM tumor cells highlighted transcriptional similarity between BM and CTCs, validating their utility as a surrogate for analyzing BM tumor cells. DE analysis also revealed 8 genes significantly upregulated and 3 genes significantly downregulated in CTCs compared to BM tumor cells, providing novel insights into genes involved in PC circulatory potential. Pathway enrichment analysis revealed genes upregulated in CTCs were associated with epithelial mesenchymal transition, interferon response, and inflammation, consistent with CTC studies on MM patients, suggesting these pathways are dysregulated earlier in the disease continuum.
Conclusions
In the largest scRNA-seq study on CTCs to date, we demonstrate the utility of CTC-based molecular profiling for prognostication of patients with early-stage disease and provide novel insights into PC circulatory potential. Additional analyses are ongoing to gain further insight into intra-patient CTC heterogeneity and define high-risk disease CTC signatures that emerge throughout the MM disease continuum.