Expression profiling of cerebrospinal fluid identifies dysregulated antiviral mechanisms in multiple sclerosis

Despite the overwhelming evidence that multiple sclerosis is an autoimmune disease, relatively little is known about the precise nature of the immune dysregulation underlying the development of the disease. Reasoning that the CSF from patients might be enriched for cells relevant in pathogenesis, we have completed a high-resolution single-cell analysis of 96 732 CSF cells collected from 33 patients with multiple sclerosis (n = 48 675) and 48 patients with other neurological diseases (n = 48 057). Completing comprehensive cell type annotation, we identified a rare population of CD8+ T cells, characterized by the upregulation of inhibitory receptors, increased in patients with multiple sclerosis. Applying a Multi-Omics Factor Analysis to these single-cell data further revealed that activity in pathways responsible for controlling inflammatory and type 1 interferon responses are altered in multiple sclerosis in both T cells and myeloid cells. We also undertook a systematic search for expression quantitative trait loci in the CSF cells. Of particular interest were two expression quantitative trait loci in CD8+ T cells that were fine mapped to multiple sclerosis susceptibility variants in the viral control genes ZC3HAV1 (rs10271373) and IFITM2 (rs1059091). Further analysis suggests that these associations likely reflect genetic effects on RNA splicing and cell-type specific gene expression respectively. Collectively, our study suggests that alterations in viral control mechanisms might be important in the development of multiple sclerosis.

Single-cell multi-omics analysis of COVID-19 patients with pre-existing autoimmune diseases shows aberrant immune responses to infection

In COVID-19, hyperinflammatory and dysregulated immune responses contribute to severity. Patients with pre-existing autoimmune conditions can therefore be at increased risk of severe COVID-19 and/or associated sequelae, yet SARS-CoV-2 infection in this group has been little studied. Here, we performed single-cell analysis of peripheral blood mononuclear cells from patients with three major autoimmune diseases (rheumatoid arthritis, psoriasis, or multiple sclerosis) during SARS-CoV-2 infection. We observed compositional differences between the autoimmune disease groups coupled with altered patterns of gene expression, transcription factor activity, and cell-cell communication that substantially shape the immune response under SARS-CoV-2 infection. While enrichment of HLA-DRlow CD14+ monocytes was observed in all three autoimmune disease groups, type-I interferon signaling as well as inflammatory T cell and monocyte responses varied widely between the three groups of patients. Our results reveal disturbed immune responses to SARS-CoV-2 in patients with pre-existing autoimmunity, highlighting important considerations for disease treatment and follow-up.

Preleukemic single-cell landscapes reveal mutation-specific mechanisms and gene programs predictive of AML patient outcomes

Acute myeloid leukemia (AML) and myeloid neoplasms develop through acquisition of somatic mutations that confer mutation-specific fitness advantages to hematopoietic stem and progenitor cells. However, our understanding of mutational effects remains limited to the resolution attainable within immunophenotypically and clinically accessible bulk cell populations. To decipher heterogeneous cellular fitness to preleukemic mutational perturbations, we performed single-cell RNA sequencing of eight different mouse models with driver mutations of myeloid malignancies, generating 269,048 single-cell profiles. Our analysis infers mutation-driven perturbations in cell abundance, cellular lineage fate, cellular metabolism, and gene expression at the continuous resolution, pinpointing cell populations with transcriptional alterations associated with differentiation bias. We further develop an 11-gene scoring system (Stem11) on the basis of preleukemic transcriptional signatures that predicts AML patient outcomes. Our results demonstrate that a single-cell-resolution deep characterization of preleukemic biology has the potential to enhance our understanding of AML heterogeneity and inform more effective risk stratification strategies.

Yolk sac cell atlas reveals multiorgan functions during human early development

The extraembryonic yolk sac (YS) ensures delivery of nutritional support and oxygen to the developing embryo but remains ill-defined in humans. We therefore assembled a comprehensive multiomic reference of the human YS from 3 to 8 postconception weeks by integrating single-cell protein and gene expression data. Beyond its recognized role as a site of hematopoiesis, we highlight roles in metabolism, coagulation, vascular development, and hematopoietic regulation. We reconstructed the emergence and decline of YS hematopoietic stem and progenitor cells from hemogenic endothelium and revealed a YS-specific accelerated route to macrophage production that seeds developing organs. The multiorgan functions of the YS are superseded as intraembryonic organs develop, effecting a multifaceted relay of vital functions as pregnancy proceeds.

Genotoxic aldehyde stress prematurely ages hematopoietic stem cells in a p53-driven manner

Aged hematopoietic stem cells (HSCs) display diminished self-renewal and a myeloid differentiation bias. However, the drivers and mechanisms that underpin this fundamental switch are not understood. HSCs produce genotoxic formaldehyde that requires protection by the detoxification enzymes ALDH2 and ADH5 and the Fanconi anemia (FA) DNA repair pathway. We find that the HSCs in young Aldh2-/-Fancd2-/- mice harbor a transcriptomic signature equivalent to aged wild-type HSCs, along with increased epigenetic age, telomere attrition, and myeloid-biased differentiation quantified by single HSC transplantation. In addition, the p53 response is vigorously activated in Aldh2-/-Fancd2-/- HSCs, while p53 deletion rescued this aged HSC phenotype. To further define the origins of the myeloid differentiation bias, we use a GFP genetic reporter to find a striking enrichment of Vwf+ myeloid and megakaryocyte-lineage-biased HSCs. These results indicate that metabolism-derived formaldehyde-DNA damage stimulates the p53 response in HSCs to drive accelerated aging.

STAT1 is essential for HSC function and maintains MHCIIhi stem cells that resist myeloablation and neoplastic expansion

Adult hematopoietic stem cells (HSCs) are predominantly quiescent and can be activated in response to acute stress such as infection or cytotoxic insults. STAT1 is a pivotal downstream mediator of interferon (IFN) signaling and is required for IFN-induced HSC proliferation, but little is known about the role of STAT1 in regulating homeostatic hematopoietic stem/progenitor cells (HSPCs). Here, we show that loss of STAT1 altered the steady state HSPC landscape, impaired HSC function in transplantation assays, delayed blood cell regeneration following myeloablation, and disrupted molecular programs that protect HSCs, including control of quiescence. Our results also reveal STAT1-dependent functional HSC heterogeneity. A previously unrecognized subset of homeostatic HSCs with elevated major histocompatibility complex class II (MHCII) expression (MHCIIhi) displayed molecular features of reduced cycling and apoptosis and was refractory to 5-fluorouracil-induced myeloablation. Conversely, MHCIIlo HSCs displayed increased megakaryocytic potential and were preferentially expanded in CALR mutant mice with thrombocytosis. Similar to mice, high MHCII expression is a feature of human HSCs residing in a deeper quiescent state. Our results therefore position STAT1 at the interface of stem cell heterogeneity and the interplay between stem cells and the adaptive immune system, areas of broad interest in the wider stem cell field.

Single-cell transcriptomics reveals the identity and regulators of human mast cell progenitors

Mast cell accumulation is a hallmark of a number of diseases, including allergic asthma and systemic mastocytosis. Immunoglobulin E-mediated crosslinking of the FcεRI receptors causes mast cell activation and contributes to disease pathogenesis. The mast cell lineage is one of the least studied among the hematopoietic cell lineages, and controversies remain about whether FcεRI expression appears during the mast cell progenitor stage or during terminal mast cell maturation. Here, we used single-cell transcriptomics analysis to reveal a temporal association between the appearance of FcεRI and the mast cell gene signature in CD34+ hematopoietic progenitors in adult peripheral blood. In agreement with these data, the FcεRI+ hematopoietic progenitors formed morphologically, phenotypically, and functionally mature mast cells in long-term culture assays. Single-cell transcriptomics analysis further revealed the expression patterns of prospective cytokine receptors regulating development of mast cell progenitors. Culture assays showed that interleukin-3 (IL-3) and IL-5 promoted disparate effects on progenitor cell proliferation and survival, respectively, whereas IL-33 caused robust FcεRI downregulation. Taken together, we showed that FcεRI expression appears at the progenitor stage of mast cell differentiation in peripheral blood. We also showed that external stimuli regulate FcεRI expression of mast cell progenitors, providing a possible explanation for the variable FcεRI expression levels during mast cell development.

Endothelial-specific Gata3 expression is required for hematopoietic stem cell generation

To generate sufficient numbers of transplantable hematopoietic stem cells (HSCs) in vitro, a detailed understanding of how this process takes place in vivo is essential. The endothelial-to-hematopoietic transition (EHT), which culminates in the production of the first HSCs, is a highly complex process during which key regulators are switched on and off at precise moments, and that is embedded into a myriad of microenvironmental signals from surrounding cells and tissues. We have previously demonstrated an HSC-supportive function for GATA3 within the sympathetic nervous system and the sub-aortic mesenchyme, but show here that it also plays a cell-intrinsic role during the EHT. It is expressed in hemogenic endothelial cells and early HSC precursors, where its expression correlates with a more quiescent state. Importantly, endothelial-specific deletion of Gata3 shows that it is functionally required for these cells to mature into HSCs, placing GATA3 at the core of the EHT regulatory network.

p57Kip2 regulates embryonic blood stem cells by controlling sympathoadrenal progenitor expansion

Hematopoietic stem cells (HSCs) are of major clinical importance, and finding methods for their in vitro generation is a prime research focus. We show here that the cell cycle inhibitor p57Kip2/Cdkn1c limits the number of emerging HSCs by restricting the size of the sympathetic nervous system (SNS) and the amount of HSC-supportive catecholamines secreted by these cells. This regulation occurs at the SNS progenitor level and is in contrast to the cell-intrinsic function of p57Kip2 in maintaining adult HSCs, highlighting profound differences in cell cycle requirements of adult HSCs compared with their embryonic counterparts. Furthermore, this effect is specific to the aorta-gonad-mesonephros (AGM) region and shows that the AGM is the main contributor to early fetal liver colonization, as early fetal liver HSC numbers are equally affected. Using a range of antagonists in vivo, we show a requirement for intact β2-adrenergic signaling for SNS-dependent HSC expansion. To gain further molecular insights, we have generated a single-cell RNA-sequencing data set of all Ngfr+ sympathoadrenal cells around the dorsal aorta to dissect their differentiation pathway. Importantly, this not only defined the relevant p57Kip2-expressing SNS progenitor stage but also revealed that some neural crest cells, upon arrival at the aorta, are able to take an alternative differentiation pathway, giving rise to a subset of ventrally restricted mesenchymal cells that express important HSC-supportive factors. Neural crest cells thus appear to contribute to the AGM HSC niche via 2 different mechanisms: SNS-mediated catecholamine secretion and HSC-supportive mesenchymal cell production.

Unique molecular and functional features of extramedullary hematopoietic stem and progenitor cell reservoirs in humans

Rare hematopoietic stem and progenitor cell (HSPC) pools outside the bone marrow (BM) contribute to blood production in stress and disease but remain ill-defined. Although nonmobilized peripheral blood (PB) is routinely sampled for clinical management, the diagnosis and monitoring potential of PB HSPCs remain untapped, as no healthy PB HSPC baseline has been reported. Here we comprehensively delineate human extramedullary HSPC compartments comparing spleen, PB, and mobilized PB to BM using single-cell RNA-sequencing and/or functional assays. We uncovered HSPC features shared by extramedullary tissues and others unique to PB. First, in contrast to actively dividing BM HSPCs, we found no evidence of substantial ongoing hematopoiesis in extramedullary tissues at steady state but report increased splenic HSPC proliferative output during stress erythropoiesis. Second, extramedullary hematopoietic stem cells/multipotent progenitors (HSCs/MPPs) from spleen, PB, and mobilized PB share a common transcriptional signature and increased abundance of lineage-primed subsets compared with BM. Third, healthy PB HSPCs display a unique bias toward erythroid-megakaryocytic differentiation. At the HSC/MPP level, this is functionally imparted by a subset of phenotypic CD71+ HSCs/MPPs, exclusively producing erythrocytes and megakaryocytes, highly abundant in PB but rare in other adult tissues. Finally, the unique erythroid-megakaryocytic-skewing of PB is perturbed with age in essential thrombocythemia and β-thalassemia. Collectively, we identify extramedullary lineage-primed HSPC reservoirs that are nonproliferative in situ and report involvement of splenic HSPCs during demand-adapted hematopoiesis. Our data also establish aberrant composition and function of circulating HSPCs as potential clinical indicators of BM dysfunction.

Coagulation factor V is a T-cell inhibitor expressed by leukocytes in COVID-19

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 in vitro. Anticoagulants that reduce FV conversion to FVa, including heparin, may have the unintended consequence of suppressing the adaptive immune system.

Simultaneous Analysis of Single-Cell Transcriptomes and Cell Surface Protein Expression of Human Hematopoietic Stem Cells and Progenitors Using the 10x Genomics Platform

The CITE-seq workflow combines conventional single-cell transcriptomic analysis with simultaneous analysis of cell surface protein expression using oligonucleotide-conjugated antibodies. This addition of immunophenotyping to mRNA data allows for a more detailed characterization of single-cell heterogeneity and can help to identify markers for the prospective isolation of transcriptionally defined novel cell subsets. Here, we describe the workflow for the preparation of human cord blood mononuclear cells and CD34+-enriched hematopoietic progenitors for the simultaneous characterization of protein and RNA using the commercially available TotalSeq™ antibodies from BioLegend and the droplet-based single-cell RNA-seq commercial platform from 10x Genomics.

Single-cell molecular profiling provides a high-resolution map of basophil and mast cell development

BACKGROUND

Basophils and mast cells contribute to the development of allergic reactions. Whereas these mature effector cells are extensively studied, the differentiation trajectories from hematopoietic progenitors to basophils and mast cells are largely uncharted at the single-cell level.

METHODS

We performed multicolor flow cytometry, high-coverage single-cell RNA sequencing analyses, and cell fate assays to chart basophil and mast cell differentiation at single-cell resolution in mouse.

RESULTS

Analysis of flow cytometry data reconstructed a detailed map of basophil and mast cell differentiation, including a bifurcation of progenitors into two specific trajectories. Molecular profiling and pseudotime ordering of the single cells revealed gene expression changes during differentiation. Cell fate assays showed that multicolor flow cytometry and transcriptional profiling successfully predict the bipotent phenotype of a previously uncharacterized population of peritoneal basophil-mast cell progenitors.

CONCLUSIONS

A combination of molecular and functional profiling of bone marrow and peritoneal cells provided a detailed road map of basophil and mast cell development. An interactive web resource was created to enable the wider research community to explore the expression dynamics for any gene of interest.

Single-cell multi-omics analysis of the immune response in COVID-19

Analysis of human blood immune cells provides insights into the coordinated response to viral infections such as severe acute respiratory syndrome coronavirus 2, which causes coronavirus disease 2019 (COVID-19). We performed single-cell transcriptome, surface proteome and T and B lymphocyte antigen receptor analyses of over 780,000 peripheral blood mononuclear cells from a cross-sectional cohort of 130 patients with varying severities of COVID-19. We identified expansion of nonclassical monocytes expressing complement transcripts (CD16+C1QA/B/C+) that sequester platelets and were predicted to replenish the alveolar macrophage pool in COVID-19. Early, uncommitted CD34+ hematopoietic stem/progenitor cells were primed toward megakaryopoiesis, accompanied by expanded megakaryocyte-committed progenitors and increased platelet activation. Clonally expanded CD8+ T cells and an increased ratio of CD8+ effector T cells to effector memory T cells characterized severe disease, while circulating follicular helper T cells accompanied mild disease. We observed a relative loss of IgA2 in symptomatic disease despite an overall expansion of plasmablasts and plasma cells. Our study highlights the coordinated immune response that contributes to COVID-19 pathogenesis and reveals discrete cellular components that can be targeted for therapy.

Single-Cell Analysis of Hematopoietic Stem Cells

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.

Two Aldehyde Clearance Systems Are Essential to Prevent Lethal Formaldehyde Accumulation in Mice and Humans

Reactive aldehydes arise as by-products of metabolism and are normally cleared by multiple families of enzymes. We find that mice lacking two aldehyde detoxifying enzymes, mitochondrial ALDH2 and cytoplasmic ADH5, have greatly shortened lifespans and develop leukemia. Hematopoiesis is disrupted profoundly, with a reduction of hematopoietic stem cells and common lymphoid progenitors causing a severely depleted acquired immune system. We show that formaldehyde is a common substrate of ALDH2 and ADH5 and establish methods to quantify elevated blood formaldehyde and formaldehyde-DNA adducts in tissues. Bone-marrow-derived progenitors actively engage DNA repair but also imprint a formaldehyde-driven mutation signature similar to aging-associated human cancer mutation signatures. Furthermore, we identify analogous genetic defects in children causing a previously uncharacterized inherited bone marrow failure and pre-leukemic syndrome. Endogenous formaldehyde clearance alone is therefore critical for hematopoiesis and in limiting mutagenesis in somatic tissues.

Interactions between lineage-associated transcription factors govern haematopoietic progenitor states

Recent advances in molecular profiling provide descriptive datasets of complex differentiation landscapes including the haematopoietic system, but the molecular mechanisms defining progenitor states and lineage choice remain ill-defined. Here, we employed a cellular model of murine multipotent haematopoietic progenitors (Hoxb8-FL) to knock out 39 transcription factors (TFs) followed by RNA-Seq analysis, to functionally define a regulatory network of 16,992 regulator/target gene links. Focussed analysis of the subnetworks regulated by the B-lymphoid TF Ebf1 and T-lymphoid TF Gata3 revealed a surprising role in common activation of an early myeloid programme. Moreover, Gata3-mediated repression of Pax5 emerges as a mechanism to prevent precocious B-lymphoid differentiation, while Hox-mediated activation of Meis1 suppresses myeloid differentiation. To aid interpretation of large transcriptomics datasets, we also report a new method that visualises likely transitions that a progenitor will undergo following regulatory network perturbations. Taken together, this study reveals how molecular network wiring helps to establish a multipotent progenitor state, with experimental approaches and analysis tools applicable to dissecting a broad range of both normal and perturbed cellular differentiation landscapes.

Interactions between lineage-associated transcription factors govern haematopoietic progenitor states

Recent advances in molecular profiling provide descriptive datasets of complex differentiation landscapes including the haematopoietic system, but the molecular mechanisms defining progenitor states and lineage choice remain ill-defined. Here, we employed a cellular model of murine multipotent haematopoietic progenitors (Hoxb8-FL) to knock out 39 transcription factors (TFs) followed by RNA-Seq analysis, to functionally define a regulatory network of 16,992 regulator/target gene links. Focussed analysis of the subnetworks regulated by the B-lymphoid TF Ebf1 and T-lymphoid TF Gata3 revealed a surprising role in common activation of an early myeloid programme. Moreover, Gata3-mediated repression of Pax5 emerges as a mechanism to prevent precocious B-lymphoid differentiation, while Hox-mediated activation of Meis1 suppresses myeloid differentiation. To aid interpretation of large transcriptomics datasets, we also report a new method that visualises likely transitions that a progenitor will undergo following regulatory network perturbations. Taken together, this study reveals how molecular network wiring helps to establish a multipotent progenitor state, with experimental approaches and analysis tools applicable to dissecting a broad range of both normal and perturbed cellular differentiation landscapes.

Manipulating niche composition limits damage to haematopoietic stem cells during Plasmodium infection

Severe infections are a major stress on haematopoiesis, where the consequences for haematopoietic stem cells (HSCs) have only recently started to emerge. HSC function critically depends on the integrity of complex bone marrow (BM) niches; however, what role the BM microenvironment plays in mediating the effects of infection on HSCs remains an open question. Here, using a murine model of malaria and combining single-cell RNA sequencing, mathematical modelling, transplantation assays and intravital microscopy, we show that haematopoiesis is reprogrammed upon infection, whereby the HSC compartment turns over substantially faster than at steady-state and HSC function is drastically affected. Interferon is found to affect both haematopoietic and mesenchymal BM cells and we specifically identify a dramatic loss of osteoblasts and alterations in endothelial cell function. Osteo-active parathyroid hormone treatment abolishes infection-triggered HSC proliferation and-coupled with reactive oxygen species quenching-enables partial rescuing of HSC function.

CHD7 and Runx1 interaction provides a braking mechanism for hematopoietic differentiation

Hematopoietic stem and progenitor cell (HSPC) formation and lineage differentiation involve gene expression programs orchestrated by transcription factors and epigenetic regulators. Genetic disruption of the chromatin remodeler chromodomain-helicase-DNA-binding protein 7 (CHD7) expanded phenotypic HSPCs, erythroid, and myeloid lineages in zebrafish and mouse embryos. CHD7 acts to suppress hematopoietic differentiation. Binding motifs for RUNX and other hematopoietic transcription factors are enriched at sites occupied by CHD7, and decreased RUNX1 occupancy correlated with loss of CHD7 localization. CHD7 physically interacts with RUNX1 and suppresses RUNX1-induced expansion of HSPCs during development through modulation of RUNX1 activity. Consequently, the RUNX1:CHD7 axis provides proper timing and function of HSPCs as they emerge during hematopoietic development or mature in adults, representing a distinct and evolutionarily conserved control mechanism to ensure accurate hematopoietic lineage differentiation.

Single-cell chromatin accessibility maps reveal regulatory programs driving early mouse organogenesis

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.

Gata3 targets Runx1 in the embryonic haematopoietic stem cell niche

Runx1 is an important haematopoietic transcription factor as stressed by its involvement in a number of haematological malignancies. Furthermore, it is a key regulator of the emergence of the first haematopoietic stem cells (HSCs) during development. The transcription factor Gata3 has also been linked to haematological disease and was shown to promote HSC production in the embryo by inducing the secretion of important niche factors. Both proteins are expressed in several different cell types within the aorta-gonads-mesonephros (AGM) region, in which the first HSCs are generated; however, a direct interaction between these two key transcription factors in the context of embryonic HSC production has not formally been demonstrated. In this current study, we have detected co-localisation of Runx1 and Gata3 in rare sub-aortic mesenchymal cells in the AGM. Furthermore, the expression of Runx1 is reduced in Gata3 -/- embryos, which also display a shift in HSC emergence. Using an AGM-derived cell line as a model for the stromal microenvironment in the AGM and performing ChIP-Seq and ChIP-on-chip experiments, we demonstrate that Runx1, together with other key niche factors, is a direct target gene of Gata3. In addition, we can pinpoint Gata3 binding to the Runx1 locus at specific enhancer elements which are active in the microenvironment. These results reveal a direct interaction between Gata3 and Runx1 in the niche that supports embryonic HSCs and highlight a dual role for Runx1 in driving the transdifferentiation of haemogenic endothelial cells into HSCs as well as in the stromal cells that support this process.

Cohesin-dependent regulation of gene expression during differentiation is lost in cohesin-mutated myeloid malignancies

Cohesin complex disruption alters gene expression, and cohesin mutations are common in myeloid neoplasia, suggesting a critical role in hematopoiesis. Here, we explore cohesin dynamics and regulation of hematopoietic stem cell homeostasis and differentiation. Cohesin binding increases at active regulatory elements only during erythroid differentiation. Prior binding of the repressive Ets transcription factor Etv6 predicts cohesin binding at these elements and Etv6 interacts with cohesin at chromatin. Depletion of cohesin severely impairs erythroid differentiation, particularly at Etv6-prebound loci, but augments self-renewal programs. Together with corroborative findings in acute myeloid leukemia and myelodysplastic syndrome patient samples, these data suggest cohesin-mediated alleviation of Etv6 repression is required for dynamic expression at critical erythroid genes during differentiation and how this may be perturbed in myeloid malignancies.

Single-cell approaches identify the molecular network driving malignant hematopoietic stem cell self-renewal

Recent advances in single-cell technologies have permitted the investigation of heterogeneous cell populations at previously unattainable resolution. Here we apply such approaches to resolve the molecular mechanisms driving disease in mouse hematopoietic stem cells (HSCs), using JAK2V617F mutant myeloproliferative neoplasms (MPNs) as a model. Single-cell gene expression and functional assays identified a subset of JAK2V617F mutant HSCs that display defective self-renewal. This defect is rescued at the single HSC level by crossing JAK2V617F mice with mice lacking TET2, the most commonly comutated gene in patients with MPN. Single-cell gene expression profiling of JAK2V617F-mutant HSCs revealed a loss of specific regulator genes, some of which were restored to normal levels in single TET2/JAK2 mutant HSCs. Of these, Bmi1 and, to a lesser extent, Pbx1 and Meis1 overexpression in JAK2-mutant HSCs could drive a disease phenotype and retain durable stem cell self-renewal in functional assays. Together, these single-cell approaches refine the molecules involved in clonal expansion of MPNs and have broad implications for deconstructing the molecular network of normal and malignant stem cells.

A single-cell hematopoietic landscape resolves 8 lineage trajectories and defects in Kit mutant mice

Hematopoietic stem and progenitor cells (HSPCs) maintain the adult blood system, and their dysregulation causes a multitude of diseases. However, the differentiation journeys toward specific hematopoietic lineages remain ill defined, and system-wide disease interpretation remains challenging. Here, we have profiled 44 802 mouse bone marrow HSPCs using single-cell RNA sequencing to provide a comprehensive transcriptional landscape with entry points to 8 different blood lineages (lymphoid, megakaryocyte, erythroid, neutrophil, monocyte, eosinophil, mast cell, and basophil progenitors). We identified a common basophil/mast cell bone marrow progenitor and characterized its molecular profile at the single-cell level. Transcriptional profiling of 13 815 HSPCs from the c-Kit mutant (W41/W41) mouse model revealed the absence of a distinct mast cell lineage entry point, together with global shifts in cell type abundance. Proliferative defects were accompanied by reduced Myc expression. Potential compensatory processes included upregulation of the integrated stress response pathway and downregulation of proapoptotic gene expression in erythroid progenitors, thus providing a template of how large-scale single-cell transcriptomic studies can bridge between molecular phenotypes and quantitative population changes.

Reconstructing blood stem cell regulatory network models from single-cell molecular profiles

Adult blood contains a mixture of mature cell types, each with specialized functions. Single hematopoietic stem cells (HSCs) have been functionally shown to generate all mature cell types for the lifetime of the organism. Differentiation of HSCs toward alternative lineages must be balanced at the population level by the fate decisions made by individual cells. Transcription factors play a key role in regulating these decisions and operate within organized regulatory programs that can be modeled as transcriptional regulatory networks. As dysregulation of single HSC fate decisions is linked to fatal malignancies such as leukemia, it is important to understand how these decisions are controlled on a cell-by-cell basis. Here we developed and applied a network inference method, exploiting the ability to infer dynamic information from single-cell snapshot expression data based on expression profiles of 48 genes in 2,167 blood stem and progenitor cells. This approach allowed us to infer transcriptional regulatory network models that recapitulated differentiation of HSCs into progenitor cell types, focusing on trajectories toward megakaryocyte-erythrocyte progenitors and lymphoid-primed multipotent progenitors. By comparing these two models, we identified and subsequently experimentally validated a difference in the regulation of nuclear factor, erythroid 2 (Nfe2) and core-binding factor, runt domain, alpha subunit 2, translocated to, 3 homolog (Cbfa2t3h) by the transcription factor Gata2. Our approach confirms known aspects of hematopoiesis, provides hypotheses about regulation of HSC differentiation, and is widely applicable to other hierarchical biological systems to uncover regulatory relationships.

Resolving early mesoderm diversification through single-cell expression profiling

In mammals, specification of the three major germ layers occurs during gastrulation, when cells ingressing through the primitive streak differentiate into the precursor cells of major organ systems. However, the molecular mechanisms underlying this process remain unclear, as numbers of gastrulating cells are very limited. In the mouse embryo at embryonic day 6.5, cells located at the junction between the extra-embryonic region and the epiblast on the posterior side of the embryo undergo an epithelial-to-mesenchymal transition and ingress through the primitive streak. Subsequently, cells migrate, either surrounding the prospective ectoderm contributing to the embryo proper, or into the extra-embryonic region to form the yolk sac, umbilical cord and placenta. Fate mapping has shown that mature tissues such as blood and heart originate from specific regions of the pre-gastrula epiblast, but the plasticity of cells within the embryo and the function of key cell-type-specific transcription factors remain unclear. Here we analyse 1,205 cells from the epiblast and nascent Flk1(+) mesoderm of gastrulating mouse embryos using single-cell RNA sequencing, representing the first transcriptome-wide in vivo view of early mesoderm formation during mammalian gastrulation. Additionally, using knockout mice, we study the function of Tal1, a key haematopoietic transcription factor, and demonstrate, contrary to previous studies performed using retrospective assays, that Tal1 knockout does not immediately bias precursor cells towards a cardiac fate.

Advancing haematopoietic stem and progenitor cell biology through single-cell profiling

Haematopoietic stem and progenitor cells (HSPCs) sit at the top of the haematopoietic hierarchy, and their fate choices need to be carefully controlled to ensure balanced production of all mature blood cell types. As cell fate decisions are made at the level of the individual cells, recent technological advances in measuring gene and protein expression in increasingly large numbers of single cells have been rapidly adopted to study both normal and pathological HSPC function. In this review we emphasise the importance of combining the correct computational models with single-cell experimental techniques, and illustrate how such integrated approaches have been used to resolve heterogeneities in populations, reconstruct lineage differentiation, identify regulatory relationships and link molecular profiling to cellular function.

An experimentally validated network of nine haematopoietic transcription factors reveals mechanisms of cell state stability

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.

Index sorting resolves heterogeneous murine hematopoietic stem cell populations

Recent advances in the cellular and molecular biology of single stem cells have uncovered significant heterogeneity in the functional properties of stem cell populations. This has prompted the development of approaches to study single cells in isolation, often performed using multiparameter flow cytometry. However, many stem cell populations are too rare to test all possible cell surface marker combinations, and virtually nothing is known about functional differences associated with varying intensities of such markers. Here we describe the use of index sorting for further resolution of the flow cytometric isolation of single murine hematopoietic stem cells (HSCs). Specifically, we associate single-cell functional assay outcomes with distinct cell surface marker expression intensities. High levels of both CD150 and EPCR associate with delayed kinetics of cell division and low levels of differentiation. Moreover, cells that do not form single HSC-derived clones appear in the 7AAD(dim) fraction, suggesting that even low levels of 7AAD staining are indicative of less healthy cell populations. These data indicate that when used in combination with single-cell functional assays, index sorting is a powerful tool for refining cell isolation strategies. This approach can be broadly applied to other single-cell systems, both to improve isolation and to acquire additional cell surface marker information.

Combined Single-Cell Functional and Gene Expression Analysis Resolves Heterogeneity within Stem Cell Populations

Heterogeneity within the self-renewal durability of adult hematopoietic stem cells (HSCs) challenges our understanding of the molecular framework underlying HSC function. Gene expression studies have been hampered by the presence of multiple HSC subtypes and contaminating non-HSCs in bulk HSC populations. To gain deeper insight into the gene expression program of murine HSCs, we combined single-cell functional assays with flow cytometric index sorting and single-cell gene expression assays. Through bioinformatic integration of these datasets, we designed an unbiased sorting strategy that separates non-HSCs away from HSCs, and single-cell transplantation experiments using the enriched population were combined with RNA-seq data to identify key molecules that associate with long-term durable self-renewal, producing a single-cell molecular dataset that is linked to functional stem cell activity. Finally, we demonstrated the broader applicability of this approach for linking key molecules with defined cellular functions in another stem cell system.

CODEX: a next-generation sequencing experiment database for the haematopoietic and embryonic stem cell communities

CODEX (http://codex.stemcells.cam.ac.uk/) is a user-friendly database for the direct access and interrogation of publicly available next-generation sequencing (NGS) data, specifically aimed at experimental biologists. In an era of multi-centre genomic dataset generation, CODEX provides a single database where these samples are collected, uniformly processed and vetted. The main drive of CODEX is to provide the wider scientific community with instant access to high-quality NGS data, which, irrespective of the publishing laboratory, is directly comparable. CODEX allows users to immediately visualize or download processed datasets, or compare user-generated data against the database's cumulative knowledge-base. CODEX contains four types of NGS experiments: transcription factor chromatin immunoprecipitation coupled to high-throughput sequencing (ChIP-Seq), histone modification ChIP-Seq, DNase-Seq and RNA-Seq. These are largely encompassed within two specialized repositories, HAEMCODE and ESCODE, which are focused on haematopoiesis and embryonic stem cell samples, respectively. To date, CODEX contains over 1000 samples, including 221 unique TFs and 93 unique cell types. CODEX therefore provides one of the most complete resources of publicly available NGS data for the direct interrogation of transcriptional programmes that regulate cellular identity and fate in the context of mammalian development, homeostasis and disease.

Key regulators control distinct transcriptional programmes in blood progenitor and mast cells

Despite major advances in the generation of genome-wide binding maps, the mechanisms by which transcription factors (TFs) regulate cell type identity have remained largely obscure. Through comparative analysis of 10 key haematopoietic TFs in both mast cells and blood progenitors, we demonstrate that the largely cell type-specific binding profiles are not opportunistic, but instead contribute to cell type-specific transcriptional control, because (i) mathematical modelling of differential binding of shared TFs can explain differential gene expression, (ii) consensus binding sites are important for cell type-specific binding and (iii) knock-down of blood stem cell regulators in mast cells reveals mast cell-specific genes as direct targets. Finally, we show that the known mast cell regulators Mitf and c-fos likely contribute to the global reorganisation of TF binding profiles. Taken together therefore, our study elucidates how key regulatory TFs contribute to transcriptional programmes in several distinct mammalian cell types.

GFI1 and GFI1B control the loss of endothelial identity of hemogenic endothelium during hematopoietic commitment

Recent studies have established that during embryonic development, hematopoietic progenitors and stem cells are generated from hemogenic endothelium precursors through a process termed endothelial to hematopoietic transition (EHT). The transcription factor RUNX1 is essential for this process, but its main downstream effectors remain largely unknown. Here, we report the identification of Gfi1 and Gfi1b as direct targets of RUNX1 and critical regulators of EHT. GFI1 and GFI1B are able to trigger, in the absence of RUNX1, the down-regulation of endothelial markers and the formation of round cells, a morphologic change characteristic of EHT. Conversely, blood progenitors in Gfi1- and Gfi1b-deficient embryos maintain the expression of endothelial genes. Moreover, those cells are not released from the yolk sac and disseminated into embryonic tissues. Taken together, our findings demonstrate a critical and specific role of the GFI1 transcription factors in the first steps of the process leading to the generation of hematopoietic progenitors from hemogenic endothelium.

The transcriptional programme controlled by Runx1 during early embryonic blood development

Transcription factors have long been recognised as powerful regulators of mammalian development yet it is largely unknown how individual key regulators operate within wider regulatory networks. Here we have used a combination of global gene expression and chromatin-immunoprecipitation approaches during the early stages of haematopoietic development to define the transcriptional programme controlled by Runx1, an essential regulator of blood cell specification. Integrated analysis of these complementary genome-wide datasets allowed us to construct a global regulatory network model, which suggested that key regulators are activated sequentially during blood specification, but will ultimately collaborate to control many haematopoietically expressed genes. Using the CD41/integrin alpha 2b gene as a model, cellular and in vivo studies showed that CD41 is controlled by both Scl/Tal1 and Runx1 in fully specified blood cells, and initiation of CD41 expression in E7.5 embryos is severely compromised in the absence of Runx1. Taken together, this study represents the first global analysis of the transcriptional programme controlled by any key haematopoietic regulator during the process of early blood cell specification. Moreover, the concept of interplay between sequentially deployed core regulators is likely to represent a design principle widely applicable to the transcriptional control of mammalian development.

Long-range DNA looping and gene expression analyses identify DEXI as an autoimmune disease candidate gene

The chromosome 16p13 region has been associated with several autoimmune diseases, including type 1 diabetes (T1D) and multiple sclerosis (MS). CLEC16A has been reported as the most likely candidate gene in the region, since it contains the most disease-associated single-nucleotide polymorphisms (SNPs), as well as an imunoreceptor tyrosine-based activation motif. However, here we report that intron 19 of CLEC16A, containing the most autoimmune disease-associated SNPs, appears to behave as a regulatory sequence, affecting the expression of a neighbouring gene, DEXI. The CLEC16A alleles that are protective from T1D and MS are associated with increased expression of DEXI, and no other genes in the region, in two independent monocyte gene expression data sets. Critically, using chromosome conformation capture (3C), we identified physical proximity between the DEXI promoter region and intron 19 of CLEC16A, separated by a loop of >150 kb. In reciprocal experiments, a 20 kb fragment of intron 19 of CLEC16A, containing SNPs associated with T1D and MS, as well as with DEXI expression, interacted with the promotor region of DEXI but not with candidate DNA fragments containing other potential causal genes in the region, including CLEC16A. Intron 19 of CLEC16A is highly enriched for transcription-factor-binding events and markers associated with enhancer activity. Taken together, these data indicate that although the causal variants in the 16p13 region lie within CLEC16A, DEXI is an unappreciated autoimmune disease candidate gene, and illustrate the power of the 3C approach in progressing from genome-wide association studies results to candidate causal genes.

RUNX1 regulates the CD34 gene in haematopoietic stem cells by mediating interactions with a distal regulatory element

The transcription factor RUNX1 is essential to establish the haematopoietic gene expression programme; however, the mechanism of how it activates transcription of haematopoietic stem cell (HSC) genes is still elusive. Here, we obtained novel insights into RUNX1 function by studying regulation of the human CD34 gene, which is expressed in HSCs. Using transgenic mice carrying human CD34 PAC constructs, we identified a novel downstream regulatory element (DRE), which is bound by RUNX1 and is necessary for human CD34 expression in long-term (LT)-HSCs. Conditional deletion of Runx1 in mice harbouring human CD34 promoter-DRE constructs abrogates human CD34 expression. We demonstrate by chromosome conformation capture assays in LT-HSCs that the DRE physically interacts with the human CD34 promoter. Targeted mutagenesis of RUNX binding sites leads to perturbation of this interaction and decreased human CD34 expression in LT-HSCs. Overall, our in vivo data provide novel evidence about the role of RUNX1 in mediating interactions between distal and proximal elements of the HSC gene CD34.

Deciphering transcriptional control mechanisms in hematopoiesis:the impact of high-throughput sequencing technologies

One of the key challenges facing biomedical research is to extract biologically meaningful information from the ever-increasing scale and complexity of datasets generated through high-throughput approaches. Hematopoiesis represents one of the most experimentally tractable mammalian organ systems and, therefore, has historically tended to be at the forefront of applying new technologies within biomedical research. The combination of massive parallel sequencing technologies with chromatin-immunoprecipitation (ChIP-Seq) permits genome-scale characterization of histone modification status and identification of the complete set of binding sites for transcription factors. Because transcription factors have long been recognized as essential regulators of cell fate choice in hematopoiesis, ChIP-Seq technology has rapidly entered the arena of modern experimental hematology. Here we review the biological insights gained from ChIP-Seq studies performed in the hematopoietic system since the earliest studies just 4 years ago. A surprisingly large number of different approaches have already been implemented to extract new biological knowledge from ChIP-Seq datasets. By focusing on successful insights from multiple different approaches, we hope to provide stimulating reading for anyone wanting to utilize ChIP-Seq technology within their particular research field.

Cell of origin in AML: susceptibility to MN1-induced transformation is regulated by the MEIS1/AbdB-like HOX protein complex

Pathways defining susceptibility of normal cells to oncogenic transformation may be valuable therapeutic targets. We characterized the cell of origin and its critical pathways in MN1-induced leukemias. Common myeloid (CMP) but not granulocyte-macrophage progenitors (GMP) could be transformed by MN1. Complementation studies of CMP-signature genes in GMPs demonstrated that MN1-leukemogenicity required the MEIS1/AbdB-like HOX-protein complex. ChIP-sequencing identified common target genes of MN1 and MEIS1 and demonstrated identical binding sites for a large proportion of their chromatin targets. Transcriptional repression of MEIS1 targets in established MN1 leukemias demonstrated antileukemic activity. As MN1 relies on but cannot activate expression of MEIS1/AbdB-like HOX proteins, transcriptional activity of these genes determines cellular susceptibility to MN1-induced transformation and may represent a promising therapeutic target.

A compendium of genome-wide hematopoietic transcription factor maps supports the identification of gene regulatory control mechanisms

OBJECTIVE

Key regulators of blood stem cell differentiation into the various mature hematopoietic lineages are commonly encoded by transcription factor genes. Elucidation of transcriptional regulatory mechanisms therefore holds great promise in advancing our understanding of both normal and malignant hematopoiesis. Recent technological advances have enabled the generation of genome-wide transcription factor binding maps using chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq). However, transcription factors operate in a combinatorial fashion suggesting that integrated analysis of genome-wide maps for multiple transcription factors will be essential to fully exploit these new genome-scale data sets.

MATERIALS AND METHODS

Here we have generated a compendium that integrates 53 ChIP-Seq studies covering 30 factors across all major hematopoietic lineages with a total of 754,380 binding peaks. We also used transgenic mouse assays to validate a newly predicted transcriptional enhancer.

RESULTS

Integrated analysis of all 53 ChIP-Seq studies demonstrated that cell-type identity exerts a larger influence on global transcription factor binding patterns than the nature of the individual transcription factors. Furthermore, regions highlighted by multifactor binding within specific gene loci overlap with known regulatory elements and also provide a useful guide for identifying novel elements, as demonstrated by transgenic analysis of a previously unrecognized enhancer in the Maml3 gene locus.

CONCLUSIONS

The ChIP-Seq compendium described here provides a valuable resource for the wider research community by accelerating the discovery of transcriptional mechanisms operating in the hematopoietic system.

Combinatorial transcriptional control in blood stem/progenitor cells: genome-wide analysis of ten major transcriptional regulators

Combinatorial transcription factor (TF) interactions control cellular phenotypes and, therefore, underpin stem cell formation, maintenance, and differentiation. Here, we report the genome-wide binding patterns and combinatorial interactions for ten key regulators of blood stem/progenitor cells (SCL/TAL1, LYL1, LMO2, GATA2, RUNX1, MEIS1, PU.1, ERG, FLI-1, and GFI1B), thus providing the most comprehensive TF data set for any adult stem/progenitor cell type to date. Genome-wide computational analysis of complex binding patterns, followed by functional validation, revealed the following: first, a previously unrecognized combinatorial interaction between a heptad of TFs (SCL, LYL1, LMO2, GATA2, RUNX1, ERG, and FLI-1). Second, we implicate direct protein-protein interactions between four key regulators (RUNX1, GATA2, SCL, and ERG) in stabilizing complex binding to DNA. Third, Runx1(+/-)::Gata2(+/-) compound heterozygous mice are not viable with severe hematopoietic defects at midgestation. Taken together, this study demonstrates the power of genome-wide analysis in generating novel functional insights into the transcriptional control of stem and progenitor cells.

Transcriptional regulation of haematopoietic transcription factors

The control of differential gene expression is central to all metazoan biology. Haematopoiesis represents one of the best understood developmental systems where multipotent blood stem cells give rise to a range of phenotypically distinct mature cell types, all characterised by their own distinctive gene expression profiles. Small combinations of lineage-determining transcription factors drive the development of specific mature lineages from multipotent precursors. Given their powerful regulatory nature, it is imperative that the expression of these lineage-determining transcription factors is under tight control, a fact underlined by the observation that their misexpression commonly leads to the development of leukaemia. Here we review recent studies on the transcriptional control of key haematopoietic transcription factors, which demonstrate that gene loci contain multiple modular regulatory regions within which specific regulatory codes can be identified, that some modular elements cooperate to mediate appropriate tissue-specific expression, and that long-range approaches will be necessary to capture all relevant regulatory elements. We also explore how changes in technology will impact on this area of research in the future.

Transcriptional regulation of Elf-1: locus-wide analysis reveals four distinct promoters, a tissue-specific enhancer, control by PU.1 and the importance of Elf-1 downregulation for erythroid maturation

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.

From genes to cells to tissues--modelling the haematopoietic system

Haematopoiesis (or blood formation) in general and haematopoietic stem cells more specifically represent some of the best studied mammalian developmental systems. Sophisticated purification protocols coupled with powerful biological assays permit functional analysis of highly purified cell populations both in vitro and in vivo. However, despite several decades of intensive research, the sheer complexity of the haematopoietic system means that many important questions remain unanswered or even unanswerable with current experimental tools. Scientists have therefore increasingly turned to modelling to tackle complexity at multiple levels ranging from networks of genes to the behaviour of cells and tissues. Early modelling attempts of gene regulatory networks have focused on core regulatory circuits but have more recently been extended to genome-wide datasets such as expression profiling and ChIP-sequencing data. Modelling of haematopoietic cells and tissues has provided insight into the importance of phenotypic heterogeneity for the differentiation of normal progenitor cells as well as a greater understanding of treatment response for particular pathologies such as chronic myeloid leukaemia. Here we will review recent progress in attempts to reconstruct segments of the haematopoietic system. A variety of modelling strategies will be covered from small-scale, protein-DNA or protein-protein interactions to large scale reconstructions. Also discussed will be examples of how stochastic modelling may be applied to multi cell systems such as those seen in normal and malignant haematopoiesis.

Levels of persistent organic pollutants in several child day care centers

The concentrations of a suite of persistent organic chemicals were measured in multiple media in 10 child day care centers located in central North Carolina. Five centers served mainly children from low-income families, as defined by the federal Women, Infants, and Children (WIC) assistance program, and five served mainly children from middle-income families. The targeted chemicals were chosen because of their probable carcinogenicity, acute or chronic toxicity, or hypothesized potential for endocrine system disruption. Targeted compounds included polycyclic aromatic hydrocarbons (PAHs), pentachloro- and nonyl-phenol, bisphenol-A, dibutyl and butylbenzyl phthalate, polychlorinated biphenyls (PCBs), organochlorine pesticides, the organophosphate pesticides diazinon and chlorpyrifos, and the herbicide 2,4-dichlorophenoxyacetic acid (2,4D). Sampled media were indoor and outdoor air, food and beverages, indoor dust, and outdoor play area soil. Concentrations of the targeted compounds were determined using a combination of extraction and analysis methods, depending on the media. Analysis was predominantly by gas chromatography/mass spectrometry (GC/MS) or gas chromatography with electron capture detection (GC/ECD). Concentrations of the targeted pollutants were low and well below the levels generally considered to be of concern as possible health hazards. Potential exposures to the target compounds were estimated from the concentrations in the various media, the children's daily time-activity schedules at day care, and the best currently available estimates of the inhalation rates (8.3 m(3)/day) and soil ingestion rates (100 mg/day) of children ages 3-5. The potential exposures for the target compounds differed depending on the compound class and the sampled media. Potential exposures through dietary ingestion were greater than those through inhalation, which were greater than those through nondietary ingestion, for the total of all PAHs, the phenols, the organophosphate pesticides, and the organochlorine pesticides. Potential exposures through dietary ingestion were greater than those through nondietary ingestion, which were greater than those through inhalation, for those PAHs that are probable human carcinogens (B2 PAH), the phthalate esters, and 2,4D. For the PCBs, exposures through inhalation were greater than those through nondietary ingestion, and exposures through dietary ingestion were smallest. Differences in targeted compound levels between the centers that serve mainly low-income clients and those that serve mainly middle-income clients were small and depended on the compound class and the medium.

Polycyclic aromatic hydrocarbon exposures of children in low-income families

Children in low-income families may have high exposures to polycyclic aromatic hydrocarbons (PAH). Such exposures could result from household proximity to heavy traffic or industrial sources, environmental tobacco smoke, contaminated house dust or soil, among others. The objectives of this study were: to establish methods for measuring total PAH exposure of children in low-income families, to estimate the PAH exposures of these children, and to estimate the relative importance of the environmental pathways for PAH exposure. Analytical methods to determine PAH in air, dust, soil, and food and to determine hydroxy-PAH in urine samples were evaluated and validated. A two-home pilot study was conducted in downtown Durham, North Carolina (NC) during February 1994. One smoker's and one nonsmoker's household, which had preschool children and income at or below the official U.S. poverty level, participated. A nine-home winter and a nine-home summer study were conducted in Durham and the NC Piedmont area during February 1995 and August 1995, respectively. A summer study in four smokers' homes was also conducted. In each of these studies, multimedia samples were collected and analyzed for PAH or hydroxy-PAH. Summary statistics, Pearson correlations, and analysis of variance were performed on the combined data from these four field studies. An effective screening method was established for recruiting low-income families. The field protocol involved measurements of three homes in 2-day periods. This protocol should be suitable for large-scale studies. The results showed that indoor PAH levels were generally higher than outdoor PAH levels. Higher indoor PAH levels were observed in the smokers' homes compared to nonsmokers' homes. Higher outdoor PAH levels were found in inner city as opposed to rural areas. The relative concentration trend for PAH in dust and soil was: house dust > entryway dust > pathway soil. The PAH concentrations in adults' food samples were generally higher than those in children's food samples. Children's potential daily doses of PAH were higher than those of adults in the same household, when intakes were normalized to body weights. Inhalation is an important pathway for children's exposure to total PAH because of the high levels of naphthalene present in both indoor and outdoor air. Dietary ingestion and nondietary ingestion pathways became more important for children's exposure to the B2 PAH (ranked as probable human carcinogens, B2 by the U.S. EPA's Integrated Risk System), most of which are of low volatility. The analysis of variance results showed that inner city participants had higher total exposure to B2 PAH than did rural participants.

Evaluation of enzyme-linked immunosorbent assay for the determination of polycyclic aromatic hydrocarbons in house dust and residential soil

Two commercially available enzyme-linked immunosorbent assays (ELISA) for total polycyclic aromatic hydrocarbon (PAH) and carcinogenic PAH (C-PAH) were evaluated. The testing procedures were refined for application to screening PAH and C-PAH in house dust and soil samples for human exposure studies. The overall method precision expressed as percent relative standard deviation (%RSD) of triplicate real world dust and soil samples was within +/- 29% (12-29%) for PAH ELISA and +/- 21% (5.9-21%) for C-PAH ELISA. Spike recoveries from real world dust/soil samples were 114 +/- 30% for phenanthrene from PAH ELISA and 120 +/- 8.2% for benzo[a]pyrene from C-PAH ELISA. The overall method accuracy for PAH and C-PAH assays cannot be assessed for multiple PAH components in dust/soil samples (which represent real-world samples), because of the assays' cross reactivities with other PAH components. Over 100 dust/soil samples from 13 North Carolina homes and 22 Arizona homes were analyzed by PAH and C-PAH assays, as well as by the conventional gas chromatography/mass spectrometry (GC/MS) method. Statistical analysis showed that dust/soil PAH data from ELISA and GC/MS methods are significantly different. In general PAH ELISA responses were higher than PAH GC/MS responses. The regression analysis showed that the linear relationship between ELISA and GC/MS measurements is not strong in the combined data. The relationship became stronger for the data from the same type of dust/soil samples. The screening performance of ELISA was evaluated based on the frequency distribution of ELISA and GC/MS data. The results indicated that the ELISA PAH and C-PAH assays cannot be used as a quantitative analytical tool for determining PAH in real-world dust/soil samples. However, the ELISA is an effective screening tool for ranking PAH concentrations in similar types of real world dust/soil samples.