A Single-Cell Transcriptional Analysis of Tumour Cells and the Immune Microenvironment during Disease Evolution in a Transgenic Mouse Model of Myeloma
Authors: Croucher Danielle C , Chesi Marta , Li Zhihua , Garbitt Victoria Marie , Sharik Meaghen E , Waller Daniel , Sebag Michael , Bergsagel P. Leif , Pugh Trevor J. , Trudel Suzanne 
Published: October 2019

Abstract
Introduction: Multiple Myeloma (MM) is consistently preceded by pre-malignant asymptomatic monoclonal gammopathies (AMG). To date, our understanding of the pathogenesis of progression to MM remains incomplete. Genetic analyses of AMG cells compared to MM-derived plasma cells (PCs) have found few differences, suggesting that progression may be mediated in part by tumour-extrinsic mechanisms. To comprehensively examine the cellular and molecular complexities of MM pathogenesis, we performed an unbiased single cell RNA-sequencing (scRNA-seq) analysis of tumour cells as well as immune cells from the tumour microenvironment (TME) derived from transgenic mice transitioning from AMG to MM.
Methods: We employed the Vk*MYC immune-competent mouse model of MM (C57BL/6/KaLwRij), which is a clinically and biologically faithful model of untreated disease that similarly progresses from AMG to MM with age. We established an age-based cohort of Vk*MYC mice to recapitulate a range of MM disease stages, generated single-cell suspensions from flushed bone marrow and subjected these cells to scRNA-seq profiling (10x Genomics).
Results: Across 12 samples profiled to date, our scRNA-seq dataset contains 82,853 high-quality cells, expressing 17,922 genes. We employed dimensionality reduction and unsupervised graph-based clustering to visualize and group transcriptionally-similar cell populations, which revealed 42 clusters. Expression of known marker genes and computed correlation scores with bulk gene expression reference datasets enabled annotation of cell types, revealing both malignant cells and non-malignant immune cell populations. We first focused on single cell T/NK profiles in our data given the emerging utility of immune checkpoint inhibitors that target these populations. Although we did not observe numerical differences in the proportion of CD8+ T cells across disease stages, analysis of immune checkpoint receptor genes revealed increased expression of Pdcd1 (PD-1) and Lag3 in CD8+ T cells from mice with disease. Co-expression of LAG3 and PD-1 proteins was also confirmed using a Vk*MYC transplantable model, with a positive correlation between disease burden (%CD138+/B220- cells) and %PD1+LAG3+ CD8+ T cells by flow cytometry. Consistent with reports of PD-1 and LAG3 co-expression on non-functional exhausted T cells, CD8+ T cells from diseased mice demonstrated elevated T cell exhaustion scores in our scRNA-seq dataset. These observations suggest that T cell exhaustion may be mediated by multiple immune checkpoint receptors during disease evolution. Although combinatorial treatment with PD-1 and LAG3 antibodies failed to induce tumour regression in mice with established disease, the addition of cyclophosphamide (Cy) to these antibodies resulted in marked improvement in survival of the mice compared to Cy alone, presumably by promoting immunogenic cell death. Studies exploring the combination of LAG3 and PD-1 antibodies as a strategy to inhibit transition from AMG to MM in the Vk*MYC mice are ongoing and will be reported. We also performed subclustering analysis of 5,228 Sdc1+ (CD138) PCs in our scRNA-seq dataset revealing 11 distinct clusters, with evidence of inter- and intra-tumoural heterogeneity across all Vk*MYC mice. Differential gene expression analysis revealed a non-malignant PC (nPC) cluster as supported by lower Myc transgene and Ccnd2 expression. Moreover, this cluster was predominantly comprised of cells from age-matched control mice or mice with earlier disease. Single-cell chromosomal copy number analysis revealed loss of Chr5 in the majority of tumour cells from MM mice, but not in the nPC cluster. Loss of Chr5 was observed in tumor subclones from all AMG mice suggesting that it is an early and potentially unifying event in Vk*MYC mice during disease progression. Further, the data support the establishment of intratumoural heterogeneity early in disease evolution.
Conclusions: Our approach of using scRNA-seq to characterize the pathogenesis of disease evolution in MM has enabled simultaneous measurement of intratumoural heterogeneity and immune cell phenotypes in the TME. In turn, this has provided insights into mechanisms that may contribute to transition from AMG to MM, including induction of T cell exhaustion and loss of mouse Chr5. Ongoing and future work aims to evaluate whether these mechanisms can be exploited therapeutically in pre-malignant AMG.

Disclosures
Sebag: Amgen Canada: Membership on an entity's Board of Directors or advisory committees; Janssen Inc.: Membership on an entity's Board of Directors or advisory committees; Celgene Canada: Membership on an entity's Board of Directors or advisory committees; Takeda Canada: Membership on an entity's Board of Directors or advisory committees. Pugh:Prosigna: Honoraria; N/A: Patents & Royalties: Hybrid-capture sequencing for determining immune cell clonality; N/A: Patents & Royalties: Combined hybrid-capture DNA sequencing for disease detection; Boehringer Ingelheim: Research Funding; Chrysalis Biomedical Advisors: Honoraria; Merck: Honoraria; DynaCare: Consultancy.