Gene expression analysis of embryonic stem cells expressing VE-cadherin (CD144) during endothelial differentiation
© Nikolova-Krstevski et al; licensee BioMed Central Ltd. 2008
Received: 24 July 2007
Accepted: 22 May 2008
Published: 22 May 2008
Endothelial differentiation occurs during normal vascular development in the developing embryo. This process is recapitulated in the adult when endothelial progenitor cells are generated in the bone marrow and can contribute to vascular repair or angiogenesis at sites of vascular injury or ischemia. The molecular mechanisms of endothelial differentiation remain incompletely understood. Novel approaches are needed to identify the factors that regulate endothelial differentiation.
Mouse embryonic stem (ES) cells were used to further define the molecular mechanisms of endothelial differentiation. By flow cytometry a population of VEGF-R2 positive cells was identified as early as 2.5 days after differentiation of ES cells, and a subset of VEGF-R2+ cells, that were CD41 positive at 3.5 days. A separate population of VEGF-R2+ stem cells expressing the endothelial-specific marker CD144 (VE-cadherin) was also identified at this same time point. Channels lined by VE-cadherin positive cells developed within the embryoid bodies (EBs) formed by differentiating ES cells. VE-cadherin and CD41 expressing cells differentiate in close proximity to each other within the EBs, supporting the concept of a common origin for cells of hematopoietic and endothelial lineages.
Microarray analysis of >45,000 transcripts was performed on RNA obtained from cells expressing VEGF-R2+, CD41+, and CD144+ and VEGF-R2-, CD41-, and CD144-. All microarray experiments were performed in duplicate using RNA obtained from independent experiments, for each subset of cells. Expression profiling confirmed the role of several genes involved in hematopoiesis, and identified several putative genes involved in endothelial differentiation.
The isolation of CD144+ cells during ES cell differentiation from embryoid bodies provides an excellent model system and method for identifying genes that are expressed during endothelial differentiation and that are distinct from hematopoiesis.
A close temporal and spatial association exists between the development of hematopoietic and endothelial stem cells during embryogenesis. This is first observed in extra-embryonic mesodermal tissues such as the yolk sac. Within the developing yolk sac primitive erythrocytes are surrounded by a layer of angioblasts [1, 2]. The close temporal and spatial association of these two lineages led to the hypothesis that they arise from a common precursor of mesodermal origin known as the hemangioblast . Further support for the existence of a common precursor comes from gene targeting disruption studies, which demonstrate that some genes are essential for the development of both lineages. For example, the SCL/tal-1 gene, previously described for its essential role during hematopoiesis, is also expressed in the developing vasculature and is required for vascular remodeling during embryogenesis [4, 5]. One of the earliest cell surface markers that defines this bipotential cell is the VEGF receptor-2 (VEGF-R2). Targeted disruption of VEGF-R2 similarly leads to defects in vasculogenesis and hematopoiesis .
One powerful tool to study the earliest steps in the differentiation of hematopoietic and endothelial cells is the use of embryonic stem cells in culture . Most studies suggest that the molecular events leading to the development of the hematopoietic and endothelial lineages in embryoid bodies (EBs) derived from ES cells in culture are very similar to those observed during embryonic development [8, 9]. During ES cell differentiation, blast colony-forming cells (BL-CFCs) can be isolated from EBs that exhibit endothelial and hematopoietic potential . BL-CFCs expressing VEGF-R2 can be identified within EBs about 2.5 days.
Several transcription factors that are required for definitive hematopoiesis have been identified (Table 2). The transcription factor Runx1 (AML1), a member of the core binding factor (CBF) family of transcription factors is required for the establishment of definitive, but not primitive hematopoiesis . The transcription factor GATA-1 is required for differentiation of the erythroid lineage . Finally, the homeobox transcription factor HoxB4 was shown to promote the potential of hematopoietic stem cells derived from ES cells that contribute to definitive hematopoiesis, promoting the development of self-renewing, long term, hematopoietic progenitors .
Considerably less is known about the transcription factors and the molecular mechanisms that regulate endothelial differentiation. The purpose of this study was to use an ES cell differentiation system to identify genes that are preferentially expressed at early and later stages of endothelial differentiation. We identified four populations of cells; cells expressing VEGF-R2 (day 2.5), CD41 expressing cells (day 3.5), cells expressing CD144 (VE-cadherin, day 3.5), and cells expressing CD144 (day 6.5). As negative control we have also isolated VEGF-R2 negative (day 2.5), CD41 negative (day 3.5), negative CD144 (VE-Cadherin, day 3.5), and negative CD144 (day 6.5). In addition, we isolated CD144 positive and CD144 negative cells from mouse embryos at embryonic day 9.5. Microarray analysis of the RNA isolated from each of these populations of cells allowed for the identification of genes expressed in these subsets of cells and suggest their potential role during the differentiation of embryonic stem cells along the endothelial lineage.
ES cell culture
CCE Embryonic stem cells (ATCC) were maintained on irradiated primary embryonic fibroblasts (Chemicon) in knockout DMEM (Invitrogen/Gibco-BRL) supplemented with 15% fetal bovine serum (Hyclone, Logan, UT), penicillin/streptomycin 1% (Invitrogen/Gibco-BRL), l-glutamine 2 mM (Invitrogen/Gibco-BRL), non-essential amino acids 0.1 mM, nucleosides 0.1 mM, 2-mercaptoethanol 0.1 mM (Sigma, St. Louis, MO), monothioglycerol (MTG) 2 mM (Sigma), and ESGRO leukemia inhibitory factor (LIF) 1000 units/ml (Sigma). ES cells were grown in LIF-containing media for 48 hours. To generate EBs, ES cells were separated from the feeder cells, re-plated at a density of 2 × 106 cells per 10 cm dish (Fisher Scientific), and grown in media without LIF to promote differentiation.
Flow Cytometry and Cell Sorting
To prepare single-cell suspensions EBs were washed with PBS and disrupted by pipetting after 20 min treatment with cell dissociation solution at 37°C (Sigma). Whole embryo cells were dissected from the decidual tissue of E9.5 pregnant 129/Sv mice and washed in PBS (Invitrogen). Dissected embryos were incubated for 60–90 minutes at 37°C in 0.1% collagenase/dispase (Sigma) and 20% fetal bovine serum (FBS, Hyclone) and subsequently dissociated into a single-cell suspension. All cells were fixed in 3:1 ethanol:glacial acetic acid and washed extensively in PBS prior to immunofluorescent staining. Cells were initially incubated with rat-monoclonal anti-VE-cadherin antibody (1:100, Pharmingen, San Diego, CA) at room temperature for 30 min, followed by labeling with anti-rat fluorescein-isothiocynate (FITC)-conjugated antibody (1:100, Jackson Laboratories). Subsequent double labeling with phycoerythrin (PE)-conjugated anti-VEGF-R2 and -CD41 antibodies (1:100, Pharmingen) was also performed at room temperature for 30 min. Flow cytometry and cell sorting was completed on the FC5000 Flow Cytometer at the Beth Israel Deaconess Medical Center Flow Cytometry Core Facility using CXP Analysis Program.
RNA Extraction and Quantitative RT-PCR
Total RNA from FACS-sorted cells was isolated using an RNeasy kit (Qiagen). Real-time PCR was performed containing SYBR green I (1:1000 dilution, Molecular Probes), forward and reverse primers (0.8 μM each), 0.8 mM dNTPs, 2 mM MgCl2, 1U Taq Polymerase (Promega), and 1U Platinum Taq Antibody (Invitrogen) in 1xPCR Buffer under the following conditions: denaturation at 95°C (3 min); 35 cycles at 94°C (1 min), at 60°C (1 min), at 72°C (1 min). GAPDH was used as an internal reference in each reaction. Amplification was followed by melting curve analysis using the program run at the step acquisition mode to verify the presence of a single amplification product in DNA. Accumulation of PCR products was monitored and determined using the Opticon Monitor (MJ Research). The threshold cycle (CT) was determined using the Opticon analysis software. Oligonucleotide Primers for Quantitative RT-PCR for each gene were:
Ikaros: Forward GCCCTATGACAGTGCCAACT, Reverse CAGCTGGTACATGGAGCTGA
AML1: Forward TGTTGGGCATTTGACTTTGA, Reverse TTACTACCGGAGGGTTGTGG
Hemoglobin Y: Forward AAGCTCCGAGCACACCCACT, Reverse AAGCTCTGAGCACACCCACT
CD41: Forward AAGCTCTGAGCACACCCACT, Reverse CTCAGCCCTTCACTCTGACC
Endothelin Receptor B: Forward CAGGAAGAAGAGCGGTATGC, Reverse CACACCTGTGTGGATTGCTC
Thrombospondin: Forward CCAAAGCCTGCAAGAAAGAC, Reverse CCTGCTTGTTGCAAACTTGA
Tie1: Forward CAGGCACAGCAGGTTGTAGA, Reverse GTGCCACCATTTTGACACTG
Tie2: Forward AAGCATGCCCATCTGGTTAC, Reverse GTAGGTAGTGGCCACCCAGA
EphrinA1: Forward CCCACATTACGAGGACGACT, Reverse CCCAAGCTAAAAGGCCTCAA
Flt-1: Forward CCAAGGCCTCCATGAAGATA, Reverse ATACTGTCAGGGGCTGGTTG
MEF2C: Forward ACGCCTGTCACCTAACATCC, Reverse AGCTCTCAAACGCCACACTT
KLF2: Forward CCAAGAGCTCGCACCTAAAG, Reverse GTGGCACTGAAAGGGTCTGT
For transcriptional profiling, the mouse genome 430 2.0 Affymetrix GeneChip, containing more than 45,000 transcripts was used. RNA for the microarray experiments was obtained in duplicates from two separately conducted experiments using the murine ES cells. Microarray analysis was conducted via the Genomics Center at the Beth Israel Deaconess Medical Center, according to previously described protocols for total RNA extraction and purification, cDNA synthesis, in vitro transcription reaction for production of biotin-labeled cRNA, hybridization of cRNA with mouse genome 430 2.0 Affymetrix gene chips, and scanning of image output files . The quality of the scanned array images was determined on the basis of background values, percent present calls, scaling factors, and 3'-5' ratio of β-actin and GAPDH using the BioConductor R packages .
Scanned array images were analyzed by dChip, as it is more robust than MAS5.0 and RMA in signal calculation. The raw probe level data was normalized using smoothing-spline invariant set method. The signal value for each transcript was summarized using perfect-match-only (PMO) based signal modeling algorithm described in dChip. The PMO based modeling algorithm yields less number of false positives as compared to the PMO-MM model (MM meaning "mismatch-match"). In this way, the signal value corresponds to the absolute level of expression of a transcript . These normalized and modeled signal values for each transcript were used for further high-level bioinformatics analysis. During the calculation of model based expression signal values, array and probe outliers are interrogated and images spike are treated as signal outliers. The outlier detection was carried out using dChip outlier detection algorithm. A chip is considered as an outlier if the probe, single or array outlier percentage exceeds a default threshold of 5%.
When comparing two groups of samples to identify genes enriched in a given phenotype, if 90% lower confidence bound (LCB) of the fold change (FC) between the two groups was above 2.0, the corresponding gene was considered to be differentially expressed . LCB is a stringent estimate of FC and has been shown to be the better ranking statistic . It has been suggested that a criterion of selecting genes that have a LCB above 2.0 most likely corresponds to genes with an "actual" fold change of at least 3 in gene expression [16, 17].
The final list of positive, phenotype-specific, differentially expressed genes was generated by removing all the genes that were also differentially expressed in the negative control population during the same transition. For example, the initial list of differentially expressed genes for VEGF-R2+/CD41+ at day 3.5 as compared to VEGF-R2+ at day 2.5 was corrected by removing the genes that are also differentially expressed in the negative control cells (VEGF-R2-/CD41- at day 3.5 compared to VEGF-R2- at day 2.5). Further corrections in the identified differentially expressed genes were made by comparing changes in the expression in negative populations (e.g. CD144-) versus positive populations (e.g. VEGF-R2+ at day 2.5) over time.
The raw MIAME compliant microarray data is available online in the GEO repository at NCBI .
Western Blot Analysis
Whole cell lysates (50 μg) of cells were boiled in Laemmli buffer and resolved on 10% SDS-PAGE acrylamide gel. Proteins were transferred on Hybond-PVDF membranes and immunoblotted using the ECL Western Blotting Detection reagents (Amersham Biosciences, UK). Polyclonal antibodies against VE-cadherin, CD41, and alpha actin were obtained from Pharmingen.
Immunohistochemistry and microscopy of embryoid bodies and the yolk sac
EBs from different time points (day 6.5 to 10.5) were washed with PBS and embedded in OCT freezing compound (Miles Scientific, Cambridge, ML) after fixation in 4% paraformaldehyde for 1 hour at 4°C and cryoprotection in 20% (wt/vol) sucrose overnight. Frozen sections (8–10 μm thick) were mounted onto gelatin-coated glass slides, air-dried overnight, permeabilized for 3 minutes with 0.5% Triton X-100, and washed two times with PBS. Slides were then processed for immunofluorescence microscopy using FITC-conjugated CD41, PE-conjugated VE-cadherin, and FITC-vWF rat antibodies (1:100, Pharmingen). 129/sv mice were bred to generate day E8.5 embryos. The yolk sac was dissected from the embryo proper and fixed in 4% paraformaldehyde, and frozen in OCT. Immunohistochemistry was performed using 10 μm sections. Fluorescent images were obtained using a Nikon Eclispse E800 microscope (Nikon, Tokyo, Japan) equipped with a 40 X/0.75 numeric aperture objective lens as well as a Zeiss AxioCam digital camera (Carl Zeiss International, Heidelberg, Germany). Images were analyzed using Openlab imaging software version 3.0.4 (Improvision, Lexington, MA).
VE-cadherin is expressed on the surface of endothelial cells, during vascular development and in the adult cells . To verify the expression of these cell surface markers at the protein level, western blot analysis was performed at different time points during ES cell differentiation. No expression of CD41 or CD144 was observed at 2.5 days (Figure 1D). However, at day 3.5 there was a significant increase in the expression of both CD144 and CD41.
Early VEGF-R2 positive ES cells express several mesodermal markers
Mesodermal Markers Up-regulated in VEGF-R2 positive cells at day 2.5 of ES Cell Differentiation*
Gene Bank accession
Lim homeobox protein 1
Even skipped homeotic gene
Comparison of CD41-postive (day 3.5) to VEGF-R2 positive cells (day 2.5)*
Gene Bank Accession number
Krueppel like factor-1
Ikaros transcription factor
Helios transcription factor
Comparison of CD144-postive (day 3.5) to VEGF-R2 positive cells (day 2.5)*
Gene Bank Accession number
Endothelin receptor type B
Transcription factor EC
Krueppel like factor 2
Forkhead box P1
Comparison of CD144-postive (day 6.5) to VEGF-R2 positive cells (day 2.5)*
Endothelin receptor type B
CD 109 antigen
Comparison of embryonic CD144-postive to CD144-negative cells (day 9.5)*
Von Willebrands Factor
Endothelial nitric oxide synthase
Forkhead box 01
Genes Enriched in all CD144-positive cells
Fold Change (vs. CD144-E9.5)
(d 6.5 vs HB)
(d 3.5 vs HB)
Genes that showed increased expression in CD144-positive cells from EBs (day 6.5) and in the embryo (day 9.5) compared to CD41-positive (day 3.5) and CD144-negative cells (day 6.5) from Ebs
(EC vs HB)
(EC vs CD144-cells)
Markers of Hematopoiesis in the CD41 expressing ES cells at day 3.5
We identified a small population of VEGF-R2 positive cells expressing the early hematopoietic marker CD41 and cells not expressing CD41 (negative control) at day 3.5. Microarray analysis was performed on RNA derived from these cells. VEGFR2+/CD41+ compared to the VEGF-R2 positive cells at day 2.5 have 565 differentially expressed transcripts. Out of these 565 transcripts, all the transcripts that were differentially expressed in the negative control cells in the same transition (VEGF-R2-/CD41- cells at day 3.5 as compared to VEGF-R2- cells at day 2.5 or VEGF-R2-/CD41- cells at day 3.5 as compared to VEGF-R2+ cells at day 2.5) were filtered out. The final list of VEGFR2+/CD41+ specific differentially expressed transcripts consists of genes associated with hematopoiesis (Table 2). It has also been suggested that in vitro ES cell differentiation systems closely recapitulates the events associated with hematopoietic development in the yolk sac, including primitive erythropoiesis, megakaryocyte, and mast cell differentiation . The first of these to occur is primitive erythropoiesis. We identified transcription factors upstream of erythropoiesis, including AML-1, GATA-1, and erythroid Krueppel-like factor (KLF1), and genes predominantly associated with the erythroid lineage (Table 2). However, Ikaros and Helios, two zinc finger transcription factors known to regulate lymphopoiesis, were also enriched in this CD41-positive population [22, 23]. Helios expression is restricted to cells of T cell origin. These results suggest that either these transcription factors exhibit additional functions at earlier stages of hematopoiesis or that initiation of other hematopoietic lineages, in addition to the erythroid lineage, occurs during ES cell differentiation.
Markers of endothelial differentiation expressed in CD144 positive cells
To evaluate whether genes of the endothelial lineage are expressed in early CD144 (VE-cadherin) positive cells during ES cell differentiation, we compared microarray data of RNA isolated from CD144 positive cells at day 3.5 to VEGF-R2 positive cells at day 2.5. To identify genes that are uniformly upregulated in the CD144 positive cells, we also performed microarray analysis of RNA samples from the negative controls (i.e CD144 negative cells at day 3.5 to VEGF-R2 negative cells at day 2.5). The CD144 positive cells at day 3.5 have 483 differentially expressed transcripts that are further corrected by removing the transcripts differentially expressed in the negative control cells in the same transition (CD144- day 3.5 vs. VEGF-R2- day 2.5) or (CD144- day 3.5 vs. VEGF-R2+ day 2.5). This final list has enrichment of selected endothelial-specific genes, including Tie1 and the endothelin receptor B (Table 3). In addition, several transcription factors were identified as being enriched at this stage, including Krueppel-like factor 2 (KLF2), the SRY-box (Sox) 11, and the Forkhead box P1 transcription factor.
To examine the specific gene expression profile at later stages of endothelial differentiation, we performed microarray analysis of CD144 positive and CD144 negative cells isolated from ES cells 6.5 days after differentiation. The CD144 positive cells at day 6.5 compared to VEGF-R2 positive cells obtained at day 2.5 showed differential expression of 1068 transcripts. Transcripts that were differentially expressed in the negative control cells (CD144- day 6.5 vs VEGF-R2- day 2.5 and CD144- day 6.5 vs VEGF-R2+ day 2.5) were removed from the list. Additional analysis was performed to eliminate remaining false positive genes (see methods). The final list has several endothelial-specific genes including PECAM-1, claudin 5, and roundabout homolog 4 (Robo 4) (Table 4).
To identify genes that are selectively enriched in endothelial cells during embryonic vascular development, we also isolated CD144 positive cells from murine embryos at day 9.5. CD144 positive cells were sorted and separated from CD144 negative cells, and microarray analysis was performed. Because the number of CD144 positive cells isolated from the embryos was much smaller, an amplification step was performed prior to the microarray analysis. The results of the microarray analysis also demonstrate an enrichment of several endothelial-specific genes including ICAM-2, VE-cadherin, PECAM-1, and endothelial nitric oxide synthase (eNOS) (Table 5). The expression of several transcription factors was also determined to be enriched in this population of cells, including SOX 17, the Ets factor Erg and Elk3, Vezf1, Foxo1, and C/EBP delta. Table 6 shows a condensed list of genes whose expression is increased in all CD144 positive cells, including the ones isolated at different stages of EB differentiation (day 3.5 and 6.5), and day 9.5 embryos. The genes upregulated in the CD144 positive cells were not detected in hematopoietic cells (CD41 positive cells at day 3.5) (Table 7).
To further define the genes that are upregulated at the angioblast stage (day 3.5 CD144 positive), and later stages of endothelial differentiation in the EBs (day 6.5 CD144 positive), and in CD144 positive cells of the embryo (day 9.5), we performed a comparative analysis of these genes (on-line Additional file 1). We identified 41 genes (section VI) that are enriched in all three populations, 139 genes that are enriched in CD144 positive cells of the EB at early and later stages (Section IV), and 119 genes (Section V) that are enriched in CD144 positive cells in later stages of the EB (day 6.5) and in CD144 positive cells of the embryo at day E9.5. Some of the genes identified in the overlapping sections of the Venn diagram are listed in Table 6 and 7. A complete list of the up and down-regulated genes are provided in an on-line supplement .
In addition to identifying genes that are enriched at various stages of differentiation, we were also interested in determining whether the genes identified have previously been shown to be linked to each other, or to other genes with regard to specific pathways. In the past few years a number of bioinformatic programs that are designed to evaluate potential links between the genes identified and potential pathways have been developed. We used the program PathwayAssist (Ariadne Genomics Inc.) to evaluate the potential of underlying existing pathways at various stages of endothelial and hematopoietic development. In this approach a natural language processing module, MedScan, is used to describe events of regulation, interaction and modification between proteins, cell processes and small molecules and the results are generated in the forms of biological pathways, gene regulation networks and protein interaction maps . Because considerably more is known regarding the genes involved in hematopoiesis, we first performed analysis of the genes enriched in the CD41 positive cells. The initial analysis was performed using transcription factors identified in the CD41 positive cells at day 3.5 in the differentiating EBs compared to the VEGF-R2 positive cells at day 2.5. When the transcription factors alone are entered into this pathway analysis, the specific factors are linked in multiple ways. For example, KLF1 and GATA1, which are known to be involved in erythropoiesis cluster together (Additional file 2) with links to references for pathways). Similarly, Ikaros, Aiolos, and Helios, zinc finger transcription factors, are involved in lymphopoiesis are linked together. AML1 (Runx1) is known to be upstream of all of the hematopoietic lineages. The pathway analysis was further expanded to include other proteins that are enriched in the CD41 positive population, including the hemoglobin genes, myeloperoxidase, c-kit, glycophorin, and other genes, in addition to the transcription factors identified. Although, significantly more complex, the pathway analysis similarly supports the linkage of these transcription factors and pathways with specific target genes along various hematopoietic lineages (on-line Additional file 3). In addition to demonstrating the utility of this kind of pathway analysis for linking transcription factors with known gene targets and pathways that have previously been defined for hematopoiesis, these results also support the potential of this approach for identifying similar pathways during the process of endothelial differentiation.
We used a similar approach to try to identify the underlying molecular mechanisms involved at various stages of endothelial differentiation. We first evaluated the genes that are upregulated in CD144 positive cells at day 3.5 (angioblast) compared to the hemangioblast. The transcription factors that are highly enriched in this population include KLF2, Sox11 and the forkhead box P1 transcription factor. Some of the potential networks and pathways were also evaluated using the PathwayAssist program (Additional files 4 and 5). Similar analyses were conducted using the Pathway Assist program for genes identified to be enriched in CD144 positive cells (ES cell day 6.5) compared to day 2.5 VEGF-R2 positive ES cells (Additional files 6 and 7) and for CD144 positive cells isolated from mouse embryos at day 9.5, which is comparable to day 6.5 in the ES cell differentiation system, compared to CD144 negative cells (Additional files 8 and 9). We recognize that there are limitations of this type of analysis, and that the results of this analysis are only suggestive of potential pathways among the identified genes and transcription factors, and these analyses are therefore presented as supplemental data. Future studies will be needed to verify specific aspects of the putative pathways identified.
Overlap between hematopoietic and endothelial lineages
The close association of endothelial and hematopoietic lineage has now been demonstrated in several studies in a wide variety of species including xenopus, mice, and humans. Our studies similarly support a very close temporal and spatial expression pattern of endothelial and hematopoietic progenitor cells in differentiating ES cells. In human embryonic tissues, a population of endothelial-like cells with endothelial and hematopoietic potential expressing VEGF-R2, VE-cadherin, but not CD45, was similarly identified . We observed weak VE-cadherin expression on CD41 positive cells in both the yolk sac and differentiating EBs, but did not observe CD41 expression on endothelial cells lining the vascular channels in the EB or yolk sac, suggesting that as CD41 expression occurs the cells bud off of the developing vasculature.
Differentiation of the hemangioblast toward the hematopoietic lineage
Several transcription factors have been shown to be critical for different stages of hematopoietic cell differentiation. Some of the critical transcription factors known to play a role early in hematopoietic cell differentiation include Runx-1 (AML-1), GATA-1, and members of the Ikaros zinc finger transcription factor family. We also observed enrichment for these factors in early CD41 positive cells within differentiating EBs based on our microarray data. Although the EB is generally considered as a model of primitive hematopoiesis, similar to what is observed in the yolk sac, the expression of many Ikaros family members that are known to regulate several different hematopoietic lineages suggest a possible broader role for the CD41 positive cells identified in developing EBs in several hematopoietic lineages .
Differentiation of the hemangioblast toward the endothelial lineage
Considerably less is known about the transcription factors and molecular mechanisms involved in regulating the differentiation of the hemangioblast along the endothelial lineage. In early (day 3.5) CD144 positive cells we observed that the increased expression of genes involved in regulating endothelial function such as Tie1 and the endothelin receptor B, compared to genes expressed in VEGF-R2 positive hemangioblast cells (day 2.5). We similarly identified a select number of transcription factors that were enriched in these cells, including KLF2, Sox11, and forkhead box (Fox) P1. Targeted disruption of FOXP1, which is expressed in the endocardium of the heart, is associated with abnormalities in cardiac development, similar to that observed with NFATc1 . The Krueppel-like factor (KLF) family members have also been shown to be involved in cellular differentiation. KLF2 is expressed in the endothelium and is upregulated in response to fluid shear stress [29, 30].
Because the early CD144 positive cells may not represent fully differentiated endothelial cells, microarray analysis was also performed on CD144 positive cells at a later stage of differentiation in the EBs (day 6.5) and isolated from the embryo (day 9.5). In addition to increased expression of several Sox family members (11 and 17), the expression of several other transcription factors including the nuclear hormone receptor COUP-TFII (NR2F1), MEF2C, and KLF2, were observed. Targeted disruption of the orphan nuclear receptor COUP-TFII results in abnormal vascular and cardiac development . One of the main downstream targets of COUP-TFII is the angiogenic factor Angiopoietin-1. COUP family members are involved in regulating the expression of BMP-4 during embryogenesis. BMP-4 promotes hematopoiesis and angiogenesis during vascular development .
In CD144 expressing cells of the embryo (day 9.5) we observed significant overlap in the endothelial-specific genes and transcription factors upregulated in the CD144 positive cells derived from the embryo and differentiating EBs at day 6.5. The transcription factors that were highly expressed and enriched in these cells included, Sox 7, Sox 17, MEF2C, and KLF2. In addition several other transcription factors were also identified in CD144 positive cells of the embryo including the ETS transcription factor family members ERG and ELK3, GATA2, NFATc1, Sox18, Forkhead O1A, and C/EBP delta. The Ets transcription factor ELK3 is highly expressed at sites of vasculogenesis and angiogenesis during early mouse development . Furthermore, down-regulation of ELK3 is associated with marked reductions in VEGF expression. Although ELK3 principally functions as a repressor, Ras activation leads to a phosphorylated form of ELK3 that is highly active . The Ets factor ERG is also enriched in endothelial cells and down regulation of ERG in endothelial cells is associated with reductions in the expression of the endothelial-specific genes . The GATA factors have similarly been shown to regulate endothelial-specific genes. GATA2 is involved in the regulation of PECAM-1, endothelin-1, and ICAM-2 [35–37]. Both the forkhead family members and NFATc1 function as transcriptional mediators of VEGF in endothelial cells [38, 39]. There are no previous reports suggesting a role for C/EBP delta in endothelial function or gene expression.
Microarray Studies of Hematopoietic (HP) and Endothelial Cell (EC) Differentiation in EBs
1. Chen D et al
2. Woll PS et al
CD31, CD34, Flk-1
3. Schenke K et al
Col IV, FGF4
4. Lu SJ et al
5. Wang H et al
6. Williamson AJ et al
7. Rolny C et al
8. Wang C. et al
9. Wu Y. et al
10. Ferguson JE et al
11. Ng YS et al
During the hemangioblast phase a number of early hematopoietic and endothelial markers are found to be co-expressed on the VEGFR-2+ cells from day 3–6 of differentiation. Genes found in these group include CD34, TAL/SCL, Runx1/AML1, GATA-2/3, Lmo2, Notch1, c-kit, CD41, HoxB4 (all of them referred to as hematopoietic commitment markers), and genes thought to be expressed on both, hematopoietic and endothelial progenitors including CD31, VE-cadherin, CXCRF, Flt-1 [[37, 38] and ], KDR , TGFα-1, Bmp7, Smad1, Wnt associated genes FrzB, β-catenin, N-cadherin, Msx2, Cyclin D1, Dkk1 , lysocardiolipin acyltransferase (Lycat)  and finally PRDM6, encoding a Krueppel-like zinc finger protein (KLF) that functions in cell-fate decisions and malignant transformation . Most of these genes are downregulated as the divergence of the two cell lines takes place, however, a number of these markers maintain steady levels of expression throughout the course of differentiation. The onset of the definitive hematopoiesis is marked by the expression of CD45, the genes involved in adult hematopoietic cell renewal (PIK3R1, MDR1, IRS1, RGS18, SUMO-1 and Wnt5) , CXCL12, Alcam , mayeloblastosis oncogene (Myb) and Hemoglobin genes .
Genes whose expression is enriched in more mature endothelial cells from day 6 and onwards include VE-cadherin, CXCR4 , Angiopoietin 2 (Ang2), Endothelin 1 (End1) , ICAM-2, neuropilin-1, endoglin, transforming growth factor genes (TGFβ-2, BMP4 and BMPR-1), Wnt associated genes (Wnt2, Wnt5a, Nkd1, cyclin D2, Frzd7, Lef1, Frzb, Sox17, Msx2, fibrinogen, Myc) , Tie-2, genes involved in angiogenic response to VEGF such as PDGFBB, PDGFR-1, Shb , ASB4 (which functions in vascular differentiation in an oxygen-dependent manner) , caveolin, Hey1, FGFR-1, and adhesion protein genes cadherin 5 and claudin 3 . The findings from our microarray studies are in accord with the previously published studies. In addition, we have identified several additional genes such as Ikaros and Helios zinc-finger transcription factors, Erg, claudin 5, and Sox17, that segregate with the hematopoietic or endothelial cell populations (see Table 1,2,3,4,5,6,7). Future studies will be needed to identify whether or not these genes play a significant role as determinants of endothelial or hematopoietic cell differentiation.
In conclusion, the results of our study support the use of embryonic stem cells to study the molecular events underlying hematopoiesis and endothelial differentiation. Our results demonstrate the expression of distinct sets of genes in populations of embryonic stem cells expressing CD41 or CD144 that are consistent with genes associated with the hematopoietic or endothelial lineages. Future studies should be directed at defining the role of selected transcription factors that are enriched in selected subsets of stem cells during the process of differentiation.
This work was supported by NIH grants HL-67219 (PO) and P01 HL76540 (PO), and AHA award EIA0740012 (PO).
- Kinder SJ, Tsang TE, Quinlan GA, Hadjantonakis AK, Nagy A, Tam PP: The orderly allocation of mesodermal cells to the extraembryonic structures and the anteroposterior axis during gastrulation of the mouse embryo. Development. 1999, 126 (21): 4691-4701.Google Scholar
- Moore MA, Metcalf D: Ontogeny of the haemopoietic system: yolk sac origin of in vivo and in vitro colony forming cells in the developing mouse embryo. Br J Haematol. 1970, 18 (3): 279-296. 10.1111/j.1365-2141.1970.tb01443.x.View ArticleGoogle Scholar
- Huber TL, Kouskoff V, Fehling HJ, Palis J, Keller G: Haemangioblast commitment is initiated in the primitive streak of the mouse embryo. Nature. 2004, 432 (7017): 625-630. 10.1038/nature03122.View ArticleGoogle Scholar
- Visvader JE, Fujiwara Y, Orkin SH: Unsuspected role for the T-cell leukemia protein SCL/tal-1 in vascular development. Genes Dev. 1998, 12 (4): 473-479. 10.1101/gad.12.4.473.View ArticleGoogle Scholar
- Endoh M, Ogawa M, Orkin S, Nishikawa S: SCL/tal-1-dependent process determines a competence to select the definitive hematopoietic lineage prior to endothelial differentiation. Embo J. 2002, 21 (24): 6700-6708. 10.1093/emboj/cdf674.View ArticleGoogle Scholar
- Shalaby F, Rossant J, Yamaguchi TP, Gertsenstein M, Wu XF, Breitman ML, Schuh AC: Failure of blood-island formation and vasculogenesis in Flk-1-deficient mice. Nature. 1995, 376 (6535): 62-66. 10.1038/376062a0.View ArticleGoogle Scholar
- Fraser ST, Ogawa M, Nishikawa S: Embryonic stem cell differentiation as a model to study hematopoietic and endothelial cell development. Methods Mol Biol. 2002, 185: 71-81.Google Scholar
- Keller G, Kennedy M, Papayannopoulou T, Wiles MV: Hematopoietic commitment during embryonic stem cell differentiation in culture. Mol Cell Biol. 1993, 13 (1): 473-486.View ArticleGoogle Scholar
- Fraser ST, Ogawa M, Yokomizo T, Ito Y, Nishikawa S: Putative intermediate precursor between hematogenic endothelial cells and blood cells in the developing embryo. Dev Growth Differ. 2003, 45 (1): 63-75. 10.1046/j.1440-169X.2003.00675.x.View ArticleGoogle Scholar
- Wang Q, Stacy T, Binder M, Marin-Padilla M, Sharpe AH, Speck NA: Disruption of the Cbfa2 gene causes necrosis and hemorrhaging in the central nervous system and blocks definitive hematopoiesis. Proc Natl Acad Sci U S A. 1996, 93 (8): 3444-3449. 10.1073/pnas.93.8.3444.View ArticleGoogle Scholar
- Pevny L, Simon MC, Robertson E, Klein WH, Tsai SF, D'Agati V, Orkin SH, Costantini F: Erythroid differentiation in chimaeric mice blocked by a targeted mutation in the gene for transcription factor GATA-1. Nature. 1991, 349 (6306): 257-260. 10.1038/349257a0.View ArticleGoogle Scholar
- Kyba M, Perlingeiro RC, Daley GQ: HoxB4 confers definitive lymphoid-myeloid engraftment potential on embryonic stem cell and yolk sac hematopoietic progenitors. Cell. 2002, 109 (1): 29-37. 10.1016/S0092-8674(02)00680-3.View ArticleGoogle Scholar
- Jones J, Otu H, Spentzos D, Kolia S, Inan M, Beecken WD, Fellbaum C, Gu X, Joseph M, Pantuck AJ, Jonas D, Libermann TA: Gene Signatures of Progression and Metastasis in Renal Cell Cancer. Clin Cancer Res. 2005, 11: 5730-5739. 10.1158/1078-0432.CCR-04-2225.View ArticleGoogle Scholar
- Ihaka RG: A language for data analysis and graphics. J Comp Grap[h Stat. 1996, 5 (3): 299-314. 10.2307/1390807.Google Scholar
- Li C, Wong WH: Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci U S A. 2001, 98 (1): 31-36. 10.1073/pnas.011404098.View ArticleGoogle Scholar
- Ramalho-Santos M, Yoon S, Matsuzaki Y, Mulligan RC, Melton DA: "Stemness": transcriptional profiling of embryonic and adult stem cells. Science. 2002, 298 (5593): 597-600. 10.1126/science.1072530.View ArticleGoogle Scholar
- Yuen T, Wurmbach E, Pfeffer RL, Ebersole BJ, Sealfon SC: Accuracy and calibration of commercial oligonucleotide and custom cDNA microarrays. Nucleic Acids Res. 2002, 30 (10): e48-10.1093/nar/30.10.e48.View ArticleGoogle Scholar
- GEO repository. [http://www.ncbi.nlm.nih.gov/geo]
- Mikkola HK, Fujiwara Y, Schlaeger TM, Traver D, Orkin SH: Expression of CD41 marks the initiation of definitive hematopoiesis in the mouse embryo. Blood. 2003, 101 (2): 508-516. 10.1182/blood-2002-06-1699.View ArticleGoogle Scholar
- Lampugnani MG, Resnati M, Raiteri M, Pigott R, Pisacane A, Houen G, Ruco LP, Dejana E: A novel endothelial-specific membrane protein is a marker of cell-cell contacts. J Cell Biol. 1992, 118 (6): 1511-1522. 10.1083/jcb.118.6.1511.View ArticleGoogle Scholar
- Keller GM: In vitro differentiation of embryonic stem cells. Curr Opin Cell Biol. 1995, 7 (6): 862-869. 10.1016/0955-0674(95)80071-9.View ArticleGoogle Scholar
- Hahm K, Cobb BS, McCarty AS, Brown KE, Klug CA, Lee R, Akashi K, Weissman IL, Fisher AG, Smale ST: Helios, a T cell-restricted Ikaros family member that quantitatively associates with Ikaros at centromeric heterochromatin. Genes Dev. 1998, 12 (6): 782-796. 10.1101/gad.12.6.782.View ArticleGoogle Scholar
- Koipally J, Georgopoulos K: A molecular dissection of the repression circuitry of Ikaros. J Biol Chem. 2002, 277 (31): 27697-27705. 10.1074/jbc.M201694200.View ArticleGoogle Scholar
- Additional files. [http://www.bidmcgenomics.org/EndothelialDifferentiation/index.html]
- Nikitin A, Egorov S, Daraselia N, Mazo I: Pathway studio--the analysis and navigation of molecular networks. Bioinformatics. 2003, 19 (16): 2155-2157. 10.1093/bioinformatics/btg290.View ArticleGoogle Scholar
- Wang L, Li L, Shojaei F, Levac K, Cerdan C, Menendez P, Martin T, Rouleau A, Bhatia M: Endothelial and hematopoietic cell fate of human embryonic stem cells originates from primitive endothelium with hemangioblastic properties. Immunity. 2004, 21 (1): 31-41. 10.1016/j.immuni.2004.06.006.View ArticleGoogle Scholar
- Wang JH, Nichogiannopoulou A, Wu L, Sun L, Sharpe AH, Bigby M, Georgopoulos K: Selective defects in the development of the fetal and adult lymphoid system in mice with an Ikaros null mutation. Immunity. 1996, 5 (6): 537-549. 10.1016/S1074-7613(00)80269-1.View ArticleGoogle Scholar
- Wang B, Weidenfeld J, Lu MM, Maika S, Kuziel WA, Morrisey EE, Tucker PW: Foxp1 regulates cardiac outflow tract, endocardial cushion morphogenesis and myocyte proliferation and maturation. Development. 2004, 131 (18): 4477-4487. 10.1242/dev.01287.View ArticleGoogle Scholar
- Huddleson JP, Srinivasan S, Ahmad N, Lingrel JB: Fluid shear stress induces endothelial KLF2 gene expression through a defined promoter region. Biol Chem. 2004, 385 (8): 723-729. 10.1515/BC.2004.088.View ArticleGoogle Scholar
- Dekker RJ, van Soest S, Fontijn RD, Salamanca S, de Groot PG, VanBavel E, Pannekoek H, Horrevoets AJ: Prolonged fluid shear stress induces a distinct set of endothelial cell genes, most specifically lung Kruppel-like factor (KLF2). Blood. 2002, 100 (5): 1689-1698. 10.1182/blood-2002-01-0046.View ArticleGoogle Scholar
- Pereira FA, Qiu Y, Zhou G, Tsai MJ, Tsai SY: The orphan nuclear receptor COUP-TFII is required for angiogenesis and heart development. Genes Dev. 1999, 13 (8): 1037-1049. 10.1101/gad.13.8.1037.View ArticleGoogle Scholar
- Li F, Lu S, Vida L, Thomson JA, Honig GR: Bone morphogenetic protein 4 induces efficient hematopoietic differentiation of rhesus monkey embryonic stem cells in vitro. Blood. 2001, 98 (2): 335-342. 10.1182/blood.V98.2.335.View ArticleGoogle Scholar
- Zheng H, Wasylyk C, Ayadi A, Abecassis J, Schalken JA, Rogatsch H, Wernert N, Maira SM, Multon MC, Wasylyk B: The transcription factor Net regulates the angiogenic switch. Genes Dev. 2003, 17 (18): 2283-2297. 10.1101/gad.272503.View ArticleGoogle Scholar
- McLaughlin F, Ludbrook VJ, Cox J, von Carlowitz I, Brown S, Randi AM: Combined genomic and antisense analysis reveals that the transcription factor Erg is implicated in endothelial cell differentiation. Blood. 2001, 98 (12): 3332-3339. 10.1182/blood.V98.12.3332.View ArticleGoogle Scholar
- Cowan PJ, Tsang D, Pedic CM, Abbott LR, Shinkel TA, d'Apice AJ, Pearse MJ: The human ICAM-2 promoter is endothelial cell-specific in vitro and in vivo and contains critical Sp1 and GATA binding sites. J Biol Chem. 1998, 273 (19): 11737-11744. 10.1074/jbc.273.19.11737.View ArticleGoogle Scholar
- Kawana M, Lee ME, Quertermous EE, Quertermous T: Cooperative interaction of GATA-2 and AP1 regulates transcription of the endothelin-1 gene. Mol Cell Biol. 1995, 15 (8): 4225-4231.View ArticleGoogle Scholar
- Gumina RJ, Kirschbaum NE, Piotrowski K, Newman PJ: Characterization of the human platelet/endothelial cell adhesion molecule-1 promoter: identification of a GATA-2 binding element required for optimal transcriptional activity. Blood. 1997, 89 (4): 1260-1269.Google Scholar
- Minami T, Horiuchi K, Miura M, Abid MR, Takabe W, Noguchi N, Kohro T, Ge X, Aburatani H, Hamakubo T, Kodama T, Aird WC: Vascular endothelial growth factor- and thrombin-induced termination factor, Down syndrome critical region-1, attenuates endothelial cell proliferation and angiogenesis. J Biol Chem. 2004, 279 (48): 50537-50554. 10.1074/jbc.M406454200.View ArticleGoogle Scholar
- Abid MR, Guo S, Minami T, Spokes KC, Ueki K, Skurk C, Walsh K, Aird WC: Vascular endothelial growth factor activates PI3K/Akt/forkhead signaling in endothelial cells. Arterioscler Thromb Vasc Biol. 2004, 24 (2): 294-300. 10.1161/01.ATV.0000110502.10593.06.View ArticleGoogle Scholar
- Chen D, Wang P, Lewis RL, Daigh CA, Ho C, Chen X, Thomson JA, Kendziorski C: A microarray analysis of the emergence of embryonic definitive hematopoiesis. Exp Hematol. 2007, 35 (9): 1344-1357. 10.1016/j.exphem.2007.06.004.View ArticleGoogle Scholar
- Woll PS, Morris JK, Painschab MS, Marcus RK, Kohn AD, Biechele TL, Moon RT, Kaufman DS: Wnt signaling promotes hemato-endothelial cell development from human embryonic stem cells. Blood. 2007Google Scholar
- Schenke-Layland K, Angelis E, Rhodes KE, Heydarkhan-Hagvall S, Mikkola HK, Maclellan WR: Collagen IV induces trophoectoderm differentiation of mouse embryonic stem cells. Stem Cells. 2007, 25 (6): 1529-1538. 10.1634/stemcells.2006-0729.View ArticleGoogle Scholar
- Lu SJ, Hipp JJ, Feng QQ, Lanza RR, Atala AA: Genechip analysis of human embryonic stem cell differentiation into hemangioblasts: An in silico dissection of mixed phenotype. Genome Biol. 2007, 8 (11): R240-10.1186/gb-2007-8-11-r240.View ArticleGoogle Scholar
- Wang H, Gilner JB, Bautch VL, Wang DZ, Wainwright BJ, Kirby SL, Patterson C: Wnt2 coordinates the commitment of mesoderm to hematopoietic, endothelial, and cardiac lineages in embryoid bodies. J Biol Chem. 2007, 282 (1): 782-791. 10.1074/jbc.M606610200.View ArticleGoogle Scholar
- Williamson AJ, Smith DL, Blinco D, Unwin RD, Pearson S, Wilson C, Miller C, Lancashire L, Lacaud G, Kouskoff V, Whetton AD: Quantitative proteomic analysis demonstrates post-transcriptional regulation of embryonic stem cell differentiation to hematopoiesis. Mol Cell Proteomics. 2007Google Scholar
- Rolny C, Lu L, Agren N, Nilsson I, Roe C, Webb GC, Welsh M: Shb promotes blood vessel formation in embryoid bodies by augmenting vascular endothelial growth factor receptor-2 and platelet-derived growth factor receptor-beta signaling. Exp Cell Res. 2005, 308 (2): 381-393. 10.1016/j.yexcr.2005.04.020.View ArticleGoogle Scholar
- Wang H, Charles PC, Wu Y, Ren R, Pi X, Moser M, Barshishat-Kupper M, Rubin JS, Perou C, Bautch V, Patterson C: Gene expression profile signatures indicate a role for Wnt signaling in endothelial commitment from embryonic stem cells. Circ Res. 2006, 98 (10): 1331-1339. 10.1161/01.RES.0000220650.26555.1d.View ArticleGoogle Scholar
- Wu Y, Ferguson JE, Wang H, Kelley R, Ren R, McDonough H, Meeker J, Charles PC, Patterson C: PRDM6 is enriched in vascular precursors during development and inhibits endothelial cell proliferation, survival, and differentiation. J Mol Cell Cardiol. 2007Google Scholar
- Ferguson JE, Wu Y, Smith K, Charles P, Powers K, Wang H, Patterson C: ASB4 is a hydroxylation substrate of FIH and promotes vascular differentiation via an oxygen-dependent mechanism. Mol Cell Biol. 2007, 27 (18): 6407-6419. 10.1128/MCB.00511-07.View ArticleGoogle Scholar
- Ng YS, Ramsauer M, Loureiro RM, D'Amore PA: Identification of genes involved in VEGF-mediated vascular morphogenesis using embryonic stem cell-derived cystic embryoid bodies. Lab Invest. 2004, 84 (9): 1209-1218. 10.1038/labinvest.3700150.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.