An immunophenotype-coupled transcriptomic atlas of human hematopoietic progenitors.

Analysis of the human hematopoietic progenitor compartment is being transformed by single-cell multimodal approaches. Cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) enables coupled surface protein and transcriptome profiling, thereby revealing genomic programs underlying progenitor states. To perform CITE-seq systematically on primary human bone marrow cells, we used titrations with 266 CITE-seq antibodies (antibody-derived tags) and machine learning to optimize a panel of 132 antibodies. Multimodal analysis resolved >80 stem, progenitor, immune, stromal and transitional cells defined by distinctive surface markers and transcriptomes. This dataset enables flow cytometry solutions for in silico-predicted cell states and identifies dozens of cell surface markers consistently detected across donors spanning race and sex. Finally, aligning annotations from this atlas, we nominate normal marrow equivalents for acute myeloid leukemia stem cell populations that differ in clinical response. This atlas serves as an advanced digital resource for hematopoietic progenitor analyses in human health and disease.

ThPOK is a critical multifaceted regulator of myeloid lineage development.

The transcription factor ThPOK (encoded by Zbtb7b) is well known for its role as a master regulator of CD4 lineage commitment in the thymus. Here, we report an unexpected and critical role of ThPOK as a multifaceted regulator of myeloid lineage commitment, differentiation and maturation. Using reporter and knockout mouse models combined with single-cell RNA-sequencing, progenitor transfer and colony assays, we show that ThPOK controls monocyte-dendritic cell versus granulocyte lineage production during homeostatic differentiation, and serves as a brake for neutrophil maturation in granulocyte lineage-specified cells through transcriptional regulation of lineage-specific transcription factors and RNA via altered messenger RNA splicing to reprogram intron retention.

Essential role of a ThPOK autoregulatory loop in the maintenance of mature CD4+ T cell identity and function.

The transcription factor ThPOK (encoded by the Zbtb7b gene) controls homeostasis and differentiation of mature helper T cells, while opposing their differentiation to CD4+ intraepithelial lymphocytes (IELs) in the intestinal mucosa. Thus CD4 IEL differentiation requires ThPOK transcriptional repression via reactivation of the ThPOK transcriptional silencer element (SilThPOK). In the present study, we describe a new autoregulatory loop whereby ThPOK binds to the SilThPOK to maintain its own long-term expression in CD4 T cells. Disruption of this loop in vivo prevents persistent ThPOK expression, leads to genome-wide changes in chromatin accessibility and derepresses the colonic regulatory T (Treg) cell gene expression signature. This promotes selective differentiation of naive CD4 T cells into GITRloPD-1loCD25lo (Triplelo) Treg cells and conversion to CD4+ IELs in the gut, thereby providing dominant protection from colitis. Hence, the ThPOK autoregulatory loop represents a key mechanism to physiologically control ThPOK expression and T cell differentiation in the gut, with potential therapeutic relevance.

Resilient anatomy and local plasticity of naive and stress haematopoiesis.

The bone marrow adjusts blood cell production to meet physiological demands in response to insults. The spatial organization of normal and stress responses are unknown owing to the lack of methods to visualize most steps of blood production. Here we develop strategies to image multipotent haematopoiesis, erythropoiesis and lymphopoiesis in mice. We combine these with imaging of myelopoiesis1 to define the anatomy of normal and stress haematopoiesis. In the steady state, across the skeleton, single stem cells and multipotent progenitors distribute through the marrow enriched near megakaryocytes. Lineage-committed progenitors are recruited to blood vessels, where they contribute to lineage-specific microanatomical structures composed of progenitors and immature cells, which function as the production sites for each major blood lineage. This overall anatomy is resilient to insults, as it was maintained after haemorrhage, systemic bacterial infection and granulocyte colony-stimulating factor (G-CSF) treatment, and during ageing. Production sites enable haematopoietic plasticity as they differentially and selectively modulate their numbers and output in response to insults. We found that stress responses are variable across the skeleton: the tibia and the sternum respond in opposite ways to G-CSF, and the skull does not increase erythropoiesis after haemorrhage. Our studies enable in situ analyses of haematopoiesis, define the anatomy of normal and stress responses, identify discrete microanatomical production sites that confer plasticity to haematopoiesis, and uncover unprecedented heterogeneity of stress responses across the skeleton.

Combinatorial Single-Cell Analyses of Granulocyte-Monocyte Progenitor Heterogeneity Reveals an Early Uni-potent Neutrophil Progenitor.

Granulocyte-monocyte progenitors (GMPs) have been previously defined for their potential to generate various myeloid progenies such as neutrophils and monocytes. Although studies have proposed lineage heterogeneity within GMPs, it is unclear if committed progenitors already exist among these progenitors and how they may behave differently during inflammation. By combining single-cell transcriptomic and proteomic analyses, we identified the early committed progenitor within the GMPs responsible for the strict production of neutrophils, which we designate as proNeu1. Our dissection of the GMP hierarchy led us to further identify a previously unknown intermediate proNeu2 population. Similar populations could be detected in human samples. proNeu1s, but not proNeu2s, selectively expanded during the early phase of sepsis at the expense of monocytes. Collectively, our findings help shape the neutrophil maturation trajectory roadmap and challenge the current definition of GMPs.

In situ mapping identifies distinct vascular niches for myelopoiesis.

In contrast to nearly all other tissues, the anatomy of cell differentiation in the bone marrow remains unknown. This is owing to a lack of strategies for examining myelopoiesis-the differentiation of myeloid progenitors into a large variety of innate immune cells-in situ in the bone marrow. Such strategies are required to understand differentiation and lineage-commitment decisions, and to define how spatial organizing cues inform tissue function. Here we develop approaches for imaging myelopoiesis in mice, and generate atlases showing the differentiation of granulocytes, monocytes and dendritic cells. The generation of granulocytes and dendritic cells-monocytes localizes to different blood-vessel structures known as sinusoids, and displays lineage-specific spatial and clonal architectures. Acute systemic infection with Listeria monocytogenes induces lineage-specific progenitor clusters to undergo increased self-renewal of progenitors, but the different lineages remain spatially separated. Monocyte-dendritic cell progenitors (MDPs) map with nonclassical monocytes and conventional dendritic cells; these localize to a subset of blood vessels expressing a major regulator of myelopoiesis, colony-stimulating factor 1 (CSF1, also known as M-CSF)1. Specific deletion of Csf1 in endothelium disrupts the architecture around MDPs and their localization to sinusoids. Subsequently, there are fewer MDPs and their ability to differentiate is reduced, leading to a loss of nonclassical monocytes and dendritic cells during both homeostasis and infection. These data indicate that local cues produced by distinct blood vessels are responsible for the spatial organization of definitive blood cell differentiation.

Granulocyte-Monocyte Progenitors and Monocyte-Dendritic Cell Progenitors Independently Produce Functionally Distinct Monocytes.

Granulocyte-monocyte progenitors (GMPs) and monocyte-dendritic cell progenitors (MDPs) produce monocytes during homeostasis and in response to increased demand during infection. Both progenitor populations are thought to derive from common myeloid progenitors (CMPs), and a hierarchical relationship (CMP-GMP-MDP-monocyte) is presumed to underlie monocyte differentiation. Here, however, we demonstrate that mouse MDPs arose from CMPs independently of GMPs, and that GMPs and MDPs produced monocytes via similar but distinct monocyte-committed progenitors. GMPs and MDPs yielded classical (Ly6Chi) monocytes with gene expression signatures that were defined by their origins and impacted their function. GMPs produced a subset of "neutrophil-like" monocytes, whereas MDPs gave rise to a subset of monocytes that yielded monocyte-derived dendritic cells. GMPs and MDPs were also independently mobilized to produce specific combinations of myeloid cell types following the injection of microbial components. Thus, the balance of GMP and MDP differentiation shapes the myeloid cell repertoire during homeostasis and following infection.

Mouse models of neutropenia reveal progenitor-stage-specific defects.

Advances in genetics and sequencing have identified a plethora of disease-associated and disease-causing genetic alterations. To determine causality between genetics and disease, accurate models for molecular dissection are required; however, the rapid expansion of transcriptional populations identified through single-cell analyses presents a major challenge for accurate comparisons between mutant and wild-type cells. Here we generate mouse models of human severe congenital neutropenia (SCN) using patient-derived mutations in the GFI1 transcription factor. To determine the effects of SCN mutations, we generated single-cell references for granulopoietic genomic states with linked epitopes1, aligned mutant cells to their wild-type equivalents and identified differentially expressed genes and epigenetic loci. We find that GFI1-target genes are altered sequentially, as cells go through successive states of differentiation. These insights facilitated the genetic rescue of granulocytic specification but not post-commitment defects in innate immune effector function, and underscore the importance of evaluating the effects of mutations and therapy within each relevant cell state.

Author Correction: Single-cell analysis of mixed-lineage states leading to a binary cell fate choice.

In this Letter, the first name of author Virendra K. Chaudhri was incorrectly spelled 'Viren'; author Meenakshi Venkatasubramanian should also be associated with 'Department of Electrical Engineering and Computer Science, University of Cincinnati, Cincinnati, Ohio 45221, USA'; authors Bruce J. Aronow, Nathan Salomonis, Harinder Singh and H. Leighton Grimes should also be associated with 'Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, Ohio 45229, USA'. The Letter has not been corrected online.

pyInfinityFlow: optimized imputation and analysis of high-dimensional flow cytometry data for millions of cells.

MOTIVATION: While conventional flow cytometry is limited to dozens of markers, new experimental and computational strategies, such as Infinity Flow, allow for the generation and imputation of hundreds of cell surface protein markers in millions of cells. Here, we describe an end-to-end analysis workflow for Infinity Flow data in Python. RESULTS: pyInfinityFlow enables the efficient analysis of millions of cells, without down-sampling, through direct integration with well-established Python packages for single-cell genomics analysis. pyInfinityFlow accurately identifies both common and extremely rare cell populations which are challenging to define from single-cell genomics studies alone. We demonstrate that this workflow can nominate novel markers to design new flow cytometry gating strategies for predicted cell populations. pyInfinityFlow can be extended to diverse cell discovery analyses with flexibility to adapt to diverse Infinity Flow experimental designs. AVAILABILITY AND IMPLEMENTATION: pyInfinityFlow is freely available in GitHub (https://github.com/KyleFerchen/pyInfinityFlow) and on PyPI (https://pypi.org/project/pyInfinityFlow/). Package documentation with tutorials on a test dataset is available by Read the Docs (pyinfinityflow.readthedocs.io). The scripts and data for reproducing the results are available at https://github.com/KyleFerchen/pyInfinityFlow/tree/main/analysis_scripts, along with the raw flow cytometry input data.

Erythroblastic islands foster granulopoiesis in parallel to terminal erythropoiesis.

The erythroblastic island (EBI), composed of a central macrophage surrounded by maturing erythroblasts, is the erythroid precursor niche. Despite numerous studies, its precise composition is still unclear. Using multispectral imaging flow cytometry, in vitro island reconstitution, and single-cell RNA sequencing of adult mouse bone marrow (BM) EBI-component cells enriched by gradient sedimentation, we present evidence that the CD11b+ cells present in the EBIs are neutrophil precursors specifically associated with BM EBI macrophages, indicating that erythro-(myelo)-blastic islands are a site for terminal granulopoiesis and erythropoiesis. We further demonstrate that the balance between these dominant and terminal differentiation programs is dynamically regulated within this BM niche by pathophysiological states that favor granulopoiesis during anemia of inflammation and favor erythropoiesis after erythropoietin stimulation. Finally, by molecular profiling, we reveal the heterogeneity of EBI macrophages by cellular indexing of transcriptome and epitope sequencing of mouse BM EBIs at baseline and after erythropoietin stimulation in vivo and provide a searchable online viewer of these data characterizing the macrophage subsets serving as hematopoietic niches. Taken together, our findings demonstrate that EBIs serve a dual role as niches for terminal erythropoiesis and granulopoiesis and the central macrophages adapt to optimize production of red blood cells or neutrophils.

HDAC11 deficiency disrupts oncogene-induced hematopoiesis in myeloproliferative neoplasms.

Protein acetylation is an important contributor to cancer initiation. Histone deacetylase 6 (HDAC6) controls JAK2 translation and protein stability and has been implicated in JAK2-driven diseases best exemplified by myeloproliferative neoplasms (MPNs). By using novel classes of highly selective HDAC inhibitors and genetically deficient mouse models, we discovered that HDAC11 rather than HDAC6 is necessary for the proliferation and survival of oncogenic JAK2-driven MPN cells and patient samples. Notably, HDAC11 is variably expressed in primitive stem cells and is expressed largely upon lineage commitment. Although Hdac11is dispensable for normal homeostatic hematopoietic stem and progenitor cell differentiation based on chimeric bone marrow reconstitution, Hdac11 deficiency significantly reduced the abnormal megakaryocyte population, improved splenic architecture, reduced fibrosis, and increased survival in the MPLW515L-MPN mouse model during primary and secondary transplantation. Therefore, inhibitors of HDAC11 are an attractive therapy for treating patients with MPN. Although JAK2 inhibitor therapy provides substantial clinical benefit in MPN patients, the identification of alternative therapeutic targets is needed to reverse MPN pathogenesis and control malignant hematopoiesis. This study establishes HDAC11 as a unique type of target molecule that has therapeutic potential in MPN.

Single-cell analysis of mixed-lineage states leading to a binary cell fate choice.

Delineating hierarchical cellular states, including rare intermediates and the networks of regulatory genes that orchestrate cell-type specification, are continuing challenges for developmental biology. Single-cell RNA sequencing is greatly accelerating such research, given its power to provide comprehensive descriptions of genomic states and their presumptive regulators. Haematopoietic multipotential progenitor cells, as well as bipotential intermediates, manifest mixed-lineage patterns of gene expression at a single-cell level. Such mixed-lineage states may reflect the molecular priming of different developmental potentials by co-expressed alternative-lineage determinants, namely transcription factors. Although a bistable gene regulatory network has been proposed to regulate the specification of either neutrophils or macrophages, the nature of the transition states manifested in vivo, and the underlying dynamics of the cell-fate determinants, have remained elusive. Here we use single-cell RNA sequencing coupled with a new analytic tool, iterative clustering and guide-gene selection, and clonogenic assays to delineate hierarchical genomic and regulatory states that culminate in neutrophil or macrophage specification in mice. We show that this analysis captured prevalent mixed-lineage intermediates that manifested concurrent expression of haematopoietic stem cell/progenitor and myeloid progenitor cell genes. It also revealed rare metastable intermediates that had collapsed the haematopoietic stem cell/progenitor gene expression programme, instead expressing low levels of the myeloid determinants, Irf8 and Gfi1 (refs 9, 10, 11, 12, 13). Genetic perturbations and chromatin immunoprecipitation followed by sequencing revealed Irf8 and Gfi1 as key components of counteracting myeloid-gene-regulatory networks. Combined loss of these two determinants 'trapped' the metastable intermediate. We propose that mixed-lineage states are obligatory during cell-fate specification, manifest differing frequencies because of their dynamic instability and are dictated by counteracting gene-regulatory networks.

cellHarmony: cell-level matching and holistic comparison of single-cell transcriptomes.

To understand the molecular pathogenesis of human disease, precision analyses to define alterations within and between disease-associated cell populations are desperately needed. Single-cell genomics represents an ideal platform to enable the identification and comparison of normal and diseased transcriptional cell populations. We created cellHarmony, an integrated solution for the unsupervised analysis, classification, and comparison of cell types from diverse single-cell RNA-Seq datasets. cellHarmony efficiently and accurately matches single-cell transcriptomes using a community-clustering and alignment strategy to compute differences in cell-type specific gene expression over potentially dozens of cell populations. Such transcriptional differences are used to automatically identify distinct and shared gene programs among cell-types and identify impacted pathways and transcriptional regulatory networks to understand the impact of perturbations at a systems level. cellHarmony is implemented as a python package and as an integrated workflow within the software AltAnalyze. We demonstrate that cellHarmony has improved or equivalent performance to alternative label projection methods, is able to identify the likely cellular origins of malignant states, stratify patients into clinical disease subtypes from identified gene programs, resolve discrete disease networks impacting specific cell-types, and illuminate therapeutic mechanisms. Thus, this approach holds tremendous promise in revealing the molecular and cellular origins of complex disease.

SKI controls MDS-associated chronic TGF-ß signaling, aberrant splicing, and stem cell fitness.

The transforming growth factor beta (TGF-ß) signaling pathway controls hematopoietic stem cell (HSC) behavior in the marrow niche; however, TGF-ß signaling becomes chronic in early-stage myelodysplastic syndrome (MDS). Although TGF-ß signaling normally induces negative feedback, in early-stage MDS, high levels of microRNA-21 (miR-21) contribute to chronic TGF-ß signaling. We found that a TGF-ß signal-correlated gene signature is sufficient to identify an MDS patient population with abnormal RNA splicing (eg, CSF3R) independent of splicing factor mutations and coincident with low HNRNPK activity. Levels of SKI messenger RNA (mRNA) encoding a TGF-ß antagonist are sufficient to identify these patients. However, MDS patients with high SKI mRNA and chronic TGF-ß signaling lack SKI protein because of miR-21 activity. To determine the impact of SKI loss, we examined murine Ski -/- HSC function. First, competitive HSC transplants revealed a profound defect in stem cell fitness (competitive disadvantage) but not specification, homing, or multilineage production. Aged recipients of Ski -/- HSCs exhibited mild phenotypes similar to phenotypes in those with macrocytic anemia. Second, blastocyst complementation revealed a dramatic block in Ski -/- hematopoiesis in the absence of transplantation. Similar to SKI-high MDS patient samples, Ski -/- HSCs strikingly upregulated TGF-ß signaling and deregulated expression of spliceosome genes (including Hnrnpk). Moreover, novel single-cell splicing analyses demonstrated that Ski -/- HSCs and high levels of SKI expression in MDS patient samples share abnormal alternative splicing of common genes (including those that encode splicing factors). We conclude that miR-21-mediated loss of SKI activates TGF-ß signaling and alternative splicing to impair the competitive advantage of normal HSCs (fitness), which could contribute to selection of early-stage MDS-genic clones.

Clonal hematopoiesis of indeterminate potential and its impact on patient trajectories after stem cell transplantation.

Clonal hematopoiesis of indeterminate potential (CHIP) is a recently identified process where older patients accumulate distinct subclones defined by recurring somatic mutations in hematopoietic stem cells. CHIP's implications for stem cell transplantation have been harder to identify due to the high degree of mutational heterogeneity that is present within the genetically distinct subclones. In order to gain a better understanding of CHIP and the impact of clonal dynamics on transplantation outcomes, we created a mathematical model of clonal competition dynamics. Our analyses highlight the importance of understanding competition intensity between healthy and mutant clones. Importantly, we highlight the risk that CHIP poses in leading to dominance of precancerous mutant clones and the risk of donor derived leukemia. Furthermore, we estimate the degree of competition intensity and bone marrow niche decline in mice during aging by using our modeling framework. Together, our work highlights the importance of better characterizing the ecological and clonal composition in hematopoietic donor populations at the time of stem cell transplantation.

The miR-23a~27a~24-2 microRNA cluster buffers transcription and signaling pathways during hematopoiesis.

MicroRNA cluster mirn23a has previously been shown to promote myeloid development at the expense of lymphoid development in overexpression and knockout mouse models. This polarization is observed early in hematopoietic development, with an increase in common lymphoid progenitors (CLPs) and a decrease in all myeloid progenitor subsets in adult bone marrow. The pool size of multipotential progenitors (MPPs) is unchanged; however, in this report we observe by flow cytometry that polarized subsets of MPPs are changed in the absence of mirn23a. Additionally, in vitro culture of MPPs and sorted MPP transplants showed that these cells have decreased myeloid and increased lymphoid potential in vitro and in vivo. We investigated the mechanism by which mirn23a regulates hematopoietic differentiation and observed that mirn23a promotes myeloid development of hematopoietic progenitors through regulation of hematopoietic transcription factors and signaling pathways. Early transcription factors that direct the commitment of MPPs to CLPs (Ikzf1, Runx1, Satb1, Bach1 and Bach2) are increased in the absence of mirn23a miRNAs as well as factors that commit the CLP to the B cell lineage (FoxO1, Ebf1, and Pax5). Mirn23a appears to buffer transcription factor levels so that they do not stochastically reach a threshold level to direct differentiation. Intriguingly, mirn23a also inversely regulates the PI3 kinase (PI3K)/Akt and BMP/Smad signaling pathways. Pharmacological inhibitor studies, coupled with dominant active/dominant negative biochemical experiments, show that both signaling pathways are critical to mirn23a's regulation of hematopoietic differentiation. Lastly, consistent with mirn23a being a physiological inhibitor of B cell development, we observed that the essential B cell transcription factor EBF1 represses expression of mirn23a. In summary, our data demonstrates that mirn23a regulates a complex array of transcription and signaling pathways to modulate adult hematopoiesis.

The cell polarity determinant CDC42 controls division symmetry to block leukemia cell differentiation.

As a central regulator of cell polarity, the activity of CDC42 GTPase is tightly controlled in maintaining normal hematopoietic stem and progenitor cell (HSC/P) functions. We found that transformation of HSC/P to acute myeloid leukemia (AML) is associated with increased CDC42 expression and activity in leukemia cells. In a mouse model of AML, the loss of Cdc42 abrogates MLL-AF9-induced AML development. Furthermore, genetic ablation of CDC42 in both murine and human MLL-AF9 (MA9) cells decreased survival and induced differentiation of the clonogenic leukemia-initiating cells. We show that MLL-AF9 leukemia cells maintain cell polarity in the context of elevated Cdc42-guanosine triphosphate activity, similar to nonmalignant, young HSC/Ps. The loss of Cdc42 resulted in a shift to depolarized AML cells that is associated with a decrease in the frequency of symmetric and asymmetric cell divisions producing daughter cells capable of self-renewal. Importantly, we demonstrate that inducible CDC42 suppression in primary human AML cells blocks leukemia progression in a xenograft model. Thus, CDC42 loss suppresses AML cell polarity and division asymmetry, and CDC42 constitutes a useful target to alter leukemia-initiating cell fate for differentiation therapy.

Temporal Expression of Bim Limits the Development of Agonist-Selected Thymocytes and Skews Their TCRß Repertoire.

CD8αα TCRαß+ intestinal intraepithelial lymphocytes play a critical role in promoting intestinal homeostasis, although mechanisms controlling their development and peripheral homeostasis remain unclear. In this study, we examined the spatiotemporal role of Bim in the thymic selection of CD8αα precursors and the fate of these cells in the periphery. We found that T cell-specific expression of Bim during early/cortical, but not late/medullary, thymic development controls the agonist selection of CD8αα precursors and limits their private TCRß repertoire. During this process, agonist-selected double-positive cells lose CD4/8 coreceptor expression and masquerade as double-negative (DN) TCRαßhi thymocytes. Although these DN thymocytes fail to re-express coreceptors after OP9-DL1 culture, they eventually mature and accumulate in the spleen where TCR and IL-15/STAT5 signaling promotes their conversion to CD8αα cells and their expression of gut-homing receptors. Adoptive transfer of splenic DN cells gives rise to CD8αα cells in the gut, establishing their precursor relationship in vivo. Interestingly, Bim does not restrict the IL-15-driven maturation of CD8αα cells that is critical for intestinal homeostasis. Thus, we found a temporal and tissue-specific role for Bim in limiting thymic agonist selection of CD8αα precursors and their TCRß repertoire, but not in the maintenance of CD8αα intraepithelial lymphocytes in the intestine.