- Open Access
Incomplete cellular reprogramming of colorectal cancer cells elicits an epithelial/mesenchymal hybrid phenotype
Journal of Biomedical Sciencevolume 25, Article number: 57 (2018)
Induced pluripotency in cancer cells by ectopic expression of pluripotency-regulating factors may be used for disease modeling of cancers. MicroRNAs (miRNAs) are negative regulators of gene expression that play important role in reprogramming somatic cells. However, studies on the miRNA expression profile and the expression patterns of the mesenchymal-epithelial transition (MET)/epithelial-mesenchymal transition (EMT) genes in induced pluripotent cancer (iPC) cells are lacking.
iPC clones were generated from two colorectal cancer (CRC) cell lines by retroviral transduction of the Yamanaka factors. The iPC clones obtained were characterized by morphology, expression of pluripotency markers and the ability to undergo in vitro tri-lineage differentiation. Genome-wide miRNA profiles of the iPC cells were obtained by microarray analysis and bioinformatics interrogation. Gene expression was done by real-time RT-PCR and immuno-staining; MET/EMT protein levels were determined by western blot analysis.
The CRC-iPC cells showed embryonic stem cell-like features and tri-lineage differentiation abilities. The spontaneously-differentiated post-iPC cells obtained were highly similar to the parental CRC cells. However, down-regulated pluripotency gene expression and failure to form teratoma indicated that the CRC-iPC cells had only attained partial pluripotency. The CRC-iPC cells shared similarities in the genome-wide miRNA expression profiles of both cancer and pluripotent embryonic stem cells. One hundred and two differentially-expressed miRNAs were identified in the CRC-iPC cells, which were predicted by bioinformatics analysis be closely involved in regulating cellular pluripotency and the expression of the MET/EMT genes, possibly via the phosphatidylinositol-3 kinases-protein kinase B (PI3K-Akt) and transforming growth factor beta (TGF-β) signaling pathways. Irregular and inconsistent expression patterns of the EMT vimentin and Snai1 and MET E-cadherin and occludin proteins were observed in the four CRC-iPC clones analyzed, which suggested an epithelial/mesenchymal hybrid phenotype in the partially reprogrammed CRC cells. MET/EMT gene expression was also generally reversed on re-differentiation, also suggesting epigenetic regulation.
Our data support the elite model for cancer cell-reprogramming in which only a selected subset of cancer may be fully reprogrammed; partial cancer cell reprogramming may also elicit an epithelial-mesenchymal mixed phenotype, and highlight opportunities and challenges in cancer cell-reprogramming.
Cellular reprogramming is a de-differentiation process that allows terminally-differentiated cells to acquire pluripotency [1, 2]. To date, numerous somatic cell-derived induced pluripotent stem cell (iPSC) lines have been established via retroviral transfection of the Yamanaka factors, OCT4, SOX2, KLF4 and c-MYC (OSKM) and by other means . Induced pluripotent cancer (iPC) cells are iPSCs of cancer origin, which also exhibit the distinctive iPSC characteristics of self-renewal and trilineage differentiation . Diminished cancer phenotypes have also been reported in the iPC cells [5, 6]; reduced in vitro and in vivo tumorigenicity  is thought to simulate the early stages of the tumorigenesis , opening the possibility that reprogrammed cancer cells may be developed into cancer disease models to recapitulate important early stages in cancer development for mechanistic studies, biomarker discovery and the development of novel therapeutic approaches to treat cancer .
In view of accumulated mutations and altered epigenetic programs in cancer cells, significant molecular and biological differences are predicted to exist between somatic and cancer- cell-derived reprogrammed iPSCs and iPCs. MicroRNAs (miRNAs) are important negative regulators of gene expression in most cellular processes , and have been implicated in regulating cell fate and maintaining pluripotency . Distinctive miRNA expression profiles have been shown to be associated with OSKM-mediated transcriptional networks and specific signaling pathways in somatic cells . However, systematic studies of miRNA expression and the biological changes elicited by the dysregulated miRNAs in reprogrammed cancer cells are still lacking.
In this study, we first established iPC clones from two colorectal cancer (CRC) cell lines. Using the CRC-iPC cells, we aimed to investigate the pluripotency features of the reprogrammed CRC cells, and possible miRNA involvement in regulating the epithelial-mesenchymal transition (EMT) and the reverse mesenchymal-epithelia transition (MET) processes, both of which are relevant to the establishment and maintenance of pluripotency in reprogrammed normal somatic and cancer cells [12,13,14,15]. Surprisingly, our data indicated that the CRC-iPC clones obtained did not attain full pluripotency, and analysis of EMT/MET protein expression further suggested that the partially reprogrammed CRC-iPC cells have acquired an epithelial/mesenchymal (E/M) hybrid phenotype. Our results support the elite model of cancer-cell reprogramming, and highlight intrinsic differences in reprogrammed normal somatic and cancer cells.
Cell lines and culture
A colorectal adenocarcinoma-derived cell line, HCT-15, and a second cell line, SK-CO-1, derived from a metastatic site (ascites) of a 65-year-old patient with colorectal adenocarcinoma, were used in this study. Both cell lines were obtained from ATCC (Manassas, VA, USA). The cell lines were cultured in DMEM high-glucose medium (HCT-15) or MEM (SK-CO-1) supplemented with 10% fetal bovine serum (FBS) (Gibco, Gaithersburg, MD, USA). Human embryonic stem cell (hESC) line H9 and iPC clones were cultured on mitomycin C-treated mouse embryonic fibroblasts (MEF) (Merck Millipore, Darmstadt, German) using a standard hESC medium. For feeder-free culture, cells were passaged to hESC-qualified matrigel (Corning, Avon, NY, USA)-coated culture vessels in mTeSR-1 (Stemcell Technologies, Vancouver, BC, Canada) with daily medium changes. All cultured cell lines were incubated in a humidified atmosphere at 37 °C with 5% CO2.
Retroviral transduction and iPC generation
pMXs retroviral constructs, each harboring one of the OCT4, SOX2, KLF4 or c-MYC genes, were amplified in 293FT cells and the supernatant was filtered through a 0.45-μm pore size PVDF filter as previously described [16, 17]. For retrovirus transduction of the CRC cell lines, the virus supernatant was added to plated CRC cells supplemented with 5 μg/ml polybrene (Merck Millipore, Darmstadt, Germany). The transfected cells were incubated for 24 h before a medium change. Upon reaching confluency, the OSKM-transduced cells were passaged to inactivated MEF. The next day, the medium was replaced with standard hESC medium and cultured using the hESC protocol until the emergence of hESC-like colonies after 21–23 days. Colonies were picked and transferred to fresh MEF feeder layer and continuously cultured in hESC medium .
The cells were fixed with 4% paraformaldehyde, incubated at room temperature for 30 min, followed by blocking for 2 h using 1% bovine serum albumin. The cells were washed twice with 1× PBS before addition of primary antibodies of the pluripotency markers, TRA-1-60, TRA-1-81, SSEA-4 or OCT4 (Stemcell Technologies) at 1:100 dilutions and incubated overnight at 4 °C. A FITC-conjugated rabbit anti-mouse antibody (Merck Millipore) was added and the mixture was further incubated for 1 h at room temperature. Nuclei were counterstained with DAPI (Gibco) and observed under an inverted fluorescent microscope.
In vitro lineage-directed and spontaneous differentiation
Putative CRC-derived induced pluripotent cancer (CRC-iPC) colonies were passaged to a 24-well plate pre-coated with hESC-qualified matrigel, and continuously cultured with osteogenic or adipogenic medium to induce mesoderm differentiation as described . The differentiation medium was changed every alternate day for 21–23 days before staining with Alizarin Red S or Oil Red O (Merck Millipore). For ectoderm-directed differentiation, the putative iPCs were cultured in DMEM/F12 medium, 10% FBS supplemented with 100 ng/ml Noggin (R&D Systems, Minneapolis, MN, USA) for 1 week. For endoderm lineage differentiation, CRC-iPCs were cultured with DMEM/F12 medium with 10% FBS supplemented with 100 ng/ml Activin A (R&D Systems). Ectoderm (MAP2) and endoderm (AFP) (Merck Millipore) markers were also used and observed by immunofluorescence staining. For in vitro spontaneous differentiation, iPCs were cultured in suspension culture in a standard hESC medium for 7 days. The embryoid bodies formed were transferred to 0.1% gelatin coated-culture dish for attachment and further differentiation in FGF-2-free hESC medium. On day 14, the attached embryoid bodies were subcultured to form post-iPC cells in DMEM/F12, 10% FBS and 1% penicillin-streptomycin .
MiRNA microarray and RNA analysis
Total RNA preparations, isolated using the miRNeasy kit (Qiagen, Redwood, CA, USA), were subjected to analysis using Agilent SurePrint Human MiRNA Microarray Release 21.0 (Agilent Technologies, Santa Clara, CA, USA), which consisted of 2549 miRNAs as annotated in miRBase 21. The dataset obtained were analyzed using the GeneSpring GX software version 13.0 (Agilent). Comparative analysis between samples was carried out using the t-test (p-values) and the Benjamini-Hochberg False Discovery Rate (FDR) correction (adjusted p-values) to remove false-positive miRNAs. The criteria set for selection of differentially-expressed miRNAs was log2(fold change) ≥ 2 or ≤ − 2 with adjusted p-values < 0.05. Hierarchical clustering was performed to identify and visualize pattern of miRNAs expression between samples. The microarray dataset has been deposited with NCBI GEO with the accession number GSE87280. Quantification of mRNA and miRNA expression levels was performed as previously described [18, 19] using primers as shown in Additional file 1: Table S1. For quantitative RT-PCR, the data were analyzed using the comparative ΔΔCt method, normalizing to the GAPDH expression level for mRNAs and that of the RNU6 gene for miRNAs. In all experiments, reactions were performed in triplicates and the data presented as mean ± standard error of the mean derived from three independent experiments. Data were analyzed by Student’s t test. P < 0.05 was considered statistically significant.
Western blot analysis
Protein lysates were harvested from the cells using 1× RIPA lysis buffer (Nacalai Tesque, Kyoto, Japan) and kept on ice for 30 min before being centrifuged at 10,000 g for 10 min at 4 °C to obtain the supernatant. Protein lysates (50 μg) were subjected to SDS-polyacrylamide gel electrophoresis and western blot analysis as described . Dilutions and the sources of the antibodies used are as follows: E-cadherin (CDH1, 1:1000, Cell Signaling Technology), occludin (OCLN, 1:1000, Merck Millipore), vimentin (VIM, 1:1000, Cell Signaling), SNAI1 (1:250, Cell Signaling) and GAPDH (1:2000, Cell Signaling). The protein bands were detected using a horseradish peroxidase-conjugated secondary antibody (1: 10,000, Abcam, Cambridge, UK) for 1 h at room temperature and visualized with the Amersham ECL Western blotting substrate (GE Healthcare), according to the manufacturer’s protocol. The mean of multiple blots is presented (Fig. 4b); the error bars represent standard errors of the mean, SEM.
Reprogrammed colorectal cancer (CRC)-induced pluripotent cancer (iPC) cells show ESC-like features and trilineage differentiation
In this study, the colorectal cancer cell lines HCT-15 and SK-CO-1 were used in cellular reprogramming via retroviral transduction of the Yamanaka factors, OCT4, SOX2, KLF4 and c-MYC (OSKM) as previously described [2, 16]. The OSKM-transduced cells formed rounded colonies at 21–23 days post-transduction (Fig. 1a). The colonies were picked and were designated as iHCT-15 and iSK-CO-1 for colonies derived from HCT-15 and SK-CO-1, respectively; different independently isolated clones were further designated with the prefix “C” followed by a clone number. Both the iHCT-15 and iSK-CO-1 colonies showed embryonic stem cell (ESC)-like morphology (Fig. 1a, leftmost panels). However, the cells tended to pile up to form multilayers in the center of the putative induced pluripotent cancer (iPC) colonies. Furthermore, the borders of the colonies were not as well-defined as the ESC colonies. The putative CRC-iPCs were stained positive for the pluripotency markers, TRA-1-60, TRA-1-81 and alkaline phosphatase, and a pluripotency marker, OCT4 (Fig. 1a) . When cultured in appropriate differentiation media, the CRC-iPC cells differentiated into osteocyte- and adipocyte-like cells of the mesoderm lineage (Fig. 1b, left panels). In the presence of Noggin and Activin A, the iPCs differentiated towards the ectoderm and endoderm lineages, respectively, as evidenced by immunofluorescence staining of the germ-layer markers, microtubule-associated protein 2 (MAP2, ectoderm) and alpha-fetoprotein (AFP, endoderm) (Fig. 1b, right panels). Using a previously published protocol , the CRC-iPCs were also spontaneously differentiated into post-iPCs via embryoid body formation (Fig. 1c, middle and right panels) . The post-iPCs grew in monolayer and were morphologically similar to the parental CRCs (Fig. 1c, left panel). The differentiated post-iPCs re-expressed the germ-layer markers, coudal type homeobox protein 2 (CDX2; ectoderm), Msh homeobox 1 (MSX1; mesoderm) and GATA binding protein-6 (GATA6; endoderm), which were initially significantly down-regulated in the iPC cells on reprogramming (Fig. 1d). The observations suggest reversion to a more primitive state when cancer cells are subjected to reprogramming.
Down-regulated but reversible expression of pluripotency genes on CRC reprogramming and re-differentiation
After OSKM transduction, it was important to distinguish between ectopic and endogenous OSKM expression. For this purpose, RT-PCR was conducted using primers specific for the transgenes using total RNA isolated from iPC cells at passage 25 or later. In each case, two independent iHCT-15 clones C1 and C5, and iSK-CO-1 clones C1 and C2 were analyzed. The results showed that expression of the transgenes was extinguished with one notable exception (Fig. 2a). In iSK-CO-1 clone C2, SOX2 and c-MYC transgene sequences were still detectable, indicating continued ectopic transgene expression. Persistent post-transduction ectopic transgene expression has been reported and interpreted as to indicate incomplete reprogramming [13, 21]. In the RT-PCR analysis, the control H9 embryonic stem cells (ESC) expressed the endogenous levels of the pluripotency genes as anticipated (Fig. 2b). It is, however, noted that the parental HCT-15 and SK-CO-1 cells were already endogenously expressing the pluripotency genes (Fig. 2b), consistent with other reports that pluripotency genes are frequently activated in cancer cell lines [22, 23].
Endogenous expression levels of the pluripotency genes were quantified by real-time RT-PCR (Fig. 2c). Post-iPC cell lines were also included in the pluripotency gene expression analysis. The results showed that, the expression of the pluripotency genes was generally down-regulated in all four CRC-iPC clones, except for SOX2, which was up-regulated in the two iHCT-15 clones. Moreover, on re-differentiation to post-iPC cells, expression of the pluripotency genes was generally up-regulated relative to the iPC clones, and often seemed to be reversed back to the parental levels. In the case of REX1, re-differentiation did not seem to have altered the expression levels. A notable exception of expression reversal was found in KLF4 in which its expression was further down-regulated on re-differentiation. Taken together, our data indicated that expression of the pluripotency genes was generally down-regulated on reprogramming, but was not totally shutdown, which raises the possibility of partial or incomplete reprogramming [13, 21]. Such down-regulated pluripotency gene expression was reversible, consistent with epigenetic regulation, and was a clear sign of alteration of the cancer phenotype at the molecular level. It is also noteworthy that the CRC-iPC clones thus derived had been maintained for up to 60 passages (data not shown), implying that the maintenance of pluripotency and self-renewal properties of the reprogrammed colorectal cancer cells may not be entirely dependent on the expression levels of the pluripotency genes.
CRC-iPC cells share similarities in the miRNA expression profiles of both cancer and pluripotent embryonic stem cells
Previous reports have shown that reprogramming alters the miRNA expression profile in induced pluripotent stem cells (iPSCs) derived from both somatic and cancer cells [10, 24]. In addition, miRNA involvement in cellular reprogramming has also been reported , indicating important miRNA roles in the induction and maintenance of pluripotency and cell fate [19, 25, 26]. However, there are few studies on genome-wide miRNA expression profiling in reprogrammed cancer cells. For this purpose, total RNA preparations of the HCT-15 and SK-CO-1 cells and two iPC clones from each cell line were subjected to miRNA microarray analysis (NCBI GEO accession number: GSE87280). A human ESC cell line, H9, was also included as a reference cell line.
Hierarchical clustering analysis of the miRNA profiles showed that the differentially-expressed miRNAs could be grouped into three major clusters based on their expression patterns (Fig. 3a; see also Additional file 2: Table S2). Cluster I miRNAs were highly or moderately expressed in ESC cells, consistent with the inclusion in this cluster the known reprogramming miR-130 , the mir-302 family (miR-302a, −302b, −302c & -302d [28, 29], and the chromosome 19 miRNA cluster (C19MC) miRNAs, which have previously been shown to modulate the early phases of reprogramming () (Fig. 3a & Additional file 2: Table S2). On the other hand, the cluster I miRNAs were not expressed, or expressed at low levels, in the parental CRC cells, and that these miRNAs remained unexpressed on reprogramming. Cluster II consists of miRNAs that were unexpressed, or expressed at low expression levels, in CRC cells, but were up-regulated in the iPC cells to levels comparable to H9 ESC. However, some cluster II miRNAs, including miR-642a and − 3162-5p, were already highly expressed in CRC, and remained consistently expressed iPC cells as were in the ESC cells (Fig. 3a & Additional file 2: Table S2). The results suggest that cluster II miRNAs are associated with pluripotency, but are already activated in CRC cells. Cluster III includes miRNAs that were already highly expressed in the parental CRC cells and expression levels remained high, or were up-regulated on reprograming in the CRC-iPC cells to levels higher than in the ESC cells, which showed only low expression levels, suggesting that these miRNAs are associated more with tumorigenesis than with the reprogramming process. Cluster III includes the let-7 and miR-200 (miR-200a,-200b,-200c) families (Fig. 3a & Additional file 2: Table S2). The let-7 family suppresses pluripotency by modulating expression of pluripotency genes, including c-MYC and Lin28 , whereas the miR-200 cluster is important in the early stage of somatic cell reprogramming by activating mesenchymal-to-epithelial transition (MET) .
Taken together, miRNA expression profiling revealed that the reprogrammed CRC-iPC cells retained some molecular signatures of the parental cells, indicating retention of cancer phenotype at the molecular level. On the other hand, despite the demonstrated pluripotency, CRC-iPC cells showed distinct dissimilarity with the pluripotent ESC cells, hinting incomplete reprogramming due to the tumor phenotype of CRC.
Reprogramming activates miRNAs targeting the TGF-β cell signaling pathway to regulate expression of EMT/MET genes
MiRNA profiling showed that the miRNA expression profiles were highly similar between the four iPC clones, irrespective of the parental sources (Fig. 3a). Using a stringent selection criteria of log2(fold change) (FC) ≥ 2.0 or ≤ − 2.0 and p-value < 0.05, a total of 102 statistically significant differentially-expressed miRNAs were thus identified, fifty of which were down-regulated and fifty-two miRNAs were up-regulated (Fig. 3b; for a full miRNA list, see Additional file 3: Table S3). The top ten most significantly up- or down-regulated miRNAs are shown in Table 1. Interestingly, three of the top ten up-regulated miRNA in the CRC-iPC cells are mapped in the highly unstable regions of the 19p13.13 and 19q13.11 deletion syndromes [32, 33], whereas the gene encoding miR-362-5p in the down-regulated group is mapped at Xp11.23 associated with the Xp11.22-Xp11.23 duplication syndrome . About half of the top 10 miRNAs have not been annotated.
To elucidate possible functional roles, putative targets of the differentially expressed miRNAs were identified by interrogating the TargetScan and MicroRNA.org databases and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were mapped, which showed that six of the top 10 KEGG pathways were related to cellular signaling (Fig. 3c), including the TGF-β pathway involved in the regulation of the EMT processes . MET is a critical step in setting up the initial phase of somatic cell reprogramming; hence, the expression status of MET and EMT genes is an important indicator of the cell state . On database interrogations, 25 of the differentially expressed miRNAs, both in the up- and down-regulated groups, were found to target transcripts of two representative EMT genes, vimentin (VIM) and Snail1 (SNAI1) and two typical MET genes, E-cadherin (CDH1) and occludin (OCLN) (Table 2). For validation, miRNAs targeting SNAI1 and CDH1 were randomly chosen for qRT-PCR analysis (Fig. 3d). The predicted expression was largely consistent with the exception that miR-194-3p was slightly up-regulated in the HCT-15 iPC clones when down-regulation was shown in microarray analysis. The minor discrepancy may be explained by the fact that in the miRNA profiling analysis (Table 1 and Additional file 2: Table S2), mean values of four iPC clones from the two parental CRC cell lines were derived, whereas individual iPC clones were subjected to qRT-PCR analysis (Fig. 3d).
Target prediction data (Table 2) showed that each of the four EMT/MET genes was targeted by multiple miRNAs that were both up- or down-regulated on reprogramming. Some miRNAs also targeted two or more genes. Notably, miR-335-3p was down-regulated on reprogramming and the miRNA was predicted to target an EMT gene SNAI1, and both the MET CDH1 and OCLN genes (Table 2). Both the miR-362-5p and -3p species, which have different sequences, were down-regulated and targeted different EMT and MET gene.
The TGF-β pathway, which has previously been shown to regulate the EMT/MET processes , is one of the top ten predicted KEGG pathways targeted by the reprogramming-induced dysregulated miRNAs (Fig. 3c) To further establish possible relationship between the miRNAs and the TGF-β pathway, targeting miRNAs were mapped to the major components of the pathways in relation to EMT regulation (Fig. 3e). In this study, bioinformatics analysis predicted that TGFβ2, a ligand of the TGFβ family , was regulated by ten miRNAs, seven of which were down-regulated and three were up-regulated on reprogramming (Fig. 3e). Furthermore, SMAD4, which complexes with phosphorylated R-SMAD further to activate or enhance transcription of downstream effector genes , was also predicted to be targeted by seven differentially expressed miRNAs. The down-regulation of these miRNAs may up-regulate TGFβ2 and SMAD4 expression, the TGF-β signaling cascade and expression of the EMT-associated genes, SNAI1, SLUG and ZEB1, resulting in a cell transition towards the mesenchymal-like state . In addition, up-regulated SNAI1 expression may also lead to repression of the epithelial E-cadherin and occludin proteins, resulting in the reversal of the MET process . Taken together, many of the cancer cell-reprogramming-induced dysregulated miRNAs are predicted to be associated with EMT/MET regulation, and subtle interplay between the up- and down-regulated miRNAs may decide the predominance of either or both the EMT and MET processes.
Dysregulated expression of EMT and MET proteins on reprogramming suggests an E/M hybrid phenotype
Effects of reprogramming on EMT/MET gene expression in CRC cells were next investigated by western blot analysis in the iPC and post-iPC cells, focusing on the same EMT proteins, vimentin and Snail1, and MET proteins, E-cadherin and occludin, examined above. Representative western blots are shown in Fig. 4a; quantification of the data from three independently performed western blots is shown in Fig. 4b. Vimentin levels were low in the parental HCT-15 and SK-CO-1 cells. On reprogramming, the protein levels were up-regulated between ~ 2.5- to ~ 8-fold in the different iPC clones of both cell lines. On differentiation into post-iPC clones, vimentin levels were generally suppressed or reverted back to levels similar to those of the parental cells. SNAI1 level was high in HCT-15 but was very low in SK-CO-1. On reprogramming, the SNAI1 levels were suppressed in the iHCT-15 clones, but were up-regulated to various extents in the two iSK-CO-1 clones. On differentiation, the SNAI1 levels were restored to approximately the parental levels in post-iHCT-15 clones, but the protein was further up-regulated in post-iSK-CO-1 cells. In the E-cadherin western blots, a 120-kD band that was detected all cell types (Fig. 4a, E-cadherin, arrowhead); however, a second 86-kD species appeared in the iPC cells, which was most likely the E-cadherin soluble ectodomain, a cleavage product of matrix metalloproteinases of the 120-kD protein [40, 41]. E-cadherin was consistently up-regulated in all iPC clones from the two different cell lines; on differentiation, the protein levels were all reverted back to a level comparable to that of the parental CRC cells. Expression of occludin was clearly detected in the HCT-15 cells, but the protein was present only in a very low level in SK-CO-1. Occludin levels were down-regulated in the iHCT-15 iPC clones, and were restored to approximately the parental levels on differentiation. On the other hand, occludin was up-regulated in the reprogrammed iSK-CO-1 cells, and differentiation further up-regulated the protein to higher levels.
Paradoxically, the western blot data revealed that on reprogramming, both the EMT and MET proteins were generally up-regulated, with exception in the SNAI1 and, to a lesser extent, occludin, in the iHCT-15 cells. Hence, the conflicting EMT and MET expression patterns suggested the possibility that reprogramming of cancer cells might have elicited an epithelial/mesenchymal (E/M) hybrid phenotype [42, 43]. The data also showed that initial differentiation of the iPC clones generally reversed expression of the EMT/MET protein, echoing epigenetic regulation, but with the exception of SNAI1 and occludin in the reprogrammed iSK-CO-1 clones. Different clones of the same cell line, and different CRC cell lines also showed different EMT/MET expression profiles, which could be explained by cellular heterogeneity and cancer stage-dependent expression modes of the EMT/MET genes in the cancer cells [44, 45].
In this work, the OSKM-reprogrammed CRC-iPC cells were able to undergo in vitro tri-lineage differentiation showed other pluripotency characteristics similar to those of embryonic stem cells. However, several lines of experimental evidence suggested that the CRC-iPC clones obtained had not undergone complete reprogramming. Firstly, the border of the iPC colonies formed was not well-defined, and the cells piled up in the center (Fig. 1a). Secondly, in some CRC-iPC clones, ectopic expression of the transgenes had not been completely extinguished (Fig. 2a). In other reports, partially-reprogrammed somatic cells, which are also called pre-iPSC, are often characterized by low levels of endogenous pluripotency genes while the exogenous transgenes may still be sustained in the pre-iPSC [13, 21], as observed in this work for the iSK-CO-1 clone 2, which failed to attenuate ectopic expression of SOX2 and c-MYC (Fig. 2a), suggesting incomplete reprogramming. However, residual transgene expression is not specific to retroviral transduction method  since OSKM-retroviral transduced sarcoma  and colorectal cancer cells  frequently showed efficient and complete silencing of the OSKM transgenes. A third line of evidence to support incomplete reprogramming is that the reprogrammed CRC-iPC cells showed distinct dissimilarity in the miRNA profiles compared with the pluripotent H9 ESC cells, while retaining similarity in the molecular signature with the parental cancer cells (Fig. 3a). Furthermore, the CRC-iPC cells failed to form teratoma in vivo in nude mice (Choo et al., unpublished data;). Indeed, others have also shown that reprogrammed cancer cells frequently fail to form teratoma, suggesting that limited pluripotency may be a common feature in reprogrammed cancer cells [26, 47]. In addition, members of the miR-302/− 367 families, which were previously used to generate iPSCs and iPCs , were not found among the 102 differentially-expressed miRNAs in the CRC-iPCs generated (Additional file 3: Table S3). On the contrary, the miRNA expression profiles indicated miRNA involvement in the regulation of EMT (Table 2 and Fig. 3e), which raises the possibility of an alternate route in cancer cell-reprogramming as discussed above.
Furthermore, the six pluripotency genes examined, including OSKM, NANOG and REX1, are already expressed in the parental CRC cells, as reported for other cancer cells [22, 23]. In the reprogrammed CRC-iPC cells, expression of these genes was generally down-regulated, which was reversible on re-differentiation (Fig. 2). Others have shown down-regulation, or limited activation, of endogenous pluripotency genes, in somatic cell reprogramming .
The finding that endogenous expression of the pluripotency genes was down-regulated on CRC reprogramming may, in fact, be linked to the epithelial or mesenchymal state of the reprogrammed cells, as proposed in the scheme shown in Fig. 5. OCT4 and SOX2 have previously been shown to inhibit expression of the EMT-inducing factor, SNAI1, which works synergistically with KLF4 to induce an epithelial phenotype by up-regulating E-cadherin . However, KLF4 is, in turn, negatively regulated by SNAI1 . Thus, the up-regulation of SNAI1 may have inhibited KLF4, leading to down-regulation of E-cadherin and, hence, the suppression of MET. However, it is unclear how down-regulation of NANOG could promote EMT. In addition, the core pluripotency genes OSKM have previously been shown to inhibit various factors of the TGF-β signaling pathway, and that OCT4 and SOX2 suppress SNAI1 . Hence, OSKM down-regulation is likely to have contributed to the TGF-β signaling pathway to up-regulate SNAI1 expression and, therefore, promoting EMT (Fig. 5). Pluripotency gene-modulated transition between MET and EMT may, therefore, play a vital role in inducing pluripotency in cells of cancer origin.
The strongest evidence to suggest partial reprogramming came from the expression analysis of the EMT/MET genes in the reprogrammed iPC cells (Fig. 4). As pointed out above, previous studies have shown that MET is an essential early step in somatic cell reprogramming [14, 51], whereas EMT counterbalances the MET effects to hinder reprogramming . In the reprogrammed CRC-iPC clones, the four EMT and MET proteins investigated were not consistently negatively or positively regulated on reprogramming to clearly confer either an epithelial or mesenchymal phenotype to the iPC cells, and the reprogramming-induced alterations in EMT and MET gene expression were generally reversible (Fig. 4). The western blot data could be interpreted to indicate that the reprogrammed cancer cells had acquired an intermediate state of an epithelial/mesenchymal (E/M) hybrid phenotype , as has previously been reported in adult epithelial, hepatic and neoplastic stem cells . Others have also developed models to suggest that such a mixed E/M phenotype is associated with multiple stemness characteristics of cancer stem cells, including tumor initiating and self-renewal abilities [42, 54, 55].
Increasing evidences are also suggesting that EMT may also be an important participant in regulating cellular reprogramming . Sequential introduction of the OSKM factors induced an early and transient response of EMT prior to the appearance of a delayed MET response in reprogrammed somatic cells . Hence, beside the proposed E/M hybrid phenotype, there may also exist a delicate and dynamic balance between EMT and MET per se to confer either an EMT or MET phenotype to reprogrammed cancer cells. In terms of miRNA regulation, the four MET/EMT genes analyzed here are also targeted by miRNAs both up- or down-regulated on reprogramming (Table 2 & Fig. 4). Hence, expression of the EMT/MET proteins in the iPC cells may further dependent on delicate interplays between multiple targeting miRNAs and other regulatory factors.
Zhang et al. (2013) have previously proposed a pluripotency hierarchy model in which the cancer cell-derived iPCs are thought to be induced to a pre-iPSC state, just a level above the multipotent mesenchymal stem cells . The pre-iPSC state of reprogrammed cancer cells is also thought to be reflected in the down-regulation of pluripotency genes . Hence, reprogramming of normal somatic cells follows a continuous stochastic model, where all transduced cells have equal probability to be transformed into a pluripotent state [57, 58]. In contrast, evidences have been presented to support the elite model for cancer-cell reprogramming in which only a selected subset of cells in the heterogeneous cancer cell population may be fully reprogrammed into iPCs . Our work here highlights challenges and in fully reprogramming cancer cells.
In this work, Yamanaka factor-reprogrammed colorectal cancer-induced pluripotent cancer cells show ESC-like features and trilineage differentiation. Down-regulated expression of pluripotency genes and dissimilar miRNA expression profiles to that of pluripotent embryonic stem cells indicate incomplete reprogramming. Reprogramming activates miRNAs that target the TGF-β signaling pathway to modulate expression of EMT/MET genes; however, dysregulated expression of EMT and MET proteins suggests an epithelial-mesenchymal hybrid phenotype, consistent with partial CRC reprogramming. Our data further highlight challenges in obtaining fully reprogrammed cancer cells likely due to the accumulated mutations and epigenetic modifications in the cancer cells. Nonetheless, the reprogrammed CRC-iPC cells offer an opportunity for the development of disease models to further elucidate the molecular mechanism in the tumorigenesis of colorectal cancer.
Coudal type homeobox protein 2
Embryonic stem cells
GATA binding protein-6
Induced pluripotent cancer
Induced pluripotent stem cell
Microtubule-associated protein 2
Mouse embryonic fibroblasts
Msh homeobox 1
OCT4, SOX2, KLF4 and c-MYC
Takahashi K, Yamanaka S. Induction of pluripotent stem cells from mouse embryonic and adult fibroblast cultures by defined factors. Cell. 2006;126(4):663–76.
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(5):861–72.
Vitale AM, Wolvetang E, Mackay-Sim A. Induced pluripotent stem cells: a new technology to study human diseases. Int J Biochem Cell Biol. 2011;43(6):843–6.
Ramos-Mejia V, Fraga MF, Menendez P. iPSCs from cancer cells: challenges and opportunities. Trends Mol Med. 2012;18:3.
Miyoshi N, Ishii H, Nagai K, Hoshino H, Mimori K, Tanaka F, Nagano H, Sekimoto M, Doki Y, Mori M. Defined factors induce reprogramming of gastrointestinal cancer cells. Proc Natl Acad Sci U S A. 2010;107(1):40–5.
Miyazaki S, Yamamoto H, Miyoshi N, Wu X, Ogawa H, Uemura M, Nishimura J, Hata T, Takemasa I, Mizushima T, et al. A Cancer reprogramming method using MicroRNAs as a novel therapeutic approach against Colon Cancer: research for reprogramming of Cancer cells by MicroRNAs. Ann Surg Oncol. 2015;22(Suppl 3):S1394–401.
Zhang X, Cruz FD, Terry M, Remotti F, Matushansky I. Terminal differentiation and loss of tumorigenicity of human cancers via pluripotency-based reprogramming. Oncogene. 2013;32(18):2249–60. 2260.e1–21
Kim J, Zaret KS. Reprogramming of human cancer cells to pluripotency for models of cancer progression. EMBO J. 2015;34(6):739–47.
Wahid F, Shehzad A, Khan T, Kim YY. MicroRNAs: synthesis, mechanism, function, and recent clinical trials. Biochim Biophys Acta. 2010;1803(11):1231–43.
Luningschror P, Hauser S, Kaltschmidt B, Kaltschmidt C. MicroRNAs in pluripotency, reprogramming and cell fate induction. Biochim Biophys Acta. 2013;1833(8):1894–903.
Tanabe S. Signaling involved in stem cell reprogramming and differentiation. World J Stem Cells. 2015;7(7):992–8.
Liu X, Sun H, Qi J, Wang L, He S, Liu J, Feng C, Chen C, Li W, Guo Y, et al. Sequential introduction of reprogramming factors reveals a time-sensitive requirement for individual factors and a sequential EMT-MET mechanism for optimal reprogramming. Nat Cell Biol. 2013;15(7):829–38.
Buganim Y, Faddah DA, Jaenisch R. Mechanisms and models of somatic cell reprogramming. Nat Rev Genet. 2013;14(6):427–39.
David L, Polo JM. Phases of reprogramming. Stem Cell Res. 2014;12(3):754–61.
Takaishi M, Tarutani M, Takeda J, Sano S. Mesenchymal to epithelial transition induced by reprogramming factors attenuates the malignancy of cancer cells. PLoS One. 2016;11(6):e0156904.
Sugii S, Kida Y, Berggren WT, Evans RM. Feeder-dependent and feeder-independent iPS cell derivation from human and mouse adipose stem cells. Nat Protoc. 2011;6(3):346–58.
Choong PF, Teh HX, Teoh HK, Ong HK, Choo KB, Sugii S, Cheong SK, Kamarul T. Heterogeneity of osteosarcoma cell lines led to variable responses in reprogramming. Int J Med Sci. 2014;11(11):1154–60.
Nguyen PN, Huang CJ, Sugii S, Cheong SK, Choo KB. Selective activation of miRNAs of the primate-specific chromosome 19 miRNA cluster (C19MC) in cancer and stem cells and possible contribution to regulation of apoptosis. J Biomed Sci. 2017a;24(1):20.
Nguyen PNN, Choo KB, Huang CJ, Sugii S, Cheong SK, Kamarul T. miR-524-5p of the primate-specific C19MC miRNA cluster targets TP53IPN1- and EMT-associated genes to regulate cellular reprogramming. Stem Cell Res Ther. 2017b;8(1):214.
Zhao W, Ji X, Zhang F, Li L, Ma L. Embryonic stem cell markers. Molecules. 2012;17(6):6196–236.
Plath K, Lowry WE. Progress in understanding reprogramming to the induced pluripotent state. Nat Rev Genet. 2011;12(4):253–65.
Schoenhals M, Kassambara A, De Vos J, Hose D, Moreaux J, Klein B. Embryonic stem cell markers expression in cancers. Biochem Biophys Res Commun. 2009;383(2):157–62.
Amini S, Fathi F, Mobalegi J, Sofimajidpour H, Ghadimi T. The expressions of stem cell markers: Oct4, Nanog, Sox2, nucleostemin, Bmi, Zfx, Tcl1, Tbx3, Dppa4, and Esrrb in bladder, colon, and prostate cancer, and certain cancer cell lines. Anat Cell Biol. 2014;47(1):1–11.
Koga C, Kobayashi S, Nagano H, Tomimaru Y, Hama N, Wada H, Kawamoto K, Eguchi H, Konno M, Ishii H, et al. Reprogramming using microRNA-302 improves drug sensitivity in hepatocellular carcinoma cells. Ann Surg Oncol. 2014;21(Suppl 4):S591–600.
Anokye-Danso F, Snitow M, Morrisey EE. How microRNAs facilitate reprogramming to pluripotency. J Cell Sci. 2012;125(Pt 18):4179–87.
Nishikawa S, Ishii H, Haraguchi N, Kano Y, Fukusumi T, Ohta K, Ozaki M, Dewi DL, Sakai D, Satoh T, et al. MicroRNA-based cancer cell reprogramming technology. Exp Ther Med. 2012;4(1):8–14.
Pfaff N, Fiedler J, Holzmann A, Schambach A, Moritz T, Cantz T, Thum T. miRNA screening reveals a new miRNA family stimulating iPS cell generation via regulation of Meox2. EMBO Rep. 2011;12(11):1153–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(10):2115–24.
Sandmaier SE, Telugu BP. MicroRNA-mediated reprogramming of somatic cells into induced pluripotent stem cells. Methods Mol Biol. 2015;1330:29–36.
Li MA, He L. microRNAs as novel regulators of stem cell pluripotency and somatic cell reprogramming. Bioessays. 2012;34(8):670–80.
Wang G, Guo X, Hong W, Liu Q, Wei T, Lu C, Gao L, Ye D, Zhou Y, Chen J, et al. Critical regulation of miR-200/ZEB2 pathway in Oct4/Sox2-induced mesenchymal-to-epithelial transition and induced pluripotent stem cell generation. Proc Natl Acad Sci. 2013;110(8):2858–63.
Malan V, Raoul O, Firth HV, Royer G, Turleau C, Bernheim A, Willatt L, Munnich A, Vekemans M, Lyonnet S, et al. 19q13.11 deletion syndrome: a novel clinically recognisable genetic condition identified by array comparative genomic hybridisation. J Med Genet. 2009;46(9):635–40.
Dolan M, Mendelsohn NJ, Pierpont ME, Schimmenti LA, Berry SA, Hirsch B. A novel microdeletion/microduplication syndrome of 19p13.13. Genet Med. 2010;12(8):503–11.
Giorda R, Bonaglia MC, Beri S, Fichera M, Novara F, Magini P, Urquhart J, Sharkey FH, Zucca C, Grasso R, et al. Complex segmental duplications mediate a recurrent dup (X) (p11.22-p11.23) associated with mental retardation, speech delay, and EEG anomalies in males and females. Am J Hum Genet. 2009;85(3):394–400.
Tan EJ, Olsson AK, Moustakas A. Reprogramming during epithelial to mesenchymal transition under the control of TGFbeta. Cell Adhes Migr. 2015;9(3):233–46.
Chen J, Han Q, Pei D. EMT and MET as paradigms for cell fate switching. J Mol Cell Biol. 2012;4(2):66–9.
Massague J. TGFbeta signalling in context. Nat Rev Mol Cell Biol. 2012;13(10):616–30.
Miyazono K. Transforming growth factor-β signaling in epithelial-mesenchymal transition and progression of cancer. Proc Jpn Acad Ser B Phys Biol Sci. 2009;85(8):314–23.
Lamouille S, Xu J, Derynck R. Molecular mechanisms of epithelial-mesenchymal transition. Nat Rev Mol Cell Biol. 2014;15(3):178–96.
Matos ML, Lapyckyj L, Rosso M, Besso MJ, Mencucci MV, Briggiler CI, Giustina S, Furlong LI, Vazquez-Levin MH. Identification of a novel human E-cadherin splice variant and assessment of its effects upon EMT-related events. J Cell Physiol. 2017;232(6):1368–86.
van Roy F, Berx G. The cell-cell adhesion molecule E-cadherin. Cell Mol Life Sci. 2008;65(23):3756–88.
Jolly MK, Boareto M, Huang B, Jia D, Lu M, Ben-Jacob E, Onuchic JN, Levine H. Implications of the hybrid epithelial/mesenchymal phenotype in metastasis. Front Oncol. 2015;5:155.
Garg M. Urothelial cancer stem cells and epithelial plasticity: current concepts and therapeutic implications in bladder cancer. Cancer Metastasis Rev. 2015;34(4):691–701.
Loboda A, Nebozhyn MV, Watters JW, Buser CA, Shaw PM, Huang PS, Van't Veer L, Tollenaar RA, Jackson DB, Agrawal D, et al. EMT is the dominant program in human colon cancer. BMC Med Genet. 2011;4:9.
Heerboth S, Housman G, Leary M, Longacre M, Byler S, Lapinska K, Willbanks A, Sarkar S. EMT and tumor metastasis. Clin Transl Med. 2015;4:6.
Montserrat N, Garreta E, Gonzalez F, Gutierrez J, Eguizabal C, Ramos V, Borros S, Izpisua Belmonte JC. Simple generation of human induced pluripotent stem cells using poly-beta-amino esters as the non-viral gene delivery system. J Biol Chem. 2011;286(14):12417–28.
Lim KL, Teoh HK, Choong PF, Teh HX, Cheong SK, Kamarul T. Reprogramming cancer cells: overview & current progress. Expert Opin Biol Ther. 2016;16(7):941–51.
Kuo CH, Ying SY. Advances in microRNA-mediated reprogramming technology. Stem Cells Int. 2012;2012:823709.
Li R, Liang J, Ni S, Zhou T, Qing X, Li H, He W, Chen J, Li F, Zhuang Q, et al. A mesenchymal-to-epithelial transition initiates and is required for the nuclear reprogramming of mouse fibroblasts. Cell Stem Cell. 2010;7(1):51–63.
De Craene B, Gilbert B, Stove C, Bruyneel E, van Roy F, Berx G. The transcription factor snail induces tumor cell invasion through modulation of the epithelial cell differentiation program. Cancer Res. 2005;65(14):6237–44.
Samavarchi-Tehrani P, Golipour A, David L, Sung HK, Beyer TA, Datti A, Woltjen K, Nagy A, Wrana JL. Functional genomics reveals a BMP-driven mesenchymal-to-epithelial transition in the initiation of somatic cell reprogramming. Cell Stem Cell. 2010;7(1):64–77.
Maherali N, Hochedlinger K. Tgfbeta signal inhibition cooperates in the induction of iPSCs and replaces Sox2 and cMyc. Curr Biol. 2009;19(20):1718–23.
Jia D, Jolly MK, Tripathi SC, Den Hollander P, Huang B, Lu M, Celiktas M, Ramirez-Peña E, Ben-Jacob E, Onuchic JN, et al. Distinguishing mechanisms underlying EMT tristability. Cancer Convergence. 2017;1(1):2–19.
Strauss R, Li ZY, Liu Y, Beyer I, Persson J, Sova P, Moller T, Pesonen S, Hemminki A, Hamerlik P, et al. Analysis of epithelial and mesenchymal markers in ovarian cancer reveals phenotypic heterogeneity and plasticity. PLoS One. 2011;6(1):e16186.
Jolly MK, Tripathi SC, Somarelli JA, Hanash SM, Levine H. Epithelial-mesenchymal plasticity: how have quantitative mathematical models helped improve our understanding? Mol Oncol. 2017;11(7):739–54.
De Chiara L, Crean J. Emerging transcriptional mechanisms in the regulation of epithelial to mesenchymal transition and cellular plasticity in the kidney. J Clin Med. 2016;5(1):6.
Hanna J, Saha K, Pando B, van Zon J, Lengner CJ, Creyghton MP, van Oudenaarden A, Jaenisch R. Direct cell reprogramming is a stochastic process amenable to acceleration. Nature. 2009;462(7273):595–601.
Yamanaka S. Elite and stochastic models for induced pluripotent stem cell generation. Nature. 2009;460(7251):49–52.
Lai J, Kong CM, Mahalingam D, Xie X, Wang X. Elite model for the generation of induced pluripotent cancer cells (iPCs). PLoS One. 2013;8(2):e56702.
We are grateful to Dr. Shigeki Sugii, Singapore BioImaging Consortium & Duke-NUS Graduate Medical School, Singapore, for his advice on cellular reprogramming and for providing the retroviral plasmids harboring the Yamanaka factors.
This work was supported by HIR-MoE Grant (UM.C/625/1/HIR/MOHE/ CHAN/03, account no. A000003–50001) to KBC, SKC and TK, and also by a UTARRF grant IPSR/RMC/UTARRF/2013 C2/C11 to MHSY, CJH and KBC.
Availability of data and materials
All data generated in this study are available from the corresponding author on reasonable request.
The miRNA expression profile of HCT-15, SK-CO-1 and four CRC-iPC cell lines has been deposited NCBI GEO with the accession number GSE87280.
Ethics approval and consent to participate
Not applicable since no patients subjects were used in the study.
Consent for publication
Not applicable. All authors read and approved the manuscript.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primer sequences of genes and miRNAs analyzed in RT-PCR or qRT-PCR. (PDF 99 kb)
Table S2. Regulated expression of miRNAs known to modulate cellular reprogramming. (PDF 38 kb)
Table S3. Full list of 102 differentially expressed miRNAs in 4 iPCs vs 2 CRCs. (PDF 220 kb)