A Single-Cell Transcriptional Analysis of Tumour Cells and the Immune Microenvironment during Disease Evolution in a Transgenic Mouse Model of Myeloma
Bibliographic record
Abstract
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.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.000 | 0.000 |
| Meta-epidemiology (narrow) | 0.000 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.001 |
| Bibliometrics | 0.001 | 0.001 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.001 | 0.000 |
| Open science | 0.000 | 0.000 |
| Research integrity | 0.000 | 0.001 |
| Insufficient payload (model declined to judge) | 0.001 | 0.001 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".