Key Points
A new Gata2 reporter indicates that all HSCs express Gata2 and corroborates findings that Gata2 is not required for generation of all HPCs.
Isolatable non–Gata2-expressing HPCs show less potency and a distinct genetic program, thus having implications for reprogramming strategies.
Abstract
The Gata2 transcription factor is a pivotal regulator of hematopoietic cell development and maintenance, highlighted by the fact that Gata2 haploinsufficiency has been identified as the cause of some familial cases of acute myelogenous leukemia/myelodysplastic syndrome and in MonoMac syndrome. Genetic deletion in mice has shown that Gata2 is pivotal to the embryonic generation of hematopoietic stem cells (HSCs) and hematopoietic progenitor cells (HPCs). It functions in the embryo during endothelial cell to hematopoietic cell transition to affect hematopoietic cluster, HPC, and HSC formation. Gata2 conditional deletion and overexpression studies show the importance of Gata2 levels in hematopoiesis, during all developmental stages. Although previous studies of cell populations phenotypically enriched in HPCs and HSCs show expression of Gata2, there has been no direct study of Gata2 expressing cells during normal hematopoiesis. In this study, we generate a Gata2Venus reporter mouse model with unperturbed Gata2 expression to examine the hematopoietic function and transcriptome of Gata2 expressing and nonexpressing cells. We show that all the HSCs are Gata2 expressing. However, not all HPCs in the aorta, vitelline and umbilical arteries, and fetal liver require or express Gata2. These Gata2-independent HPCs exhibit a different functional output and genetic program, including Ras and cyclic AMP response element-binding protein pathways and other Gata factors, compared with Gata2-dependent HPCs. Our results, indicating that Gata2 is of major importance in programming toward HSC fate but not in all cells with HPC fate, have implications for current reprogramming strategies.
Introduction
Gata2 is one of the “heptad” transcription factors that acts on regulatory regions of hematopoietic genes.1 It is upregulated in vivo in Ly6aGFP+ cells undergoing endothelial-to-hematopoietic cell transition (EHT), a process by which definitive hematopoietic progenitors (HPCs) and hematopoietic stem cells (HSCs) are generated in the embryo.2,3 As one of the major regulators of HPC and HSC generation, germline deficiency of Gata2 results in embryonic lethality between embryonic day (E)10 and E10.5 and an anemic phenotype, with a decreased number of primitive and definitive HPCs in the yolk sac (YS) and in Gata2−/− embryonic stem (ES) cell hematopoietic differentiation cultures.4-6 Chimeric embryo generation with Gata2−/− ES cells revealed defective production of all hematopoietic lineages.5 The E10.5 lethality of Gata2−/− embryos precludes the study of HSC generation in the aorta-gonad-mesonephros (AGM) region, the first site of de novo HSC production. Gata2+/− embryos contain greatly reduced number of HSCs in the AGM region.7,8 Gata2 haploinsufficiency perturbs adult HSC homeostasis in mice9 and, in humans, leads to MonoMac syndrome,10 which is associated with sporadic myelodysplasia and myeloid leukemia. Also, rearrangement of the remote Gata2 enhancer drives acute myeloid leukemogenesis by activating Evi1 expression.11,12 Overexpression studies also reveal that levels of Gata2 expression are important for its hematopoietic function.13-15 In situ hybridization studies localize Gata2 expression to aortic endothelial cells, intra-aortic hematopoietic cluster cells, placenta (PL), and fetal liver (FL) in the midgestation mouse.16-18 Conditional knockout of Gata2 or Gata2 regulatory elements in vascular endothelial cells indicates that Gata2 is essential for hematopoietic cluster formation and HSC generation.7,19,20 Gata2 plays a role in the emergence of cKit-expressing hematopoietic cells from the endothelium.7 Later, as shown in VavCre conditional knockout mice, Gata2 is essential for HSC maintenance,7 thus demonstrating a role for Gata2 as previously recognized in bone marrow LSK HSCs.21
To date, the correlation between Gata2 and hematopoietic cell generation in the embryo has been made in the absence of prospective isolation of viable Gata2-expressing cells.16 Although some hematopoietic cells remain in the embryo in the absence of Gata2,5-8 the identity of these cells is unknown. In this study, to further understand the requirement for Gata2 in normal hematopoietic development, we create and use a mouse model in which a fluorescent reporter for Gata2 (IRES-Venus knock-in gene) does not affect the normal level or function of Gata2. We demonstrate that all long-term repopulating HSCs and a large percentage of HPCs in the midgestation mouse embryo are Venus+. We isolate and characterize a Venus− HPC population that corresponds to the HPCs found in Gata2-null embryos. Gata2-independent hematopoietic progenitors are functionally less complex and do not follow the same genetic program as Gata2-dependent HPCs.
Materials and methods
Gata2Venus ES cells and mice
Generation of the Gata2-Venus mouse model is described in the supplemental Methods, available on the Blood Web site. In short, an IRES-Venus fragment and a loxP-PGK-Puro-loxP fragment were inserted in the Gata2 3′ untranslated region (UTR). IB10 ES cells were transfected and puromycin selected, and 360 clones were polymerase chain reaction (PCR) screened for Gata2Venus (right arm junction, 2292 bp). Correct integration was verified by Southern blot (left arm) for 2 clones with normal karyotype. Founders were identified by Venus PCR. First-generation G2V offspring were crossed with CAG-Cre mice22 and backcrossed (>10 generations) with C57BL/6.
Mice and embryo production
Gata2+/− mice,5 Ly5.1 (6-8 weeks) and C57BL/6 mice were obtained/maintained (Harlan or locally) and genotyped by PCR (supplemental Methods). Day of plug discovery is E0. Embryos were staged by somite pair (sp): E9.5 = 16 to 28 sp, E10 = 28 to 40 sp, early E10 = 28 to 34 sp, E10.5 = 35 to 40 sp, and E11 = 40 to 50 sp. All mouse experimentation was performed under the UK Animals Scientific Procedures Act 1986 Project License 70/8076 and NL Ethics Committee approval and performed in compliance with Standards for Care and Use of Laboratory Animals.
Immunostaining
Whole-mount conceptuses were stained and imaged23 ; cryosections and cells for flow cytometry were stained8 using anti-CD34-biotin (1:50; BD), anti-Gata3 (1:10 KT122, 111207H09; Absea), anti-Gata4 (1:50 H-112, sc-9053; SantaCruz), and anti-green fluorescent protein antibodies. For flow cytometry,8 cells were stained with anti-CD31 (390; BD), anti-CD34 (RAM34; BD), anti-cKit (2B8; BD), anti-CD41 (MWReg3; SantaCruz), anti-Sca1 (D7; Ebiosciences), and anti-CD16/32 (2.4G2; BD) antibodies and Hoechst 33258 (BD) and analyzed (FACSARIAIII; SORP).
Hematopoietic assays
Venus-sorted E9 to E11 AGM, vitelline+umbilical arteries (VA+UA), PL, and YS and E10 FL or Gata2-deficient E9 to E10 AGM, VA+UA, and PL were seeded in 3.6 mL methylcellulose (1 mL per dish; M3434; Stem Cell Tech) for 10 to 12 days.24 Colonies were counted, isolated, and washed, and Venus expression was examined (FACSARIAIII; Fortessa). Sorted G2V E11 AGM (Ly5.1/Ly5.2) cells were transplanted24 into 9.5-Gy irradiated (C57BL/6X129) F1 recipients (Ly5.2/Ly5.2) together with 2 × 105 spleen cells from the recipient strain. Peripheral blood (PB) donor chimerism was determined by Venus PCR and/or fluorescence-activated cell sorter (FACS) at 4 months after transplantation and scored positive if PB donor chimerism was >10%.
RNA analyses
Detailed RNA procedures are provided in supplemental Methods. The Gene Expression Omnibus data accession number is GSE76254. Briefly, RNA was isolated from E10.5 AGM CD31 and cKit sorted cells with the mirVana miRNA Kit (Ambion), and quality/quantity was assessed using the 2100 Bioanalyzer (Agilent; RNA Nano/Pico chip). For RNA sequencing, samples were prepared with the SMARTER protocol and analyzed on an Illumina HiSeq200 system.25-27 For quantitative reverse transcription-PCR (qRT-PCR), SuperScriptII/III ReverseTranscriptase (LifeTechnologies) was used for first-strand cDNA synthesis. Primers are specified in supplemental Table 1.
Results
Generation and validation of a novel Gata2 reporter mouse model
Previously, analysis of Gata2-expressing cells has been limited to a reporter mouse model that results in Gata2 haploinsufficiency.16 Our approach allows for the expression of the reporter within the Gata2 genomic locus without affecting the levels of Gata2 expression or protein function. This is particularly important because Gata2 haploinsufficiency greatly reduces the number of HS/PCs generated during development.7,8,19,20,28 Briefly, an internal ribosome entry site sequence (IRES) followed by the Venus fluorochrome gene was recombined into the Gata2 3′UTR (Figure 1A) in ES cells. The resulting Gata2Venus (G2V) mice bred normally and showed no overt growth or hematopoietic defects.
To determine whether Venus reporter expression parallels that of Gata2, G2V bone marrow (BM) cells were sorted into Venus-expressing (Venus+) and nonexpressing (Venus−) fractions (Figure 1B). qRT-PCR for Gata2 and Venus transcripts (Gata2V/+ BM) demonstrated that only Venus+ cells express Venus and Gata2 mRNA (Figure 1C). Western blot analysis revealed that equivalent amounts of Gata2 protein were present in Gata2+/+ and Gata2V/V adult BM cells (data not shown). FACS analysis showed that Gata2v/v BM LSK frequency (388/105 cells) are comparable to wild-type (WT) BM LSK frequency (378/105 cells). Importantly, the results of competitive limiting dilution transplantation analyses of Gata2V/V and WT BM cells demonstrate that HSCs are qualitatively and quantitatively normal in this mouse model (Figure 1D). Thus, Venus expression correctly reports Gata2 expression without interfering with its normal expression levels or function.
Gata2 is expressed in emerging aortic hematopoietic cluster cells and other embryonic hematopoietic tissues
Venus expression was examined in midgestation G2V hematopoietic tissues. Flow cytometry revealed that E9 to E11 AGM, YS, PL, and FL contained Venus+ cells (Figure 1E). At E9, 6.28 ± 0.47% of viable YS cells and 1.82 ± 0.31% of AGM cells are Venus+. At E10.5 (when the first HSCs are generated), 3.27 ± 0.52% of AGM cells are Venus+, and this increases to 7.86 ± 1.1% at E11. Table 1 shows the frequencies of Gata2-expressing cells.29,30
Tissue . | Stage . | Number of experiments, embryos analyzed . | % Venus+ cells/tissue . |
---|---|---|---|
AGM | E9, 16-25 sp | n = 4, 22 | 1.82 ± 0.31 |
E10, 28-36 sp | n = 4, 29 | 3.27 ± 0.52 | |
E11, 43-49 sp | n = 4, 25 | 7.86 ± 1.1 | |
FL | E9 | nd | nd |
E10, 28-36 sp | n = 4, 10 | 13.89 ± 0.7 | |
E11, 43-49 sp | n = 4, 19 | 19.27 ± 2.14 | |
YS | E9, 16-25 sp | n = 5, 8 | 6.28 ± 0.47 |
E10, 28-36 sp | n = 6, 9 | 6.17 ± 0.86 | |
E11, 42-46 sp | n = 1, 3 | 5.25 ± 0.59 | |
PL | E9, 17-23 sp | n = 3, 4 | 10.91 ± 0.49 |
E10, 28-35 sp | n = 3, 6 | 10.8 ± 1.92 | |
E11, 42-46 sp | n = 1, 3 | 10.01 ± 4.64 |
Tissue . | Stage . | Number of experiments, embryos analyzed . | % Venus+ cells/tissue . |
---|---|---|---|
AGM | E9, 16-25 sp | n = 4, 22 | 1.82 ± 0.31 |
E10, 28-36 sp | n = 4, 29 | 3.27 ± 0.52 | |
E11, 43-49 sp | n = 4, 25 | 7.86 ± 1.1 | |
FL | E9 | nd | nd |
E10, 28-36 sp | n = 4, 10 | 13.89 ± 0.7 | |
E11, 43-49 sp | n = 4, 19 | 19.27 ± 2.14 | |
YS | E9, 16-25 sp | n = 5, 8 | 6.28 ± 0.47 |
E10, 28-36 sp | n = 6, 9 | 6.17 ± 0.86 | |
E11, 42-46 sp | n = 1, 3 | 5.25 ± 0.59 | |
PL | E9, 17-23 sp | n = 3, 4 | 10.91 ± 0.49 |
E10, 28-35 sp | n = 3, 6 | 10.8 ± 1.92 | |
E11, 42-46 sp | n = 1, 3 | 10.01 ± 4.64 |
The frequency of Venus+ cells within the viable cell fraction of embryonic tissues per embryo is presented. FACS analysis of single cell suspensions of dissected embryonic tissues was performed to define the percentage of cells expressing Venus. AGM contains part of the vitelline and umbilical arteries. Our data for the total number of cells in each tissue (data not shown) correlated with published data for YS and E9.5 AGM29 and for FL.30 The data represent mean ± SEM of 3 to 6 independent experiments, with the exception of E11 YS and PL data (data represent the mean ± standard deviation of 1 experiment. E, embryonic day; n, number of independent experiments, number of individual embryos analyzed; nd, not done; sp, somite pairs.
Whole-mount images of E10 and E11 G2V embryos immunostained with anti-CD31 antibody (marks all endothelial cells and hematopoietic cluster cells) shows Venus+ cells along the aorta (DA). Venus+ cells are also observed in cells of the neural tube (NT), olfactory bulb (OB), and FL (Figure 2A-C). In the E10.5 AGM region (4,6 diamidino-2-phenylindole [DAPI] and CD31 stained [blue and red, respectively]), Venus expression is found in endothelial and hematopoietic cluster cells mainly on the ventral side of the DA and in the urogenital (UG) region (Figure 2D-E). In the FL, Venus-expressing cells are found in a punctate distribution pattern (Figure 2D,F). At E9, Venus is expressed in some of the CD34+ (red) endothelial cells of the paired aorta (Figure 2G, arrowhead) and also in some of the endothelial and hematopoietic cluster cells of the UA (Figure 2H). Venus continues to be expressed at E11 in some aortic endothelial cells and emerging/other hematopoietic cluster cells (Figure 2I, arrowhead). The E10.5 YS shows Venus expression in some of the CD31+ (red) endothelial cells. Overall, Venus expression is similar to what has been previously documented for Gata2 in situ hybridization analysis.31,32 Thus, our model allows for the prospective identification, isolation, and characterization of Gata2-expressing cells during normal development.
All HSCs, but not all HPCs, express Gata2
To test for HSC activity, E11 AGM Venus+ and Venus− cells were transplanted into irradiated adult recipients. All long-term repopulating HSCs were found in the Gata2Venus-expressing fraction (Figure 3A). Nine of 19 recipients receiving Venus+ cells were engrafted (15-71%), whereas none of the 14 Venus− recipients showed donor-derived hematopoietic cells. These HSCs were multilineage repopulating (supplemental Figure 1) and self-renewing (8 repopulated of 12 transplanted with 3 × 106 BM cells from primary repopulated mice; n = 4).
The relationship between Gata2 expression and HPC function was also examined. E9 and E10 AGM Venus+ and Venus− cells were plated in the colony forming unit–culture (CFU-C) assay. High enrichment of HPCs was found in the Venus+ fractions (Table 2). Surprisingly, HPCs were found also in the Venus- fraction, although there were very few. At E9, the Venus+ and Venus− fractions, respectively, yielded 8.0 ± 2.1 and 0.4 ± 0.4 CFU-C per AGM. CFU-C numbers increased in the Venus+ cell fraction at E10 (69.0 ± 7.1) and the Venus− fraction increased to 8.0 ± 2.3 CFU per AGM. However, bi- and multipotent progenitors were found only in the Venus+ fraction (Figure 3B).
Tissue . | Stage . | Number of experiments, embryos analyzed . | CFU-C/tissue/sorted cell fraction . | |
---|---|---|---|---|
Venus− . | Venus+ . | |||
AGM | E9, 20-23sp | n = 2, 5 | 0.4 ± 0.4 | 8.0 ± 2.1 |
E10, 32-35sp | n = 2, 6 | 8.0 ± 2.3 | 69.0 ± 7.1 | |
VA+UA | E9, 20-23sp | n = 2, 5 | 0.8 ± 0.5 | 2.4 ± 1.0 |
E10, 32-35sp | n = 2, 6 | 7.0 ± 1.5 | 197.7 ± 36.9 | |
YS | E9, 20-23sp | n = 2, 5 | 270.0 ± 69.8 | 1020 ± 137.3 |
E10, 32-35sp | n = 2, 5 | 130.0 ± 25.5 | 1252.0 ± 156.0 | |
PL | E9, 20-23sp | n = 2, 5 | 4.0 ± 1.3 | 22.0 ± 2.5 |
E10, 32-35sp | n = 2, 6 | 14.5 ± 5.3 | 251.7 ± 32.4 | |
FL | E10, 32-35sp | n = 2, 6 | 20.5 ± 7.6 | 248.7 ± 75.2 |
Tissue . | Stage . | Number of experiments, embryos analyzed . | CFU-C/tissue/sorted cell fraction . | |
---|---|---|---|---|
Venus− . | Venus+ . | |||
AGM | E9, 20-23sp | n = 2, 5 | 0.4 ± 0.4 | 8.0 ± 2.1 |
E10, 32-35sp | n = 2, 6 | 8.0 ± 2.3 | 69.0 ± 7.1 | |
VA+UA | E9, 20-23sp | n = 2, 5 | 0.8 ± 0.5 | 2.4 ± 1.0 |
E10, 32-35sp | n = 2, 6 | 7.0 ± 1.5 | 197.7 ± 36.9 | |
YS | E9, 20-23sp | n = 2, 5 | 270.0 ± 69.8 | 1020 ± 137.3 |
E10, 32-35sp | n = 2, 5 | 130.0 ± 25.5 | 1252.0 ± 156.0 | |
PL | E9, 20-23sp | n = 2, 5 | 4.0 ± 1.3 | 22.0 ± 2.5 |
E10, 32-35sp | n = 2, 6 | 14.5 ± 5.3 | 251.7 ± 32.4 | |
FL | E10, 32-35sp | n = 2, 6 | 20.5 ± 7.6 | 248.7 ± 75.2 |
Number of total CFU-C (mean ± SEM) per tissue per G2V embryo for the sorted Venus− and Venus+ cell fractions at E9 and E10.
HPC activity was also examined in the Venus+ and Venus− fractions of other hematopoietic tissues (Table 2). E9 and E10 VA+UA (Figure 3C), YS (Figure 3D), and PL (Figure 3E), and E10 FL (Figure 3F) contained progenitors in both fractions. Most HPCs were Venus+. The greatest number of CFU-C arising from Venus− cells was found in the E9 YS (270.0 ± 69.8 CFU-C/YS). These data indicate that some HPCs are not expressing Gata2. BFU-E, CFU-G, and CFU-M were the predominant colony types in both fractions, and in contrast to the Venus+ fractions, the Venus− fractions of VA+UA, YS, FL, and PL yielded few or no CFU-GEMM. Thus, all AGM HSCs express Gata2, Gata2 expression is associated with immature HPCs, but not all HPCs are Gata2 expressing.
Some HPCs and vascular cluster cells are formed in the absence of Gata2
Because the Venus− fractions of midgestation G2V hematopoietic tissues contain CFU-C, we tested whether such hematopoietic progenitors are present in Gata2-deleted embryos. CFU-Cs were detected in the Gata2−/− E9 AGM, E10 AGM (Figure 4A), and E10 VA+UA (Figure 4B), although significantly fewer compared with WT (Table 3). Gata2+/− tissues also contained fewer CFU-Cs compared with WT. The E9 Gata2−/− YS contained the most CFU-C (64.4 ± 12.2; Figure 4C). In VEC-Cre:Gataf/f embryos, E10 PL showed significantly decreased CFU-C numbers (Figure 4D), as did E10 AGM and YS7 compared with WT. The CFU-C remaining in Gata2−/− embryos are predominantly CFU-G and CFU-M. Very few Gata2−/− CFU-GM and no CFU-GEMM were observed. These data support and validate our findings in G2V embryos that not all HPCs are Gata2-expressing, Gata2-independent progenitors exist in each of the early hematopoietic tissues, and the Gata2-expressing cell fraction is more enriched in multipotent progenitors.
Tissue . | Stage . | Genotype . | |||||
---|---|---|---|---|---|---|---|
WT . | Gata2+/− . | Gata2−/− . | |||||
AGM | E9, 20-23sp | 8.5 ± 2.3 | n = 1, 4 | 2.25 ± 0.9* | n = 1, 4 | 3.0 ± 1.0 | n = 1, 2 |
E10, 28-34sp | 62.4 ± 20.8 | n = 3, 7 | 24.0 ± 7.4* | n = 3, 13 | 16.7 ± 6.8* | n = 3, 9 | |
VA+UA | E9, 20-23sp | 7.0 ± 0.6 | n = 1, 4 | 2.25 ± 0.6* | n = 1, 4 | 5.0 ± 2.0 | n = 1, 2 |
E10, 28-34sp | 241.1 ± 67.8 | n = 3, 7 | 78.7 ± 12.3** | n = 3, 13 | 12.0 ± 2.5** | n = 3, 9 | |
YS | E9, 20-23sp | 772.5 ± 85.3 | n = 1, 4 | 500.0 ± 71.5* | n = 1, 4 | 64.4 ± 12.2** | n = 1, 2 |
E10, 28-34sp | 918.6 ± 147.9 | n = 3, 7 | 584.6 ± 89.0** | n = 3, 9 | 25.7 ± 8.0** | n = 2, 3 | |
Gata2f/+or Gata2f/f | VEC-Cre:Gata2f/+ | VEC-Cre:Gata2−/− | |||||
FL | E10, 30-34sp | 115.8 ± 35.0 | n = 2, 5 | 83.0 ± 20.8 | n = 2, 6 | 22.2 ± 6.7* | n = 2, 6 |
PL | E10, 30-34sp | 406.0 ± 134.0 | n = 2, 5 | 81.0 ± 25.0* | n = 2, 6 | 12.0 ± 4.0* | n = 2, 6 |
Tissue . | Stage . | Genotype . | |||||
---|---|---|---|---|---|---|---|
WT . | Gata2+/− . | Gata2−/− . | |||||
AGM | E9, 20-23sp | 8.5 ± 2.3 | n = 1, 4 | 2.25 ± 0.9* | n = 1, 4 | 3.0 ± 1.0 | n = 1, 2 |
E10, 28-34sp | 62.4 ± 20.8 | n = 3, 7 | 24.0 ± 7.4* | n = 3, 13 | 16.7 ± 6.8* | n = 3, 9 | |
VA+UA | E9, 20-23sp | 7.0 ± 0.6 | n = 1, 4 | 2.25 ± 0.6* | n = 1, 4 | 5.0 ± 2.0 | n = 1, 2 |
E10, 28-34sp | 241.1 ± 67.8 | n = 3, 7 | 78.7 ± 12.3** | n = 3, 13 | 12.0 ± 2.5** | n = 3, 9 | |
YS | E9, 20-23sp | 772.5 ± 85.3 | n = 1, 4 | 500.0 ± 71.5* | n = 1, 4 | 64.4 ± 12.2** | n = 1, 2 |
E10, 28-34sp | 918.6 ± 147.9 | n = 3, 7 | 584.6 ± 89.0** | n = 3, 9 | 25.7 ± 8.0** | n = 2, 3 | |
Gata2f/+or Gata2f/f | VEC-Cre:Gata2f/+ | VEC-Cre:Gata2−/− | |||||
FL | E10, 30-34sp | 115.8 ± 35.0 | n = 2, 5 | 83.0 ± 20.8 | n = 2, 6 | 22.2 ± 6.7* | n = 2, 6 |
PL | E10, 30-34sp | 406.0 ± 134.0 | n = 2, 5 | 81.0 ± 25.0* | n = 2, 6 | 12.0 ± 4.0* | n = 2, 6 |
Number of total CFU-C (mean ± SEM) per tissue shown for WT, Gata2 germline, and conditional knockout embryos at E9 and E10.
P < .05.
P < .01.
Because hematopoietic clusters appear in the VA and UA prior to appearance in the aorta, and are larger than in the AGM,33 we further examined these vessels. Whole-mount microscopic analysis demonstrates that clusters form in the absence of Gata2. The number and size of cKit+ hematopoietic clusters in early E10 Gata2+/− and Gata2−/− VA+UA are decreased compared with WT (Figure 4E). The number of cKit+ cells decreases 20-fold in the E10 Gata2−/− VA+UA (Figure 4F) in correspondence to the decrease in VA+UA CFU-C (Figure 4B), suggesting that these emerging cKit+ hematopoietic cluster cells are part of the cohort of Gata2-independent HPCs.
Alternative genetic program is expressed in Venus− hematopoietic cells
The molecular basis for the functional differences observed in Gata2-dependent and -independent HPCs was examined by RNA sequencing. As most CD31+Venus− HPCs showed cKit intermediate expression, we compared this population to CD31+Venus+cKitint HPCs (Figure 5A). Gene set enrichment analysis on genes sorted by log ratio of Venus+ vs Venus− FPKMs revealed that genes in the Ras signaling pathway were significantly enriched in the Venus+ compared with the Venus− fraction (Figure 5B). Genes upregulated by Ras were enriched in the Venus+ fraction, and highly upregulated genes included Kras, Grb2 (Ras adaptor), and Sos1 and Sos2 (RasGEF activators) (Figure 5C). Genes downregulated by Ras were enriched in the Venus− fraction. RasGAP gene (renders Ras inactive) Rasa2 was highly upregulated in the Venus− fraction, whereas Rasa1 and Rasa3 were highly upregulated in the Venus+ fraction. Rasa4 and NF1 were expressed to similar levels. Also, Venus+ HPCs showed increased levels of CREB and CBP expression compared with Venus− HPCs and express protein kinase A catalytic subunit genes, suggesting that Venus+ HPCs have the potential to activate CREB target genes. Gata2 has CREB response element consensus sites (−3 kb, −300 bp upstream transcription start site), suggesting that it is a downstream target.34,35 As Gata2 is a Notch target,18 a two- to fourfold higher expression of Notch1 and Notch4 was found in the Venus+ fraction (Figure 5D). Moreover, Snw1 and Maml1 (transcriptional coactivators in the Notch pathway that interact with Notch) were upregulated (2- and 30-fold, respectively) in Gata2-expressing HPCs.
Because Venus− HPCs are mainly restricted in their differentiation potential to the macrophage and granulocytic lineages, we evaluated their similarity to YS-derived erythromyeloid progenitors (EMPs) that give rise to tissue-resident macrophages. Flow cytometric analysis for EMP markers36 showed that 3.89% of E10 YS and 0.72% of E10 AGM cells were EMPs (Sca1−cKit+CD41+CD16/32+). The majority of EMPs were Venus− (74% in YS, 85% in AGM; Figure 5F). At E11, the frequency of EMPs in the E11 YS and AGM decreased to 2.11% and 0.11%, respectively (Figure 5F), with 60% of YS and 77% of AGM EMPs now being Venus+. Published transcriptome data on mouse YS EMPs show the low expression of several chemokine receptors/ligands (Cx3cr1, Cx3cl1, Ccl2, Ccr1, Ccl9, and Ccr7).37 The expression of these genes was low or absent in Venus− AGM cells compared with Venus+ cells (Figure 5F). Also, Cxcr4 (highly expressed in EMPs) was highly expressed in Venus− AGM cells compared with Venus+ cells. These results suggest that the Venus− population shares similarities to EMPs at the transcription level.
Analysis of FPKMs for heptad transcription factors previously described as expressed in AGM HSCs and HPCs1-3 showed expression in both the Venus+ and Venus− AGM fractions (data not shown). Also, other Gata factors were expressed in both fractions. In the mouse, Gata1, 2, and 3 are hematopoietic transcription factors, whereas the Gata4, 5, and 6 factors are not directly related to hematopoiesis. qRT-PCR performed on E10.5 AGM CD31+cKit+ cells (Figure 6A) confirmed Gata3 expression by both the Venus+ and Venus− fractions, and Gata4 was significantly higher in the Venus− fraction. Gata1, 5, and 6 transcripts were low/not detected. Immunostaining of E10.5 AGM (WT) showed Gata3-expressing cells in the mesenchyme underlying the ventral aspect of the aorta, aortic endothelial cells, and some cells emerging from the aortic wall (Figure 6B) in agreement with Fitch.38 G2V E10.5 AGM confirmed that some aortic endothelial cells coexpress Gata3 and Gata2 (Figure 6C). Gata4 expression was also found in aortic endothelial (CD34+) cells, but it did not overlap with Gata2 expression (Figure 6D). Together, these results suggest that Gata3 and/or Gata4 may provide some function in Gata2-independent hematopoietic cells.
Discussion
In this study, we prospectively enriched and characterized Gata2-dependent and -independent HPC subsets from our novel Gata2Venus reporter mouse. Molecular analyses, together with the fact that some vascular hematopoietic cluster cells and HPCs persist in the absence of Gata2 expression, suggest that an alternative genetic program exists for the production of HPCs. The transcriptome differences observed between Venus+ and Venus− HPCs may offer possibilities for pathway modifications to achieve the programming complexities necessary for the generation/function of normal definitive HPCs and provide insights into the factors involved in myeloid leukemogenesis.
Gata2 expression in the developing hematopoietic system
We showed the temporal and quantitatively coordinate transcription of Venus and Gata2 in our G2V mouse model. The strategy used39 eliminates expression level and protein alterations that affect HP/SC development. In G2V embryos, we showed that the cells with the most robust and complex hematopoietic potential (all HSCs and most HPCs) are Gata2 expressing. Imaging and FACS analyses of G2V embryos confirm that Gata2 is expressed in all hematopoietic sites during midgestation and that the numbers of Gata2-expressing cells reflect the developmental and temporal hematopoietic changes occurring in each site. At E9, Gata2-expressing cells are found predominantly in the YS, which at this time produces the highest numbers of the hematopoietic progenitors (EMP) in the conceptus. Slightly later as hematopoiesis begins in the AGM and FL, the numbers of Gata2-expressing cells also increase. The highest numbers of CD31+cKit+ cluster cells are found in the aorta, VA, and UA at E10.5, as quantitated by whole-mount embryo imaging.33 Most, but not all, hematopoietic cluster cells express Gata2, and Gata2 expression may be downregulated as HPCs differentiate. However, we found some hematopoietic cluster cells and HPCs in the E10 Gata2−/− vasculature, confirming the existence of Gata2-independent HPCs.5
Importantly, Gata2 is expressed in the endothelial cells of the DA. Already at E8.5, endothelial cells lining the paired dorsal aortae express Gata2, and it continues to be expressed in the E10.5 aorta when HSCs are generated, thus highlighting an involvement of Gata2 in the hemogenic program of endothelial cells. Data in VE-cadherin conditional Gata2-deficient mice and other models7,19,28,40 strongly support the notion that Gata2 is required in hemogenic endothelium for the emergence of HSCs, as does the morpholino knockdown of Gata2b in zebrafish.41
Gata2 and the relationship with hematopoietic function
Prospective isolation and in vivo transplantation showed that all HSCs are Gata2 expressing. In contrast, some HPCs are present in the Venus− cell fractions of G2V hematopoietic tissues and Gata2−/− hematopoietic tissues. In both cases, the HPCs are restricted in their differentiation potential to predominantly the macrophage and granulocytic lineages. Currently, the EMP population is of high interest as a novel hematopoietic cell subset providing tissue resident macrophages.37,42-44 Our FACS data revealed that EMPs are mainly in the Venus− cell population of E10 YS and AGM and increased in the Venus+ population at E11. Chemokine receptor/ligand gene sets obtained from a study on EMP/microglia transcriptome comparisons allowed us to find similarities in chemokine receptor/ligand expression between EMPs and the Venus− HPC fraction.
Despite prevalence of EMPs in the Venus− cell population in E10 YS and AGM, definitive progenitors are largely Venus+. The coexistence of these HPC subsets highlights the fact that there is more diversity in the types of progenitors generated in the embryo than was previously appreciated. In support of this are recent data from ES cell hematopoietic differentiation cultures suggesting that there are 2 different hemogenic endothelial cell subsets45 and the fact that, in vivo, the AGM, VA/UA, YS, PL, and head are all hemogenic tissues.36,46-49
The highest number of Venus− HPCs was found in E9 and E10 YS (270.0 ± 69.8 and 130.0 ± 25.5 CFU-C, respectively) compared with other tissues (PL, AGM, VA+UA). It is clear that Gata2 has an important role in EHT in the hemogenic endothelial cell compartment before or during the generation/emergence of hematopoietic cells, as evidenced by the decrease (but not absence) in the hematopoietic cluster cells in Gata2−/− aorta, VA, and UA. However, it is as yet unclear at what frequency EHT occurs in the YS, thus raising the possibility that Gata2-independent HPCs arise differently than Gata2-dependent HPCs (perhaps directly from hemangioblasts50 ).
We found differences in the number of CFUs from E9 YS Venus- cells (270.0 ± 69) and Gata2−/− cells (64.4 ± 12.2; Figures 3D and 4C; Tables 2 and 3). The fourfold lower CFU number is likely related to observations (ours and others) that colonies from Gata2−/− embryos, YS explants, and ES cell differentiations were smaller/less proliferative than WT colonies, due to the complete absence of Gata2.6,7 Venus− cells are not defective for Gata2, and the resulting colonies are normal in size. Whereas at the time of sorting they did not express Gata2, Gata2 expression could initiate after seeding Venus− HPCs in methylcellulose, and cells thus undergo normal proliferation/differentiation. To test whether Venus− HPCs can convert to Venus+ cells, we analyzed Venus expression in colonies derived from sorted YS fractions after 10 days of differentiation (Figure 7). Venus expression was found in colonies derived from both fractions, indicating that a portion of Venus− cells start to express Gata2 during formation of a hematopoietic colony. Interestingly, colonies derived from Venus+ cells showed a Gr1+ and Mac1+ phenotype, whereas Venus−-derived colony cells were only Mac1+. This demonstrates that Gata2 is not necessary for a subset of HPCs and that Gata2 promotes more complex hematopoietic function in other progenitor subsets.
Gata3/Gata4 redundancy in Gata2-independent progenitors
The expression of Gata3 and Gata4 in Venus− AGM HPCs and aortic endothelial cells is intriguing and highlights the potential redundancy of Gata transcription factors in hematopoietic cell generation. Gata2 and Gata3 can partially rescue the erythroid phenotype in Gata1-deficient mice,51,52 and although recently it was suggested that Gata3 is redundant in HSCs,53 others clearly show that it regulates HSC cell cycle entry54 and self-renewal.55 The fact that Gata3 and Gata4 are expressed in Gata2-nonexpressing enriched HPCs suggests that they may function in this early progenitor subset. Gata3-deficient embryos show decreased numbers of FL HP/SCs.56 Gata3 affects HSC development non-cell autonomously by activating the expression of Th (tyrosine hydroxylase) and hence, catecholamine production in the ventro-lateral cells of the sympathetic nervous system underlying the embryonic aorta.38 HSC production was rescued when catecholamines were administered to the pregnant dams. These investigators also found that some aortic endothelial cells express a Gata3LacZ reporter, leaving open the possibility of a direct and overlapping role for Gata3 in some HPCs.
Much less is known concerning Gata4 in hematopoietic development. In zebrafish, there is a close relationship between anterior hemangioblasts and cardiac precursors.57 Together with Gata5 and Gata6, Gata4 specifies these 2 anterior mesoderm derivatives. In mice, Gata4 is a key component of the cardiac developmental program, with close associations between cardiac, vascular, and hematopoietic lineages.58-60 Moreover, a subset of mouse endocardial and YS endothelial cells express cardiac markers, possess hemogenic potential, and give rise to transient definitive erythroid/myeloid progenitors.61 Our results suggest the Gata4 aortic endothelial cells and Venus− Gata4-expressing HPCs may be derivatives of mesodermal cells with a genetic program that retains cardiac-vascular-hematopoietic potential and can produce HPCs. Further examination of double reporter and deficient mice should reveal the overlapping and/or redundant roles of these Gata factors.
Gata2 as a pivotal regulator of complex hematopoietic function
RNA sequence comparisons of the functionally distinct Venus+ and Venus− HPC subsets revealed a strong upregulation of Kras and Ras pathway genes in Venus+ HPCs. This pathway is particularly important in cell differentiation, acting as a molecular switch to relay extracellular growth signals.62 Kras mutations confer a competitive-repopulating advantage to BM HSCs in transplantations and initiate leukemia in mice.63 In humans, Kras mutations (together with other cooperating gene mutations) are prevalent in patients with various forms of myelomonocytic and myeloid leukemia.64 Interactions between oncogenic Ras and Gata2 have been proposed.65 The normal function of Kras has not yet been explored fully. However, conditional deletion of Kras by Vav-Cre or Mx1-Cre does not affect HSCs or the adult hematopoietic system.66 However, chimeric mice, produced by Kras−/− ES cell blastocyst injection, show no contribution of Kras−/− cells to the hematopoietic system, suggesting that Kras may be important during the embryonic development of the hematopoietic system but not after its generation.
The low expression of Notch1, Notch4, and coactivators in the Venus− compared with the Venus+ HPC fraction supports the fact that early hematopoietic cells are generated independent of this signaling pathway or implies that these are differentiated cells that have turned off Notch signaling.67,68 Others have shown that Notch1 deletion impairs the development of HSCs and angiogenesis,69 but not YS primitive or definitive hematopoiesis. Moreover, Gata2 expression in the aortic endothelium is lost when Jagged1 (ligand) is deleted.70 Our data demonstrate a direct relationship for Notch and Gata2 expression, strongly supporting a pivotal role for this pathway in the generation of functionally complex hematopoietic cells. In the absence of Notch signaling, less complex HPCs emerge in the AGM or are immigrants from the YS.18,69,70 In addition, our observed upregulated expression of some CREB genes in Venus+ HPCs supports the involvement of these regulators in definitive hematopoietic cell generation.71 Our Gata2Venus model, in combination with recently reported Gata2 distal enhancer-Evi1 mouse model, will allow for a direct examination of the cells relevant to leukemogenesis.12
In conclusion, we enriched, localized, and characterized Gata2-dependent and -independent subsets of hematopoietic progenitors in Gata2Venus embryos. The combination of this reporter with other reporter and knockout models will lead to a better understanding of the role of Gata2 (and other factors) in the development and function of multipotential HP/SCs in health, leukemogenesis, and reprogramming.
The online version of this article contains a data supplement.
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
Acknowledgments
The authors thank Prof Jim Palis for critical comments on this manuscript; Dr Dorota Kurek for providing anti-Gata3 and Gata4 antibodies; Drs Dorota Kurek and Mihaela Crisan for immunostaining support; Dr Siska Driegen for ES cell culture support; and Dr Derk ten Berge for providing Wnt for ES cell cultures.
This work was supported by ZonMw (Dutch Medical Research Council Grant 911-09-036), FES Netherlands Institute for Regenerative Medicine (101675), the National Institutes of Health, National Institute of Diabetes and Digestive and Kidney Diseases (RO37 DK54077), ZonMw TOP (91211068), European Research Council Advanced Grant (341096), and Landsteiner Foundation for Blood Transfusion Research (1344).
Authorship
Contribution: P.K., E.d.P., C.E., M.-L.K., M.J., C.S.V., and T.Y. performed research; P.S.K. analyzed RNAseq data; R.v.d.L. performed/analyzed flow cytometric data; D.M. provided reagents; P.K., E.d.P., C.E., D.M., and E.D. designed experiments and analyzed and interpreted the data; and P.K., E.D., E.d.P., and C.E. wrote the manuscript.
Conflict-of-interest disclosure: The authors declare no competing financial interests.
Correspondence: Elaine Dzierzak, University of Edinburgh, Centre for Inflammation Research, Queens Medical Research Institute, 47 Little France Crescent, Edinburgh EH16 4TJ, United Kingdom; e-mail: e.dzierzak@erasmusmc.nl or elaine.dzierzak@ed.ac.uk.
References
Author notes
P.K., E.d.P., and C.E. contributed equally to this work.