Abstract:Understanding the molecular signatures of human hematopoietic stem/progenitor cells (HSCPs) and gene network changes leading to differentiation induction are of utmost importance to achieve therapeutic hematopoietic stem cell (HSC) expansion. The identification of surface markers has improved the prospective isolation of human bone marrow HSCs over the last decades, but the isolated population still remains very heterogeneous. Single cell sequencing approaches have the potential to unbiased assort stem cells along their differentiation stages and trajectories, and to decipher molecular markers in differentiation hierarchies. However, whole transcriptome single cell sequencing often suffers from inadequate detection of low-expressed genes with increased drop-outs. To address this limitation and to provide sensitive mRNA quantification in human HSPCs, we developed a selected targeted gene panel (600 genes) in combination with 50 protein markers to cover most relevant genes/proteins of early HSPC differentiation and cell cycle status. We employed single cell CITE-Seq (cellular indexing of transcriptomes and epitopes by sequencing) based multi-omic profiling to >60.000 high-quality FACS-enriched human bone marrow HSPCs (CD34 + and CD34 +CD38 - cells) from 15 healthy donors representing young, mid and old age groups. Here, we demonstrate the most comprehensive molecular map of early human HSC differentiation and provide evidence of early branching points into lineage trajectories across different age groups using Slingshot pseudotime analysis. Our data demonstrate that the HSC branched into two main lineages, with the first branching point giving rise to megakaryocytic and erythroid progenitors (MKP/ERP), followed by lymphoid and myeloid progenitors (LMP/MDP). As expected, the most immature HSCs were defined by high expression of CD34, CD90, CRHBP, HOPX, HLF, MLLT3. However, we observed a number of hitherto undescribed HSC markers. The MKP/ERP lineage was defined by GATA1, TFRC, CD36, MPL and VWF expression. Furthermore, we found CD274, the Programmed Cell Death receptor ligand 1, to be selectively expressed in the MKP cluster. Re-clustering and trajectory analysis of the most immature CD34 +CD38 -/low cells confirmed early lineage commitment, suggesting that an MKP/ERP development does not include myeloid or lymphoid progenitors. A similar pattern was observed for all age groups. Interestingly, we observed increased expression of novel marker genes ( HLA-E, DLK1, ADGRG6, ADGRG1, MMRN1) and proteins (CD137 and CD273) in the most immature HSCs. Besides a gradual upregulation of CD38 along differentiation, differentiated cells showed lower CD45 levels. Furthermore, MKP/ERP committed cells showed a downregulation of SPINK2. TradeSeq analysis defined the continuous single cell expression changes in genes and markers upon early differentiation and lineage commitment, allowing the detection of decision-making gene networks. Our data adds further evidence to the revised model of hematopoiesis that challenges the dogma of a common myeloid progenitor giving rise to erythroid, megakaryocytic and myeloid cells. Furthermore, we see a structured hierarchical organization with a sequential loss of lineage potential and not an early commitment of uni-lineage restricted precursors. We report on the upregulation of ADGRG6, ADGRG1, DLK1, HLA-E, MMRN1 CD137 and CD273 in the most immature HSCs. Molecular and functional characterization of CD273 hi and low HSCs are presented in an adjacent report by our group at this meeting (Schmachtel et al.). In conclusion, we provide a valuable resource for the continuous regulation of genes and surface proteins along early differentiation of very immature human HSCs.

Single-Cell Trajectory Detection Uncovers Progression and Regulatory Coordination in Human B Cell Development

Tviblindi algorithm identifies branching developmental trajectories of human B cell development

Beyond pseudotime: Following T-cell maturation in single-cell RNAseq time series

Tviblindi algorithm identifies branching developmental trajectories of human B-cell development and describes abnormalities in RAG-1 and WAS patients

A Single-Cell Hematopoietic Landscape Resolves 8 Lineage Trajectories and Defects in Kit Mutant Mice.

Concepts and Limitations for Learning Developmental Trajectories from Single Cell Genomics.

Harnessing spatial transcriptomics to study fate decisions and lineage specification in T cell development

Decoding the Regulatory Network of Early Blood Development from Single-Cell Gene Expression Measurements.

Pseudocell Tracer-A method for inferring dynamic trajectories using scRNAseq and its application to B cells undergoing immunoglobulin class switch recombination

Mapping the human DC lineage through the integration of high-dimensional techniques

Deconstructing Complexity: A Computational Topology Approach to Trajectory Inference in the Human Thymus with

Deciphering cell states and genealogies of human hematopoiesis

Single-cell RNA-seq analysis revealed the stemness of a specific cluster of B cells in acute lymphoblastic leukemia progression

Single-cell analysis reveals multi-faceted features of B cell development together with age-associated B cell subpopulations

High-resolution mapping of cell-cycle dynamics during T-cell development and regeneration in vivo

A landscape of haematopoiesis informed by cellular barcoding and agent-based modelling

High-resolution mapping of cell-cycle dynamics during steady-state T-cell development and regeneration in vivo

A single-cell multi-omic and spatial atlas of B cell lymphomas reveals differentiation drives intratumor heterogeneity

Advances and challenges in investigating B-cells via single-cell transcriptomics

Identification and Isolation of Hematopoietic Stem Cells and Progenitors with Discrete Developmental Gene Expression Programs

Early Molecular Steps in Human Bone Marrow Stem Cell Differentiation Trajectories