Decoding lymphomyeloid divergence and immune hyporesponsiveness in G-CSF-primed human bone marrow by single-cell RNA-seq

Granulocyte colony-stimulating factor (G-CSF) has been widely used to mobilize bone marrow hematopoietic stem/progenitor cells for transplantation in the treatment of hematological malignancies for decades. Additionally, G-CSF is also accepted as an essential mediator in immune regulation, leading to reduced graft-versus-host disease following transplantation. Despite the important clinical roles of G-CSF, a comprehensive, unbiased, and high-resolution survey into the cellular and molecular ecosystem of the human G-CSF-primed bone marrow (G-BM) is lacking so far. Here, we employed single-cell RNA sequencing to profile hematopoietic cells in human bone marrow from two healthy donors before and after 5-day G-CSF administration. Through unbiased bioinformatics analysis, our data systematically showed the alterations in the transcriptional landscape of hematopoietic cells in G-BM, and revealed that G-CSF-induced myeloid-biased differentiation initiated from the stage of lymphoid-primed multipotent progenitors. We also illustrated the cellular and molecular basis of hyporesponsiveness of T cells and natural killer (NK) cells caused by G-CSF stimulation, including the potential direct mechanisms and indirect regulations mediated by ligand-receptor interactions. Taken together, our data extend the understanding of lymphomyeloid divergence and potential mechanisms involved in hyporesponsiveness of T and NK cells in human G-BM, which might provide basis for optimization of stem cell transplantation in hematological malignancy treatment.