This article has Open Peer Review reports available.
In vitro characterization of CD133lo cancer stem cells in Retinoblastoma Y79 cell line
© The Author(s). 2017
Received: 6 May 2017
Accepted: 1 November 2017
Published: 21 November 2017
Retinoblastoma (Rb), the most common childhood intraocular malignant tumor, is reported to have cancer stem cells (CSCs) similar to other tumors. Our previous investigation in primary tumors identified the small sized cells with low CD133 (Prominin-1) and high CD44 (Hyaluronic acid receptor) expression to be putative Rb CSCs using flow cytometry (FSClo/SSClo/CD133lo/CD44hi). With this preliminary data, we have now utilized a comprehensive approach of in vitro characterization of Y79 Rb cell line following CSC enrichment using CD133 surface marker and subsequent validation to confirm the functional properties of CSCs.
The cultured Rb Y79 cells were evaluated for surface markers by flow cytometry and CD133 sorted cells (CD133lo/CD133hi) were compared for CSC characteristics by size/percentage, cell cycle assay, colony formation assay, differentiation, Matrigel transwell invasion assay, cytotoxicity assay, gene expression using microarray and validation by semi-quantitative PCR.
Rb Y79 cell line shared the profile (CD133, CD90, CXCR4 and ABCB1) of primary tumors except for CD44 expression. The CD133lo cells (16.1 ± 0.2%) were FSClo/SSClo, predominantly within the G0/G1 phase, formed larger and higher number of colonies with ability to differentiate to CD133hi cells, exhibited increased invasive potential in a matrigel transwell assay (p < 0.05) and were resistant to Carboplatin treatment (p < 0.001) as compared to CD133hi cells. The CD133lo cells showed higher expression of several embryonic stem cell genes (HOXB2, HOXA9, SALL1, NANOG, OCT4, LEFTY), stem cells/progenitor genes (MSI2, BMI1, PROX1, ABCB1, ABCB5, ABCG2), and metastasis related gene- MACC1, when compared to the CD133hi cells.
This study validates the observation from our earlier primary tumor study that CSC properties in Rb Y79 cell line are endowed within the CD133lo population, evident by their characteristics- i.e. small sized, dormant in nature, increased colony forming ability, differentiation to CD133hi cells, higher invasiveness potential, drug resistance and primitive gene expression pattern. These findings provide a proof of concept for methodological characterization of the retinoblastoma CSCs with future implications for improved diagnostic and treatment strategies.
Retinoblastoma is the most common paediatric ocular malignant tumor occurring in 1 of every 15,000–20,000 live births [1, 2]. This tumor is caused due to inactivation of both the alleles of Retinoblastoma (RB1) gene resulting in the defective pRB protein. RB1 is a major tumor suppressor gene that is involved in cell cycle progression, DNA replication and terminal differentiation . Loss of pRB activity in the retinal progenitor cells leads to impaired cell cycle, uncontrolled cell proliferation and tumor progression. In addition to RB1 as the rate-limiting step for tumor initiation, there are multiple genes (oncogenes and tumor suppressor genes) that undergo mutations, such as MYCN gain, loss of 16q, etc., thereby promoting tumorigenesis [4, 5]. Recent studies have shown that there are cases of unilateral Retinoblastoma that are devoid of Rb mutations and these tumors have distinct histological and genomic landscapes (e.g. high MYCN expression) that facilitate aggressive tumor formation similar to that seen in RB1(−/−) tumors [4, 6]. The hypothesis of cancer stem cells (CSCs), which is now synonymous with tumor initiating cells (TICs) and stem-like cancer cells (SLCCs), originated first from blood related cancers, wherein a small fraction of the tumor cells were reported to be responsible for tumor formation and were attributed with properties of normal stem cells such as quiescence, proliferation, and drug resistance . The salient features of both CSCs and normal stem cells are their potency for self-renewal and forming a cellular hierarchy within the tumor/normal tissue. Additionally, both stem cells and CSCs have the ability to differentiate and migrate . In paediatric brain cancers, tumor derived progenitors form neurospheres that can be passaged at clonal density and are able to self-renew. These cells express several genes characteristic of neural and other stem cells including CD133, NESTIN, SOX2, MSI1, BMI1, MELK, OCT4, etc. [9, 10].
In retinoblastoma, several studies have reported the presence of stem cells in both primary tumors as well as cell lines, using a few properties attributed to cancer stem cell phenotype, ability to actively extrude drugs, slow cycling, clone formation post nutrient starvation, etc. [11–16]. Seigel and co-workers showed the presence of Hoechst dye exclusion, Bromodeoxyuridine (BrdU) label retaining cells, Aldehyde Dehydrogenase 1 (ALDH1) and stem cell markers in the total population of Rb cell lines, primary tumors and simian virus-40 luteinising hormone β sub-unit Large T-antigen (SV40 βLH-T-Ag) mouse tumors [12, 13]. In cultured primary tumors cells, expression of a few retinal development related genes and in vivo tumorigenicity was demonstrated by Zhong et al., thereby hinting at a presence of stem-like cancer cells within Rb tumors . Our group provided evidence of putative stem-like cells in primary Rb tumor cells using a bi-parameter model by flow cytometry with a phenotype of low CD133, high CD44 expression and small sized cells (FSClo/SSClo/CD133lo/CD44hi) expressing progenitor cell markers (PROX1 and SYX1A) . Similar to our study on primary Rb, we observed the two parametric distribution of Rb Y79 cell line with a small subset of cells exhibiting low forward and side scatter profile with low CD133 expression (FSClo/SSClo/CD133lo). It is interesting to note that other than the studies on normal developing retina and Rb deficient retinal cells, which suggests that CD133 expression is low in progenitors and high in differentiated photoreceptors [19, 20], the Rb tumor studies suggest that CD133 expression is specific to cancer stem cells [15, 17, 21]. This prompted us to investigate the surface marker profile of Y79 cell line and compare the functional properties of CSCs  within the CD133 sorted populations.
Retinoblastoma Y79 cells (Riken: RCB1645 Y79 - a generous gift from Dr. S. Krishnakumar, Sankara Nethralaya, Chennai, India) were revived and cultured in Roswell Park Memorial Institute-1640 (RPMI-1640) media supplemented with 10% Fetal Bovine Serum (FBS), antibiotics and L-Glutamine (Gibco™, ThermoFisher Scientific). The cell line was authenticated at the time of purchase from Riken. Regular mycoplasma tests have been carried out and the experiments were conducted in compliance with good laboratory practices (GLP). Media was changed every 3 days and the cells were sub-cultured following observation of cell confluency of about 70%. Enrichment and characterization of putative cancer stem cells were carried out in the cultured Y79 cells and then sorted using the surface marker, CD133. In vitro functional characterization of the sorted subsets was carried out to assess the cell cycle status, clone forming ability, differentiation, invasion, chemoresistance and gene expression signature.
Flow Cytometry analysis and sorting
One million Y79 cells were stained by incubating with directly labelled primary antibodies (CD133-Miltenyi Biotech, CD44, CD90, CXCR4-Ebioscience and ABCB1-Abcam) for 45 min at 4 °C. The antibodies were standardized by varying their dilutions and checking the expression percentage. The cells were then washed thrice with wash buffer to remove excess antibody and run in the BD LSRFortessa™ flow cytometry analyser and the analysis was done using FACSDiva™ software version 6.2. Appropriate controls were used for the experiments. A total of 20,000 to 50,000 events were acquired for analysis. The cells were gated based on size, granularity and doublet discrimination as described previously . In brief, Y79 cells were first selected based on forward and side scatter, and the doublets were excluded using the doublet discrimination plots. The negative control (unstained cells) is used to set the laser voltages, sorting gates and establish the Allophycocyanin (APC) expression profile for the population. The labeled cells were then run through the cytometer and the two populations (CD133hi and lo) are collected in tubes with medium containing 2X antibiotic solution. The post-sort purity and viability was determined. The sorted cells were then used for CSC characterization.
Magnetic activated cell sorting
The Y79 cells in their growth phase were sorted using the CD133 Microbead Kit according to the manufacturer’s protocol (Miltenyi Biotec Inc., Auburn, CA). Briefly, the cells were centrifuged at 300 g for 10 min and resuspended in 300 μL of buffer per 107cells. 100 μL of FcR Blocking Reagent was added and mixed well. To the cell suspension, 100 μL of CD133 Microbeads were added, mixed well and incubated for 30 min at 4 °C. The cells were washed twice with the buffer and resuspended in 500 μL of buffer. Magnetic separation was carried out in LS columns using the MiniMACS™ separator and the two populations were collected in labelled tubes. The sorted cells were then assessed for CD133 expression using flow cytometry for determining the sorting purity and further experiments were carried out following viability count using Trypan blue.
Cell cycle analysis
Sorted CD133hi, CD133lo and unsorted total Y79 cell populations were pelleted by centrifugation and resuspended in PBS with 50 μg/ml propidium iodide for cell cycle analysis. After incubation on ice for 30 min, cell populations were treated with 0.25 mg/ml RNaseA for 45 min at 37 °C to remove RNA. Cells were analysed by flow cytometry at an excitation wavelength of 488 nm and the cell cycle histogram was assessed using the BD FACSDiva software.
Soft agar Colony formation assay and differentiation
The sorted CD133hi and CD133lo Y79 cells were grown in agarose as single cells to assess their colony forming potential. Briefly, a base coat of 0.8% agarose was added into the wells of a 24-well plate and further covered with cell suspension (1000 cells/well in 0.48% agarose). Plates were incubated for 2 weeks following which the resulting colonies were then fixed with 3.7% paraformaldehyde and stained with crystal violet. The images of the colonies were taken at 1.5X and 4X magnification and analysed using ImageJ and OpenCFU software. For colony forming efficiency (CFE) analysis, colonies greater than 50 cells were counted and percentage of colonies were calculated. The morphology of the colonies was assessed for characteristics of holoclones, meroclones and paraclones. The CD133lo clones were then expanded and the expression of CD133 was checked for four passages in vitro.
Matrigel Transwell invasion assay
The Matrigel transwell invasion assays were performed using Corning transwell 24-well inserts with 8 μm pore size as per manufacturer’s instructions. Briefly, 104 cells each of CD133hi and CD133lo subsets were serum-starved overnight and plated into 100 μl serum-free medium onto the inner chamber of the transwell plate, which was previously coated with Matrigel (1:50) (BD Biosciences). The bottom well was then filled with 600 μl of media with 10% serum. The cells were incubated for 24 h at 37 °C and 5% CO2 level. Following incubation, the media was removed from the plate and the non-invasive cells were scraped off from the upper side of the insert and the cells in the lower side of the membrane were fixed in 3.7% Paraformaldehyde (PFA), washed, and stained using crystal violet. Experiments were performed in triplicate transwell and the invaded cells were quantified by counting the average number of cells per 20X field of view and 10 fields per chamber.
The sorted Y79 cells (CD133hi and CD133lo) were assessed for chemoresistance against Carboplatin (Alkem Pharmaceuticals) using MTT assay. Briefly, 5000 cells/90 μL media each of the populations were seeded in a 96 well plate and incubated overnight at 37 °C and 5% CO2 level. Carboplatin was added at varying concentrations (1 μM–100 μM) at a final volume of 10 μL in the wells and incubated for 48 h. Following incubation, 20 μL of 5 mg/mL MTT reagent was added to each well and incubated for 3 h. The Formazan crystals formed were dissolved in 100 μL of Dimethyl sulfoxide (DMSO) and the absorbance was recorded at 595 nm using an ELISA plate reader. The percentage of viability was calculated compared to the controls for each of the population and drug concentration.
Gene expression microarray and pathway analysis
Microarray was performed in duplicates using human whole genome (4x44K) cDNA arrays (Agilent technologies, USA). Labelling reactions were performed using 500 ng of RNA from CD133hi and CD133lo populations. Labelling of the probes was carried out using the Low RNA input linear amplification kit (Agilent technologies, USA) where total RNA was first converted to cDNA using T7-oligo d(T) primers. From this cDNA, labelled cRNA was generated via an in vitro transcription reaction using T7 RNA polymerase and Cy3 (for CD133hi population) or Cy5 (CD133lo population) labelled CTP respectively. Probes with higher labelling efficiency (specific activity ≥8 pmol Cy3 or Cy5/ng cRNA) were selected for competitive hybridization as per the manufacturer’s instructions. 825 ng each of Cy5 and Cy3 labelled cRNAs from CD133lo and CD133hi populations were mixed and added to hybridization buffer and placed on the array. Hybridization was done in a chamber (Agilent technologies, USA) for 17 h at 65 °C with gentle rotation. The slide was scanned and image was analysed using feature extraction tool version 220.127.116.11 (Agilent technologies, USA) and data was analysed using GeneSpring version 10 (Agilent technologies, USA). Lowess algorithm was used to normalize the data. Fold change was calculated based on ratio of Cy5/Cy3 intensities and genes with fold change ≥ + 1.5 or ≤ − 1.5 where considered differentially regulated for which with p-values were assessed. Further, we performed functional enrichment analysis to identify enriched biological processes and pathways in our differentially regulated genes using DAVID bioinformatics resources (version 6.8) and KEGG pathway database [23, 24]. Biological processes and pathways with p-value ≤0.05 where considered significantly enriched. Additionally, we validated the up-regulated genes involved in pathways using Polymerase Chain Reaction (PCR).
Semi quantitative PCR
Primer sequences for the genes used in semi-quantitative PCR
The quantitative data were stated as Mean ± SEM, and GraphPad Prism (GraphPad Software, La Jolla, CA) was used for unpaired Student’s t-test and ANOVA with Bonferroni’s Post-hoc tests. The representative images were analysed using ImageJ software. The experiments were repeated at least thrice with biological replicates and p < 0.05 was considered for statistical significant difference between the groups.
Phenotypic characterization of Y79 cells and CD133 cell sorting
CD133lo cells are predominantly in the resting phase of cell cycle
The cell cycle status was assessed in the two populations using flow cytometry to compare the dormant and proliferative compartments. Cell cycle analysis revealed that the majority of the CD133lo population segregated with resting phase i.e. G0/G1 (83.3 ± 4.1%). On the contrary, 81.1 ± 4.1% of CD133hi cells were in S/G2/M phase as shown in Fig. 1l-n, suggesting that CD133hi cells are mitotically active when compared to dormant CD133lo cells.
CD133lo cells are clonal in nature and differentiate to CD133hi cells in vitro
CD133lo cells display increased invasive ability when compared to CD133hi cells
CD133lo cells exhibit increased chemoresistance to Carboplatin treatment
MTT assay was carried out to assess cytotoxicity of Carboplatin on CD133 sorted Y79 cells following 48 h exposure. CD133lo cells were observed to be more resistant to Carboplatin treatment when compared to CD133hi cells following 48 h exposure (p < 0.001) indicating the phenomenon of chemoresistance in this population (Fig. 3b).
Gene expression microarray and pathway analysis of CD133lo vs CD133hi cells
Deregulated genes in CD133lo CSCs when compared to CD133hi non-CSCs. List of top 30 genes up-regulated in CD133lo cells
Homo sapiens P25 mRNA
Homo sapiens Kruppel like factor 17
hypothetical protein MGC5566
cDNA FLJ43186 fis, clone FCBBF3022767.
phosphodiesterase 4D interacting protein [Homo sapiens (human)]
Human DNA damage repair and recombination protein RAD52 pseudogene mRNA, partial cds
vanin 3 [Homo sapiens (human)]
zinc finger DHHC-type containing 15 [Homo sapiens (human)]
RST23879 Athersys RAGE Library Homo sapiens cDNA
chromosome 1 open reading frame 131 [Homo sapiens (human)]
prostaglandin-endoperoxide synthase 2 [Homo sapiens (human)]
family with sequence similarity 115, member C pseudogene [Homo sapiens (human)]
uncharacterized protein LOC100041774
GTPase IMAP family member 1 [Homo sapiens]
alternative protein PPP1R14C [Homo sapiens]
MDS1 and EVI1 complex locus [Homo sapiens (human)]
CD69 molecule [Homo sapiens (human)]
Homo sapiens mRNA for hSSH-2, complete cds. [AB072358]
MORN repeat containing 5 [Homo sapiens (human)]
ens|cDNA FLJ37906 fis, clone COLON2004318 [Source:UniProtKB/TrEMBL;Acc:Q8N9A9] [ENST00000442408]
ring finger protein 175 [Homo sapiens (human)]
synaptotagmin 4 [Homo sapiens (human)]
serine/threonine kinase 32B [Homo sapiens (human)]
centrosomal protein 126
immunoglobulin heavy variable 3–66
homeobox A11 [Homo sapiens (human)]
activating transcription factor 7 interacting protein 2 [Homo sapiens (human)]
MACC1, MET transcriptional regulator [Homo sapiens (human)]
Deregulated genes in CD133lo CSCs when compared to CD133hi non-CSCs. List of top 30 genes down-regulated in CD133lo cells
follicular lymphoma variant translocation 1
coiled-coil glutamate rich protein 1 [Homo sapiens (human)]
thyroid stimulating hormone receptor [Homo sapiens (human)]
AU146536 HEMBB1 Homo sapiens cDNA clone HEMBB1000770
cadherin 6 [Homo sapiens (human)]
1.00E + 00
G elongation factor mitochondrial 1 [Homo sapiens (human)]
1.00E + 00
distal-less homeobox 2 [Homo sapiens (human)]
receptor tyrosine kinase like orphan receptor 1 [Homo sapiens (human)]
transcription factor EC [Homo sapiens (human)]
melanogenesis associated transcription factor [Homo sapiens (human)]
kininogen 1 [Homo sapiens (human)]
sodium voltage-gated channel alpha subunit 9 [Homo sapiens (human)]
1.00E + 00
1.00E + 00
leiomodin 3 [Homo sapiens (human)]
1.00E + 00
fibroblast growth factor 5 [Homo sapiens (human)]
coiled-coil domain containing 68 [Homo sapiens (human)]
chromosome X open reading frame 57 [Homo sapiens (human)]
Deregulated genes in CD133lo CSCs when compared to CD133hi non-CSCs. Embryonic and Neural stem cell markers up-regulated in CD133lo cells
Genes highly expressed in CD133lo cells
spalt like transcription factor 1
ATP binding cassette subfamily B member 1
ATP binding cassette subfamily B member 5
left-right determination factor
musashi RNA binding protein 2
BMI1 proto-oncogene, polycomb ring finger
Retinoblastoma is a small round cell tumor that comprises of rapidly dividing tumor cells, arising from the retina, within extensive areas of ischemic necrosis as they outgrow their own blood supply . This study attempts to characterize the stem-like cancer cells in the well-established Rb Y79 cell line using the surface marker widely used in other tumors for CSC isolation- CD133 (Prominin). The study confirms that Y79 Rb cell line harbours stem-like cancer cells endowed within the CD133lo enriched population. This was demonstrated by CSC properties such as exclusive surface marker phenotype, size and percentage, slow cycling/dormancy, clone forming ability and differentiation, invasiveness, chemoresistance and primitive gene expression markers.
While in several solid tumors , CD133hi cells are heralded as the cells with CSC properties, there are also contrasting reports of CD133lo subset being capable of exhibiting tumorigenicity and clone forming ability in tumors such as Glioma , Glioblastoma , Colon cancer , etc. Sun and co-workers reported that in neural stem cells, the CD133 negative population are clonogenic and slow cycling in nature . In the developing retina, CD133, which is a conserved surface glycoprotein, is reported to be acquired upon differentiation. In support of this observation is the study by Lakowski et al. who reported that the neuroblastic layer of developing photoreceptors showed lower expression of CD133 when compared to the outer segments of the photoreceptors that had intense expression . Also, Xu and their group showed that in Rb deficient retinal cells, CD133 expression was strong in maturing photoreceptors and weak in the retinal progenitor cell population . Studies have shown that CD133 (Prominin) is crucial for photoreceptor outer segment morphogenesis and the mutations within the gene is associated with several retinal dystrophies which ascertains its pivotal role in the visual cycle . It therefore appears logical to speculate that the CD133 is a marker acquired after cellular differentiation in Retinoblastoma, and the putative CSCs would lack the expression similar to progenitors, since they are believed to arise from the undifferentiated retinal cells.
One of the challenging issues in characterizing CSCs is the sorting/enrichment strategies which are crucial for validating the functional attributes of stem-like cancer cells .With regard to Rb primary tumors, our earlier work highlighted, using a bimodal pattern, that FSClo/SSClo CD133lo/CD44hi cells show more primitive markers as compared to CD133hi cells . In contrast to our observation of CD133lo cells as the Rb CSC subset, Hu et al. using the principle of generating stem-like cancer cells by serum-free culture in WERI-Rb cell line, documented higher percentage of CD133hi cells within the clones when compared to cultured WERI-Rb cells . Based on the evidence of varying expression of CD133 expression in different stages of embryonic stem cell development and its differentiation into neural lineage , it is logical to speculate that there could be a range of expression within tumors of different lineages. Similarly, it could also be affected by modifications induced in primary tumors and the culture conditions of cell lines, which needs to be explored further. This interesting and contrasting expression possibly reflects the evolving concept that the CSC properties of tumor cells enhance or decrease in a gradual fashion and is possibly tissue specific. The surface markers CD90 and CXCR4 are within the range reported by our earlier study in patient samples and correspond to 79.7 ± 1.3% and 14.4 ± 0.5% respectively. This study confirms the presence of CSCs in Y79 cell line comparable to the surface markers observed in the primary tumors , however with one major difference which is the absence of CD44 expression in all cells. Lack of expression of CD44 (Hyaluronic acid receptor) could possibly be attributed to culture conditions and absence of hyaluronic acid, which is abundantly present in the vitreous fluid . This observation is also supported by Ma et al., who showed that long-term serum-free cultures of neurospheres from primary Retinoblastoma showed increased expression of CD44 marker in addition to CD133 when compared to the in vitro differentiated cells . The evidence from retinal studies and increasing expression of CD133 in long-term cultures of Rb supports our observation of CD133 as a marker of differentiation. A thorough comparison of CSC functional properties between the CD133 enriched populations, that was lacking in earlier reports, is highlighted in this study which ascertains the robustness of CD133lo phenotype as a CSC marker.
Another important characteristic of CSCs appears to be the dormancy of cells and segregation within the G0/G1 phase of the cell cycle indicating that they are inherently slow cycling. It has been reported that CSCs are quiescent/slow cycling in nature unless triggered by selection pressure . While the bulk of the tumor population remains in a constant balance between the G1/G2-S-M phases, CSCs are believed to remain in a G0/G1 phase . This attribute makes the CSCs resistant to conventional therapies that target the rapidly dividing tumor cells. In ovarian cancer, Kusumbe et al. showed that the quiescent fraction of the tumor exhibited stem cell activity and ability to revert to a state of self-renewal and differentiation . Our study also demonstrated that CD133lo cells were observed to be in the G0/G1 phase implicating that these are the slow cycling and dormant cells.
The ability to form large colonies within a short period is another hallmark of stem-like cancer cells. Using the soft agar colony forming assay, the study documented that the CD133lo cells were capable of forming larger colonies in soft agar when compared to smaller colonies of the CD133hi subset in 2 weeks of culture. Siegel et al. has shown the evidence for neurosphere forming ability of Y79 cells in culture for 5 days . However, clone forming ability of CD133 sorted populations is one of the novel features of this study. The experiments revealed CD133lo cells formed increased number of large sized colonies as against CD133hi cells in an anchorage independent colony formation assay. This is in concordance with the defined characteristics of holoclones and paraclones in human epidermal keratinocytes by Barrandon and Green and the more recent study on prostate cancer cell clones by Beaver et al. [37, 38]. The CD133lo cells form large clusters of small sized cells with a smooth outline representative of holoclones, whereas CD133hi cells form much smaller colonies with irregular margins depicting paraclones. CD133lo colonies exhibited larger colony area than the CD133hi clusters after 2 weeks hinting at their ability to self-renew and proliferate. The CD133lo population with higher clonogenicity further supports the CSC nature of these cells. The clones could be expanded to passage 5 and showed gradual increase in CD133hi cells with increasing passages suggesting that these cells are capable of self-renewal as well as differentiation to CD133hi cells. In support of this observation is a study by Wang et al. whose group reported that glioma cells that were CD133 negative were capable of forming tumors in nude rats and differentiating into CD133 positive cells . The CSCs generated CD133hi cells upon multiple passages thereby confirming both self-renewal and differentiation. The results from these functional studies of CD133lo vs CD133hi Y79 cells strongly hint that CD133 is a marker of differentiated/mature Rb tumor cells. Validation of this by in vivo tumorigenicity assays would be a valuable study in future.
The invasive potential of these CSCs were evaluated using a matrigel transwell system and the serum starved CD133lo cells exhibited higher invasiveness (p < 0.05) when compared to their counterparts. These observations are in agreement with other tumors such as breast cancer, colon cancer, ovarian cancer, etc., in which the CSCs were shown to have increased metastatic potential which concurs with our findings [39–41].
The CD133lo cells also exhibited chemoresistance towards Carboplatin at higher doses (p < 0.001), which caused cytotoxicity in the CD133hi population at the end of 48 h exposure. CSCs have shown to have increased chemoresistance which helps them to overcome the therapeutic killing by being able to efflux drugs and resist cell death [42, 43]. In the CD133lo population, it was noted that these cells were not only capable of resisting cell death but also exhibited high proliferation in presence of Carboplatin at the end of a 48 h exposure. On the contrary, CD133hi cells showed a profile similar to the total Y79 population with decreasing viability upon exposure to increasing drug concentrations. It is possible to extrapolate that this could be one of the reasons why in some of the cases of Retinoblastoma, the enucleated eyeballs from patients who have received chemotherapy, showed evidence of viable tumor cells with mitotic activity .
Gene expression studies in Y79 cell line revealed several pathways deregulated in the putative CSCs with significant up-regulation of stem cell genes and genes involved in proliferation, chemoresistance and metastasis when compared to the non-CSC population. In support of the evidence in primary tumor cells , where CD44hiCD133lo cells exhibited the progenitor cell markers, the CD133lo cells showed similar marker profile except for CD44 expression. Microarray data also shows that there was significant higher expression of embryonic and neural stem cell markers in the stem-like cancer cells of Y79 cell line. Stem and progenitor cell genes such as OCT4, NANOG, BMI1, ABCG2 and PROX1 were highly expressed in CD133lo population. This is in concordance with the gene expression findings by few groups on total and ABCG2 enriched Rb Y79 cells, and WERI-Rb cells [13, 15, 16]. This study also identified an important regulator of metastasis, Metastasis-associated in colon cancer 1 (MACC1) gene to be expressed in Y79 cell line. To the best of our knowledge, this is the first report of MACC1 overexpression in Retinoblastoma and specifically within the CD133lo subset of Y79 cell line. This gene has been recently identified as an important factor in tumor cell proliferation, invasion and metastasis in several cancers such as Colon, HCC, Breast, Ovarian, etc. via the c-MET/HGF signalling pathway and is being evaluated as a potential therapeutic target [45, 46].
This study also has its set of limitations, the foremost being lack of a positive co-marker for CSCs. Even though the functional studies validate the CD133lo Y79 cells being endowed with CSC properties, which are in tune with the findings in primary tumors; positive co-markers would add more strength to the data. The multidrug resistance markers that were identified in the microarray study (ABCB1, ABCB5) could be used for co-localization studies in future. Though quiescence is one of the hallmarks of CSCs, documentation of the same is a challenge. Cell cycle analysis does not segregate the G0/G1 phases or the different population of cells within the G0 phase of cell cycle, which could also contain differentiated, senescent and stem cells. The study partly addressed it using surface phenotype and clone forming ability, however, segregating the different populations in G0/G1 phase would add value to the study. The CFE measurement was performed for the colonies formed at a 2-week period. Longer incubation period and characterization of the primary clones could further ascertain the clone forming ability and stem cell phenotype. However, the cumulative evidence of all the investigations point towards the CD133lo cells being endowed with CSC properties in Rb Y79 cell line which is a novel finding of this study. We believe that similar CSC characterization studies on other Rb cell lines and primary cells would add more strength to the hypothesis.
This study confirms that the CD133lo subset of Rb Y79 cell line is endowed with the characteristics of cancer stem cells as demonstrated by the in vitro functional assays and is in agreement with our previous findings in primary Rb tumors. This insight of tumor heterogeneity and hierarchy is quite vital for designing better strategies towards a successful therapeutic outcome.
Analysis of variance
Colony forming efficiency
hepatocyte growth factor receptor/hepatocyte growth factor
Cancer stem cells
Cyanine 3/5 dye
Database for annotation, visualization and integrated discovery
Enzyme-linked immunosorbent assay
Fetal Bovine Serum
Good laboratory practices
Kyoto encyclopaedia of genes and genomes
Polymerase chain reaction
Roswell Park Memorial Institute-1640
Standard error of the mean
Stem-like cancer cells
- SV40 βLH-T-Ag:
simian virus-40 luteinising hormone β sub-unit Large T-antigen
Transforming growth factor beta
Tumor initiating cells
Homo sapiens ATP binding cassette subfamily B member 1
Homo sapiens ATP binding cassette subfamily B member 15
Homo sapiens ATP binding cassette subfamily G member 2
Homo sapiens actin beta
Homo sapiens aldehyde dehydrogenase 1 family member A1
Homo sapiens BMI1 proto-oncogene, polycomb ring finger
Homo sapiens prominin 1
Homo sapiens CD44 molecule
Homo sapiens C-X-C motif chemokine receptor 4
Homo sapiens homeobox A9
Homo sapiens homeobox B2
Homo sapiens left-right determination factor
Homo sapiens metastasis associated in colon cancer 1
Homo sapiens maternal embryonic leucine zipper kinase
Homo sapiens musashi RNA binding protein 1
Homo sapiens musashi RNA binding protein 2
Homo sapiens proto-oncogene, bHLH transcription factor
Homo sapiens Nanog homeobox
Homo sapiens nestin
Homo sapiens POU class 5 homeobox 1
Homo sapiens prospero homeobox 1
Homo sapiens RB transcriptional corepressor 1
Homo sapiens spalt like transcription factor 1
Homo sapiens snail family transcriptional repressor 2
Homo sapiens SRY (sex determining region Y)-box 2
Homo sapiens Syntaxin 1A
The authors acknowledge the central instrumentation facilities of L.V. Prasad. Eye Institute and University of Hyderabad, Hyderabad, India.
This study was supported by academic research grants from the Indian Council of Medical Research (ICMR), Hyderabad Eye Research Foundation (HERF), DST-PURSE (University of Hyderabad) and UPE-II (University of Hyderabad). RMN is supported by UGC, Govt. of India (Senior Research Fellowship) and MMSB by CSIR, Govt. of India. The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
All data generated or analysed during this study are included in this published article. The microarray datasets generated and analysed during the current study are available from the corresponding author on reasonable request.
The study was conceptualized by GKV, experiments were designed by RMN, MMSB and GKV. Acquisition of data were accomplished by RMN and MMSB. Microarray experiments and expression data analysis were done by MMSB, IK, RKR and KP. Interpretation of data was done by RMN, MMSB and GKV. All the authors contributed to the preparation, reading and approval of the final manuscript.
Ethics approval and consent to participate
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Deegan WF. Emerging strategies for the treatment of retinoblastoma. Curr Opin Ophthalmol. 2003;14(5):291–5.View ArticlePubMedGoogle Scholar
- DiCiommo D, Gallie BL, Bremner R. Retinoblastoma: the disease, gene and protein provide critical leads to understand cancer. Semin Cancer Biol. 2000;10(4):255–69.View ArticlePubMedGoogle Scholar
- Sachdeva UM, O'Brien JM. Understanding pRb: toward the necessary development of targeted treatments for retinoblastoma. J Clin Invest. 2012;122(2):425–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Theriault BL, Dimaras H, Gallie BL, Corson TW. The genomic landscape of retinoblastoma: a review. Clin Experiment Ophthalmol. 2014;42(1):33–52.View ArticlePubMedGoogle Scholar
- Corson TW, Gallie BL. One hit, two hits, three hits, more? Genomic changes in the development of retinoblastoma. Genes Chromosomes Cancer. 2007;46(7):617–34.View ArticlePubMedGoogle Scholar
- Rushlow DE, Mol BM, Kennett JY, Yee S, Pajovic S, Theriault BL, Prigoda-Lee NL, Spencer C, Dimaras H, Corson TW, et al. Characterisation of retinoblastomas without RB1 mutations: genomic, gene expression, and clinical studies. Lancet Oncol. 2013;14(4):327–34.View ArticlePubMedGoogle Scholar
- Clevers H. The cancer stem cell: premises, promises and challenges. Nat Med. 2011;17(3):313–9.View ArticlePubMedGoogle Scholar
- Magee JA, Piskounova E, Morrison SJ. Cancer stem cells: impact, heterogeneity, and uncertainty. Cancer Cell. 2012;21(3):283–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Singh SK, Clarke ID, Terasaki M, Bonn VE, Hawkins C, Squire J, Dirks PB. Identification of a cancer stem cell in human brain tumors. Cancer Res. 2003;63(18):5821–8.PubMedGoogle Scholar
- Galli R, Binda E, Orfanelli U, Cipelletti B, Gritti A, De Vitis S, Fiocco R, Foroni C, Dimeco F, Vescovi A. Isolation and characterization of tumorigenic, stem-like neural precursors from human glioblastoma. Cancer Res. 2004;64(19):7011–21.View ArticlePubMedGoogle Scholar
- Krishnakumar S, Mallikarjuna K, Desai N, Muthialu A, Venkatesan N, Sundaram A, Khetan V, Shanmugam MP. Multidrug resistant proteins: P-glycoprotein and lung resistance protein expression in retinoblastoma. Br J Ophthalmol. 2004;88(12):1521–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Seigel GM, Campbell LM, Narayan M, Gonzalez-Fernandez F. Cancer stem cell characteristics in retinoblastoma. Mol Vis. 2005;11:729–37.PubMedGoogle Scholar
- Seigel GM, Hackam AS, Ganguly A, Mandell LM, Gonzalez-Fernandez F. Human embryonic and neuronal stem cell markers in retinoblastoma. Mol Vis. 2007;13:823–32.PubMedPubMed CentralGoogle Scholar
- Silva AK, Yi H, Hayes SH, Seigel GM, Hackam AS. Lithium chloride regulates the proliferation of stem-like cells in retinoblastoma cell lines: a potential role for the canonical Wnt signaling pathway. Mol Vis. 2010;16:36–45.PubMedPubMed CentralGoogle Scholar
- Hu H, Deng F, Liu Y, Chen M, Zhang X, Sun X, Dong Z, Liu X, Ge J. Characterization and retinal neuron differentiation of WERI-Rb1 cancer stem cells. Mol Vis. 2012;18:2388–97.PubMedPubMed CentralGoogle Scholar
- Seigel GM, Chen-Tsai RY, Diaz R, Cassidy L. Comparative analysis of ABCG2+ and ABCG2-retinoblastoma cells in three-dimensional culture. Oncocytology J. 2014;4:1.View ArticleGoogle Scholar
- Zhong X, Li Y, Peng F, Huang B, Lin J, Zhang W, Zheng J, Jiang R, Song G, Ge J. Identification of tumorigenic retinal stem-like cells in human solid retinoblastomas. Int J Cancer. 2007;121(10):2125–31.View ArticlePubMedGoogle Scholar
- Balla MM, Vemuganti GK, Kannabiran C, Honavar SG, Murthy R. Phenotypic characterization of retinoblastoma for the presence of putative cancer stem-like cell markers by flow cytometry. Invest Ophthalmol Vis Sci. 2009;50(4):1506–14.View ArticlePubMedGoogle Scholar
- Lakowski J, Han YT, Pearson RA, Gonzalez-Cordero A, West EL, Gualdoni S, Barber AC, Hubank M, Ali RR, Sowden JC. Effective transplantation of photoreceptor precursor cells selected via cell surface antigen expression. Stem Cells. 2011;29(9):1391–404.PubMedPubMed CentralGoogle Scholar
- Xu XL, Singh HP, Wang L, Qi DL, Poulos BK, Abramson DH, Jhanwar SC, Cobrinik D. Rb suppresses human cone-precursor-derived retinoblastoma tumours. Nature. 2014;514(7522):385–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Ma B, Lei X, Guan Y, Mou LS, Yuan YF, Yue H, Lu Y, Xu GT, Qian J. Maintenance of retinal cancer stem cell-like properties through long-term serum-free culture from human retinoblastoma. Oncol Rep. 2011;26(1):135–43.PubMedGoogle Scholar
- Duan JJ, Qiu W, Xu SL, Wang B, Ye XZ, Ping YF, Zhang X, Bian XW, Yu SC. Strategies for isolating and enriching cancer stem cells: well begun is half done. Stem Cells Dev. 2013;22(16):2221–39.View ArticlePubMedPubMed CentralGoogle Scholar
- Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):P3.View ArticlePubMedGoogle Scholar
- Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–61.View ArticlePubMedGoogle Scholar
- Eagle RC Jr. The pathology of ocular cancer. Eye (Lond). 2013;27(2):128–36.View ArticleGoogle Scholar
- Grosse-Gehling P, Fargeas CA, Dittfeld C, Garbe Y, Alison MR, Corbeil D, Kunz-Schughart LA. CD133 as a biomarker for putative cancer stem cells in solid tumours: limitations, problems and challenges. J Pathol. 2013;229(3):355–78.View ArticlePubMedGoogle Scholar
- Wang J, Sakariassen PO, Tsinkalovsky O, Immervoll H, Boe SO, Svendsen A, Prestegarden L, Rosland G, Thorsen F, Stuhr L, et al. CD133 negative glioma cells form tumors in nude rats and give rise to CD133 positive cells. Int J Cancer. 2008;122(4):761–8.View ArticlePubMedGoogle Scholar
- Beier D, Hau P, Proescholdt M, Lohmeier A, Wischhusen J, Oefner PJ, Aigner L, Brawanski A, Bogdahn U, Beier CP. CD133(+) and CD133(−) glioblastoma-derived cancer stem cells show differential growth characteristics and molecular profiles. Cancer Res. 2007;67(9):4010–5.View ArticlePubMedGoogle Scholar
- Shmelkov SV, Butler JM, Hooper AT, Hormigo A, Kushner J, Milde T, St Clair R, Baljevic M, White I, Jin DK, et al. CD133 expression is not restricted to stem cells, and both CD133+ and CD133- metastatic colon cancer cells initiate tumors. J Clin Invest. 2008;118(6):2111–20.PubMedPubMed CentralGoogle Scholar
- Sun Y, Kong W, Falk A, Hu J, Zhou L, Pollard S, Smith A. CD133 (Prominin) negative human neural stem cells are clonogenic and tripotent. PLoS One. 2009;4(5):e5498.View ArticlePubMedPubMed CentralGoogle Scholar
- Gurudev N, Florek M, Corbeil D, Knust E. Prominent role of prominin in the retina. Adv Exp Med Biol. 2013;777:55–71.View ArticlePubMedGoogle Scholar
- Sundberg M, Jansson L, Ketolainen J, Pihlajamaki H, Suuronen R, Skottman H, Inzunza J, Hovatta O, Narkilahti S. CD marker expression profiles of human embryonic stem cells and their neural derivatives, determined using flow-cytometric analysis, reveal a novel CD marker for exclusion of pluripotent stem cells. Stem Cell Res. 2009;2(2):113–24.View ArticlePubMedGoogle Scholar
- Laurie N, Mohan A, McEvoy J, Reed D, Zhang J, Schweers B, Ajioka I, Valentine V, Johnson D, Ellison D, et al. Changes in retinoblastoma cell adhesion associated with optic nerve invasion. Mol Cell Biol. 2009;29(23):6268–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Kreso A, Dick JE. Evolution of the cancer stem cell model. Cell Stem Cell. 2014;14(3):275–91.View ArticlePubMedGoogle Scholar
- Moore N, Lyle S. Quiescent, slow-cycling stem cell populations in cancer: a review of the evidence and discussion of significance. J Oncol. 2011;2011. doi:10.1155/2011/396076.
- Kusumbe AP, Bapat SA. Cancer stem cells and aneuploid populations within developing tumors are the major determinants of tumor dormancy. Cancer Res. 2009;69(24):9245–53.View ArticlePubMedGoogle Scholar
- Barrandon Y, Green H. Three clonal types of keratinocyte with different capacities for multiplication. Proc Natl Acad Sci U S A. 1987;84(8):2302–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Beaver CM, Ahmed A, Masters JR. Clonogenicity: holoclones and meroclones contain stem cells. PLoS One. 2014;9(2):e89834.View ArticlePubMedPubMed CentralGoogle Scholar
- Charafe-Jauffret E, Ginestier C, Iovino F, Wicinski J, Cervera N, Finetti P, Hur MH, Diebel ME, Monville F, Dutcher J, et al. Breast cancer cell lines contain functional cancer stem cells with metastatic capacity and a distinct molecular signature. Cancer Res. 2009;69(4):1302–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Leng Z, Tao K, Xia Q, Tan J, Yue Z, Chen J, Xi H, Li J, Zheng H. Kruppel-like factor 4 acts as an oncogene in colon cancer stem cell-enriched spheroid cells. PLoS One. 2013;8(2):e56082.View ArticlePubMedPubMed CentralGoogle Scholar
- Meng E, Long B, Sullivan P, McClellan S, Finan MA, Reed E, Shevde L, Rocconi RP. CD44+/CD24- ovarian cancer cells demonstrate cancer stem cell properties and correlate to survival. Clin Exp Metastasis. 2012;29(8):939–48.View ArticlePubMedGoogle Scholar
- Alisi A, Cho WC, Locatelli F, Fruci D. Multidrug resistance and cancer stem cells in neuroblastoma and hepatoblastoma. Int J Mol Sci. 2013;14(12):24706–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Khan IN, Al-Karim S, Bora RS, Chaudhary AG, Saini KS. Cancer stem cells: a challenging paradigm for designing targeted drug therapies. Drug Discov Today. 2015;20(10):1205-16.Google Scholar
- Gupta R, Vemuganti GK, Reddy VA, Honavar SG. Histopathologic risk factors in retinoblastoma in India. Arch Pathol Lab Med. 2009;133(8):1210–4.PubMedGoogle Scholar
- Stein U, Walther W, Arlt F, Schwabe H, Smith J, Fichtner I, Birchmeier W, Schlag PM. MACC1, a newly identified key regulator of HGF-MET signaling, predicts colon cancer metastasis. Nat Med. 2009;15(1):59–67.View ArticlePubMedGoogle Scholar
- Wang G, Fu Z, Li D. MACC1 overexpression and survival in solid tumors: a meta-analysis. Tumour Biol. 2015;36(2):1055–65.View ArticlePubMedGoogle Scholar