Example 1

Temporal modulation of canonical Wnt signaling has been shown to be sufficient to generate functional cardiomyocytes at high yield and purity from numerous hPSC lines (Lian, X. et al. (2012) Proc. Natl. Acad. Sci. USA 109:E1848-1857; Lian, X. et al. (2013) Nat. Protoc. 8:162-175). In this approach, Wnt/β-catenin signaling first is activated in the hPSCs, followed by an incubation period, followed by inhibition of Wnt/β-catenin signaling. In the originally published protocol, Wnt/β-catenin signaling activation was achieved by incubation with the Gsk3 inhibitor CHIR99021 (GSK-3 α, IC50=10 nM; GSK-3, β IC50=6.7 nM) and Wnt/β-catenin signaling inhibition was achieved by incubation with the Porcn inhibitor IWP2 (IC50=27 nM). Because we used Gsk3 inhibitor and Wnt production inhibitor for cardiac differentiation, this protocol was termed GiWi protocol. To improve the efficiency of the original protocol and reduce the potential side effects of the small molecules used in the original protocol, a second generation protocol was developed that uses another set of small molecules with higher inhibition potency. In this second generation GiWi protocol, Wnt/β-catenin signaling activation was achieved by incubation with the Gsk3 inhibitor CHIR98014 (CAS 556813-39-9; commercially available from, e.g., Selleckchem) (GSK-3 α, IC50=0.65 nM; GSK-3, β IC50=0.58 nM) and Wnt/β-catenin signaling inhibition was achieved by incubation with the Porcn inhibitor Wnt-059 (CAS 1243243-89-1; commercially available from, e.g., Selleckchem or Tocris) (IC50=74 pM). The Gsk3 inhibitor CHIR98014 was used to promote cardiac mesodermal differentiation, whereas the Porcn inhibitor Wnt-059 was used to enhance ventricular progenitor differentiation from mesoderm cells.

For cardiomyocyte differentiation via the use of these small molecules, hPSCs were maintained on Matrigel (BD Biosciences) coated plates (Corning) in E8 medium (described in Chen, G. et al. (2011) Nature Methods, 8:424-429; commercially available; STEMCELL Technologies) or mTeSR1 medium (commercially available; STEMCELL Technologies). Suitable hPSCs include induced pluripotent stem cells (iPSCs) such as 19-11-1, 19-9-7 or 6-9-9 cells (Yu, J. et al. (2009) Science, 324:797-801) and human embryonic stem cells (hESCs), such as ES03 (WiCell Research Institute) and H9 cells (Thomson, J. A. et al. (1998) Science, 282:1145-1147).

hPSCs maintained on a Matrigel-coated surface in mTeSR1 medium were dissociated into single cells with Accutase (Life Technologies) at 37° C. for 5 minutes and then seeded onto a Matrigel-coated cell culture dish at 100,000-200,000 cells/cm2 in mTeSR1 medium supplemented with 5 μM ROCK inhibitor Y-27632 (Selleckchem)(day −2) for 24 hours. Cells were then cultured in mTeSR1, changed daily. At day 0, cells were then treated with 1 μM Gsk3 inhibitor CHIR98014 (Selleckchem) for 24 hours (day 0 to day 1) in RPMI/B27-ins (500 ml RPMI with 10 ml B27 supplement without insulin). The medium was then changed to the corresponding medium containing 2 μM the Porcn inhibitor Wnt-059 (Selleckchem) at day 3, which was then removed during the medium change on day 5. Cells were maintained in RPMI/B27 (stock solution: 500 ml RMPI medium+10 ml B27 supplement) starting from day 7, with the medium changed every three days. This exemplary culturing protocol for generating cardiomyogenic progenitor cells is illustrated schematically in the drawing.

Flow cytometry and immunostaining were preformed to examine the expression of particular lineage markers. After 24 hour treatment with CHIR-98014, more than 99% of the hPSCs expressed the mesoderm marker Brachyury. Three days after treatment with CHIR-98014, more than 95% of differentiated cells expressed Mesp1, which marks the cardiac mesoderm. The culture protocol not only allowed the cells to synchronously differentiate into the cardiac mesodermal lineage, but also reproducibly generated more than 90% of ventricular myocytes after 14 days of differentiation, as determined by cTnT flow cytometry and electrophysiology analysis.

To further assess cardiac differentiation of the hPSCs over time, Western blot analysis was performed on days 0-7 and d11 to examine the expression of Isl1 and Nkx2.5 (cardiomyogenic progenitor markers) and cTnI (a cardiac myocyte marker). Cells were lysed in M-PER Mammalian Protein Extraction Reagent (Pierce) in the presence of Halt Protease and Phosphatase Inhibitor Cocktail (Pierce). Proteins were separated by 10% Tris-Glycine SDS/PAGE (Invitrogen) under denaturing conditions and transferred to a nitrocellulose membrane. After blocking with 5% dried milk in TBST, the membrane was incubated with primary antibody overnight at 4° C. The membrane was then washed, incubated with an anti-mouse/rabbit peroxidase-conjugated secondary antibody at room temperature for 1 hour, and developed by SuperSignal chemiluminescence (Pierce). During cardiac differentiation of hPSCs, Isl1 expression started on day 4 and increased to its maximum expression on day 6, whereas NK×2.5 only started to express on day 6 and reached its maximum expression after day 10. Cardiomyoctes (cTnI+ cells) were not induced until day 11 of differentiation.

In addition, immunostaining of the day 6 cells was performed for Isl1 expression. Cells were fixed with 4% formaldehyde for 15 minutes at room temperature and then stained with primary (anti-Isl1) and secondary antibodies in PBS plus 0.4% Triton X-100 and 5% non-fat dry milk (Bio-Rad). Nuclei were stained with Gold Anti-fade Reagent with DAPI (Invitrogen). An epifluorescence microscope (Leica DM IRB) with a QImaging® Retiga 4000R camera was used for imaging analysis. The results showed substantial numbers of Isl1+ cells.

Flow cytometry analysis of day 6 cells for Isl1 expression also was performed. Cells were dissociated into single cells with Accutase for 10 minutes and then fixed with 1% paraformaldehyde for 20 minutes at room temperature and stained with primary and secondary antibodies in PBS 0.1% Triton X-100 and 0.5% BSA. Data were collected on a FACSCaliber flow cytometer (Beckton Dickinson) and analyzed using FloJo. The results showed that more than 95% of cells expressed Isl1 at this stage.

In summary, this example provides a protocol for human ventricular progenitor generation (HVPG protocol) that allows for the large-scale production of billions of Isl1+ human HPVs efficiently within 6 days.

Example 2

To profile the transcriptional changes that occur during the cardiac differentiation process at a genome-scale level, RNA sequencing (RNA-seq) was performed at different time points following differentiation to build cardiac development transcriptional landscapes. We performed RNA-seq experiments on day 0 to day 7 samples, as well as day 19 and day 35 samples (two independent biological replicates per time point). Two batches of RNA-seq (100 bp and 50 bp read length) were performed using the illumine Hiseq 2000 platform. In total, 20 samples were examined. Bowtie and Tophat were used to map our reads into a reference human genome (hg19) and we calculate each gene expression (annotation of the genes according to Refseq) using RPKM method (Reads per kilobase transcript per million reads). Differentiation of hPSCs to cardiomyocytes involves five major cell types: pluripotent stem cells (day 0), mesoderm progenitors (day 1 to day 2), cardiac mesoderm cells (day 3 to day 4), heart field progenitors (day 5, day 6 and day 7), and cardiomyocytes (day 10 after).

Molecular mRNA analysis of cardiac differentiation from hPSCs using the HVPG protocol revealed dynamic changes in gene expression, with down-regulation of the pluripotency markers OCT4, NANOG and SOX2 during differentiation. Induction of the primitive streak-like genes T and MIXL1 occurred within the first 24 hours following CHIR-98014 addition, and was followed by upregulation of the cardiac mesodermal marker MESP1 on day 2 and day 3. Expression of the cardiac muscle markers TNNT2, TNNC1, MYL2, MYL7, MYH6, MYH7 and IRX4 was detected at later stage of differentiation (after day 10).

By this analysis, genes enriched at each differentiation stage, including mesoderm cells, cardiac progenitors and cardiomyocytes, were identified. Mesoderm cells, which are related to day 1 differentiated cells, express brachyury. We identified potential surface markers for mesoderm cells, including: FZD10, CD48, CD1D, CD8B, IL15RA, TNFRSF1B, TNFSF13, ICOSLG, SEMA7A, SLC3A2, SDC1, HLA-A. Through similar analysis, we also identified surface markers for cardiac mesoderm mesp1 positive cells, including: CXCR4, ANPEP, ITGA5, TNFRSF9, FZD2, CD1D, CD177, ACVRL1, ICAM1, L1CAM, NGFR, ABCG2, FZD7, TNFRSF13C, TNFRSF1B.

Consistent with western blot analysis, ISL1 mRNA was expressed as early as day 4 and peaked on day 5, one day before its protein expression reached its peak. On day 5 of differentiation (the cardiac progenitor stage, isl1 mRNA expression maximum on day 5, isl1 protein expression maximum on day 6), the day 5 enriched genes were compared with an anti-CD antibody array (a panel of 350 known CD antibodies) and a number of potential cell-surface protein markers were identified. We identified many cell-surface proteins expressed at this stage, including: FZD4, JAG1, PDGFRA, LIFR (CD118), TNFSF9, FGFR3.

The cell surface protein Jagged 1 (JAG1), Frizzled 4 (FZD4), LIFR (CD118) and FGFR3 were selected for further analysis. Jagged 1 expression was further studied as described below and in Example 3. Frizzled 4 expression was further studied as described in Example 4. LIFR (CD118) and FGFR3 expression was further studied as described in Example 5.

Firstly, the expression of Isl1 and Jag1 was profiled using the double staining flow cytometry technique. Flow cytometric analysis was carried out essentially as described in Example 1, using anti-Isl1 and anti-Jag1 antibodies for double staining. Jagged 1 expression was found to trace the expression of Islet 1 and on day 6 of differentiation, all of the Islet 1 positive cells also expressed Jagged 1, and vice versa. Because of the co-expression pattern of these two markers, a Jagged 1 antibody was used to enrich the 94.1% Islet 1+ cells differentiated population to 99.8% purity of Islet1+Jagged1+ cells.

It also was confirmed that Islet 1 is an earlier developmental gene than the Nkx2.5 gene using double immunostaining of ISL1 and NKX2.5 expression in HVPs. The purified HVPs uniformly express the ISL1 gene, but at this stage, only a few of the cells started to express Nkx2.5.

Furthermore, immunostaining with both anti-Isl1 and anti-Jag 1 was performed, essentially as described in Example 1, on week 4 human fetal heart tissue, neonatal heart tissue and 8-year old heart tissue. The results revealed that in the in vivo fetal heart, all of the Islet 1 positive cells also expressed Jagged 1. However, the neonatal heart and 8-year old heart did not express Islet 1 or Jagged 1. In the ventricle of week 4 human fetal heart, cardiac Troponin T (cTnT) staining revealed visible sarcomere structures. In addition, over 50% of ventricular cells in the week 4 fetal heart expressed both Islet1 and Jagged1, which was markedly decreased during subsequent maturation, with the loss of expression of both Islet1 and Jagged1 in the ventricular muscle cells of the human neonatal hearts.

The above-described experiments demonstrate that Jagged 1 is a cell surface marker for Islet 1 positive cardiomyogenic progenitor cells.

Example 3

To characterize the clonal differentiation potential of Isl1+Jag1+ cells, cardiomyogenic progenitor cells were generated by the culturing protocol described in Example 1, and one single Isl1+Jag1+ cell was seeded into one well of a Matrigel-coated 48-well plate. Cells were purified with antibody of Jag1 and then one single cell was seeded into one well. The single cells were then cultured for 3 weeks in Cardiac Progenitor Culture (CPC) medium (advanced DMEM/F12 supplemented with 2.5 mM GlutaMAX, 100 μg/ml Vitamin C, 20% Knockout Serum Replacement).

Immunostaining of the 3-week differentiation cell population was then performed with three antibodies: cardiac troponin I (cTn1) for cardiomyocytes, CD144 (VE-cadherin) for endothelial cells and smooth muscle actin (SMA) for smooth muscle cells. The results showed that the single cell-cultured, Isl1+Jag1+ cells gave rise to cTnI positive and SMA positive cells, but not VE-cadherin positive endothelial cells, indicating these generated Islet1+ cells are heart muscle progenitors that have limited differentiation potential to endothelial lineages. Purified Islet1+Jagged1+ cells differentiated with the HVPG protocol from human induced pluripotent stem cells (iPSC 19-9-11 line) also showed similar in vitro differentiation potential and predominantly differentiate to cTnI+SMA+ cells, but not VE-cadherin+ cells. Over the course of several weeks, the cells expressed the ventricular specific marker MLC2v, indicating that the initial ISL1+ subset was already committed to the ventricular cell fate. Because of the limited vascular differentiation potential of Islet1+ cells generated using the HVPG protocol, these generated Islet1+ cells might represent a distinct progenitor population from the previously reported KDR+ population (Yang, L. et al. (2008) Nature 453:524-528) or multipotent ISL1+ cells (Bu, L. et al. (2009) Nature 460:113-117; Moretti, A. et al. (2006) Cell 127:1151-1165), which can give rise to all three lineages of cardiovascular cells.

These results demonstrated that the Isl1+Jag1+ cardiomyogenic progenitor cells can be successfully cultured in vitro from a single cell to a significantly expanded cell population (1×109 cells or greater) that contains all three types of cardiac lineage cells, with a predominance of cardiomyocytes. Furthermore, these cells can be cultured in vitro for extended periods of time, for at least 2-3 weeks, and even for months (e.g., six months or more). Since the cardiomyogenic progenitor cells gradually differentiate into cardiomyocytes, which do not proliferate, a culture period of approximately 2-3 weeks is preferred.

Example 4

As described in Example 2, Frizzled 4 (FZD4) was identified by RNA-seq analysis as being expressed in cardiac progenitor cells. Thus, to confirm FZD4 as a cell surface marker of cardiac progenitor cells, FZD4 expression was assessed during cardiac differentiation via Western blot analysis. The results demonstrated that FZD4 was not expressed in pluripotent stem cells and the first 3 days differentiated cells. However, FZD4 started to express on day 4 and maximize its expression on day 5 of expression.

In order to quantify the co-expression pattern of FZD4 and Isl1 at the single cell level, FACS analysis was performed. On day 5 of differentiation, more than 83% of cells express both isl1 and FZD4, demonstrating that FZD4 is a cell surface marker for Isl1 positive cells during cardiac progenitor differentiation using the GiWi protocol.

In order to confirm that both JAG1 and FZD4 were indeed co-expressed with ISL1 on the human ventricular progenitor cells, triple immunofluorescence analysis of day 6 differentiated cells from hPSCs was performed with antibodies to Islet 1, Jagged 1 and Frizzled 4. The triple staining experiment demonstrated that Isl1+ cells expressed both Jagged 1 and Frizzled 4.

Example 5

As described in Example 2, LIFR(CD118) and FGFR3 were identified by RNA-seq analysis as being expressed in cardiac progenitor cells. In this example, expression of these additional cell surface markers for the human ventricular progenitor cells was confirmed by flow cytometry analysis. Human ventricular progenitor (HVP) cells were generated as described in Example 1 or 11 and day 6 cells were analyzed by standard flow cytometry. A double staining flow cytometry experiment using anti-Islet 1 and anti-Leukemia Inhibitory Factor Receptor (LIFR) antibodies was performed. The results demonstrate that the HVP cells co-express Islet 1 and LIFR, thereby confirming that LIFR is a cell surface marker for the HVP cells. Furthermore, flow cytometry experiments were performed comparing the expression of LIFR and Fibroblast Growth Factor Receptor 3 (FGFR3) on day 6 HVP cells to undifferentiated embryonic stem (ES) cells. The results demonstrate that LIFR and FGFR3 are both highly enriched for expression on the HVP cells, thereby confirming that LIFR and FGFR3 are both cell surface markers for the HVP cells.

Example 6

The ES03 human embryonic stem cell (hESC) line (obtained from WiCell Research Institute) expresses green fluorescent protein (GFP) driven by the cardiac-specific cTnT promoter. ES03 cells were used to generate Isl1+Jag1+ cardiomyogenic progenitor cells using the culturing protocol described in Example 1. The Isl1+Jag1+ cardiomyogenic progenitor cells were transplanted into the hearts of severe combined immunodeficient (SCID) beige mice to document their developmental potential in vivo.

Briefly, Isl1+Jag1+ cells were injected (1,000,000 cells per recipient) directly into the left ventricular wall of NOD/SCID-gamma mice in an open-chest procedure. Hearts were harvested 2-3 weeks post surgery, fixed in 1% PFA and sectioned at 10 μm (n=12). Histological analyses of the hearts of the transplanted mice revealed the presence of GFP+ donor cells, detected by epifluorescence and by staining with an anti-GFP antibody, demonstrating that the Isl1+Jag1+ cardiomyogenic progenitor cells were capable of differentiating into cardiomyocytes when transplanted in vivo.

The Isl1+Jag1+ cardiomyogenic progenitor cells were also transplanted directly into infarcted hearts of SCID beige mice (“injured mice”), as compared to similarly transplanted normal mice. When analyzed two weeks later, injured mice transplanted with the Isl1+Jag1+ cardiomyogenic progenitor cells had a larger graft size than the normal mice similarly transplanted, demonstrating the cardiomyocyte regeneration capacity of the Isl1+Jag1+ cardiomyogenic progenitor cells in vivo.

Example 12

In this example, genes in the angiogenic family that are expressed in human ventricular progenitor cells (HVPs) were identified. HVPs were generated as described in Examples 1 or 11 and RNA sequencing (RNA-seq) was performed at different time points following differentiation as described in Example 2. Cluster analysis of gene expression profiles at different time points during HVP differentiation identified stage-specific signature genes. These genes were clustered hierarchically on the basis of the similarity of their expression profiles. First, genes showing expression in four different categories were identified: (i) cell surface expression; (ii) co-expression with Islet 1; (iii) high expression on day 5 of differentiation; and (iv) high d5/d0 ratio. This analysis confirmed the cell surface markers for HVPs of: JAG1, FZD4, FGFR3, LIFR (CD118) and TNFSF9. Next, from this same population of HVPs that identified the cell surface markers, gene ontogeny searches were performed to identify angiogenic family genes that were expressed in this population of HVPs, to thereby identify a gene fingerprint profile that identifies genes critical for cell engraftment.

Statistically, Pearson's correlation with Isl1 expression was used to identify those angiogenic genes whose expression in the HVPs best correlated with Isl1 expression. Table 1 below lists the angiogenic genes that correlate with Isl1 expression with a Pearson's correlation of 0.50 or higher.

TABLE 1
Angiogenic genes expressed in HVPs with a Pearson
Correlation with Isl1 Expression of 0.50 or greater
Pearson's
Correlation with
GeneAngiogenic genes (GO:0001525)Isl1 Expression
FGF10fibroblast growth factor 100.98
PRKD1protein kinase D10.95
CCBE1collagen and calcium binding EGF domains 10.94
PDGFRAplatelet-derived growth factor receptor, alpha polypeptide0.94
EPHB2EPH receptor B20.92
GATA2GATA binding protein 20.92
NTRK1neurotrophic tyrosine kinase, receptor, type 10.92
PTGISprostaglandin I2 (prostacyclin) synthase0.87
BMPERBMP binding endothelial regulator0.85
BMP4bone morphogenetic protein 40.84
C1GALT1core 1 synthase, glycoprotein-N-acetylgalactosamine 3-0.84
beta-galactosyltransferase 1
MEIS1Meis homeobox 10.83
TBX1T-box 10.83
PKNOX1PBX/knotted 1 homeobox 10.83
ID1inhibitor of DNA binding 1, dominant negative helix-loop-0.82
helix protein
TCF21transcription factor 210.82
HEY1hes-related family bHLH transcription factor with YRPW0.80
motif 1
HOXB3homeobox B30.78
JAG1jagged 10.75
HGFhepatocyte growth factor (hepapoietin A; scatter factor)0.74
IL6interleukin 60.74
GHRLghrelin/obestatin prepropeptide0.73
IHHindian hedgehog0.70
SRPK2SRSF protein kinase 20.70
GATA6GATA binding protein 60.69
HAND1heart and neural crest derivatives expressed 10.69
AMOTangiomotin0.69
NRP2neuropilin 20.65
PTENphosphatase and tensin homolog0.65
SEMA3Esema domain, immunoglobulin domain (Ig), short basic0.64
domain, secreted, (semaphorin) 3E
APOLD1apolipoprotein L domain containing 10.62
SETD2SET domain containing 20.62
DAB2IPDAB2 interacting protein0.61
KDRkinase insert domain receptor0.60
PGFplacental growth factor0.60
EMP2epithelial membrane protein 20.59
TAL1T-cell acute lymphocytic leukemia 10.58
ACVR1activin A receptor, type I0.58
HIPK2homeodomain interacting protein kinase 20.56
CSPG4chondroitin sulfate proteoglycan 40.55
TNFAIP3tumor necrosis factor, alpha-induced protein 30.55
NRP1neuropilin 10.55
NFATC4nuclear factor of activated T-cells, cytoplasmic,0.54
calcineurin-dependent 4
CDC42cell division cycle 420.54
ANGPTL4angiopoietin-like 40.53
BCAS3breast carcinoma amplified sequence 30.53
HIPK1homeodomain interacting protein kinase 10.53
NRXN3neurexin 30.52
FZD5frizzled class receptor 50.52
HHEXhematopoietically expressed homeobox0.50
Table 2 below lists the angiogenic genes that correlate with Isl1 expression with a Pearson's correlation of 0.49-0.00.

TABLE 2
Angiogenic genes expressed in HVPs with a Pearson
Correlation with Isl1 Expression of 0.49 to 0.00
Pearson's
Correlation with
GeneAngiogenic genes (GO:0001525)Isl1 Expression
ACVRL1activin A receptor type II-like 10.49
ENPEPglutamyl aminopeptidase (aminopeptidase A)0.49
EFNA1ephrin-A10.49
CHRNA7cholinergic receptor, nicotinic, alpha 7 (neuronal)0.49
TMEM100transmembrane protein 1000.48
NOS3nitric oxide synthase 3 (endothelial cell)0.47
LEF1lymphoid enhancer-binding factor 10.47
NRXN1neurexin 10.46
EPHB3EPH receptor B30.44
ROCK1Rho-associated, coiled-coil containing protein kinase 10.42
NF1neurofibromin 10.42
CYSLTR2cysteinyl leukotriene receptor 20.42
FGFR2fibroblast growth factor receptor 20.41
GATA4GATA binding protein 40.40
FMNL3formin-like 30.40
C3complement component 30.40
WASF2WAS protein family, member 20.40
CALCRLcalcitonin receptor-like0.39
HIF1Ahypoxia inducible factor 1, alpha subunit (basic helix-loop-0.39
helix transcription factor)
VEGFAvascular endothelial growth factor A0.39
KRIT1KRIT1, ankyrin repeat containing0.39
CDH13cadherin 130.39
COL18A1collagen, type XVIII, alpha 10.39
STK4serine/threonine kinase 40.38
C5complement component 50.38
HDAC7histone deacetylase 70.38
ANGPT2angiopoietin 20.38
PLCG1phospholipase C, gamma 10.37
EDNRAendothelin receptor type A0.35
TGFB2transforming growth factor, beta 20.35
HAND2heart and neural crest derivatives expressed 20.35
CD34CD34 molecule0.35
BTG1B-cell translocation gene 1, anti-proliferative0.34
TGFBR1transforming growth factor, beta receptor 10.33
FGFR1fibroblast growth factor receptor 10.33
FN1fibronectin 10.31
TWIST1twist family bHLH transcription factor 10.31
ELK3ELK3, ETS-domain protein (SRF accessory protein 2)0.30
THSD7Athrombospondin, type I, domain containing 7A0.30
RGCCregulator of cell cycle0.30
PLCD1phospholipase C, delta 10.29
SPARCsecreted protein, acidic, cysteine-rich (osteonectin)0.29
TBX20T-box 200.28
PIK3CAphosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic0.27
subunit alpha
MMRN2multimerin 20.27
FOXO4forkhead box O40.26
RAMP2receptor (G protein-coupled) activity modifying protein 20.25
FLT1fms-related tyrosine kinase 10.25
ADRB2adrenoceptor beta 2, surface0.25
SLC12A6solute carrier family 12 (potassium/chloride transporter),0.25
member 6
ADMadrenomedullin0.25
NPPBnatriuretic peptide B0.24
SPINK5serine peptidase inhibitor, Kazal type 50.24
MAPK14mitogen-activated protein kinase 140.24
MMP2matrix metallopeptidase 20.24
PTPRMprotein tyrosine phosphatase, receptor type, M0.23
OVOL2ovo-like zinc finger 20.23
CTNNB1catenin (cadherin-associated protein), beta 1, 88 kDa0.22
OTULINOTU deubiquitinase with linear linkage specificity0.21
B4GALT1UDP-Gal:betaGlcNAc beta 1,4-galactosyltransferase,0.21
polypeptide 1
PDGFRBplatelet-derived growth factor receptor, beta polypeptide0.20
F3coagulation factor III (thromboplastin, tissue factor)0.20
PRKCAprotein kinase C, alpha0.20
LRP5low density lipoprotein receptor-related protein 50.20
MAP3K7mitogen-activated protein kinase kinase kinase 70.20
NRCAMneuronal cell adhesion molecule0.19
MAP2K5mitogen-activated protein kinase kinase 50.18
S1PR1sphingosine-1-phosphate receptor 10.18
NFATC3nuclear factor of activated T-cells, cytoplasmic,0.18
calcineurin-dependent 3
TSPAN12tetraspanin 120.18
LAMA5laminin, alpha 50.17
LOXL2lysyl oxidase-like 20.17
ANGPT1angiopoietin 10.17
GTF2Igeneral transcription factor IIi0.16
E2F8E2F transcription factor 80.16
PDE3Bphosphodiesterase 3B, cGMP-inhibited0.15
SHBSrc homology 2 domain containing adaptor protein B0.14
MYH9myosin, heavy chain 9, non-muscle0.14
FZD8frizzled class receptor 80.14
NOVnephroblastoma overexpressed0.14
SH2D2ASH2 domain containing 2A0.14
FGF8fibroblast growth factor 8 (androgen-induced)0.13
TIE1tyrosine kinase with immunoglobulin-like and EGF-like0.13
domains 1
EGLN1egl-9 family hypoxia-inducible factor 10.12
RORARAR-related orphan receptor A0.11
MFGE8milk fat globule-EGF factor 8 protein0.11
ARHGAP24Rho GTPase activating protein 240.10
ITGA5integrin, alpha 5 (fibronectin receptor, alpha polypeptide)0.10
PARVAparvin, alpha0.10
ADIPOR2adiponectin receptor 20.09
NPR1natriuretic peptide receptor 10.09
ITGB1integrin, beta 1 (fibronectin receptor, beta polypeptide,0.09
antigen CD29 includes MDF2, MSK12)
HIF3Ahypoxia inducible factor 3, alpha subunit0.08
EPAS1endothelial PAS domain protein 10.08
FOXC2forkhead box C20.07
ANXA2annexin A20.06
RBM15RNA binding motif protein 150.06
PITX2paired-like homeodomain 20.06
FOXC1forkhead box C10.06
SRFserum response factor0.06
ECSCRendothelial cell surface expressed chemotaxis and0.05
apoptosis regulator
SOX17SRY (sex determining region Y)-box 170.04
HDAC5histone deacetylase 50.04
LRG1leucine-rich alpha-2-glycoprotein 10.04
ADAM8ADAM metallopeptidase domain 80.03
UBP1upstream binding protein 1 (LBP-1a)0.02
VASH1vasohibin 10.02
ANXA3annexin A30.01
RRASrelated RAS viral (r-ras) oncogene homolog0.01
TYMPthymidine phosphorylase0.01
PRCPprolylcarboxypeptidase (angiotensinase C)0.01
SEMA5Asema domain, seven thrombospondin repeats (type 1 and0.00
type 1-like), transmembrane domain (TM) and short
cytoplasmic domain, (semaphorin) 5A
GREM1gremlin 1, DAN family BMP antagonist0.00

Angiogenic genes whose expression negatively correlated with Isl1 expression in the HVPs were also identified. Table 3 below lists the angiogenic genes that negatively correlate with Isl1 expression with a Pearson's correlation of −0.50 or less.

TABLE 3
Angiogenic genes expressed in HVPs with a Pearson
Correlation with Isl1 Expression of −0.50 or less
Pearson's
Correlation with
GeneAngiogenic genes (GO:0001525)Isl1 Expression
ETS1v-ets avian erythroblastosis virus E26 oncogene homolog 1−0.50
BAXBCL2-associated X protein−0.50
XBP1X-box binding protein 1−0.52
TDGF1teratocarcinoma-derived growth factor 1−0.53
C5AR1complement component 5a receptor 1−0.53
EPHA1EPH receptor A1−0.53
HS6ST1heparan sulfate 6-O-sulfotransferase 1−0.56
SHC1SHC (Src homology 2 domain containing) transforming−0.56
protein 1
SP100SP100 nuclear antigen−0.58
JAM3junctional adhesion molecule 3−0.58
CASP8caspase 8, apoptosis-related cysteine peptidase−0.60
FLT4fms-related tyrosine kinase 4−0.60
SFRP2secreted frizzled-related protein 2−0.61
HPSEheparanase−0.61
BAK1BCL2-antagonist/killer 1−0.65
GPX1glutathione peroxidase 1−0.65
VAV3vav 3 guanine nucleotide exchange factor−0.70
VAV2vav 2 guanine nucleotide exchange factor−0.72
EGFepidermal growth factor−0.72
ADAM15ADAM metallopeptidase domain 15−0.73
AGGF1angiogenic factor with G patch and FHA domains 1−0.76
Table 4 below lists the angiogenic genes that negatively correlate with Isl1 expression with a Pearson's correlation of −0.01 to −0.49

TABLE 4
Angiogenic genes expressed in HVPs with a Pearson Correlation
with Isl1 Expression of −0.01 to −0.49
Pearson's
Correlation with
GeneAngiogenic genes (GO:0001525)Isl1 Expression
EIF2AK3eukaryotic translation initiation factor 2-alpha kinase 3−0.01
ROCK2Rho-associated, coiled-coil containing protein kinase 2−0.01
WNT5Awingless-type MMTV integration site family, member 5A−0.02
NR4A1nuclear receptor subfamily 4, group A, member 1−0.02
CYP1B1cytochrome P450, family 1, subfamily B, polypeptide 1−0.02
PTK2protein tyrosine kinase 2−0.03
SFRP1secreted frizzled-related protein 1−0.04
STAT1signal transducer and activator of transcription 1, 91 kDa−0.04
ITGAVintegrin, alpha V−0.04
EPHB4EPH receptor B4−0.05
CYR61cysteine-rich, angiogenic inducer, 61−0.05
TEKTEK tyrosine kinase, endothelial−0.06
COL15A1collagen, type XV, alpha 1−0.06
COL4A1collagen, type IV, alpha 1−0.07
ANGangiogenin, ribonuclease, RNase A family, 5−0.07
HSPB1heat shock 27 kDa protein 1−0.07
PLXND1plexin D1−0.08
HSPG2heparan sulfate proteoglycan 2−0.09
VEGFCvascular endothelial growth factor C−0.09
SYNJ2BPsynaptojanin 2 binding protein−0.09
THBS1thrombospondin 1−0.09
CTGFconnective tissue growth factor−0.10
ITGB3integrin, beta 3 (platelet glycoprotein IIIa, antigen CD61)−0.12
AAMPangio-associated, migratory cell protein−0.12
GJA5gap junction protein, alpha 5, 40 kDa−0.12
PRKCBprotein kinase C, beta−0.13
EGR3early growth response 3−0.13
JMJD6jumonji domain containing 6−0.13
TGFBItransforming growth factor, beta-induced, 68 kDa−0.14
SIRT1sirtuin 1−0.14
ANGPTL3angiopoietin-like 3−0.14
ACKR3atypical chemokine receptor 3−0.14
SAT1spermidine/spermine N1-acetyltransferase 1−0.15
VEGFBvascular endothelial growth factor B−0.16
UTS2urotensin 2−0.16
JUNjun proto-oncogene−0.16
TNFSF12tumor necrosis factor (ligand) superfamily, member 12−0.16
EGFL7EGF-like-domain, multiple 7−0.17
MED1mediator complex subunit 1−0.17
SLIT2slit guidance ligand 2−0.17
SERPINF1serpin peptidase inhibitor, clade F (alpha-2 antiplasmin,−0.18
pigment epithelium derived factor), member 1
NOTCH3notch 3−0.18
FGF9fibroblast growth factor 9−0.19
DLL4delta-like 4 (Drosophila)−0.19
CCL2chemokine (C-C motif) ligand 2−0.19
MMP14matrix metallopeptidase 14 (membrane-inserted)−0.19
TMPRSS6transmembrane protease, serine 6−0.19
EPGNepithelial mitogen−0.20
RBPJrecombination signal binding protein for immunoglobulin−0.20
kappa J region
COL4A2collagen, type IV, alpha 2−0.20
PRKD2protein kinase D2−0.20
ALOX12arachidonate 12-lipoxygenase−0.21
RNH1ribonuclease/angiogenin inhibitor 1−0.21
APOHapolipoprotein H (beta-2-glycoprotein I)−0.21
CHI3L1chitinase 3-like 1 (cartilage glycoprotein-39)−0.21
ESM1endothelial cell-specific molecule 1−0.22
PTGS2prostaglandin-endoperoxide synthase 2 (prostaglandin G/H−0.22
synthase and cyclooxygenase)
ANPEPalanyl (membrane) aminopeptidase−0.22
LEMD3LEM domain containing 3−0.22
UTS2Rurotensin 2 receptor−0.22
CIB1calcium and integrin binding 1 (calmyrin)−0.22
ITGB1BP1integrin beta 1 binding protein 1−0.22
AQP1aquaporin 1 (Colton blood group)−0.22
IL18interleukin 18−0.22
EPHA2EPH receptor A2−0.22
EPHB1EPH receptor B1−0.22
AGTangiotensinogen (serpin peptidase inhibitor, clade A,−0.22
member 8)
PLAUplasminogen activator, urokinase−0.22
VEZF1vascular endothelial zinc finger 1−0.23
SPHK1sphingosine kinase 1−0.23
SRPX2sushi-repeat containing protein, X-linked 2−0.23
PDCL3phosducin-like 3−0.23
COL8A1collagen, type VIII, alpha 1−0.24
HDAC9histone deacetylase 9−0.24
CTSHcathepsin H−0.24
EDN1endothelin 1−0.24
CXCL8chemokine (C-X-C motif) ligand 8−0.24
ECM1extracellular matrix protein 1−0.24
BRCA1breast cancer 1, early onset−0.24
EFNB2ephrin-B2−0.25
SERPINE1serpin peptidase inhibitor, clade E (nexin, plasminogen−0.25
activator inhibitor type 1), member 1
SASH1SAM and SH3 domain containing 1−0.25
WNT7Bwingless-type MMTV integration site family, member 7B−0.25
RAMP1receptor (G protein-coupled) activity modifying protein 1−0.26
SCG2secretogranin II−0.26
COL8A2collagen, type VIII, alpha 2−0.26
SULF1sulfatase 1−0.26
CLIC4chloride intracellular channel 4−0.26
FGF1fibroblast growth factor 1 (acidic)−0.27
NODALnodal growth differentiation factor−0.27
RASIP1Ras interacting protein 1−0.28
RLN2relaxin 2−0.28
POFUT1protein O-fucosyltransferase 1−0.28
FGF18fibroblast growth factor 18−0.28
AIMP1aminoacyl tRNA synthetase complex-interacting−0.28
multifunctional protein 1
TGFBR2transforming growth factor, beta receptor II (70/80 kDa)−0.28
RHOBras homolog family member B−0.28
GBX2gastrulation brain homeobox 2−0.28
ENPP2ectonucleotide pyrophosphatase/phosphodiesterase 2−0.29
MAPK7mitogen-activated protein kinase 7−0.30
PROK2prokineticin 2−0.30
E2F7E2F transcription factor 7−0.30
ERAP1endoplasmic reticulum aminopeptidase 1−0.31
MTDHmetadherin−0.31
KLF5Kruppel-like factor 5 (intestinal)−0.31
DICER1dicer 1, ribonuclease type III−0.32
LECT1leukocyte cell derived chemotaxin 1−0.32
CX3CL1chemokine (C-X3-C motif) ligand 1−0.32
PTK2Bprotein tyrosine kinase 2 beta−0.33
SEMA4Asema domain, immunoglobulin domain (Ig),−0.34
transmembrane domain (TM) and short cytoplasmic
domain, (semaphorin) 4A
ARHGAP22Rho GTPase activating protein 22−0.34
RSPO3R-spondin 3−0.34
KLF4Kruppel-like factor 4 (gut)−0.34
ROBO1roundabout guidance receptor 1−0.34
GPLD1glycosylphosphatidylinositol specific phospholipase D1−0.35
NUS1NUS1 dehydrodolichyl diphosphate synthase subunit−0.35
NRARPNOTCH-regulated ankyrin repeat protein−0.35
PDCD10programmed cell death 10−0.36
PF4platelet factor 4−0.36
PRKXprotein kinase, X-linked−0.36
PMLpromyelocytic leukemia−0.36
ATP5BATP synthase, H+ transporting, mitochondrial F1 complex,−0.36
beta polypeptide
TNFRSF12Atumor necrosis factor receptor superfamily, member 12A−0.36
ENGendoglin−0.37
THY1Thy-1 cell surface antigen−0.37
FGF2fibroblast growth factor 2 (basic)−0.37
CXCL12chemokine (C-X-C motif) ligand 12−0.37
CAV1caveolin 1, caveolae protein, 22 kDa−0.38
PDGFAplatelet-derived growth factor alpha polypeptide−0.38
PNPLA6patatin-like phospholipase domain containing 6−0.38
PLCD3phospholipase C, delta 3−0.38
DDAH1dimethylarginine dimethylaminohydrolase 1−0.39
GNA13guanine nucleotide binding protein (G protein), alpha 13−0.39
ADM2adrenomedullin 2−0.39
HMOX1heme oxygenase 1−0.40
MCAMmelanoma cell adhesion molecule−0.41
RAPGEF3Rap guanine nucleotide exchange factor (GEF) 3−0.41
TNFAIP2tumor necrosis factor, alpha-induced protein 2−0.41
HTATIP2HIV-1 Tat interactive protein 2, 30 kDa−0.42
NCLnucleolin−0.42
ERBB2erb-b2 receptor tyrosine kinase 2−0.43
NAA15N(alpha)-acetyltransferase 15, NatA auxiliary subunit−0.43
ATPIF1ATPase inhibitory factor 1−0.43
THBS4thrombospondin 4−0.43
SYKspleen tyrosine kinase−0.44
LIFleukemia inhibitory factor−0.44
THBS2thrombospondin 2−0.44
PPP1R16Bprotein phosphatase 1, regulatory subunit 16B−0.44
NOTCH1notch 1−0.44
RUNX1runt-related transcription factor 1−0.45
PDCD6programmed cell death 6−0.45
VASH2vasohibin 2−0.45
GPIglucose-6-phosphate isomerase−0.46
ZC3H12Azinc finger CCCH-type containing 12A−0.46
WARStryptophanyl-tRNA synthetase−0.46
HYAL1hyaluronoglucosaminidase 1−0.47
PIK3CBphosphatidylinositol-4,5-bisphosphate 3-kinase, catalytic−0.47
subunit beta
TNMDtenomodulin−0.49

Example 13

In this example, genes in the extracellular matrix family that are expressed in human ventricular progenitor cells (HVPs) were identified. HVPs were generated as described in Examples 1 or 11 and RNA sequencing (RNA-seq) was performed at different time points following differentiation as described in Example 2. Cluster analysis of gene expression profiles at different time points during HVP differentiation identified stage-specific signature genes. These genes were clustered hierarchically on the basis of the similarity of their expression profiles. First, genes showing expression in four different categories were identified: (i) cell surface expression; (ii) co-expression with Islet 1; (iii) high expression on day 5 of differentiation; and (iv) high d5/d0 ratio. This analysis confirmed the cell surface markers for HVPs of: JAG1, FZD4, FGFR3, LIFR (CD118) and TNFSF9. Next, from this same population of HVPs that identified the cell surface markers, gene ontogeny searches were performed to identify extracellular matrix family genes that were expressed in this population of HVPs, to thereby identify a gene fingerprint profile that identifies genes critical for cell engraftment.

Statistically, Pearson's correlation with Isl1 expression was used to identify those extracellular matrix genes whose expression in the HVPs best correlated with Isl1 expression. Table 5 below lists the extracellular matrix genes that correlate with Isl1 expression with a Pearson's correlation of 0.50 or higher.

TABLE 5
Extracellular matrix genes expressed in HVPs with a Pearson
Correlation with Isl1 Expression of 0.50 or greater
Pearson's
Correlation with
GeneExtracellular matrix genes (GO:0031012)Isl1 Expression
FGF10fibroblast growth factor 100.98
SMOC1SPARC related modular calcium binding 10.97
CCBE1collagen and calcium binding EGF domains 10.94
COL6A6collagen, type VI, alpha 60.89
ADAMTS12ADAM metallopeptidase with thrombospondin type 10.85
motif, 12
COL19A1collagen, type XIX, alpha 10.85
LAMA1laminin, alpha 10.85
BMP4bone morphogenetic protein 40.84
FBLN7fibulin 70.81
FBLN2fibulin 20.81
NDNFneuron-derived neurotrophic factor0.80
HTRA1HtrA serine peptidase 10.80
HAPLN1hyaluronan and proteoglycan link protein 10.79
EMILIN1elastin microfibril interfacer 10.79
SPOCK3sparc/osteonectin, cwcv and kazal-like domains0.76
proteoglycan (testican) 3
PODNL1podocan-like 10.73
IHHindian hedgehog0.70
ACANaggrecan0.69
NID2nidogen 2 (osteonidogen)0.69
COL4A6collagen, type IV, alpha 60.68
LAMC1laminin, gamma 1 (formerly LAMB2)0.65
FMODfibromodulin0.65
MUC4mucin 4, cell surface associated0.64
EMID1EMI domain containing 10.62
HMCN1hemicentin 10.61
NID1nidogen 10.60
VCANversican0.58
CILP2cartilage intermediate layer protein 20.57
SOD3superoxide dismutase 3, extracellular0.56
ADAMTS3ADAM metallopeptidase with thrombospondin type 10.54
motif, 3
ZP3zona pellucida glycoprotein 3 (sperm receptor)0.54
ANGPTL4angiopoietin-like 40.53
CRTAC1cartilage acidic protein 10.52
LTBP4latent transforming growth factor beta binding protein 40.50
FREM1FRAS1 related extracellular matrix 10.50
Table 6 below lists the extracellular matrix genes that correlate with Isl1 expression with a Pearson's correlation of 0.49-0.00.

TABLE 6
Extracellular matrix genes expressed in HVPs with a Pearson
Correlation with Isl1 Expression of 0.49 to 0.00
Pearson's
Correlation with
GeneExtracellular matrix genes (GO:0031012)Isl1 Expression
SSC5Dscavenger receptor cysteine rich family, 5 domains0.49
GPC6glypican 60.49
COL1A1collagen, type I, alpha 10.49
ADAMTSL3ADAMTS-like 30.48
FLRT3fibronectin leucine rich transmembrane protein 30.48
FBLN1fibulin 10.48
ADAMTS9ADAM metallopeptidase with thrombospondin type 10.48
motif, 9
COL27A1collagen, type XXVII, alpha 10.47
RELNreelin0.46
COL9A2collagen, type IX, alpha 20.46
EFEMP2EGF containing fibulin-like extracellular matrix protein 20.45
AGRNagrin0.44
PCOLCEprocollagen C-endopeptidase enhancer0.44
NTN4netrin 40.44
CD248CD248 molecule, endosialin0.44
TGFB1transforming growth factor, beta 10.43
ADAMTS2ADAM metallopeptidase with thrombospondin type 10.43
motif, 2
CTHRC1collagen triple helix repeat containing 10.42
FGFR2fibroblast growth factor receptor 20.41
APOEapolipoprotein E0.41
MMP11matrix metallopeptidase 110.41
MMP15matrix metallopeptidase 15 (membrane-inserted)0.41
PODNpodocan0.39
VEGFAvascular endothelial growth factor A0.39
COL18A1collagen, type XVIII, alpha 10.39
GLG1golgi glycoprotein 10.39
GPC2glypican 20.37
DAG1dystroglycan 1 (dystrophin-associated glycoprotein 1)0.35
TGFB2transforming growth factor, beta 20.35
PRELPproline/arginine-rich end leucine-rich repeat protein0.35
CHADchondroadherin0.33
COL2A1collagen, type II, alpha 10.33
FN1fibronectin 10.31
SMC3structural maintenance of chromosomes 30.31
COL4A5collagen, type IV, alpha 50.30
FBN3fibrillin 30.30
MMP23Bmatrix metallopeptidase 23B0.30
CCDC80coiled-coil domain containing 800.29
SPARCsecreted protein, acidic, cysteine-rich (osteonectin)0.29
TNXBtenascin XB0.28
COL6A2collagen, type VI, alpha 20.28
ADAMTS13ADAM metallopeptidase with thrombospondin type 10.28
motif, 13
LOXL1lysyl oxidase-like 10.28
HAPLN2hyaluronan and proteoglycan link protein 20.28
TNCtenascin C0.28
ENTPD2ectonucleoside triphosphate diphosphohydrolase 20.28
TGFB3transforming growth factor, beta 30.28
MFAP4microfibrillar-associated protein 40.27
VWFvon Willebrand factor0.27
WNT2wingless-type MMTV integration site family member 20.27
MMRN2multimerin 20.27
SPON1spondin 1, extracellular matrix protein0.26
ADAMTS1ADAM metallopeptidase with thrombospondin type 10.26
motif, 1
F2coagulation factor II (thrombin)0.26
FLRT2fibronectin leucine rich transmembrane protein 20.25
MMP2matrix metallopeptidase 20.24
COL26A1collagen, type XXVI, alpha 10.24
CASKcalcium/calmodulin-dependent serine protein kinase0.24
(MAGUK family)
NTN3netrin 30.23
SLC1A3solute carrier family 1 (glial high affinity glutamate0.22
transporter), member 3
F3coagulation factor III (thromboplastin, tissue factor)0.20
ADAMTS6ADAM metallopeptidase with thrombospondin type 10.20
motif, 6
COL5A2collagen, type V, alpha 20.19
ERBB2IPerbb2 interacting protein0.18
LAMB1laminin, beta 10.18
COLQcollagen-like tail subunit (single strand of homotrimer) of0.18
asymmetric acetylcholinesterase
LAMA5laminin, alpha 50.17
LOXL2lysyl oxidase-like 20.17
WNT11wingless-type MMTV integration site family, member 110.17
LAMB2laminin, beta 2 (laminin S)0.17
COL5A1collagen, type V, alpha 10.17
AEBP1AE binding protein 10.17
COL9A3collagen, type IX, alpha 30.16
CTSDcathepsin D0.16
COL21A1collagen, type XXI, alpha 10.16
EGFLAMEGF-like, fibronectin type III and laminin G domains0.16
FBN2fibrillin 20.15
NAV2neuron navigator 20.15
EMILIN2elastin microfibril interfacer 20.14
WNT9Bwingless-type MMTV integration site family, member 9B0.14
NOVnephroblastoma overexpressed0.14
CHL1cell adhesion molecule L1-like0.13
DLG1discs, large homolog 1 (Drosophila)0.11
MFGE8milk fat globule-EGF factor 8 protein0.11
TIMP1TIMP metallopeptidase inhibitor 10.11
CST3cystatin C0.10
APLP1amyloid beta (A4) precursor-like protein 10.10
PRTN3proteinase 30.10
ADAMTS10ADAM metallopeptidase with thrombospondin type 10.09
motif, 10
ILKintegrin-linked kinase0.09
FRAS1Fraser extracellular matrix complex subunit 10.09
ANXA2P2annexin A2 pseudogene 20.08
SMOC2SPARC related modular calcium binding 20.07
ANXA2annexin A20.06
ODAModontogenic, ameloblast asssociated0.06
FREM2FRAS1 related extracellular matrix protein 20.05
HAPLN3hyaluronan and proteoglycan link protein 30.05
GPC3glypican 30.03
LGALS1lectin, galactoside-binding, soluble, 10.02
ADAMTS8ADAM metallopeptidase with thrombospondin type 10.02
motif, 8
LUMlumican0.01
HSP90B1heat shock protein 90 kDa beta (Grp94), member 10.00
HAPLN4hyaluronan and proteoglycan link protein 40.00
MATN2matrilin 20.00

Extracellular matrix genes whose expression negatively correlated with Isl1 expression in the HVPs were also identified. Table 7 below lists the extracellular matrix genes that negatively correlate with Isl1 expression with a Pearson's correlation of −0.50 or less.

TABLE 7
Extracellular matrix genes expressed in HVPs with a Pearson
Correlation with Isl1 Expression of −0.50 or less
Pearson's
Correlation with
GeneExtracellular matrix genes (GO:0031012)Isl1 Expression
FKBP1AFK506 binding protein 1A, 12 kDa−0.51
CLUclusterin−0.52
TFPI2tissue factor pathway inhibitor 2−0.52
PLSCR1phospholipid scramblase 1−0.53
FBLN5fibulin 5−0.53
VWA1von Willebrand factor A domain containing 1−0.54
ADAMTS16ADAM metallopeptidase with thrombospondin type 1−0.55
motif, 16
MMP25matrix metallopeptidase 25−0.55
SFRP2secreted frizzled-related protein 2−0.61
SOD1superoxide dismutase 1, soluble−0.68
Table 8 below lists the extracellular matrix genes that negatively correlate with Isl1 expression with a Pearson's correlation of −0.01 to −0.49.

TABLE 8
Extracellular matrix genes expressed in HVPs with a Pearson
Correlation with Isl1 Expression of −0.01 to −0.49
Pearson's
Correlation with
GeneExtracellular matrix genes (GO:0031012)Isl1 Expression
PAPLNpapilin, proteoglycan-like sulfated glycoprotein−0.01
SOSTsclerostin−0.01
CDONcell adhesion associated, oncogene regulated−0.02
HMCN2hemicentin 2−0.02
WNT5Awingless-type MMTV integration site family, member 5A−0.02
PCSK6proprotein convertase subtilisin/kexin type 6−0.02
GSTO1glutathione S-transferase omega 1−0.02
LTBP1latent transforming growth factor beta binding protein 1−0.03
KAZALD1Kazal-type serine peptidase inhibitor domain 1−0.03
LTBP2latent transforming growth factor beta binding protein 2−0.03
SFRP1secreted frizzled-related protein 1−0.04
ADAM11ADAM metallopeptidase domain 11−0.05
COL6A1collagen, type VI, alpha 1−0.05
COL22A1collagen, type XXII, alpha 1−0.05
CYR61cysteine-rich, angiogenic inducer, 61−0.05
ELNelastin−0.06
COL9A1collagen, type IX, alpha 1−0.06
VTNvitronectin−0.06
COL15A1collagen, type XV, alpha 1−0.06
COL4A1collagen, type IV, alpha 1−0.07
ANGangiogenin, ribonuclease, RNase A family, 5−0.07
HSPG2heparan sulfate proteoglycan 2−0.09
CRIP2cysteine-rich protein 2−0.09
CD151CD151 molecule (Raph blood group)−0.09
THBS1thrombospondin 1−0.09
ADAMTS4ADAM metallopeptidase with thrombospondin type 1−0.09
motif, 4
CTGFconnective tissue growth factor−0.10
CRISPLD2cysteine-rich secretory protein LCCL domain containing 2−0.10
BMP7bone morphogenetic protein 7−0.11
COL6A3collagen, type VI, alpha 3−0.11
COL3A1collagen, type III, alpha 1−0.11
COL14A1collagen, type XIV, alpha 1−0.11
MATN3matrilin 3−0.11
CPZcarboxypeptidase Z−0.11
BMP1bone morphogenetic protein 1−0.11
WISP1WNT1 inducible signaling pathway protein 1−0.12
ADAMTS18ADAM metallopeptidase with thrombospondin type 1−0.12
motif, 18
COL7A1collagen, type VII, alpha 1−0.12
IGFBP7insulin-like growth factor binding protein 7−0.12
COCHcochlin−0.13
ADAMTS5ADAM metallopeptidase with thrombospondin type 1−0.13
motif, 5
COL11A2collagen, type XI, alpha 2−0.13
TGFBItransforming growth factor, beta-induced, 68 kDa−0.14
COL16A1collagen, type XVI, alpha 1−0.14
ACHEacetylcholinesterase (Yt blood group)−0.14
THSD4thrombospondin, type I, domain containing 4−0.15
DGCR6DiGeorge syndrome critical region gene 6−0.15
TGFB1I1transforming growth factor beta 1 induced transcript 1−0.15
ADAMTSL1ADAMTS-like 1−0.15
SERPINA1serpin peptidase inhibitor, clade A (alpha-1 antiproteinase,−0.16
antitrypsin), member 1
MAMDC2MAM domain containing 2−0.16
LAMA4laminin, alpha 4−0.17
LTBP3latent transforming growth factor beta binding protein 3−0.17
EGFL7EGF-like-domain, multiple 7−0.17
NPNTnephronectin−0.17
SERPINF1serpin peptidase inhibitor, clade F (alpha-2 antiplasmin,−0.18
pigment epithelium derived factor), member 1
ABI3BPABI family, member 3 (NESH) binding protein−0.18
SERPINE2serpin peptidase inhibitor, clade E (nexin, plasminogen−0.18
activator inhibitor type 1), member 2
WNT6wingless-type MMTV integration site family, member 6−0.19
TIMP3TIMP metallopeptidase inhibitor 3−0.19
SNCAsynuclein, alpha (non A4 component of amyloid precursor)−0.19
PKMpyruvate kinase, muscle−0.19
FGF9fibroblast growth factor 9−0.19
VITvitrin−0.19
WNT1wingless-type MMTV integration site family, member 1−0.19
LAMC3laminin, gamma 3−0.19
MMP14matrix metallopeptidase 14 (membrane-inserted)−0.19
PXDNperoxidasin−0.19
HNRNPMheterogeneous nuclear ribonucleoprotein M−0.19
FBN1fibrillin 1−0.20
ASPNasporin−0.20
ADAMTSL5ADAMTS-like 5−0.20
SPON2spondin 2, extracellular matrix protein−0.20
COL1A2collagen, type I, alpha 2−0.20
BGNbiglycan−0.20
COL4A2collagen, type IV, alpha 2−0.20
ADAMTSL4ADAMTS-like 4−0.21
APOHapolipoprotein H (beta-2-glycoprotein I)−0.21
CHI3L1chitinase 3-like 1 (cartilage glycoprotein-39)−0.21
ADAMTS7ADAM metallopeptidase with thrombospondin type 1−0.22
motif, 7
CALRcalreticulin−0.22
MMP9matrix metallopeptidase 9−0.22
MMP24matrix metallopeptidase 24 (membrane-inserted)−0.22
SPOCK2sparc/osteonectin, cwcv and kazal-like domains−0.22
proteoglycan (testican) 2
COL11A1collagen, type XI, alpha 1−0.23
MMP7matrix metallopeptidase 7−0.23
MMP16matrix metallopeptidase 16 (membrane-inserted)−0.23
MFAP2microfibrillar-associated protein 2−0.23
POSTNperiostin, osteoblast specific factor−0.24
COL8A1collagen, type VIII, alpha 1−0.24
WNT2Bwingless-type MMTV integration site family, member 2B−0.24
DCNdecorin−0.24
EGFL6EGF-like-domain, multiple 6−0.24
MMP10matrix metallopeptidase 10−0.24
MGPmatrix Gla protein−0.24
ECM1extracellular matrix protein 1−0.24
SERPINE1serpin peptidase inhibitor, clade E (nexin, plasminogen−0.25
activator inhibitor type 1), member 1
MMP1matrix metallopeptidase 1−0.25
WNT10Awingless-type MMTV integration site family, member 10A−0.25
B4GALT7xylosylprotein beta 1,4-galactosyltransferase, polypeptide 7−0.25
COL12A1collagen, type XII, alpha 1−0.25
LAMA3laminin, alpha 3−0.25
LAMA2laminin, alpha 2−0.25
LAMB3laminin, beta 3−0.25
WNT7Bwingless-type MMTV integration site family, member 7B−0.25
FLRT1fibronectin leucine rich transmembrane protein 1−0.25
ADAMTS15ADAM metallopeptidase with thrombospondin type 1−0.26
motif, 15
COL8A2collagen, type VIII, alpha 2−0.26
MFAP1microfibrillar-associated protein 1−0.26
TINAGL1tubulointerstitial nephritis antigen-like 1−0.26
FGF1fibroblast growth factor 1 (acidic)−0.27
OLFML2Aolfactomedin-like 2A−0.27
CPA6carboxypeptidase A6−0.27
COL17A1collagen, type XVII, alpha 1−0.27
SPARCL1SPARC-like 1 (hevin)−0.27
MFAP5microfibrillar associated protein 5−0.27
COL4A4collagen, type IV, alpha 4−0.28
WNT8Bwingless-type MMTV integration site family, member 8B−0.28
ADAMTS19ADAM metallopeptidase with thrombospondin type 1−0.29
motif, 19
CRTAPcartilage associated protein−0.29
WNT5Bwingless-type MMTV integration site family, member 5B−0.30
WNT3wingless-type MMTV integration site family, member 3−0.30
UCMAupper zone of growth plate and cartilage matrix associated−0.30
GPC1glypican 1−0.30
TIMP2TIMP metallopeptidase inhibitor 2−0.30
ALPLalkaline phosphatase, liver/bone/kidney−0.30
LECT1leukocyte cell derived chemotaxin 1−0.32
GPC4glypican 4−0.32
SPOCK1sparc/osteonectin, cwcv and kazal-like domains−0.32
proteoglycan (testican) 1
HSD17B12hydroxysteroid (17-beta) dehydrogenase 12−0.32
LGALS3lectin, galactoside-binding, soluble, 3−0.33
EMILIN3elastin microfibril interfacer 3−0.34
GFOD2glucose-fructose oxidoreductase domain containing 2−0.34
VWC2von Willebrand factor C domain containing 2−0.34
SERAC1serine active site containing 1−0.34
WNT8Awingless-type MMTV integration site family, member 8A−0.34
LMCD1LIM and cysteine-rich domains 1−0.34
CPXM2carboxypeptidase X (M14 family), member 2−0.34
ADAMTS14ADAM metallopeptidase with thrombospondin type 1−0.34
motif, 14
GPLD1glycosylphosphatidylinositol specific phospholipase D1−0.35
FGFBP3fibroblast growth factor binding protein 3−0.35
BCANbrevican−0.35
ITGB4integrin, beta 4−0.35
LGALS3BPlectin, galactoside-binding, soluble, 3 binding protein−0.36
LPLlipoprotein lipase−0.38
LAD1ladinin 1−0.39
WNT3Awingless-type MMTV integration site family, member 3A−0.39
TGFBR3transforming growth factor, beta receptor III−0.39
DSTdystonin−0.40
WNT10Bwingless-type MMTV integration site family, member 10B−0.40
LEFTY2left-right determination factor 2−0.41
TNFRSF11Btumor necrosis factor receptor superfamily, member 11b−0.41
WNT9Awingless-type MMTV integration site family, member 9A−0.41
TIMP4TIMP metallopeptidase inhibitor 4−0.42
WNT4wingless-type MMTV integration site family, member 4−0.42
NCANneurocan−0.42
ADAMTS20ADAM metallopeptidase with thrombospondin type 1−0.43
motif, 20
ITGA6integrin, alpha 6−0.43
LOXlysyl oxidase−0.43
THBS4thrombospondin 4−0.43
THBS2thrombospondin 2−0.44
ADAMTSL2ADAMTS-like 2−0.44
ENTPD1ectonucleoside triphosphate diphosphohydrolase 1−0.45
RUNX1runt-related transcription factor 1−0.45
VWA2von Willebrand factor A domain containing 2−0.45
RELL2RELT-like 2−0.46
PTPRZ1protein tyrosine phosphatase, receptor-type, Z polypeptide 1−0.46
LAMC2laminin, gamma 2−0.46
DSTdystonin−0.40
WNT10Bwingless-type MMTV integration site family, member 10B−0.40
LEFTY2left-right determination factor 2−0.41
TNFRSF11Btumor necrosis factor receptor superfamily, member 11b−0.41
WNT9Awingless-type MMTV integration site family, member 9A−0.41
TIMP4TIMP metallopeptidase inhibitor 4−0.42

Example 14

In this example, the gene expression profile was determined for Islet 1 negative cells within the Day 6 HVP population to further characterize a subpopulation of cells within the Day 6 population that do not express the necessary markers to qualify as engraftable HVPs. Day 6 HVP populations were generated as described in Examples 1 or 11 and RNA sequencing (RNA-seq) was performed following differentiation as described in Example 2. Cells that were Islet 1 negative (Isl1−) were further analyzed with respect to their gene expression profile. Genes expressed in the Isl1− cells with an average RNA copy number of 2000 or higher are shown below in Table 9.

TABLE 9
Gene Expression Profile of Day 6 Islet 1 Negative Cells
GeneSample #1Sample #2Avg. RNA Copy #
ACTB122882812620207
MTRNR2L224511977417142.5
MALAT1141631809216127.5
EEF1A1116631245612059.5
KRT888841408711485.5
MTRNR2L812836768810262
KRT185215105527883.5
FN14900105817740.5
MTRNR2L1855057197134.5
TTN314996016375
GAPDH490763915649
YWHAZ534954145381.5
MTRNR2L9567338504761.5
RPL3334659114628.5
AHNAK672721974462
KCNQ1OT1583530644449.5
TUBB487039364403
SLC2A3331143013806
FTL348440213752.5
HSP90B1417327783475.5
KRT19350232023352
HSPA8345529033179
MYL6189843753136.5
RPLP0231939223120.5
BSG251935933056
COL3A153126953003.5
TPM1293830592998.5
VCAN256334222992.5
ENO1244935352992
RPL4261933282973.5
ACTG1268732532970
MTRNR2L10340924872948
HMGN2268431532918.5
PRTG259429802787
TPI1241831132765.5
HMGB1257728802728.5
VIM262127042662.5
ATP5B300022192609.5
HSP90AB1273524192577
RPL7213228962514
CBX5279922192509
MYL7161433822498
SERPINH1254723272437
HNRNPK287819322405
SRRM2275820462402
PODXL368311122397.5
EEF2257921192349
SPARC302616452335.5
ACTC143741522294.5
HUWE1258319772280
COL1A235449412242.5
LINC00506296514962230.5
HSPA5207823562217
MDK222321442183.5
HNRNPC229220742183
HSP90AA1222021382179
RGS5218021502165
LAMC1275715652161
APLNR86832462057
UGDH-AS1263314572045
RPS3A160123992000
Accordingly, the data shown in Table 9 provides a gene expression profile for Islet 1 negative, non-engraftable cells within a Day 6 HVP population that are not suitable for transplantation and thus are to be selected against when choosing cells for transplantation and engraftment.

Free full text: Click here