Search

SummaryHematopoietic stem cells (HSCs) develop within a short time window from the hemogenic endothelium in the aorta- gonads-and mesonephros (AGM) region during embryonic development. The first HSCs reside within Intra-aortic hematopoietic clusters (IAHC) along with hematopoietic progenitors (HPC). The signalling mechanisms that divert HSCs from HPCs are unknown. Notch signaling is essential for arterial specification, IAHC formation and HSC activity, but current studies on how Notch drives these different fates are inconsistent. To determine the role of Notch in the specification of hemogenic endothelium, HSC and/or HPCs, we extensively analyzed Notch dynamics in the period of HSC generation. We defined the expression pattern of Notch signalling molecules at the gene and protein level and established a molecular mechanism that reconcile previous studies demonstrating the loss of HSC activity in NOTCH1, JAG1 and RBPJ null mutants, the enhanced HSC generation by blocking specific Notch activities or the abrogation of emerging HSCs by high Notch activation. We now demonstrate that Notch activity is highest in a subset of Gfi1+ hemogenic endothelial cells and is gradually lost with HSC maturation. We uncover that the HSC phenotype is maintained through loss of Notch activity due to increasing levels of NOTCH1 and JAG1 interactions on the surface of the same cell (cis) that renders the NOTCH1 receptor from being activated. Forcing activation of the NOTCH1 receptor in IAHC cells activates a hematopoietic differentiation program and supports acis-inhibitory function for JAG1 and NOTCH1. Furthermore, we demonstrate that thiscis-inhibitory interaction is enabled by RADICAL FRINGE (RFNG), a glycosyltransferase that enhances the affinity of NOTCH1 to JAG1 incis. Finally, our results indicate that NOTCH1-JAG1cis-inhibition is necessary for preserving the HSC phenotype in the hematopoietic clusters of the aorta.

Early organogenesis represents a key step in animal development, where pluripotent cells divide and diversify to initiate formation of all major organs. Here, we used scRNA-Seq to profile over 300,000 single cell transcriptomes sampled in 6 hour intervals from mouse embryos between E8.5 and E9.5. Combining this dataset with our previous E6.5 to E8.5 atlas resulted in a densely sampled time course of over 400,000 cells from early gastrulation to organogenesis. Computational lineage reconstruction at full organismal scale identified complex waves of blood and endothelial development, including a new molecular programme for somite-derived endothelium. To assess developmental fates across the primitive streak, we dissected the E7.5 primitive streak into four adjacent regions, performed scRNA- Seq and predicted cell fates computationally. We next defined early developmental state/fate relationships experimentally by a combination of orthotopic grafting, microscopic analysis of graft contribution as well as scRNA-Seq to transcriptionally determine cell fates of the grafted primitive streak regions after 24h ofin vitroembryo culture. Experimentally determined fate outcomes were in good agreement with the fates predicted computationally, thus demonstrating how classical grafting experiments can be revisited to establish high-resolution cell state/fate relationships. Such interdisciplinary approaches will benefit future studies in both developmental biology as well as guide thein vitroproduction of cells for organ regeneration and repair.

Hematopoietic aging is marked by a loss of regenerative capacity and skewed differentiation from hematopoietic stem cells (HSC) leading to dysfunctional blood production. Signals from the bone marrow (BM) niche dynamically tailor hematopoiesis, but the effect of aging on the niche microenvironment and the contribution of the aging niche to blood aging still remains unclear. Here, we characterize the inflammatory milieu in the aged marrow cavity that drives both stromal and hematopoietic remodeling. We find decreased numbers and functionality of osteogenic mesenchymal stromal cells (MSC) at the endosteum and expansion of pro-inflammatory perisinusoidal MSCs with deterioration of sinusoidal endothelium in the central marrow, which together create a degraded and inflamed old niche. Molecular mapping at single cell resolution confirms disruption of cell identities and enrichment of inflammatory response genes in niche populations. Niche inflammation, in turn, drives chronic activation of emergency myelopoiesis pathways in old HSCs and multipotent progenitors (MPP), which promotes myeloid differentiation at the expense of lymphoid and erythroid commitment and hinders hematopoietic regeneration. Remarkably, niche deterioration, HSC dysfunction and defective hematopoietic regeneration, can be improved by blocking inflammatory IL-1 signaling. Our results demonstrate that targeting niche inflammation is a tractable strategy to restore blood production during aging.

T1 Figure 1ConclusionsOverall, this largest paediatric single cell COVID-19 study to date showed significant differences in response to SARS-CoV-2 between children and adults, reflecting the changes of the immune landscape over developmental time, which in children are dominated by naive and innate responses.

Clotting Factor V (FV) is primarily synthesized in the liver and when cleaved by thrombin forms pro-coagulant Factor Va (FVa). Using whole blood RNAseq and scRNAseq of peripheral blood mononuclear cells, we find that FV mRNA is expressed in leukocytes, and identify neutrophils, monocytes, and T regulatory cells as sources of increased FV in hospitalized patients with COVID-19. Proteomic analysis confirms increased FV in circulating neutrophils in severe COVID-19, and immunofluorescence microscopy identifies FV in lung-infiltrating leukocytes in COVID-19 lung disease. Increased leukocyte FV expression in severe disease correlates with T-cell lymphopenia. Both plasma-derived and a cleavage resistant recombinant FV, but not thrombin cleaved FVa, suppress T-cell proliferation

Common genetic variants modulate the cellular response to viruses and are implicated in a range of immune pathologies, including infectious and autoimmune diseases. The transcriptional antiviral response is known to vary between infected cells from a single individual, yet how genetic variants across individuals modulate the antiviral response (and its cell-to-cell variability) is not well understood. Here, we triggered the antiviral response in human fibroblasts from 68 healthy donors, and profiled tens of thousands of cells using single-cell RNA-seq. We developed GASPACHO (GAuSsian Processes for Association mapping leveraging Cell HeterOgeneity), the first statistical approach designed to identify dynamic eQTLs across a transcriptional trajectory of cell populations, without aggregating single-cell data into pseudo-bulk. This allows us to uncover the underlying architecture and variability of antiviral response across responding cells, and to identify more than two thousands eQTLs modulating the dynamic changes during this response. Many of these eQTLs colocalise with risk loci identified in GWAS of infectious and autoimmune diseases. As a case study, we focus on a COVID-19 susceptibility locus, colocalised with the antiviral OAS1 splicing QTL. We validated it in blood cells from a patient cohort and in the infected nasal cells of a patient with the risk allele, demonstrating the utility of GASPACHO to fine-map and functionally characterise a genetic locus. In summary, our novel analytical approach provides a new framework for delineation of the genetic variants that shape a wide spectrum of transcriptional responses at single-cell resolution.

The study of hematopoiesis has been revolutionized in recent years by the application of single-cell RNA sequencing technologies. The technique coupled with rapidly developing bioinformatic analysis has provided great insight into the cell type compositions of many populations previously defined by their cell surface phenotype. Moreover, transcriptomic information enables the identification of individual molecules and pathways which define novel cell populations and their transitions including cell lineage decisions. Combining single-cell transcriptional profiling with molecular perturbations allows functional analysis of individual factors in gene regulatory networks and better understanding of the earliest stages of malignant transformation. In this chapter we describe a comprehensive protocol for scRNA-Seq analysis of the mouse bone marrow, using both plate-based (low throughput) and droplet-based (high throughput) methods. The protocol includes instructions for sample preparation, an antibody panel for flow cytometric purification of hematopoietic progenitors with index sorting for plate-based analysis or in bulk for droplet-based methods. The plate-based protocol described in this chapter is a combination of the Smart-Seq2 and mcSCRB-Seq protocols, optimized in our laboratory. It utilizes off-the-shelf reagents for cDNA preparation, is amenable to automation using a liquid handler, and takes 4 days from preparation of the cells for sorting to producing a sequencing-ready library. The droplet-based method (using for instance the 10× Genomics platform) relies on the manufacturer's user guide and commercial reagents, and takes 3 days from isolation of the cells to the production of a library ready for sequencing.

While a substantial proportion of adults infected with SARS-CoV-2 progress to develop severe disease, children rarely manifest respiratory complications. Therefore, understanding differences in the local and systemic response to SARS-CoV-2 infection between children and adults may provide important clues about the pathogenesis of SARS-CoV-2 infection. To address this, we first generated a healthy reference multi-omics single cell data set from children (n=30) in whom we have profiled triple matched samples: nasal and tracheal brushings and PBMCs, where we track the developmental changes for 42 airway and 31 blood cell populations from infancy, through childhood to adolescence. This has revealed the presence of naive B and T lymphocytes in neonates and infants with a unique gene expression signature bearing hallmarks of innate immunity. We then contrast the healthy reference with equivalent data from severe paediatric and adult COVID-19 patients (total n=27), from the same three types of samples: upper and lower airways and blood. We found striking differences: children with COVID-19 as opposed to adults had a higher proportion of innate lymphoid and non-clonally expanded naive T cells in peripheral blood, and a limited interferon-response signature. In the airway epithelium, we found the highest viral load in goblet and ciliated cells and describe a novel inflammatory epithelial cell population. These cells represent a transitional regenerative state between secretory and ciliated cells; they were found in healthy children and were enriched in paediatric and adult COVID-19 patients. Epithelial cells display an antiviral and neutrophil-recruiting gene signature that is weaker in severe paediatricversusadult COVID-19. Our matched blood and airway samples allowed us to study the spatial dynamics of infection. Lastly, we provide a user-friendly interface for this data1as a highly granular reference for the study of immune responses in airways and blood in children.

Hematopoietic aging is marked by a loss of regenerative capacity and skewed differentiation from hematopoietic stem cells (HSC) leading to dysfunctional blood production. Signals from the bone marrow (BM) niche dynamically tailor hematopoiesis, but the effect of aging on the niche microenvironment and the contribution of the aging niche to blood aging still remains unclear. Here, we characterize the inflammatory milieu in the aged marrow cavity that drives both stromal and hematopoietic remodeling. We find decreased numbers and functionality of osteogenic mesenchymal stromal cells (MSC) at the endosteum and expansion of pro-inflammatory perisinusoidal MSCs with deterioration of sinusoidal endothelium in the central marrow, which together create a degraded and inflamed old niche. Molecular mapping at single cell resolution confirms disruption of cell identities and enrichment of inflammatory response genes in niche populations. Niche inflammation, in turn, drives chronic activation of emergency myelopoiesis pathways in old HSCs and multipotent progenitors (MPP), which promotes myeloid differentiation at the expense of lymphoid and erythroid commitment and hinders hematopoietic regeneration. Remarkably, niche deterioration, HSC dysfunction and defective hematopoietic regeneration, can be improved by blocking inflammatory IL-1 signaling. Our results demonstrate that targeting niche inflammation is a tractable strategy to restore blood production during aging.

The COVID-19 pandemic, caused by SARS coronavirus 2 (SARS-CoV-2), has resulted in excess morbidity and mortality as well as economic decline. To characterise the systemic host immune response to SARS-CoV-2, we performed single-cell RNA-sequencing coupled with analysis of cell surface proteins, providing molecular profiling of over 800,000 peripheral blood mononuclear cells from a cohort of 130 patients with COVID-19. Our cohort, from three UK centres, spans the spectrum of clinical presentations and disease severities ranging from asymptomatic to critical. Three control groups were included: healthy volunteers, patients suffering from a non-COVID-19 severe respiratory illness and healthy individuals administered with intravenous lipopolysaccharide to model an acute inflammatory response. Full single cell transcriptomes coupled with quantification of 188 cell surface proteins, and T and B lymphocyte antigen receptor repertoires have provided several insights into COVID-19: 1. a new non-classical monocyte state that sequesters platelets and replenishes the alveolar macrophage pool; 2. platelet activation accompanied by early priming towards megakaryopoiesis in immature haematopoietic stem/progenitor cells and expansion of megakaryocyte-primed progenitors; 3. increased clonally expanded CD8+ effector:effector memory T cells, and proliferating CD4+ and CD8+ T cells in patients with more severe disease; and 4. relative increase of IgA plasmablasts in asymptomatic stages that switches to expansion of IgG plasmablasts and plasma cells, accompanied with higher incidence of BCR sharing, as disease severity increases. All data and analysis results are available for interrogation and data mining through an intuitive web portal. Together, these data detail the cellular processes present in peripheral blood during an acute immune response to COVID-19, and serve as a template for multi-omic single cell data integration across multiple centers to rapidly build powerful resources to help combat diseases such as COVID-19.

Severe Covid-19 is associated with elevated plasma Factor V (FV) and increased risk of thromboembolism. We report that neutrophils, T regulatory cells (Tregs), and monocytes from patients with severe Covid-19 express FV, and expression correlates with T cell lymphopenia. In vitro full length FV, but not FV activated by thrombin cleavage, suppresses T cell proliferation. Increased and prolonged FV expression by cells of the innate and adaptive immune systems may contribute to lymphopenia in severe Covid-19. Activation by thrombin destroys the immunosuppressive properties of FV. Anticoagulation in Covid-19 patients may have the unintended consequence of suppressing the adaptive immune system.

Chimeric antigen receptor (CAR) T-cell adoptive therapy is set to transform the treatment of a rapidly expanding range of malignancies. Although the activation process of normal T cells is well characterized, comparatively little is known about the activation of cells via the CAR. Here we have used flow cytometry together with single-cell transcriptome profiling to characterize the starting material (peripheral blood mononuclear cells) and CAR therapeutic products of 3 healthy donors in the presence and absence of antigen-specific stimulation. Analysis of 53,191 single-cell transcriptomes showed APRIL-based CAR products to contain several subpopulations of cells, with cellular composition reproducible from donor to donor, and all major cellular subsets compatible with CAR expression. Only 50% of CAR-expressing cells displayed transcriptional changes upon CAR-specific antigen exposure. The resulting molecular signature for CAR T-cell activation provides a rich resource for future dissection of underlying mechanisms. Targeted data interrogation also revealed that a small proportion of antigen-responding CAR-expressing cells displayed an exhaustion signature, with both known markers and genes not previously associated with T-cell exhaustion. Comprehensive single-cell transcriptomic analysis thus represents a powerful way to guide the assessment and optimization of clinical-grade CAR-T-cells, and inform future research into the underlying molecular processes.

The ability of epithelial cells to rewire their cell fate program beyond their physiological repertoire has become a new paradigm in stem cell biology. This plasticity leaves behind the concept of strict stem cell hierarchies, opening up new exciting questions about its limits and underlying regulation. Here we developed a heterotypic 3D culture system to study the mechanisms modulating changes in the identity of adult esophageal epithelial cells. We demonstrate that, when exposed to the foreign stroma of adult skin, esophageal cells transition towards hair follicle identity and architecture. Heterotypic transplantation experiments recapitulated this cell fate conversion processin vivo. Single-cell RNA sequencing and histological analysis, capturing the temporality of this process, reveal that most esophageal cells switching towards skin identity remain in an intermediate state marked by a transient regenerative profile and a particularly strong hypoxic signature. Inhibition of HIF1a establishes the central role of this pathway in regulating epithelial cell plasticity, driving cells away from their transition state in favor of cell fate conversion.

The study of hematopoiesis has been revolutionized in recent years by the application of single-cell RNA sequencing technologies. The technique coupled with rapidly developing bioinformatic analysis has provided great insight into the cell type compositions of many populations previously defined by their cell surface phenotype. Moreover, transcriptomic information enables the identification of individual molecules and pathways which define novel cell populations and their transitions including cell lineage decisions. Combining single-cell transcriptional profiling with molecular perturbations allows functional analysis of individual factors in gene regulatory networks and better understanding of the earliest stages of malignant transformation. In this chapter we describe a comprehensive protocol for scRNA-Seq analysis of the mouse bone marrow, using both plate-based (low throughput) and droplet-based (high throughput) methods. The protocol includes instructions for sample preparation, an antibody panel for flow cytometric purification of hematopoietic progenitors with index sorting for plate-based analysis or in bulk for droplet-based methods. The plate-based protocol described in this chapter is a combination of the Smart-Seq2 and mcSCRB-Seq protocols, optimized in our laboratory. It utilizes off-the-shelf reagents for cDNA preparation, is amenable to automation using a liquid handler, and takes 4 days from preparation of the cells for sorting to producing a sequencing-ready library. The droplet-based method (using for instance the 10× Genomics platform) relies on the manufacturer's user guide and commercial reagents, and takes 3 days from isolation of the cells to the production of a library ready for sequencing.

Summary Whereas hundreds of cells in the mouse embryonic aorta transdifferentiate to hematopoietic cells, only very few establish hematopoietic stem cell (HSC) identity at a single time point. The Gata2 transcription factor is essential for HSC generation and function. In contrast to surface-marker-based cell isolation, Gata2-based enrichment provides a direct link to the internal HSC regulatory network. Here, we use iterations of index-sorting of Gata2-expressing intra-aortic hematopoietic cluster (IAHC) cells, single-cell transcriptomics, and functional analyses to connect HSC identity to specific gene expression. Gata2-expressing IAHC cells separate into 5 major transcriptomic clusters. Iterative analyses reveal refined CD31, cKit, and CD27 phenotypic parameters that associate specific molecular profiles in one cluster with distinct HSC and multipotent progenitor function. Thus, by iterations of single-cell approaches, we identify the transcriptome of the first functional HSCs as they emerge in the mouse embryo and localize them to aortic clusters containing 1–2 cells., Graphical Abstract, Highlights • Single-cell iterations capture the transcriptome of the first functional HSCs in mouse • Gata2 connects the “inner” HSC regulatory network with single-cell function • Specific Gata2, cKit, and CD27 levels define all HSCs in embryonic aortic clusters • HSCs emerge from endothelium as single cells within aortic clusters of 1–2 cells, Vink and colleagues capture the transcriptome of the first functional HSCs in mouse by single-cell RNA-seq, index-sorting, and in vivo and invitro hematopoietic analyses. The HSC transcriptome is unique compared to HPCs, and heterogeneous expression of pivotal genes suggests that establishment of functional hematopoietic fate during cell emergence from embryonic aortic endothelium is stochastic.

During mouse embryonic development, pluripotent cells rapidly divide and diversify, yet the regulatory programs that define the cell repertoire for each organ remain ill-defined. To delineate comprehensive chromatin landscapes during early organogenesis, we mapped chromatin accessibility in 19,453 single nuclei from mouse embryos at 8.25 days post-fertilization. Identification of cell-type-specific regions of open chromatin pinpointed two TAL1-bound endothelial enhancers, which we validated using transgenic mouse assays. Integrated gene expression and transcription factor motif enrichment analyses highlighted cell-type-specific transcriptional regulators. Subsequent in vivo experiments in zebrafish revealed a role for the ETS factor FEV in endothelial identity downstream of ETV2 (Etsrp in zebrafish). Concerted in vivo validation experiments in mouse and zebrafish thus illustrate how single-cell open chromatin maps, representative of a mammalian embryo, provide access to the regulatory blueprint for mammalian organogenesis.

Hematopoietic stem cells (HSCs) develop from the hemogenic endothelium in cluster structures that protrude into the embryonic aortic lumen. Although much is known about the molecular characteristics of the developing hematopoietic cells, we lack a complete understanding of their origin and the three-dimensional organization of the niche. Here, we use advanced live imaging techniques of organotypic slice cultures, clonal analysis, and mathematical modeling to show the two-step process of intra-aortic hematopoietic cluster (IACH) formation. First, a hemogenic progenitor buds up from the endothelium and undergoes division forming the monoclonal core of the IAHC. Next, surrounding hemogenic cells are recruited into the IAHC, increasing their size and heterogeneity. We identified the Notch ligand Dll4 as a negative regulator of the recruitment phase of IAHC. Blocking of Dll4 promotes the entrance of new hemogenic Gfi1+ cells into the IAHC and increases the number of cells that acquire HSC activity. Mathematical modeling based on our data provides estimation of the cluster lifetime and the average recruitment time of hemogenic cells to the cluster under physiologic and Dll4-inhibited conditions. We thank Thomas Graf (CRG), Minhong Yan (Genentech), Georges Lacaud (CRUK Manchester Institut) and Ralf Adams (Max Planck Institute for Molecular Medicine, Münster) for important mice models and reagents. Thanks to Catherine Robin, Timo Zimmermann and, Raul Gómez-Riera for imaging advice. Thanks to the PRBB facilities: animal facility, flow cytometry, CNAGCRG genomics, confocal microscopy for critical technical assistance. Thanks to Marta Garrido and Luis Galán for experimental support and all the members of the lab for helpful critical discussions. This research was funded by the Ministerio de Economía y Competitividad (SAF2016 -75613-R) and Generalitat de Catalunya, Agència de Gestió d’Ajuds Universitaris i de Recerca (AGAUR) (2017 SGR 135) and PERIS-SLT002/16/00299 to AB and ERC AdG 341096 to ED. CP was a recipient of Juan de la Cierva fellowship (FJCI2014-19870), RT is a recipient of Beatriu de Pinos (2016 BP 00021), JG is a recipient of PERIS (SLT002/16/00070).

By profiling the transcriptomes of individual cells, single-cell RNA sequencing provides unparalleled resolution to study cellular heterogeneity. However, this comes at the cost of high technical noise, including cell-specific biases in capture efficiency and library generation. One strategy for removing these biases is to add a constant amount of spike-in RNA to each cell and to scale the observed expression values so that the coverage of spike-in transcripts is constant across cells. This approach has previously been criticized as its accuracy depends on the precise addition of spike-in RNA to each sample. Here, we perform mixture experiments using two different sets of spike-in RNA to quantify the variance in the amount of spike-in RNA added to each well in a plate-based protocol. We also obtain an upper bound on the variance due to differences in behavior between the two spike-in sets. We demonstrate that both factors are small contributors to the total technical variance and have only minor effects on downstream analyses, such as detection of highly variable genes and clustering. Our results suggest that scaling normalization using spike-in transcripts is reliable enough for routine use in single-cell RNA sequencing data analyses.

Although many aspects of blood production are well understood, the spatial organization of myeloid differentiation in the bone marrow remains unknown. Here we use imaging to track granulocyte/macrophage progenitor (GMP) behaviour in mice during emergency and leukaemic myelopoiesis. In the steady state, we find individual GMPs scattered throughout the bone marrow. During regeneration, we observe expanding GMP patches forming defined GMP clusters, which, in turn, locally differentiate into granulocytes. The timed release of important bone marrow niche signals (SCF, IL-1β, G-CSF, TGFβ and CXCL4) and activation of an inducible $\textit{Irf8}$ and β-catenin progenitor self-renewal network control the transient formation of regenerating GMP clusters. In leukaemia, we show that GMP clusters are constantly produced owing to persistent activation of the self-renewal network and a lack of termination cytokines that normally restore haematopoietic stem-cell quiescence. Our results uncover a previously unrecognized dynamic behaviour of GMPs $\textit{in situ}$, which tunes emergency myelopoiesis and is hijacked in leukaemia.

Leukaemogenic mutations commonly disrupt cellular differentiation and/or enhance proliferation, thus perturbing the regulatory programs that control self-renewal and differentiation of stem and progenitor cells. Translocations involving the Mll1 (Kmt2a) gene generate powerful oncogenic fusion proteins, predominantly affecting infant and paediatric AML and ALL patients. The early stages of leukaemogenic transformation are typically inaccessible from human patients and conventional mouse models. Here, we take advantage of cells conditionally blocked at the multipotent haematopoietic progenitor stage to develop a MLL-r model capturing early cellular and molecular consequences of MLL-ENL expression based on a clear clonal relationship between parental and leukaemic cells. Through a combination of scRNA-seq, ATAC-seq and genome-scale CRISPR-Cas9 screening, we identify pathways and genes likely to drive the early phases of leukaemogenesis. Finally, we demonstrate the broad utility of using matched parental and transformed cells for small molecule inhibitor studies by validating both previously known and other potential therapeutic targets., The oncogene MLL is frequently translocated in leukemia, resulting in oncogenic fusion proteins. Here, the authors report a temporally controlled mouse model of MLL-ENL driven leukemia AND identify therapeutic targets associated with early MLL-ENL driven leukaemogenesis.

Hematopoietic stem cells (HSCs) develop from the hemogenic endothelium in cluster structures that protrude into the embryonic aortic lumen. Although much is known about the molecular characteristics of the developing hematopoietic cells, we lack a complete understanding of their origin and the three-dimensional organization of the niche. Here we use advanced live imaging techniques of organotypic slice cultures, clonal analysis, and mathematical modelling to show the two-step process of intra-aortic hematopoietic cluster (IACH) formation. First, a hemogenic progenitor buds up from the endothelium and undergoes division forming the monoclonal core of the IAHC. Next, surrounding hemogenic cells are recruited into the IAHC, increasing their size and heterogeneity. We identified the Notch ligand Dll4 as a negative regulator of the recruitment phase of IAHC. Blocking of Dll4 promotes the entrance of new hemogenic Gfi1+ cells into the IAHC and increases the number of cells that acquire HSC activity. Mathematical modelling based on our data provides estimation of the cluster lifetime and the average recruitment time of hemogenic cells to the cluster under physiologic and Dll4-inhibited conditions.

Summary In vertebrates, GATA2 is a master regulator of hematopoiesis and is expressed throughout embryo development and in adult life. Although the essential role of GATA2 in mouse hematopoiesis is well established, its involvement during early human hematopoietic development is not clear. By combining time-controlled overexpression of GATA2 with genetic knockout experiments, we found that GATA2, at the mesoderm specification stage, promotes the generation of hemogenic endothelial progenitors and their further differentiation to hematopoietic progenitor cells, and negatively regulates cardiac differentiation. Surprisingly, genome-wide transcriptional and chromatin immunoprecipitation analysis showed that GATA2 bound to regulatory regions, and repressed the expression of cardiac development-related genes. Moreover, genes important for hematopoietic differentiation were upregulated by GATA2 in a mostly indirect manner. Collectively, our data reveal a hitherto unrecognized role of GATA2 as a repressor of cardiac fates, and highlight the importance of coordinating the specification and repression of alternative cell fates., Graphical Abstract, Highlights • GATA2 promotes hemogenic emergence during human iPSC differentiation • GATA2 enhances hematopoietic differentiation of human iPSCs • GATA2 acts as a direct repressor of cardiac fates during mesoderm specification, Giorgetti A. and colleagues demonstrate that GATA2 induction in early mesodermal cells leads to the robust generation of hemogenic endothelial and hematopoietic cells from human iPSCs. They show the ability of GATA2 to instruct mesodermal precursors toward a hematopoietic cell fate and concurrently inhibit cardiac fates. This study expands our understanding of the regulatory networks that control early human hematopoiesis.

Across the animal kingdom, gastrulation represents a key developmental event during which embryonic pluripotent cells diversify into lineage-specific precursors that will generate the adult organism. Here we report the transcriptional profiles of 116,312 single cells from mouse embryos collected at nine sequential time points ranging from 6.5 to 8.5 days post-fertilization. We construct a molecular map of cellular differentiation from pluripotency towards all major embryonic lineages, and explore the complex events involved in the convergence of visceral and primitive streak-derived endoderm. Furthermore, we use single-cell profiling to show that Tal1−/− chimeric embryos display defects in early mesoderm diversification, and we thus demonstrate how combining temporal and transcriptional information can illuminate gene function. Together, this comprehensive delineation of mammalian cell differentiation trajectories in vivo represents a baseline for understanding the effects of gene mutations during development, as well as a roadmap for the optimization of in vitro differentiation protocols for regenerative medicine. Single-cell profiling is used to create a molecular-level atlas of cell differentiation trajectories during gastrulation and early organogenesis in the mouse.

Cell fate and behavior are regulated by dynamic intrinsic and extrinsic processes. These processes affect single cells to produce a fully integrated organism with distinct functional cell types and tissues. Experimental hematology provides insight into the independent iterations of blood cell development (primitive, definitive progenitor and definitive stem) in the different embryonic sites/times prior to adult hematopoietic system formation. Hematopoietic stem cell (HSC) generation occurs in a short developmental time window from a natural transdifferentiation of aortic endothelial-to-hematopoietic cells (EHT). The Gata2 transcription factor is a dosage-dependent intrinsic regulator of EHT controlled by signals from the extrinsic developmental niche. Hematopoietic progenitor cells (HPC) are also generated by EHT, although at an earlier stage. So, what mechanistically determines HPC or HSC fate during EHT? Vital time-lapse imaging reveals highly dynamic Gata2 expression changes in cells emerging during EHT. We used Gata2 levels to examine the “inner” HSC transcriptional regulatory network. Single cell iterative RNAseq, index-sorting and in vivo/vitro hematopoietic analyses show that the transcriptome of the first aortic HSCs is unique compared to HPCs. The observed heterogeneity and quantitative expression differences of pivotal EHT/HSC genes suggest that establishment of functional hematopoietic fate is a stochastic process. Additionally, dynamic changes occur in the aortic niche. Transient cell-cell interactions occur between motile yolk sac HPC-derived macrophages and emerging aortic hematopoietic cells to affect hematopoietic output/function. The extrinsic signals from this new pro-inflammatory macrophage subset, together with intrinsic dynamic Gata2 expression changes in endothelial cells could assist strategies for ex vivo HSC generation.

SummaryDespite great advances in understanding the mechanisms underlying blood production, lineage specification at the level of multipotent progenitors (MPPs) remains poorly understood. Here, we show that MPP2 and MPP3 are distinct myeloid-biased MPP subsets that work together with lymphoid-primed MPP4 cells to control blood production. We find that all MPPs are produced in parallel by hematopoietic stem cells (HSCs), but with different kinetics and at variable levels depending on hematopoietic demands. We also show that the normally rare myeloid-biased MPPs are transiently overproduced by HSCs in regenerating conditions, hence supporting myeloid amplification to rebuild the hematopoietic system. This shift is accompanied by a reduction in self-renewal activity in regenerating HSCs and reprogramming of MPP4 fate toward the myeloid lineage. Our results support a dynamic model of blood development in which HSCs convey lineage specification through independent production of distinct lineage-biased MPP subsets that, in turn, support lineage expansion and differentiation.

During gastrulation, cell types from all three germ layers are specified and the basic body plan is established 1 . However, molecular analysis of this key developmental stage has been hampered by limited cell numbers and a paucity of markers. Single-cell RNA sequencing circumvents these problems, but has so far been limited to specific organ systems 2 . Here, we report single-cell transcriptomic characterization of >20,000 cells immediately following gastrulation at E8.25 of mouse development. We identify 20 major cell types, which frequently contain substructure, including three distinct signatures in early foregut cells. Pseudo-space ordering of somitic progenitor cells identifies dynamic waves of transcription and candidate regulators, which are validated by molecular characterization of spatially resolved regions of the embryo. Within the endothelial population, cells that transition from haemogenic endothelial to erythro-myeloid progenitors specifically express Alox5 and its co-factor Alox5ap, which control leukotriene production. Functional assays using mouse embryonic stem cells demonstrate that leukotrienes promote haematopoietic progenitor cell generation. Thus, this comprehensive single-cell map can be exploited to reveal previously unrecognized pathways that contribute to tissue development.

By profiling the transcriptomes of individual cells, single-cell RNA sequencing provides unparalleled resolution to study cellular heterogeneity. However, this comes at the cost of high technical noise, including cell-specific biases in capture efficiency and library generation. One strategy for removing these biases is to add a constant amount of spike-in RNA to each cell, and to scale the observed expression values so that the coverage of spike-in RNA is constant across cells. This approach has previously been criticized as its accuracy depends on the precise addition of spike-in RNA to each sample, and on similarities in behaviour (e.g., capture efficiency) between the spike-in and endogenous transcripts. Here, we perform mixture experiments using two different sets of spike-in RNA to quantify the variance in the amount of spike-in RNA added to each well in a plate-based protocol. We also obtain an upper bound on the variance due to differences in behaviour between the two spike-in sets. We demonstrate that both factors are small contributors to the total technical variance and have only minor effects on downstream analyses such as detection of highly variable genes and clustering. Our results suggest that spike-in normalization is reliable enough for routine use in single-cell RNA sequencing data analyses.

Objectives: To identify sub-populations of intestinal preproglucagon-expressing (PPG) cells producing Glucagon-like Peptide-1, and their associated expression profiles of sensory receptors, thereby enabling the discovery of therapeutic strategies that target these cell populations for the treatment of diabetes and obesity. Methods: We performed single cell RNA sequencing of PPG-cells purified by flow cytometry from the upper small intestine of 3 GLU-Venus mice. Cells from 2 mice were sequenced at low depth, and from the third mouse at high depth. High quality sequencing data from 234 PPG-cells were used to identify clusters by tSNE analysis. qPCR was performed to compare the longitudinal and crypt/villus locations of cluster-specific genes. Immunofluorescence and mass spectrometry were used to confirm protein expression. Results: PPG-cells formed 3 major clusters: a group with typical characteristics of classical L-cells, including high expression of Gcg and Pyy (comprising 51% of all PPG-cells); a cell type overlapping with Gip-expressing K-cells (14%); and a unique cluster expressing Tph1 and Pzp that was predominantly located in proximal small intestine villi and co-produced 5-HT (35%). Expression of G-protein coupled receptors differed between clusters, suggesting the cell types are differentially regulated, and would be differentially targetable. Conclusions: Our findings support the emerging concept that many enteroendocrine cell populations are highly overlapping, with individual cells producing a range of peptides previously assigned to distinct cell types. Different receptor expression profiles across the clusters highlight potential drug targets to increase gut hormone secretion for the treatment of diabetes and obesity., Research in the FR/FMG lab is funded by a Wellcome joint investigator award (106262/Z/14/Z and 106263/Z/14/Z) and a joint MRC programme within the Metabolic Diseases Unit (MRC_MC_UU_12012/3). Single cell collection and analysis was supported through MRC Clinical Research Infrastructure funds for the Cambridge Single Cell Analysis Clinical Core Facility. Work in the Göttgens laboratory is supported by grants from the Wellcome Trust, Bloodwise, Cancer Research UK, NIDDK and core support grants by the Wellcome Trust to the Wellcome Trust-MRC Cambridge Stem Cell Institute.

Hematopoietic stem cells (HSC) are responsible for life-long maintenance and regeneration of the adult vertebrate blood system. The first HSCs arise from intra-aortic hematopoietic cluster cells (IAHC) in the mouse aorta at embryonic day (E)10.5 through a transdifferentiation process called endothelial-to-hematopoietic-transition (EHT). Previously we showed that the heptad factor Gata2, one of a group of 7 (heptad) transcription factors pivotal to definitive hematopoiesis, is required for HSC generation and survival. Gata2-/- embryos suffer lethality at E10.5 and Gata2+/- HSCs are qualitatively defective. Utilizing our Gata2Venus reporter mouse for time-lapse imaging of EHT in vivo, we found rapid pulsatile Gata2 expression level changes in single transitioning cells, implicating transcriptional instability during establishment of hematopoietic fate. Flow cytometric analysis of IAHCs revealed cells with varying levels of Venus (Vmedium, Vhigh). Hematopoietic progenitors (HPC) were highly enriched and HSCs exclusively found in the Vmed fraction. From these results, we hypothesized that fate choice is based on stochastic expression (levels and combinations) of the pivotal heptad factors. Here we test this by single cell transcriptomics and functional assays. Single-cell RNA sequencing of index-sorted Vmed cells shows surprising heterogeneity in heptad gene expression and has facilitated further HSC enrichment, defining intensity levels of known markers and providing additional markers. These new markers were functionally validated and their expression visualized by immunohistochemistry. Single-cell HPC assays also allowed for enrichment of immature HPCs. Further sorting and sequencing reiterations will refine the HSC transcriptomic profile and provide information on whether transcriptome heterogeneity correlates with and plays a role in generating both HSCs and HPCs.

Small molecule inhibition of the BET family of proteins, which bind acetylated lysines within histones, has been shown to have a marked therapeutic benefit in pre-clinical models of mixed lineage leukemia (MLL) fusion protein-driven leukemias. Here, we report that I-BET151, a highly specific BET family bromodomain inhibitor, leads to growth inhibition in a human erythroleukemic (HEL) cell line as well as in erythroid precursors isolated from polycythemia vera patients. One of the genes most highly downregulated by I-BET151 was LMO2, an important oncogenic regulator of hematopoietic stem cell development and erythropoiesis. We previously reported that LMO2 transcription is dependent upon Janus kinase 2 (JAK2) kinase activity in HEL cells. Here, we show that the transcriptional changes induced by a JAK2 inhibitor (TG101209) and I-BET151 in HEL cells are significantly over-lapping, suggesting a common pathway of action. We generated JAK2 inhibitor resistant HEL cells and showed that these retain sensitivity to I-BET151. These data highlight I-BET151 as a potential alternative treatment against myeloproliferative neoplasms driven by constitutively active JAK2 kinase.

Cellular decision-making is mediated by a complex interplay of external stimuli with the intracellular environment, in particular transcription factor regulatory networks. Here we have determined the expression of a network of 18 key haematopoietic transcription factors in 597 single primary blood stem and progenitor cells isolated from mouse bone marrow. We demonstrate that different stem/progenitor populations are characterized by distinctive transcription factor expression states, and through comprehensive bioinformatic analysis reveal positively and negatively correlated transcription factor pairings, including previously unrecognized relationships between Gata2, Gfi1 and Gfi1b. Validation using transcriptional and transgenic assays confirmed direct regulatory interactions consistent with a regulatory triad in immature blood stem cells, where Gata2 may function to modulate cross-inhibition between Gfi1 and Gfi1b. Single-cell expression profiling therefore identifies network states and allows reconstruction of network hierarchies involved in controlling stem cell fate choices, and provides a blueprint for studying both normal development and human disease. © 2013 Macmillan Publishers Limited. All rights reserved.

LMO1 is a transcriptional regulator and a T-acute lymphoblastic leukaemia (T-ALL) oncogene. Although first identified in association with a chromosomal translocation in T-ALL, the ectopic expression of LMO1 occurs far more frequently in the absence of any known mutation involving its locus. Given that LMO1 is barely expressed in any haematopoietic lineage, and activation of transcriptional drivers in leukaemic cells is not well described, we investigated the regulation of this gene in normal haematopoietic and leukaemic cells. We show that LMO1 has two promoters that drive reporter gene expression in transgenic mice to neural tissues known to express endogenous LMO1. The LMO1 promoters display bivalent histone marks in multiple blood lineages including T-cells, and a 3' flanking region at LMO1 +57 contains a transcriptional enhancer that is active in developing blood cells in transgenic mouse embryos. The LMO1 promoters become activated in T-ALL together with the 3' enhancer, which is bound in primary T-ALL cells by SCL/TAL1 and GATA3. Taken together, our results show that LMO1 is poised for expression in normal progenitors, where activation of SCL/TAL1 together with a breakdown of epigenetic repression of LMO1 regulatory elements induces ectopic LMO1 expression that contributes to the development and maintenance of T-ALL.

The ETS transcription factor Elf-4 is an important regulator of hematopoietic stem cell (HSC) and T cell homeostasis. To gain insights into the transcriptional circuitry within which Elf-4 operates, we used comparative sequence analysis coupled with chromatin immunoprecipitation (ChIP) with microarray technology (ChIP-chip) assays for specific chromatin marks to identify three promoters and two enhancers active in hematopoietic and endothelial cell lines. Comprehensive functional validation of each of these regulatory regions in transgenic mouse embryos identified a tissue-specific enhancer (−10E) that displayed activity in fetal liver, dorsal aorta, vitelline vessels, yolk sac, and heart. Integration of a ChIP-sequencing (ChIP-Seq) data set for 10 key stem cell transcription factors showed Pu.1, Fli-1, and Erg were bound to the −10E element, and mutation of three highly conserved ETS sites within the enhancer abolished its activity. Finally, the transcriptional repressor Gfi1b was found to bind to and repress one of the Elf-4 promoters (−30P), and we show that this repression of Elf-4 is important for the maturation of primary fetal liver erythroid cells. Taken together, our results provide a comprehensive overview of the transcriptional control of Elf-4 within the hematopoietic system and, thus, integrate Elf-4 into the wider transcriptional regulatory networks that govern hematopoietic development.

The C/EBPα transcription factor plays a pivotal role in myeloid differentiation and E2F-mediated cell cycle regulation. Although CEBPA mutations are common in acute myeloid leukaemia (AML), little is known regarding pre-leukemic alterations caused by mutated CEBPA. Here, we investigated early events involved in pre-leukemic transformation driven by CEBPA N321D in the LMPP-like cell line Hoxb8-FL (Redecke et al., Nat Methods 2013), which can be maintained in vitro as a self-renewing LMPP population using Flt3L and estradiol, as well as differentiated both in vitro and in vivo into myeloid and lymphoid cell types. Hoxb8-FL cells were retrovirally transduced with Empty Vector (EV), wild-type CEBPA (CEBPA WT) or its N321D mutant form (CEBPA N321D). CEBPA WT-transduced cells showed increased expression of cd11b and SIRPα and downregulation of c-kit, suggesting that wild-type CEBPA was sufficient to promote differentiation even under LMPP growth conditions. Interestingly, we did not observe the same phenotype in CEBPA N321D-transduced cells. Upon withdrawal of estradiol, both EV and CEBPA WT-transduced cells differentiated rapidly into a conventional dendritic cell (cDC) phenotype by day 7 and died within 12 days. By contrast, CEBPA N321D-transduced cells continued to grow for in excess of 56 days, with an initial cDC phenotype but by day 30 demonstrating a plasmacytoid dendritic cell precursor phenotype. CEBPA N321D-transduced cells were morphologically distinct from EV-transduced cells. To test leukemogenic potential in vivo, we performed transplantation experiments in lethally irradiated mice. Serial monitoring of peripheral blood demonstrated that Hoxb8-FL derived cells had disappeared by 4 weeks, and did not reappear. However, at 6 months CEBPA N321D-transduced cells could still be detected in bone marrow in contrast to EV-transduced cells but without any leukemic phenotype. To identify early events involved in pre-leukemic transformation, the differentiation profiles of EV, CEBPA WT and CEBPA N321D-transduced cells were examined with single cell RNA-seq (scRNA-seq). 576 single cells were taken from 3 biological replicates at days 0 and 5 post-differentiation, and analysed using the Automated Single-Cell Analysis Pipeline (Gardeux et al., Bioinformatics 2017). Visualisation by t-SNE (Fig 1) demonstrated: (i) CEBPA WT-transduced cells formed a distinct cluster at day 0 before withdrawal of estradiol; (ii) CEBPA N321D-transduced cells separated from EV and CEBPA WT-transduced cells after 5 days of differentiation, (iii) two subpopulations could be identified within the CEBPA N321D-transduced cells at day 5, with a cluster of five CEBPA N321D-transduced single cells distributed amongst or very close to the day 0 non-differentiated cells. Differential expression analysis identified 224 genes upregulated and 633 genes downregulated specifically in the CEBPA N321D-transduced cells when compared to EV cells after 5 days of differentiation. This gene expression signature revealed that CEBPA N321D-transduced cells switched on a HSC/MEP/CMP transcriptional program and switched off a myeloid dendritic cell program. Finally, in order to further dissect the effect of the N321D mutation, the binding profile of endogenous and CEBPA N321D was compared by ChIP-seq before and after 5 days of differentiation. Integration with scRNA-seq data identified 160 genes specifically downregulated in CEBPA N321D-transduced cells which were associated with the binding of the mutant protein. This list of genes included genes previously implicated in dendritic cell differentiation (such as NOTCH2, JAK2), as well as a number of genes not previously implicated in the evolution of AML, representing potentially novel therapeutic targets. Disclosures No relevant conflicts of interest to declare.

Ets transcription factors play important roles during the development and maintenance of the haematopoietic system. One such factor, Elf-1 (E74-like factor 1) controls the expression of multiple essential haematopoietic regulators including Scl/Tal1, Lmo2 and PU.1. However, to integrate Elf-1 into the wider regulatory hierarchies controlling haematopoietic development and differentiation, regulatory elements as well as upstream regulators of Elf-1 need to be identified. Here, we have used locus-wide comparative genomic analysis coupled with chromatin immunoprecipitation (ChIP-chip) assays which resulted in the identification of five distinct regulatory regions directing expression of Elf-1. Further, ChIP-chip assays followed by functional validation demonstrated that the key haematopoietic transcription factor PU.1 is a major upstream regulator of Elf-1. Finally, overexpression studies in a well-characterized erythroid differentiation assay from primary murine fetal liver cells demonstrated that Elf-1 downregulation is necessary for terminal erythroid differentiation. Given the known activation of PU.1 by Elf-1 and our newly identified reciprocal activation of Elf-1 by PU.1, identification of an inhibitory role for Elf-1 has significant implications for our understanding of how PU.1 controls myeloid–erythroid differentiation. Our findings therefore not only represent the first report of Elf-1 regulation but also enhance our understanding of the wider regulatory networks that control haematopoiesis.

The basic helix-loop-helix transcription factor Scl/Tal1 controls the development and subsequent differentiation of hematopoietic stem cells (HSCs). However, because few Scl target genes have been validated to date, the underlying mechanisms have remained largely unknown. In this study, we have used ChIP-Seq technology (coupling chromatin immunoprecipitation with deep sequencing) to generate a genome-wide catalog of Scl-binding events in a stem/progenitor cell line, followed by validation using primary fetal liver cells and comprehensive transgenic mouse assays. Transgenic analysis provided in vivo validation of multiple new direct Scl target genes and allowed us to reconstruct an in vivo validated network consisting of 17 factors and their respective regulatory elements. By coupling ChIP-Seq in model cell lines with in vivo transgenic validation and sophisticated bioinformatic analysis, we have identified a widely applicable strategy for the reconstruction of stem cell regulatory networks in which biologic material is otherwise limiting. Moreover, in addition to revealing multiple previously unrecognized links to known HSC regulators, as well as novel links to genes not previously implicated in HSC function, comprehensive transgenic analysis of regulatory elements provided substantial new insights into the transcriptional control of several important hematopoietic regulators, including Cbfa2t3h/Eto2, Cebpe, Nfe2, Zfpm1/Fog1, Erg, Mafk, Gfi1b, and Myb.

The human interleukin-3 (IL-3) and granulocyte-macrophage colony-stimulating-factor (GM-CSF, or CSF2) gene cluster arose by duplication of an ancestral gene. Although just 10 kb apart and responsive to the same signals, the IL-3 and GM-CSF genes are nevertheless regulated independently by separate, tissue-specific enhancers. To understand the differential regulation of the IL-3 and GM-CSF genes we have investigated a cluster of three ubiquitous DNase I-hypersensitive sites (DHSs) located between the two genes. We found that each site contains a conserved CTCF consensus sequence, binds CTCF, and recruits the cohesin subunit Rad21 in vivo. The positioning of these sites relative to the IL-3 and GM-CSF genes and their respective enhancers is conserved between human and mouse, suggesting a functional role in the organization of the locus. We found that these sites effectively block functional interactions between the GM-CSF enhancer and either the IL-3 or the GM-CSF promoter in reporter gene assays. These data argue that the regulation of the IL-3 and the GM-CSF promoters depends on the positions of their enhancers relative to the conserved CTCF/cohesin-binding sites. We suggest that one important role of these sites is to enable the independent regulation of the IL-3 and GM-CSF genes.

Transcription factor (TF) networks determine cell-type identity by establishing and maintaining lineage-specific expression profiles, yet reconstruction of mammalian regulatory network models has been hampered by a lack of comprehensive functional validation of regulatory interactions. Here, we report comprehensive ChIP-Seq, transgenic and reporter gene experimental data that have allowed us to construct an experimentally validated regulatory network model for haematopoietic stem/progenitor cells (HSPCs). Model simulation coupled with subsequent experimental validation using single cell expression profiling revealed potential mechanisms for cell state stabilisation, and also how a leukaemogenic TF fusion protein perturbs key HSPC regulators. The approach presented here should help to improve our understanding of both normal physiological and disease processes. DOI: http://dx.doi.org/10.7554/eLife.11469.001, eLife digest Blood stem cells and blood progenitor cells replenish a person’s entire blood system throughout their life and are crucial for survival. The stem cells have the potential to become any type of blood cell – including white blood cells and red blood cells – while the progenitor cells are slightly more restricted in the types of blood cell they can become. It is important to understand how the balance of cell types is maintained because, in cancers of the blood (also known as leukaemias), this organisation is lost and some cells proliferate abnormally. Almost all of a person’s cells will contain the same genetic information, but different cell types arise when different genes are switched on or off. The genes encoding proteins called transcription factors are particularly important because the proteins can control – either by activating or repressing – many other genes. Importantly, some of these genes will encode other transcription factors, meaning that these proteins essentially work together in networks. Schütte et al. have now combined extensive biochemical experiments with computational modelling to study some of the transcription factors that define blood stem cells and blood progenitor cells in mice. Firstly, nine transcription factors, which were already known to be important in blood stem cells, were thoroughly studied in mouse cells that could be grown in the laboratory. These experiments provided an overall view of which other genes these transcription factors control. Additional targeted investigations of the nine transcription factors then revealed how these proteins act in combination to activate or repress their respective activities. With this information, Schütte et al. built a computational model, which accurately reproduced how real mouse blood stem and progenitor cells behave when, for example, a transcription factor is deleted. Furthermore, the model could also predict what happens in single cells if the amounts of the transcription factors change. Lastly, Schütte et al. studied a common type of leukaemia. The model showed that the mutations that occur in this cancer change the finely tuned balance of the nine transcription factors; this may explain why leukaemia cells behave abnormally. In future these models could be extended to more transcription factors and other cell types and cancers. DOI: http://dx.doi.org/10.7554/eLife.11469.002

Comprehensive study of transcriptional control processes will be required to enhance our understanding of both normal and malignant hematopoiesis. Modern sequencing technologies have revolutionized our ability to generate genome-scale expression and histone modification profiles, transcription factor (TF)-binding maps, and also comprehensive chromatin-looping information. Many of these technologies, however, require large numbers of cells, and therefore cannot be applied to rare hematopoietic stem/progenitor cell (HSPC) populations. The stem cell factor-dependent multipotent progenitor cell line HPC-7 represents a well-recognized cell line model for HSPCs. Here we report genome-wide maps for 17 TFs, 3 histone modifications, DNase I hypersensitive sites, and high-resolution promoter-enhancer interactomes in HPC-7 cells. Integrated analysis of these complementary data sets revealed TF occupancy patterns of genomic regions involved in promoter-anchored loops. Moreover, preferential associations between pairs of TFs bound at either ends of chromatin loops led to the identification of 4 previously unrecognized protein-protein interactions between key blood stem cell regulators. All HPC-7 data sets are freely available both through standard repositories and a user-friendly Web interface. Together with previously generated genome-wide data sets, this study integrates HPC-7 data into a genomic resource on par with ENCODE tier 1 cell lines and, importantly, is the only current model with comprehensive genome-scale data that is relevant to HSPC biology.