$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Characterization of undifferentiated hiPSCs
To assess the phenotype of the hiPSCs, the analysis of colony/cell morphology, the determination of PSC-specific markers, and the examination of gene expression and alkaline phosphate activity should be performed. Undifferentiated hiPSCs should be round, with large nucleoli and without abundant cytoplasm. The majority of the colonies should be characterized by a flat and tightly-packed morphology, indicative of an undifferentiated phenotype (Figure 2A and Figure 2B). Additionally, more than 80% of the colonies should be positive for alkaline phosphatase activity staining (Figure 2C).
About 80% of the cells should be positive for classic pluripotency-related markers, such as Oct4, SSEA3, SSEA4, and Tra1-60 (Figure 2D-H), as shown by immunocytochemistry and flow cytometry, whilst percentages of nestin+ and β-III-Tubulin+ cells should be significantly low (about 8% and 3%, respectively, as shown in Figure 2H). These results should be reproducible over passages.

Figure 2. Characterization of Undifferentiated IMR90-hiPSCs. (A and B) Representative phase-contrast images (10X and 20X magnifications) of undifferentiated IMR90-hiPSC colonies. (C) representative images of alkaline phosphatase-stained colonies (4X magnification). (D-F) Representative immunocytochemical images of (D) Oct4 (red), (E) SSEA3 (green), and (F) Tra1-60 (green). (G) Representative dot plot of SSEA1 (CD15) and SSEA4 staining, analyzed by flow cytometry. (H) The bar graph shows the percentages of Oct4+ (~75 - 80%), Tra1-60+ (~75 - 80%), SSEA3+ (~75 - 80%), nestin+ (~10 - 15%), and β-III-tubulin+ (~3 - 7%) cells, counterstained with DAPI and quantified by HCI, with a mean of 3 to 5 biological replicates ± the S.E.M. (graph modified from Reference 6). Please click here to view a larger version of this figure.
Assessment of pluripotency via EB formation
HiPSCs are pluripotent, meaning that they express three germ layer-related genes under appropriate conditions. To assess hiPSC pluripotency, it is possible to apply a common approach based on spontaneous EB formation, which induces the formation of the three germ layers14. Analyses of germ layer-specific genes must indicate a time-dependent increase of endoderm (α-fetoprotein (AFP) and Cytokeratin 18 (KRT18)), ectoderm (Nestin, SRY-box 1 (Sox1), and paired box 6 (Pax6)), and mesoderm (natriuretic peptide A (NPPA) and Brachyury-T) -related gene expression (Figure 3A and Figure 3B); see the Table of Materials.

Figure 3. Assessment of Pluripotency by EB formation. (A) Representative phase-contrast image of EBs at day 1. (B) The bar graph shows qPCR analyses of mesodermal (NPPA and brachyury), ectodermal (Pax6, Sox1, and nestin), and endodermal (AFP and KRT18) genes, normalized to the reference genes, β-actin and GAPDH, and calibrated to undifferentiated cells (Day 0). This is the ΔΔCt method, with a mean of 5 independent analyses ± the S.E.M. * p < 0,05, ** p < 0.01, *** p < 0.001; graph modified from Reference 6. Please click here to view a larger version of this figure.
Induction of neuronal and glial differentiation
IMR90-iPSCs can be differentiated into mixed cultures of post-mitotic neurons and glial cells following the steps summarized in Figure 1 and in the protocol sections. 5-8 days after plating the EBs onto standard matrix- or laminin-coated dishes in the presence of complete NRI medium, rosette-like structures should start becoming visible (Figure 4A). Rosettes are characterized by the presence of nestin+ cells (neuronal precursors, ~90%), with few β-III-tubulin+ cells (committed neuronal cells, ~5 - 10%), the latter generally localized mainly at the periphery of the rosettes (Figure 4B, rosettes on Day 12).
Upon rosette dissociation and replating onto laminin- or standard matrix-coated dishes or plates in the presence of complete ND medium, cells start undergoing differentiation into mixed cultures of neurons and glia, progressively forming clusters of neuronal cell bodies connected by bundles of neurites (Figure 4C and Figure 4D). Similar results should be obtained when analyzing neuronal populations obtained by expanding rosette-derived NSCs and differentiating them into neurons and glia. NSCs expanded from rosettes should be nestin+ (Figure 4E, inset showing nestin+ cells).
After 21 days of differentiation, the cells should be positive for β-III-tubulin; NF200; Tau; and MAP2, the late marker of dendrites (Figure 4D, 4F, and Figure 4H), with at least 10 - 15% of the cells positive for the glial fibrillary acidic protein (GFAP), an astroglial marker (Figure 4G and Figure 4H). Moreover, ~20 - 30% of the cells should retain the expression of nestin after differentiation (Figure 4H). It is important to consider that the percentage of each cell type (i.e., neuron, astrocyte, nestin+ cells) may vary over passages, and user-dependent variability may be observed.
By analyzing specific neuronal subpopulations, GABAergic neurons represent ~15 - 20% of the total cells, dopaminergic neurons ~13 - 20%, and glutamatergic neurons ~35 - 42%, as shown by immunostaining for gamma-aminobutyric acid (GABA), tyrosine hydroxylase (TH), and vesicular glutamate transporter 1 (VGlut1), respectively (see the representative quantification in Figure 4H). The induction of differentiation can also be assessed by the analysis of pluripotency-related markers (e.g., Oct4, Tra1-60, and SSEA3), which should be significantly downregulated in differentiated cells versus undifferentiated hiPSCs (not shown, see Reference 6). This can also be confirmed through the analysis of gene expression by qPCR, which should indicate a decrease of Oct4 and Nanog and the upregulation of neuronal genes, such as neural cell adhesion molecule 1 (NCAM1) and microtubule-associated protein 2 (MAP2); the presynaptic gene, synaptophysin (SYP); and the post-synaptic gene, microtubule-associated protein tau (MAPT), as shown in Figure 4I. Furthermore, dopaminergic (TH and NR4A1), noradrenergic (PHOX2A and PHOX2B), glutamatergic (NARG2, GRIA1 and GAP43), GABAergic (GABRA1 and GABRA3), motor neurons (ISL1 and LHX3), and cholinergic (SLC5A7 and SLC18A13) -related genes result in upregulated neuronal cells compared to undifferentiated cells (Figure 4J).
The analysis of spontaneous electrical activity, by means of MEA, is a valuable readout to assess the functionality of the neuronal network in differentiated hiPSCs. At the end of the differentiation period, neuronal derivatives are generally characterized by a mean firing rate (MFR) of at least 60 spikes/min (see the representative raster plot in Figure 4K). However, bursts are not observed.

Figure 4. Differentiation of IMR90-hiPSCs into Mixed Cultures of Neurons and Glia. (A and B) Representative images of rosettes after 7 DIV (A) and after 12 DIV (B), stained for nestin (green) and β-III-tubulin (red)). (C and D) Representative images of differentiated cells after 22 DIV (C) and 28 DIV (D), stained for β-III-tubulin (red) and NF200 (green)). (E) Representative image of NSCs derived from rosette dissociation and expansion (inset shows nestin+ cells, red). (F and G) Representative images of neuronal cells (F, stained for NF200 (red) and Tau (green)) and glial cells (G, stained for GFAP (red)) differentiated from NSCs (after 21 DIV). (H) Quantification of nestin, MAP2, GFAP, gamma-aminobutyric acid (GABA), vesicular glutamate transporter 1 (VGlut1), and tyrosine hydroxylase (TH) -positive cells by HCI, comparing IMR90-hiPSC-derivatives and cells differentiated from IMR90-hiPSC-derived NSCs (graph modified from Reference 7). (I and J) Bar graphs showing qPCR analyses of pluripotency genes (Oct4 and Nanog) and neuronal genes (NCAM1, MAP2, SYP, and MAPT) (I) and of dopaminergic (TH and NR4A1), noradrenergic (PHOX2A and PHOX2B), glutamatergic (NARG2, GRIA1, and GAP43), GABAergic (GABRA1 and GABRA3), motor neuron (ISL1 and LHX3), and cholinergic (SLC5A7 and SLC18A13) -related genes (J). All analyses were normalized to the reference genes, β-actin and GAPDH, and calibrated to undifferentiated cells (green bars). This is the ΔΔCt method, with a mean of 5 independent analyses ± the S.E.M. * p < 0,05, ** p < 0.01, *** p < 0.001. The graphs in I and J were modified from Reference 6. (K) Representative raster plot of IMR90-NSC-derived neurons (the recording was carried out for a minimum of 600 s; the vertical bars represent single spikes). Please click here to view a larger version of this figure.
A specific signature of neuronal markers is upregulated in differentiated IMR90-hiPSCs
In the new toxicity testing paradigm, it is essential to define the molecular and cellular events occurring within a cell following exposure to a given toxicant. Therefore, it is relevant to characterize which signaling pathways are activated and/or upregulated within the cellular model under investigation.
Commercially available arrays for the analyses of protein kinase gene expression can be used to compare undifferentiated hiPSCs versus differentiated cells. Differentiated IMR90-hiPSCs undergo the upregulation of genes involved in the control of neurotrophin receptors, axon guidance regulation, neurite outgrowth modulation, glial-derived neurotrophic factor (GDNF) receptors, bone morphogenetic protein (BMP)/TGF-beta pathway, and platelet-derived growth factor (PDGF) receptors (Figure 5A and Figure 5B).
The analysis of RPPA shows the upregulation of a specific neuronal signature in differentiated IMR90-hiPSCs. In particular, Erk/CREB, Akt/PDK1/mTOR, and Notch1 signaling pathways are activated upon differentiation (Figure 5C).

Figure 5. Differentiated Neuronal and Glial Cells Exhibit the Activation of Neuron-related Pathways. (A and B) Bar graphs report the qPCR analyses of neuron-related kinases (A) and other kinase-related genes (B). Gene expression data were normalized to the reference genes 18S and GAPDH (provided in the array) and calibrated to undifferentiated cells. For these analyses, a gene was considered significantly upregulated when its expression was at least 2x higher than in undifferentiated cells (2-ΔΔCt≥ 2); mean of 3 independent analyses ± the S.E.M.). (C) The bar graph shows the absolute protein quantifications by means of RPPA analysis, comparing differentiated (red bars) and undifferentiated cells (green bars). The proteins belonging to the same signaling pathway cascades are clustered together as follows: CREB/Erk/Akt, PDK1/mTOR/Erk/Akt, Notch1, Shh and Wnt, SAPK/JNK, and BMP/SMAD. Mean ± the S.E.M. of 4 independent analyses. * p < 0,05, ** p < 0.01, *** p < 0.001; graph modified from Reference 6. Please click here to view a larger version of this figure.
IMR90-hiPSC-derived neuronal/glial cultures can be used to assess the effects of rotenone
Rotenone, an inhibitor of complex I of the mitochondrial respiratory chain, is known to cause oxidative stress by triggering the activation of the Nrf2 pathway. Under quiescent conditions, Nrf2 is anchored in the cytoplasm by Keap1 (Kelch-like ECH-associated protein 1), an Nrf2 repressor, which facilitates Nrf2 ubiquitination and proteolysis15. Upon the induction of oxidative stress, Nrf2 translocates into the nucleus and activates the expression of Nrf2-ARE target genes16.
IMR90-hiPSC-derived neurons and glial cells can be used to assess the effects of rotenone on Nrf2 activation by exposing the cells to different concentrations of rotenone (e.g., 1, 10, and 100 nM) for 24 h. These concentrations were established according to previous studies17,18.
At these concentrations and times of exposure, rotenone did not result in cytotoxicity, as shown by the quantification of live DAPI+ cell nuclei (Figure 6A). Rotenone induced Nrf2 nuclear translocation, especially after exposing the cells to 100 nM rotenone (Figure 6B and Figure 6E). At the same concentration, a significant decrease in cytoplasmic Keap1 was observed (Figure 6C), together with an increase in both NAD(P)H quinone oxidoreductase 1 (NQO1) and Sulfiredoxin 1 (SRXN1), two Nrf2-target enzymes19,20 (Figure 6D).

Figure 6. Effects of Rotenone on Nrf2 Nuclear Translocation, Keap1, SRXN1, and NQO1 protein levels. (A) Quantification of live DAPI+ cells (i.e., non-pyknotic nuclei) upon 24 h treatment with 1, 10, 100 nM rotenone and normalized to untreated cells (Control, Ctr). (B) Nrf2 protein nuclear translocation (i.e., nuclear/cytoplasmic ratios) after 24 h of exposure to rotenone, assessed by measurements of fluorescence intensity using HCI analysis. (C) Quantification of cytoplasmic Keap1 protein levels upon rotenone treatment, assessed by HCI analysis. (D) Quantification of NAD(P)H quinone oxidoreductase 1 (NQO1) and Sulfiredoxin 1 (SRXN1) by means of immunofluorescence and HCI after 24 h of treatment with rotenone. (E) Representative images of Nrf2 protein localization (green). All values are shown as the mean ± S.E.M. of 3 biological replicates. * p < 0,05, ** p < 0.01; figure modified from Reference 7. Please click here to view a larger version of this figure.
At these concentrations and times of treatment, rotenone also elicited a concentration-dependent increase of astroglial (GFAP+) cell percentage (Figure 7A and Figure 7B), without affecting the proportions of NSCs (nestin+) and neurons (MAP2+) (Figure 7B). By looking at the proportions of specific neuronal subtypes, rotenone treatment (10 nM and 100 nM) significantly decreased the number of dopaminergic neurons (TH+) (Figure 7C and D), whilst the percentages of GABAergic (GABA+) and glutamatergic (VGlut1+) neurons did not change (Figure 7D). Analogously, previous in vivo and in vitro studies have described a rotenone-dependent and selective dopaminergic neuronal cell death21,22,23.

Figure 7. Effects of Rotenone on Glial Cells and Dopaminergic Neurons. (A) Representative pictures of GFAP+ cells (red), with 40X magnification in the insets, untreated or treated with 100 nM rotenone for 24 h. (B) Quantification of nestin+, MAP2+, and GFAP+ cells, normalized to untreated cells (control, Ctr). (C) Representative pictures of dopaminergic TH+ neurons (green), with 40X magnification in the insets, untreated or treated with 100 nM rotenone for 24 h. (D) Quantification of GABA+, VGlut1+, and TH+ neuronal cells, normalized to untreated cells (Ctr). All values are shown as the mean ± S.E.M. of 3 biological replicates. * p < 0,05, ** p < 0.01, *** p < 0.001; figure modified from Reference 7. Please click here to view a larger version of this figure.
Statistical significance was assessed by one-way ANOVA with Dunnett's Multiple Comparison Test as a post-test (comparing all columns versus the control column)24 or by two-tailed unpaired or paired t-test according to the type of analysis. All data represent the average of at least three biological replicates ± the standard error of the mean (S.E.M.). An asterisk over a bar indicates a significant difference with the control group. * p < 0.05, ** p < 0.01, *** p < 0.001.
Table 1.
Note on dissociated rosette plating density: if rosette fragments do not look completely dissociated, in order to reach a cell plating density of about 15.000 cells/cm2, dissociated rosettes fragments deriving from about 50 EBs/1 x 60mm-dish can be resuspended in 50 mL of complete NRI medium and plated as follows (depending on the plate format):
| Multiwell plates / MEA | Growth area (cm2/well) | Volume of cell suspension to plate per well (or MEA chip) | Maximum number of plates that can be plated (with 50 ml of cell suspension) |
| 96 wells | 0.3 | 100 ul | 5 |
| 48 wells | 0.7 | 220 ul | 4 |
| 24 wells | 2 | 625 ul | 3 |
| 12 wells | 4 | 1.25 ml | 3 |
| 6 wells | 10 | 3.125 ml | 2 |
| Single well MEA chip | 3.5 | 1.1 ml | 45 |
Table 2: Acceptance Criteria
| Marker / Antibody | Percentage (on DAPI+ (live) cells) after 28 DIV |
| B-III-tubulin (Tuj1) | 35-45% |
| MAP2 | 50-60% |
| NF200 | 45-55% |
| GFAP | 10-25% |
| Nestin | 15-25% |