Functional Analysis of Human Hematopoietic Stem Cell Gene Expression Using Zebra

干细胞之家 - 中国干细胞行业门户第一站 › 干细胞之家论坛 › 干细胞文献资源库 › 《科学公立图书馆生物学》 › Functional Analysis of Human Hematopoietic Stem Cell ...



免疫细胞治疗专区	欢迎关注干细胞微信公众号

12 3 4 5 6 7 8 9 10 ... 24 下一页

返回列表

查看: 308077\|回复: 231	go Functional Analysis of Human Hematopoietic Stem Cell Gene Expression Using Zebra

春天的风筝

新手上路

Rank: 1

积分: 0
威望: 0
包包: 0

发表于 2009-4-23 08:50 |显示全部帖子

1 Department of Medicine, Division of Hematology, Oncology, and Transplantation, and Stem Cell Institute, University of Minnesota, Minneapolis, Minnesota, United States of America,2 Genetics, Cell Biology, and Development and Arnold and Mabel Beckman Center for Transposon Research, University of Minnesota, Minneapolis, Minnesota, United States of America,3 Cancer Center Bioinformatics Division, University of Minnesota, Minneapolis, Minnesota, United States of America

Although several reports have characterized the hematopoietic stem cell (HSC) transcriptome, the roles of HSC-specific genes in hematopoiesis remain elusive. To identify candidate regulators of HSC fate decisions, we compared the transcriptome of human umbilical cord blood and bone marrow CD34 CD33CD38Rholoc-kit  cells, enriched for hematopoietic stem/progenitor cells with CD34 CD33CD38Rhohi cells, enriched in committed progenitors. We identified 277 differentially expressed transcripts conserved in these ontogenically distinct cell sources. We next performed a morpholino antisense oligonucleotide (MO)-based functional screen in zebrafish to determine the hematopoietic function of 61 genes that had no previously known function in HSC biology and for which a likely zebrafish ortholog could be identified. MO knock down of 14/61 (23%) of the differentially expressed transcripts resulted in hematopoietic defects in developing zebrafish embryos, as demonstrated by altered levels of circulating blood cells at 30 and 48 h postfertilization and subsequently confirmed by quantitative RT-PCR for erythroid-specific hbae1 and myeloid-specific lcp1 transcripts. Recapitulating the knockdown phenotype using a second MO of independent sequence, absence of the phenotype using a mismatched MO sequence, and rescue of the phenotype by cDNA-based overexpression of the targeted transcript for zebrafish spry4 confirmed the specificity of MO targeting in this system. Further characterization of the spry4-deficient zebrafish embryos demonstrated that hematopoietic defects were not due to more widespread defects in the mesodermal development, and therefore represented primary defects in HSC specification, proliferation, and/or differentiation. Overall, this high-throughput screen for the functional validation of differentially expressed genes using a zebrafish model of hematopoiesis represents a major step toward obtaining meaningful information from global gene profiling of HSCs.! O6 H9 x! b2 u% L+ |9 @
. V5 O' [/ N7 u* u) n
Introduction- y3 x! x0 F7 I# u

Hematopoiesis is the process by which hematopoietic stem cells (HSCs) give rise to all hematopoietic lineages during the lifetime of an individual. To sustain lifelong hematopoiesis, HSCs must self-renew to maintain or expand the HSC pool [1], and they must differentiate to form committed hematopoietic progenitor cells (HPCs) that progressively lose self-renewal potential and become increasingly restricted in their lineage potential. A combination of extrinsic and intrinsic signals are thought to converge to regulate HSC differentiation versus self-renewal decisions, but the molecular mechanisms that regulate these processes are poorly understood [2].
2 I" _4 I' l7 ?& z/ b9 F2 J: c
A multitude of cytokines have been cloned that affect HSCs and HPCs; however, to date none of these, alone or in combination, can induce the symmetrical, self-renewing HSC cell division in vitro that is required for HSC expansion. Recently, several novel regulators of HSC fate decisions have been identified. For instance, overexpression of Hoxb4 results in expansion of murine and human HSCs with an increased competitive repopulation potential [3–5]; novel extrinsic regulators implicated in self-renewal of HSCs include Notch [6], Wnt [7,8], and the morphogens, SHH (sonic hedgehog) [9] and BMP4 (bone morphogenetic protein 4) [9]. While the discovery of these novel regulators provides credence to the hypothesis that extrinsic and intrinsic signals can influence HSC fate, a more global gene and/or protein expression analysis of human HSC should provide additional insight into pathways that support HSC self-renewal.% M# D& {' l- W; I9 e
; i3 |; ^1 k$ p5 W* N
Our current understanding of the expressed gene profile of HSCs comes primarily from murine HSCs that can be purified to near homogeneity [10–14]. The inability to purify human HSCs to similar degrees of homogeneity makes study of the transcriptome of human HSCs more difficult. Human HSCs and HPCs are CD34 positive, while cells that engraft in severe combined immunodeficiency (SCID) mice are enriched in the CD34 Lineage(Lin)CD38 fraction [15]. As fewer than 1/500 CD34 LinCD38 cells can repopulate SCID mice [15], the expressed gene profile of CD34 LinCD38 cells is likely only partially enriched for HSC-specific genes [12,16]. We previously demonstrated that the rhodamine (Rho) 123 and c-kit  subpopulation of CD34 LinCD38 cells (Rholo) cells are highly enriched for primitive HPCs with myeloid–lymphoid initiating cell (ML–IC) capacity relative to CD34 CD38CD33Rhohi (Rhohi) cells [17]. Thus, such selection separates CD34 LinCD38 cells into HSC-enriched and HSC-depleted populations.' @8 S0 }. _, a8 ~

We hypothesized that comparison of the transcriptome of Rholo and Rhohi cells from umbilical cord blood (UCB) and bone marrow (BM) should identify conserved genes and gene pathways that define the human HSC. Because of the inherent limitations of using gene expression data to infer biological gene function, we also assessed the hematopoietic role of these genes in a high-throughput in vivo functional genomics screen in the zebrafish. Using this strategy we have not only identified a series of genes that may represent novel regulators of human HSC fate decisions, but this work also represents, to our knowledge, the first example of a functional genetic screening strategy that is a critical step toward obtaining biologically relevant functional data from global gene-profiling studies.

Results/Discussion7 H! b( y- k# j/ [" d+ R; x

ML–ICs Are Highly Enriched in Rholo Compared to Rhohi Cells6 W" u7 E% x4 t- a" _
# R" z8 _2 V7 J# @' z  p: b
The study of human HSCs has been limited since the CD34 LinCD38 fraction of hematopoietic cells, commonly used as an HSC-enriched population, contains fewer than 0.2% SCID-repopulating cells [15], suggesting considerable heterogeneity. We have shown that ML–IC s, single hematopoietic cells that can generate several daughter cells capable of reinitiating long-term myeloid and long-term lymphoid cultures, are highly enriched by selecting the Rholo fraction of CD34 LinCD38 cells. While the Rholo population still only contains 15%–25% ML–ICs and therefore remains heterogeneous, the enrichment factor is 5- to 10-fold greater than CD34 LinCD38 cells [17]. Similar to our previous studies, the ML–IC frequency was greater than or equal to 10-fold higher in UCB Rholo compared to Rhohi cells (Figure S1).
% k# P4 j( g. }/ S1 p) u
Genes Differentially Expressed between Rholo and Rhohi Cells from Both UCB and BM7 X+ j7 Q6 t% D; `2 ^

We hypothesized that comparing genes differentially expressed between Rholo and Rhohi cells from ontogenically distinct sources would identify conserved genes and gene pathways that govern self-renewal and differentiation of human HSCs. The experimental design used is illustrated in Figure 1. We defined differentially expressed probe sets as those with p 6 U* X. n4 O/ O2 n" ?

UCB or adult BM CD34 CD33CD38Rholoc-kit  (Rholo; stem cell enriched) and CD34 CD33CD38Rhohi (Rhohi; stem cell depleted) cell populations were sorted for subsequent global gene expression profiling. Total RNA was isolated from Rholo and Rhohi cell populations prior to linear amplification and labeling for hybridization to the Affymetrix HG-U133 GeneChip set (approximately 45,000 probe sets) and subsequent data analysis.

We identified 2,707 and 4,667 probe sets differentially expressed between Rholo and Rhohi cells from UCB and BM, respectively (see Tables S1 and S2). The fidelity of our microarray results was confirmed using quantitative RT-PCR (Q-RT-PCR) (see Figure S2). We focused our further analysis on 277 unique transcripts, represented by 304 probe sets that were differentially expressed between Rholo and Rhohi cells from both UCB and BM with a fold change greater than 1.5 in either UCB or BM (Table S3).. R3 z) f# G; Z% ^
8 V1 X0 B+ @9 c, ]
Among the conserved genes enriched in Rholo cells, many have been implicated in early hematopoiesis, including CDKN1A, a cell cycle regulator required for maintenance of murine HSCs [18], and ABCB1, the ABC-transporter family member responsible for the Rholo phenotype [17]. Several transcription factors (TFs) known to play a role in early hematopoiesis or leukemogenesis were also identified, including HLF, involved in leukemogenic chromosomal translocations [19] and EVI1, a TF associated with myeloid leukemias [20]. Other TFs without a known role in hematopoiesis were also more highly expressed in Rholo cells, including HMGA2, a high-mobility group gene, and the zinc finger TFs ZNF165, ZNF331, and KLF5. All Rholo-enriched genes are listed in Table S3. As demonstrated in previous HSC gene-profiling studies [12–14], more than 40% of genes enriched in Rholo cells lack a functional annotation, are hypothetical proteins, or are expressed sequence tags, and thus may represent currently uncharacterized regulators of HSC fate decisions (Figure S3).1 ?5 t! R/ k. Z- J. H6 l! O+ X- `

Some genes with well-established roles in HSC self-renewal and early differentiation are not present in the Rholo enriched gene list. However, most of these were differentially expressed in both datasets, but differences did not reach statistical significance. For instance, LMO2 [21] and GATA2 [22], known to be involved in HSC development and self-renewal, were expressed significantly higher in BM Rholo than in Rhohi cells. Although similar trends were seen in UCB Rholo cells, these differences were not statistically significant (Table S4). Conversely, HOXB4 [4] expression was significantly higher in UCB Rholo than in Rhohi cells, but this difference was not statistically significant in BM. Although our stringent criteria for differential expression likely contribute to the omission of some genes that might be differentially expressed, another explanation might be that expression of these genes is maintained when cells differentiate from a Rholo to a Rhohi stage. The latter is consistent with the fact that most known HSC-associated genes were expressed at much higher levels than the normalized average microarray expression level in Rholo and in Rhohi cells from both UCB and BM (Table S4).9 V5 O! n2 G) g

Conserved genes enriched in the Rhohi cells included LEF1, an effector of Wnt signaling expressed in pre-B and T cells [23], and NOTCH2, involved in hematopoietic differentiation cell-fate decision [24]. Several TFs known to play a role in hematopoietic cell differentiation were more highly expressed in Rhohi compared to Rholo cells, including HELLS and MAFB [25,26]. Additional TFs with no known role in hematopoietic development were also enriched in Rhohi cells, such as the zinc finger homeobox gene, ZFHX1B and the polycomb genes, EZH2 and SUZ12, the latter required for germ cell development [27]. Globin (Hb) gene family members were also more highly expressed in UCB and BM Rhohi than in Rholo cells. Consistent with the ontogenic expression patterns of fetal versus adult Hb genes, Hb污 genes were more highly expressed in perinatal UCB Rhohi cells, while Hb汕 genes were more highly expressed in adult BM Rhohi cells. Additional genes enriched in Rhohi cells are listed in Table S3.

Because of functional redundancy among gene families, we examined the data for common differentially expressed gene family members. The Id family of transcriptional repressors [28] was enriched in the Rholo fraction, but was represented by different family members in UCB (ID4) and BM (ID1, ID2, and ID3). Similarly, various H1 and H2 histone genes were enriched in the Rholo fraction in both datasets, but were represented by distinct family members.
6 I8 E# W, _8 k6 s3 ~
We also evaluated whether common differentially expressed genes were concentrated on specific chromosomes. We found that genes were not only concentrated on certain chromosomes, but at specific g-band addresses. Of the genes enriched in Rholo cells, 9% reside at 6p21, a region involved in recurrent chromosomal translocations in myeloid [29] and lymphoid [30] leukemias, and home to the PIM1 oncogene [31]. Six members of the H2B and one member of the H1 histone family, as well as CDKN1A, more highly expressed in Rholo than Rhohi cells, reside at 6p21. The remaining Rholo-enriched genes at 6p21 consist of six class II major histocompatibility complex (MHC) family members and a putative testis-specific zinc finger TF, ZNF165. H1 and H2 histone gene family members [11,13], class II MHC antigens [12,13], and CDKN1A [13] were also found among the genes identified in studies characterizing the transcriptome of murine HSC. The differential expression of such a large number of genes located at this chromosomal address suggests that, like CDKN1A, other genes located at 6p21 with as yet unknown hematopoietic function may play a role in HSC proliferation or differentiation.
, p  I; ^( F8 {- H( J" I: a) g" b
We also compared the genes expressed more highly in Rholoversus Rhohi cells with published gene expression data. Comparison with the study by Ivanova et al. [12] that compared human CD34 LinCD38 with CD34 LinCD38  cells, yielded only seven genes in common: ARMCX2, CRYGD, HLF, KIAA1102, RBPMS, SLCO3A1, and SSBP2. The lack of overlap may not be that surprising as Rholo and Rhohi cells are subpopulations of the CD34 LinCD38 population used by Ivanova et al. Comparison of genes expressed more highly in Rholo versus Rhohi cells with genes expressed more highly in murine side population/KLS/CD34 compared to total BM cells published by Ramalho-Santos et al. [13] identified 16 likely orthologs and 38 common gene family members (Table S5), suggesting that HSC-specific genes are conserved across species.

In vivo functional genomics screen in zebrafish Because gene profiling per se does not prove functional importance, we developed an in vivo functional genomics screen in zebrafish (Figure 2). The zebrafish, Danio rerio, is an ideal organism for high-throughput genetic screens [32] as organogenesis is highly conserved from zebrafish to man [33]. There is abundant evidence that hematopoiesis in zebrafish occurs via a highly conserved genetic program. As in mammals, hematopoiesis in zebrafish occurs via specification of mesoderm to a hemangioblast stage that subsequently commits to either HSC or angioblasts [34], and genes and signals involved in specification (BMP signaling) and commitment (vascular endothelial growth factor signaling; flk1, also known as kdr; lmo2; scl, also known as tal1; gata2; gata1) are conserved from fish to man [35]. This high degree of homology in the genetic control of zebrafish and human hematopoietic development makes genetic screens in zebrafish a powerful tool to elucidate the role of genes in hematopoiesis. Additionally, rapid reverse genetic screens can be accomplished using morpholino antisense oligonucleotides (MOs) to knock down gene expression in the developing zebrafish embryo [36].2 j6 u' P& i7 f: j+ M- y4 y& H  z

The hematopoietic function of differentially expressed candidate genes was determined by injecting MOs into one- to two-cell embryos from gata1

sRed Tg zebrafish to disrupt gene expression. Injected embryos were scored at 30–48 hpf for the presence of DsRed blood cells by fluorescence microscopy. Subsequently, MO-targeted embryos with gross hematopoietic defects were analyzed for the expression of the early hematopoietic markers gata1 and scl by whole-mount in situ hybridization, and the late hematopoietic markers hbae1 and lcp1 by Q-RT-PCR.( [! J3 ]4 t. \. J
' O: p% l" b% p% A$ Z/ ^
From the 277 unique transcripts that were differentially expressed between Rholo and Rhohi cells of both UCB and BM (Table S3), we eliminated genes with known function in hematopoiesis, MHC genes, histones, and genes that are known to play a role in glucose and protein metabolism and RNA and DNA synthesis, resulting in a final list of 158 genes. Of these, we identified a putative zebrafish ortholog for 86, and designed MOs against 61 (Table S6). The 61 MOs were injected in one- to two-cell zebrafish embryos and assayed for effects on blood development. Initial dosing studies identified 16/61 MOs that reduce blood cell production without confounding toxicities (Table 1). The 16 MOs induced a blood defect in more than 70% of embryos in two or more independent injections. Blood defects identified by gata1

sRed transgenic (Tg) fluorescence microscopy were confirmed by Q-RT-PCR of the erythroid-specific hbae1 and myeloid-specific lcp1 transcripts in MO-targeted embryos compared to uninjected controls in three or more independent experiments of n = 5 embryos per experiment. A greater than 2-fold reduction in both erythroid and myeloid gene expression levels were seen for five of seven MOs analyzed (see Table 1; Figures 3 and 4). Thus, the addition of Q-RT-PCR to the screening process provides an independent confirmation and quantitation of the observed phenotypes, thereby limiting the false-positive rate, while maintaining the high-throughput nature of the screen. Additionally, the observed reduction of both erythroid and myeloid gene expression following knockdown of candidate genes is consistent with their presumed roles in HSC fate decisions prior to specification of the common myeloid progenitor. The validity of the Q-RT-PCR analysis was corroborated by analysis of hbae1 and lcp1 transcript levels in gata1 MO-targeted embryos, in which there was a virtually complete loss of hbae1 expression and an almost 2-fold increase in myeloid-specific lcp1 transcripts (Figure 4) consistent with the published expression patterns for these genes following loss of gata1 expression [37]. Analysis of vascular development by injecting MOs into fli1:enhanced green fluorescent protein (EGFP) Tg embryos, which express EGFP in endothelial cells, demonstrated no major abnormalities in vascular morphogenesis or remodeling that precludes the circulation of the remaining blood cells, indicating that the hematopoietic defect is not secondary to a vascular phenotype (Figure 4; Table S7).

(A) The observed frequency of hematopoietic defects in gata1

sRed Tg zebrafish embryos are indicated for two independent spry4-targeted MOs (MO1 and MO2; black bars), a four-base mismatched control MO (MM MO; no bar indicates a 0% frequency), a low-dose injection of MO1 and MO2 individually and in combination (gray bars), and MO2 coinjected with a human SPRY1 DNA expression vector or a GFP control (white bars) (error bars = standard deviation of the mean, p

(B) Representative pictures of the phenotypes seen with spry4 MO1 at 48 h using fli1:EGFP/gata1

sRed double Tg zebrafish embryos that have EGFP vascular endothelial cells and DsRed erythroid cells. The embryos display a more drastic reduction in DsRed blood cells with some blood pooling and pericardial edema at higher MO doses ( MO) without major vasculature defects. The gata1

sRed and fli1:EGFP Tg images are representative of three experiments of n > 40 embryos each./ T% a6 \9 L: k* \3 f; z- L) o

(C) Q-RT-PCR analysis of hbae1 and lcp1 transcripts in spry4 MO1, control mismatch MO and a gata1 MO-injected zebrafish at 48 hpf (p * x4 o* h( a  u1 v

(D) Injection of 30 pg of human SPRY1 DNA expression vector resulted in an expansion of DsRed  hematopoietic cells in the posterior ICM in greater than 50% of successfully injected embryos at 32 hpf (three experiments of n = 30), and a representative bright-field image is pictured with the fluorescence micrograph overlaid (left). Q-RT-PCR analysis of hbae1 and lcp1 transcripts in SPRY1 overexpressing embryos relative to uninjected clutchmate controls (three experiments of n = 5 embryos) (right).

The greater than 20% frequency of blood defects seen in our screen compares very favorably with the 0.5%–1% frequency of hematopoietic phenotypes seen by ethylnitrosourea mutagenesis screens that mutate genes in a near random fashion [38] and the 4% of hematopoietic phenotypes seen in a morpholino-based functional screen of the zebrafish secretome [39] (S. C. Ekker, unpublished data). The high incidence of blood defects also demonstrates that the candidate genes identified by comparing the transcriptome of Rholo and Rhohi cells represent genes with important roles in HSC biology. We were unable to identify a candidate ortholog in zebrafish for 72/158 of the differentially expressed human genes. This may be because currently only one quarter of the zebrafish genome is high-quality finished sequence. Hence, some genes with important roles in hematopoiesis may have been untested. However, from a sample of ten genes that lacked a zebrafish match, only one has a likely ortholog in the Fugu or Medaka sequencing projects, and thus incomplete genome coverage provides only a partial explanation. Alternative possibilities include a reduced level of primary sequence conservation between functional orthologs that would have been missed using our simple comparative genomics criteria, or that a number of genes are not conserved between fish and man, and thus might be less important for the conserved processes of hematopoietic self-renewal and differentiation. Therefore, the relatively high incidence of blood defects in conserved genes may in part reflect ※evolutionary filtering§ in our screen. Consistent with this hypothesis, all zebrafish genes whose mutation resulted in a visible embryonic phenotype identified using a retroviral insertion strategy have a likely human ortholog [40].

Although the frequency of blood defects is high, the screen is not as well suited for the identification of knockdown phenotypes that result in increased HSC proliferation and/or differentiation. We and others have successfully shown that dramatically increased hematopoietic development can be modeled in zebrafish, as is the case for knockdown of the BMP-antagonist chd (chordin) and the corresponding dino mutant [41,42]; however, more modest increases in hematopoietic cell production or skewing of lineage differentiation may be undetectable. Although these caveats may lead to underestimation of the true frequency of genes with a role in early hematopoietic development and differentiation, the screening procedure used has proven effective for extracting functional information from a global gene expression profiling dataset in a high-throughput manner.2 _) C  v4 A3 `: P5 _6 O  [

Of note, viable knockout mice exist for 4/14 genes identified in our functional screen (Sparc, Irak3, Ccr7, and Prkch) [43–46]. One could argue that if a viable knockout mouse exists, the gene of interest may not be important in hematopoiesis. However, lack of an overt hematopoietic phenotype does not preclude a role of a gene in HSC self-renewal and differentiation, as this may only be detectable under conditions where the hematopoietic system is stressed or in transplantation experiments. For example, Hoxb4/ mice develop normally and present with only subtle differences in spleen and BM cellularity [47]. However, the proliferative response of HSC in vitro and in vivo is decreased, consistent with the observation that overexpression of Hoxb4 supports expansion of competitive repopulating units and SCID-repopulating cells [4,5].+ x7 [( A; i. D" w9 s
) M& X* x: M& B& g3 p4 E0 D
Characterization of Sprouty Family Members in Zebrafish Hematopoiesis

To further verify the hematopoietic role of genes identified by gene array, we performed a more extensive evaluation of the zebrafish targeted with an MO against spry4 (spry4 morphant or spry4MO). Although human SPRY1 was differentially expressed, an MO against zebrafish spry4 was used, as it is expressed in the region of the lateral plate mesoderm, the first site of zebrafish hematopoiesis [48], and it was the full-length zebrafish Sprouty gene with the greatest protein homology to human SPRY1. Recently the partial sequence of a potential zebrafish spry1 ortholog was predicted by Ensembl's gene prediction software based on genomic sequence information. However, the single exon that was predicted does not contain an ATG start codon, or a conserved splice donor or acceptor site, and the putative zebrafish spry1 sequence only partially covers the human SPRY1 gene. Moreover, the genomic location of Ensembl's putative zebrafish spry1 is currently not known, and therefore it is not possible to use syntenic relationships to determine the most likely zebrafish ortholog for human SPRY1. At present, there is not sufficient sequence data available to design gain- or loss-of-function experiments for the putative zebrafish spry1, thus precluding an analysis of hematopoietic function in the zebrafish model. Therefore, spry4 is currently the best full-length, MO-targetable candidate ortholog for human SPRY1, and based on our results, at the very least zebrafish spry4 and SPRY1 share a conserved function in hematopoiesis.9 \" t( a3 O1 k9 s1 G# M

To confirm the specificity of MO targeting in the spry4MO, a second spry4 MO of independent sequence and a four-base mismatched spry4 MO were injected into zebrafish embryos. Injection of the independent spry4 MOs induced a hematopoietic phenotype in more than 65% of injected embryos, while the four-base mismatched MO did not induce any phenotypic changes (Figure 4). Q-RT-PCR for hbae1 and lcp1 transcripts also did not show changes in expression in four-base mismatched MO embryos. Furthermore, the two independent spry4 MOs acted synergistically when coinjected (Figure 4). In addition to the blood phenotype, a slight facial outgrowth was seen at a low frequency. Also, a weak dorsalization phenotype was seen at 2-fold higher MO doses based on a decreased somite size.3 v9 M8 u0 S: A0 o( q- Q) ?
  F# c- b2 d6 w' s1 s9 R
To rule out the possibility that the hematopoietic defect observed in the spry4MO was secondary to a vascular defect, we injected spry4 MO into fli1:EGFP Tg zebrafish. While the resulting embryos exhibited minor defects in cardinal vein remodeling and morphogenesis of intersegmental vessels in the posterior tail, there were no major defects in vascular development (Figure 4). The integrity of the vascular network in spry4MO fish was further demonstrated by the unimpeded circulation of the remaining DsRed  blood cells in gata1

sRed Tg spry4MO (data not shown).
' s4 g' E a/ ]( C
Thisse et al. [35] have shown that overexpression of zebrafish spry4 mRNA leads to an expansion of the posterior intermediate cell mass (ICM) [48]. We overexpressed human SPRY1 in gata1

sRed Tg zebrafish embryos, and observed a similar dose-dependent expansion of DsRed  blood cells in the posterior ICM (Figure 4). Q-RT-PCR analysis of embryos overexpressing SPRY1 revealed a 1.7- and 1.4-fold increase hbae1 and lcp1 mRNA levels, respectively, confirming the observed expansion of blood cells in the ICM. Coinjection of human SPRY1 cDNA with spry4 MOs ameliorated the spry4MO phenotype, yet another confirmation of the specificity of MO targeting (Figure 4). The similar hematopoietic phenotypes observed following the overexpression of either human SPRY1 or zebrafish spry4 indicate that these genes encode proteins with a similar functional potential in blood development. In addition to the blood phenotype, the overexpression often caused a reduction and/or curve in the posterior tail (Figure 4) not seen in the overexpression of zebrafish spry4 [48].& g) q' x4 ~( J2 B5 B

Characterization of hematopoietic gene expression in the spry4MO by whole-mount in situ hybridization revealed a reduction in scl expression at four somites (8/15), and virtually no scl (12/20) or gata1 (18/25) expression at 20 somites (Figure 5), consistent with a defect in mesodermal commitment to HSCs and/or HSC proliferation and differentiation. The few hematopoietic cells that are present in the morphant are hemoglobinized based on o-dianisidine staining (data not shown). To determine whether the hematopoietic defects observed in spry4MO were the result of a defect in mesoderm specification during development, we performed whole-mount in situ hybridization for the vasculature-specific flk1 and muscle-specific myod transcripts. At 10 h postfertilization (hpf) there was a slight defect in myod expression (5/10), while myod expression at 26 hpf was comparable to wild type (21/21) (Figure 5). This suggests that there were no major defects in mesodermal commitment in the spry4 morphants. Considering the nearly absent expression of the early hemangioblast and HSC cell marker scl, these results may suggest that the defect induced in the spry4MO occurs prior to HSC specification. However, the normal expression pattern of other mesodermal genes, such as flk1 and myod at 26 hpf in the spry4MO (Figure 5), indicate that the hematopoietic phenotype is not merely a consequence of defective specification of mesoderm. Finally, the hematopoietic gene cmyb, a presumed marker of definitive HSC in zebrafish [49], was absent at 38 hpf (10/12) (data not shown), suggesting that the spry4MO embryos are devoid of definitive HSC.% B3 I" l5 g" ^% D7 K: Z2 i, s

(A) Expression of the early hematopoietic marker scl was reduced at the four-somite stage (arrow) and virtually absent at the 20-somite stage in spry4MO embryos compared to uninjected controls.( w! j5 ~4 a7 \& M: O( V

(B) The more mature hematopoietic marker gata1 was also reduced in spry4MO compared to uninjected zebrafish embryos at 20 somites.% A9 X& y2 Q: F0 y

(C and D) In contrast to hematopoietic genes, the mesodermal-specific flk1 (vasculature) transcripts were expressed at similar levels to uninjected controls at 26 hpf, and myod (muscle) had a slight defect (arrowhead) at 10 hpf, while expression was similar to controls at 26 hpf. All images are representative of greater than or equal to two experiments of n ≡ 8 embryos.: w" E3 w$ z& e- u* N7 h% N- s  i

In vertebrates, Sprouty family members act as antagonists for fibroblast growth factor (FGF), vascular endothelial growth factor, and epidermal growth factor signaling, and they may be involved in feedback regulation, as Sprouty gene expression is induced by activation of these signaling pathways [50]. Sprouty genes antagonize receptor tyrosine kinase signaling at the level of the Ras/Raf/mitogen-activated protein kinase pathway; however, they also can serve as positive regulators of these pathways in some cell types [50]. Therefore, our current hypothesis is that SPRY1 may affect HSC by modulating FGF-mediated, perhaps in combination with other receptor tyrosine kinase-mediated, signaling. In fact, three of the 14 genes that induce a hematopoietic defect in the zebrafish screen, SPRY1,MAFB and SPARC, are all involved in FGF signaling [50–52], thus suggesting a role for FGF in hematopoiesis. Studies are ongoing to confirm a role of SPRY1 in mammalian hematopoiesis, by testing the effect of overexpression of SPRY1 on the repopulating ability of human and murine HSCs. Similar definitive studies in zebrafish and subsequent confirmation in mammalian HSC models are underway for genes that were functionally validated in our zebrafish screen. We believe that the sequential genetic screen in zebrafish followed by confirmation in mammalian models as described here will establish a hematopoietic function for genes identified by gene array analysis in a high-throughput and efficient manner.
/ [* f' E; x$ T8 D7 `0 h5 R* l# O
Materials and Methods4 k: v! `  N# J
6 P8 f; N, E- `' ^0 A% P
Isolation of Rholo and Rhohi cell populations from UCB and BM/ a' c0 k  i) Q! }/ b% f) |
5 M$ W2 ~8 i  F/ ^+ P8 V$ n
Human UCB from full-term delivered infants and BM from healthy donors were obtained after informed consent in accordance with guidelines approved by the University of Minnesota Committee on the Use of Human Subjects in Research. Each biologically distinct replicate was comprised of one to four donors for UCB and individual donors for BM samples. CD34 CD38CD33Rholoc-kit  and CD34 CD38CD33Rhohi fractions were selected by sequential Ficoll-Hypaque separation, MACS column depletion, and fluorescence-activated cell sorting as previously described [17]. Postsort analysis demonstrated that sorted populations contained fewer than 1%–2% contaminating cells from the opposing population (see Figure 1).4 Q& Q2 `/ U; r# T; s1 `. j/ D7 N
! y. o* S+ G) [8 p* x( f
Determination of ML–IC frequencies

ML–IC frequencies for UCB samples (n = 3) were determined as previously described [17]. An ML–IC was defined as a single cell that gave rise to at least one LTC–IC and one NK–IC. Results are presented as ML–IC frequency ± standard deviation of the mean.  v0 E4 x, q7 O  N1 c
0 J  @0 s. k! V$ ^6 D
Processing of RNA samples and oligonucleotide microarray analysis

Total cellular RNA was isolated from UCB (n = 5) and BM (n = 4) Rholo and Rhohi cells using the PicoPure RNA Isolation Kit (Arcturus, Mountain View, California, United States) per the manufacturer's instructions. Seven to 10,000 Rholo and Rhohi cells were sorted directly into 100 μl extraction buffer (XB) provided with the PicoPure RNA Isolation Kit prior to RNA isolation. Labeled cRNA was generated by one round of IVT-based, linear amplification using the RiboAmp OA RNA Amplification Kit (Arcturus) followed by labeling with the Enzo Bioarray HighYield RNA Transcript Labeling Kit (Enzo Life Sciences, Farmingdale, New York, United States) according to manufacturer's instructions. Samples were hybridized to Affymetrix HG-U133 A and B chips (Affymetrix Inc., Santa Clara, California, United States) were washed, and scanned at the University of Minnesota Affymetrix Microarray Core Facility as described in the Affymetrix GeneChip Expression Analysis Technical Manual.
6 R7 f- \: w& w) {
Oligonucleotide microarray data analysis

Affymetrix HG-U133 GeneChips were processed using GeneData Refiner software (GeneData, San Francisco, California, United States) to assess overall quality. Feature intensities for each chip were condensed into a single intensity value per gene using the Affymetrix Statistical Algorithm (MAS 5.0) with 而 = 0.015, 汐1 = 0.04, 汐2 = 0.06, and a scaling factor of 500. Expression data was analyzed using GeneData's Expressionist and Microsoft Excel (Microsoft, Redmond, Washington, United States). Differential gene expression for comparison of Rholo versus Rhohi cells was defined by using a paired Student's t-test with a threshold of p
& h% h7 f3 j1 U  X8 b. N) C
Microarray Q-RT-PCR confirmation$ D; {( X. D+ l$ E" g4 _/ H
3 m) c( ^7 X1 Q' h% u
Labeled cRNA was reverse transcribed to generate cDNA using SuperScript II Reverse Transcriptase (Invitrogen, Carlsbad, California, United States) according to manufacturer's instructions. Q-RT-PCR was performed using the ABI Prism 7000 Sequence Detection System (Applied Biosystems, Foster City, California, United States). Briefly, 3 ng of cDNA was amplified by 40 cycles of a two-step PCR reaction (95 °C for 15 s denaturation and 60 °C for 1 min. annealing/elongation) containing 100 nM gene-specific primers (Table S8) and the SYBR Green PCR Master Mix (Applied Biosystems). Gene expression was normalized using human ACTB (also 汕-actin) expression levels. Results are presented as percent expression relative to uninjected clutchmate control embryos ± standard deviation of the mean.* p+ y7 g% }1 v& v" S
0 @1 I; a% e$ J
High-throughput loss-of-function genetic screen in zebrafish
1 ]9 V, y/ i9 \
The likely zebrafish orthologs of differentially expressed human genes were identified using Ensembl's gene homology prediction program (http://www.ensembl.org, build Zv3) in combination with comparison of the human protein sequence to the Institute for Genomic Research zebrafish EST database (http://www.tigr.org, release 13–15). The criteria for a likely ortholog was ≡ 40% amino acid identity over the entire length of the protein or ≡ 50% if the fish protein sequence was only partial. MO (Gene Tools, Philomath, Oregon, United States) sequences were designed complimentary to the region of translational initiation of the zebrafish orthologs in order to inhibit protein translation (Table S6) and injected into one- to two-cell zebrafish embryos as previously described [36]. MOs were injected at 1.5, 3, 4.5, and 6 ng into 50–70 embryos derived from mating Tg zebrafish containing a fluorescent DsRed reporter cassette driven by the gata1 promoter (gata1

sRed Tg) or a fluorescent EGFP reporter cassette driven by the fli1 promoter (fli1:EGFP Tg). MO-targeted (morphant) zebrafish were evaluated for defects in hematopoietic and vascular development compared to uninjected controls from the same embryo clutch at 30 and 48 hpf, using fluorescence microscopy for visualization of DsRed blood cells and EGFP vasculature.' F5 I. F$ n! R. O+ r3 z S$ ~* x' @

Whole zebrafish Q-RT-PCR confirmation
6 C* I) y/ o% o+ ^9 a& [& a3 w
Total RNA was isolated from five MO-injected zebrafish embryos at 48 hpf with hematopoietic defects, determined based on gata1

sRed hematopoietic cell production, or five uninjected clutchmate controls using the RNeasy Mini Kit (Qiagen, Valencia, California, United States) according to the manufacturer's instructions. Total RNA was incubated with DNaseI (Invitrogen) to digest contaminating genomic DNA, and reverse transcribed to generate cDNA using SuperScript II Reverse Transcriptase (Invitrogen) according to the manufacturer's protocol. Q-RT-PCR was performed using the ABI Prism 7000 Sequence Detection System (Applied Biosystems). Briefly, 1/20 of the total cDNA from five zebrafish embryos was amplified by 40 cycles of a two-step PCR reaction (95 °C for 15 s denaturation and 60 °C for 1 min. annealing/elongation) containing 100 nM gene-specific primers (Table S8) and the SYBR Green PCR Master Mix (Applied Biosystems). Gene expression was normalized using a zebrafish gapd expression level. Results are presented as percentage expression relative to uninjected clutchmate control embryos ± standard deviation of the mean.1 K8 \6 Z& }6 Z3 A) S( Z) V
: ?6 A" _0 A! a
Whole-mount in situ hybridization
/ E) h8 e/ C; L& Q5 v$ v4 w
The scl, myod, flk1, gata1, and cmyb riboprobes were generated and whole-mount in situ hybridization of zebrafish embryos was conducted as previously described [41].; I7 B. Z" {% B0 c

Overexpression of human and zebrafish Sprouty gene family members

The XhoI and KpnI fragment of the SPRY1 open reading frame (ORF) (Open Biosystems, Huntsville, Alabama, United States) was cloned into pENTR1A (Invitrogen) and subsequently transferred into a modified pFRM2.1 zebrafish expression vector using the Gateway cloning system (Invitrogen) to create the pFRM2.1–SPRY1 vector. pFRM2.1–SPRY1 was coinjected with pFRM2.1–EGFP at a 5:1 ratio into the yolk/cell interface of one-cell gata1

sRed Tg zebrafish embryos as described for MO injections. Defects in hematopoietic development of EGFP  embryos were analyzed by comparison to embryos from the same clutch injected with pFRM2.1–EGFP alone, using fluorescence microscopy to visualize DsRed  blood cells.
2 s) j, W# o# ~% d
Supporting Information# M- _1 d- z4 X8 x0 @( Y/ B" {
0 J4 |* y( G! n; C7 n
ML–ICs are highly enriched in Rholo (white bar) compared to Rhohi (black bar) cells from UCB (n = 3, p

(3.6 MB TIF).

The Rholo to Rhohi fold change for DLK1, ABCB1, BMP6, HELLS, CDC25A, MAFB, and S100A8 was determined using Affymetrix GeneChip analysis (gray bars) and Q-RT-PCR (black bars) to confirm the fidelity of the microarray results for (A) adult BM (n ≡ 2) and (B) UCB (n ≡ 3).
5 h$ X' T% H: ]: a8 H
(7.2 MB TIF).

Percentages of Gene Ontology classifications of the genes differentially expressed between Rholo and Rhohi cells from both UCB and adult BM (p  1.5 in either UCB or BM).- ]: M. ~" C; M0 U" Q

(2.7 MB TIF).

(437 KB XLS).$ l- o4 ?+ \5 ~$ w' e3 _0 \6 B

(738 KB XLS).

(75 KB XLS).2 a1 i% i* T6 C6 o6 f9 {, c
6 }/ s8 o5 a7 x7 @2 j
(17 KB XLS).) C8 M0 T# p2 G9 R
. ^* V: @" B1 V% ~' s! b0 _6 ^; B7 I2 U
(22 KB XLS).

(20 KB XLS).
3 w9 }$ N* y  _9 R% a
(19 KB XLS).
, s% t/ D9 j) U4 J# q, a8 s: Y
(16 KB XLS).& ^/ Y# ~1 @8 O' C

Accession Numbers- _8 _& g* d1 M! ^& Y( i

The National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/) Entrez Gene accession numbers for the following genes are: ABCB1 (5243), ARMCX2 (9823), BMP4 (652), C12orf2 (11228), CCR7 (1236), Ccr7 (12775), CDKN1A (1026), chd (30161), cmyb (30519), CRYGD (1421), EVI1 (2122), EZH2 (2146), FLJ14917 (84947), flk1 (also known as kdr) (58106), FOXM1 (2305), gata1 (30481), GATA2 (2624), gata2 (30480), hbae1 (30597), HDHD2 (84064), HELLS (3070), HLF (3131), HMGA2 (8091), Hoxb4 (15412), HOXB4 (3214), HSPC039 (51124), IRAK3 (11213), Irak3 (73914), KIAA1102 (22998), KLF5 (688), lcp1 (30583), LEF1 (51176), LMO2 (4005), lmo2 (30332), MAFB (9935), MGC15875 (85007), MRPS6 (64968), myod (30513), NOTCH2 (4853), PIM1 (5292), PRKCH (5583), Prkch (18755), RBPMS (11030), scl (also known as tal1) (30766), SHH (6469), SLC40A1 (30061), SLCO3A1 (28232), SNX5 (27131), SPARC (6678), Sparc (20692), spry4 (114437), SPRY1 (10252), SSBP2 (23635), SUZ12 (23512), ZFHX1B (9839), ZNF165 (7718), and ZNF331 (55422).

The microarray data in this manuscript have been deposited in the GEO database (http://www.ncbi.nlm.nih.gov/geo/, and have been assigned the accession number GSE2666.

Acknowledgments; G3 J$ n5 R  G, r2 ^  L  W$ s

We thank Brent Bill, Aubrey Nielsen, and Shridhar Sivasubbu for zebrafish injections, Greg Veltri and Joel Sederstrom for cell sorting at the University of Minnesota Flow Cytometry Facility, Dr. Leonard Zon for the gata1:DsRed Tg zebrafish line, and Drs. Stephanie Salesse and David Largaespada for critical review of the manuscript. This research was supported by Medical Scientist Training Program (GM008244) and National Cancer Institute (NCI) (CA009138) training grants to CEE, a MinnCResT training fellowship to MAP sponsored by the National Institute of Dental and Craniofacial Research (DE07288), NCI Program Project Grant (CA065493) to CMV, and National Institute of General Medical Sciences R01 (GM63904) to SCE.# q; ?: z) ?" X$ E- S
1 s: ?3 x8 ^, ^$ U- w4 c2 s  ?6 f
Author contributions. CEE, EMM, SCE, and CMV conceived and designed the experiments. CEE, EMM, CMF, and TFW performed the experiments. CEE, EMM, CMF, TFW, SMG, SCE, and CMV analyzed the data. MAP and SMG contributed reagents/materials/analysis tools. CEE, EMM, SCE, and CMV wrote the paper.
  N/ j5 a6 l3 C
References
/ o) P0 s& l2 T' p8 b
Morrison SJ, Uchida N, Weissman IL (1995) The biology of hematopoietic stem cells. Annu Rev Cell Dev Biol 11:35–71.( r, ?& ^+ E5 ~4 j, ]
# y" s3 W3 [; q, F3 X+ u
Krause DS (2002) Regulation of hematopoietic stem cell fate. Oncogene 21:3262–3269.% V7 J, y# _' m/ C

Sauvageau G, Thorsteinsdottir U, Eaves CJ, Lawrence HJ, Largman C (1995) Overexpression of HOXB4 in hematopoietic cells causes the selective expansion of more primitive populations in vitro and in vivo. Genes Dev 9:1753–1765.- `/ [! L$ v/ Z' _+ H9 c" W

Antonchuk J, Sauvageau G, Humphries RK (2001) HOXB4 overexpression mediates very rapid stem cell regeneration and competitive hematopoietic repopulation. Exp Hematol 29:1125–1134.4 H  A4 I- F5 Z2 X* ~/ Z# i0 L, M
3 g( f: l6 }4 D) X
Buske C, Feuring-Buske M, Abramovich C, Spiekermann K, Eaves CJ (2002) Deregulated expression of HOXB4 enhances the primitive growth activity of human hematopoietic cells. Blood 100:862–868.
4 x: ~- i- {3 ^/ H& E. R* @- t
Varnum-Finney B, Xu L, Brashem-Stein C, Nourigat C, Flowers D (2000) Pluripotent, cytokine-dependent, hematopoietic stem cells are immortalized by constitutive Notch1 signaling. Nat Med 6:1278–1281.

Reya T, Duncan AW, Ailles L, Domen J, Scherer DC (2003) A role for Wnt signalling in self-renewal of haematopoietic stem cells. Nature 423:409–414.
7 w# b$ _& j2 N9 E' L+ o
Murdoch B, Chadwick K, Martin M, Shojaei F, Shah KV (2003) Wnt-5A augments repopulating capacity and primitive hematopoietic development of human blood stem cells in vivo. Proc Natl Acad Sci U S A 100:3422–3427.
- U& w" s- g% N1 T4 k
Bhardwaj G, Murdoch B, Wu D, Baker DP, Williams KP (2001) Sonic hedgehog induces the proliferation of primitive human hematopoietic cells via BMP regulation. Nat Immunol 2:172–180.

Terskikh AV, Miyamoto T, Chang C, Diatchenko L, Weissman IL (2003) Gene expression analysis of purified hematopoietic stem cells and committed progenitors. Blood 102:94–101.
8 f8 s5 y4 s* b4 s& m, ^4 B; U' n
Park IK, He Y, Lin F, Laerum OD, Tian Q (2002) Differential gene expression profiling of adult murine hematopoietic stem cells. Blood 99:488–498.
7 x0 o5 y& [& H9 x3 V
Ivanova NB, Dimos JT, Schaniel C, Hackney JA, Moore KA (2002) A stem cell molecular signature. Science 298:601–604.8 ~6 ^2 L% x. l% G
' F1 v; f) @; n8 e+ \
Ramalho-Santos M, Yoon S, Matsuzaki Y, Mulligan RC, Melton DA (2002) ※Stemness§: Transcriptional profiling of embryonic and adult stem cells. Science 298:597–600.$ p/ K( B4 y$ w
! C; m; W( k+ b
Fortunel NO, Otu HH, Ng HH, Chen J, Mu X (2003) Comment on ※＆Stemness': Transcriptional profiling of embryonic and adult stem cells§ and ※a stem cell molecular signature.§. Science 302:393.

Bhatia M, Wang JC, Kapp U, Bonnet D, Dick JE (1997) Purification of primitive human hematopoietic cells capable of repopulating immune-deficient mice. Proc Natl Acad Sci U S A 94:5320–5325.+ p1 Y4 u- P) h; P8 ~, u5 q

Shojaei F, Gallacher L, Bhatia M (2004) Differential gene expression of human stem progenitor cells derived from early stages of in utero human hematopoiesis. Blood 103:2530–2540.
# [# ]' ?8 x9 a8 y; b( q# o0 a
Liu H, Verfaillie CM (2002) Myeloid–lymphoid initiating cells (ML–IC) are highly enriched in the rhodamine-c-kit( )CD33()CD38() fraction of umbilical cord CD34( ) cells. Exp Hematol 30:582–589.2 ]3 M; x8 V  Z8 V% `% c& b
1 o( E+ f+ e4 W4 t, U$ h
Cheng T, Rodrigues N, Shen H, Yang Y, Dombkowski D (2000) Hematopoietic stem cell quiescence maintained by p21cip1/waf1. Science 287:1804–1808.. Z! Z8 S$ t* u2 T- Y; @% ]7 A
+ @# _& t: d- P3 d) i& Z
Seidel MG, Look AT (2001) E2A-HLF usurps control of evolutionarily conserved survival pathways. Oncogene 20:5718–5725.

Buonamici S, Chakraborty S, Senyuk V, Nucifora G (2003) The role of EVI1 in normal and leukemic cells. Blood Cells Mol Dis 31:206–212." Y: Z* T7 u& I& X8 C: m

Warren AJ, Colledge WH, Carlton MB, Evans MJ, Smith AJ (1994) The oncogenic cysteine-rich LIM domain protein rbtn2 is essential for erythroid development. Cell 78:45–57.: `9 G# ]* E! ]( W6 k. y6 r

Zhou Y, Lim KC, Onodera K, Takahashi S, Ohta J (1998) Rescue of the embryonic lethal hematopoietic defect reveals a critical role for GATA-2 in urogenital development. EMBO J 17:6689–6700.

Travis A, Amsterdam A, Belanger C, Grosschedl R (1991) LEF-1, a gene encoding a lymphoid-specific protein with an HMG domain, regulates T-cell receptor alpha enhancer function [corrected]. Genes Dev 5:880–894.
7 B9 ~; ^4 }" I! w( ^
Bigas A, Martin DI, Milner LA (1998) Notch1 and Notch2 inhibit myeloid differentiation in response to different cytokines. Mol Cell Biol 18:2324–2333.

Hashizume H, Hamalainen H, Sun Q, Sucharczuk A, Lahesmaa R (2003) Downregulation of mafB expression in T-helper cells during early differentiation in vitro. Scand J Immunol 57:28–34.

Geiman TM, Muegge K (2000) Lsh, an SNF2/helicase family member, is required for proliferation of mature T lymphocytes. Proc Natl Acad Sci U S A 97:4772–4777.- S* H; s; E7 |3 S( j# f

Birve A, Sengupta AK, Beuchle D, Larsson J, Kennison JA (2001) Su(z)12 a novel Drosophila Polycomb group gene that is conserved in vertebrates and plants. Development 128:3371–3379.
4 A6 Z* q& m2 v4 D4 U, O3 t
Ruzinova MB, Benezra R (2003) Id proteins in development, cell cycle and cancer. Trends Cell Biol 13:410–418.3 {6 z7 A# v# f# b0 q5 u
* D: K* b; p1 m! T1 v& W' f& C) z1 @
Hoyle CF, Sherrington P, Hayhoe FG (1988) Translocation (3;6)(q21;p21) in acute myeloid leukemia with abnormal thrombopoiesis and basophilia. Cancer Genet Cytogenet 30:261–267.

Sonoki T, Harder L, Horsman DE, Karran L, Taniguchi I (2001) Cyclin D3 is a target gene of t(6;14)(p21.1;q32.3) of mature B-cell malignancies. Blood 98:2837–2844.
# \) j3 T+ a  ]
Nagarajan L, Louie E, Tsujimoto Y, ar-Rushdi A, Huebner K (1986) Localization of the human pim oncogene (PIM) to a region of chromosome 6 involved in translocations in acute leukemias. Proc Natl Acad Sci U S A 83:2556–2560.

Ekker SC (2000) Morphants: A new systematic vertebrate functional genomics approach. Yeast 17:302–306." ]5 Y! G% Z: E

Ackermann GE, Paw BH (2003) Zebrafish: A genetic model for vertebrate organogenesis and human disorders. Front Biosci 8:d1227–1253.1 E) _% L2 n/ ?

Liao W, Ho CY, Yan YL, Postlethwait J, Stainier DY (2000) Hhex and scl function in parallel to regulate early endothelial and blood differentiation in zebrafish. Development 127:4303–4313.
* s, b! A" ~4 X3 x2 ~1 {# o
Thisse C, Zon LI (2002) Organogenesis〞Heart and blood formation from the zebrafish point of view. Science 295:457–462.
  X" x8 M0 h, l. g
Nasevicius A, Ekker SC (2000) Effective targeted gene ※knockdown§ in zebrafish. Nat Genet 26:216–220.' K# ?% L+ S7 E7 B: G* d( w9 M

Galloway JL, Wingert RA, Thisse C, Thisse B, Zon LI (2005) Loss of gata1 but not gata2 converts erythropoiesis to myelopoiesis in zebrafish embryos. Dev Cell 8:109–116.  Y: _8 [6 O" m7 [& Q  `
1 C6 \& t/ h3 F- p
Haffter P, Granato M, Brand M, Mullins MC, Hammerschmidt M (1996) The identification of genes with unique and essential functions in the development of the zebrafish, Danio rerio Development 123:1–36.
7 G& K1 k& [, p+ g- U% Y
Pickart MA, Sivasubbu S, Nielsen AL, Shriram S, King RA (2004) Functional genomics tools for the analysis of zebrafish pigment. Pigment Cell Res 17:461–470.  I  C  G) g2 G, }- \& M
; p* ~4 H" s% @- X" a2 Q
Golling G, Amsterdam A, Sun Z, Antonelli M, Maldonado E (2002) Insertional mutagenesis in zebrafish rapidly identifies genes essential for early vertebrate development. Nat Genet 31:135–140.

Leung AY, Mendenhall EM, Kwan TT, Liang R, Eckfeldt C (2005) Characterization of expanded intermediate cell mass in zebrafish chordin morphant embryos. Dev Biol 277:235–254.
# h" P- r* |5 F3 e5 v, O* w
Hammerschmidt M, Pelegri F, Mullins MC, Kane DA, van Eeden FJ (1996) Dino and mercedes, two genes regulating dorsal development in the zebrafish embryo. Development 123:95–102.& j  [  M  `+ Y" i- Y) ?

Kobayashi K, Hernandez LD, Galan JE, Janeway CA Jr, Medzhitov R (2002) IRAK-M is a negative regulator of toll-like receptor signaling. Cell 110:191–202.
& t( q- z/ z% R1 V
Forster R, Schubel A, Breitfeld D, Kremmer E, Renner-Muller I (1999) CCR7 coordinates the primary immune response by establishing functional microenvironments in secondary lymphoid organs. Cell 99:23–33.! \  d7 K4 g5 l* N" N  k* ?
5 n: U( y2 w5 Q' P: o
Gilmour DT, Lyon GJ, Carlton MB, Sanes JR, Cunningham JM (1998) Mice deficient for the secreted glycoprotein SPARC/osteonectin/BM40 develop normally but show severe age-onset cataract formation and disruption of the lens. EMBO J 17:1860–1870.0 u% [, H, {  U1 n% }: L
7 m+ F' }+ C" w: m, M; [
Chida K, Hara T, Hirai T, Konishi C, Nakamura K (2003) Disruption of protein kinase Ceta results in impairment of wound healing and enhancement of tumor formation in mouse skin carcinogenesis. Cancer Res 63:2404–2408.
. o3 \- f- G% E5 P! K. N
Brun AC, Bjornsson JM, Magnusson M, Larsson N, Leveen P (2004) Hoxb4-deficient mice undergo normal hematopoietic development but exhibit a mild proliferation defect in hematopoietic stem cells. Blood 103:4126–4133.; E  ~$ O7 b9 ^& i
, {7 f6 O+ T" L1 w
Furthauer M, Reifers F, Brand M, Thisse B, Thisse C (2001) Sprouty4 acts in vivo as a feedback-induced antagonist of FGF signaling in zebrafish. Development 128:2175–2186.5 U! x4 e5 o6 g

Thompson MA, Ransom DG, Pratt SJ, MacLennan H, Kieran MW (1998) The cloche and spadetail genes differentially affect hematopoiesis and vasculogenesis. Dev Biol 197:248–269.- ^  `6 c( u& C! L( Y: s8 D: u( V

Kim HJ, Bar-Sagi D (2004) Modulation of signalling by Sprouty A developing story. Nat Rev Mol Cell Biol 5:441–450.$ I4 D2 M0 [2 l+ n" _! S9 j# d; u: b
' U! @& m4 N5 Q, X! z' p
Brekken RA, Sage EH (2001) SPARC, a matricellular protein: At the crossroads of cell-matrix communication. Matrix Biol 19:816–827.$ K& d6 T- U+ |+ g  y

Marin F, Charnay P (2000) Hindbrain patterning: FGFs regulate Krox20 and mafB/kr expression in the otic/preotic region. Development 127:4925–4935.(Craig E. Eckfeldt, Eric M)