Transcriptional, Protein-Level and Functional Profiling of Human Fetal Liver (FL)-Derived Hematopoietic Stem Cells (HSCs) at Single Cell Resolution

Intro: The complex and tightly regulated process of human hematopoietic development culminates in the production of hematopoietic stem cells (HSCs), which subsequently acquire functional competence and undergo expansion in the fetal liver (FL). The establishment of a high-resolution molecular signature of FL HSCs provides insights into HSC biology with potential utility in the purification and expansion of engraftable HSCs ex vivo and the generation of HSCs from pluripotent cell sources. To profile HSCs at this developmental stage, we performed CITE-Seq, a technique that combines single cell RNA sequencing (scRNAseq) and cell surface marker interrogation using oligo-tagged antibodies to simultaneously map transcriptional and protein-level expression in individual cells. To connect expression profiles with functional engraftment, we have coupled this with transplantation assays in immunocompromised mice. Methods: In these studies, three populations of human FL cells were used: CD34- cells, CD34+ cells and CD34+ cells furtherenrichedby expression of GPI-80, a marker tightly linked to engraftment potential, to explicitly identify HSCs capable of long-term engraftment. These populations were stained with a panel of oligo-tagged antibodies, processed via the 10X Genomics platform, and sequenced (26,407 total cells). Results: Transplantation experiments using the same sorted fractions that were assayed by CITE-seq revealed superior engraftment potential of the GPI-80+ fraction, and thus enrichment for bona fide HSCs at the functional level. This functional signature coincided with enrichment for known HSC markers such as ITGA6 (CD49f), PROCR (EPCR), CD164, MLLT3, HLF, CLEC9A and HMGA2 at the transcriptional level. As such, by profiling >7000 GPI-80+ cells, we have achieved unprecedented resolution of the engraftable HSC compartment within the FL. Combined analysis of all captured FL fractions accurately recapitulated the hematopoietic landscape of the FL at this developmental stage, representing the expected hematopoietic lineages and cell types. To gain further insight into FL HSCs, we next focused on the CD34+ HSC/progenitor compartment where we tracked cluster dynamics upon functional HSC enrichment between the CD34+ bulk and GPI-80+ sample. We noted a prominent (4-fold) increase in a cluster marked by enrichment for genes including RGCC, LMNA, VIM, ID1 and ID3 as well as components of the AHR pathway, suggesting that this expression profile strongly correlat