Identification of phenotypically, functionally, and anatomically distinct stromal niche populations in human bone marrow based on single-cell RNA sequencing
Abstract
Hematopoiesis is regulated by the bone marrow (BM) stroma. However, cellular identities and functions of the different BM stromal elements in humans remain poorly defined. Based on single-cell RNA sequencing (scRNAseq), we systematically characterized the human non-hematopoietic BM stromal compartment and we investigated stromal cell regulation principles based on the RNA velocity analysis using scVelo and studied the interactions between the human BM stromal cells and hematopoietic cells based on ligand-receptor (LR) expression using CellPhoneDB. scRNAseq led to the identification of six transcriptionally and functionally distinct stromal cell populations. Stromal cell differentiation hierarchy was recapitulated based on RNA velocity analysis and in vitro proliferation capacities and differentiation potentials. Potential key factors that might govern the transition from stem and progenitor cells to fate-committed cells were identified. In situ localization analysis demonstrated that different stromal cells were localized in different niches in the bone marrow. In silico cell-cell communication analysis further predicted that different stromal cell types might regulate hematopoiesis through distinct mechanisms. These findings provide the basis for a comprehensive understanding of the cellular complexity of the human BM microenvironment and the intricate stroma-hematopoiesis crosstalk mechanisms, thus refining our current view on human hematopoietic niche organization.
Data availability
The scRNA-seq matrix data generated in this study have been deposited in the GEO database (GSE190965).
-
Transcriptomic profiling of human bone marrow non-hematopoietic cellsNCBI Gene Expression Omnibus, GSE190965.
Article and author information
Author details
Funding
Swedish Cancer Foundation (20-1163PjF 01H)
- Stefan Scheding
Swedish Childhood Cancer Foundation (PR2018-0078)
- Stefan Scheding
Swedish Childhood Cancer Foundation (PR2021-0065)
- Stefan Scheding
Swedish Bloodcancer Association (N.A.)
- Stefan Scheding
Swedish Research Council (N.A.)
- Stefan Scheding
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
Human subjects: Human bone marrow (BM) cells were collected at the Hematology Department, Skåne University Hospital Lund, Sweden, from consenting healthy donors. The use of human samples was approved by the Regional Ethics Review Board in Lund, Sweden.
Copyright
© 2023, Li et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,478
- views
-
- 703
- downloads
-
- 20
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Stem Cells and Regenerative Medicine
Haploinsufficiency for GATA6 is associated with congenital heart disease (CHD) with variable comorbidity of pancreatic or diaphragm defects, although the etiology of disease is not well understood. Here, we used cardiac directed differentiation from human embryonic stem cells (hESCs) as a platform to study GATA6 function during early cardiogenesis. GATA6 loss-of-function hESCs had a profound impairment in cardiac progenitor cell (CPC) specification and cardiomyocyte (CM) generation due to early defects during the mesendoderm and lateral mesoderm patterning stages. Profiling by RNA-seq and CUT&RUN identified genes of the WNT and BMP programs regulated by GATA6 during early mesoderm patterning. Furthermore, interactome analysis detected GATA6 binding with developmental transcription factors and chromatin remodelers, suggesting cooperative regulation of cardiac lineage gene accessibility. We show that modulating WNT and BMP inputs during the first 48 hr of cardiac differentiation is sufficient to partially rescue CPC and CM defects in GATA6 heterozygous and homozygous mutant hESCs. This study provides evidence of the regulatory functions for GATA6 directing human precardiac mesoderm patterning during the earliest stages of cardiogenesis to further our understanding of haploinsufficiency causing CHD and the co-occurrence of cardiac and other organ defects caused by human GATA6 mutations.
-
- Cell Biology
- Stem Cells and Regenerative Medicine
Conjugated linoleic acids (CLAs) can serve as a nutritional intervention to regulate quality, function, and fat infiltration in skeletal muscles, but the specific cytological mechanisms remain unknown. Here, we applied single-nucleus RNA-sequencing (snRNA-seq) to characterize the cytological mechanism of CLAs regulates fat infiltration in skeletal muscles based on pig models. We investigated the regulatory effects of CLAs on cell populations and molecular characteristics in pig muscles and found CLAs could promote the transformation of fast glycolytic myofibers into slow oxidative myofibers. We also observed three subpopulations including SCD+/DGAT2+, FABP5+/SIAH1+, and PDE4D+/PDE7B+ subclusters in adipocytes and CLAs could increase the percentage of SCD+/DGAT2+ adipocytes. RNA velocity analysis showed FABP5+/SIAH1+ and PDE4D+/PDE7B+ adipocytes could differentiate into SCD+/DGAT2+ adipocytes. We further verified the differentiated trajectory of mature adipocytes and identified PDE4D+/PDE7B+ adipocytes could differentiate into SCD+/DGAT2+ and FABP5+/SIAH1+ adipocytes by using high intramuscular fat (IMF) content Laiwu pig models. The cell-cell communication analysis identified the interaction network between adipocytes and other subclusters such as fibro/adipogenic progenitors (FAPs). Pseudotemporal trajectory analysis and RNA velocity analysis also showed FAPs could differentiate into PDE4D+/PDE7B+ preadipocytes and we discovered the differentiated trajectory of preadipocytes into mature adipocytes. Besides, we found CLAs could promote FAPs differentiate into SCD+/DGAT2+ adipocytes via inhibiting c-Jun N-terminal kinase (JNK) signaling pathway in vitro. This study provides a foundation for regulating fat infiltration in skeletal muscles by using nutritional strategies and provides potential opportunities to serve pig as an animal model to study human fat infiltrated diseases.