Efficient differentiation of hPSCs towards the three germ layers
Primed human pluripotent stem cells (hPSCs) represent in vitro the human epiblast of the pre-gastrulating embryo and are a great model to recapitulate early human embryogenesis (Figure 1A,B). They self-renew and retain pluripotency -the ability to differentiate towards all the cell types of the future organism. Here, we direct PSC differentiation towards lineage-specific paths upon primed pluripotency exit through a chemical and growth factor approach (Figure 1B,C). The protocols described in Figure 1C have been implemented based on previous lineage specification methods using hPSCs, including those generating the primitive streak stage6,21,27, lateral and paraxial mesoderm6,21,27, definitive endoderm5,21,28, and neuroectoderm7,21.
Human pluripotent stem cells are characterized by the expression of core pluripotency factors, including OCT4 (Figure 1D). Upon reaching 70% confluency, removal of pluripotent medium conditions and modulation of specific signaling cascades trigger a downregulation of POU5F1 and NANOG mRNA levels, and the protein it codes, pluripotency factor OCT4, towards differentiation to any lineage tested (Figure 1D–F). Strong activation of Wnt/β-catenin signaling through GSK3 inhibition robustly induces a primitive streak-like cell state, characterized by the upregulation of Brachyury (TBXT) at both the mRNA and protein levels, with ~97% of cells quantified as TBXT+ (Figure 1D,E,G). Activation of WNT, BMP, and Nodal/Activin signaling to specify a posterior primitive streak-like cell, followed next day by inhibition of Wnt and TGFβ signaling while maintaining BMP activation, drives the formation of lateral mesoderm (LM)-like cells, characterized by high HAND1 protein and mRNA expression (Figure 1D,E,G). On the other hand, activation of Wnt signaling and Activin-A in the absence of BMP signaling drives the formation of an anterior primitive streak-like state. Subsequent inhibition of the BMP cascade while maintaining Wnt activity drives the formation of paraxial mesoderm (PM)-like cells expressing TBX6 protein and TBX6 mRNA levels (Figure 1D,E,G). Sustained activation of Nodal/Activin-A after WNT activity peak at the PS stage induces differentiation of hPSCs towards definitive endoderm (DE), characterized by high mRNA and protein levels of SOX17 (Figure 1D,E,G). Last, maintenance during several days of hPSCs under a double TGFβ and BMP inhibition regime triggered the formation of neuroectodermal cells expressing PAX6 at the protein (~65% positive cells) and mRNA levels, which corresponds in vivo to the epithelial layer of the anterior epiblast at the onset of gastrulation (Figure 1A). Of note, the protocols described yield a robust and homogeneous differentiation towards post-implantation lineages (>92% positive for lineage-specific markers), except for neuroectodermal differentiation (~65% PAX6+ cells) (Figure 1D,E,G), with an efficiency similar to other studies5,6,7,27,28. Additionally, we provide a list of reference cell fate markers that characterize the lineages described here (Table 2).
Mapping genome and chromosome (in)stability in pluripotency
Alterations in chromosome structure and number (aneuploidy) represent a hallmark of tumorigenesis, which often impair cellular differentiation of hPSCs29,30,31,32,33, representing a hurdle for their application in cell therapy and regenerative medicine. Given the underlying chromosomal instability that can arise during long-term hPSC culture, driven by the acquisition and selection of recurrent chromosomal abnormalities and by endogenously elevated replicative stress22,32,34. Quantification of chromosomal integrity and DNA damage remains critical: as an experimental checkpoint to validate the health of the cells prior to starting differentiation protocols, and ii) for those interested in studying the mechanisms of genome integrity regulation in stem cells. Here, the utility of a combined chromosomal instability and DNA damage assessment workflow for distinguishing different forms of genome instability in hPSCs (Figure 2A) is demonstrated. As proof-of-concept perturbations, cells were exposed to a panel of replication stress-inducing agents and mitotic inhibitors (Figure 2B–O, Table 3).
Compared to somatic cells, hPSCs display fast cell cycle dynamics characterized by a shortened G1 phase and a high proportion of cells in S-phase, which makes them susceptible to replication stress (RS)22,35,36. Replication stress can arise from diverse mechanisms, including replication fork slowing, nucleotide depletion, and replication-transcription conflicts, ultimately leading to DNA damage. To model distinct forms of replication stress, cells were treated with aphidicolin (Aphi; 50 nM), an inhibitor of replicative DNA polymerases; hydroxyurea (HU; 1–3 mM), which depletes intracellular deoxyribonucleotide pools through inhibition of ribonucleotide reductase; or camptothecin (CPT; 1–2.5 µM), a topoisomerase I inhibitor that generates severe replication-associated DNA lesions (Figure 2B, Table 3). Low doses of aphidicolin and HU induced a modest increase in double-stranded DNA breaks, as measured by γH2AX (Figure 2C,D). In contrast, treatment with CPT produced substantially greater γH2AX accumulation, indicating enhanced replication-associated DNA damage (Figure 2C,D). Similarly, phospho-RPA levels were most strongly and consistently elevated following CPT treatment. HU-treated cells in S-phase displayed a bimodal response, with some cells showing a mild increase in phospho-RPA, while others remained unperturbed, but at the same threshold as aphidicolin-treated cells. Overall, the phospho-RPA levels in the HU-treated condition were comparable to those observed following aphidicolin treatment. (Figure 2C,D). Together, these results demonstrate that the workflow can quantitatively distinguish distinct degrees and types of replication stress, highlighting its utility for evaluating genome stability in hPSCs and their response to genotoxic stress.
Replication stress can give rise to incomplete DNA replication at vulnerable genomic regions, allowing cells that have escaped DNA repair mechanisms to enter mitosis with under-replicated DNA21,37. During chromosome segregation in mitosis (M-phase), unreplicated DNA can generate ultrafine bridges (UFBs), which are visualized as BLM+ threads connecting chromatids during segregation38,39. Treatment with aphidicolin at either 50 or 200 nM for 16 h increased the percentage of mitotic hPSCs displaying BLM-positive UFBs relative to the DMSO-treated control (Figure 2E–G), highlighting the sensitivity of this developmental stage to genomic instability. The largest mean increase, approximately threefold relative to control, was observed following treatment with 50 nM aphidicolin. These results demonstrate that mild replication stress can persist into mitosis as unresolved DNA structures in hPSCs. Intriguingly, we have recently shown that long-term culture of hPSCs also increases the rate of UFBs in mitosis20, linking basal replication stress with mitotic defects (Figure 2A).
Unrepaired genomic instability often causes mitotic catastrophe at the onset of cell division. Chromosome missegregation represents a phenotype characterized by the presence of “lagging” chromosomes lying between the separating chromatids during anaphase, which frequently occurs in cells with unresolved damage in S-phase. These events can lead to incorrect segregation of chromosomes into the daughter cells after mitosis, which can give rise to an imbalance in the structure/number of chromosomes (aneuploidy), which is frequent during long-term pluripotent stem cell culture20,21. Primed hPSCs display a basal rate of chromosome missegregation of ~5%–9%21(Figure 2H–J), which drastically increases to 30%–70% upon MPS1 checkpoint spindle inhibition via reversine in a concentration-dependent manner (Figure 2H–J). As a consequence, many of these lagging chromosomes cause nuclear envelope rupture and reside in the cytoplasm, forming micronuclei (MN) after mitosis. Micronuclei are often considered a hallmark and driver of chromosomal instability, since the presence of DNA outside the nucleus triggers an immune cellular response that can promote tumor progression and further replicative stress40,41. In pluripotent conditions, MN are scarce in interphase hPSCs (Figure 2K,L). While a low dose of reversine treatment (0.25 µM) increases the presence of MN to ~5% of the cells, 0.5 µM reversine triggers ~40% of hPSCs in interphase with micronuclei (Figure 2K,L). In conclusion, mitotic instability can be easily measured through immunofluorescence and static imaging by quantifying the percentage of anaphases displaying lagging chromosomes and the percentage of interphase cells with MN.
To functionally confirm gains and losses of chromosomes after genomic stress, we perform karyotype or Giemsa staining. This technique allows the spread of the compact chromosomes of metaphase-arrested cells for direct visual quantification after staining, confirming whether each cell is euploid (46 chromosomes) or aneuploid (imbalance in this number). Pre-treatment with reversine rendered a high percentage of aneuploid cells (~40%) compared to the control untreated (~20%) (Figure 2M–O). In conclusion, an array of methods is provided to map genomic stress and chromosomal instability across the human pluripotent stem cell cycle.
Troubleshooting and identification of suboptimal outcomes
Representative suboptimal outcomes are shown in Figure 3, which helps to troubleshoot the protocols described above. The two differentiated hPSC cultures shown in Figure 3A display high cell densities and overlaps between adjacent DAPI-positive nuclei. Following thresholding and watershed processing, groups of overlapping nuclei were merged into large, segmented objects rather than being identified as individual nuclei (Figure 3A, right pictures). In other regions, individual nuclei were fragmented into multiple nuclei or were not accurately segmented. Consequently, the resulting regions of interest did not reliably correspond to individual cells and were not suitable for quantitative measurement of nuclear fluorescence intensity.
When this occurs, the threshold, watershed, and particle-size parameters should first be reviewed, and the segmentation overlay should be checked to confirm that each region of interest corresponds to one intact nucleus. The same analysis criteria should be applied consistently across experimental conditions to avoid quantification bias. To reduce segmentation failure, the initial seeding density should be optimized to prevent excessive confluence, and representative fields containing distinguishable nuclei should be acquired without selectively excluding regions based on the experimental phenotype. If individual nuclei cannot be separated reliably after parameter adjustment, the image should be excluded from quantitative analysis, or cells should be manually segmented. In this scenario, more powerful software, including Ilastik or CellPose, may be of great help to avoid manual segmentation42,43.
Figure 3B shows representative suboptimal Giemsa-stained metaphase preparations. Instead of forming a single spread containing spatially separated chromosomes, the chromosomes remain clustered and overlap each other. Some fields also contain multiple chromosome clusters, cellular debris, and poorly defined or incompletely dispersed metaphase preparations. These features prevent reliable identification and counting of individual chromosomes and may result from inadequate hypotonic swelling, incomplete cellular rupture, excessive cell density on the slide, insufficient dropping force, or inconsistent dropping height. Additionally, it is not rare to observe metaphase spreads with condensation defects (~5%–10% in basal conditions culturing hPSCs, Figure 3C). We and others22 have graded these metaphases where chromosomes are highly condensed, making it impossible to quantify the chromosome number.
To improve metaphase-spread quality, the duration of hypotonic treatment should be kept consistent, the fixed-cell suspension should be resuspended gently but completely, and the cell concentration should be adjusted before slide preparation. The slide should be sufficiently cold and wet, and the suspension should be dropped at a consistent height and force. Only metaphase spreads containing clearly distinguishable, sufficiently condensed chromosomes within an imaging field should be included in chromosome-number quantification. Dense chromosome clusters, overlapping spreads, incomplete cellular rupture, and fields in which chromosome loss or mixing between adjacent spreads is suspected should be excluded.

Figure 1. In vitro differentiation of human pluripotent stem cells into early embryonic lineages and validation of cell fate identity. (A) Schematic representation of early human embryonic development until the gastrulation stage, highlighting the developmental context of the in vitro lineages derived in this study. Inspiration for the creation of this scheme comes from Shahbazi et al44. (B) Overview of the differentiation of hPSCs to generate primitive streak (PS), lateral mesoderm (LM), paraxial mesoderm (PM), definitive endoderm (DE), and neuroectoderm (NE)-like cells. (C) Experimental timeline and culture conditions used for human lineage specification. (D) Representative immunofluorescence experiment characterizing lineage-specific identity acquisition following differentiation. Nuclei were counterstained with DAPI. Primitive streak (PS)-like cells were identified by TBXT, lateral mesoderm (LM) by HAND1, paraxial mesoderm (PM) by TBX6, definitive endoderm (DE) by SOX17, and neuroectoderm (NE) by PAX6. Scale bars: 50 µm. (E) Quantification of lineage-associated marker-positive cells following differentiation. Primitive-streak-like, lateral-mesoderm-like, paraxial-mesoderm-like, definitive-endoderm-like, and neuroectoderm-like cells were identified by TBXT, HAND1, TBX6, SOX17, and PAX6 expression, respectively. The percentage of marker-positive cells was calculated relative to the total number of DAPI-positive nuclei in each field. Bars represent the mean ± SD of 10 imaging fields with n>1000 cells analyzed per lineage in one representative differentiation experiment. No inferential statistical comparisons were performed. (F) Representative experiment of a qPCR of hPSCs undergoing differentiation. Expression values were normalized to β-actin and plotted as relative expression (2^-ΔCt). Data are mean ± SD of n = 3 technical replicates. (G) Representative experiment of a qPCR of hPSCs undergoing differentiation. Expression values were normalized to β-actin and plotted as log₂ fold change (log₂FC) relative to the pluripotent state. Data are mean ± SD of n = 3 technical replicates. Figures 1A–C have been created using BioRender and are used under the terms of the applicable license. Please click here to view a larger version of this figure.

Figure 2. Mapping DNA damage and mitotic defects in hPSCs. (A) Schematic overview of the genomic and chromosomal instability phenotypes associated with each phase of the human cell cycle, including single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) breaks, ultrafine bridges (UFBs), chromosome missegregation, micronuclei formation, and aneuploidy. (B) Experimental workflow to assess DNA damage responses to replication stress during S-phase using Aphidicolin, hydroxyurea (HU), and camptothecin (CPT). (C) Representative immunofluorescence images of untreated and drug-treated hPSCs stained with DAPI, EdU, γH2AX, and phospho-RPA (pRPA). EdU-positive nuclei identify cells in S phase, whereas γH2AX marks DNA double-strand breaks. pRPA serves as a marker of DNA damage associated with single-stranded DNA and DNA replication stress, and can also be detected at sites of double-strand break repair. Scale bars: 50 µm. (D) Quantification of γH2AX and pRPA mean fluorescence intensity (MFI) in EdU-positive nuclei following treatment with HU, aphidicolin, or CPT, as shown in C. Data are from a representative biological experiment with >300 EdU-positive nuclei analyzed per condition. The number of EdU-positive nuclei analyzed were: control (n = 326), HU 1 mM (n = 314), HU 3 mM (n = 304), aphidicolin 50 nM (n = 377), CPT 1 µM (n = 426), and CPT 2.5 µM (n = 504). Values are shown in arbitrary units (a.u.). (E) Experimental workflow used to assess ultrafine bridge formation following aphidicolin treatment. (F) Representative immunofluorescence images of anaphase cells without and with BLM+ ultrafine bridges. Cells were stained with DAPI and BLM. Scale bars, 10 µm. (G) Quantification of the percentage of anaphase cells containing BLM-positive ultrafine bridges after treatment with aphidicolin. Data are mean ± SD of n = 3 independent biological replicates, with > 50 anaphases quantified per condition in each experiment; each dot represents one independent experiment. Statistical significance was assessed using repeated-measures one-way ANOVA followed by Dunnett’s multiple comparisons test, comparing each aphidicolin-treated condition with the DMSO control. *p < 0.05 and **p < 0.01. (H) Experimental workflow used to induce chromosome segregation defects and micronuclei formation. (I) Representative immunofluorescence images of control and reversine-treated cells stained for DAPI, Lamin B1, α-tubulin, and CREST. The arrow indicates a chromosome missegregation event. Scale bars = 10 µm. (J) Quantification of chromosome missegregation following reversine treatment and 7 h recovery. Data are mean ± SD of n = 3 independent biological replicates, with > 100 anaphases quantified per condition in each experiment; each dot represents one independent experiment. Statistical significance was assessed using one-way ANOVA followed by Dunnett’s multiple comparisons test. ****p < 0.0001. (K) Quantification of interphase cells displaying Lamin B1-positive micronuclei following reversine treatment. Data are mean ± SD of n > 3 independent biological replicates, with > 100 interphase cells quantified per condition in each experiment. Statistical significance was assessed using Welch’s ANOVA followed by Dunnett’s T3 multiple comparisons test. *p < 0.05, ***p < 0.001. (L) Representative immunofluorescence images of micronuclei stained for DAPI and Lamin B1. Arrows indicate Lamin B1-positive micronuclei. Scale bars: 10 µm. (M) Experimental workflow used for metaphase spread preparation and karyotype analysis. (N) Representative metaphase chromosome spreads showing cells with normal chromosome number, 46 chromosomes, and abnormal chromosome number, 43 chromosomes. Scale bars: 10 µm. (O) Quantification of cells with normal and abnormal chromosome number following the indicated treatments (n = 50 metaphase spreads per condition combined from n = 3 different biological experiments/3 independent coverslips). Figure 2A has been created using BioRender and is used under the terms of the applicable license. Please click here to view a larger version of this figure.

Figure 3: Troubleshooting and identification of suboptimal outcomes. (A) Representative unsuccessful nuclear segmentation of differentiated definitive endoderm cells. The high density and overlap of DAPI-positive nuclei caused adjacent nuclei to be merged into large objects in the segmentation, preventing reliable single-nucleus fluorescence quantification. Scale bar: 50 µm. (B) Representative suboptimal Giemsa-stained metaphase preparations showing clustered and overlapping chromosomes, incomplete chromosome spreading, and poorly separated metaphase spreads. These preparations were unsuitable for chromosome-number quantification because individual chromosomes could not be distinguished and counted confidently. Scale bar: 10 µm. (C) (Left) Representative metaphase spreads after Giemsa staining clustered in four groups (g1–g4) based on ascending condensation defects; (Right) Quantification of the left phenotypes (n = 1 independent biological replicate, with n > 100 metaphase spreads quantified). Scale bar: 100 µm. Please click here to view a larger version of this figure.
| Media Name | Basal Condition | Growth Factors & Inhibitors |
| SFM | Stem Flex supplemented with 1% P/S | - |
| PS | Advanced RPMI supplemented with 1X L-Glutamine and 1% P/S | 5 µM CHIR99021 |
| LM | Advanced RPMI supplemented with 1X L-Glutamine and 1% P/S | LM-A Media: |
| 6 µM CHIR99021 |
| 30 ng/ml Activin A |
| 20 ng/ml FGF2 |
| 100 nM PIK90 |
| 40 ng/ml BMP4 |
| LM-B Media: |
| 1 µM A8301 |
| 20 ng/ml BMP4 |
| 1 µM C59 |
| PM | Advanced RPMI supplemented with 1X L-Glutamine and 1% P/S | PM-A Media: |
| 4 µM CHIR99021 |
| 30 ng/ml Activin A |
| 20 ng/ml FGF2 |
| 100 nM PIK90 |
| PM-B Media: |
| 3 µM CHIR99021 |
| 1 µM A8301 |
| 20 ng/ml FGF2 |
| 250 nM LDN193189 |
| DE | Advanced RPMI supplemented with 1X L-Glutamine and 1% P/S | DE-A Media: |
| 2 µM CHIR99021 |
| 100 ng/ml Activin A |
| 50 nM PIK90 |
| DE-B Media: |
| 250 nM LDN193189 |
| 100 ng/ml Activin A |
| NE | E6 supplemented with 1X L-Glutamine and 1% P/S | 500 nM LDN193189 |
| 10 µM SB431542 |
Table 1: Media compositions used in this study for differentiation of hiPSCs into early embryonic lineages. The ingredients of the media for each early embryonic lineage, with the respective concentrations, are indicated in that table.
| Cell State | Core markers for cell fate validation |
| Primded pluripotency | OCT4, NANOG, SOX2, TRA-1-60, SSEA4, LIN28A, DNMT3B, CDH1 |
| Primitive streak (PS) | TBXT, MIXL1, EOMES, SNAI1, CDX2, GSC |
| Lateral mesoderm (LM) | HAND1, FOXF, ISL1, KDR/VEGFR2, GATA4 |
| Paraxial mesoderm (PM) | TBX6, MSGN1, MEOX1, PAX3, FOXC1 |
| Definitive endoderm (DE) | SOX17, FOXA2, CXCR4, GATA4, GATA6, HNF4A |
| Neuroectoderm (NE) | PAX6, SOX1, SOX2, NESTIN, OTX2, ZIC1, ZEB2 |
Table 2: Marker panels for validation of early human cell lineages. Early human embryonic lineages and their markers for cell fate validation are indicated in that table. The markers listed in this table are not necessarily restricted to a single lineage when considered individually. Cell identity should therefore be checked using the combined expression of multiple positive markers, together with loss of pluripotency-associated markers and, where necessary, the absence of markers associated with alternative lineages. Marker interpretation should also consider the differentiation stage and signaling conditions. For instance, CDX2 is listed as a primitive streak (PS) marker, but it is not unique to this developmental lineage. CDX2 is also expressed in nascent posterior (lateral) mesoderm, extraembryonic mesoderm, as well as in amnion and trophoblast-like cells, depending on the species and the developmental time and context.
| Compound | Concentration | Mechanism/Funciton | Cell Cycle Phase Affected |
| Aphidicolin | 50–200 nm | DNA polymerase inhibitor. It leads to slowed DNA synthesis and replication fork stalling. | S phase/early S phase |
| Hydroxyurea (HU) | 1 mM, 3 mM | Inhibits ribonucleotide reductase, reducing dNTP pools and slowing/stalling DNA replication forks. | S phase |
| Camptothecin (CPT) | 1 µM, 2.5 µM | Inhibits topoisomerase I by trapping Topoisomerase I–DNA cleavage complexes. Collision of replication forks with these trapped complexes converts the lesions into replication-associated DNA damage, including single-strand breaks (SSBs) and double-strand breaks (DSBs). | S phase |
| Reversine | 0.1 µM, 0.25 µM, 0.5 µM | ATP-competitive drug that binds to ATP pocket of monopolar spindle 1 kinase (MPS1) and thus inhibiting the spindle assembly checkpoint (SAC). Blockage of SAC causes improper chromosome missegregation and premature dividing which leads to aneuploidy. | Mitosis |
Table 3: Compounds used in this study to induce replication stress and DNA damage in hiPSCs. Compounds used for inducing DNA errors are listed alongside their brief respective mechanism, concentrations, and target cell cycle phases.
Supplementary Table 1: Primer sequences. The primer sequences used in this study are listed in this table. Please click here to download this file.
Supplementary File 1: hPSC culture and image analysis approach. The procedure to freeze, thaw, and passage hPSCs and the image analysis approaches are described here. Please click here to download this file.