Skip to main content

Single-cell RNA sequencing identifies the properties of myelodysplastic syndrome stem cells

Letter to editor,

Myelodysplastic syndromes (MDS) are heterogeneous clonal diseases characterized by cytopenia caused by ineffective hematopoiesis and high risk of transformation into acute myeloid leukemia (AML). At present, the pathogenesis of MDS has not been elucidated. MDS is a group of stem cell diseases. The abnormal proliferation and blockade in differentiation of hematopoietic stem cells (HSCs) result in cytopenia and leukemic transformation. HSCs architectures in MDS can also predict therapeutic reaction [1]. Therefore, the high-resolution analysis of HSCs is of great significance. We investigated the properties of MDS stem cells by single-cell RNA sequencing (scRNA-seq). Lineage negative (Lin) cells from bone marrow aspirates of 5 patients with MDS and 2 patients with secondary AML (sAML) were sorted out for scRNA-seq (Additional file 1: Additional methods, Additional file 2: Table S1). The scRNA-seq data of bone marrow mononuclear cells from 17 healthy donors (HDs) were downloaded from the Gene Expression Omnibus database (GSE120221) [2].

A total of 65,509 Lin hematopoietic cells were analyzed. After dimension reduction, clustering, and visualization, we identified HSC/multipotent progenitor (MPP) populations and conducted an in-depth analysis. We investigated the properties of MDS stem cells by analyzing the differentially expressed genes (DEGs) between patients and HDs. We found that genes associated with neutrophil granule, such as MPO, AZU1, DEFA3, had elevated expression in MDS patients compared with HDs (Fig. 1A and 1B). This is consistent with the tendency of predominantly myeloid differentiation trajectories of HSCs in MDS. We identified that the down-regulated DEGs in HSCs/MPPs were enriched in ribosome, translation, mRNA catabolic process, Th17 cell differentiation and antigen processing and presentation in MDS patients (Fig. 1C), which supported the opinion that MDS stem cells were in a relatively static state and had impaired immune function. The progression of MDS to AML usually bases on acquired mutations [3]. Our results suggested that abnormal proliferation, RNA metabolism and ribosome biogenesis also exist in MDS stem cells during leukemic transformation (Fig. 1D - F).

Fig. 1
figure 1

Analyses of aberrantly expressed genes and functional enrichment in myelodysplastic syndromes (MDS) stem cells by single-cell RNA sequencing

A Volcano plot of differentially expressed genes (DEGs) in hematopoietic stem cells (HSCs)/multipotent progenitors (MPPs) between MDS patients and healthy donors (HDs). B Functional enrichment bar chart of up-regulated DEGs in HSCs/MPPs of MDS patients compared with HDs. Red: Kyoto Encyclopedia of Genes and Genomes (KEGG) gene sets, blue: Biological process (BP). C Functional enrichment bar chart of down-regulated DEGs in HSCs/MPPs of MDS patients compared with HDs. D Volcano plot of DEGs in HSCs/MPPs between secondary acute myeloid leukemia (sAML) and MDS patients. Red plot: up-regulated in sAML patients, blue plot: up-regulated in MDS patients. E Functional enrichment bar chart of up-regulated DEGs in HSCs/MPPs of sAML patients compared with MDS patients. F Functional enrichment bar chart of down-regulated DEGs in HSCs/MPPs of sAML patients compared with MDS patients. G Gene sets enrichment analysis (GSEA) showed decreased expression of ribosomal genes in HSCs/MPPs of MDS patients compared with HDs. H The violin plots of expression of ribosomal genes in HSCs/MPPs of HDs and patients with lower-risk MDS (LR-MDS), higher-risk MDS (HR-MDS) and sAML. I The Cancer Genome Atlas (TCGA) database analysis showed the expression of several ribosomal genes in AML patients. J Quantitative real-time PCR (qPCR) analyses of RPL31 mRNA in CD34+ hematopoietic stem and progenitor cells from bone marrow of controls (n = 40), LR-MDS (n = 40), HR-MDS (n = 40) and sAML patients (n = 25). Significant difference is analyzed using analysis of variance. ***, P < 0.001. K Analyses of correlation between RPL31 mRNA and hemoglobin (HB), platelet (PLT) in MDS patients (n = 80)

It is noteworthy that ribosomal genes were widely down-regulated in HSCs/MPPs of MDS/sAML patients compared with those of HDs (Fig. 1G and H), which suggested that abnormal ribosome biogenesis may be involved in the pathogenesis of MDS and leukemic transformation. The Cancer Genome Atlas (TCGA) database analysis also showed decreased expression of some ribosomal genes in AML (Fig. 1I). The results of quantitative real-time PCR verification showed that the expression of RPL31 and RPL21 mRNA in CD34+ cells of patients with LR-MDS (n = 40), HR-MDS (n = 40) and sAML (n = 25) were significantly lower than those of control (n = 40) (P < 0.001) (Fig. 1J, Additional file 3: Fig. S1A) and positively correlated with the levels of hemoglobin and platelet (Fig. 1K, Additional file 3: Fig. S1B).

Hematopoietic homeostasis depends on the balance between self-renewal and differentiation of HSCs. Increasing genome data showed that abnormal expression of some crucial regulatory genes may break the homeostasis and eventually lead to hematological malignant diseases. Our results revealed MDS stem cells have abnormal proliferation, ribosome biogenesis, RNA metabolism and impaired immune function. The expression of ribosomal genes was down-regulated in MDS stem cells. Abnormal ribosome biogenesis and regulation of translation was frequent in myeloid diseases. Gene mutations of ribosomal protein could lead to Diamond-Blackfan Anemia, Schwachman Diamond Syndrome, Congenital Dyskeratosis and some other diseases, all of which have defects in hematopoiesis. In addition, the dyserythropoiesis in MDS with isolated del(5q) is related to the loss of heterozygosity of RPS14 gene. Haploinsufficiency of RPS14 leads to activation of innate immune system and p53 pathway, excessive apoptosis of erythroblasts and megaloblastic anemia. Study by Saha et al.[4] showed that the expression of ribosomal genes in leukemia stem cells (LSCs) was down-regulated compared with normal HSCs, which was consistent with the results in our study. The biosynthesis of proteins in HSCs is lower than that in committed progenitor cells and differentiated cells [5]. LSCs may have similar mechanisms to survive and resist to chemotherapy. Therefore, abnormal ribosome biogenesis is an important feature during development and progression of MDS and may be also related to therapeutic resistance.

Availability of data and materials

The raw data reported in this study are deposited in the NCBI Sequence Read Archive under bioproject No.PRJNA720840.

Abbreviations

DEG:

differentially expressed gene

HD:

healthy donor

HSC:

hematopoietic stem cell

HR-MDS:

higher-risk myelodysplastic syndrome

LSC:

leukemia stem cell

Lin :

Lineage negative

LR-MDS:

lower-risk myelodysplastic syndrome

MPP:

multipotent progenitor

MDS:

myelodysplastic syndrome

sAML:

secondary acute myeloid leukemia

scRNA-seq:

single-cell RNA sequencing

TCGA:

The Cancer Genome Atlas

References

  1. Ganan-Gomez I, Yang H, Ma F, et al. Stem cell architecture drives myelodysplastic syndrome progression and predicts response to venetoclax-based therapy. Nat Med. 2022;28(3):557–67.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Oetjen KA, Lindblad KE, Goswami M, et al. Human bone marrow assessment by single-cell RNA sequencing, mass cytometry, and flow cytometry. JCI Insight. 2018;3(23):e124928.

    Article  PubMed Central  Google Scholar 

  3. Miles LA. Bowman RL, Merlinsky TR, et al Single-cell mutation analysis of clonal evolution in myeloid malignancies. Nature. 2020;587(7834):477–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Saha S, Murmu KC, Biswas M, et al. Transcriptomic analysis identifies RNA Binding proteins as putative regulators of Myelopoiesis and Leukemia. Front Oncol. 2019;9:692.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Signer RA, Qi L, Zhao Z, et al. The rate of protein synthesis in hematopoietic stem cells is limited partly by 4E-BPs. Genes Dev. 2016;30(15):1698–703.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank Christopher S Hourigan for providing single-cell RNA sequencing data of HDs.

Funding

This research was funded by the National Natural Science Foundation of China (No. 81970116, 81800123), the Zhao Yicheng (Youth Incubation) Fund of Tianjin Medical University General Hospital (No. ZYYFY2017026), Tianjin Key Medical Discipline(Specialty) Construction Project.

Author information

Authors and Affiliations

Authors

Contributions

YL conceived the study, designed and performed research, analyzed data and wrote the manuscript. HN and NS designed and performed research, performed bioinformatics and statistical analyses. WZ and LL supervised the study, analyzed and interpreted data. HW, RF and ZS conceived and supervised the study, interpreted the data and revised the manuscript. All authors discussed the results and approved the final version of the manuscript.

Corresponding authors

Correspondence to Huaquan Wang, Rong Fu or Zonghong Shao.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of Tianjin Medical University General Hospital (NO. IRB2021-WZ-051, April 2021). Informed consent was obtained from all subjects involved in the study.

Consent for publication

The authors affirm that human research participants provided informed consent for publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Supplementary Material 2

Supplementary Material 3

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.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, Y., Niu, H., Song, N. et al. Single-cell RNA sequencing identifies the properties of myelodysplastic syndrome stem cells. J Transl Med 20, 499 (2022). https://doi.org/10.1186/s12967-022-03709-9

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-022-03709-9