- Open Access
Single-cell transcriptome analysis reveals stem cell-like subsets in the progression of Waldenström’s macroglobulinemia
Experimental Hematology & Oncology volume 12, Article number: 18 (2023)
Waldenström’s macroglobulinemia (WM) is an uncommon lymphoproliferative disorder, and the precise cellular landscape and the mechanisms of progression from IgM monoclonal gammopathy of undetermined significance (MGUS) to WM remain unclear. We performed single-cell RNA sequencing of CD19 + and CD19-CD38 + cells from healthy donors, IgM MGUS and WM patients. We found that samples from IgM MGUS and WM patients were composed of fewer early B-cell subsets and more T cells and NK cells than those from healthy controls. Compared with those of IgM MGUS patients, mature B cells of WM patients showed upregulation of HES1, GADD45B, NEAT1, DUSP22, RGS1, RGS16, and PIM1. We also identified a subpopulation of CD3 + CD19 + cells in IgM MGUS and WM patients, and trajectory analysis suggested that this subpopulation might be a stem cell-like subset. Further targeted gene sequencing of CD3 + CD19 + and CD3-CD19 + cells proved that MYD88 might be the early events in tumorigenesis according to variant allele fraction analysis. Additional subclonal hits such as CXCR4 and MAP2K1 mutations could be acquired during tumor progression. CXCL signaling, CCL signaling, IL2 signaling and TGFβ signaling pathways were involved in communication between CD3 + CD19 + cells and other immune cells. Our findings reveal the composition of CD38 + immune microenvironment together with B cells and plasma cells in IgM MGUS and WM patients, and provide comprehensive insights into mechanisms of progression from IgM MGUS to WM. The rare CD3 + CD19 + cells might be cells with “stemness” feature.
To the editor,
Waldenström’s macroglobulinemia (WM) is an uncommon lymphoproliferative disorder characterized by monoclonal immunoglobulin M protein production and bone marrow infiltration by lymphoplasmacstetytoid cells . It remains an incurable disease. IgM monoclonal gammopathy of undetermined significance (MGUS) is the premalignant condition and is associated with the risk of developing WM . MYD88L265P is a recurrent mutation found in more than 90% of WM patients, and this genetic event is detectable in at least half of IgM MGUS patients . However, the precise cellular landscape and the mechanisms of progression from IgM MGUS to WM remain unclear. Single-cell RNA sequencing (scRNA-seq) provides us a powerful approach to explore tumor heterogeneity and identify rare malignant cells . Anti-CD38 monoclonal antibodies have been introduced into the therapeutic arsenal for plasma cell diseases . Therefore, we analyzed CD38 + immune microenvironment together with B cells and plasma cells to illustrate the first cellular landscape of IgM MGUS and WM in single cell resolution.
We isolated bone marrow CD19 + and CD19-CD38 + cells from WM patients(n = 3), IgM MGUS patients(n = 3) and healthy donors(n = 3) and employed the 10 × Genomics platform to perform single-cell transcriptomic sequencing. After quality control, our dataset contained a total of 73,024 cells. We used uniform manifold approximation and projection (UMAP) to visualize the cell superclusters and identified clusters including B cells (CD19, MS4A1, and CD79A), plasma cells (SDC1, CD38, and PRDM1), CD3 + CD20 + CD19 + cells, T cells (CD3D, CD8A, and CD4) and NK cells (NCAM1, GNLY, and NKG7) based on the expression of canonical lineage markers and cluster-specific markers (Fig. 1A and Additional file 2: Figure S1A, B). We found that healthy donors had the highest percentage of B cells (72.7% in healthy donors, 19.0% in IgM MGUS patients, 60.5% in WM patients; p < 0.001) and that IgM MGUS and WM patients had more NK and T cells than healthy controls (0.06% and 17.7%, respectively, in healthy donors; 24.1% and 20.5% in IgM MGUS patients; 8.5% and 24.0% in WM patients; p < 0.001) (Fig. 1B). We also noticed that CD3 + CD20 + (CD19 +) cells were present in both IgM MGUS and WM patients (Additional file 2: Figure S1C–E).
First, we investigated the cell clusters correlated with WM pathology including B cells and plasma cells. We identified pre-B cells, pro-B cells and mature B cells based on the progenitor cell marker CD34, pre-BCR genes IGLL1 and VPREB1, and B-lymphoid lineage surface markers MME, CD38, FCER2 and CD27. We compared the B-cell gene expression profiles after removing the early B cells (Additional file 2: Figure S1F). Compared with those of IgM MGUS patients, mature B cells of WM patients showed upregulation of HES1, GADD45B, NEAT1, DUSP22, RGS1, RGS16, and PIM1 (Fig. 1C and Additional file 3: Table S1). Light chain restriction and higher copy number variants (CNV) levels were found in B cells of WM patients (Fig. 1C, D). SCENIC analysis showed that FOS, FOSB, EGR1 and JUNB were downregulated in B cells from WM and IgM MGUS patients (Fig. 1E). REL, a proto-oncogene promoting the survival and proliferation of B lymphocytes, was downregulated in B cells of WM patients, and IRF8 and ELF1 were upregulated in B cells of IgM MGUS patients (Fig. 1E). Nevertheless, plasma cells of WM patients, IgM MGUS patients and healthy donors were colocalized in one cluster, which might indicate their similar gene expression patterns. The significant differentially expressed genes (DEGs) across the groups are highlighted in scatter plots (Fig. 1F) and showed in Additional file 4: Table S2. In particular, we noticed that the expression levels of part of genes decrease in IgM MGUS patients and then rise again when the disease progresses to WM. We assumed that the regulation might be correlated with the progression from IgM MGUS to WM.
Next, we further determined the characteristics of the rare CD3 + CD19 + cells. Trajectory analysis revealed that CD3 + CD20 + (CD19 +) cells and pre-B cells were located at the same end of the branched structure and were then directed toward mature B cell and plasma cell fates, respectively (Fig. 2A). This inferred developmental trajectory suggested that CD3 + CD20 + (CD19 +) cells may be tumor stem cell-like subset. Using flow cytometry, we validated that CD3 was expressed in a fraction of CD19 + cells (mean: 3.69%, range: 0.28%-15.10%) of bone marrow samples from 8 WM patients (Fig. 2B, Additional file 2: Figure S2A and Additional file 5: Table S3). We compared the gene expression of CD3 + CD20 + CD19 + cells between WM patients and IgM MGUS patients (Additional file 6: Table S4). Genes including IGKC, JCHAIN, AREG, RGS1, IGHM were significantly upregulated in CD3 + CD20 + CD19 + cells of WM patients. In addition, targeted gene sequencing (TGS) was performed on CD3 + CD19 + and CD3-CD19 + cells from 6 WM patients (Additional file 7: Table S5). MYD88 were proved to be the most frequently mutated gene in this cell subset (Fig. 2C). And we inferred that MYD88 mutation might be the early events in tumorigenesis by variant allele fraction (VAF) analysis (Fig. 2D, E). Additional subclonal hits, such as CXCR4 and MAP2K1 mutations, could be acquired during tumor progression.
We observed higher proportions of T lymphocytes and NK cells in WM patients; thus, CellChat was employed to explore the cell–cell communication between immune cells and CD3 + CD20 + (CD19 +) cells in WM patients. We observed that the CXCR4-CXCL12 ligand-receptor pair was enriched in the interaction between CD3 + CD20 + cells and B cells (Additional file 2: Figure. S2B, C). The CCL signaling pathway were highly expressed in CD8 + T cell/NK-cell—CD3 + CD20 + cell interactions (Additional file 2: Figure. S2B, C). In addition, IL2 signaling, as well as TGF β signaling, were also involved in communication networks between CD3 + CD20 + cells and T/NK cells (Additional file 2: Figure. S2B, C).
In conclusion, our study provides comprehensive insights into mechanisms of progression from IgM MGUS to WM. We identified the rare CD3 + CD19 + cell subpopulation in WM patients. It will be interesting to explore novel therapeutic strategies targeting rare potential cells with “stemness” in future.
Monoclonal gammopathy of undetermined significance
Single-cell RNA sequencing
Uniform manifold approximation and projection
Targeted gene sequencing
Differentially expressed genes
Copy number variants
Variant allele fraction
Morton LM, Turner JJ, Cerhan JR, Linet MS, Treseler PA, Clarke CA, et al. Proposed classification of lymphoid neoplasms for epidemiologic research from the pathology working group of the international lymphoma epidemiology consortium (InterLymph). Blood. 2007;110(2):695–708.
Kyle RA, Therneau TM, Rajkumar SV, Remstein ED, Offord JR, Larson DR, et al. Long-term follow-up of IgM monoclonal gammopathy of undetermined significance. Blood. 2003;102(10):3759–64.
Jiménez C, Sebastián E, Chillón MC, Giraldo P, Mariano Hernández J, Escalante F, et al. MYD88 L265P is a marker highly characteristic of, but not restricted to Waldenström’s macroglobulinemia. Leukemia. 2013;27(8):1722–8.
Peng J, Sun BF, Chen CY, Zhou JY, Chen YS, Chen H, et al. Single-cell RNA-seq highlights intra-tumoral heterogeneity and malignant progression in pancreatic ductal adenocarcinoma. Cell Res. 2019;29(9):725–38.
Nooka AK, Kaufman JL, Hofmeister CC, Joseph NS, Heffner TL, Gupta VA, et al. Daratumumab in multiple myeloma. Cancer. 2019;125(14):2364–82.
Database Resources of the National Genomics Data Center. China national center for bioinformation in 2022. Nucleic Acids Res. 2022;50(D1):D27-d38.
Chen T, Chen X, Zhang S, Zhu J, Tang B, Wang A, et al. The genome sequence archive family: toward explosive data growth and diverse data types. Genom Proteom Bioinform. 2021;12:279.
We appreciate Springer Nature Author Services for English language editing.
This work was supported by institutional research funding provided by National Key R&D Program of China, Ministry of Science and Technology of the People’s Republic of China (2022YFC2304605), the National High Level Hospital Clinical Research Funding (2022-PUMCH-B-046, 2022-PUMCH-A-193 and 2022-PUMCH-A-080).
Ethics approval and consent to participate
Written Informed consent was obtained from all patients, and the study was approved by the Peking Union Medical College Hospital Ethics Committee. The current study was performed in accordance with the ethical standards of the 1964 Declaration of Helsinki and its later amendments.
Consent for publication
All authors have read and approved the submitted manuscript.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Demographics of patients. Table S2. Summary and quality of sequencing. Table S3. Gene list of targeted gene sequencing panel
Additional file 2: Figure S1.
Cellular landscape of patients with WM and IgM MGUS. A. UMAP plot of merged single-cell transcriptomes. Cells are colored according to sample origin. B. Dot plot of feature genes expression in each cluster. C. UMAP plot showing clusters identified in WM patients. D. UMAP plot showing clusters identified in healthy donors. E. UMAP plot showing clusters identified in IgM MGUS patients. F. UMAP plot of B cell sub-clusters in patients and healthy donors. Figure S2. Identification of CD3+CD19+ cells and cell-cell communication. A. Flow cytometric analysis showing the population of CD3+CD19+ cells in six WM patients. B. Cell-cell communication inferred by Cellchat. C.Relative contribution of each ligand-receptor pair.
Additional file 3.
Differentially expressed genes of mature B cells.
Additional file 4.
Differentially expressed genes of plasma cells.
Additional file 5: Table S3.
Percentage of CD3+CD19+ cells in 8 WM samples.
Additional file 6.
Differentially expressed genes of CD3+CD19+ cells between WM and IgM MGUS patients.
Additional file 7
: Table S5. Mutation in the 6 WM patients
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Qiu, Y., Wang, Xs., Yao, Y. et al. Single-cell transcriptome analysis reveals stem cell-like subsets in the progression of Waldenström’s macroglobulinemia. Exp Hematol Oncol 12, 18 (2023). https://doi.org/10.1186/s40164-023-00382-6