- Open Access
Using low-risk factors to generate non-integrated human induced pluripotent stem cells from urine-derived cells
Stem Cell Research & Therapy volume 8, Article number: 245 (2017)
Because the lack of an induced pluripotent stem cell (iPSC) induction system with optimal safety and efficiency limits the application of these cells, development of such a system is important.
To create such an induction system, we screened a variety of reprogrammed plasmid combinations and multiple compounds and then verified the system’s feasibility using urine cells from different individuals. We also compared large-scale iPSC chromosomal variations and expression of genes associated with genomic stability between this system and the traditional episomal system using karyotype and quantitative reverse transcription polymerase chain reaction analyses.
We developed a high-efficiency episomal system, the 6F/BM1-4C system, lacking tumorigenic factors for human urine-derived cell (hUC) reprogramming. This system includes six low-risk factors (6F), Oct4, Glis1, Klf4, Sox2, L-Myc, and the miR-302 cluster. Transfected hUCs were treated with four compounds (4C), inhibitor of lysine-demethylase1, methyl ethyl ketone, glycogen synthase kinase 3 beta, and histone deacetylase, within a short time period. Comparative analysis revealed significantly decreased chromosomal variation in iPSCs and significantly increased Sirt1 expression compared with iPSCs induced using the traditional episomal system.
The 6F/BM1-4C system effectively induces reprogramming of urine cells in samples obtained from different individuals. iPSCs induced using the 6F/BM1-4C system are more stable at the cytogenetic level and have potential value for clinical application.
Advancements in induced pluripotent stem cell (iPSC) technology have provided great opportunities for regenerative medicine and tumor immunotherapy [1,2,3]. However, human iPSCs (hiPSCs) are primarily induced using retroviral or lentiviral vectors carrying reprogramming factors [4, 5], and exogenous DNA fragments can randomly insert into genomic DNA and induce cell transformation, thereby preventing the clinical application of iPSCs. To date, many non-integrating methods have been generated, including mRNA and protein transfection [6, 7], Sendai virus (SeV) , piggyback (PB) transposons , and episomal vectors . Nonetheless, mRNA and protein transfections are associated with high preparation costs and low induction efficiency, and a risk of transformation is associated with the retention of SeV RNA in the first passages of iPSC lines [8, 11] and PB transposons. Although episomal induction systems can avoid these problems and are widely used for reprogramming, most of these systems utilize at least one tumorigenic factor such as c-Myc, SV40-LT, and p53 inhibitors, including p53 RNA interference (RNAi) or small molecule inhibitors [10,13,14,15,16,, 12–17]. As induction efficiency varies when the same method is used to reprogram different types of somatic cells or when different methods are applied to reprogram the same type of somatic cells , it is crucial to optimize the induction method for each type of cell. Human urine-derived cells (hUCs) are ideal donors for iPSC generation: their isolation is simple and nontraumatic. In addition, these cells are easy to expand in vitro and can be used as the main source for iPSCs; thus, their use is cost-effective and universal [18,19,20]. In the present study, an episomal vector was used for hUC reprogramming [14, 17].
The nontransformative MYC family protein L-Myc can be replaced with c-Myc to induce iPSCs [15, 21]. Glis1, which is enriched in oocytes, can also replace c-Myc in the classical OSKM (Oct4, Sox2, Klf4, and c-Myc) system; chimeric mice generated from iPSCs induced using OSK and Glis1 have longer survival times than those generated from iPSCs induced by OSKM . Moreover, there is a positive correlation between chimeric mouse mortality and mouse tumor mortality , suggesting that Glis1 is a safety factor for iPSC generation. The miR-302 family, which is specifically expressed in embryonic stem cells (ESCs), can partially or completely replace reprogramming factors and increase reprogramming efficiency [14, 23, 24]. Furthermore, the miR-302 family activates Ink4a and Arf to suppress the tumorigenesis of human pluripotency stem cells by targeting the oncogene Bmi1 , and Arf/p53 pathway activation suppresses somatic cell reprogramming [15, 16]. Therefore, miR-302 s are typically important factors that promote somatic cell reprogramming, but targeted factors that inhibit reprogramming exist in some signal pathways. Several studies to date have suggested that long noncoding RNAs (lncRNAs) regulate development and tumorigenesis; for example, long intergenic noncoding RNA, regulator of reprogramming (lincRNA-ROR) regulates the self-renewal and pluripotency of human ESCs (hESCs) and the reprogramming of hiPSCs [26, 27]. In this study, we applied hUCs as donor cells to induce iPSCs using low-risk factors, and then we screened a combination of low-risk reprogramming factors, including Oct4, Glis1, Klf4, Sox2, L-Myc, and the miR-302 cluster.
To improve non-integrated reprogramming efficiency, we optimized our culture system and observed iPSC induction with high efficiency when four compounds (Parnate, PD0325901, CHIR99021, and sodium butyrate) were added to the medium for no more than 4 days. To analyze iPSC safety, a karyotype analysis was performed, and the results showed that significantly lower iPSC chromosomal variation was induced when using this system than when using episomal systems containing SV40-LT and c-Myc.
hUCs were collected according to methods reported previously [14, 18]. Briefly, 100–1000 ml of urine was collected from donors, centrifuged at 1010 × g for 5 minutes, and washed with phosphate-buffered saline (PBS). The cells were maintained in 24-well plates coated with 0.1% gelatin (ES-006-B; Millipore, Germany) in RM1 medium (50% Renal Epithelial Cell Growth Medium (REGM) (CC-3190; Lonza, USA) and 44% Dulbecco's Modified Eagle Medium (DMEM) (SH30022; HyClone, USA) supplemented with 5% fetal bovine serum (FBS) (P30-3302; PAN Biotech, Germany), 0.5% nonessential amino acids (NEAA) (11140050; Gibco, USA), 0.5% GlutaMax (35050-061; Gibco, USA)) and 1 × Primocin (ant-pm-2; InvivoGen, USA); 0.25% trypsin-EDTA (25200072; Gibco, USA) was used for dissociation of primary hUCs. RM1 or RM2 (82% DMEM (SH30022; HyClone, USA) supplemented with 5% FBS, 1% human keratinocyte growth supplement (HKGS) (S-001-5; Gibco, USA), 1% NEAA, and 1% GlutaMax) was used for hUC culture.
The HN4 hESC line was obtained from the Chinese Academy of Sciences, and both HN4 and hiPSCs were maintained in the hESC medium BioCISO (BC-PM0001; BIOCARE Biotech, China) in plates coated with Matrigel (354277; Corning, USA).
pCEP4 (V04450; Invitrogen, USA) was digested using the restriction enzymes NruI and SalI and ligated with synthesized multiple cloning site (MCS) oligonucleotides to obtain the plasmid pE2.1. The EF1α promoter (NruI, NheI), BGH-PA element (BamHI, PmeI), EF1α promoter (PmeI, BglII), and BGH-PA element (PacI, SalI) were cloned into pE2.1 to obtain the plasmid pE3.1. The process chart for pE3.1 plasmid construction is shown in Additional file 1: Figure S1a. To obtain the plasmid pE3.2, the EF1α promoter (Pme I, Bgl II) was replaced with the cytomegalovirus (CMV) promoter (PmeI, BglII). The sequences of Oct4, Sox2, and Klf4 were subcloned from OKSIM (Plasmid 24603; Addgene, USA). The Glis1 sequence was subcloned from pMXs-Glis1 (Plasmid 30166; Addgene, USA) and the L-Myc sequence from pMXs-Hu-L-Myc (Plasmid 30166; Addgene, USA). The Oct4-P2A-Glis1, Klf4-P2A-Sox2, and Oct4-P2A-L-Myc sequences were obtained using overlap PCR. The sequence of the miR-302 cluster was cloned from genomic DNA. The exon sequences of lincRNA-ROR were amplified using genomic DNA and synthesized DNA, and overlap PCR was used to obtain the complete lincRNA-ROR sequence. These DNA sequences were cloned into the plasmids pE3.1 and pE3.2, respectively, to generate the plasmids pE3.1-OL--KS, pE3.1-OG--KS, pE3.1-Oct4--Klf4, pE3.1-Glis1--LINC-ROR, pE3.1-L-Myc--hmiR-302 cluster, pE3.1-Oct4--Sox2, and pE3.2-L-Myc--hmiR-302 cluster. Information regarding the factors, primer sequences, and MCS is shown in Additional file 2: Table S1.
For hUC16 reprogramming, 1.5 × 106 hUC16 cells were transfected with plasmids (Additional file 3: Table S2) using the T-020 program of a Lonza Nucleofector 2b Device and a Basic Epithelial Cells Nucleofector Kit (VPI-1005; Lonza, USA). Transfected hUC16 cells were seeded into six-well plates coated with 0.1% gelatin and cultured using RM2 medium. After 24 h, the cells were dissociated using 0.25% trypsin-EDTA (25200-072; Gibco, USA), and 2 × 104 cells were seeded into 12-well plates coated with Matrigel. To induce iPSCs, the medium was changed to BioCISO-BM1 medium (BC-BM001; BIOCARE Biotech, China) containing 4i (A83-01 (0.5 μM, BC-SMC-A01-10; BIOCARE Biotech, China), Thiazovivin (0.5 μM, BC-SMC-T01-10; BIOCARE Biotech, China), CHIR99021 (3 μM, BC-SMC-C01-10; BIOCARE Biotech, China), and PD03254901 (0.5 μM, BC-SMC-P01-10; BIOCARE Biotech, China)) after 24 h. The medium was then changed to BioCISO on day 15. Alkaline phosphatase (AP) staining was performed on day 18, and the induction efficiency was calculated according to the formula:
For reprogramming using the 6F/BM1-4C system, approximately 2.8 × 106–3.5 × 106 hUCs were transfected with 4.0 μg pE3.1-OG--KS and 2.8 μg pE3.1-L-Myc--hmiR-302 cluster using the same nucleofector method. The transfected hUCs were placed in plates coated with Matrigel and cultured with RM1 medium. On day 3 after nucleofector addition, the medium was changed to BioCISO-BM1 medium containing 2 μM Parnate (also known as tranylcypromine hydrochloride, 1986-47-6; Curegenix, China). The medium was then changed to BioCISO-BM1 medium containing 2 μM Parnate, 0.25 mM sodium butyrate (NaB) (303410-100G; Sigma, USA), 3 μM CHIR99021, and 0.5 μM PD03254901 on day 5, to BioCISO-BM1 on day 7, and to BioCISO on day 17. iPSC colonies were collected or stained with AP on day 19. The induction efficiency was calculated according to the formula:
Compounds used in the present study also included dimethyloxaloylglycine (DMOG) (0.1 μM, D1070; Frontier Scientific, USA), PS48 (5 μM, 1180676-32-7; Curegenix, China), SC-79 (0.5 μM, 4635; Tocris, USA), forskolin (5 μM, 66575-29-9; Curegenix, China), and 3-deazaneplanocin A (DZNEP) (0.05 μM, 4703; Tocris, USA).
For reprogramming using the 4F2L-6C system, 3.0 × 106 hUCs were transfected with 4.0 μg pEP4-E02S-ET2K and 2.8 μg pCEP4-M2L using the same nucleofector method. The transfected hUCs were placed in plates coated with Matrigel and cultured with RM1 medium. The medium was changed to BioCISO-BM1 medium containing 2 μM Parnate on day 3 after nucleofector addition, and then changed to BioCISO-BM1 medium containing 2 μM Parnate, 0.25 mM NaB, 3 μM CHIR99021, 0.5 μM PD03254901, 0.5 μM A83-01, and 0.5 μM Thiazovivin on day 5. The medium was changed to BioCISO-BM1 on day 7 and to BioCISO on day 17. iPSC colonies were collected on day 34.
AP staining, non-integrated PCR analysis, flow cytometry analysis, immunofluorescence analysis, bisulfate sequencing, in-vitro embryoid body (EB) differentiation assays, and in-vivo teratoma formation were conducted according to previous methods [14, 17]. Briefly, AP staining was carried out using nitroblue tetrazolium (NBT) (N104908-1 g; Aladdin, China) and 5-bromo-4-chloro-3-indolyl phosphate (BCIP) (BIMB1018; J&K Chemical Technology, China). PCR was applied to analyze the integration of exogenous reprogramming factors and the episomal backbone with the primers presented in Additional file 2: Table S1. Flow cytometry and immunofluorescence analyses were employed to examine human pluripotency markers with the following antibodies: anti-Oct4 (130-105-606; Miltenyi Biotec, Germany), anti-SSEA4 (sc-21704; Santa Cruz, USA), anti-Tra-1-60 (sc-21705; Santa Cruz, USA), anti-Tra-1-81 (sc-21706; Santa Cruz, USA), anti-IgM-PE (sc-3768; Santa Cruz, USA), and anti-IgG3-PE (sc-3767; Santa Cruz, USA). Bisulfate sequencing was used to determine Oct4 and Nanog promoter methylation with the primers presented in Additional file 2: Table S1. The PCR product was cloned into the pMD18-T vector and subsequently sequenced. For the in-vitro EB differentiation assay, cells were scraped from plates after dissociation using BioC-PDE1 (BC-PDE1; BIOCARE Biotech, China) and cultured in six-well suspension culture plates (657185; Greiner, Germany) with BioCISO-EB1 medium (BC-EB001; BIOCARE Biotech, China) for 7 days to obtain EBs. The EBs were then cultured in six-well culture plates (657160; Greiner, Germany) coated with Matrigel for 7–14 days. The cells were collected, and expression profiles of marker genes in the three germ layers were determined using quantitative reverse transcription PCR (qRT-PCR) with primers obtained from BIOCARE Biotech (BioEB-primer). For in-vivo teratoma formation, iPSCs were cultured to approximately 85% confluence; after 10–15 minutes of dissociation using BioC-PDE1, cells were scraped from the plates. The hind-limb muscle and forelimb subcutaneous muscle of 6-week-old NOD/SCID mice were injected with 130 μl BioCISO culture medium, 70 μl Matrigel, and cell suspensions. The formation of teratomas could be observed after 6–8 weeks; when the teratomas reached a certain size, they were removed and fixed with 4% paraformaldehyde. The tissues were embedded with paraffin, sectioned, stained with hematoxylin and eosin, and analyzed under a microscope. The procedures were performed according to IACUC (Institutional Animal Care and Use Committee; YS-YFStudy060-20160315).
Sample preparation for karyotyping was conducted as described previously . Cells were treated with 50 ng/ml colchicine (Xy008; Xiangya Gene Technology, China) for 16 h, and the Ikaros karyotyping system was used to analyze karyotypes. The aneuploid evaluation is shown in Additional file 4: Table S3.
Quantitative reverse transcription polymerase chain reaction
Total RNA was isolated using RNAiso Plus (TaKaRa), and M-MLV Reverse Transcriptase (TaKaRa) was used to synthetize cDNA. Specific stem-loop primers and random primers were used for reverse transcription of microRNAs and mRNAs into cDNA, respectively. mRNA and miRNA expression levels were determined using SYBR Premix Ex Taq™ (TaKaRa). Reactions were performed in triplicate using a LightCycler 480II/96 system (Roche, Switzerland). mRNA expression was normalized to GAPDH, and microRNA (miRNA) expression was normalized to U6 small nuclear RNA (snRNA). The primers are presented in Additional file 2: Table S1.
Western blot analysis
Radioimmunoprecipitation assay (RIPA) buffer (CW2333; Cwbiotech, China) supplemented with a protease inhibitor cocktail (PI003; BOCAI Technology, China) and phenylmethylsulfonyl fluoride (PMSF; Dingguo Changshen Biotech, China) was used to isolate cellular proteins. Equivalent amounts of protein were separated by sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to polyvinylidene difluoride (PVDF) membranes. The membranes were incubated with specific primary antibodies against Oct4 (2840; Cell Signaling Technology, USA), Glis1 (SAB2700289; Sigma, USA), Klf4 (ab72543; Abcam, UK), Sox2 (3579; Cell Signaling Technology, USA), L-Myc (sc-790; Santa Cruz Biotechnology, USA), and GAPDH (KC-5G4; KangChen Biotech, China), followed by horseradish peroxidase-conjugated secondary antibodies: goat anti-Rabbit IgG (ZB-2301; ZsBio, China) and anti-mouse IgG-HRP (IH-0031; Dingguo Changshen Biotech, China). Bands were visualized using enhanced chemiluminescence (ECL) (34087; Thermo, USA).
GeneChip Human Transcriptome Array 2.0 (Affymetrix HTA 2.0, USA) was utilized to determine the gene expression profiles of human ESCs, iPSCs, and hUCs. The experiments were conducted according to the manufacturer’s instructions.
SPSS 18.0 was used to perform statistical analysis. The results are presented as the mean ± standard deviation (SD) of at least three repeated individual experiments for each group. Statistical differences were examined using Student’s t test. For analysis of the chromosome abnormality rate, a four-table chi-square test was applied. P < 0.05 was considered statistically significant.
Microarray data for human ESs, iPSCs, and hUCs have been submitted to Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) under accession number GSE85885.
Screening low-risk reprogramming factors using hUCs
Since Yamanaka used OSKM to induce reprogramming, many genes and non-RNAs that improve reprogramming efficiency have been reported. To induce reprogramming in this study, we employed low-risk factors, including L-Myc, Glis1, lincRNA-ROR, and the miR-302 cluster, and randomly combined them with Oct4, Sox2, and Klf4 in the Epstein–Barr virus-encoded nuclear antigen-1 (EBNA)-oriP episomal vector (Additional file 1: Figure S1a, S1b, S1c). The miR-302 family can increase reprogramming efficiency by replacing reprogramming factors [14, 23, 24]. Furthermore, miR-302 family members target the oncogene Bmi1 and suppress tumorigenesis, which can inhibit somatic cell reprogramming [15, 16, 25]. Therefore, miR-302 s have a positive impact on somatic cell reprogramming, but in some pathways these family members can inhibit reprogramming by indirectly activating targets that inhibit reprogramming. Because the expression levels of the miR-302 cluster must be precisely regulated, we used different promoters to exogenously express the miR-302 cluster and screened for optimal expression (Additional file 1: Figure S1a, S1b).
To evaluate the best plasmid combination for reprogramming, we used hUC16 cells constructed in our laboratory that showed high proliferation (Additional file 5: Figure S2a) and transfected these cells with different plasmid combinations (Additional file 5: Figure S2b), followed by AP staining to identify iPSCs after 18 days of nucleofection. Three groups of cells harboring the reprogramming factors Oct4, Glis1, Klf4, Sox2, L-Myc, lincRNA-ROR, and the miR-302 cluster with high AP-positive scores were selected for further analysis (Additional file 5: Figure S2c, S2d). The initial screen was performed in cells (UC16) with high proliferative ability and strong anti-stress capacity, and the best three combinations were further tested using other hUCs to ensure reliability. To induce reprogramming, we transfected the three plasmid combinations into hUCs and cultured in hESC basal medium (BioCISO-BM1) with the addition of four small inhibitors that have been widely used for reprogramming [12, 29]: 4i, the TGF-β/Activin/Nodal receptor inhibitor A83-01 (0.5 μM); the MEK inhibitor PD0325901 (0.5 μM); the GSK3β inhibitor CHIR99021 (3 μM); and the ROCK inhibitor Thiazovivin (0.5 μM) (Fig. 1a, b). As revealed by AP staining, the combination termed 6F, which includes Oct4, Glis1, Klf4, Sox2, L-Myc, and the miR-302 cluster initially expressed from the CMV promoter (Fig. 1c, d), showed the highest reprogramming efficiency at 19 days post nucleofection. Thus, we selected 6F, which does not contain high-risk tumorigenic factors such as c-Myc, SV40-LT, and p53 inhibitors, as the reprogramming induction combination for hUCs and found that it successfully induced hUCs into iPSCs.
Optimization of compounds in the 6F combination system
Different cell lineages exhibit different gene expression profiles, and ideal reprogramming factor combinations and induction conditions depend on the cell type [12, 30, 31]. To determine whether 4i is the best induction condition for the 6F combination, we first examined the effects of compounds that are reported to regulate reprogramming in our system, including many compounds involved in signaling pathways. Some inhibitors or activators were observed to be unsuitable for our system; for example, the TGF-β/Activin/Nodal receptor inhibitor A83-01 and the ROCK inhibitor Thiazovivin did not promote programming (Additional file 6: Figure S3a–S3f). In addition, other compounds [31, 32], such as the PDK1 activator PS48, adenylyl cyclase activator forskolin, histone methyltransferase inhibitor DZNEP, HIF-α prolyl hydroxylase inhibitor DMOG, and Akt activator SC-79, did not promote or inhibit programming (Additional file 6: Figure S3g–S3l). Next, we examined the effects of combinations of compounds that promote reprogramming on the reprogramming efficiency. These compounds included the MEK inhibitor PD0325901 (PD, 0.5 μM), the GSK3β inhibitor CHIR99021 (CHIR, 3 μM), the histone deacetylase inhibitor NaB (0.25 mM), and the lysine-specific demethylase 1 inhibitor Parnate (Par, 2 μM). We observed that when Par was used in 6F combination reprogramming for an extended time interval, many cells in the cell layer shrank and were floating; however, the reprogramming efficiency increased when Par was added for a short time period (Additional file 6: Figure S3m–S3q). Thus, compound combinations, including Par, were used for no longer than 4 days. AP staining showed the highest reprogramming efficiency and the best iPSC quality when Par was added to the BioCISO-BM1 medium on days 3–4 and when PD, CHIR, NaB, and Par, termed 4C, were added on days 5–6 (Fig. 2a–c). Therefore, we selected BM1-4C as the optimum induction culture condition for the 6F combination induction system and named it 6F/BM1-4C.
Reprogramming hUCs from different human sources using the 6F/BM1-4C system
When hUCs were isolated from different individuals or from the same individual at different time points, the cells exhibited different morphologies when cultured, suggesting that these cells were of different types, which is consistent with previous reports [17,34,, 33–35]. To confirm that the 6F/BM1-4C system is suitable for reprogramming different types of hUCs, we reprogrammed seven groups of hUCs isolated from different individuals (Fig. 3a, b, Table 1). At approximately 19 days post nucleofection, we collected iPSCs for further purification and expansion culture (Fig. 3c). The remaining iPSCs were subjected to AP staining, which demonstrated that the 6F/BM1-4C system resulted in effective reprogramming of all hUCs into iPSCs. AP staining also revealed an AP-positive rate that varied between 0.00021 and 0.0741% for iPSCs reprogrammed from different types of hUCs (Fig. 3d, Table 1). Furthermore, PCR analysis demonstrated a lack of genomic integration of the exogenous gene sequence in 13 of 14 iPSC colonies (Fig. 3e, Additional file 7: Figure S4a). Moreover, karyotype analysis revealed that 13 iPSC colonies had normal chromosome numbers and G-band distributions (Fig. 4a, Additional file 7: Figure S4b). Flow cytometry expression and immunofluorescence analyses revealed that iPSCs induced using the 6F/BM1-4C system express the hESC-specific markers Oct4, SSEA4, Tra-1-60, and Tra-1-81, suggesting that 6F/BM1-4C-iPSCs possess the molecular characteristics of hESCs (Fig. 4b, c, Additional file 8: Figure S5a, S5b). In addition, bisulfite sequencing PCR analysis indicated that the endogenous pluripotency genes Oct4 and Nanog were activated and that their promoters were demethylated in 6F/BM1-4C-iPSCs, similar to what occurs in hESCs (Fig. 4d, Additional file 8: Figure S5c). The generated iPSCs were also able to differentiate into derivatives of all three germ layers, as determined using an in-vitro EB differentiation assay and an in-vivo teratoma formation assay (Fig. 4e, f, Additional file 8: Figure S5d, S5e). The gene expression profile of iPSCs was similar to that of hESCs and differed from that of hUCs, which was determined using Affymetrix gene microarray HTA 2.0. (Fig. 4g, Additional file 8: Figure S5f). The induction efficiencies of the different sources of hUCs as well as the iPSC characteristics are presented in Table 1. Together, these results indicate that the 6F/BM1-4C system has high reliability and versatility for reprogramming hUCs into iPSCs.
The 6F/BM1-4C system is safer than traditional episomal induction systems
We observed that before day 5 following nucleofection, the observed cell masses had high nuclear–cytoplasmic ratios, a phenotype that was similar to the lentiviral preprogramming process. From days 9 to 11, the colony mass underwent massive death, and the surviving cells grew slowly to maturity. However, because early induction using the traditional episomal induction system (termed 4F2L-6C), which includes Oct4, Sox2, Klf4, c-Myc, SV40-LT, and Lin28 , is a gradual process, the cells aggregated slowly and displayed little massive death (Fig. 5a, Additional file 9: Figure S6a, S6b); these results suggest that the two induction systems differ regarding iPSC generation. Therefore, we examined the safety of iPSCs generated using the 6F/BM1-4C system.
Evaluation of the application of iPSCs is based on safety. Schlaeger et al.  analyzed chromosomal variation among iPSCs induced using retrovirus, mRNA transfection, SeV, episomal, and lentivirus systems and observed large differences using these various methods. For example, the aneuploidy rate of iPSCs induced using the mRNA transfection method was only 2.3%, whereas that of iPSCs induced episomally was as high as 11.5%. However, these methods employ different reprogramming genes, and it is difficult to determine whether differences in chromosomal variation reflect the genes or the methods used for induction. In our study, we utilized an episomal vector to reprogram the same batch of hUCs isolated from the same donor using the 6F/BM1-4C and 4F2L-6C systems, and at least 60 iPSC colonies (65 for 6F/BM1-4C and 64 for 4F2L-6C) were selected for karyotype analysis (the specific criteria are summarized in Methods). We observed a significantly lower chromosomal abnormality rate for 6F/BM1-4C-iPSCs than for 4F2L-6C-iPSCs (P = 0.017, χ2 test; Fig. 5b, Table 2). We also determined the expression profiles of genes associated with genomic stability, such as Sirt1, p53, and CHK1, and found significantly high Sirt1 expression in iPSCs induced using the 6F/BM1-4C system (Fig. 5c). These data show that iPSCs induced using the 6F/BM1-4C system are safer than those induced using the traditional episomal induction method, which would be beneficial for clinical applications.
Reflecting their high efficiency and controllable cost, episomal plasmid-carried reprogramming factors are the most widely used approach for obtaining non-integrated iPSCs. Most episomal induction methods employ at least one oncogene or tumorigenic molecule, such as c-Myc, SV40-LT, p53 short hairpin RNA (shRNA), and the p53 small molecule inhibitor cyclin pifithrin-α [10,13,14,15,16,, 12–17], and iPSCs induced using previous methods cannot be used in clinical applications. Furthermore, induction methods that do not include tumorigenic factors are essential. In the present study, we constructed a low-risk 6F/BM1-4C reprogramming system, in which we eliminated the tumorigenic factors used in traditional episomal reprogramming systems, such as c-Myc, SV40-LT, and p53 inhibitor, and included Oct4, Glis1, Klf4, Sox2, L-Myc, the miR-302 cluster and four compounds, and then treated cells for no longer than 48 h to efficiently generate iPSCs from hUCs. This system also successfully converted hUCs from different sources into iPSCs and showed good reproducibility. Analyzing a large number of iPSCs by karyotype analysis, the 6F/BM1-4C-hiPSCs we generated exhibited fewer chromosome abnormalities compared with traditional 4F2L-6C-hiPSCs. In addition, expression of Sirt1, the NAD-dependent deacetylase necessary for maintaining iPSC genomic stability , in 6F/BM1-4C-iPSCs was high compared with iPSCs induced using the 4F2L-6C system, suggesting that 6F/BM1-4C-iPSC chromosome inheritance is more stable. Moreover, the presented method has a low cost, and the use of episomal plasmids makes this system suitable for clinical non-integrated iPSC preparation.
To obtain large-scale amounts of clinical-grade iPSCs, a reprogramming method with good reproducibility, non-tumorigenic reprogramming factors, and cost-effectiveness is needed; xeno-free components and a medium for primary somatic cell isolation to iPSC generation are also necessary. Besides, although the 6F/BM1-4C reprogramming system has relatively high reprogramming reproducibility, due to the heterogeneity of hUCs [18, 33, 34], it is difficult to accurately and separately perform multiplication culture in vitro to obtain a variety of different types of cells that meet the required number of experiments; hUCs from different donors and different batches of cells also show a wide range of induction efficiencies in the 6F/BM1-4C reprogramming system (Table 1). Therefore, except for a xeno-free induction reprogramming system, in the future the best reprogramming system should be screened for different types of hUCs or general suitability for a variety hUCs; in addition, the particular type of hUCs that is more easily reprogrammed or the particular type of hUCs that is more suitable for the 6F/BM1-4C system should be screened to find a more specific cultivation method for a particular type of hUCs, so that we can obtain a high-efficiency reprogramming system for screening high-quality clinical-grade iPSCs from a large number of iPSCs. Furthermore, when we used the xeno-free hESC E8 medium  to induce hUC reprogramming based on episomal vectors, we found it to be unsuitable after the addition of certain compounds, with deformed cells all dying (data not shown). Xeno-free extracellular matrices such as Vitronectin exhibit poor maintenance of iPSC self-renewal capacity . Conversely, Laminin521 maintains iPSC self-renewal capacity, but it is extremely expensive  and thus is not suitable for large-scale production of clinical-grade iPSCs. Accordingly, the selection of appropriate cell culture materials remains essential for further clinical applications.
We developed a safe method based on an episomal vector for inducing iPSCs from hUCs. This method does not involve the use of tumorigenic factors, such as c-Myc, SV40-LT, and p53 inhibitor. Karyotype analysis revealed that the chromosomal variation that occurred during iPSC generation in the present study was significantly low compared with the traditional method. Such low variability is critical for clinical applications of iPSCs.
Induced pluripotent stem cell
Human urine-derived cell
Human embryonic stem cell
Multiple cloning site
PD, CHIR, NaB, and Par
episomal-induced system containing six reprogramming factors (Oct4, Sox2, Klf4, c-Myc, SV40-LT, and Lin28) and six compounds (PD, CHIR, NaB, Par, thi, and A83-01)
PD, CHIR, NaB, Par, thi, and A83-01
episomal-induced system containing six reprogramming factors (Oct4, Glis1, Klf4, Sox2, L-Myc, and miR-302 cluster) and four compounds (PD, CHIR, NaB, and Par)
Oct4, Sox2, Klf4, and c-Myc
Oct4, Glis1, Klf4, Sox2, L-Myc, and the miR-302 cluster
Oct4, Sox2, Klf4, c-Myc, SV40-LT, and Lin28
Ando M, Nishimura T, Yamazaki S, Yamaguchi T, Kawana-Tachikawa A, Hayama T, Nakauchi Y, Ando J, Ota Y, Takahashi S, Nishimura K, Ohtaka M, Nakanishi M, Miles JJ, Burrows SR, Brenner MK, Nakauchi H. A safeguard system for induced pluripotent stem cell-derived rejuvenated T cell therapy. Stem Cell Rep. 2015;5:597–608.
Vizcardo R, Masuda K, Yamada D, Ikawa T, Shimizu K, Fujii S, Koseki H, Kawamoto H. Regeneration of human tumor antigen-specific T cells from iPSCs derived from mature CD8(+) T cells. Cell Stem Cell. 2013;12:31–6.
Srivastava D, DeWitt N. In vivo cellular reprogramming: the next generation. Cell. 2016;166:1386–96.
Yu J, Vodyanik MA, Smuga-Otto K, Antosiewicz-Bourget J, Frane JL, Tian S, Nie J, Jonsdottir GA, Ruotti V, Stewart R, Slukvin II, Thomson JA. Induced pluripotent stem cell lines derived from human somatic cells. Science. 2007;318:1917–20.
Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, Yamanaka S. Induction of pluripotent stem cells from adult human fibroblasts by defined factors. Cell. 2007;131:861–72.
Warren L, Manos PD, Ahfeldt T, Loh YH, Li H, Lau F, Ebina W, Mandal PK, Smith ZD, Meissner A, Daley GQ, Brack AS, Collins JJ, Cowan C, Schlaeger TM, Rossi DJ. Highly efficient reprogramming to pluripotency and directed differentiation of human cells with synthetic modified mRNA. Cell Stem Cell. 2010;7:618–30.
Zhou H, Wu S, Joo JY, Zhu S, Han DW, Lin T, Trauger S, Bien G, Yao S, Zhu Y, Siuzdak G, Scholer HR, Duan L, Ding S. Generation of induced pluripotent stem cells using recombinant proteins. Cell Stem Cell. 2009;4:381–4.
Fusaki N, Ban H, Nishiyama A, Saeki K, Hasegawa M. Efficient induction of transgene-free human pluripotent stem cells using a vector based on Sendai virus, an RNA virus that does not integrate into the host genome. Proc Jpn Acad Ser B Phys Biol Sci. 2009;85:348–62.
Woltjen K, Michael IP, Mohseni P, Desai R, Mileikovsky M, Hamalainen R, Cowling R, Wang W, Liu P, Gertsenstein M, Kaji K, Sung HK, Nagy A. piggyBac transposition reprograms fibroblasts to induced pluripotent stem cells. Nature. 2009;458:766–70.
Yu J, Hu K, Smuga-Otto K, et al. Human induced pluripotent stem cells free of vector and transgene sequences. Science. 2009;324:797-801.
Okita K, Ichisaka T, Yamanaka S. Generation of germline-competent induced pluripotent stem cells. Nature. 2007;448:313–7.
Yu J, Chau KF, Vodyanik MA, Jiang J, Jiang Y. Efficient feeder-free episomal reprogramming with small molecules. PLoS One. 2011;6:e17557.
Valamehr B, Robinson M, Abujarour R, Rezner B, Vranceanu F, Le T, Medcalf A, Lee TT, Fitch M, Robbins D, Flynn P. Platform for induction and maintenance of transgene-free hiPSCs resembling ground state pluripotent stem cells. Stem Cell Reports. 2014;2:366–81.
Xue Y, Cai X, Wang L, Liao B, Zhang H, Shan Y, Chen Q, Zhou T, Li X, Hou J, Chen S, Luo R, Qin D, Pei D, Pan G. Generating a non-integrating human induced pluripotent stem cell bank from urine-derived cells. PLoS One. 2013;8:e70573.
Okita K, Matsumura Y, Sato Y, Okada A, Morizane A, Okamoto S, Hong H, Nakagawa M, Tanabe K, Tezuka K, Shibata T, Kunisada T, Takahashi M, Takahashi J, Saji H, Yamanaka S. A more efficient method to generate integration-free human iPS cells. Nat Methods. 2011;8:409–12.
Okita K, Yamakawa T, Matsumura Y, Sato Y, Amano N, Watanabe A, Goshima N, Yamanaka S. An efficient nonviral method to generate integration-free human-induced pluripotent stem cells from cord blood and peripheral blood cells. Stem Cells. 2013;31:458–66.
Li D, Wang L, Hou J, Shen Q, Chen Q, Wang X, Du J, Cai X, Shan Y, Zhang T, Zhou T, Shi X, Li Y, Zhang H, Pan G. Optimized approaches for generation of integration-free iPSCs from human urine-derived cells with small molecules and autologous feeder. Stem Cell Rep. 2016;6:717–28.
Zhou T, Benda C, Duzinger S, Huang Y, Li X, Li Y, Guo X, Cao G, Chen S, Hao L, Chan YC, Ng KM, Ho JC, Wieser M, Wu J, Redl H, Tse HF, Grillari J, Grillari-Voglauer R, Pei D, Esteban MA. Generation of induced pluripotent stem cells from urine. Clin J Am Soc Nephrol. 2011;22:1221–8.
Wang L, Wang L, Huang W, Su H, Xue Y, Su Z, Liao B, Wang H, Bao X, Qin D, He J, Wu W, So KF, Pan G, Pei D. Generation of integration-free neural progenitor cells from cells in human urine. Nat Methods. 2013;10:84–9.
Cheng L, Hu W, Qiu B, Zhao J, Yu Y, Guan W, Wang M, Yang W, Pei G. Generation of neural progenitor cells by chemical cocktails and hypoxia. Cell Res. 2014;24:665–79.
Nakagawa M, Takizawa N, Narita M, Ichisaka T, Yamanaka S. Promotion of direct reprogramming by transformation-deficient Myc. Proc Natl Acad Sci U S A. 2010;107:14152–7.
Maekawa M, Yamaguchi K, Nakamura T, Shibukawa R, Kodanaka I, Ichisaka T, Kawamura Y, Mochizuki H, Goshima N, Yamanaka S. Direct reprogramming of somatic cells is promoted by maternal transcription factor Glis1. Nature. 2011;474:225–9.
Lin SL, Chang DC, Chang-Lin S, Lin CH, Wu DT, Chen DT, Ying SY. Mir-302 reprograms human skin cancer cells into a pluripotent ES-cell-like state. RNA. 2008;14:2115–24.
Lin SL, Chang DC, Lin CH, Ying SY, Leu D, Wu DT. Regulation of somatic cell reprogramming through inducible mir-302 expression. Nucleic Acids Res. 2011;39:1054–65.
Lin SL, Chang DC, Ying SY, Leu D, Wu DT. MicroRNA miR-302 inhibits the tumorigenecity of human pluripotent stem cells by coordinate suppression of the CDK2 and CDK4/6 cell cycle pathways. Cancer Res. 2010;70:9473–82.
Yan X, Lv Y, Ma J, et al. The development of extrahepatic portacaval shunt device based on magnetic compression technique through the interventional procedure. Chin J Med Instrum. 2013;37:421-22, 36.
Loewer S, Cabili MN, Guttman M, Loh YH, Thomas K, Park IH, Garber M, Curran M, Onder T, Agarwal S, Manos PD, Datta S, Lander ES, Schlaeger TM, Daley GQ, Rinn JL. Large intergenic non-coding RNA-RoR modulates reprogramming of human induced pluripotent stem cells. Nat Genet. 2010;42:1113–7.
Meisner LF, Johnson JA. Protocols for cytogenetic studies of human embryonic stem cells. Methods. 2008;45:133–41.
Xu Y, Zhu X, Hahm HS, Wei W, Hao E, Hayek A, Ding S. Revealing a core signaling regulatory mechanism for pluripotent stem cell survival and self-renewal by small molecules. Proc Natl Acad Sci U S A. 2010;107:8129–34.
Li W, Zhou H, Abujarour R, Zhu S, Young Joo J, Lin T, Hao E, Scholer HR, Hayek A, Ding S. Generation of human-induced pluripotent stem cells in the absence of exogenous Sox2. Stem Cells. 2009;27:2992–3000.
Zhu S, Li W, Zhou H, Wei W, Ambasudhan R, Lin T, Kim J, Zhang K, Ding S. Reprogramming of human primary somatic cells by OCT4 and chemical compounds. Cell Stem Cell. 2010;7:651–5.
Hou P, Li Y, Zhang X, Liu C, Guan J, Li H, Zhao T, Ye J, Yang W, Liu K, Ge J, Xu J, Zhang Q, Zhao Y, Deng H. Pluripotent stem cells induced from mouse somatic cells by small-molecule compounds. Science. 2013;341:651–4.
Rahmoune H, Thompson PW, Ward JM, Smith CD, Hong G, Brown J. Glucose transporters in human renal proximal tubular cells isolated from the urine of patients with non-insulin-dependent diabetes. Diabetes. 2005;54:3427–34.
Dorrenhaus A, Muller JI, Golka K, Jedrusik P, Schulze H, Follmann W. Cultures of exfoliated epithelial cells from different locations of the human urinary tract and the renal tubular system. Arch Toxicol. 2000;74:618–26.
Zhou T, Benda C, Dunzinger S, Huang Y, Ho JC, Yang J, Wang Y, Zhang Y, Zhuang Q, Li Y, Bao X, Tse HF, Grillari J, Grillari-Voglauer R, Pei D, Esteban MA. Generation of human induced pluripotent stem cells from urine samples. Nat Protoc. 2012;7:2080–9.
Schlaeger TM, Daheron L, Brickler TR, Entwisle S, Chan K, Cianci A, DeVine A, Ettenger A, Fitzgerald K, Godfrey M, Gupta D, McPherson J, Malwadkar P, Gupta M, Bell B, Doi A, Jung N, Li X, Lynes MS, Brookes E, Cherry AB, Demirbas D, Tsankov AM, Zon LI, Rubin LL, Feinberg AP, Meissner A, Cowan CA, Daley GQ. A comparison of non-integrating reprogramming methods. Nat Biotechnol. 2015;33:58–63.
De Bonis ML, Ortega S, Blasco MA. SIRT1 is necessary for proficient telomere elongation and genomic stability of induced pluripotent stem cells. Stem Cell Rep. 2014;2:690–706.
Beers J, Gulbranson DR, George N, Siniscalchi LI, Jones J, Thomson JA, Chen G. Passaging and colony expansion of human pluripotent stem cells by enzyme-free dissociation in chemically defined culture conditions. Nat Protoc. 2012;7:2029–40.
Rodin S, Antonsson L, Hovatta O, Tryggvason K. Monolayer culturing and cloning of human pluripotent stem cells on laminin-521-based matrices under xeno-free and chemically defined conditions. Nat Protoc. 2014;9:2354–68.
The authors thank American Journal Experts for English editing.
Availability of data and materials
Supporting data are available from the corresponding author upon reasonable request.
Ethics approval and consent to participate
hUCs were obtained from subjects who provided written informed consent for the use of their cells for iPSC generation, and the experiments involving human subjects and animal research were reviewed and approved by the Institutional Review Board at Guangzhou Biocare Institute of Cancer (BioCareIRB201304a#).
Consent for 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.
showing expression of factors from episomal vectors. a pE3.1 plasmid construction process chart (upper). Schematic representation of seven constructed episomal vectors. pEF1α EF1α promoter, pCMV CMV promoter (below). b Quantitative real-time PCR assay for Oct4, Glis1, Klf4, Sox2, L-Myc, linc-RoR, miR-367, miR-302a, miR-302b, miR-302c, and miR-302d. c Western blot assay for Oct4, Glis1, Klf4, Sox2, and L-Myc carried on episomal vectors. GAPDH was used as the loading control. (PDF 171 kb)
Presenting information for primers or functional fragments used in the present study. (XLSX 17 kb)
presenting plasmid combinations used for screening low-risk factors. (XLSX 9 kb)
presenting the method of karyotype analysis. (XLSX 10 kb)
showing use of hUC16 cells to screen for low-risk factors for iPSC generation. a hUC16 morphology. b Strategy to screen low-risk factor combinations using hUC16 cells. c AP staining for iPSC generation using different factor combinations. d Numbers of AP-positive colonies. (PDF 137 kb)
showing the effects of multiple compounds in the 6F combination system. a Strategy to optimize six-factor combinations using A83-01. b AP staining for iPSCs induced using A83-01. c Number of AP-positive colonies induced using A83-01. P(B) = 0.002. d Strategy to optimize six-factor combinations using Thiazovivin (thi). e AP staining of iPSCs induced using thi. f Number of AP-positive colonies induced using thi. P(D) = 0.000. g Strategy to optimize six-factor combinations using forskolin, PS48, and sc-79. h AP staining for iPSCs induced using forskolin, PS48, and sc-79. i Number of AP-positive colonies induced using forskolin, PS48, and sc-79. j Strategy to optimize six-factor combinations using DMOG and DZNEP. k AP staining for iPSCs induced using DMOG and DZNEP. l Number of AP-positive colonies induced using DMOG and DZNEP. P(C) = 0.041. m Strategy to optimize six-factor combination using Parnate in the early induction stage. n AP staining for iPSCs induced using Parnate in the early induction stage. Arrow indicates cell edge hemming. o Strategy to optimize six-factor combination treated with Parnate for a short time. p AP staining for iPSCs induced using Parnate for a short time. q Number of AP-positive colonies. P(C) = 0.04. Error bars indicate mean ± SD. *P < 0.05, **P < 0.01, ***P < 0.001. Scale bars, 100 μm. (PDF 227 kb)
showing non-integrating analysis and karyotype assays of iPSCs induced with the 6F/BM1-4C system. a Non-integrating analysis of genomic DNA in iPSCs. Representative lanes: 1, H2O; 2, pE3.1-OG--KS and pE3.2-L-Myc--hmiR-302 cluster; 3, UC5; 4, UC5, pE3.1-OG--KS, and pE3.2-L-Myc--hmiR-302 cluster; 5, UC5iPSC1; 6, UC5iPSC2; 7, UC6; 8, UC6, pE3.1-OG--KS, and pE3.2-L-Myc--hmiR-302 cluster; 9, UC6iPSC1; 10, UC6iPSC2; 11, UC7; 12, UC7, pE3.1-OG--KS, and pE3.2-L-Myc--hmiR-302 cluster; 13, UC7iPSC1; 14, UC7iPSC2. OriP in lane 9 exhibited integration. b Karyotype analysis of iPSCs induced from several hUCs. (PDF 340 kb)
showing pluripotent characterization of iPSCs induced using the 6F/BM1-4C system. a Flow cytometry for expression profiles of the hESC markers OCT4, SSEA4, Tra-1-60, and Tra-1-81. b Bisulfite sequencing assay for the methylation status of the Oct4 and Nanog promoters in iPSCs. Color codes indicate the proportion of methylation. y axis shows individual CpGs analyzed. x axis shows different cells. c Immunofluorescence assay for expression profiles of hESC markers. d Quantitative real-time PCR assay for expression profiles of marker genes of the three germ layers. e Hematoxylin and eosin staining of sections of iPSC-generated teratomas. f Scatter plots comparing global gene expression patterns between HN4 hESCs and UC1 iPSCs and between UC2 cells and UC2 iPSCs. Highlighted are the pluripotency factors Oct4, Sox2, Nanog, and miR-302a. Error bars indicate mean ± SD. *P < 0.05, **P < 0.01, ***P < 0.001. Scale bars, 100 μm. (PDF 375 kb)
showing morphology changes during iPSC generation using the 6F/BM1-4C system. a Morphology altered during iPSC generation using the 6F/BM1-4C system. b Schematic of episomal vectors used in the 4F2L-6C system. pEF1α EF1α promoter, pCMV CMV promoter. Scale bars, 100 μm. (PDF 158 kb)