- Open Access
The role of insulin as a key regulator of seeding, proliferation, and mRNA transcription of human pluripotent stem cells
Stem Cell Research & Therapyvolume 10, Article number: 228 (2019)
Human-induced pluripotent stem cells (hiPSCs) show a great promise as a renewable source of cells with broad biomedical applications. Since insulin has been used in the maintenance of hiPSCs, in this study we explored the role of insulin in culture of these cells.
We report conditions for insulin starvation and stimulation of hiPSCs. Crystal violet staining was used to study the adhesion and proliferation of hiPSCs. Apoptosis and cell cycle assays were performed through flow cytometry. Protein arrays were used to confirm phosphorylation targets, and mRNA sequencing was used to evaluate the effect of transcriptome.
Insulin improved the seeding and proliferation of hiPSCs. We also observed an altered cell cycle profile and increase in apoptosis in hiPSCs in the absence of insulin. Furthermore, we confirmed phosphorylation of key components of insulin signaling pathway in the presence of insulin and demonstrated the significant effect of insulin on regulation of the mRNA transcriptome of hiPSCs.
Insulin is a major regulator of seeding, proliferation, phosphorylation and mRNA transcriptome in hiPSCs. Collectively, our work furthers our understanding of human pluripotency and paves the way for future studies that use hiPSCs for modeling genetic ailments affecting insulin signaling pathways.
Human-induced pluripotent stem cells (hiPSCs) can be expanded continuously and can differentiate into all lineages of the human body. Due to these characteristics, hiPSCs provide a renewable source of cells for developmental studies, drug screening, and regenerative medicine . Additionally, undifferentiated hiPSCs can recapitulate the aspects of genetic diseases provided that the affected cellular and molecular pathways are active in the pluripotent cells [2,3,4]. Since the initial description of human-induced pluripotent stem cells [5, 6], there has been a growing interest in identification of essential components of their culture. Studying these key components and the respective signaling pathways improves our understanding of pluripotency and also demonstrates the potentials and limitations of hiPSCs in modeling of diseases associated with cell signaling.
Transforming growth factor beta 1 (TGFB1) and fibroblast growth factor 2 (FGF2) are two of the components that have been identified that seem to be essential in the maintenance of human pluripotent stem cells. Another essential component is insulin, a peptide hormone that regulates an intricate network of signaling pathways and cellular responses  and has been identified as one of the key ingredients in culture of human pluripotent stem cells [8,9,10]. However, the role of insulin in the culture of hiPSCs is not fully understood. In this study, we report conditions for insulin starvation of hiPSCs for up to 3 days and use this setting to explore the effect of insulin on several cellular and molecular functions that have been associated with insulin signaling as well as explore the effect of insulin on global mRNA transcription in hiPSCs.
Culture of hiPSCs
hiPSC lines were derived from PBMCs using integration-free Sendai virus to deliver the Yamanaka factors as reported earlier . These cell lines were initially expanded on matrigel in the presence of mTeSR1 (StemCell Technologies) and passaged using 0.5 mM EDTA. All cells were then thawed and adapted to TeSR-E8 for at least 2 weeks prior to the experiments. The medium was changed daily starting within the initial 3 days of culture. To prepare an insulin-free hPSC medium, mTeSR-E5-defined medium (StemCell Technologies) was supplemented with FGF2 (bFGF, PeproTec, final concentration of 100 μg/l) and TGFB1 (PeproTec, final concentration of 2 μg/l). Cells were washed 4 times with DMEM-F12 (Life Technologies) prior to insulin starvation experiments. When insulin treated groups were needed for the experiments, the medium was further supplemented with insulin (insulin, human recombinant, zinc solution; Life Technologies) to the final concentration of 19.4 mg/l. During maintenance in TeSR-E8 and experiments, the cultures were fed double medium from day 4 of culture onwards, unless the cells were used for overnight experiments on day 4.
Of the four hiPSC lines initially used in the preliminary experiments (Additional file 4: Table S1), one line showed sporadic patches of differentiating colonies (data not shown). While the vast majority of the cells maintained their pluripotent morphology in the line, we excluded this line from our experiments to prevent the chance that traces of differentiated cells interfere with our data analysis.
Primary and secondary antibodies used for pluripotency were obtained from a Pluripotent Stem Cell 4-Marker Immunocytochemistry Kit (Life Technologies). DAPI was used for counter staining of the nuclei. Secondary antibody stainings were used as controls (data not shown).
Crystal violet staining
Crystal violet staining was used for visualization of the cells in quantification of cell adhesion and cell expansion. Samples were fixed with 4% PFA in PBS for 5 min and were stained with 0.05% crystal violet solution for 30 min; plates were rinsed up to 4 times with water, dried, and were used for imaging. The staining in the wells was quantified using microscopic imaging of the wells for seeding experiments, or direct images of the wells for expansion experiments. The percentage of surface area covered by cells (stained regions) was measured using photoshop and was compared between groups.
Cell cycle assay
Cells were harvested with accutase (Millipore) or Stempro Accutase (Life Technologies), fixed with 70% ethanol and stained through incubation of cells in PBS-0.1% Triton-X solution supplemented with 0.2 mg/ml RNase A and 0.02 mg/ml propidium iodide (PI) at room temperature for 30 min. Prepared samples were analyzed through flow cytometry. FlowJo software (Treestar, Inc) was used for modeling and analysis of the cell cycle.
Cells were harvested with accutase. Supernatant and rinse solutions were collected separately and were pooled with the harvested cells prior to staining. Samples were stained using dead cell apoptosis kit with Annexin V FITC and PI (Life Technologies) and analyzed through flow cytometry.
Human Phospho-RTK Antibody Array, Human Phospho-MAPK Antibody Array, Human Phospho-Kinase Antibody Array, and Human Apoptosis Antibody Array Kits were obtained from R&D system and used according to the provided instructions. Briefly, lysates were harvested after treatment of cells in the described experimental conditions. The cells were rinsed once with PBS (with Ca/Mg) and lysed with inhibitors supplemented with lysis buffer. Concentration of protein lysates were subsequently determined by BCA assay. The membranes were incubated overnight with protein lysate at 4 °C on a rocking platform before incubation with a secondary antibody or directly with a detection reagent based on the array manual provided with the kit.
Li-Cor Odyssey imager and Image Studio software (Li-Cor Biosciences) were used for imaging and quantification of protein array results. Image Studio software was used to quantify the signal of the array dots. The data were normalized in each dataset through division of signal value obtained from groups treated without insulin to their respective insulin-treated groups.
RNA extraction and RNA sequencing
RNA extraction and RNA sequencing were performed as described earlier (Carcamo-Orive et al., 2017). Sequencing reads were first quality checked by FastQC (v 0.11.8). Reads above Q30 were further mapped and reads per gene were counted by STAR (v2.7). Read count normalization and differential analyses were performed using Deseq2 R package (v1.22.2) with custom codes. Pathway enrichments were performed using Gene Set Enrichment Analysis (GSEA) with MSigDB 6.2 database. For 24 h, all significantly (p_adj < 0.05) differential genes were included for the enrichment analysis, while for 72 h, the significantly (p_adj < 5E10–4) differential genes were included to obtain a similar number of significant genes as compared to those under 24 h for the enrichment analysis.
The effect of insulin on seeding of hiPSCs
To investigate the effect of insulin on seeding of hiPSC aggregates, colonies were harvested by EDTA and were seeded with or without insulin. The attachment of the hiPSC aggregates was evaluated microscopically after crystal violet staining of plates 1 day after seeding. The presence of insulin significantly improved the seeding of hiPSC aggregates (Fig. 1a, b).
To determine the effect of insulin on seeding of dissociated hiPSCs, the colonies were digested with Accutase and passed through the 40-μM cell strainer and were used for crystal violet staining. The presence of Y-27632, which is an inhibitor of rho-associated kinase (ROCK), prevents apoptosis in dissociated hiPSCs . Our results demonstrate that the presence of both Y-27632 and insulin are essential for attachment of dissociated hiPSCs (Fig. 1c, d). Inclusion of either component allowed very limited attachment of cells while no attached cells were observed when both reagents were absent.
The effect of insulin on pluripotent morphology, expression of pluripotency markers, and expansion of hiPSC colonies
To compare the effect of insulin on pluripotent morphology and expression of pluripotency markers, 3 lines of hiPSCs (Additional file 4: Table S1) were maintained in culture for 3–4 days and were subsequently cultured with or without insulin for 3 additional days. These cells maintained pluripotent morphology (Fig. 2a) and expression of key pluripotency markers SOX2, POU5F1, and TRA-1-60 after 3 days of insulin starvation (Fig. 2b).
The surface area of the tissue culture plate covered by hiPSC colonies in the presence of insulin was significantly larger than the surface area covered in the absence of insulin (Fig. 2c, d). These observations were in line with number of cells in each condition, determined by manual cell count (Fig. 2e) and indicate that the presence of insulin significantly improves the proliferation of hiPSCs.
Insulin changes the cell cycle profile and inhibits apoptosis in hiPSCs
To investigate the effect of insulin on cell cycle, the cells were cultured with or without insulin for 3 days and then used for flow cytometry-based cell cycle assay using propidium iodide (PI) staining. There was a significant decrease in the percentage of the cells at the G1 phase is the absence of insulin (Fig. 3a, b).
After hiPSCs were cultured overnight (16–18 h) with and without insulin, there was a significant increase in detected apoptotic cells in the absence of insulin in a flow cytometry-based apoptosis assay (Fig. 3c, d).
Protein array analysis of apoptotic proteins after 16 h of insulin starvation demonstrates more than a 3-fold increase in levels of cleaved caspase-3 while levels of pro-caspase-3 remain at similar levels (Fig. 3e, f). This robust increase in levels of cleaved caspase-3 is accompanied by reduction in XIAP, TNFRSF10B proteins, and serine 392 phosphorylated TP53.
Verification of insulin-mediated protein phosphorylation in hiPSCs
We used protein arrays to monitor phosphorylation of a set of tyrosine kinase receptors and phosphokinases (Fig. 4). We confirmed that the insulin receptor (INSR) and insulin-like growth factor 1 receptor (IGF1R) are highly phosphorylated in the presence of insulin. Furthermore, our results indicate that phosphorylation of AKT and WNK1 and GSK3 are regulated by insulin stimulation in hiPSCs.
Insulin regulates transcriptome of pluripotent stem cells
To gain an understanding of the transcriptional changes induced by insulin in hiPSCs, we looked at the effect of insulin on gene expression at 24 h and 72 h.
An average of 62.7 million reads were obtained for each sample, and the average genome mapping rate was 90% (Additional file 1: Figure S1 and Additional file 2: Figure S2). The number of identified genes is shown in Additional file 4: Table S2.
The prominent effect of insulin on global mRNA transcription of the cells is demonstrated by the t-SNE plot (Fig. 5a) and ternary plot (Fig. 5b). The t-SNE plot clearly indicates that the presence of insulin has a stronger global effect on transcription than internal variations across the cells at 72 h. This is manifested by closer clustering of samples based on insulin treatment rather than the individual lines at 72 h time point. A ternary plot shows that among the three major variables (time, cell line, or treatment), insulin treatment explains most of the variances at a majority of the mapped genes.
We identified 1161 differentially expressed genes at 24 h of hiPSC culture with and without insulin and 6753 differentially expressed genes at 72 h, with q value < 0.05. The 20 top differentially expressed genes at 24 h and 72 h time points are listed in Additional file 4: Tables S3 and S4. The relationship between log2 fold change and p values is shown as volcano plots for each time point with differentially expressed genes shown in red (Fig. 5c). The heatmaps of differentially expressed genes at 24 h and 72 h are shown in Fig. 5d.
We compared the effect of insulin on gene expression at 24 h and 72 h. There is a strong correlation between gene expression at 24 h and 72 h suggesting that the effect of insulin on transcriptome is generally consistent between 24 h and 72 with more prominent effect at 72 h (Fig. 6a). This consistency is also demonstrated by overlap of genes differentially expressed at 24 h with 72 h (Fig. 6b).
To better understand the affected biological processes using identified differentially expressed genes, we performed Gene Ontology (GO) term enrichment for biological processes at each time point.
The top 25 terms that were associated with differentially expressed transcripts at 24 h are shown in Table 1. The top 25 biological Gene Ontology (GO) processes terms that were associated with significantly differentially expressed transcripts at 72 h are shown in Additional file 4: Table S5. From the list of top 20 GO processes at 24 h, we selected the high-ranked terms associated with our results and heatmaps of differentially expressed genes within each term were generated (Additional file 3: Figure S3).
We compared differentially expressed genes at 24 h and 72 h with KEGG pathway database to investigate the effect of insulin on expression of known genes involved in pathways regulating the pluripotency (Additional file 4: Table S6) and insulin signaling pathway (Additional file 4: Table S7).
We investigated the effect of insulin on seeding of hiPSC aggregates or dissociated hiPSCs. Our results with dissociated hiPSCs demonstrate that insulin exerts its effect on their seeding independent of Y-27632. This observation is in line with an earlier report indicating that insulin is essential for survival of dissociated human pluripotent stem cells . We expand on this by also demonstrating the significant effect of insulin on seeding of hiPSC aggregates without use of either chemical inhibitors blebbistatin or Y-27632.
Our data indicate that insulin significantly improves the expansion of hiPSC colonies and cellular proliferation. These results are broadly in line with earlier reports , although we observe less drastic effect after insulin starvation in our experiments on the expansion of the cells, possibly because Chen et al. used dissociated cells rather than hiPSC aggregates in their experiments. In our work, we also investigated the effect of insulin in expansion of hiPSC colonies without using of blebbistatin or Y-27632 and separated the seeding experiments from expansion experiments.
Specific inhibition of IGF1R with shRNA in human embryonic stem cells (hESCs) has been reported to lead to their competitive disadvantage in proliferation , demonstrating the importance of IGF1R in expansion of human pluripotent stem cells. In addition, hiPSCs derived from insulin-resistant patients with mutations in their insulin receptor have demonstrated lower proliferation compared to their normal counterparts in the presence of insulin , suggesting that insulin signaling plays a regulatory role in proliferation of hiPSCs. Phosphorylation of IGF1R and INSR in the presence of conditioned medium containing high levels of insulin  or after insulin stimulation  has also been reported. Consistent with these studies, we observed robust phosphorylation of both the INSR as well as IGF1R. In addition, our data appear to demonstrate a more robust phosphorylation of these receptors, potentially due to different starvation conditions or higher concentrations of insulin used.
Cell division and cell death are two main contributory factors to cellular proliferation. We observed a significant increase in apoptosis by hiPSCs in the absence of insulin. This observation was consistent with the detection of higher levels of cleaved caspase-3. Caspase-3 plays a key role in execution stage of apoptosis. Interestingly, cleavage of pro-caspase-3 has been reported after inhibition of AKT phosphorylation in hiPSCs . In addition, we detected lower levels of TNFRSF10B after insulin starvation. TNFRSF10B has been reported to play a dual role in cell survival .
After insulin starvation, we detected lower levels of serine 392 phosphorylated TP53 and higher level of XIAP in an apparent contrast to the usual role attributed to them. Serine 392 phosphorylated TP53 has been reported to contribute to genotoxic stress and apoptosis . XIAP has also been identified as an anti-apoptotic molecule; nevertheless, upregulation of XIAP in response to environmental stress stimuli to induce apoptosis resistance has been reported in a cellular model . These observations can be potentially attributed to the coping mechanism in surviving cells in response to the activated apoptotic pathways.
Phosphorylation of AKT is a critical node in insulin signaling . AKTs play a central role in pluripotency and it has been demonstrated that overexpression of AKT improves the maintenance of pluripotent morphology in a primate pluripotent stem cell model . Our data confirms that AKTs including AKT1 and AKT2 are main targets of insulin signaling in hiPSCs, consistent with earlier reports . We also demonstrate that phosphorylation of GSK3B is a target of insulin stimulation in hiPSCs. Inhibition of phosphorylation of both AKT and GSK3B has been demonstrated to increase apoptosis in hiPSCs [13, 19].
Despite phosphorylation of AKT and GSK3B in the presence of insulin, we still detected a basal level of phosphorylation of these target proteins in the absence of insulin. This observation is in line with earlier reports of contribution of other growth factors including endogenously secreted ELABELA  and high concentrations of FGF  in the medium which contribute independently to activation of AKT signaling pathway in hiPSCs.
With regard to cell cycle profile of hiPSC populations after 3 days of insulin starvation, we observed a significant reduction in percentage of cells in the G1 phase. This change in cell cycle profile is similar to reports regarding inhibition of phosphorylation of either GSK3B (which we here demonstrate to be a target of insulin signaling in hiPSCs) or receptor tyrosine kinase ERBB2 in human pluripotent stem cells. In both studies, inhibition of phosphorylation led to decrease in G1 phase accompanying an increase in apoptosis or decrease in proliferation in these cells [10, 13]. This reduction in percentage of cells in the G1 phase alongside a decrease in the proliferation of the cells suggests a cell cycle delay or arrest at G2/M transition. Our data also demonstrates that WNK1 is strongly phosphorylated upon insulin stimulation in hiPSCs. WNK1 has been identified as a downstream target of phosphorylation by AKT/PKB after insulin and IGF-1 stimulation in other cell lines [21, 22]. The cellular significance of this phosphorylation in the context of hiPSC biology would be an interesting avenue for follow-up studies.
To our knowledge, we are the first group to demonstrate the robust regulatory effect of insulin on mRNA transcriptome of hiPSCs. Insulin regulates expression of wide array of genomic targets . Using our starvation protocol, we explored the effect of insulin on gene expression in pluripotent stem cells after 24 h and 72 h. Our data also indicate that the effect of insulin on the transcriptome at 24 h and 72 h seem to be generally in the same direction while becoming more robust overtime.
At the 24 h time point and as expected, the top 25 identified terms include terms consistent with the role of insulin in cell signaling and transduction. Furthermore, regulation of both cell proliferation and cell death are also among the top 25 identified terms supporting the role of insulin in cell signaling and cellular proliferation of hiPSCs. All these terms remain among the top 100 terms at 72 h.
We compared the differentially expressed genes at 24 h with identified genes in insulin signaling pathway (KEGG pathway). Among the 9 identified genes, 8 are also differentially expressed at 72 h. Among these 8, there are classic components of insulin signaling pathway including INSR, IRS1, IRS2, and CRK. These are expressed at lower levels in the presence of insulin, while signaling suppressor SOCS2 is expressed at higher levels in the presence of insulin. This observation is in line with a negative regulation feedback of components insulin signaling pathway in response to activated insulin signaling or increased sensitivity of hiPSCs cultured without insulin.
We also compared the differentially expressed genes at 24 h with identified genes involved in regulating pluripotency (KEGG pathway). All of the 7 identified differentially expressed genes are also differentially expressed at 72 h. From these genes, FZD2 and DVL3 belong to Wnt signaling pathway and APC2 is a regulator of the Wnt pathway. Expressions of FZD2 and DVL3 are lower while expression of APC2 is higher in the presence of insulin. These observations suggest that insulin might exert a regulatory effect on Wnt signaling pathway in hiPSCs. The role of Wnt signaling in human pluripotent cells remains controversial . An interesting identified gene is LIF, with a higher expression in the presence of insulin. LIF plays a critical role in maintenance of naïve pluripotent stem cells which represent an earlier developmental stage of pluripotency than hiPSCs (primed stage) . LIF has been shown to increase the proliferation in other human cells [26, 27] and is a potential inducer of proliferation in the presence of insulin. Another interesting target is BMP4. Expression of BMP4 appears to be significantly lower in the presence of insulin at both time points. Elevated levels of BMP4 can trigger differentiation of human pluripotent stem cells into mesodermal lineages [28, 29] although differentiation process is suppressed by cell contact . These observations suggest that insulin might maintain hiPSCs at a higher pluripotency stage and less susceptible to differentiation.
In summary, we confirmed that insulin plays a significant role in regulation of proliferation, apoptosis, protein phosphorylation, and mRNA transcription in hiPSCs. Our findings will improve our understanding of the role of insulin in hiPSC biology and are reasonably expected to contribute to better disease modeling experiments affecting the insulin signaling pathways using hiPSCs. We are now focusing on understanding the global transcriptional differences between hiPSCs of insulin-resistant and insulin-sensitive individuals in the presence of insulin.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. RNA-Sequencing data will be available through GEO database.
APC regulator of WNT signaling pathway 2
Bone morphogenetic protein 4
Disheveled segment polarity protein 3
Fibroblast growth factor 2
Frizzled class receptor 2
Human embryonic stem cells
Human-induced pluripotent stem cells
Human phosphokinase kinase
Insulin-like growth factor 1 receptor
Insulin receptor substrate
Leukemia inhibitory factor
Mitogen-activated protein kinase
Receptor tyrosine kinase
Suppressor of cytokine signaling 2
Transforming growth factor beta 1
Takahashi K, Yamanaka S. Induced pluripotent stem cells in medicine and biology. Development. 2013;140(12):2457–61.
Iovino S, Burkart AM, Kriauciunas K, Warren L, Hughes KJ, Molla M, et al. Genetic insulin resistance is a potent regulator of gene expression and proliferation in human iPS cells. Diabetes. 2014;63(12):4130–42.
Lemonnier T, Blanchard S, Toli D, Roy E, Bigou S, Froissart R, et al. Modeling neuronal defects associated with a lysosomal disorder using patient-derived induced pluripotent stem cells. Hum Mol Genet. 2011;20(18):3653–66.
Quarto N, Leonard B, Li S, Marchand M, Anderson E, Behr B, et al. Skeletogenic phenotype of human Marfan embryonic stem cells faithfully phenocopied by patient-specific induced-pluripotent stem cells. Proc Natl Acad Sci U S A. 2012;109(1):215–20.
Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, et al. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007;131(5):861–72.
Yu J, Vodyanik MA, Smuga-Otto K, Antosiewicz-Bourget J, Frane JL, Tian S, et al. Induced pluripotent stem cell lines derived from human somatic cells. Science. 2007;318(5858):1917–20.
Taniguchi CM, Emanuelli B, Kahn CR. Critical nodes in signalling pathways: insights into insulin action. Nat Rev Mol Cell Biol. 2006;7(2):85–96.
Chen G, Gulbranson DR, Hou Z, Bolin JM, Ruotti V, Probasco MD, et al. Chemically defined conditions for human iPSC derivation and culture. Nat Methods. 2011;8(5):424–9.
Ludwig TE, Levenstein ME, Jones JM, Berggren WT, Mitchen ER, Frane JL, et al. Derivation of human embryonic stem cells in defined conditions. Nat Biotechnol. 2006;24(2):185–7.
Wang L, Schulz TC, Sherrer ES, Dauphin DS, Shin S, Nelson AM, et al. Self-renewal of human embryonic stem cells requires insulin-like growth factor-1 receptor and ERBB2 receptor signaling. Blood. 2007;110(12):4111–9.
Carcamo-Orive I, Hoffman GE, Cundiff P, Beckmann ND, D'Souza SL, Knowles JW, et al. Analysis of transcriptional variability in a large human iPSC library reveals genetic and non-genetic determinants of heterogeneity. Cell Stem Cell. 2017;20(4):518–32 e9.
Watanabe K, Ueno M, Kamiya D, Nishiyama A, Matsumura M, Wataya T, et al. A ROCK inhibitor permits survival of dissociated human embryonic stem cells. Nat Biotechnol. 2007;25(6):681–6.
Romorini L, Garate X, Neiman G, Luzzani C, Furmento VA, Guberman AS, et al. AKT/GSK3beta signaling pathway is critically involved in human pluripotent stem cell survival. Sci Rep. 2016;6:35660.
Shlyakhtina Y, Pavet V, Gronemeyer H. Dual role of DR5 in death and survival signaling leads to TRAIL resistance in cancer cells. Cell Death Dis. 2017;8(8):e3025.
Castrogiovanni C, Waterschoot B, De Backer O, Dumont P. Serine 392 phosphorylation modulates p53 mitochondrial translocation and transcription-independent apoptosis. Cell Death Differ. 2018;25(1):190–203.
Holcik M, Yeh C, Korneluk RG, Chow T. Translational upregulation of X-linked inhibitor of apoptosis (XIAP) increases resistance to radiation induced cell death. Oncogene. 2000;19(36):4174–7.
Watanabe S, Umehara H, Murayama K, Okabe M, Kimura T, Nakano T. Activation of Akt signaling is sufficient to maintain pluripotency in mouse and primate embryonic stem cells. Oncogene. 2006;25(19):2697–707.
Ho L, Tan SY, Wee S, Wu Y, Tan SJ, Ramakrishna NB, et al. ELABELA is an endogenous growth factor that sustains hESC self-renewal via the PI3K/AKT pathway. Cell Stem Cell. 2015;17(4):435–47.
Hossini AM, Quast AS, Plotz M, Grauel K, Exner T, Kuchler J, et al. PI3K/AKT signaling pathway is essential for survival of induced pluripotent stem cells. PLoS One. 2016;11(5):e0154770.
Dalton S. Signaling networks in human pluripotent stem cells. Curr Opin Cell Biol. 2013;25(2):241–6.
Jiang ZY, Zhou QL, Holik J, Patel S, Leszyk J, Coleman K, et al. Identification of WNK1 as a substrate of Akt/protein kinase B and a negative regulator of insulin-stimulated mitogenesis in 3T3-L1 cells. J Biol Chem. 2005;280(22):21622–8.
Vitari AC, Deak M, Collins BJ, Morrice N, Prescott AR, Phelan A, et al. WNK1, the kinase mutated in an inherited high-blood-pressure syndrome, is a novel PKB (protein kinase B)/Akt substrate. Biochem J. 2004;378(Pt 1:257–68.
O'Brien RM, Streeper RS, Ayala JE, Stadelmaier BT, Hornbuckle LA. Insulin-regulated gene expression. Biochem Soc Trans. 2001;29(Pt 4:552–8.
Sokol SY. Maintaining embryonic stem cell pluripotency with Wnt signaling. Development. 2011;138(20):4341–50.
Onishi K, Zandstra PW. LIF signaling in stem cells and development. Development. 2015;142(13):2230–6.
Salehnia M, Fayazi M, Ehsani S. Leukemia inhibitory factor increases the proliferation of human endometrial stromal cells and expression of genes related to pluripotency. Int J Reprod Biomed (Yazd). 2017;15(4):209–16.
Broholm C, Laye MJ, Brandt C, Vadalasetty R, Pilegaard H, Pedersen BK, et al. LIF is a contraction-induced myokine stimulating human myocyte proliferation. J Appl Physiol (1985). 2011;111(1):251–9.
Chadwick K, Wang L, Li L, Menendez P, Murdoch B, Rouleau A, et al. Cytokines and BMP-4 promote hematopoietic differentiation of human embryonic stem cells. Blood. 2003;102(3):906–15.
Richter A, Valdimarsdottir L, Hrafnkelsdottir HE, Runarsson JF, Omarsdottir AR, Ward-van Oostwaard D, et al. BMP4 promotes EMT and mesodermal commitment in human embryonic stem cells via SLUG and MSX2. Stem Cells. 2014;32(3):636–48.
Tashiro S, Le MNT, Kusama Y, Nakatani E, Suga M, Furue MK, et al. High cell density suppresses BMP4-induced differentiation of human pluripotent stem cells to produce macroscopic spatial patterning in a unidirectional perfusion culture chamber. J Biosci Bioeng. 2018;126(3):379–88.
We would like to thank Paul Burridge, PhD, Jared Churko, PhD, and Joseph C. Wu, MD, PhD, for the defined medium for culture of hiPSCs, used for our preliminary experiments prior to the reported experiments. Flow cytometry analysis for this project was in part done on instruments in the Stanford Shared FACS Facility.
This work was supported by NIH grants U01HL107388 and the Diabetes Genomics and Analysis Core of the Stanford Diabetes Research Center (P30DK116074).
Ethics approval and consent to participate
Study is approved by Stanford Institutional Review Board.
Consent for publication
All subjects gave written informed consent for study participation and publication.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. RNA sequencing Data QC, related to results section regarding RNA sequencing. (A) Average read quality of each sequencing data. The x-axis shows the average read Q-score and the y-axis indicates the number of reads. (B) Base quality along the reads of the sequencing data. The x-axis shows the position in the read is plotted and the y-axis shows the Q-score is plotted. Median value Q-score are shown by the red lines. The mean value Q-score are represented by dark blue line. The inter-quartile range is shown by the boxplot, while the 10% and 90% points are represented by the whiskers. (TIFF 6593 kb)
Figure S2. Mapping results, related to the “Results” section regarding RNA sequencing. Samples’ summary of mapping results of the reads. In each sample reads could be classified into: mapped (aligned to reference genome), high abundance or out-mapped (e.g., rRNA, mtRNA, polyA, polyC) and unmapped (reads which did not align to anything). (TIFF 6593 kb)
Figure S3. Differentially expressed genes in selected GO pathways related to 24 h of hiPSC culture with or without insulin, related to Table 1 regarding the top 25 significant GO (biological process) terms. (TIFF 6593 kb)
Table S1. List of hiPSC lines used. Table S2. Number of identified genes in each sample, related to results section regarding RNA sequencing. Table S3. The 20 most significantly differentially expressed genes at 24 h of hiPSC culture with or without insulin, related to results section regarding RNA sequencing. The list is ordered by adjusted p values. Average read counts of samples with no insulin treatment are shown in the “baseMean” column. Fold change is the log2 fold change of normalized read count per gene between groups. Followed are the lfcSE (log fold change Standard Error) and stat (Wald statistics). Also listed are p values and adjusted p values corrected for multiple testing using the Benjamini-Hochberg False Discovery Rate (FDR) approach. Table S4. The 20 most significantly differentially expressed genes at 72 h of hiPSC culture with or without insulin, related to results section regarding RNA sequencing. The list is ordered by adjusted p values. Average read counts of samples with no insulin treatment are shown in the “baseMean” column. Fold change is the log2 fold change of normalized read count per gene between groups. Followed are the lfcSE (log fold change Standard Error) and stat (Wald statistics). Also listed are p values and adjusted p values corrected for multiple testing using the Benjamini-Hochberg False Discovery Rate (FDR) approach. Table S5. The top 100 significant GO (biological process) terms associated with differentially expressed transcripts after 72 h of hiPSC culture with or without insulin, related to results sections regarding RNA sequencing. Table S6. Differentially expressed genes from culture of hiPSCs with or without insulin that have been identified as a part of KEGG signaling pathways regulating pluripotency of stem cells (human). Table S7. Differentially expressed genes from culture of hiPSCs with or without insulin that have been identified as a part of KEGG insulin signaling pathway (human). (DOCX 74 kb)