Skip to main content

Analysis of multi-lineage gene expression dynamics during primordial germ cell induction from human induced pluripotent stem cells

Abstract

Background

In mammals, specification of primordial germ cells (PGCs) is established in the early post­implantation embryo. The bone morphogenetic protein (BMP)-SMAD and WNT3-β-catenin signaling initiate the gene regulatory network for PGC specification. The activation of SOX17-BLIMP1 axis is critical for human PGC program. Moreover, EpCAM and INTEGRINα6 were identified as surface markers of human PGC-like cells (PGCLCs) recently. However, the signaling mechanism for PGC specification in non­rodent mammals remains to be clarified.

Methods

We differentiated human induced pluripotent stem cells (hiPSCs) into PGCLCs in vitro in response to Activin A and BMP4. The percentage of EpCAM/INTEGRINα6 double-positive cells (PGCLCs) was analyzed by flow cytometry. The expression of PGC genes was evaluated by qRT-PCR and immunofluorescence. The expression dynamic of multi-lineage genes during the differentiation process was evaluated by qRT-PCR.

Results

Under the stimulation for PGCLC induction, the embryoids derived from hiPSCs initiated significant upregulation of the early PGC genes (BLIMP1, TFAP2C, and NANOS3), but maintained low or no levels of DPPA3 and late PGC genes (DAZL and DDX4). The percentage of EpCAM/INTEGRINα6 double-positive PGCLCs reached the highest at day 6 of induction. After pre-induction, the incipient mesoderm-like cells (iMeLCs) upregulated most of the mesoderm genes (EOMES, T, MSXI, RUNX2, and MIXL1). The differentiating embryoids showed high levels of key pluripotency genes, OCT4 and NANOG, but became negative for SOX2. In contrast to iMeLCs, the differentiating embryoids downregulated mesoderm genes RUNX2 and EOMES, and ectoderm gene PAX6, but increased the expression of endoderm gene SOX17.

Conclusions

During PGCLC induction process in vitro, the differentiating embryoids not only activated the PGC-related genes, but also displayed complex regulation of pluripotency genes and multi-lineage genes. These results would be meaningful for future research investigating the regulation of human early germ line development.

Introduction

The germ cell lineage is the source of totipotency and transmits genetic and epigenetic information across generations. In mammals, primordial germ cells (PGCs) are the founder cells of the germ cell lineage and emerge at early stage of post-implantation development. The PGCs undergo a number of complex developmental events, including repression of somatic programs, (re)acquisition of potential pluripotency, and epigenetic reprogramming, and finally initiate oogenesis or spermatogenesis to form gametes [1]. Currently, the mechanism of PGC specification has been extensively investigated in mice, which provides insight into mammalian development. In mice, the bone morphogenetic protein (BMP) and WNT signals from extra-embryonic tissues induce the expression of several PGC fate regulators in a few germ line competent cells [2,3,4]. Particularly, a tripartite transcription factor network of PR domain zinc finger protein 1 (PRDM1, also known as BLIMP1), PRDM14, and transcription factor AP-2 gamma (TFAP2C) represses the somatic fate and promotes the mouse PGC specification [5].

By contrast, the knowledge about human PGC specification is relatively limited because of the technical and ethical obstacles to obtain germ line cells from early human embryos. Notably, in the last decade, generation of germ cells from pluripotent stem cells (embryonic stem cells, ESCs, and induced pluripotent stem cells, iPSCs) has provided a surrogate model for germ line development in vitro. In mice, ESCs and iPSCs with naive pluripotency could be induced into primordial germ cell-like cells (PGCLCs) through epiblast-like cells (EpiLCs) [6]. However, human ESCs (hESCs) and iPSCs (hiPSCs) are considered to show a primed pluripotency with limited potential for germ cell fate and respond poorly to the method used for mouse PGCLC induction [7,8,9]. Accordingly, Gafni et al. established defined conditions to generate human naive pluripotent stem cells from primed hESCs and hiPSCs [10]. Moreover, Irie et al. developed an approach for human PGCLC (hPGCLCs) specification from the ground state hESCs and hiPSCs [7]. Strikingly, hPGCLCs are robustly induced in vitro from hiPSCs in a primed pluripotent state through incipient mesoderm-like cells (iMeLCs) [11]. Based on these differentiation models, several key regulators of human PGC (hPGC) fate as well as the regulation network they formed were clarified, including EOMES and SOX17 [7, 12].

The mammalian germline is set aside from somatic lineages in early post-implantation embryos [1]. During the in vitro hPGCLC specification process, only some of the pluripotent stem cells respond to the induction signals, and the other cells still spontaneously differentiated into somatic lineages. Here, we set out to differentiate hiPSCs into hPGCLCs in vitro with the iMeLC method and determine the expression dynamics of multi-lineage genes, in order to uncover new clues for cell fate decision.

Materials and methods

Culture of hiPSCs

Human fibroblasts were isolated from foreskin of three volunteers at Wuhan Tongji Reproductive Medicine Hospital with written informed consent. The fibroblasts were reprogrammed with the Yamanaka KOSM (KLF4, OCT4, SOX2, and c-MYC) transcriptional factors using the lentivirus vectors. The derived hiPSCs were maintained in mTeSR1 medium (Stem Cell Technologies) on Matrigel (Corning)-coated dishes. The medium was changed every day. Cells were passaged every 3 to 5 days using Accutase (Life Technologies). For single-cell dissociation, the cells were treated with 1 to 1 mixture of TrypLE Select (Life Technologies) and 0.5 mM EDTA/PBS. Ten micrometers of a ROCK inhibitor (Y-27632, TOCRIS bioscience) was added for 24 h after passaging.

Induction of hPGCLCs

For pre-induction, hiPSCs were dissociated with 0.5 mM EDTA/PBS, and 3× 105 cells per well were plated on Matrigel-coated 12-well plates in GK15 medium (G-MEM [Thermo Fisher] supplemented with 15% KSR [Thermo Fisher], 0.1 mM NEAA [Thermo Fisher], 2 mM L-glutamine [Thermo Fisher], 1 mM sodium pyruvate [Thermo Fisher], 0.1 mM 2-mercaptoethanol [Sigma]) containing 3 μM CHIR (Selleck Chemicals), 50 ng/ml Activin A (PEPRO TECH), and 10 μM ROCK inhibitor (Y-27632, TOCRIS bioscience). After 2 days of pre-induction, the cells were dissociated with Accutase (Thermo Fisher) and plated into ultra-low cell attachment U-bottom 96-well plates (Corning) at a density of 2000–4000 cells per well to form embryoid bodies in 200 μl of GK15 medium containing 200 ng/ml BMP4 (R&D Systems), 20 ng/ml human LIF (R&D Systems), 100 ng/ml SCF (R&D Systems), 50 ng/ml EGF (R&D Systems), and 10 μM ROCK inhibitor (Y-27632, TOCRIS bioscience). H1 hESC was used as a control for PGC induction.

Flow cytometry

The floating embryoid bodies were dissociated with 0.05%Trypsin-EDTA/PBS for 15 min at 37 °C. After washing with PBS, the cell suspension was filtered by cell strainer to remove cell clumps and then subjected to centrifugation. Then, the dissociated cells were stained with PE-conjugated anti-human EpCAM (eBioscience) and FITC-conjugated anti-human INTEGRINα6 (eBioscience). The stained cells were resuspended in PBS and analyzed with a flow cytometer (Beckman, DxFLEX).

RNA extraction and quantitative RT-PCR

Total RNA was extracted using DirectZol RNA mini-prep (Zymo research) according to manufacturer’s instructions. Reverse transcription was performed using the RevertAid First Strand cDNA synthesis kit (Thermo Fisher Scientific). The quantitative RT-PCR was performed using SYBR Premix Ex Taq II (Takara). All gene expression analyses were performed with samples from three independent differentiation experiments. Values normalized to GAPDH are shown. Primers are listed in Table 1.

Table 1 Primers used in this study

Immunofluorescence

Cells were fixed in 4% paraformaldehyde/PBS for 15 min at room temperature and permeabilized with 1% Triton X-100/PBS (Sigma) for 10 min. After blocking with 10% donkey serum in PBS (Jackson ImmunoResearch) for 1 h at room temperature, cells were incubated with primary antibodies overnight at 4 °C, followed by incubation with appropriate fluorescent-conjugated secondary antibodies for 1 h at room temperature the next day. Primary antibodies were as follows: OCT4 (Abcam), SOX2 (Abcam), SSEA4 (Abcam), TFAP2C (Santa Cruz), PRDM14 (Abcam), and NANOS3 (Abcam). Secondary antibodies were as follows: AlexaFluor 488 conjugated donkey anti-rabbit IgG and AlexaFluor 594 conjugated donkey anti-mouse IgG (all Life Technologies). The nuclei were counterstained with DAPI (Thermo Fisher Scientific). The cells were observed with a Zeiss inverted confocal microscope.

Statistical analyses

All of the data were presented as mean ± standard deviation (SD). Statistical analyses were performed using the Student’s t test or one-way analysis of variance (ANOVA) with SPSS 17.0 software. A p value less than 0.05 was considered to be statistically significant.

Results

Differentiation of hiPSCs into PGCLCs in vitro

We established hiPSC lines from dermal fibroblasts of three male volunteers using a method published before [13]. The hiPSCs displayed typical hESC morphology and were positive for pluripotency markers, including OCT4, SOX2, and SSEA (Fig. 1a, b). After pre-induction, the hiPSCs were differentiated into flat iMeLCs with distinct cell borders, which were also positive for pluripotency markers (Fig. 1a, c). For PGCLC induction, the differentiating cells were maintained under a floating culture condition and aggregated to form embryoids (Fig. 1a). We analyzed the expression of PGC genes, including BLIMP1, TFAP2C, NANOS3, DPPA3, DDX4, and DAZL, during the differentiation process. The hiPSCs and hESCs showed no or low expression of the PGC genes, which remained at low levels after pre-induction, except for BLIMP1and TFAP2C, exhibiting slightly increased levels at iMeLC stage (Fig. 2a and Supplemental fig. 1 A). Under the stimulation for PGCLC induction, the embryoids differentiated from hiPSCs initiated significant upregulation of the early PGC genes (BLIMP1, TFAP2C, and NANOS3, p < 0.05), but maintained low or no levels of DPPA3 and late PGC genes (DAZL and DDX4) (Fig. 2a). However, we observed that the embryoids from hESCs increased the expression of NANOS3 and DPPA3 from day 4 of PGC induction (Supplemental fig. 1 A). Immunofluorescence analysis demonstrated that the expression of TFAP2C coincided with PRDM14, OCT4, and NANOS3 in day 4 embryoids (Fig. 2b). These findings indicate that the hiPSCs are capable of early PGC fate following PGCLC induction in vitro.

Fig. 1
figure 1

Characterization of hiPSC, iMeLCs, and embryoids. a Typical morphology of hiPSCs, iMeLCs, and embryoids. Scale bar, 200 μm. bd Immunofluorescence analysis for pluripotency markers of hiPSCs (b) and iMeLCs (c). The uterine tumor cells, AN3CA, were used as negative control (d). Scale bar, 100 μm

Fig. 2
figure 2

Expression analysis of germ cell specific genes during hPGCLC specification process. a Expression dynamics of germ cell specific genes during hPGCLC specification process, including iPSCs, iMeLCs, and the whole floating embryoids at day 2, 4, 6, and 8 of induction, as measured by qRT-PCR. Relative expression levels are shown with normalization to housekeeping gene GAPDH. Error bars indicate mean ± SD of three independent experiments. Red squares indicate values for embryoids exposed to cytokine stimulation; black triangles indicate values for embryoids formed spontaneously without cytokine stimulation. *p < 0.05 vs. the iPSC groups. n.d., not detected. b Immunofluorescence analysis for the expression of TFAP2C, PRDM14, OCT4, and NANOS3 in day 4 embryoids derived from hiPSCs. Scale bars, 100 μm

Efficiency of PGCLC induction in vitro

Recently, EpCAM and INTEGRINα6 were identified as surface markers of hPGCLCs in embryoids [11]. We performed fluorescent-activated cell sorting (FACS) analysis for the differentiated cells by these two markers. The results showed that the percentage of EpCAM/INTEGRINα6 double-positive cells in iMeLCs (day 0) was up to 46~49% without obvious segregation, but the germ cell fate was not activated at this point and these cells were not PGCLCs. As early as day 2 of PGCLC induction, the differentiating cells in embryoids began to segregate into two different populations. The EpCAM/INTEGRINα6 double-positive proportion was increased progressively until day 6, resulting in around 39~44% of double-positive putative hPGCLCs, and the two subpopulations became more distinct and persisted until day 8 of induction, albeit with reduced ratios at day 8 (Fig. 3 and Supplemental fig. 1 B).

Fig. 3
figure 3

Efficiency of PGCLC induction from iPSCs in vitro. a FACS analysis by EpCAM and INTEGRINα6 expression of cells during hPGCLC induction (until day 8) from hiPSCs by BMP4, LIF, SCF, and EGF (left) or by no cytokines (right). Boxed areas indicate EpCAM/INTEGRINα6 double-positive cells with their percentages. b Percentage of cells double positive for EpCAM/INTEGRINα6, only positive for EpCAM and only positive for INTEGRINα6 in day 2, 4, 6 and 8 floating embryoids determined by FACS, respectively. Error bars indicate mean ± SD of three independent experiments. *p < 0.001 vs. the day 2 groups

Multi-lineage gene expression dynamics during hPGCLC induction process

We also detected the expression dynamics of pluripotency and embryonic lineage genes during the PGCLC induction process by qPCR (Fig. 4). Overall, the embryoids derived from all lines upregulated or remained similar levels of key pluripotency genes, OCT4 and NANOG, but became negative for SOX2. Compared to undifferentiated cells, the differentiating cells exhibited high levels of genes associated with naive pluripotency such as ZFP42, KLF4, and TFCP2L1, whereas the expression of PRDM14, which plays an important role in the mouse PGC specification [3], did not change evidently after induction. For the embryonic lineage genes, we observed that hiPSCs showed low or no expression of genes involved in endoderm (GATA6, SOX17, and FOXA2), mesoderm (EOMES, T, MSXI, RUNX2, and MIXL1), and ectoderm (PAX6) development. After pre-induction, the differentiating cells upregulated most of the mesoderm genes (EOMES, T, MSXI, RUNX2, and MIXL1) and one ectoderm gene (PAX6) to some extent, but still showed a lack of endoderm gene expression. These results were consistent with the designated identity of iMeLCs by Sasaki et al. [11]. In contrast to iMeLCs, the differentiating embryoids downregulated mesoderm genes RUNX2 and EOMES and ectoderm gene PAX6. In mice, WNT/BMP signaling stimulates the upregulation of T expression, a key mesodermal factor for the onset of mPGC specification [14]. Here, the expression of T exhibited modest upregulation in PGCLCs, although it may not play a major role in hPGCLC specification [12]. Recently, it has been reported that WNT signaling activates EOMES at iMeLC stage to induce the expression of SOX17, a critical driver of hPGC fate, at PGCLC stage [7, 12]. Particularly, we noted that the hPGCLCs decreased the expression of EOMES and increased the expression of SOX17, compared with iMeLCs. Additionally, immunofluorescence analyses of the day 4 embryoids validated the expression of SOX17 and the repression of SOX2 in PGCLCs (Fig. 5 and Supplemental fig. 1 C), which were identical to the observation in human embryonic PGCs [15, 16].

Fig. 4
figure 4

Multi-lineage gene expression analysis during hPGCLC specification process. Relative expression levels of genes in hiPSCs, iMeLCs, and the whole floating embryoids were measured by qRT-PCR and shown with normalization to housekeeping gene GAPDH. Error bars indicate mean ± SD of three independent experiments. p < 0.05 vs. the iPSC-1 group, *p < 0.05 vs. the iPSC-2 group, $p < 0.05 vs. the iPSC-3 group, #p < 0.05 vs. the H1 ESC group. n.d., not detected

Fig. 5
figure 5

Expression analysis of SOX17 and SOX2 in PGC specification. a Immunofluorescence for the co-expression of SOX17 with SOX2 in day 4 embryoids derived from hiPSCs. Scale bars, 100 μm. b Quantification of the SOX17 positive cells in day 4 embryoids. Error bars indicate mean ± SD of five independent experiments

Discussion

Here, we differentiated hiPSCs into hPGCLCs in vitro and analyzed the expression dynamics of multi-lineage genes during the differentiation process. While the PGC fate regulator genes were activated, the embryonic lineage genes in the embryoids present complex expression patterns after induction, providing insights for future studies on human early embryo development.

From fertilization to blastocyst formation, mouse and human pre-implantation embryo development are morphologically similar; however, timing and molecular details of important developmental events occurring at this stage may differ, like zygotic genome activation (ZGA), X-inactivation, signaling requirement for pluripotency, and blastocyst lineage specification [17,18,19,20]. From implantation to gastrulation, human and mouse embryos become structurally distinct, and the germ layer formation takes place at this stage [18]. In mice, PGC specification occurs at the posterior epiblast before gastrulation [2], which is supported by the complicated signaling interactions between embryonic and extra-embryonic tissues in the egg cylinder [21]. However, less is known about the origin and specification mechanism of hPGCs in the early post-implantation human embryos. It is demonstrated that BMP signaling is very likely to be conserved for PGC specification in mammals, including human [7, 11, 22].

Until recently, several studies reconstituted hPGC specification from hESCs in vitro and defined the signaling and transcriptional programs for human germ cell specification in vitro. In brief, WNT signaling (ACTIVIN A) activates EOMES to induce the expression of SOX17, which is a critical regulator of hPGC fate and is upstream of BLIMP1 [7, 11, 12]. SOX17 upregulates BLIMP1 and, potentially, endoderm genes. TFAP2C is initially activated independently from SOX17 through the BMP signaling. SOX17 and TFAP2C work together to establish the hPGCLC transcriptional program in competent cells upon induction by BMP signaling. The expression of downstream BLIMP1 represses SOX17-induced endoderm genes, BMP- and WNT-induced mesoderm genes, as well as other somatic genes [7]. Hypothetically, this signaling and transcription model for hPGC specification would occur in the nascent amnion in early post-implantation embryos.

Additionally, PRDM14 is expressed in the differentiating cells at a relatively low level and the expression of T is increased after BMP induction, whereas the role of PRDM14 and T in hPGC specification remains to be clarified. Notably, SOX2 is suppressed as soon as the induction starts. BMP signaling and BLIMP1 expression may contribute to the rapid downregulation of SOX2 during hPGC specification [23, 24]. Moreover, the embryoids derived from hiPSCs show different regulation for the three germ layer genes, suggesting that only parts of the hiPSCs respond to the BMP signaling to initiate the germ line specification; the other cells may still get into somatic cell fate.

Conclusions

In conclusion, we demonstrated that during PGCLC induction process in vitro, the differentiating embryoids not only activated the PGC-related genes, but also displayed complex regulation of pluripotency genes and multi-lineage genes. Nevertheless, further studies on human and non-human primate embryo development, especially for high-throughput studies, are needed to explore the complicated regulation networks of hPGC specification.

Availability of data and materials

All relevant data are available from the authors upon reasonable request.

Abbreviations

PGCs:

Primordial germ cells

BMP:

Bone morphogenetic protein

iPSCs:

Induced pluripotent stem cells

PGCLCs:

PGC-like cells

PRDM1:

PR domain zinc finger protein 1

TFAP2C:

Transcription factor AP-2 gamma

ESCs:

Embryonic stem cells

EpiLCs:

Epiblast-like cells

iMeLCs:

Incipient mesoderm-like cells

ZGA:

Zygotic genome activation

FACS:

Fluorescent-activated cell sorting

References

  1. Kurimoto K, Saitou M. Mechanism and reconstitution in vitro of germ cell development in mammals. Cold Spring Harb Symp Quant Biol. 2015;80:147–54.

    Article  Google Scholar 

  2. Ohinata Y, Payer B, O’Carroll D, Ancelin K, Ono Y, Sano M, Barton SC, Obukhanych T, Nussenzweig M, Tarakhovsky A, et al. Blimp1 is a critical determinant of the germ cell lineage in mice. Nature. 2005;436(7048):207–13.

    Article  CAS  Google Scholar 

  3. Yamaji M, Seki Y, Kurimoto K, Yabuta Y, Yuasa M, Shigeta M, Yamanaka K, Ohinata Y, Saitou M. Critical function of Prdm14 for the establishment of the germ cell lineage in mice. Nat Genet. 2008;40(8):1016.

    Article  CAS  Google Scholar 

  4. Weber S, Eckert D, Nettersheim D, Gillis AJ, Schafer S, Kuckenberg P, Ehlermann J, Werling U, Biermann K, Looijenga LH, et al. Critical function of AP-2 gamma/TCFAP2C in mouse embryonic germ cell maintenance. Biol Reprod. 2010;82(1):214–23.

    Article  CAS  Google Scholar 

  5. Magnusdottir E, Dietmann S, Murakami K, Gunesdogan U, Tang F, Bao S, Diamanti E, Lao K, Gottgens B, Azim SM. A tripartite transcription factor network regulates primordial germ cell specification in mice. Nat Cell Biol. 2013;15(8):905–15.

    Article  CAS  Google Scholar 

  6. Hayashi K, Ohta H, Kurimoto K, Aramaki S, Saitou M. Reconstitution of the mouse germ cell specification pathway in culture by pluripotent stem cells. Cell. 2011;146(4):519–32.

    Article  CAS  Google Scholar 

  7. Irie N, Weinberger L, Tang WW, Kobayashi T, Viukov S, Manor YS, Dietmann S, Hanna JH, Surani MA. SOX17 is a critical specifier of human primordial germ cell fate. Cell. 2015;160(1–2):253–68.

    Article  CAS  Google Scholar 

  8. Nichols J, Smith A. Naive and primed pluripotent states. Cell Stem Cell. 2009;4(6):487–92.

    Article  CAS  Google Scholar 

  9. Theunissen TW, Powell BE, Wang H, Mitalipova M, Faddah DA, Reddy J, Fan ZP, Maetzel D, Ganz K, Shi L. Systematic identification of culture conditions for induction and maintenance of naive human pluripotency. Cell Stem Cell. 2014;15(4):471–87.

    Article  CAS  Google Scholar 

  10. Gafni O, Weinberger L, Mansour AA, Manor YS, Chomsky E, Ben-Yosef D, Kalma Y, Viukov S, Maza I, Zviran A. Derivation of novel human ground state naive pluripotent stem cells. Nature. 2013;504(7479):282.

    Article  CAS  Google Scholar 

  11. Sasaki K, Yokobayashi S, Nakamura T, Okamoto I, Yabuta Y, Kurimoto K, Ohta H, Moritoki Y, Iwatani C, Tsuchiya H, et al. Robust in vitro induction of human germ cell fate from pluripotent stem cells. Cell Stem Cell. 2015;17(2):178–94.

    Article  CAS  Google Scholar 

  12. Kojima Y, Sasaki K, Yokobayashi S, Sakai Y, Nakamura T, Yabuta Y, Nakaki F, Nagaoka S, Woltjen K, Hotta A, et al. Evolutionarily distinctive transcriptional and signaling programs drive human germ cell lineage specification from pluripotent stem cells. Cell Stem Cell. 2017;21(4):517–32 e515.

    Article  CAS  Google Scholar 

  13. Fang F, Li Z, Zhao Q, Ye Z, Zhao K, Gu X, Shen S, Reijo Pera RA, Xiong C. Derivation of human induced pluripotent stem cell line (HUSTi001-A) from a 40-year-old patient with idiopathic infertility. Stem Cell Res. 2018;27:136–9.

    Article  CAS  Google Scholar 

  14. Aramaki S, Hayashi K, Kurimoto K, Ohta H, Yabuta Y, Iwanari H, Mochizuki Y, Hamakubo T, Kato Y, Shirahige K. A mesodermal factor, T, specifies mouse germ cell fate by directly activating germline determinants. Dev Cell. 2013;27(5):516–29.

    Article  CAS  Google Scholar 

  15. Perrett RM, Turnpenny L, Eckert JJ, O'shea M, Sonne SB, Cameron IT, Wilson DI, Meyts ER-D, Hanley NA. The early human germ cell lineage does not express SOX2 during in vivo development or upon in vitro culture. Biol Reprod. 2008;78(5):852–8.

    Article  CAS  Google Scholar 

  16. De Jong J, Stoop H, Gillis A, van Gurp R, van de Geijn GJ, de Boer M, Hersmus R, Saunders P, Anderson R, Oosterhuis J. Differential expression of SOX17 and SOX2 in germ cells and stem cells has biological and clinical implications. J Pathol. 2008;215(1):21–30.

    Article  Google Scholar 

  17. Rossant J, Tam PPL. New insights into early human development: lessons for stem cell derivation and differentiation. Cell Stem Cell. 2017;20(1):18–28.

    Article  CAS  Google Scholar 

  18. Tang WW, Kobayashi T, Irie N, Dietmann S, Surani MA. Specification and epigenetic programming of the human germ line. Nat Rev Genet. 2016;17(10):585.

    Article  CAS  Google Scholar 

  19. Blakeley P, Fogarty NM, del Valle I, Wamaitha SE, Hu TX, Elder K, Snell P, Christie L, Robson P, Niakan KK. Defining the three cell lineages of the human blastocyst by single-cell RNA-seq. Development. 2015;142(18):3151–65.

    Article  CAS  Google Scholar 

  20. Niakan KK, Eggan K. Analysis of human embryos from zygote to blastocyst reveals distinct gene expression patterns relative to the mouse. Dev Biol. 2013;375(1):54–64.

    Article  CAS  Google Scholar 

  21. Ohinata Y, Ohta H, Shigeta M, Yamanaka K, Wakayama T, Saitou M. A signaling principle for the specification of the germ cell lineage in mice. Cell. 2009;137(3):571–84.

    Article  CAS  Google Scholar 

  22. Kee K, Gonsalves JM, Clark AT, Pera RA. Bone morphogenetic proteins induce germ cell differentiation from human embryonic stem cells. Stem Cells Dev. 2006;15(6):831–7.

    Article  CAS  Google Scholar 

  23. Rao J, Pfeiffer MJ, Frank S, Adachi K, Piccini I, Quaranta R, Arauzo-Bravo M, Schwarz J, Schade D, Leidel S, et al. Stepwise clearance of repressive roadblocks drives cardiac induction in human ESCs. Cell Stem Cell. 2016;18(3):341–53.

    Article  CAS  Google Scholar 

  24. Lin IY, Chiu FL, Yeang CH, Chen HF, Chuang CY, Yang SY, Hou PS, Sintupisut N, Ho HN, Kuo HC, et al. Suppression of the SOX2 neural effector gene by PRDM1 promotes human germ cell fate in embryonic stem cells. Stem Cell Reports. 2014;2(2):189–204.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (81901484) and the Health commission of Hubei Province scientific research project (WJ2019H263).

Author information

Authors and Affiliations

Authors

Contributions

The study was conceived and designed by FF, ZL and KN. FF and ZL conducted the majority of experiments. FF and KN analyzed the data. QZ contributed to the cell culture and the immunofluorescence staining. FF and KN wrote the manuscript with the help from CX. The study was supervised by CX and KN. The final version was approved by all the authors.

Corresponding author

Correspondence to Ke Ni.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethics Committee of Tongji Medical College, Huazhong University of Science and Technology (S096).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S1.

Differentiation of hESCs into PGCLCs in vitro. (A) Expression dynamics of germ cell specific genes during hPGCLC specification process, including hESCs, iMeLCs, and the whole floating embryoids at day 2, 4, 6 and 8 of induction, as measured by qRT-PCR. Relative expression levels are shown with normalization to housekeeping gene GAPDH. Error bars indicate mean ± SD of three independent experiments. Red squares indicate values for embryoids exposed to cytokine stimulation; black triangles indicate values for embryoids formed spontaneously without cytokine stimulation. *p < 0.05 vs. the hESC groups. n.d., not detected. (B) FACS analysis by EpCAM and INTEGRINα6 expression of cells during hPGCLC induction (until day 8) from hESCs by BMP4, LIF, SCF, and EGF (left) or by no cytokines (right). Boxed areas indicate EpCAM /INTEGRINα6 double positive cells with their percentages. (C) Immunofluorescence for the co-expression of SOX17 with SOX2 in day 4 embryoids derived from hiPSCs. Scale bars, 100 μm.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fang, F., Li, Z., Zhao, Q. et al. Analysis of multi-lineage gene expression dynamics during primordial germ cell induction from human induced pluripotent stem cells. Stem Cell Res Ther 11, 100 (2020). https://doi.org/10.1186/s13287-020-01620-y

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13287-020-01620-y

Keywords