Computing cell state discriminates the aberrant hematopoiesis and activated microenvironment in Myelodysplastic syndrome (MDS) through a single cell genomic study

Abstract

Abstract Background Myelodysplastic syndrome (MDS) is a complicated hematopoietic malignancy characterized by bone marrow (BM) dysplasia with symptoms like anemia, neutropenia, or thrombocytopenia. MDS exhibits considerable heterogeneity in prognosis, with approximately 30% of patients progressing to acute myeloid leukemia (AML). Single cell RNA-sequencing (scRNA-seq) is a new and powerful technique to profile disease landscapes. However, the current available scRNA-seq datasets for MDS are only focused on CD34+ hematopoietic progenitor cells. We argue that using entire BM cell for MDS studies probably will be more informative for understanding the pathophysiology of MDS. Methods Five MDS patients and four healthy donors were enrolled in the study. Unsorted cells from BM aspiration were collected for scRNA-seq analysis to profile overall alteration in hematopoiesis. Results Standard scRNA-seq analysis of unsorted BM cells successfully profiles deficient hematopoiesis in all five MDS patients, with three classified as high-risk and two as low-risk. While no significant increase in mutation burden was observed, high-risk MDS patients exhibited T-cell activation and abnormal myelogenesis at the stages between hematopoietic stem and progenitor cells (HSPC) and granulocyte–macrophage progenitors (GMP). Transcriptional factor analysis on the aberrant myelogenesis suggests that the epigenetic regulator chromatin structural protein-encoding gene HMGA1 is highly activated in the high-risk MDS group and moderately activated in the low-risk MDS group. Perturbation of HMGA1 by CellOracle simulated deficient hematopoiesis in mouse Lineage-negative (Lin-) BM cells. Projecting MDS and AML cells on a BM cell reference by our newly developed MarcoPolo pipeline intuitively visualizes a connection for myeloid leukemia development and abnormalities of hematopoietic hierarchy, indicating that it is technically feasible to integrate all diseased bone marrow cells on a common reference map even when the size of the cohort reaches to 1,000 patients or more. Conclusion Through scRNA-seq analysis on unsorted cells from BM aspiration samples of MDS patients, this study systematically profiled the development abnormalities in hematopoiesis, heterogeneity of risk, and T-cell microenvironment at the single cell level.

Keywords

• Myelodysplastic syndrome
• Single cell omics
• Computational hematopoiesis
• Bone marrow cell reference
• HMGA1