hPSCs differentiation into CD144+hematoendothelial progenitors
hPSC-EBs were generated from one hiPSC line, A29, and two hESC lines, SA01 and H1, following an adapted version of the Lu et al. protocol (47). Briefly, mesoderm was induced during the first 48 hours of hPSC-EB formation followed by hematoendothelial specification step from the mesoderm during the next 36 hours in a serum-free medium (Fig. 1A). A progressively loss of pluripotency during hPSC-EB formation was confirmed by a decreased in the mRNA levels of two pluripotency markers, POU5F1 (OCT3/4) and NANOG, for all hPSCs (Fig. 1B-C). This loss of pluripotency coincided in time with a peak in the expression of the mesodermal marker TBXT at 48h (Fig. 1D) followed by a significant increase in the expression of KDR (CD309) at 84h of hPSC-EB formation (Fig. 1E). No significant differences in this kinetics were observed among the hPSC lines.
Since KDR (CD309) has been described as an important marker of endothelial and hematopoietic precursors (48), we hypothesized that the increased in KDR expression may mark the apparition of hematoendothelial progenitors within 84h-hPSC-EBs. Thus, we analyzed the expression of other endothelial and hematopoietic markers in 84h-hPSC-EB by flow cytometry (Fig. 1F). Between 15–25% of cells within 84h-hPSC-EBs expressed hematoendothelial markers as CD309, CD34 and CD143, and the endothelial marker CD144, independently on the hPSC line. Fewer cells (≈ 10%) expressed the endothelial marker CD31 and almost none expressed the hematopoietic markers CD43, CD45 and CD41. Interestingly, the CD309+ population did not contain exclusively hematoendothelial progenitors since only 60% of these CD309+ cells coexpressed CD34, CD143 and CD144 (supplementary Fig. 1). Instead, the CD144+ population seemed a more homogeneous population composed by more than 90% of CD34+CD143+CD309+ cells (Fig. 1G). CD144 is a classical endothelial marker expressed by both hemogenic and non-hemogenic endothelium. Thus, the CD144+ population was isolated from 84h-hPSC-EBs as a possible source of both ECs and HCs cells. As expected, hPSC-EB-CD144+ expressed significantly higher levels of two master transcription factors required for endothelial and hematopoietic specification, ETV2 and RUNX1, compared to the non-sorted 84h-hPSC-EBs and hPSC-EB-CD144− (Fig. 1H-I). No significant differences among the hPSC lines were observed. These results suggest that hPSC-EB-CD144+ may contain not only ECs but also hemogenic ECs.
hPSC-EB-CD144+ generates functional and mature hPSC-ECs
We, then, first sought to determine whether hPSC-EB-CD144+ include proliferative endothelial progenitors. hPSC-EB-CD144+ were seeded at one cell per well in a 96 well plate. Small colonies were observed at day 7 and kept growing until day 15 confirming the existence of proliferative ECs within hPSC-EB-CD144+ (Fig. 2A-B). Similar results were obtained when fresh and frozen hPSC-EB-CD144+ were cultured at a normal cell density. hPSC-EB-CD144+ differentiate into cobblestone-like cells that proliferate until reaching confluence (Fig. 2C). hPSC-ECs were composed of more than 90% cells expressing CD144, CD31, CD34 and vWF (Fig. 2D-E). Most of these cells also expressed CD309 and CD143. We did not observe expression of the hematopoietic markers CD43, CD45 and CD41 (Fig. 2D). hPSC-ECs could be maintained in culture up to 6 passages without modifying their endothelial phenotype (Fig. 2C-D). However, the cumulative population doublings (CPD) decreased during time suggesting hPSC-ECs do not proliferate indefinitely (Fig. 2F). It is worth to note that the CPD was dependent on the hPSC line, being generally smaller for hPSC-ECs coming from A29 compared to hPSC-ECs derived from SA01 and H1.
Next, hPSC-EB-CD144+ were cultured in the presence or absence of an inhibitor of Notch pathway (DAPT). Notch signaling pathway plays a pivotal role in arterial but not in venous differentiation (49). hPSC-EB-CD144+ proliferated and formed a confluent homogeneous layer of hPSC-ECs in both conditions, independently on the hPSC line of origin (supplementary Fig. 2). Their endothelial identity was confirmed by flow cytometry analysis showing that more than 95% of cells expressed CD309, CD144, CD31, CD34 and CD143 independently on the presence or absence of DAPT. No expression of the hematopoietic markers CD43 was detected (supplementary Fig. 2). Similarly, RTqPCR analysis revealed no differences in the expression of the panendothelial marker CDH5, suggesting than inhibition of the Notch pathway does not affect the differentiation of hPSC-EB-CD144+ into hPSC-ECs. No differences in the expression of the arterial marker NRP-1 were observed between DAPT-treated and untreated ECs (Fig. 2G-I). However, DAPT-treated ECs downregulated the arterial marker EFNB2. Some heterogeneity between hPSC lines was observed regarding the expression of venous markers. While hPSC-EC from A29 cell line cultured in the presence of DAPT showed a significantly higher expression of NRP2, EPHB4 and NR2F2 compared to the control, hPSC-EC from SA01 and H1 cell lines only showed upregulation of NRP2 and NR2F2, and NR2F2, respectively (Fig. 2G-I). These results suggest that the hPSC-EB-CD144+ arterial/venous specification may be guided in vitro under the appropriate culture conditions although their response to environmental factors might be influenced by the hPSC line from which they are derived.
We also investigated whether hPSC-EB-CD144+ were able to differentiate into fully mature ECs in vitro. To this, we compared the in vitro functionality of hPSC-ECs to human umbilical cord blood endothelial progenitors, ECFCs. Ac-LDL was observed in the cytoplasm of hPSC-ECs derived from A29, SA01 and H1 hPSCs like ECFCs (Fig. 3A) after 4h of incubation with a fluorochrome-conjugated ac-LDL, suggesting hPSC-ECs are able to endocyte LDL. In addition, more than 90% of hPSC-ECs and ECFC expressed high levels of the adhesion molecule ICAM-1 in response to TNF-α. In contrast, less than 1% of untreated cells expressed similar levels of ICAM-1 suggesting hPSC-ECs may be activated by proinflammatory cytokines (Fig. 3B).
All hPSC-ECs formed a tubular network comparable to the one formed by ECFCs after seeding on Matrigel (Fig. 3C, 3C’ and supplementary Fig. 2). No significant differences in the number of master segments and meshes, neither in the total segment length nor in the branching length were observed among hPSC lines and ECFCs. However, ECFCs formed networks with a slightly reduced number of segments, nodes and branching intervals compared to hPSC-ECs (Fig. 3C’).
Production of nitric oxide (NO) by the endothelial nitric oxide synthase (eNOS) is known to play an important role in endothelial homeostasis, inflammation, and vasodilation. We detected eNOS in 72 to 86% of hPSC-ECs and 84.5% of ECFCs (Fig. 3D). We further analyzed NO production by using the probe DAF-FM diacetate. This permeable molecule is de-acetylated inside the cell and reacts with NO to form fluorescent DAF-FM that can be detectable by flow cytometry. In the absence of the probe, no fluorescence was detected for any of the hPSC-ECs nor ECFCs. However, after incubation with the probe, we detected the fluorescent DAF-FM in 83 to 92% of hPSC-ECs and 88% of ECFCs (Fig. 3D’). Slightly increase of DAF-FM fluorescence was detected when ECs were incubated with LPS (known to induce the NO production) and with a S-nitrosothiol molecule (SNAP) as a positive control.
Finally, we tested the in vivo angiogenic potential of hPSC-ECs in a mouse dorsal chamber model where an ischemic injury has been induced. Since not apparent differences in the functionality of the three hPSC-ECs were observed in vitro, we performed the in vivo experiments with mcherry-ECs derived from SA01 cell line. mcherry-ECs were detected for 4 weeks in the mice. Real time in vivo imaging showed that mcherry-ECs migrated and localized around the proximities of mouse blood vessels at day 5 independently on the induction or not of an ischemia. Intravital microscopy analysis demonstrated hPSC-ECs appear to line up around blood vessels and capillaries (Fig. 3E). Quantification of the blood vessels show no significant differences between the analyzed conditions, although the number of new blood vessels was slightly higher in mice injected with mcherry-ECs (Fig. 3F, supplementary Fig. 3). The immunofluorescent analysis post-mortem confirmed the survival of human CD144+ hPSC-ECs in both ischemic and non-ischemic mice and their localization at the wall of the blood vessels. No teratomas or any other complications were observed at the dorsal chambers. Thus, SA01-ECs survive and participate to the formation of new vessels in vivo (Fig. 3G).
hPSC-EB-CD144+ contain hemogenic hPSC-ECs
To confirm the existence of hemogenic ECs within the hPSC-EB-CD144+ population, fresh and frozen hPSC-EB-CD144+ and hPSC-EB-CD144− cells were seeded on a serum-free methylcellulose enriched with hematopoietic cytokines. After 4 to 6 days, hPSC-BCs were observed for the three-hPSC lines only in the CD144+ fraction (Fig. 4A). Flow cytometry analysis showed that more than 90% of hPSC-BCs were positive for the hematopoietic marker CD43 (Fig. 4B). Analyses of other hematopoietic and endothelial markers revealed heterogeneity within hPSC-BCs. Between 33 to 62% of hPSC-BCs expressed CD41 and 33 to 56% of hPSC-BCs expressed CD45 (Fig. 4B). In addition, a subpopulation of hPSC-BCs coexpressed the hematoendothelial markers CD34 and CD143 (11–46%). Interestingly, only a small fraction of hPSC-BCs retained the expression of CD144 and around 5 to 25% of hPSC-BCs co-expressed CD144 and CD45 markers (Fig. 4C) suggesting that hPSC-EB-CD144+ may undergo endothelial-to-hematopoietic transition (EHT) to give rise to HCs. Even though hPSC-BCs were observed from all the three-hPSC lines, significant differences in the expression of hematopoietic and endothelial markers were observed among the hPSC-BCs. Notably, hPSC-BCs from H1 expressed significantly more CD34, CD143, CD45 and CD41, suggesting a difference in hematopoietic potential among the hPSC cell lines (Fig. 4B-C).
We further analyzed the expression of key hematopoietic transcription factors in hPSC-BCs compared to hPSC-EB-CD144+ and hPSC-ECs RTqPCR. Consistently with their common origin and the role of SCL and GATA2 in both endothelial and hematopoietic differentiation, SCL and GATA2 mRNAs were detected in hPSC-EB-CD144+, hPSC-ECs and hPSC-BCs independently on the hPSC-line (Fig. 4D-E). However, the expression of SCL was significantly higher in hPSC-BCs for A29 and SA01 (Fig. 4D). hPSC-BCs expressed significantly higher levels of the RUNX1 isoform c and GATA1, transcription factors implicated in HSC emergence and erythropoiesis, respectively (Fig. 4F-G). In addition, hPSC-BCs expressed reduced levels of HOXA3 compared to hPSC-EB-CD144+ and hPSC-ECs, a transcription factor key for arterial and hemogenic endothelium determination (Fig. 4H).
Given the hematopoietic profile of hPSC-BCs, we next assessed hPSC-BC hematopoietic potential in vitro. hPSC-BC derived from hiPSCs and hESCs lines gave rise to myeloid and erythroid colonies CFU-GM, CFU-GEMM and BFU-E with an efficiency of 79 to 117 CFU per 105 seeded hPSC-BCs (Fig. 5A-B). From these, most of the colonies accounted for CFU-GM (60% − 87%), around 22–38% corresponded to BFU-E and only 0.1–4% were CFU-GEMM, confirming the heterogeneity previously described by flow cytometry (Fig. 5C). Although similar percentages were found for the three hPSC lines, it is worth to note that hPSC-BCs from SA01 generated slightly more CFU-GEMM (4.3 ± 2.4%) than A29 cell line (0.1 ± 0.4%) and H1 (0.8 ± 0.9%) (Fig. 5C). Since hemoglobin subunits are differentially expressed during the embryonic, fetal, and adult erythropoiesis, we analyzed by RTqPCR the expression of the embryonic (ε), fetal (γ) and adult (β) hemoglobin chains from hPSC-derived CFU-GEMM and BFU-E colonies and we compared them to the embryonic human liver transcripts. CFU-GEMM and BFU-E derived from hPSCs expressed significantly higher levels of ε hemoglobin (HBE1) (10–103 -fold change) compared to ELs (Fig. 5D). In contrast, all the hPSC-derived hematopoietic colonies expressed significantly less γ (HBG1) (0.025-0.12-fold change) and β (HBB) hemoglobin (7x10− 5-1x10− 4-fold change) than human ELs (Fig. 5D). In other words, hPSC-derived CFU-GEMM and BFU-E expressed around 102 times more ε than β hemoglobin (Fig. 5E). These results suggests that our protocol yields principally embryonic erythrocytes and to lesser extent, definitive erythrocytes.
To confirm the existence of bona fide HSCs/HPs in hPSC-EB-CD144+, we analyzed the in vivo hematopoietic reconstitution potential of hPSC-EB-CD144+ in a xenograft model of immunodeficient NSG mice. When hPSC-EB-CD144+ from H1 were intravenously injected, we did not observe any hematopoietic engraftment in the transplanted mice while human cord blood CD34+ cells engrafted efficiently (data not shown). To circumvent any caveat due to homing defects, we injected hPSC-EB-CD144+ cells directly into the right femur of sublethally irradiated NSG mice by intra-BM transplantation (IBMT). Mice were monitored for human reconstitution by peripheral blood analyses (at 6 weeks, see Table 3) and were euthanatized at 12 weeks post-transplant. BM was, then, analyzed for human CD45 expression by FACS analysis and the presence of human DNA in the BM of the right and left femur of transplanted mice was also analyzed by qPCR. Six weeks after transplantation, we observed human CD45+ cells in the peripheral blood of mice injected with CD34+ cells (1% and 9%) and with hPSC-EB-CD144+ from H1 (between 0.5% and 3.6%). We did not detect any human CD45+ in the control mice injected with PBS (Table 3). After 3 months of engraftment, human specific CD45 was detected by flow cytometry only in the right femur in 5 out of 5 mice injected with hPSC-EB-CD144+ from H1 (between 0,1% and 0,6%) (Table 3). In the mice injected with human umbilical cord blood CD34+ cells as control, human specific CD45 was detected in both right (9% and 10%) and left femurs (1% and 5%) (Table 3). These results were further confirmed by qPCR where CD45 human specific DNA was detected in the right femur for hPSC-EB-CD144+ and in the right and left femur for hCD34+ cells from human cord blood. One left femur out of five was also positive for human DNA, however it was negative for CD45 in the flow cytometry analysis. Mice injected with PBS did not express human specific CD45 neither contained human DNA in the femurs (Fig. 5F). No developed lymph nodes or tumors were observed in any of the examined mice. These results suggest hPSC-EB-CD144+ have a hematopoietic reconstitution potential only by IBMT however, this is considerably lower compared to human CD34+ isolated from umbilical cord blood.
Table 3
Hematopoietic reconstitution in NSG mice.
Sample | Human CD34 | hPSC-EB-CD144+ | PBS |
#1 | #2 | #1 | #2 | #3 | #4 | #5 | #1 | #2 |
%hCD45 in the PB | 9% | 1% | 1.5% | 0.6% | 0.5% | 0.9% | 3.6% | 0% | 0% |
%hCD45 in the BM | RF | 10% | 9% | 0.1% | 0.3% | 0.1% | 0.1% | 0.6% | 0% | 0% |
LF | 5% | 1% | 0% | 0% | 0% | 0% | 0% | 0% | 0% |
Human DNA Fold change | RF | 1 | 0.662 | 0.007 | 0.210 | 0.013 | 0.0004 | 0.051 | ND | ND |
LF | 0.175 | 0.008 | ND | ND | 0.0002 | ND | ND | ND | ND |
Percentage of human CD45+ cells detected by flow cytometry in peripheral blood (PB) and in the bone marrow (BM) from right (RF) and left femur (LF) of NSG mice injected with either human CD34+ cord blood cells, hPSC-EB-CD144+ cells or PBS. Relative quantification of human DNA (based on the amplification of human CD45 gene) in BM of mice injected with human cells or PBS. ND: non detected.
Transcriptomic profile of CD144+-EBs, hPSC-ECs and hPSC-BCs
To investigate whether hPSC-EB-CD144+, hPSC-ECs and hPSC-BCs may correspond to their physiological equivalent, we performed a transcriptome analysis to compare hPSC-derived cell populations to undifferentiated hPSC, to differentiated-hPSC themselves and to human EL cell populations. In human, EL is considered as a major organ of hematopoiesis during development. From day 30 to the 30th week of human gestation, mature ECs (CD144+CD34+CD45−) coexist with diverse hematopoietic populations resulted from the migration, proliferation, and maturation of HCs, namely pre-HSCs (CD144+CD34+CD45Low), HSCs and HPs (CD144−CD34+CD45Low) and mature HCs (CD144−CD45High) (14). Thus, EL allow us to compare hPSCs-derived cells to ECs and HCs, including intermediate pre-HSCs no longer found in adults but also HSCs/HPs and mature HCs, all with an identical genetic background.
We first confirmed the identity of the samples by a principal component analysis (PCA). As expected, undifferentiated-hPSCs and hPSC-derived populations were grouped in two independent groups (Fig. 6A). In addition, hPSC-derived cells did not cluster with EL samples, yet they were spatially closer to EL than to hPSCs (Fig. 6A). Independent PCA within each group (EL and hPSC-derived populations) clustered samples according to their endothelial and hematopoietic identity. As shown in Fig. 6B, EL sample distribution along the first component of EL-PCA showed a clear difference between EL-ECs and EL-mature HCs, located at the opposite extremes of the axis, while EL-pre-HSCs and EL-HSCs/HPs are in between these two populations, as result of the EHT giving place to HCs. A similar analysis within hPSC-derived cells clustered hPSC-EB-CD144+ and hPSC-ECs in a separate group from hPSC-BC-CD144+CD45+ sorted and unsorted hPSC-BCs (Fig. 6C). Supervised analysis performed to identify transcripts specifically up- or down-regulated in hPSC-derived cell populations confirmed these results. Indeed, when hPSC-ECs and hPSC-BC-CD144+CD45+ were compared to hPSC-EBs, EC gene signatures were identified specifically in hPSC-ECs, whereas hemogenic and hematopoietic cell gene signatures were found in hPSC-BC-CD144+CD45+, confirming distinct endothelial, hemogenic endothelial and hematopoietic transcriptional profile in hPSC-derived populations (Fig. 6D).
To confirm these results, we performed a Gene Set Enrichment Analysis (GSEA) between each hPSC-derived cell population and undifferentiated hPSCs. hPSC-EB-CD144+ expressed a discrete list of endothelial and hematopoietic-related terms (see Fig. 6E and supplementary table 5), including Embryonic Hemopoiesis and Endothelium Development ones. Endothelial terms were mainly enriched in hPSC-ECs and hPSC-EB-CD144+ compared to hPSC-BC-CD144+45+ and hPSC-BCs (Fig. 6E-G and supplementary table 5). In addition, hematopoietic terms including Embryonic_Hemopoiesis were enriched in hPSC-EB-CD144+, hPSC-BCs and hPSC-BC-CD144+45+ but not in hPSC-ECs (Fig. 6E-G and supplementary table 5). From all the hPSCs-derived populations, hPSC-ECs contained the larger group of endothelial-related terms, including the following Blood vessel morphogenesis, Endothelial Proliferation, Lymphangiogenesis and Endothelial Barrier (see Fig. 6F and supplementary table 5). hPSC-ECs display also few terms involved in HSCs/HPs differentiation and term related to erythroid, myeloid, and megakaryocytic hematopoietic lineages. In contrast to hPSC-EC, hPSC-BCs and hPSC-BC-CD144+CD45+ contained a larger list of hematopoietic-related terms involved principally in the differentiation of HSCs/HPs and the maturation of HCS from erythroid, myeloid, megakaryocytes but also terms related to lymphoid lineages (see Fig. 6G and supplementary table 5). It is worth to mention that the terms related to leukocytes, lymphocytes and immune system were largely more enriched in hPSC-BC-CD144+CD45+ than in unsorted hPSC-BC (Fig. 6G). These results suggest that hPSC-BC and more specifically hPSC-BC-CD144+CD45+ have the potential to generate multilineage hematopoiesis contrary to hPSC-EB-CD144+ that were only enriched in terms involved in erythroid, myeloid, and megakaryocytic lineages (see supplementary table 5).
Once confirmed the identity of hPSC-derived populations, we performed an unsupervised hierarchical clustering between the hPSC-derived populations and the EL populations. As illustrated in Fig. 6H, two main clusters were obtained coinciding with the origin of the samples and not with their identity (e.g. endothelial or hematopoietic). An exception was observed for EL-HSCs/HPs, which were clustered together with non-sorted hPSC-BCs and hPSC-BC-CD144+CD45+, suggesting these two hPSC-derived populations are closed to physiological EL-HSCs/HPs (Fig. 6H).
We then performed a GSEA between hPSC-derived populations and their closest EL-equivalent population. Compared to EL-pre-HSCs, hPSC-EB-CD144+ displayed a stronger endothelial signature. However, they were also enriched in a discrete list of terms related to EL-HSCs/HPs differentiation. This hematopoietic potential was further confirmed by GSEA between hPSC-EB-CD144+ and EL-ECs. Interestingly, we did not observe enriched terms related to mature hematopoietic lineages (supplementary table 6).
When hPSC-ECs were compared to EL-ECs, hPSC-ECs exhibited a more mature endothelial phenotype compared to EL-ECs since they were enriched on terms associated to endothelial function so as platelet aggregation or myeloid-leukocyte migration but also on terms related to HSCs/HPs differentiation (supplementary table 6).
Despite comparison between hPSC-derived populations and undifferentiated hPSC showed that hPSC-BCs have a hematopoietic signature, GSEA between hPSC-BCs and EL-HSCs/HPs revealed that hPSC-BCs were enriched in endothelial terms when compared to human EL-HSCs/HPs hematopoietic populations. Furthermore, hPSC-BCs were not enriched in hematopoietic multilineage-related terms when compared to EL-HSCs/HPs. In contrast, hPSC-BC-CD144+-CD45+ exhibit a long list of terms related to multilineage hematopoietic progenitors when compared to EL-HSCs/HPs. Regarding their endothelial phenotype, hPSC-BC-CD144+-CD45+ expressed more endothelial-related terms than EL-HSCs/HPs but surprisingly less than EL-pre-HSC (supplementary table 6). These results suggest hPSC-BC-CD144+-CD45+ are more engaged towards hematopoietic lineages and have a largest multilineage potential compared to hPSC-BC and may correspond to a population in between EL-pre-HSCs and EL-HSCs/HSCs.
Thus, these transcriptomic analyses suggest that hPSC differentiation into endothelial and hematopoietic cells recapitulates the human hematopoietic developmental processes through an EHT leading first to primitive hematopoiesis progenitors and thereafter to definitive hematopoietic progenitors (Fig. 6I).