Acquisition of cancer stem cell properties in osteosarcoma cells by defined factors
Stem Cell Research & Therapy volume 11, Article number: 429 (2020)
Cancer stem cells (CSCs) are considered to be responsible for tumor initiation, formation, and poor prognosis of cancer patients. However, the rarity of CSCs in clinical samples makes it difficult to elucidate characteristics of CSCs, especially in osteosarcoma (OS). The aim of this study is to verify whether it is possible to generate CSC-like cells by transducing defined factors into an OS cell line.
We retrovirally transduced the Octamer-binding transcription factor 3/4 (OCT3/4), Kruppel-like factor 4 (KLF4), and SRY-box transcription factor 2 (SOX2) genes into the MG-63 human OS cell line (MG-OKS). Parental and GFP-transduced MG-63 cells were used as negative control. We assessed the properties of the generated cells in vitro and in vivo. Multiple comparisons among groups were made using a one-way analysis of variance (ANOVA) followed by post hoc testing with Tukey’s procedure.
MG-OKS cells in vitro exhibited the significantly increased mRNA expression levels of CSC markers (CD24, CD26, and CD133), decreased cell growth, increased chemoresistance and cell migration, and enhanced sphere formation. Notably, MG-OKS cells cultured under osteogenic differentiation conditions showed strongly positive staining for both Alizarin Red S and alkaline phosphatase, indicating osteogenesis of the cells. Gene ontology analysis of microarray data revealed significant upregulation of epidermal-related genes. Tumors derived from MG-OKS cells in vivo were significantly larger than those from other cells in μCT analysis, and immunohistochemical staining showed that Ki-67, osteocalcin, and HIF-1α-positive cells were more frequently detected in the MG-OKS-derived tumors.
In this study, we successfully generated OS CSC-like cells with significantly enhanced CSC properties following transduction of defined factors.
Osteosarcoma (OS) is the most common type of primary malignant bone tumor in children and young adults. OS arises predominantly in the long bones, especially the distal femur, followed by the proximal tibia, and the proximal humerus, with lungs being the most common site of metastasis of the disease . Although the progress of neoadjuvant and adjuvant chemotherapies along with limb-sparing surgery has improved both the prognosis of patients with OS and their limb function, the 5-year survival rate has plateaued at around 70% for patients with non-metastatic lesions, and the rate for the patients with metastatic or recurrent or chemotherapy-resistant lesions has remained below 20% [2, 3]. Therefore, it is important to elucidate the molecular mechanisms underlying the progression, chemoresistance, and metastasis of OS .
To a certain extent, the poor prognosis of patients with OS could be ascribed to the chemoresistance to treatment. Chemoresistance in OS appears to be associated with various mechanisms, including removal of intracellular drug accumulation , enhanced DNA repair , perturbations in signal transduction pathways , autophagy-related chemoresistance , and microRNA (miRNA) dysregulation . Recent studies have suggested that cancer stem cells (CSCs) should play crucial roles in cancer progression and relate to the chemoresistance and metastatic capabilities of cancers . However, it is known that CSCs constitute only a small percentage of tumor cells; therefore, targeting CSCs for treatment or for research, especially in rare cancers, such as OS, presents many challenges.
In previous reports, several methods have been developed to isolate subpopulations with enriched stem cell properties from OS [11, 12]. However, as CSCs comprise only a small population in cancer tissues, sampling limitations remain a major obstacle in the research of CSCs. To this day, an effective (quantitatively and qualitatively) isolation method for CSCs from OS has not been established.
To overcome this challenge, we previously developed a novel method to generate CSC-like cells. We retrovirally transduced a set of defined factors, Octamer-binding transcription factor 3/4 (OCT3/4), Kruppel-like factor 4 (KLF4), and SRY-box transcription factor 2 (SOX2) into human colon  and lung  cancer cells, followed by culturing with a conventional serum-containing medium, but not a human embryonic stem cell medium. Cancer cells transduced with the three factors showed significantly enhanced CSC properties in terms of the expression of marker genes, formation of spheres, chemoresistance, and in vivo tumorigenicity, and were capable of forming tumors that were similar to human colon and lung cancer tissues in terms of both their structural and immunohistological patterns [13, 14]. Moreover, we identified glycogen synthase kinase 3 beta (GSK3B) and interleukin-6 (IL-6) as novel potential therapeutic targets in colon and lung cancers, respectively, using organoids, 3D structures derived from the CSC-like cells [13, 14]. However, this technique has not previously been applied in OS cells, and many aspects of the mechanism remain unexplained.
In the current study, we transduced the three factors into OS cells, resulting in the induction of CSC properties, in terms of the elevated expression of CSC markers, slow cell proliferation, higher chemoresistance, enhanced sphere formation, increased migration and osteogenic ability in vitro, and enhanced tumorigenicity in vivo. Moreover, we performed microarray analysis and identified the differentially expressed genes between the OS CSC-like cells and the control cells, with the goal of identifying the principal genes associated with the molecular mechanisms of progression, metastasis, and chemoresistance, as well as revealing novel effective therapeutic targets in human OS.
The human OS cell lines (MG-63 and NOS-1) were used in this study (RIKEN BRC through the National Bio-Resource Project of the MEXT, Ibaraki, Japan). Cells were cultured in Dulbecco’s modified Eagle medium (Sigma-Aldrich, St Louis, MO, USA) supplemented with 10% fetal bovine serum (FBS; Sigma-Aldrich), 100 U/mL penicillin (Sigma-Aldrich), and 100 μg/mL streptomycin (Sigma-Aldrich) and were maintained in a humidified atmosphere with 5% CO2 at 37 °C. Plat-A packaging cells (Cosmo Bio Co., LTD, Tokyo, Japan) were used for the production of the retrovirus. In Plat-A cultures, 1 mg/mL of puromycin (Life Technologies, Grand Island, NY, USA) and 10 mg/mL of blasticidin (Life Technologies) were added.
For gene transduction, we employed a modified version of the retroviral infection method previously described [15, 16]. We used the polycistronic retroviral vector designed to encode OCT3/4, KLF4, and SOX2 in the pMX retroviral vector (pMX-OKS) (Supplemental Fig. S1A). Concomitantly, pMX with GFP (pMX-GFP) was used as a control vector. To produce retroviral particles, Plat-A cells (in DMEM with 10% FBS without antibiotics) were transfected with the retroviral vector (pMX-OKS or pMX-GFP) using the FuGene HD transfection reagent (Promega, Madison, WI, USA) following the manufacturer’s instructions. The medium was replaced at 24 h after transduction, and the retrovirus-containing supernatant was harvested at 48 h after transduction. The supernatant was filtered through a 0.45-μm pore-size syringe filter (Sartorius Stedim Biotech, Goettingen, Germany). Infection of both cell lines with the retrovirus was conducted in the presence of 4 μg/mL polybrene (Nacalai Tesque, Kyoto, Japan) for 24 h. Non-transfected cells (MG-parental and NOS-parental), as well as cells transfected with GFP (MG-GFP and NOS-GFP), were used as control. All gene transduction procedures were performed in accordance with the National Institutes of Health Guidelines, and the study protocol was approved by the Kobe University Institutional Committee (Permission no. 30-18).
RNA isolation and real-time quantitative reverse-transcription polymerase chain reaction (RT-qPCR)
Total RNA was extracted from cultured cells and tumor tissues using an RNeasy mini kit (Qiagen, Valencia, CA, USA), and qPCR reactions were performed with the SYBR Green master mix reagent (Applied Biosystems, Foster City, CA, USA) on the ABI prism 7500 sequence-detection system (Applied Biosystems) according to the manufacturer’s instructions. Relative mRNA expressions of transduced genes (OCT3/4, KLF4, and SOX2), previously reported markers (CD24, CD26, and CD133), chemoresistance-related gene (ABCB1), and osteogenic differentiation-related genes (osteocalcin, BMP2, BMP4, and BMP6) were assessed. The details are given in the Supplemental Data.
Cell proliferation and cell migration assay
We assessed cell proliferative activities by performing cell counting using a hemocytometer for measurement of total cell count and WST-8 assays using a colorimetric Cell Counting Kit-8 (CCK-8; Dojindo Inc., Kumamoto, Japan), whereby a formazan dye color intensity is directly proportional to the viable number of cells. The migration ability of cells was evaluated using a scratch wound healing assay . The details are given in the Supplemental Data.
Analysis of chemoresistance to doxorubicin
To assess the chemoresistance to doxorubicin (DOX) in cells with and without gene transduction, the viability of cells under exposure to DOX was measured by performing the WST-8 assays using CCK-8 (Dojindo Inc.). A total of 5 × 103 cells at 10 days after retroviral transduction were seeded in 96-well plates, and after 24 h, the medium was replaced with DMEM containing 0 (as control), 0.3, or 30 μM of DOX (Sigma-Aldrich). The WST-8 assays were performed after 48 h of incubation, and the relative viability of cells was calculated. The expression level of the ATP-binding cassette subfamily B member 1 (ABCB1; also known as MDR1) gene in the DOX-treated cells was also examined by qPCR analysis.
Induction of osteogenic differentiation and evaluation of osteogenic ability in vitro
Osteogenic differentiation was evaluated by alkaline phosphatase (ALP) staining and by assessment of calcium in deposited minerals using Alizarin Red S staining. A total of 1 × 105 cells at 7 days after retroviral transduction were seeded in 24-well plates, following which, the medium was replaced with the osteogenesis culture kit (Cosmo Bio). After 3 weeks of incubation, staining of ALP and Alizarin Red S was performed using an ALP staining kit (Cosmo Bio) and a Calcified Nodule Staining kit (Cosmo Bio) following the manufacturer’s instructions, respectively. Concomitantly, the mRNA expressions of osteocalcin and bone morphogenetic protein (BMP) family members (BMP2, BMP4, and BMP6), which stimulate osteoblast differentiation [18, 19], were also assessed in the cells by qPCR analysis.
Sphere formation assay
Cells (1 × 105 cells/well) at 10 days after retroviral transduction were seeded to Ultra-Low Attachment Surface 6 well plates (Corning Inc., Corning, NY, USA) in a serum-free DMEM medium containing 10 ng/mL basic fibroblast growth factor (bFGF; Wako, Osaka, Japan), 10 mg/mL human insulin (CSTI, Miyagi, Japan), 100 mg/mL human transferrin (Roche, Basel, Switzerland), and 100 mg/mL bovine serum albumin (BSA; Nacalai Tesque) and incubated at 37 °C in a 5% CO2 incubator for 10 days . Tumor spheroids were manually counted under an inverted phase contrast microscope (BZ-X710 Microscope and BZ-X Viewer, BZ-X Analyzer imaging system, Keyence, Osaka, Japan). All morphometric studies were performed by two examiners blinded to treatment conditions.
Immunoblot analysis was performed to assess the expression of epithelial-mesenchymal transition (EMT) markers (E-cadherin and vimentin) and previously reported marker genes of CSCs (CD24, CD26, and CD133). The details are given in the Supplemental Data.
Male BALB/c nude mice (5 weeks old) were purchased from CLEA Japan Inc. (Tokyo, Japan) and maintained in a facility under specific pathogen-free conditions. Mice were fed with pathogen-free laboratory chow and allowed free access to autoclaved water in an air-conditioned room under a 12-h light/dark cycle. A total of 2 × 106 cells (MG-parental n = 5, MG-OKS n = 6, or MG-GFP n = 6) in 200 μL of serum-free PBS were subcutaneously injected into the dorsal flank of mice. Consecutively, 8 weeks after cell transplantation, the volume of generated tumors was calculated by the formula 0.5 × (length) × (width)2, and all tumors were excised. All morphometric studies were performed by 2 examiners blinded to treatment conditions. The μCT Scanner (R_mCT; Rigaku Mechatronics, Tokyo, Japan) was used for microarchitectural analysis. Parameters used for the scans were as follows: FOV 30 mm, 90 kV tube voltage, and 160 μA tube current.
Hematoxylin-eosin staining and immunofluorescence staining on frozen sections
Mice were euthanized for histological analyses at 8 weeks after cell transplantation. For analyses, 10-μm-thick frozen sections were cut, washed 5 times with PBS, and stained with hematoxylin-eosin. Sections were microscopically evaluated to confirm the presence of tumor cells. To evaluate the proliferation ability, hypoxia environment, and the osteogenic markers in the xenograft tumor tissues, immunofluorescence staining was performed under a fluorescence microscope (BZ-X710 Microscope and BZ-X Viewer, BZ-X Analyzer imaging system, Keyence) with an anti-human Ki-67 rabbit polyclonal antibody (1:50; Novus Bio, Littleton, CO, USA, catalog number: NB500-170), an anti-human hypoxia-inducible factor 1α (HIF-1α) rabbit polyclonal antibody (1:100; Abcam, Cambridge, MA, USA, catalog number: ab82832), and an anti-human osteocalcin mouse monoclonal antibody (1:100; Takara Bio, Shiga, Japan, catalog number: M184), respectively. The secondary antibodies used were the Alexa-Fluor 594-conjugated chicken anti-rabbit IgG (1:100; Invitrogen, Carlsbad, CA, USA, catalog number: A-21442) for both Ki-67 and HIF-1α and the Alexa-Fluor 594-conjugated goat anti-mouse antibody (1:100; Molecular Probes, Eugene, OR, USA, catalog number: R37121) for osteocalcin, while nuclear staining was performed using the 4′,6-diamino-2-phenylindole (DAPI) solution (1:5000; Sigma-Aldrich). The positivities were quantified using ImageJ software (National Institute of Health, Bethesda, MD, USA) in four randomly selected fields. All morphometric studies were performed by two examiners blinded to treatment conditions.
Total RNA was extracted from the cells 11 days after gene transduction using an RNeasy mini kit (Qiagen). For RNA quality control, concentration and purity were checked with a NanoDrop 1000 (Thermo Fisher Scientific, Waltham, MA, USA) and integrity was evaluated with a TapeStation system (Agilent Technologies, Santa Clara, CA, USA). Gene expression profiling in the MG-parental, the MG-OKS, and the MG-GFP cells was performed using the SuperPrint G3 Human Gene Expression 8 × 60 K Ver3 Microarray (Agilent Technologies) following the manufacturer’s protocol.
The data were analyzed using the GeneSpring 13.1.1 software program (Agilent Technologies) as follows. Threshold raw signals were set to 1.0, and signal intensities were used with log base 2 transformation. The 75th percentile normalization method was used as the normalized algorithm (http://genespringsupport.com/faq/normalization). The default flag setting was used. The number of detected probes in at least one of three cell groups was 36,323, and the probes were used for further analysis. Gene ontology (GO) analysis of the obtained microarray data was also performed using GeneSpring GX 13.1.1 software (Agilent Technologies). Z-scores were calculated by subtracting the overall averaged gene intensity from the normalized intensity of each gene and by dividing this result by the standard deviation (SD) of all measured intensities, according to the following equation: Z-score = (intensity − mean intensity)/SD .
Microarray data have been deposited in NCBI GEO under accession number GSE143556 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE143556).
Obtained results were statistically analyzed using a software package (Graph Pad Prism™, MDF Software, La Jolla, CA, USA). All values were expressed as mean ± SEM (standard error of the mean). Multiple comparisons among groups were made using a one-way analysis of variance (ANOVA) followed by post hoc testing with Tukey’s procedure. A probability value (P) of less than 0.05 was considered to denote statistical significance.
Transduction of defined factors enhanced CSC properties in OS cells
The efficiency of retroviral infection was monitored with the expression of GFP in control cells (MG-GFP) using fluorescence microscopy; a strong expression of GFP was accordingly observed (Supplemental Fig. S1B). qPCR analyses revealed that the total transcript levels of OCT3/4, KLF4, and SOX2 were significantly elevated in the MG-OKS cells (P < 0.05, n = 3) (Supplemental Fig. S1C). To assess the stem cell property of the gene-transfected and parental cells, we evaluated the expression levels of previously reported marker genes of CSCs in several cancers, such as CD24 , CD26, and CD133  by qPCR analysis. The mRNA expressions of all the genes were significantly increased in the MG-OKS cells compared with those of other cell populations (P < 0.05, n = 3) (Fig. 1a). We also evaluated alterations in the markers at the protein level by immunoblot analyses (Supplemental Fig. S1D). Regarding the morphology of the cells, an elongated cell body with invasive processes was observed in the MG-OKS cells, whereas the MG-parental and the MG-GFP cells predominantly consisted of spindle-shaped cells (Fig. 1b). In addition, the EMT, a biological process that assumes stem-like features, was also assessed by immunoblotting. Accordingly, the MG-OKS cells were characterized by a decreased expression of E-cadherin and an increased expression of vimentin (Fig. 1c).
Decreased proliferation and elevated migration were observed in the MG-OKS cells
The number of the MG-OKS cells was lower than that of both the MG-parental and the MG-GFP cells at 22 days after gene transduction, but this was not statistically significant (Fig. 2a). Likewise, cell proliferation assays also revealed that the proliferation of cells was significantly decreased in the MG-OKS cells compared with other cells at 72 h after seeding (P < 0.05, n = 3) (Fig. 2b). In wound healing assays, the speed of wound healing was obviously increased in the MG-OKS group (Fig. 2c) and was significantly elevated in the MG-OKS cells compared with other cells (P < 0.05, n = 5) (Fig. 2d).
Transduction of the factors enhanced chemoresistance to DOX in OS cells
To examine the effect of the transduction of OCT3/4, KLF4, and SOX2 genes on the chemoresistance to DOX in MG-63 OS cells, we compared the viability of cells after treatment with DOX using the WST-8 assays. There was no significant difference in cell viability following treatment with a low concentration of DOX (0.3 μM) (P = 0.3, n = 3) (Fig. 2e). However, when a higher concentration of DOX (30 μM) was administered, the viability of MG-OKS cells was approximately 20% higher than that in both the MG-parental and MG-GFP cells (P < 0.05, n = 4) (Fig. 2e). In addition, the mRNA expression level of ATP-binding cassette subfamily B member 1 (ABCB1; also known as MDR1) gene, which encodes the membrane drug transporter P-glycoprotein and has been reported as a well-described mechanism of resistance to drugs , was approximately 3 times higher in the MG-OKS cells compared with other cells (ABCB1: MG-parental set to 1, MG-OKS 3.4 ± 0.3, MG-GFP 1.1 ± 0.3, P < 0.05; n = 3) (Fig. 2f).
Confirmation of the reproducibility with another osteosarcoma cell line
We confirmed the reproducibility of our results using another osteosarcoma cell line (NOS-1). Several in vitro experiments, including measurement of the mRNA expression of markers (CD24, CD26, and CD133), cell morphology, cell proliferation rate, migration ability, and chemoresistance to DOX, were performed. We obtained similar results as in MG-63 cells. The mRNA expression of CD24 was significantly increased in the NOS-OKS cells compared to that in other cell populations (P < 0.05, n = 3). The mRNA expression of CD26 and CD133 were significantly increased in NOS-OKS cells compared with NOS-parental cells (P < 0.05, n = 4) (Fig. S2A). For morphological features, NOS-OKS cells showed similar changes as in MG-OKS cells (Fig. S2B). Cell proliferation assays revealed significantly decreased proliferation in NOS-OKS cells (P < 0.05, n = 4) (Fig. S3A); in wound healing assays, the rate of wound healing was significantly increased in NOS-OKS cells compared with other cells (P < 0.05, n = 5) (Fig. S3B, C). Chemoresistance to DOX and the mRNA expression level of ABCB1 was significantly higher than those in both NOS-parental and NOS-GFP cells, which were similar to those in the experiments using MG-63 (P < 0.05, n = 3) (Fig. S3D, E).
Osteogenic ability was elevated in the MG-OKS cells
Osteosarcoma has been generally defined as a high-grade tumor of malignant mesenchymal cells with osteoid formation . Hence, we assessed the ability of osteogenic differentiation in the cells. After culturing with the osteogenesis culture kit for 3 weeks, strong staining of both ALP and Alizarin Red S was observed in the MG-OKS cells but not in the other cells (Fig. 3a). Concomitantly, the mRNA expressions of both osteocalcin, a marker of mature osteoblasts , and BMPs (BMP-2, BMP-4, and BMP-6), which are known to stimulate osteoblast differentiation [18, 19], were significantly elevated in the MG-OKS cells compared with other cells (osteocalcin: MG-parental set to 1, MG-OKS 6.7 ± 3.0, MG-GFP 1.0 ± 0.05, P < 0.05, n = 5; BMP2: MG-parental set to 1, MG-OKS 23.3 ± 9.9, MG-GFP 1.1 ± 0.2, P < 0.05, n = 5; BMP4: MG-parental set to 1, MG-OKS 3.5 ± 0.3, MG-GFP 1.1 ± 0.1, P < 0.05, n = 3; BMP6: MG-parental set to 1, MG-OKS 4.0 ± 0.6, MG-GFP 0.8 ± 0.1, P < 0.05, n = 3) (Fig. 3b).
The sphere formation ability was elevated in the MG-OKS cells
As previously reported, it is known that CSC has the ability to form sphere in anchorage-independent and serum-starved conditions [25, 26]. To examine the sphere formation ability of these cells, we performed sphere formation assay (Fig. 4a). In the MG-OKS cultures, the number of spheres was significantly increased compared to that of other cells (Fig. 4b).
In vivo OS tumorigenicity was accelerated in the MG-OKS cells
To examine their in vivo tumorigenicity, we subcutaneously transplanted parental and gene-transduced OS cells into the dorsal area of immunodeficient nude mice. In mice with MG-OKS xenografts, 100% engraftment and a significantly larger volume of tumors were observed at 8 weeks after cell transplantation compared with those carrying xenografts derived from MG-parental and MG-GFP (P < 0.05) (MG-parental, 0.1 ± 0.08 cm3, MG-OKS 0.8 ± 0.1 cm3, MG-GFP 0.04 ± 0.02 cm3, P < 0.001; n = 5 in MG-parental, n = 6 in MG-OKS, and MG-GFP) (Fig. 4c, d). Immunofluorescence staining revealed a strong staining of Ki-67, a marker of proliferative potential, in tumors derived from MG-OKS cells (Fig. 5a). We also observed increased positive staining of both osteocalcin (Fig. 5b) and HIF-1α, which permits hypoxic tumor cells to upregulate proteins that promote their survival and increase their aggressiveness , in tumors derived from MG-OKS cells (Fig. 5c).
Microarray and data analyses
We used microarray analysis to identify changes in gene expression pattern among the MG-parental, MG-OKS, and MG-GFP cells. The effects were visualized after filtering by applying the gene list to a scatter plot with a fold-change level of 10 as a threshold (genes upregulated, > 10-fold change) (Fig. 6a). In the MG-OKS cells, 269 and 302 genes were upregulated above a 10-fold change compared with the MG-parental and the MG-GFP cells, respectively. Among these upregulated genes, 219 genes were commonly upregulated as compared with both the MG-parental and the MG-GFP cells (Fig. 6b). Of these 219 genes, the top 10 differentially expressed genes are shown in Table 1. In contrast, 75 and 86 genes showed over 10-fold downregulation in the MG-OKS cells compared with MG-parental and MG-GFP cells, respectively (Fig. S4A). Among these downregulated genes, only 1 gene was commonly downregulated (Fig. S4B). The gene was C3orf49, which function had been still unclear; therefore, we focused on the upregulated 219 genes. GO analysis showed that the 219 upregulated genes were mainly enriched in “biological processes (BP),” including epidermis development (GO: 0008544), skin development (GO: 0043588), epidermal cell differentiation (GO: 0009913), keratinocyte differentiation (GO: 0030216), keratinization (GO: 0031424), animal organ development (GO: 0048513), tissue development (GO: 0009888), and followed in “cellular component (CC),” including cornified envelope (GO: 0001533), extracellular region (GO: 0005576), and extracellular region part (GO: 0044421) (Fig. 6c).
In our series of experiments, we have succeeded in obtaining for the first time CSC-like cells from a human OS cell line by transducing them with the OCT3/4, KLF4, and SOX2 genes. Recently, CSCs, a subpopulation of tumor cells, have been characterized and emerged as a major topic of interest in the cancer research field . Isolation and subsequent studies of CSCs from various types of cancers indicated CSCs as one of the crucial causes of conventional treatment failure. In consequence, targeting CSCs should be a promising perspective for the development of more effective anticancer therapies. Especially, CSCs are attracting a great deal of attention in the research subjects of rare cancers, such as OS. CSCs are considered to play a key role for the poor prognosis of patients with OS due to therapeutic resistance, and the possible reason for the development of recurrence or metastasis . However, difficulties in acquiring proper research samples by conventional methods have interfered with the elucidation of OS CSC biology and the development of novel therapies that could target OS CSCs. Several previous studies have attempted to isolate OS CSCs from clinical specimens [13,14,15,16,17]. First, CSCs were collected by using their self-renewal capacity to form spherical colonies, named “sarcospheres,” under serum-free conditions [25, 26]. The formation of sarcospheres was shown to be further improved by cultivating under the hypoxic conditions of the tumor microenvironment . Then, OS CSCs were isolated by sorting cells according to the expression of previous reported markers associated with cancer stem cells, such as CD133 or CD117, in combination with Stro-1 . Other methods used to isolate OS CSCs included the identification of a “side population” of cells able to efflux Hoechst dyes , with/without the expression of surface markers, such as CD117 or STRO-1 ; the sorting of cells with high ALDH1 activity; and the tracking of subpopulations of OCT-4-positive cells . In this study, we successfully generated OS CSC-like cells by transducing defined factors into human OS cell lines. Unlike the previous methods, our novel technique could make it feasible to obtain OS CSCs more effectively and abundantly.
CSCs are predominantly present in the dormant state of the cell cycle and undergo proliferation in response to physiological cell stimuli; the quiescent state is an important factor in chemoresistance . Our data revealed a slow cell proliferation rate in vitro and high tumorigenicity in vivo (larger tumor volume with more Ki67+ tumor cells) in MG-OKS xenografts. Subcutaneous transplantation of the cells may have led to cell cycle re-entry from the dormant state as external stimuli. Dembinski et al. reported that a subpopulation of slow cycling cells in the pancreas adenocarcinoma cell lines exhibited increased tumorigenic and invasive potentials in vivo ; our results are consistent with their findings. We assessed altered expression of genes, related to slow cell growth in MG-OKS cells using microarray data. Our data revealed increased expression of the cyclin-dependent kinase (CDK) inhibitor 1A (CDKN1A) and retinoblastoma 1 (RB1) and decreased expressions of CDK4, cyclin D (CCND), proliferating cell nuclear antigen (PCNA), and MYC (Fig. 6d). CDKN1A inhibits the cell cycle during G1 phase . CCND binds to and activates CDK4, and then the CCND-CDK4 complex inhibits members of the RB protein family including RB1 and regulates the cell cycle during G1/S transition . Therefore, slow cell growth in MG-OKS cells was likely caused by inhibition of the cell cycle at the G1 phase via increased CDKN1A expression. Decreased expression of PCNA which accumulate in the G1 phase of the cell cycle, reaching maximum expression in S phase, may also inhibit the G1/S transition . MYC is one of the most commonly activated oncogenes implicated in the pathogenesis of human cancers and induces tumorigenesis by evading multiple tumor-suppressing checkpoint mechanisms including proliferative arrest . Our results showed decreased MYC expression, which likely contribute to the slow cell proliferation of MG-OKS cells in vitro.
As previously reported, isolated OS CSCs have been associated with the metastasis and chemoresistance of the disease [29, 32, 40]. And, this increased chemoresistance of CSCs has been associated with their increased ability for DNA repair, their frequent quiescent state of low proliferation rate , suppression of the apoptotic signaling, and increased levels of lysosomal activity due to the overexpression of vacuolar ATPases . Additionally, chemoresistance of CSCs might be strongly related to a gain in the drug efflux capability due to the overexpression of the ATP-binding cassette (ABC) family transporters, especially ABCG2  or the ABCB1 multidrug resistance P-glycoprotein [44,45,46]. In the current study, our generated OS CSC-like cells showed enhanced chemoresistance to doxorubicin, the most important anticancer drug in the chemotherapy of OS , associated with an increased mRNA expression of ABCB1. Therefore, it is suggested that OS CSCs have specific characteristics making them more resistant to chemotherapies. The present results also indicated that the generated OS CSC-like cells showed elevated ability for cell migration, which is associated with both the progression and metastasis of OS cells.
Osteosarcoma is a neoplasm derived from the primitive bone-forming mesenchyme. Pathologic feature is characterized by the production of osteoid and new bone by spindle-shaped tumor cells. However, the MG-63 cell line, a human OS-derived cell line, has impaired the capability to synthesize a correct extracellular bone matrix [48, 49] and has a low ALP activity . Interestingly, in the current study, the CSC-like cells generated from MG-63 cells following transduction of defined factors were able to induce osteogenic differentiation in vitro, whereas bone formation was not observed in vivo. There was a report that bone formation was not observed following simple implantation of an OS cell line into muscle tissue of mice; however, formation of the bone was noted when cells were pretreated with dexamethasone before implantation . One of the possible explanations for our in vivo results should be that bone formation did not occur because steroids derived from nude mice could not act on human-derived OS cells. Another report suggested that the γ-carboxylation of osteocalcin by γ-carboxylase increased its affinity for hydroxyapatite, the mineral component of the bone extracellular matrix . The absence of bone formation despite the increased expression of osteocalcin both in vitro and in vivo suggested a problem with this carboxylation process.
The EMT is a biological process including various phases through which epithelial cells lose their characteristics with intercellular adhesion and proliferative capacity while acquiring mesenchymal characteristics with increased migratory and invasive potencies. This process has been reported in several cancers of epithelial origin and is important to metastatic tumor cell progression and chemoresistance . In contrast, the precise role of EMT-related processes of tumor mesenchymal origin, such as OS, is still mostly unknown . There are reports that activation of the EMT and reactivation from cellular dormancy could cause tumor metastasis and/or relapse [54, 55]. In this study, we reported that the MG-OKS cells showed a decreased expression of E-cadherin with an increased expression of vimentin, suggesting activation of EMT. Yang et al. reported that hypoxia or HIF-1α overexpression induces various EMT regulators, such as SNAI1, SNAI2, and TWIST, and that TWIST plays critical roles in the acquisition of EMT and metastatic phenotypes including an increased migration ability . Our microarray analysis results are consistent with this previous report, as MG-OKS cells showed greater migration ability with increased expression of HIF-1α and TWIST and decreased expression of E-cadherin (Fig. 6e).
We also surveyed the differentially expressed genes between MG-OKS and other cells by microarray analysis, and numerous genes showing significant differential expression were identified. Interestingly, although OS is a tumor of mesenchymal origin, upregulated genes in MG-OKS cells included epithelial-related genes, such as small proline-rich protein 2A (SPRR2A), SPRR1A, keratin 6A (KRT6A), and keratin 6B (KRT6B). There have been several reports on the relationships of keratin and SPRRs with cancer . Zachary et al. reported that urothelial cell carcinomas expressed high levels of the basal keratins (KRT6A, KRT6B, KRT6C, KRT14, and KRT16), which are not normally expressed in the urothelium, and expression of these genes was an indicator of poor prognosis . Kim et al. reported that SPRR3 promoted the proliferation of breast cancer cells by enhancing p53 degradation via the AKT and MAPK pathways . Vos et al. reported that a subset of genes, such as SPRR2A, KRT6A, S100A2, and others, which are related and form an effective barrier against external (chemical and physical) stimuli, were expressed in bronchial epithelial cells and keratinocytes . Our microarray data revealed similar changes in these genes (Fig. 6f). However, few studies have reported on the relationship between stromal cell-derived OS and epithelial markers, and how these associated genes forming a barrier in the epithelial system affect OS CSC-like cells remains unclear. So further study of these relationships might lead to the discovery of prognostic markers and new therapeutic strategies to target OS CSCs.
There have been several limitations in this study. First, to avoid heterogeneity, we designed a polycistronic retroviral vector in which the frames of OCT3/4, KLF4, and SOX2 were linked via T2A sequences. Although these genes were individually reported to be correlated with malignant behavior and poor prognosis in OS [33, 60, 61], it has been unclear which gene contributes the most toward gaining of CSC properties. To solve this problem, it will be necessary to use 3 viral vectors independently carrying OCT3/4, KLF4, or SOX2 for the generation of CSCs in future studies. Second, it is important to mention that generated OS CSC-like cells might be heterogeneous, with different subpopulations carrying different genetic alterations existing within a tumor. Moreover, these subpopulations might be highly dynamic and able of processes of dedifferentiation and phenotype switching which may render the CSCs resistant to a specific CSC therapy . In this regard, to overcome this limitation, we need to find appropriate culture conditions to maintain CSC properties and future therapies should combine different treatments to target both non-CSCs and CSCs. Finally, similar to normal stem cells, specific signaling or tumor niches might have responsibility for the regulation of OS CSCs . OS CSC niche is reported to be composed of three types of niches within the bone microenvironment where these signaling pathways are particularly active: the perivascular niche, the hypoxic niche, and the endosteal niche . In addition, Wnt and Notch signaling pathways have also been reported to play a role in the progression of OS . In OS cell lines, inhibition of the Wnt and Notch pathways enhanced the sensitization to chemotherapy . Besides this, activation of the Wnt/β-catenin signaling has been reported in the OS CSC population and has been related to the overexpression of SOX2 and cancerogenesis . Therefore, there is a need for developing and testing therapeutic strategies directed against specific signaling or tumor niches, or both. Although it is true that we need to examine CSC-specific treatments that would target multiple pathways altered in different subsets of CSCs within the tumor in future studies, our novel method could be useful in shedding light on the undisclosed molecular mechanisms and facilitating development of new therapeutic approaches for the management and treatment of human OS.
To the best of our knowledge, this is the first study to generate CSC-like cells in terms of the elevated expression of CSC markers, slow cell proliferation, higher chemoresistance, increased migration, and osteogenic abilities in vitro, as well as tumorigenicity in vivo, from an established OS cell line by exogenous expression of OCT3/4, KLF4, and SOX2 genes. Although further investigations are needed to elucidate the properties of these CSC-like cells, this method could be useful in solving the quantitative problem of CSCs and in elucidating the molecular mechanisms and identifying novel effective therapeutic targets in human OS.
Availability of data and materials
Microarray data have been deposited in NCBI GEO under accession number GSE143556 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE143556).
ATP-binding cassette subfamily B member 1
Bone morphogenetic protein
Basic fibroblast growth factor
Bovine serum albumin
Cell Counting Kit-8
Cyclin-dependent kinase 1A
Cancer stem cell
Gap junction beta-4 protein
Glycogen synthase kinase 3 beta
Hypoxia-inducible factor 1α
Insulin-like growth factor 2
Kruppel-like factor 4
Neutrophil cytosol factor 2
Octamer-binding transcription factor 3/4
Proliferating cell nuclear antigen
Real-time quantitative reverse-transcription polymerase chain reaction
SRY-box transcription factor 2
Small proline-rich protein
Wittig JC, Bickels J, Priebat D, et al. Osteosarcoma: a multidisciplinary approach to diagnosis and treatment [in eng]. Am Fam Physician. 2002;65(6):1123–32.
Ando K, Heymann MF, Stresing V, et al. Current therapeutic strategies and novel approaches in osteosarcoma [in eng]. Cancers. 2013;5(2):591–616.
ESMO/European Sarcoma Network Working Group. Bone sarcomas: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up [in eng]. Ann Oncol. 2014;25(Suppl 3):iii113–23.
He H, Ni J, Huang J. Molecular mechanisms of chemoresistance in osteosarcoma (review) [in eng]. Oncol Lett. 2014;7(5):1352–62.
Park YB, Kim HS, Oh JH, et al. The co-expression of p53 protein and P-glycoprotein is correlated to a poor prognosis in osteosarcoma [in eng]. Int Orthop. 2001;24(6):307–10.
Wang D, Luo M, Kelley MR. Human apurinic endonuclease 1 (APE1) expression and prognostic significance in osteosarcoma: enhanced sensitivity of osteosarcoma to DNA damaging agents using silencing RNA APE1 expression inhibition [in eng]. Mol Cancer Ther. 2004;3(6):679–86.
Gazitt Y, Kolaparthi V, Moncada K, et al. Targeted therapy of human osteosarcoma with 17AAG or rapamycin: characterization of induced apoptosis and inhibition of mTOR and Akt/MAPK/Wnt pathways [in eng]. Int J Oncol. 2009;34(2):551–61.
Lambert LA, Qiao N, Hunt KK, et al. Autophagy: a novel mechanism of synergistic cytotoxicity between doxorubicin and roscovitine in a sarcoma model [in eng]. Cancer Res. 2008;68(19):7966–74.
Song B, Wang Y, Xi Y, et al. Mechanism of chemoresistance mediated by miR-140 in human osteosarcoma and colon cancer cells [in eng]. Oncogene. 2009;28(46):4065–74.
Shiozawa Y, Nie B, Pienta KJ, et al. Cancer stem cells and their role in metastasis [in eng]. Pharmacol Ther. 2013;138(2):285–93.
Dela Cruz FS. Cancer stem cells in pediatric sarcomas [in eng]. Front Oncol. 2013;3:168.
Tirino V, Desiderio V, Paino F, et al. Cancer stem cells in solid tumors: an overview and new approaches for their isolation and characterization [in eng]. FASEB J. 2013;27(1):13–24.
Oshima N, Yamada Y, Nagayama S, et al. Induction of cancer stem cell properties in colon cancer cells by defined factors [in eng]. PLoS One. 2014;9(7):e101735.
Ogawa H, Koyanagi-Aoi M, Otani K, et al. Interleukin-6 blockade attenuates lung cancer tissue construction integrated by cancer stem cells. Sci Rep. 2017;7(1):12317.
Takahashi K, Yamanaka S. Induction of pluripotent stem cells from mouse embryonic and adult fibroblast cultures by defined factors [in eng]. Cell. 2006;126(4):663–76.
Takahashi K, Tanabe K, Ohnuki M, et al. Induction of pluripotent stem cells from adult human fibroblasts by defined factors [in eng]. Cell. 2007;131(5):861–72.
Xie Y, Sun W, Deng Z, et al. MiR-302b suppresses osteosarcoma cell migration and invasion by targeting Runx2 [in eng]. Sci Rep. 2017;7(1):13388.
Hee CK, Nicoll SB. Endogenous bone morphogenetic proteins mediate 1alpha, 25-dihydroxyvitamin D (3)-induced expression of osteoblast differentiation markers in human dermal fibroblasts [in eng]. J Orthop Res. 2009;27(2):162–8.
Matsubara T, Kida K, Yamaguchi A, et al. BMP2 regulates Osterix through Msx2 and Runx2 during osteoblast differentiation [in eng]. J Biol Chem. 2008;283(43):29119–25.
Cheadle C, Vawter MP, Freed WJ, et al. Analysis of microarray data using Z score transformation. J Mol Diagn. 2003;5(2):73–81.
Salaria S, Means A, Revetta F, et al. Expression of CD24, a stem cell marker, in pancreatic and small intestinal neuroendocrine tumors [in eng]. Am J Clin Pathol. 2015;144(4):642–8.
Wu C, Xie Y, Gao F, et al. Lgr5 expression as stem cell marker in human gastric gland and its relatedness with other putative cancer stem cell markers [in eng]. Gene. 2013;525(1):18–25.
Vaidyanathan A, Sawers L, Gannon AL, et al. ABCB1 (MDR1) induction defines a common resistance mechanism in paclitaxel- and olaparib-resistant ovarian cancer cells [in eng]. Br J Cancer. 2016;115(4):431–41.
Kondo A, Otsuka T, Kato K, et al. AMP-activated protein kinase regulates thyroid hormone-stimulated osteocalcin synthesis in osteoblasts [in eng]. Int J Mol Med. 2013;31(6):1457–62.
Gibbs CP, Kukekov VG, Reith JD, et al. Stem-like cells in bone sarcomas: implications for tumorigenesis [in eng]. Neoplasia (New York, NY). 2005;7(11):967–76.
Fujii H, Honoki K, Tsujiuchi T, et al. Sphere-forming stem-like cell populations with drug resistance in human sarcoma cell lines [in eng]. Int J Oncol. 2009;34(5):1381–6.
Roncuzzi L, Pancotti F, Baldini N. Involvement of HIF-1alpha activation in the doxorubicin resistance of human osteosarcoma cells [in eng]. Oncol Rep. 2014;32(1):389–94.
Eun K, Ham SW, Kim H. Cancer stem cell heterogeneity: origin and new perspectives on CSC targeting [in eng]. BMB Rep. 2017;50(3):117–25.
Salerno M, Avnet S, Bonuccelli G, et al. Sphere-forming cell subsets with cancer stem cell properties in human musculoskeletal sarcomas [in eng]. Int J Oncol. 2013;43(1):95–102.
Tirino V, Desiderio V, d'Aquino R, et al. Detection and characterization of CD133+ cancer stem cells in human solid tumours [in eng]. PLoS One. 2008;3(10):e3469.
Honoki K, Fujii H, Kubo A, et al. Possible involvement of stem-like populations with elevated ALDH1 in sarcomas for chemotherapeutic drug resistance [in eng]. Oncol Rep. 2010;24(2):501–5.
Adhikari AS, Agarwal N, Wood BM, et al. CD117 and Stro-1 identify osteosarcoma tumor-initiating cells associated with metastasis and drug resistance [in eng]. Cancer Res. 2010;70(11):4602–12.
Levings PP, McGarry SV, Currie TP, et al. Expression of an exogenous human Oct-4 promoter identifies tumor-initiating cells in osteosarcoma [in eng]. Cancer Res. 2009;69(14):5648–55.
Chen W, Dong J, Haiech J, et al. Cancer stem cell quiescence and plasticity as major challenges in cancer therapy [in eng]. Stem Cells Int. 2016;2016:1740936.
Dembinski JL, Krauss S. Characterization and functional analysis of a slow cycling stem cell-like subpopulation in pancreas adenocarcinoma [in eng]. Clin Exp Metastasis. 2009;26(7):611–23.
Engeland K. Cell cycle arrest through indirect transcriptional repression by p53: I have a DREAM. Cell Death Differ. 2018;25(1):114–32.
Bonelli P, Tuccillo FM, Borrelli A, et al. CDK/CCN and CDKI alterations for cancer prognosis and therapeutic predictivity. Biomed Res Int. 2014;2014:361020.
Guzińska-Ustymowicz K, Stepień E, Kemona A. MCM-2, Ki-67 and PCNA protein expressions in pT3G2 colorectal cancer indicated lymph node involvement [in eng]. Anticancer Res. 2008;28(1b):451–7.
Gabay M, Li Y, Felsher DW. MYC activation is a hallmark of cancer initiation and maintenance [in eng]. Cold Spring Harb Perspect Med. 2014;4(6):a014241.
Wang L, Park P, Zhang H, et al. Prospective identification of tumorigenic osteosarcoma cancer stem cells in OS99-1 cells based on high aldehyde dehydrogenase activity [in eng]. Int J Cancer. 2011;128(2):294–303.
Nunes T, Hamdan D. Targeting cancer stem cells to overcome chemoresistance. Int J Mol Sci. 2018;19(12):4036.
Salerno M, Avnet S, Bonuccelli G, et al. Impairment of lysosomal activity as a therapeutic modality targeting cancer stem cells of embryonal rhabdomyosarcoma cell line RD [in eng]. PLoS One. 2014;9(10):e110340.
Bleau AM, Hambardzumyan D, Ozawa T, et al. PTEN/PI3K/Akt pathway regulates the side population phenotype and ABCG2 activity in glioma tumor stem-like cells [in eng]. Cell Stem Cell. 2009;4(3):226–35.
Martins-Neves SR, Lopes AO, do Carmo A, et al. Therapeutic implications of an enriched cancer stem-like cell population in a human osteosarcoma cell line [in eng]. BMC Cancer. 2012;12:139.
Dean M, Fojo T, Bates S. Tumour stem cells and drug resistance [in eng]. Nat Rev Cancer. 2005;5(4):275–84.
Yang M, Yan M, Zhang R, et al. Side population cells isolated from human osteosarcoma are enriched with tumor-initiating cells [in eng]. Cancer Sci. 2011;102(10):1774–81.
Zhang Y, Yang J, Zhao N, et al. Progress in the chemotherapeutic treatment of osteosarcoma [in eng]. Oncol Lett. 2018;16(5):6228–37.
Valenti MT, Zanatta M, Donatelli L, et al. Ascorbic acid induces either differentiation or apoptosis in MG-63 osteosarcoma lineage [in eng]. Anticancer Res. 2014;34(4):1617–27.
Pautke C, Schieker M, Tischer T, et al. Characterization of osteosarcoma cell lines MG-63, Saos-2 and U-2 OS in comparison to human osteoblasts [in eng]. Anticancer Res. 2004;24(6):3743–8.
Fournier B, Price PA. Characterization of a new human osteosarcoma cell line OHS-4 [in eng]. J Cell Biol. 1991;114(3):577–83.
Hara A, Ikeda T, Nomura S, et al. In vivo implantation of human osteosarcoma cells in nude mice induces bones with human-derived osteoblasts and mouse-derived osteocytes [in eng]. Lab Invest. 1996;75(5):707–17.
Hoang QQ, Sicheri F, Howard AJ, et al. Bone recognition mechanism of porcine osteocalcin from crystal structure [in eng]. Nature. 2003;425(6961):977–80.
Sannino G, Marchetto A, Kirchner T, et al. Epithelial-to-mesenchymal and mesenchymal-to-epithelial transition in mesenchymal tumors: a paradox in sarcomas? [in eng]. Cancer Res. 2017;77(17):4556–61.
Yadav AS, Pandey PR, Butti R, et al. The biology and therapeutic implications of tumor dormancy and reactivation [in eng]. Front Oncol. 2018;8:72.
Weidenfeld K, Barkan D. EMT and stemness in tumor dormancy and outgrowth: are they intertwined processes? [in eng]. Front Oncol. 2018;8:381.
Yang M-H, Wu M-Z, Chiou S-H, et al. Direct regulation of TWIST by HIF-1α promotes metastasis. Nat Cell Biol. 2008;10(3):295–305.
Hoggarth ZE, Osowski DB, Slusser-Nore A, et al. Enrichment of genes associated with squamous differentiation in cancer initiating cells isolated from urothelial cells transformed by the environmental toxicant arsenite [in eng]. Toxicol Appl Pharmacol. 2019;374:41–52.
Kim JC, Yu JH, Cho YK, et al. Expression of SPRR3 is associated with tumor cell proliferation in less advanced stages of breast cancer [in eng]. Breast Cancer Res Treat. 2012;133(3):909–16.
Vos JB, Datson NA, van Kampen AH, et al. A molecular signature of epithelial host defense: comparative gene expression analysis of cultured bronchial epithelial cells and keratinocytes. BMC Genomics. 2006;7(1):9.
Basu-Roy U, Seo E, Ramanathapuram L, et al. Sox2 maintains self renewal of tumor-initiating cells in osteosarcomas [in eng]. Oncogene. 2012;31(18):2270–82.
Zhang L, Zhang L, Xia X, et al. Kruppel-like factor 4 promotes human osteosarcoma growth and metastasis via regulating CRYAB expression [in eng]. Oncotarget. 2016;7(21):30990–1000.
Vermeulen L, de Sousa e Melo F, Richel DJ et al. The developing cancer stem-cell model: clinical challenges and opportunities [in eng]. Lancet Oncol 2012;13(2):e83–e89.
Plaks V, Kong N, Werb Z. The cancer stem cell niche: how essential is the niche in regulating stemness of tumor cells? [in eng]. Cell Stem Cell. 2015;16(3):225–38.
Alfranca A, Martinez-Cruzado L, Tornin J, et al. Bone microenvironment signals in osteosarcoma development [in eng]. Cell Mol Life Sci. 2015;72(16):3097–113.
McQueen P, Ghaffar S, Guo Y, et al. The Wnt signaling pathway: implications for therapy in osteosarcoma [in eng]. Expert Rev Anticancer Ther. 2011;11(8):1223–32.
Ma Y, Ren Y, Han EQ, et al. Inhibition of the Wnt-beta-catenin and Notch signaling pathways sensitizes osteosarcoma cells to chemotherapy [in eng]. Biochem Biophys Res Commun. 2013;431(2):274–9.
Martins-Neves SR, Corver WE, Paiva-Oliveira DI, et al. Osteosarcoma stem cells have active Wnt/beta-catenin and overexpress SOX2 and KLF4 [in eng]. J Cell Physiol. 2016;231(4):876–86.
We thank Ms. Minako Nagata, Ms. Maya Yasuda, and Ms. Kyoko Tanaka for their expert technical assistance. We would like to thank Editage (www.editage.com) for the English language editing.
This work was supported by Grants-in-Aid for Scientific Research (C) (no. 18 K09031) from the Japan Society for the Promotion of Science.
Ethics approval and consent to participate
All animal procedures were performed in accordance with the Japanese Physiological Society Guidelines for the Care and Use of Laboratory Animals, and the study protocol was approved by the Kobe University Animal Experimentation Regulations (Permission no. P180617).
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.
About this article
Cite this article
Fujiwara, S., Kawamoto, T., Kawakami, Y. et al. Acquisition of cancer stem cell properties in osteosarcoma cells by defined factors. Stem Cell Res Ther 11, 429 (2020). https://doi.org/10.1186/s13287-020-01944-9
- Cancer stem cells