Relevance of iPSC-derived human PGC-like cells at the surface of embryoid bodies to prechemotaxis migrating PGCs
Pluripotent stem cell-derived human primordial germ cell-like cells (hPGCLCs) provide important opportunities to study primordialgerm cells (PGCs). We robustly produced CD38+ hPGCLCs [∼43% ofFACS-sorted embryoid body (EB) cells] from primed-state inducedpluripotent stem cells (iPSCs) after a 72-hour transient incubation in the four chemical inhibitors (4i)-naïve reprogramming medium and showed transcriptional consistency of our hPGCLCs withhPGCLCs generated in previous studies using various and distinct protocols. Both CD38+ hPGCLCs and CD38− EB cells significantly expressed PRDM1 and TFAP2C, although PRDM1 mRNA in CD38− cells lacked the 3′-UTR harboring miRNA binding sites regulating mRNA stability. Genes up-regulated in hPGCLCs were enriched forcell migration genes, and their promoters were enriched for the binding motifs of TFAP2 (which was identified in promoters of T, NANOS3, and SOX17) and the RREB-1 cell adhesion regulator. In EBs, hPGCLCs were identified exclusively in the outermost surface monolayer as dispersed cells or cell aggregates with strong and specific expression of POU5F1/OCT4 protein. Time-lapse live cell imaging revealed active migration of hPGCLCs on Matrigel. Whereas all hPGCLCs strongly expressed the CXCR4 chemotaxis receptor, its ligand CXCL12/SDF1 was not significantly expressed in the whole EBs. Exposure of hPGCLCs to CXCL12/SDF1 induced cell migration genes and antiapoptosis genes.
Thus, our study shows that tran- scriptionally consistent hPGCLCs can be readily produced from hiPSCs after transition of their pluripotency from the primed state using various methods and that hPGCLCs resemble the early-stage PGCs randomly migrating in the midline region of human embryos before initiation of the CXCL12/SDF1-guided chemotaxis.Pluripotency is the cellular capability to differentiate into primordial germ cells (PGCs) and three embryonic germ layers (1). Studies using mouse models have revealed that plu- ripotency in vivo is only a transient state gained by cells of the inner cell mass (ICM) and epiblast during the early stages of embryonic development. However, mouse ICM cells and epiblast cells can be maintained indefinitely without losing their pluri- potency in cell culture as ES cells (ESCs) and epiblast stem cells, respectively, and engineered epigenetic reprogramming of differ- entiated cells can gain pluripotency to generate induced pluripo- tent stem cells (iPSCs). Whereas human pluripotent stem cells (hPSCs) maintained in the mTeSR1 medium on Matrigel has a primed-state pluripotency (1), Gafni et al. (2) showed that it can be converted into a naïve state when grown in the naïve human stem cell medium (NHSM). Although transgene-free conversion of primed hPSCs to naïve state was subsequently achieved under several different cell culture conditions, features of the achieved naïve pluripotency varied significantly (1, 3), and long-term maintenance of pluripotent stem cells (PSCs) in a strongly naïvestate may increase risk of chromosomal anomalies (1, 4).Although PSCs can produce cells resembling PGCs [primor- dial germ cell-like cells (PGCLCs)] in cell culture (5–14), direct conversion of mouse PSCs to PGCLCs is inefficient (6). Hayashi et al. (6) generated cells resembling early epiblast [epiblast-likecells (EpiLCs)] from mouse PSCs and showed their strong ability to produce PGCLCs as cell aggregates in the presence of bone morphogenetic protein 4 (Bmp4).
In the absence of Bmp4, ro- bust PGCLC production can also be achieved by overexpression of three transcription factors (TFs)—namely Prdm1, Prdm14,and Tfap2c (7, 14)—or by a single TF Nanog (8) from exogenousvectors. The high competency of EpiLCs for PGCLC productionis transient, peaking at day 2 of EpiLC induction but largely lost at day 3 (6, 10). In a previous study (10), we showed robust and global demethylation of genomic DNA (gDNA) involving era- sure of DNA methylation at CpG dinucleotides in the imprinting control regions.In the first report of efficient human primordial germ cell-likecell (hPGCLC) production from PSCs, Irie et al. (15) maintained human induced pluripotent stem cells (hiPSCs) in an ERK- independent naïve pluripotency state for at least 2 wk in amodified NHSM containing LIF, FGF2, and TGF-β1 and four chemical inhibitors (4i) of MEK, p38MAPK, GSK3β, and JNK. Aggregates of hPSCs maintained in the 4i medium robustlydifferentiated into hPGCLCs, which strongly express NANOS3, TNAP, and CD38 (15). Efficiency of hPGCLC production from the conventional primed pluripotency hiPSCs was very low (<5%) (15), which was later confirmed by Sasaki et al. (16). Including Irie et al. (15), at least four laboratories reported ef-ficient hPGCLC production from PSCs (9, 15–18) but with sig- nificant variations in the intermediate cell cultures and markerantigens used for FACS enrichment of hPGCLCs (9, 16, 17) (Table S1).In this study, we show that a 72-h exposure of the primed pluripotency hiPSCs in the 4i medium is sufficient for a robust production of CD38+ hPGCLCs. In contrast to mouse germ cell development, induction of PRDM1 or TFAP2C did not seem to be the primary determinant of hiPSC differentiation to hPGCLCs, agreeing with observations made by Irie et al. (15). hiPSC differentiation to hPGCLCs in embryoid bodies (EBs) was associated with enriched induction of genes involved in cell migration, and most hPGCLCs were observed at the outermost surface monolayer of EBs. Live cell imaging revealed actively migrating hPGCLCs forming cellular protrusions. All hPGCLCs expressed the CXCR4 chemotaxis receptor, whereas its ligand CXCL12/SDF1 was not significantly expressed in any cells in EBs. Exposure of hPGCLCs to CXCL12/SDF1 induced genes involved in cell migration or antiapoptosis. These results suggestthat hPGCLCs in EBs resemble early-stage PGCs randomly migrating in the midline region of human embryos before initi- ation of their directional migration (i.e., chemotaxis) toward genital ridges under the CXCR4-CXCL12 signaling.
Results
Production of CD38+ hPGCLCs from Short-Term 4i Reprogrammed hiPSCs. Since previous studies showed robust production of hPGCLCs via various precursor cell cultures (Table S1), wespeculated that the primed pluripotency state may specifically prevent hPSCs from germline differentiation, while various de- grees of deviation from it might be more permissive. To examine this possibility, we exposed primed pluripotency hiPSCs (clone A4; 46 + XY diploid) to the 4i medium for total 72 h (48-h exposure as monolayer cultures followed by 24-h exposure as EBs) and attempted to produce hPGCLCs using the protocol described by Irie et al. (15) (Fig. 1A). Colonies of A4 iPSCs grown in the mTeSR1 medium on feeder cells showed typical morphology of primed pluripotency hPSCs (Fig. 1B, Left), whereas colonies grown for 48 h in the 4i medium already showed a densely packed, sharp-edged appearance (Fig. 1B, Center), which agrees with the reported characteristics of naïve hPSCs (2). After 72 h of growth in the 4i medium, expression of DNMT3A and DNMT3B in A4 iPSCs dramatically decreased, whereas expression of TET1 or DNMT1 was unchanged or only modestly suppressed, respectively (Fig. 1C, Left). Strong sup- pression of DNMT3B expression was obvious after only 24-hincubation in the 4i medium (Fig. 1C, Right). The rapid changes in expression of the DNMT genes agreed with the characteristics of hiPSCs grown in the NHSM observed by Gafni et al. (2). These results support that a 72-h incubation in the 4i medium was sufficient to significantly deviate the pluripotency of hiPSCs out of the standard primed state maintained in the mTeSR1 medium. Because both the NHSM and 4i media commonlycontain 1 μM PD0325901 ERK1/2 inhibitor (2, 15), we speculate that the pluripotency of hiPSCs after a 72-h 4i culture may beshifted toward the ERK-independent naïve state, although exact characteristics are unknown.
After a 48-h culture in the 4i medium, we casted hiPSCs into the AggreWell microwells for rapid formation of EBs using the spin EB method (19). EBs were formed in the 4i medium within 24 h and then incubated in the PGCLC medium for 5–8 d (Fig. 1B, Right) followed by single-cell dissociation and FACS en- richment of hPGCLCs as CD38+ cells. FACS contour plotsshowed a readily discerned, distinct cell population of CD38+ cells, which corresponded to ∼6 and ∼43% of FACS-sorted cells after 5 and 8 d of EB culture in the hPGCLC production me- dium, respectively (Fig. 1D), whereas CD38+ mRNA was scarcely detected in day 0 EB (Fig. 1C). Real-time qPCR of the total EB cell lysates revealed that expression of PRDM1, an earlyPGC marker gene that encodes a TF required for PGC specifi- cation (15, 20, 21), was induced by day 2 and reached maximalstrength by day 4 (Fig. 1E). Expression of three other early PGC marker genes, DPPA3, NANOS3, and CD38, strongly increased from day 4 to 5, whereas their expression at day 0 was scarcely detectable (Fig. 1E). Expression of DAZL, a marker of naïve pluripotency as well as late-stage germline development (22), was rapidly suppressed in EBs, agreeing with previously pub- lished studies using other methods of hPGCLC production (16, 20). These results indicate that a 72-h reprogramming of primed pluripotency hiPSCs in the 4i medium is sufficient for subsequent robust production of CD38+ hPGCLCs.Transcriptomal Profiling of CD38+ hPGCLCs Produced via 72-h Pre-EB 4i Reprogramming.
To obtain insight into mechanisms of germline differentiation of hiPSCs, we determined transcriptomes of the CD38+ hPGCLCs and CD38− cells in day 5 EBs and their pre- cursors. Unsupervised hierarchical clustering of transcriptomes (Fig. 2A and Fig. S1 A and B) identified genes significantly expressed in both CD38+ and CD38− day 5 EB cells but only weakly expressed in their precursors (cluster 1; 81 genes), genes expressed predominantly in CD38+ cells over other cell types (cluster 2; 130 genes), and genes suppressed only in CD38− cells (cluster 3; 153 genes). Positions of these clusters are indicated on the right in Fig. 2A, and their enlarged images with full lists of genes are shown in Fig. S1C. A heat map of marker genes se- lected from Fig. 2A is also shown as Fig. 2B.Cluster 1 genes included PRDM1 and naïve pluripotency/ICM markers TCL1B and TFCP2L1 (16, 20). TFAP2C and another ICM marker KLF4 (16, 20) were also expressed in both CD38+ and CD38− EB cells relatively strongly, although varying degrees of weaker expression of TFAP2C were observed with some of the precursor hPSC cells (Fig. 2 A and B). The weak expression ofTFCP2L1 and KLF4 in both the primed hiPSCs and the short- term 4i reprogrammed hiPSCs was consistent with the previously reported characteristics of the primed pluripotency hPSCs and the ERK-independent naïve pluripotency hiPSCs, respectively (2, 3, 15), indicating that robust production of CD38+ hPGCLCs does not require strong expression of these markers of naïve pluripotency in the pre-EB precursor cultures. Instead, a de- viation from the primed pluripotency achieved after a 72-h in- cubation in the 4i medium seems sufficient. However, the ICM markers TCL1B, TFCP2L1, and KLF4 were strongly expressedin CD38− cells in day 5 EBs, suggesting a certain degree of commonality in the gene regulation network between ICM and EB cells incubated in the hPGCLCs. Cluster 2 genes includedknown markers of human primordial germ cells (hPGCs)/ hPGCLCs, such as NANOS3, SOX15, SOX17, and CD38 (15, 22,23), as well as mesodermal markers T, HOXD1, and WNT3 and a naïve pluripotency marker TCL1A (16).
Cluster 3 genes were enriched with markers shared by hPSCs and hPGCs/hPGCLCs, such as ALPL, KIT, NANOG, UTF1, and POU5F1 (15, 16). Manypluripotency markers (but not hPGC/hPGCLC markers), such asLIN28A, SALL4, DPPA4, PRDM14, ZIC3, and SOX2, werestrongly expressed in the pluripotent precursor cells (including day0 EB cells) but only weakly expressed in day 5 EBs. Tran- scriptomal principal component analysis (PCA) (Fig. 2C) showed reproducibility of our protocol, although EB cell differentiation to nongermline lineages seemed somewhat stochastic.The marker gene expression profile of our CD38+ hPGCLCs was nearly identical to those of the CD38+/TNAP+ hPGCLCs produced by Irie et al. (15) from hiPSCs grown in the 4i medium for at least 2 wk, confirming that the CD38+ cells produced in this study via a 72-h 4i reprogramming were bona fide hPGCLCs. Choi et al. (4) recently reported (together with J.O. and T.S.) that there is a significant risk of chromosomal aneuploidy and anomalous DNA demethylation in mouse PSCs after prolonged maintenance (≥20 passages) in the ERK1/2-independent naïvepluripotency state; our modification of hPGCLC productionprotocol using the ERK1/2-inhibiting 4i medium for only a limited period of time may be beneficial to reduce the risk of genetic mutations and epimutations, although the effect of prolonged ERK1/2-repressing culture on hPSCs remains topolymorphism of PRDM1 mRNA and/or let-7 miRNA in hiPSC differentiation to hPGCLCs remains to be determined.Gene Ontology (GO) analysis of the differentially expressed genes (DEGs) between CD38+ and CD38− EB cells [1,445 genes; false discovery rate (FDR) < 5%; more than fourfold change]revealed very strong enrichment of genes involved in cell migra- tion (Fig. 2E, Left) and plasma membrane activities (Fig. 2E, Right), both of which included formation of cell projections.
Most of the DEGs supporting these enrichments were expressed more strongly in CD38+ hPGCLCs than CD38− EB cells. These resultssuggest that the germline differentiation of EB cells may affecttheir ability to change their shape for migration, agreeing with the reported highly migratory feature of hPGCs (29).We attempted to determine whether CD38+ hPGCLCs in day 8 EBs reflect a more advanced stage of germline development than day 5 hPGCLCs (Fig. S3). Transcriptomal PCA (Fig. S3C) and hierarchical clustering of 3,500 DEGs with the lowest FDRs (Fig. S3D) revealed a clear distinction between days 5 and 8 CD38+ hPGCLCs. However, GO analysis of these DEGs did not reveal any meaningful enrichment. Expression of marker genes for pluripotency or early-stage hPGCs/hPGCLCs was unchanged between days 5 and 8, except for a slight decrease in LIN28A expression (Fig. S3E), and expression of late-stage hPGC marker genes, such as DAZL or DDX4, was not significant. These results suggest that hPGCLCs at day 8 EB culture did not reflect more advanced germline differentiation than at day 5. Interestingly, 152 DEGs that were expressed more strongly in day 8 hPGCLCs than day 5 hPGCLCs (Fig. S3 D, bar F and F) were significantly enriched with processed pseudogenes derived from ribosomal protein (RP) genes (eight RPS and three RPL pseudogenes; P < 0.0001 by permutation test) (Fig. S3F, asterisks). Since expres- sion of processed RP pseudogenes shows strong tissue specificity, presumably reflecting distinct epigenetic organizations (30), their differential expression in hPGCLCs might reflect epigenetic changes between days 5 and 8 EB cultures.To examine the reproducibility of the transcriptomal profiles across independent hiPSC clones, we generated two more hiPSC clones “A5” and “B6” from human neonatal foreskin fibroblasts using the same method by which clone A4 was generated. All three hiPSC clones had normal diploid karyotype (46, XY) (Fig. S4 A and C) and strongly expressed TRA-1–60, TRA-1–81, POU5F1, and NANOG.
These hiPSC clones were subjected to hPGCLC production [RNA-sequencing (RNA-seq) data nor- malization (Fig. S4 D and E), hierarchical clustering (Fig. S4 F, all DEGs and G, selected marker genes), and transcriptomal PCA (Fig. S4H)], and CD38+ hPGCLCs were isolated from day7 EBs. Unsupervised clustering (Fig. S4 F and G) and PCA (Fig.significantly expressed PRDM1 and TFAP2C, which are early- stage hPGC/hPGCLC markers (15), although expression of PRDM1 was weaker in CD38− cells than CD38+ cells (Fig. 2 A and B and Fig. S1 C1 and D). However, the mRNA transcriptsfor PRDM1 expressed in CD38− cells lacked most of the 3′-UTR, whereas CD38+ hPGCLCs expressed transcripts with the full-length 3′-UTR (24, 25) (Fig. 2D and Fig. S2B). Such 3′-UTR length polymorphism seemed to be specific to PRDM1 among genes shown in Fig. 2B (Fig. S2B). The 3′-UTR of PRDM1 mRNA missing in CD38− cells contained binding sites for miRNA let-7, miR-9, and miR-BHRF1-2 (Fig. 2D and Fig. S2B), whichplay significant roles in PRDM1 protein degradation in normal or malignant human B lymphocytes (26–28). West et al. (21) reported that in vitro differentiation of mouse ESCs to PGCLCs is dependent on suppression of let-7 miRNA by Lin28 and thatoverexpression of Prdm1 rescued differentiation of lin28-deficient ESCs to PGCLCs. In our hPGCLC production system, both LIN28A and LIN28B were significantly expressed in CD38+ and CD38− cells (Fig. S1D). The possible importance of the 3′-UTRcompared CD38+ hPGCLC and CD38− EB cell transcriptomes at days 5 (Fig. 2B), 7 (Fig. S4G), and 8 (Fig. S3E) and observed a decrease in PRDM1 and TFAP2C expression in CD38− cells over time, whereas the expression of these genes was maintained in CD38+ hPGCLCs.
These results suggest that sustained expres- sion of PRDM1 and TFAP2C in EB cells may be necessary for hiPSC differentiation to hPGCLCs.Transcriptomal Consistency of hPGCLCs Produced Using Various Protocols and FACS Markers. To obtain insight into biological consistency of hPGCLCs produced using different protocols (Table S1), we compared transcriptomes of our hPGCLCs with those reported by Irie et al. (15) and Sasaki et al. (16) as well as human fetal gonocytes (23), gestation week 7 hPGCs, gonadal somatic cells, and the TCam-2 human seminoma cell line (15). Hierarchical clustering involving 46 selected pluripotency andgermline marker genes clearly separated the subjects into six clusters corresponding to cell types—namely PSCs, CD38− EB cells, hPGCLCs, advanced germline cells (AGCs), gonadal somatic cells,and TCam-2 (Fig. 3A; color-coded protocols and cell types are shown in Fig. 3B). The PSC cluster cells strongly expressed pluri- potency markers (Fig. 3A, rectangle 1). Our CD38− EB cells strongly expressed ICM/germline markers LEF1, DPPA3, KLF2, KLF4, TFC2L1, TCL1B, and TFAP2C (Fig. 3A, rectangle 2),whereas expression of other germline markers, such as CD38, NANOS3, or SOX17, was low. PRDM1 was significantly expressed in CD38− EB cells, although less strongly than hPGCLCs. The hPGCLC cluster cells strongly expressed known hPGC/hPGCLC markers (Fig. 3A, rectangle 3). The especially strong expression of five markers—namely T, CD38, PRDM1, NANOS3, and SOX17—in all hPGCLCs produced by three independent laboratories suggests their practical usefulness as an hPGCLC signature, whereas ex- pression of other hPGC/hPGCLC markers, such as TCL1A, KIT,DND1, SOX15, or GATA4, was either heterogeneous or relatively low. Expression of late-stage germline markers, including DAZL or SYCP3, was generally weak in hPGCLCs but strong in week 7 em- bryonic hPGCs, AGCs, and TCam-2 cells (Fig. 3A, rectangle 4).
PCA with these 46 marker genes also separated the subjects into five clusters along cell types (Fig. 3C and Fig. S3G). These results show significant similarity of the CD38+ hPGCLCs produced in this study with those generated in previous reports.Hierarchical clustering with 131 genes most strongly expressed in our CD38+ hPGCLCs (Fig. S1C2) also roughly separated the subjects along with the cell types (Fig. S5). In this heat map, a set of 34 genes was strongly expressed in hPGCLCs (Fig. S5 A, blue bar and B, enlargement), and 9 of them—namely HOXD1, ALOX5, WNT2, VSTM2B, PLDB1, IRX6, KANK3, TMEM123,and TMPRSS11E (RP11-387A1.6 is identical to the second exon of HOXD1)—were specific to hPGCLC. Other genes, includinggermline markers NANOS3, SOX17, TCL1A, CD38, T, and WNT3, were also expressed in 7-wk hPGCs, week 16–16.5 gon- ocytes, or TCam-2. Several known hPGC markers, such as PRDM1, TFAP2C, or KIT, were not included in this heat map, because they were also significantly expressed in CD38− EB cells and/or pluripotent precursors. Strong expression of NANOS3, SOX17, TCL1A, T, and CD38 was also observed in hPGCLCs produced and enriched as KIT+ cells by von Meyenn et al. (9) (Fig. S6 and Table S1). Expression of UTF1, DND1, or SOX15 was weak in their hPGCLCs (Fig. S6B) as well as hPGCLCsproduced by other laboratories (Fig. 3A, rectangle 3). Significant induction of NANOS3 and SOX17 was also observed by Sugawa et al. (17) with their hPGCLCs enriched as KIT+/TRA-1–81+ cells (Table S1). Taken together, these results suggest that hPGCLCs produced via different pre-EB precursors and/orFACS enrichment marker antigens have comparable gene expression characteristics.Enrichment of TF Binding Motifs in Promoters of DEGs Between CD38+ hPGCLCs and CD38− EB Cells.
To obtain insight into the mecha- nisms of hiPSC differentiation to CD38+ hPGCLCs, we exam-ined enrichment of TF binding motifs in promoter sequences of DEGs between CD38+ and CD38− cells in day 5 EBs using the Transcription Factor (TRANSFAC) database and the F-match algorithm (31). We identified 13 TF binding motifs that were significantly enriched in promoters of 537 genes expressed more strongly in CD38+ hPGCLCs than CD38− EB cells (Fig. 4), whereas no enriched motifs were found with genes down-regulated in CD38+ hPGCLCs. Strikingly, 1 of these 13 enriched TF motifs was the optimal 9-nt TFAP2C binding sequence GCCTGAGGC(32). Although other members of the TFAP2 TFs can also bind to this sequence, expression of TFAP2C was strongest among them in CD38+ hPGCLCs and their precursors (Fig. S1E). The list of 44 genes with promoters that harbor predicted TFAP2 binding consensus motifs (Table S2) included SOX17, NANOS3, and T, which are robustly expressed in hPGCLC (Figs. 2B and 4B, lo- cations of TFAP2 binding motifs), suggesting the functional sig- nificance of TFAP2 for hPGCLC induction. However, our list showed only a marginal overlap with a previously reported list of predicted Tfap2c-regulated genes in mouse PGCs (33), sharing asfew as three genes—namely NANOS3, PRUNE2, and GSTM1. Among other TFs shown in Fig. 4, RREB-1 is also noteworthy,because it is required for formation of cell protrusions and lamellipodia as well as the dynamic cell–cell interactions (34).Localization of CD38+ hPGCLCs at Surface of EBs. Since iPSC dif- ferentiation to hPGCLCs was associated with enriched expres- sion of cell migration genes (Fig. 2E), we determined locations ofhPGCLCs in EBs by immunohistochemistry.
Because an anti-CD38 antibody suitable for identification of hPGCLCs in formaldehyde- fixed, paraffin-embedded (FFPE) slides was not available, we used POU5F1 as an alternative hPGCLC marker antigen (35). Although POU5F1 was expressed in hPSCs, its expression in day 7 EBs was specific to CD38+ hPGCLCs (Fig. 5A). Unexpectedly, we observed POU5F1+ hPGCLCs almost exclusively in the outermost surface monolayer of EB cells (Fig. 5B and Fig. S7B). The majority of the POU5F1+ cells (>90%) were distributed at the EB surface without aggregation (Fig. 5B, Left), although aggregated clusters of POU5F1+ cells were occasionally observed(Fig. 5B, Right). Significantly biased localization of POU5F1+ hPGCLCs in the outermost surface monolayer of EBs was sta- tistically confirmed (Fig. 5C). A weak POU5F1 signal, readily suppressed by low concentrations of a POU5F1 blocking peptide (Fig. S7A), was also detected predominantly at the subsurface area of EBs (Fig. S7B). We also attempted to determine quan- titative aspects of expression of PRDM1 protein by immuno- fluorescence and immunohistochemistry but were unsuccessful because of technical difficulties.Because the DEGs up-regulated in CD38+ hPGCLCs rather than in CD38− EB cells were enriched for cell migration genes and because the majority of hPGCLCs at the EB surface were observed without forming cell aggregates, we speculated thathPGCLCs may have highly migratory characteristics just as hu- man embryonic PGCs do (29). To obtain evidence of hPGCLC migration, we attempted to obtain time-lapse images of CD38+ cells by immunofluorescence staining of live cells. Whole EBs at day 6 were unsuitable for live fluorescence imaging because of strong autofluorescence. However, when EBs were inoculated on Matrigel-coated dishes to spread overnight, cells that emerged from the EB mass and formed a monolayer were largely free of autofluorescence. Taking advantage of this observation, we stained live cells surrounding the immobilized EB masses with a fluorescence dye-conjugated anti-CD38 antibody and obtained antibody-dependent fluorescence images matching cell contours and cellular protrusions (Fig. S7C). Two live fluorescence images taken at a 60-min time interval showed significant changes in CD38+ cell shape and location (in Fig. 5D, fluorescence colors were digitally adjusted to green and orange for images taken be- fore and after the interval, respectively), agreeing with the highly migratory nature of human embryonic PGCs (29). Our attempts to obtain video microscopic tracings of live-stained CD38+ cells were unsuccessful, apparently because of the apparent sensitivity of CD38+ cells to changes in the cell culture environment.Expression of the CXCR4 Chemotactic Cytokine Receptor in CD38+ PGCLCs and Transcriptional Effects of CXCL12.
In mouse embryos,PGCs migrate laterally from the midline regions toward genital ridges following a gradient of chemoattractant ligand CXCL12/ SDF1 secreted by gonadal somatic cell precursors and its plasma membrane receptor CXCR4 expressed by PGCs (36). Consistent with these observations, CXCR4 was strongly expressed in hPGCLCs and in the TCam-2 seminoma cells, while CXCL12 expression was strong in human embryonic gonadal somatic cells (Fig. 5E). Neither CD38+ hPGCLCs nor CD38− EB cellsexpressed significant amounts of CXCL12 mRNA. FACS con-tour plots (Fig. 5F) showed that the EB cell population stained only with either anti-CD38 or anti-CXCR4 antibody completely disappeared when cells were double-stained with these two antibodies, while a new CD38+/CXCR4+ double-positive cell population emerged. These results indicated that all CD38+ hPGCLCs expressed CXCR4 and that all CXCR4-positive cells in EBs were CD38+ hPGCLCs.We next attempted to examine the transcriptional effects of exogenous CXCL12 on gene expression in hPGCLCs (based on three replicates of RNA-seq). CXCL12-activated CXCR4 signaling in human CD34+ hematopoietic stem cells (37) or Jurkat T cells (38) requires Rho GTPase-activated kinases 1/2 (ROCK1/2). However, a potent ROCK1/2 inhibitor Y27632 is included in our hPGCLC medium as well as the media used by Irie et al. (15), Sasaki et al. (16), and Sugawa et al. (17) but not that used by von Meyenn et al. (9) (Table S1). When we removed Y27632 from hPGCLC medium starting at day 5 EB culture, CD38+ cells isolated at day 7 still strongly expressed NANOS3, NANOG, and SOX17 as well as CD38 (Fig. 6A), indicating that Y27632 is dispensable for hPGCLC production on and after day5. Under this Y27362-free condition, addition of CXCL12 to the hPGCLC medium on and after day 5 resulted in suppression ofCXCR4 expression (Fig. 6A), consistent with an observation made with human CD34+ hematopoietic stem cells (37) and possibly reflecting a negative feedback regulation. Expression of PTGER3 (encoding prostaglandin E2 receptor isoform 3) and BCL2 (encoding an antiapoptotic protein) was dramatically in- duced after exposure to CXCL12, whereas expression of CASP8 (encoding caspase 8, the initiator caspase required for death receptors-induced apoptosis) was suppressed (Fig. 6A).
These changes in expression of apoptosis-regulating genes suggest that CXCL12 may reduce CD38+ cell sensitivity to cytokine-induced apoptosis.To predict the transcriptomal impact of CXCL12 on CD38+hPGCLCs, DEGs were identified from the three replicates of RNA-seq data (1,190 up-regulated and 735 down-regulated in the presence of CXCL12 during days 5–8 of EB culture in the PGCLC medium lacking Y27632) and subjected to GO analysis using the DAVID server. Indeed, GO analysis of CXCL12-induced DEGs in CD38+ hPGCLCs showed significant enrich- ment of genes involved in negative regulation of apoptosis as well as cell migration (Fig. 6B). However, CXCL12-suppressed DEGs were enriched for genes involved in various aspects of mitotic nuclear division (Fig. 6C). These results suggest that CXCR4 activation in hPGCLCs may promote cell migration coupled with suppression of apoptosis and mitotic nuclear di- vision. Additional studies will be needed to examine in detail the biological functions of CXCR4 in hPGCLCs and their relevance to human embryonic PGCs in vivo.
Discussion
Our study provides evidence that hPGCLCs produced using different protocols (Table S1), including ours involving only transient exposure of primed pluripotency hiPSCs to the 4i me- dium (Fig. 1), share concordant transcriptomal characteristics, despite a variety of pre-EB precursors and FACS enrichment strategies (Fig. 3). Based on the transcriptomal profiling along with steps of hPGCLC production (Fig. 2 A and B), we speculate that specific induction of cluster 2 genes may be important for germline differentiation of EB cells, while persistent expression of cluster 3 genes may prevent differentiation to somatic line- ages. Irie et al. (15) reported that human ES cells (hESCs) maintained in the 4i medium for more than 2 wk showed a slight up-regulation of T and HAND1 compared with the conventional hESCs, speculating that 4i-reprogrammed hESCs might have primitive streak-like features, which could explain the compe- tence to directly give rise to hPGCLCs. Our data also show a slight but significant increase in expression of T in the short-term 4i reprogrammed hiPSCs compared with the primed-state hiPSCs (Fig. 2B). Future studies will be needed to evaluate the relevance of the short-term 4i reprogrammed cells with the primitive streak-like state. KJF4, TFAP2C, and cluster 1 genes, with importance in specification of hPGC/hPGCLCs that has been shown or predicted, were expressed in both CD38+ and CD38− EB cells. Expression of KLF4, TCL1B, and TFCP2L1, which are markers of ICM-like strong naïve pluripotency (16, 20), may imply a certain degree of commonality of gene networks between EB cells exposed to the hPGCLC medium and ICMs.
The strong expression of hPGC/hPGCLC markers, such as NANOS3 (15), indicates that our CD38+ cells isolated from day 5 EBs are hPGCLCs. However, weak but significant expression of NANOS3 and SOX17 in CD38− EB cells at day 5 (Fig. 2B) might reflect their undergoing commitment to germline differentiation. Kobayashi et al. (18) recently showed that high levels of SOX17 are needed for robust induction of hPGCLCs. While SOX17 is essential for hPGCLC induction (15), it is not a specific germline maker but also labels mesendodermal cells. Additional studies will be necessary to determine whether the predicted hPGCLC precursor cells weakly expressing SOX17 and/or NANOS3 in the CD38− cells in day 5 EBs can later differentiate to hPGCLCs.
We also observed that the PRDM1 mRNA transcripts expressed in CD38− EB cells lacked the 3′-UTR harboring miRNA binding sites required for regulation of mRNA stability (26–28) (Fig. 2D and Fig. S2). Whether the 3′-UTR polymorphism of PRDM1 mRNA and/or let-7 miRNA have significance in hiPSC differ-entiation to hPGCLCs remains to be determined. Specification of mouse PGCs requires coordinated actions of three TFs—namely Prdm1, Prdm14, and Tfap2C—which have been described as a tripartite TF network (5, 13, 14). In contrast, specification of hPGCLCs is dependent on SOX17, which induces
PRDM1 (15) and endoderm genes. PRDM1 plays critical roles in specification of hPGCs/hPGCLCs by suppressing SOX17-induced expression of endodermal genes as well as BMP- or WNT-induced expression of mesodermal genes (13, 18). Whereas TFAP2C mRNA was strongly expressed in both CD38+ hPGCLCs and
CD38− EB cells (Figs. 2 A and B and 3 and Figs. S1, S3, and S4), promoter sequences of DEGs up-regulated in hPGCLCs were significantly enriched with the TFAP2 binding motif (Fig. 4A).
Interestingly, these possible TFAP2-regulated DEGs up-regulated in hPGCLCs (Table S2) included SOX17, NANOS3, and T (Fig. 4B). These results support the importance of TFAP2C in hiPSC differentiation to hPGCLCs, although its expression is not re- stricted to hPGCLCs.GO analysis of DEGs up-regulated in CD38+ hPGCLCs compared with CD38− EB cells revealed significant enrichment of genes involved in cellular activities relevant to migration, including cell projection morphogenesis, or regulation of loco- motion (Fig. 2E). Promoter sequences of these DEGs were enriched with the binding motif of RREB-1 (Fig. 4A), a TF required for migration and spreading of human MCF- 10A cells (34). Defects of its Drosophila homolog gene hindsignt caused excess accumulation of intercellular adhesion molecules and suppression of cell spreading in ovarian follicular epithe- lium, indicating its importance in reducing cell–cell adhesion when interconnected groups of cells undergo dynamic changes in cell shape (34). Based on these roles of hindgut/RREB-1, we speculate that human RREB-1 may be important during mi- gration of PGCs, involving transient and reversible adhesion of PGCs to other types of cells as well as to marked changes in cell shape (Fig. 5D).We detected hPGCLCs in EBs almost exclusively at their outermost surface monolayer (Fig. 5 A and C and Fig. S7 A and B). Live cell imaging revealed active migration of hPGCLCs (Fig. 5D). In mouse EBs it has been shown that only the out- ermost surface monolayer will differentiate to primitive endo- derm (PrE) and noteworthy visceral endoderm (VE) (39, 40). This is because in vivo specification of mouse PGCs from epi- blast requires BMP2 secreted from VE to activate the tripartite TF network (5). Before implantation, mouse PrE deposits basement membrane (BM) to its boundary with ICM, and this BM is apparently required for PrE differentiation to VE (40). Similarly, several BM proteins are deposited beneath the out- ermost monolayer of mouse EBs (41, 42).
Mouse PSCs partially differentiated toward PrE emerge in the standard naïve pluri- potency culture, and on EB formation, such PrE-primed PSCs move toward the outside because of differential adhesion (39, 40). Exclusive localization of hPGCLCs at the outermost monolayer of EBs might occur for similar reasons. The distri- bution of hPGCLCs in wide areas of EB surface without ag- gregation may reflect their random migration on the subsurface BM. Compared with the conventional methods for forming loose cell aggregates using low-attachment U- or V-shaped wells (Table S1), our use of microwells for the spin method EB for- mation (19) may increase initial density of hiPSCs, resulting in tighter intercellular adhesion, which in turn, affects directed differentiation of EBs (40).Our GO analysis of DEGs between CD38+ hPGCLCs and CD38− EB cells in day 5 EBs revealed significant enrichment of cell migration genes (Fig. 2E). CXCR4 was strongly expressed in all hPGCLCs (Figs. 5 E and F and 6A), whereas expression of
CXCL12 in whole EB cells was weak (Fig. 5E). Expression of c-KIT, another plasma membrane receptor involved in PGC mi- gration and survival, was also strong in CD38+ hPGCLCs (Fig. 2 A and B and Figs. S3 D and E and S4 F and G). Importantly, our hPGCLC production medium contained stem cell factor (SCF), the c-KIT ligand, which is also the case for most other hPGCLC production protocols (Table S1). Removal of SCF from the hPGCLC production medium resulted in significant dismantling of EBs with no yield of hPGCLCs.
In mouse embryos, CXCL12 is produced at the genital ridges and diffused into the surrounding mesenchyme, guiding chemotactic migration of CXCR4+ PGCs (36, 43) (Fig. S8). In contrast, SCF is expressed by somatic cells throughout the path of PGC migration to support survival and the random migratory activities of PGCs (43–45) (Fig. S8). Supported by SCF in the medium, hPGCLCs may survive and randomly migrate on the surface of EBs. We, therefore, speculate that hPGCLCs at the surface of EBs may resemble the prechemotactic embryonic hPGCs in the midline regions (Fig. S8). The reported low degree of gDNA demethylation in hPGCLCs agrees with the notion that hPGCLCs resemble midline hPGCs (9, 15–17, 20).Exposure of hPGCLCs at the surface of EBs to CXCL12 in the absence of ROCK1/2 inhibitor induced cell migration genes and antiapoptotic genes (Fig. 6 A and B), while expression of CASP8 was suppressed (Fig. 6A). These results may imply that the midline hPGCs may be prone to the extrinsic pathway apo- ptosis until the exposure to CXCL12 (Fig. S8). Exposure of hPGCLCs to CSCL12 also induced PTGER3 (Fig. 6A). Notably, Ptger3 activation has been reported to be required for expression of both Cxcr4 and Cxcl12 by dendritic cells in mouse lymph nodes (46), and a nine-gene biomarker predicting metastatic prognosis of human head and neck squamous cell carcinomas with 90% accuracy included CXCR4, CXCL12, and PTGER3 (47). It is thus interesting to speculate that there may be func- tional interactions between prostaglandin E2 signaling and CXCL12 signaling in hPGCs/hPGCLCs.
In summary, we have shown robust production of hPGCLCs after a 72-h culture of primed pluripotency hiPSCs in the 4i medium as well as transcriptomal consistency of hPGCLCs using different protocols and/or FACS enrichment antigens. Most hPGCLCs strongly expressed CXCR4 and were localized at the outermost surface monolayer of EBs. Exposure of hPGCLCs to CXCL12 induced migratory and antiapoptotic genes. Thus, our results provide evidence that hPGCLCs produced using different protocols have comparable biological characteristics resembling prechemotactic, embryonic midline hPGCs.Footprint-free hiPSC clones were generated from commercially available, deidentified neonatal foreskin dermal fibroblasts and maintained in mTeSR1 on Matrigel. hPGCLCs were generated from hiPSCs using a protocol based on the study of Irie et al. (15), with modifications involving 72-h ex- posure of primed pluripotency hiPSCs to the 4i medium and the use of AggreWell microwells for the spin EB technique of EB formation (19). CD38+ and CD38− cells were FACS-enriched from single-cell suspension using FACS.Total RNA was subjected to TaqMan qPCR or RNA-seq using Illumina Next-Seq 500, and deep sequencing data were analyzed using the STAR aligner and R/Bioconductor packages Rsubread, edgeR, prcomp, and plot3D. Pro- moter enrichment analysis of TF binding motifs was performed using the TRANSFAC database. Immunohistochemistry of FFPE serial sections was performed using standard methods. Time-lapse immunofluorescence live imaging was performed on EBs spread on Matrigel using CD38 inhibitor 1 fluorescence dye- conjugated anti-CD38 antibodies and an inverted fluorescence microscope.