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.

PUF60-related developmental disorder: A case series and phenotypic analysis of 10 additional patients with monoallelic PUF60 variants.

PUF60-related developmental disorder (also referred to as Verheij syndrome), resulting from haploinsufficiency of PUF60, is associated with multiple congenital anomalies affecting a wide range of body systems. These anomalies include ophthalmic coloboma, and congenital anomalies of the heart, kidney, and musculoskeletal system. Behavioral and intellectual difficulties are also observed. While less common than other features associated with PUF60-related developmental disorder, for instance hearing impairment and short stature, identification of specific anomalies such as ophthalmic coloboma can aid with diagnostic identification given the limited spectrum of genes linked with this feature. We describe 10 patients with PUF60 gene variants, bringing the total number reported in the literature, to varying levels of details, to 56 patients. Patients were recruited both via locally based exome sequencing from international sites and from the DDD study in the United Kingdom. Eight of the variants reported were novel PUF60 variants. The addition of a further patient with a reported c449-457del variant to the existing literature highlights this as a recurrent variant. One variant was inherited from an affected parent. This is the first example in the literature of an inherited variant resulting in PUF60-related developmental disorder. Two patients (20%) were reported to have a renal anomaly consistent with 22% of cases in previously reported literature. Two patients received specialist endocrine treatment. More commonly observed were clinical features such as: cardiac anomalies (40%), ocular abnormalities (70%), intellectual disability (60%), and skeletal abnormalities (80%). Facial features did not demonstrate a recognizable gestalt. Of note, but remaining of unclear causality, we describe a single pediatric patient with pineoblastoma. We recommend that stature and pubertal progress should be monitored in PUF60-related developmental disorder with a low threshold for endocrine investigations as hormone therapy may be indicated. Our study reports an inherited case with PUF60-related developmental disorder which has important genetic counseling implications for families.

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.

A primer on single-cell genomics in myeloid biology.

PURPOSE OF REVIEW: Understanding the fast-moving field of single-cell technologies, as applied to myeloid biology, requires an appreciation of basic molecular, informatics, and biological concepts. Here, we highlight both key and recent articles to illustrate basic concepts for those new to molecular single-cell analyses in myeloid hematology. RECENT FINDINGS: Recent studies apply single-cell omics to discover novel cell populations, construct relationships between cell populations, reconfigure the organization of hematopoiesis, and study hematopoietic lineage tree and fate choices. Accompanying development of technologies, new informatic tools have emerged, providing exciting new insights. SUMMARY: Hematopoietic stem and progenitor cells are regulated by complex intrinsic and extrinsic factors to produce blood cell types. In this review, we discuss recent advances in single-cell omics to profile these cells, methods to infer cell type identify, and trajectories from molecular omics data to ultimately derive new insights into hematopoietic stem and progenitor cell biology. We further discuss future applications of these technologies to understand hematopoietic cell interactions, function, and development. The goal is to offer a comprehensive overview of current single-cell technologies and their impact on our understanding of myeloid cell development for those new to single-cell analyses.

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.

Volunteers' experiences of providing telephone-based breast-feeding peer support in the RUBY randomised controlled trial.

OBJECTIVE: The Ringing Up About Breastfeeding earlY (RUBY) randomised controlled trial (RCT) found that a telephone-based peer volunteer support intervention increased breast-feeding duration in a setting with high breast-feeding initiation. This sub-study of the RUBY RCT describes the motivation, preparation and experiences of volunteers who provided the peer support intervention. DESIGN: An online survey was completed by 154 (67 %) volunteers after ceasing volunteering. SETTING: Volunteers provided peer support to primiparous women (n 574) who birthed at one of three public hospitals in Melbourne, Australia, between February 2013 and December 2015. PARTICIPANTS: Volunteers (n 230) had themselves breastfed for at least 6 months and received 4 h of training for the role. RESULTS: The median number of mothers supported was two (range 1-11), and two-thirds of respondents supported at least one mother for 6 months. Volunteers were motivated by a strong desire to support new mothers to establish and continue breast-feeding. Most (93 %) considered the training session adequate. The majority (60 %) reported following the call schedule 'most of the time', but many commented that 'it depends on the mother'. Overall, 84 % of volunteers were satisfied with the role and reported that the experience was enjoyable (85 %) and worthwhile (90 %). Volunteers agreed that telephone support for breast-feeding was valued by women (88 %) and that the programme would be effective in helping women to breastfeed (93 %). CONCLUSIONS: These findings are important for those developing similar peer support programmes in which recruiting volunteers and developing training requirements are an integral and recurrent part of volunteer management.