SummaryHematopoietic stem cells give rise to all blood cells in a differentiation process that involves widespread epigenome remodeling. Here we present genome-wide reference maps of the associated DNA methylation dynamics. We used a meta-epigenomic approach that combines DNA methylation profiles across many small pools of cells and performed single-cell methylome sequencing to assess cell-to-cell heterogeneity. The resulting dataset identified characteristic differences between HSCs derived from fetal liver, cord blood, bone marrow, and peripheral blood. We also observed lineage-specific DNA methylation between myeloid and lymphoid progenitors, characterized immature multi-lymphoid progenitors, and detected progressive DNA methylation differences in maturing megakaryocytes. We linked these patterns to gene expression, histone modifications, and chromatin accessibility, and we used machine learning to derive a model of human hematopoietic differentiation directly from DNA methylation data. Our results contribute to a better understanding of human hematopoietic stem cell differentiation and provide a framework for studying blood-linked diseases.

Capturing where and how multipotency is lost is crucial to understand how blood formation is controlled. Blood lineage specification is currently thought to occur downstream of multipotent haematopoietic stem cells (HSC). Here we show that, in human, the first lineage restriction events occur within the CD19−CD34+CD38−CD45RA−CD49f+CD90+ (49f+) HSC compartment to generate myelo-lymphoid committed cells with no erythroid differentiation capacity. At single-cell resolution, we observe a continuous but polarised organisation of the 49f+ compartment, where transcriptional programmes and lineage potential progressively change along a gradient of opposing cell surface expression of CLEC9A and CD34. CLEC9AhiCD34lo cells contain long-term repopulating multipotent HSCs with slow quiescence exit kinetics, whereas CLEC9AloCD34hi cells are restricted to myelo-lymphoid differentiation and display infrequent but durable repopulation capacity. We thus propose that human HSCs gradually transition to a discrete lymphoid-primed state, distinct from lymphoid-primed multipotent progenitors, representing the earliest entry point into lymphoid commitment.